Nomogram development and validation to predict Ki-67 expression of hepatocellular carcinoma derived from Gd-EOB-DTPA-enhanced MRI combined with T1 mapping

Objective: As an important biomarker to reflect tumor cell proliferation and tumor aggressiveness, Ki-67 is closely related to the high early recurrence rate and poor prognosis, and pretreatment evaluation of Ki-67 expression possibly provides a more accurate prognosis assessment and more better treatment plan. We aimed to develop a nomogram based on gadolinium ethoxybenzyl diethylenetriamine pentaacetic acid (Gd-EOB-DTPA)-enhanced magnetic resonance imaging (MRI) combined with T1 mapping to predict Ki-67 expression in hepatocellular carcinoma (HCC).

Methods: This two-center study retrospectively enrolled 148 consecutive patients who underwent preoperative Gd-EOB-DTPA-enhanced MRI T1 mapping and surgically confirmed HCC from July 2019 to December 2020. The correlation between quantitative parameters from T1 mapping, ADC, and Ki-67 was explored. Three cohorts were constructed: a training cohort (n = 73) and an internal validation cohort (n = 31) from Shunde Hospital of Southern Medical University, and an external validation cohort (n = 44) from the Sixth Affiliated Hospital, South China University of Technology. The clinical variables and MRI qualitative and quantitative parameters associational with Ki-67 expression were analyzed by univariate and multivariate logistic regression analyses. A nomogram was developed based on these associated with Ki-67 expression in the training cohort and validated in the internal and external validation cohorts.

Results: T1rt-Pre and T1rt-20min were strongly positively correlated with Ki-67 (r = 0.627, r = 0.607, P < 0.001); the apparent diffusion coefficient value was moderately negatively correlated with Ki-67 (r = -0.401, P < 0.001). Predictors of Ki-67 expression included in the nomogram were peritumoral enhancement, peritumoral hypointensity, T1rt-20min, and tumor margin, while arterial phase hyperenhancement (APHE) was not a significant predictor even included in the regression model. The nomograms achieved good concordance indices in predicting Ki-67 expression in the training and two validation cohorts (0.919, 0.925, 0.850), respectively.

Conclusions: T1rt-Pre and T1rt-20min had a strong positive correlation with the Ki-67 expression in HCC, and Gd-EOB-DTPA enhanced MRI combined with T1 mapping-based nomogram effectively predicts high Ki-67 expression in HCC.

Introduction

Hepatocellular carcinoma (HCC), as the fifth most common malignancy and the fourth leading cause of cancer-related deaths worldwide, is the most common primary malignant tumor of the liver with increasing incidence (1, 2). The main treatments for HCC include surgical resection, liver transplantation, and transarterial chemoembolization, among others. Surgical resection is recognized as an early radical treatment method, but the recurrence rate of 5 years after surgical resection is as high as 60%–70% (3–5). High early recurrence rate is an important factor affecting the long-term survival and poor prognosis of patients with HCC. The immunohistochemical marker Ki-67 is a nuclear antigen related to cell proliferation activity, and it is a common indicator that reflects the level of cell proliferation. Ki-67 has been proposed as the most valuable independent predictor for evaluating early recurrence and poor prognosis of surgically resected HCC in recent studies (6–8). Ki-67 detection relies on pathological examination. Needle biopsy is a common method to obtain pathological tissue, but it is an invasive examination and has certain disadvantages, such as poor patient compliance, possible surgical risks, and needle tract transfer (9, 10). With the formation of a multidisciplinary and multimethod comprehensive treatment model for liver cancer, a non-invasive preoperative method to predict Ki-67 status is significant in the treatment and prognostic management of patients.

Gadolinium-ethoxybenzyl-diethylenetriamine pentaacetic acid (Gd-EOB-DTPA) is a liver-specific contrast agent, which has dual properties extracellular and hepatobiliary that can be taken up by normal liver cells and is increasingly used in the diagnosis and evaluation of liver diseases. T1 mapping is a non-invasive method of quantitatively analyzing the T1 value of tissue reflecting the T1 relaxation time. T1 mapping combined with GD-EOB-DTPA can provide more accurate and objective magnetic resonance imaging (MRI) quantitative images with functional information. Previous studies have shown that T1 mapping can effectively predict the postoperative histopathological grade and recurrence status of HCC (11, 12) and effectively assess liver fibrosis (13) and liver function (14, 15). These provide a basis for the non-invasive preoperative prediction of Ki-67 expression. To our knowledge, few studies had reported the application of T1 mapping combined with GD-EOB-DTPA for the evaluation of Ki-67 expression in HCC.

Therefore, the purpose of this study was to investigate the correlation between quantitative parameters from T1 mapping, apparent diffusion coefficient (ADC), and Ki-67 expression. We aimed to develop a nomogram based on Gd-EOB-DTPA-enhanced MRI combined with T1 mapping to non-invasive preoperatively predict Ki-67 expression in HCC.

Materials and methods

Patients

This two-center retrospective study was approved by the ethics committee of each participating hospital and was conducted in accordance with the Declaration of Helsinki. The informed consent had a waiver of informed consent. The entire workflow of this retrospective study is shown in Figure 1. A total of 240 patients with surgically and pathologically confirmed HCC from two hospital centers (Shunde Hospital of Southern Medical University, and the Sixth Affiliated Hospital, South China University of Technology) were collected between July 2019 and December 2020. Finally, 148 patients with preoperative Gd-EOB-DTPA-enhanced MRI T1 mapping and clinical data were included in the final analysis, according to the inclusion and exclusion criteria. Three cohorts were constructed: a training cohort (n = 73) and an internal validation cohort (n = 31) from Shunde Hospital of Southern Medical University, and an external validation cohort (n = 44) from the Sixth Affiliated Hospital, South China University of Technology. The inclusion criteria were as follows: (a) pathologically confirmed solitary HCC; (b) underwent hepatectomy; (c) and received preoperative Gd-EOB-DTPA-enhanced MRI and the interval period between the MRI examination and operation in less than 2 weeks. The exclusion criteria were as follows: (a) received previous treatment; (b) incomplete clinical or pathological information; (c) incomplete T1 mapping image data (pre-enhanced and 20 min after enhancement); (d) incomplete MR images or poor quality with obvious artifact; (e) and multiple HCC (≥2).

FIGURE 1

Figure 1 Flowchart of patient selection in the study.

One hundred four patients in the internal cohort were followed up. We followed up patients after surgical resection every 3–6 months by contrast-enhanced ultrasound, CT, or MRI. We defined early recurrence (ER) as recurrence within 1 year after surgery.

Histopathological examination

All tissue specimens were sectioned and stained with hematoxylin eosin, and Ki-67 immunohistochemical staining was performed in the pathology department from two hospitals, respectively. The histopathological and immunohistochemical marker Ki-67 was blinded performed by two experienced pathologists from two hospitals (10, 12, 14, and 16 years’ experience, respectively). In case of disagreements, a third senior pathologist was consulted.

Histological differentiation was assessed by using the Edmondson–Steiner grading system. If different tumor grades coexist within a tumor, the diagnosis is made with the highest grade. Moreover, Edmondson–Steiner grades were divided into two groups; grades 1 and 2 were defined as high differentiation, and grades 3 and 4 were defined as low differentiation (16). Immunohistochemical staining was checked for Ki-67 expression (Beijing Zhongshan Golden Bridge Biotechnology Company, Beijing, China). Interpretation of Ki-67 by immunohistochemistry was as follows: take five FOVs under a high-power microscope (×400), count 100 cells for each FOV, and count the positively stained cells. Ki-67 was identified as the percentage of positive cells to the total cells, and the average value was used and was classified as low Ki-67 expression (≤25%) or high Ki-67 expression (>25%) according to previous studies (6, 17–19) (Figure 2).

FIGURE 2

Figure 2 Representative staining patterns of Ki-67 (×400). (A, B): Low and high Ki-67 expression (approximately 10% and 50%).

MRI examination

All MRI examinations from two hospitals were performed on a 3.0-T MR scanner (MAGNETOM Skyra or MAGNETOM Verio, Siemens Healthcare, Germany). All patients should be fasted for ≥6 h before the examination, and the patients should be trained in breathing before the scan. The examination range was from the upper edge of the liver to the lower edge of the liver. The main scanning sequence includes T1WI-volumetric interpolated breath-hold examination (VIBE), T1WI-TWIST-VIBE sequence, T2WI-BLADE, T2WI-HASTE (half Fourier single-shot turbo spin echo), diffusion-weighted imaging (DWI), and T1 mapping. The parameters are detailed in Table 1.

TABLE 1

Table 1 MRI scan sequence and parameters.

GD-EOB-DTPA (Primovist, Bayer Schering Pharma, Berlin, German) was used for enhanced scanning, the dose was 0.1 ml/kg, the flow rate was 1.0 ml/s, and the tube was flushed with 30 ml of physiological saline. The T1WI-VIBE or T1WI-TWIST-VIBE sequence was used for multiphase (two to five phases) arterial phase scanning 10–30 s after spraying, and the T1WI-VIBE sequence was used for portal vein and balance phase scanning 60 and 150 s after spraying. Hepatobiliary imaging after enhancement was performed 20 min later. T1 mapping pre-enhancement was performed before enhancement, and T1 mapping after enhancement was performed 20 min later.

Imaging analysis

The patients’ image data were exported from the PACS in the DICOM format, and the RadiAnt DICOM Viewer 2020.2 (https://www.radiantviewer.com) software was used for image reading. All MRI quantitative and qualitative features were independently assessed by two abdominal radiologists from Shunde Hospital of Southern Medical University (8 and 15 years’ experience, respectively) who were blinded to the patients’ clinical and pathological information. If there is disagreement on the reevaluated image, it was resolved by consensus.

The following qualitative features were evaluated (Figure 3): (1) tumor margin, including smooth (round or oval with smooth margin) and non-smooth margins (protrusion, depression, or irregular); (2) tumor capsule, defined as the peripheral rim of hyperintensity in the equilibrium phase and classified as complete capsule, and incomplete or non-capsular; (3) mosaic structure, defined as the presence of randomly distributed internal nodules or compartments in the T2-weighted images; (4) arterial phase hyperenhancement (APHE) was defined as enhancement in the arterial phase unequivocally greater in part or in whole than the liver (no non-rim APHE excluded); (5) arterial rim enhancement, defined as the arterial phase enhancement which is most pronounced in the observation periphery; (6) peritumoral enhancement, defined as peri-observation enhancement in the late arterial phase or early portal phase; (7) peritumoral hypointensity, defined as the wedge-shaped area or irregular around the tumor with hypointensity in the HBP (slightly higher signal intensity than the tumor and lower than the surrounding normal liver parenchyma); and (8) satellite nodules, defined as tumors ≤2 cm in size and located ≤2 cm from the main tumor.

FIGURE 3

Figure 3 Qualitative features of hepatocellular carcinoma. Assessment of qualitative features including tumor margin, tumor capsule, mosaic structure, APHE, washout, arterial rim enhancement, peritumoral enhancement, peritumoral hypointensity, and satellite nodules.

The quantitative characteristics included tumor size, T1 relaxation time obtained from T1 mapping pre-enhancement and 20 min after enhancement (T1rt-Pre, T1rt-20min), and ADC value. Selection of the region of interest (ROI) was conducted as follows: (1) select the largest level of measurement as much as possible, and keep the same level of each sequence as much as possible and ROI should be placed as far as possible in the area of obvious enhancement of the lesion; if the tumor is larger than 5 cm, the method of averaging multiple ROIs is used (avoid hemorrhage, cystic area, and blood vessels; Figure 4). (2) The ROIs were not smaller than 1.0 cm² (100–150 pixels); the same lesion was measured three times with the same ROI, and then average amounts were calculated. The average values of the measurements of two radiologists were used for further analyses. (3) The following data were measured: T1 relaxation time of pre-enhancement (T1rt-Pre) and T1 relaxation time 20 min after enhancement (T1rt-20min) on the T1 mapping; calculate the reduction rate of T1 relaxation time (rrT1rt): rrT1rt-20min = (T1rt-Pre-T1rt-20min)/T1rt-Pre; ADC value on the ADC map; and tumor diameter on the maximum diameter of the tumor measured in the transverse or coronal view.

FIGURE 4

Figure 4 (A-D) were ROI measurement methods of T1 relaxation time obtained from pre-enhanced T1 mapping and 20 min after enhancement and the ADC map.

Statistical analysis

Quantitative variables were expressed as mean and standard deviation (SD) (if normally distributed) and median and range (if non-normally distributed), while categorical data were presented as counts and proportions. The independent-sample t-test or Mann–Whitney U-test were used to assess the continuous variables, while the chi-squared test or Fisher’s exact test was used to analyze the categorical variables. Test–retest reliability was assessed using the intraclass correlation coefficient (ICC) for quantitative measurements and kappa for MRI features (>0.75 was considered to represent good agreement). The Spearman correlation analysis method was used to analyze the relationship between quantitative parameters from T1 mapping and ADC and Ki-67 expression in the whole cohort.

Patients from the internal cohort were randomly divided into the training and validation cohorts at a ratio of 7:3. The MRI features and clinical data were compared respectively between the training vs. validation cohorts in the internal cohort and between the Ki-67 high expression group and the low expression group in the internal cohort and the external validation cohort. The variables with statistical differences were determined by univariate logistic regression analysis. Then, the above-identified variables were further selected by stepwise regression based on the Akaike information criterion to construct a multivariate logistic regression model for Ki-67 expression. The nomogram was constructed based on the corresponding final logistic model prediction for these features and was separately validated with the same validation dataset and external dataset. We used Harrell’s C-index and calibration curves in the training and validation cohorts to test the performance of the model. We also calculated the area under the characteristic operating curves and 95% confidence intervals (CIs), accuracy, sensitivity, specificity, positive predictive value, and negative predictive value. Kaplan–Meier analysis was used to compare early recurrence between high Ki-67 expression groups and low Ki-67 expression groups. Statistical significance was considered as a two-sided p value of less than 0.05. The above statistical analysis was performed by the R software version 3.2.3 (Bell Laboratories; https://cran.r-project.org/bin/windows/base/old/3.2.3).

Results

Correlation analysis of quantitative parameters from T1 mapping, ADC and Ki-67 expression, and histological differentiation

T1rt-Pre and T1rt-20min are strongly positively correlated with Ki-67 (r = 0.627, r = 0.607, P < 0.001); the ADC value was moderately negatively correlated with Ki-67 (r = -0.401, P < 0.001). rrT1rt-20min was not significantly associated with Ki-67 (r = -0.150, P = 0.069). The correlation between T1rt-Pre and T1rt-20min with Ki-67 expression in different histological differentiation patients are shown in Figure 5 (all P < 0.05).

FIGURE 5

Figure 5 Correlation between the quantitative parameters of T1 mapping with Ki-67 expression in different histological differentiation patients (*<0.05; ***<0.001).

Baseline characteristics

From July 2019 to December 2020, 104 consecutive patients were initially enrolled from the internal cohort, with 73 patients assigned to the training cohort and 31 to the validation cohort (Table 2). There was no significant difference between the two cohorts in any of the listed variables in Table 2 (all P > 0.05). Among the internal cohort, 50 patients had a high Ki-67 expression (48.1%) and 54 showed a low Ki-67 expression (51.9%). Regarding quantitative data, patients with a high Ki-67 expression had higher alpha-fetoprotein (AFP) levels, larger tumor size, higher T1rt-pre, higher T1rt-20min, and lower ADC (all P < 0.05, Table 3). Regarding qualitative data, patients with high Ki-67 expression more frequently showed APHE, a non-smooth margin, an incomplete or non-complete tumor capsule, a mosaic structure, arterial rim enhancement, peritumoral enhancement, and peritumoral hypointensity (all P < 0.05, Table 3). Among the external validation cohort, 20 patients had a high Ki-67 expression (45.5%) and 24 showed a low Ki-67 expression (54.5%). Patients with a high Ki-67 expression had higher albumin levels, larger tumor size, higher T1rt-pre, higher T1rt-20min, and lower ADC (all P < 0.05, Table 3). Regarding qualitative data, patients with high Ki-67 expression more frequently showed a non-smooth margin, mosaic structure, arterial rim enhancement, peritumoral enhancement, and peritumoral hypointensity (all P < 0.05, Table 3). There were no statistically significant differences in other MR features and clinical indicators between the two groups (P > 0.05). The average concordance rate of evaluating MRI features between the two radiologists was 0.877 (95% CI, 0.855–0.899). MRI images of typical cases in the high and low Ki-67 expression groups are shown in Figure 6.

TABLE 2

Table 2 Baseline characteristics between patients in the internal training and validation cohorts.

TABLE 3

Table 3 Baseline characteristics of patients in the Ki-67 high and low expression groups.

FIGURE 6

Figure 6 Representative images of each phase on Gd-EOB-DTPA-enhanced MRI and T1 mapping 20 min after enhancement in patients with high or low Ki-67 expression. a1–4: A 48-year-old man with Edmondson–Steiner grade II HCC and low expression of Ki-67 (10%). (a1) APHE; (a2), washout and complete capsule; (a3), smooth tumor margin; (a4), T1rt-20 min = 509.92 ms. b1–4: A 58-year-old man with Edmondson–Steiner grade III HCC and high expression of Ki-67 (30%). (b1) no APHE; (b2), peritumoral enhancement; (b3), peritumoral hypointensity and non-smooth tumor margin; (b4), T1rt-20 min = 1,062.86 ms. c1–4: A 59-year-old man with Edmondson–Steiner grade II HCC and high expression of Ki-67 (40%). (c1), rim enhancement and peritumoral enhancement; (c2), washout; (c3), non-smooth tumor margin; (c4), T1rt-20min = 971.84 ms.

Patients in the internal cohort were followed up until recurrence or at the end point of this study (1 February 2022). The median follow-up time of the patients was 18.96 (range 17.07–20.86) months. The overall early recurrence rate was 34.62% (36/104) and in particular 23.32 (21.44–25.19) months for those with low Ki-67 and 13.84 (11.13–16.56) months for those with high Ki-67 (log-rank test, P<0.001; Figure 7).

FIGURE 7

Figure 7 Disease-free survival of patients with high Ki-67 expression and low Ki-67 expression.

Performance of predictors and the combined nomogram for predicting high Ki-67 expression

Using multivariable analyses, we found significant independent predictors of Ki-67 expression, including peritumoral enhancement, peritumoral hypointensity, T1rt-20min, and tumor margin, while APHE was marginally significant predictors of Ki-67 expression (Table 4).

TABLE 4

Table 4 Multivariable analyses of preoperative Gd-EOB-DTPA enhanced MRI features in prediction of Ki-67 expression.

For Ki-67 expression, a nomogram was established in the training cohort based on five imaging features: APHE, peritumoral enhancement, peritumoral hypointensity, T1rt-20min, and tumor margin. The nomograms and calibration curves are presented in Figure 8, which showed that the calibration curves for Ki-67 expression in the three cohorts all well-matched with standard lines. The C-index of the nomogram for Ki-67 expression prediction was 0.919 (95% CI, 0.858–0.970) in the training cohort, 0.925 (95% CI, 0.821–1.000) in the validation cohort, and 0.850 (95% CI, 0.736–0.952) in the external validation cohort (Table 5).

FIGURE 8

Figure 8 Nomogram with calibration curves for predicting Ki-67 expression (internal cohort and external validation cohort).

TABLE 5

Table 5 Predictive performance of preoperative Gd-EOB-DTPA-enhanced MRI features in Ki-67 expression prediction.

Discussion

Gd-EOB-DTPA can be absorbed by the liver cells through the organic anion transporter polypeptide (OATP) pathway, and it provides structural information of lesions and carries functional information in the hepatobiliary phase, due to its dual extracellular and hepatobiliary properties (20). T1 mapping technology is a non-invasive quantitative analysis method of the tissue T1 value, which has good repeatability and strong operability. Gd-EOB-DTPA-enhanced MR combined T1 mapping is believed to provide abundant diagnostic information. In our study, we retrospectively collected data of 148 consecutive patients with HCC, and we investigated the correlation with T1 mapping, ADC quantification parameters, and Ki-67 expression. Our data showed that T1rt-Pre and T1rt-20min were strongly positively correlated with Ki-67 (r = 0.627, r = 0.607, P < 0.001); the ADC value was moderately negatively correlated with Ki-67(r = -0.401, P < 0.001), as seen in previous studies (20). However, the ADC values were less correlated with Ki-67 than T1 mapping; the possible explanation could be that the ADC value can be affected by technical parameters, such as limited image quality with poor signal-to-noise ratio and low spatial resolution, motion and air artifacts, and misregistration artifacts on the ADC map (21). T1 mapping reflects the fixed characteristics of the tissue, is not limited by the scan sequence parameters, and is positively correlated with the concentration of the gadolinium contrast agent in the tissue. Therefore, the T1 mapping value is more accurate, which has also been confirmed in the liver function and liver fibrosis assessment (13–15). The author believes that when the proliferation of tumor cells is more active, the tumor cells are arranged more closely per unit volume, which results in a larger T1rt-pre value; when the proliferation of tumor cells is more active, the proportion of normal liver cells contained in the tumor is lower, and the absorption of Gd-EOB-DTPA is also reduced, resulting in a larger value of T1rt-20min. There is a certain correlation between Ki-67 expression and histopathological differentiation (22, 23). Is the T1 mapping value affected by histological differentiation? Therefore, in our study, the correlation between the T1rt-pre, T1rt-20min, and Ki-67 expressions in the different histological differentiations was further obtained. No matter in the low differentiation group or in the high differentiation group, T1rt-Pre and T1rt-20min had statistical difference between high and low expression of Ki-67 (P < 0.05). It can be seen that T1 mapping is not affected by the degree of histopathological differentiation when predicting high and low expressions of Ki-67.

The expression of Ki-67 in cancer has been intensively studied, and most studies have shown that Ki-67 is associated with the metabolic, genetic, or clinical-pathological features of HCC (5, 19, 20). Aktas’s study showed that Ki-67 was one of the independent prognostic factors of recurrence on patients who underwent liver transplant for HCC (24). Cao et al. revealed that HCC with high Ki-67 expression was more aggressive, and its recurrence-free survival and postoperative overall survival were significantly lower than those of HCC with low Ki-67 expression (6). With the rise of multidisciplinary and multimodal comprehensive treatment of HCC, the high risk of recurrence in HCC patients requires adjuvant therapy and careful follow-up. Hence, a non-invasive preoperative method to predict the Ki-67 status is needed to guide individualized HCC treatment and postoperative surveillance in clinical practice. Many studies have attempted to determine preoperative predictors of HCC with poor prognosis using Ki-67 (6–8), microvascular invasion (MVI) (25, 26), cytokeratin 19 (27, 28), or microvascular density (20). Although the results of these relevant studies have confirmed the relationship between these indicators and prognosis, no consensus has been reached. In terms of the relative molecular mechanism, except for Ki-67 which can be a marker of cellular proliferative activity, few studies have described other mechanisms. Currently, it is more acceptable to use Ki-67 as an important biomarker to reflect tumor cell proliferation and to assess prognosis in patients with HCC and other malignant tumors. However, there is no unified standard for the high and low expression groups of Ki-67 yet. In some studies, 10% (29–31) was used as the cutoff value of Ki-67 expression, but in others, 14% (32) or 35% (7) or 50% (20) was used, and we used 25% (6, 17–19) as the cutoff value because we believe that a cutoff that is too low or too high will lead to selectivity bias for inclusion criteria due to the under- or overestimation of Ki-67 expression. At the same time, our data showed a difference in the overall early recurrence rate between the two groups (P < 0.001). The high Ki-67 expression group (>25%) showed a tendency to recur easily within 1 year after surgery.

In our study, MR qualitative characteristics such as APHE, tumor margin, tumor capsule, mosaic structure, arterial rim enhancement, peritumoral hypointensity, and peritumoral enhancement were statistically different between the high Ki-67 group and low Ki-67 group in the internal cohort (P < 0.05). However, among the above qualitative characteristics, there was no statistical difference in APHE and tumor capsule between the high Ki-67 expression group and low Ki-67 expression group in the external validation cohort. Among these meaningful MR characteristics in our study, peritumoral hypointensity, peritumoral enhancement, tumor capsule, and tumor margin mainly reflected the relatively dense cell structure and infiltrative growth types of the tumor to the surrounding structures. Peritumoral hypointensity and peritumoral enhancement probably relate to the known hypothesis of invasion of the surrounding structures especially the minute portal branch occlusion and the hemodynamic changes existing in compensatory arterial hyperperfusion and decreased portal flow, which may lead to altered expression of OATP or MRP2 receptors (33, 34). However, the statistically different MR characteristics developed in our study were not included in the regression model, such as mosaic structure, arterial rim enhancement, and tumor capsule. Mosaic structure could be related to the inhomogeneity of the tumor, while arterial rim enhancement pattern is rarely observed in HCCs and is more commonly observed in intrahepatic cholangiocarcinomas or metastases (27). It is noteworthy that APHE was statistically different between the two groups (P < 0.05) and included in the nomogram, while it was not a significant predictor of high Ki-67 expression. During hepatocarcinogenesis, changes in the tumor’s blood supply often lead to APHE. However, in some benign lesions, such as atypical hemangioma, it can also present with the typical enhancement pattern of APHE (35). This study did not conduct further analysis of the APHE subtypes, but Cunha’s study (36) showed that there are still differences in observers’ perceptions of APHE subtypes (37), and no unified understanding has been formed. Quantitative assessment can directly reflect histological features. In our study, T1rt-pre and T1rt-20min in the high Ki-67 group were higher than those in the low Ki-67 group (P < 0.001), while the ADC value in the high Ki-67 group was lower than that in the low Ki-67 group (P < 0.001). However, the ADC value was not included in the regression model.

Finally, the significant independent predictors of Ki-67 expression included peritumoral enhancement, peritumoral hypointensity, T1rt-20min, and tumor margin. Although APHE was included in the regression, it was not an independent predictor. For Ki-67 expression, a nomogram was established based on five imaging features: APHE, peritumoral enhancement, peritumoral hypointensity, T1rt-20min, and tumor margin in the training cohort. The combined nomogram yielded an incremental performance in predicting the Ki-67 expression of HCC; the C-index was 0.919 (95% CI, 0.858–0.970) in the training, 0.925 (95% CI, 0.821–1.000) in the validation cohort, and 0.850 (95% CI, 0.736–0.952) in the external validation group.

In summary, the author believes that gadoxetate disodium-enhanced MRI combined with T1 mapping quantitative technology can better evaluate the expression of Ki-67. However, there are several limitations in this study: (1) This is a retrospective study with a relatively small sample size, and there may be selection bias. (2) At present, there is currently no standardized Ki-67 expression level threshold in HCC and there are differences in various studies that we defined 25% as the cutoff value. Therefore, it is necessary to study the differences of different groups with a large sample of multiple centers and compare the diagnostic efficacy of different groups. (3) ROI is drawn manually, and there may be measurement errors. (4) Multiple liver cancers (>2) were excluded from this study. Because this study is a retrospective study, it is impossible to guarantee the one-to-one correspondence between the imaging images of multiple lesions and the pathological images. However, multiple liver cancers are not uncommon in clinical practice. A prospective study was conducted in conjunction with colleagues from hepatobiliary surgery and pathology departments.

Conclusion

In conclusion, this study shows that the T1 relaxation time measured by Gd-EOB-DTPA-enhanced MRI T1 mapping has a strong positive correlation with Ki-67 expression in HCC, and our established nomogram has good predictive performance for a non-invasive preoperative prediction of Ki-67 expression in HCC.

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

This study was reviewed and approved by the ethics committee of each participating hospital. The patients/participants provided their written informed consent to participate in this study.

Author contributions

Acquisiton of the data of the external validation group of our work: JF. All the other authors listed have made substantial, direct and intellectual contribution to the article and approved the submitted version.

Funding

The study was supported by the grants of the science and technology planning project of Foshan (2020001005216) and Guangdong Medical Science and Technology Research Fund (A2021483).

Acknowledgments

The authors thank the patients and their families for their participation in the study.

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.

References

1. Caines A, Selim R, Salgia R. The changing global epidemiology of hepatocellular carcinoma. Clin Liver Dis (2020) 24(4):535–47. doi: 10.1016/j.cld.2020.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol (2019) 16(10):589–604. doi: 10.1038/s41575-019-0186-y

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Renzulli M, Brocchi S, Cucchetti A, Federico M, Cristina M, Camilla S, et al. Can current preoperative imaging be used to detect microvascular invasion of hepatocellular carcinoma? Radiology (2016) 279(2):432–42. doi: 10.1148/radiol.2015150998

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Chen B, Wu JX, Cheng SH, Wang LM, Rong WQ, Wu F, et al. Phase 2 study of adjuvant radiotherapy following narrow-margin hepatectomy in patients with HCC. Hepatology (2021) 74(5):2595–604. doi: 10.1002/hep.31993

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kitamura K, Hatano E, Higashi T, Narita M, Seo S, Nakamoto Y, et al. Proliferative activity in hepatocellular carcinoma is closely correlated with glucose metabolism but not angiogenesis. J Hepatol (2011) 55(4):846–57. doi: 10.1016/j.jhep.2011.01.038

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Cao Y, Ke R, Wang S, Zhu X, Chen J, Huang C, et al. DNA Topoisomerase IIalpha and Ki67 are prognostic factors in patients with hepatocellular carcinoma. Oncol Lett (2017) 13(6):4109–16. doi: 10.3892/ol.2017.5999

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Li HH, Qi LN, Ma L, Chen ZS, Xiang BD, Li LQ. Effect of KI-67 positive cellular index on prognosis after hepatectomy in Barcelona clinic liver cancer stage a and b hepatocellular carcinoma with microvascular invasion. Onco Targets Ther (2018) 11:4747–54. doi: 10.2147/OTT.S165244

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Zhang X, Wu Z, Peng Y, Li D, Jiang Y, Pan F, et al. Correlationship between Ki67, VEGF, and p53 and hepatocellular carcinoma recurrence in liver transplant patients. BioMed Res Int (2021) 2021:6651397. doi: 10.1155/2021/6651397

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Tian G, Kong D, Jiang T, Li L. Complications after percutaneous ultrasound-guided liver biopsy: A systematic review and meta-analysis of a population of more than 12,000 patients from 51 cohort studies. J Ultrasound Med (2020) 39(7):1355–65. doi: 10.1002/jum.15229

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Dube JP, Azzi Z, Semionov A, Sayegh K, Kosiuk J, Presscco J. Imaging of post transthoracic needle biopsy complications. Can Assoc Radiol J (2019) 70(2):156–63. doi: 10.1016/j.carj.2018.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Peng Z, Jiang M, Cai H, Chan T, Dong Z, Luo Y, et al. Gd-EOB-DTPA-enhanced magnetic resonance imaging combined with T1 mapping predicts the degree of differentiation in hepatocellular carcinoma. BMC Cancer (2016) 16:625. doi: 10.1186/s12885-016-2607-4

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Wang WT, Zhu S, Ding Y, Yang L, Chen CZ, Ye QH, et al. T1 mapping on gadoxetic acid-enhanced MR imaging predicts recurrence of hepatocellular carcinoma after hepatectomy. Eur J Radiol (2018) 103:25–31. doi: 10.1016/j.ejrad.2018.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Li Q, Chen T, Shi N, Ye W, Yuan M, Shi Y. Quantitative evaluation of hepatic fibrosis by fibro scan and gd-EOB-DTPA-enhanced T1 mapping magnetic resonance imaging in chronic hepatitis b. Abdom Radiol (NY) (2022) 47(2):684–92. doi: 10.1007/s00261-021-03300-8

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Haimerl M, Fuhrmann I, Poelsterl S, Fellner C, Nickel MD, Weigand K, et al. Gd-EOB-DTPA-enhanced T1 relaxometry for assessment of liver function determined by real-time (13)C-methacetin breath test. Eur Radiol (2018) 28(9):3591–600. doi: 10.1007/s00330-018-5337-y

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Li J, Cao B, Bi X, Chen W, Wang L, Du Z, et al. Evaluation of liver function in patients with chronic hepatitis b using gd-EOB-DTPA-enhanced T1 mapping at different acquisition time points: A feasibility study. Radiol Med (2021) 126(9):1149–58. doi: 10.1007/s11547-021-01382-4

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tang M, Zhou Q, Huang M, Sun K, Wu T, Li X, et al. Nomogram development and validation to predict hepatocellular carcinoma tumor behavior by preoperative gadoxetic acid-enhanced MRI. Eur Radiol (2021) 31(11):8615–27. doi: 10.1007/s00330-021-07941-7

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Bai K, Cao Y, Huang Q, Jiang Y, Lv L. Prognostic value of Ki67 expression for patients with surgically resected hepatocellular carcinoma: Perspectives from a high incidence area. Clin Lab (2017) 63(2):355–64. doi: 10.7754/Clin.Lab.2016.160638

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yao Z, Dong Y, Wu G, Zhang Q, Yang D, Yu J, 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

19. Yang CK, Yu TD, Han CY, Qin W, Liao XW, Yu L, et al. Genome-wide association study of MKI67 expression and its clinical implications in HBV-related hepatocellular carcinoma in southern China. Cell Physiol Biochem (2017) 42(4):1342–57. doi: 10.1159/000478963

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chen Y, Qin X, Long L, Zhang L, Huang Z, Jiang Z, et al. Diagnostic value of gd-EOB-DTPA-Enhanced MRI for the expression of Ki67 and microvascular density in hepatocellular carcinoma. J Magn Reson Imaging (2020) 51(6):1755–63. doi: 10.1002/jmri.26974

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jing M, Cao Y, Zhang P, Zhang B, Lin X, Deng L, et al. The benefit of apparent diffusion coefficient in evaluating the invasiveness of hepatocellular carcinoma. Front Oncol (2021) 11:719480. doi: 10.3389/fonc.2021.719480

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ding Y, Rao S-X, Zhu T, Chen CZ, Li RC, Zeng MS. Liver fibrosis staging using T1 mapping on gadoxetic acid-enhanced MRI compared with DW imaging. Clin Radiol (2015) 70(10):1096–103. doi: 10.1016/j.crad.2015.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ko GH, Go SI, Lee WS, Lee JH, Jeong SH, Lee YJ, et al. Prognostic impact of ki-67 in patients with gastric cancer-the importance of depth of invasion and histologic differentiation. Med (Baltimore) (2017) 96(25):e7181. doi: 10.1097/MD.0000000000007181

CrossRef Full Text | Google Scholar

24. Quaglia A, McStay M, Stoeber K, Loddo M, Caplin M, Fanshawe T, et al. Novel markers of cell kinetics to evaluate progression from cirrhosis to hepatocellular carcinoma. Liver Int (2006) 26(4):424–32. doi: 10.1111/j.1478-3231.2006.01242.x

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Aktas S, Karakayali H, Moray G, Ozdemir H, Haberal M. Effects of risk factors and ki-67 on rates of recurrence on patients who have undergone liver transplant for hepatocellular carcinoma. Transplant Proc (2011) 43(10):3807–12. doi: 10.1016/j.transproceed.2011.09.067

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Erstad DJ, Tanabe KK. Prognostic and therapeutic implications of microvascular invasion in hepatocellular carcinoma. Ann Surg Oncol (2019) 26(5):1474–93. doi: 10.1245/s10434-019-07227-9

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wang Y, Lu L, Liu H, Ding Y, Song N, Lin F, et al. Effect of microvascular invasion on the prognosis in hepatocellular carcinoma and analysis of related risk factors: A two-center study. Front Surg (2021) 8:733343. doi: 10.3389/fsurg.2021.733343

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Choi SY, Kim SH, Park CK, Min JH, Lee JE, Choi YH, et al. Imaging features of gadoxetic acid-enhanced and diffusion-weighted MR imaging for identifying cytokeratin 19-positive hepatocellular carcinoma: A retrospective observational study. Radiology (2018) 286(3):897–908. doi: 10.1148/radiol.2017162846

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Yang F, Wan Y, Xu L, Wu Y, Shen X, Wang J, et al. MRI-Radiomics prediction for cytokeratin 19-positive hepatocellular carcinoma: A multicenter study. Front Oncol (2021) 11:672126. doi: 10.3389/fonc.2021.672126

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Wu H, Han X, Wang Z, Mo L, Liu W, Guo Y, et al. Prediction of the ki-67 marker index in hepatocellular carcinoma based on CT radiomics features. Phys Med Biol (2020) 65(23):235048. doi: 10.1088/1361-6560/abac9c

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Hu XX, Yang ZX, Liang HY, Ding Y, Grimm R, Fu CX, et al. Whole-tumor MRI histogram analyses of hepatocellular carcinoma: Correlations with ki-67 labeling index. J Magn Reson Imaging (2017) 46(2):383–92. doi: 10.1002/jmri.25555

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Li Y, Chen J, Weng S, Sun H, Yan C, Xu X, et al. Small hepatocellular carcinoma: using MRI to predict histological grade and ki-67 expression. Clin Radiol (2019) 74(8):653.e1–e9. doi: 10.1016/j.crad.2019.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Fan Y, Yu Y, Wang X, Hu M, Hu C. Radiomic analysis of gd-EOB-DTPA-enhanced MRI predicts ki-67 expression in hepatocellular carcinoma. BMC Med Imaging (2021) 21(1):100. doi: 10.1186/s12880-021-00633-0

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kim KA, Kim MJ, Jeon HM, Kim KS, Choi JS, Ahn SH, et al. Prediction of microvascular invasion of hepatocellular carcinoma: usefulness of peritumoral hypointensity seen on gadoxetate disodium-enhanced hepatobiliary phase images. J Magn Reson Imaging (2012) 35(3):629–34. doi: 10.1002/jmri.22876

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Hong SB, Choi SH, Kim SY, Shim JH, Lee SS, Byun JH, et al. MRI Features for predicting microvascular invasion of hepatocellular carcinoma: A systematic review and meta-analysis. Liver Cancer (2021) 10(2):94–106. doi: 10.1159/000513704

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Wei Y, Ye Z, Yuan Y, Huang Z, Wei X, Zhang T, et al. A new diagnostic criterion with gadoxetic acid-enhanced MRI may improve the diagnostic performance for hepatocellular carcinoma. Liver Cancer (2020) 9(4):414–25. doi: 10.1159/000505696

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Cunha GM, Hasenstab KA, Delgado T, Ichikawa S, Lee MH, Medina PM, et al. Multi-arterial phase MRI depicts inconsistent arterial phase hyperenhancement (APHE) subtypes in liver observations of patients at risk for hepatocellular carcinoma. Eur Radiol (2021) 31(10):7594–604. doi: 10.1007/s00330-021-07924-8

PubMed Abstract | CrossRef Full Text | Google Scholar