ORIGINAL RESEARCH article

Front. Oncol., 20 January 2022

Sec. Cancer Imaging and Image-directed Interventions

Volume 11 - 2021 | https://doi.org/10.3389/fonc.2021.805545

The Feasibility of Combining ADC Value With Texture Analysis of T2WI, DWI and CE-T1WI to Preoperatively Predict the Expression Levels of Ki-67 and p53 of Endometrial Carcinoma

  • 1. Department of Radiology, Bengbu Medical College, Bengbu, China

  • 2. Department of Radiology, The First Affiliated Hospital of the University of Science and Technology of China, Anhui Provincial Cancer Hospital, Hefei, China

Abstract

Purpose:

To evaluate the feasibility of apparent diffusion coefficient (ADC) value combined with texture analysis (TA) in preoperatively predicting the expression levels of Ki-67 and p53 in endometrial carcinoma (EC) patients.

Methods:

Clinical, pathological and MRI findings of 110 EC patients were analyzed retrospectively. The expression levels of Ki-67 and p53 in EC tissues were detected by immunohistochemistry. ADC value was calculated, and three-dimensional (3D) texture features were measured on T2-weighted images (T2WI), diffusion-weighted images (DWI), and contrast-enhanced T1-weighted images (CE-T1WI). The univariate and multivariate logistic regression and cross-validations were used for the selection of texture features. The receiver operating characteristic (ROC) curve was performed to estimate the diagnostic efficiency of prediction model by the area under the curve (AUC) in the training and validation cohorts.

Results:

Significant differences of the ADC values were found in predicting Ki-67 and p53 (P=0.039, P=0.007). The AUC of the ADC value in predicting the expression levels of Ki-67 and p53 were 0.698, 0.853 and 0.626, 0.702 in the training and validation cohorts. The AUC of the TA model based on T2WI, DWI, CE-T1WI, and ADC value combined with T2WI + DWI + CE-T1WI in the training and validation cohorts for predicting the expression of Ki-67 were 0.741, 0.765, 0.733, 0.922 and 0.688, 0.691, 0.651, 0.938, respectively, and for predicting the expression of p53 were 0.763, 0.805, 0.781, 0.901 and 0.796, 0.713, 0.657, 0.922, respectively.

Conclusion:

ADC values combined with TA are beneficial for predicting the expression levels of Ki-67 and p53 in EC patients before surgery, and they provide higher auxiliary diagnostic values for clinical application.

Introduction

Endometrial carcinoma (EC) is one of the most common malignancies of the female reproductive system worldwide (1), and the morbidity and mortality of EC have been rising with a trend towards a younger age (2). It has been suggested that the occurrence and development of EC are related not only to estrogen levels but also to the proliferation and apoptosis of tumor cells (3). A necessary condition for normal functioning of the body is to maintain the dynamic balance between cell proliferation and apoptosis (4, 5). If this balance is broken, it will promote the occurrence of tumors. Ki-67 and p53 are closely related to the proliferation and apoptosis of tumor cells. Ki-67 is a marker of cell proliferation, is mainly expressed in the nucleus of proliferating cells, and is used to evaluate the proliferation state of tumor cells (6). Higher values of Ki-67 indicate increased malignancy and invasiveness of tumours (7). p53 is an important tumor suppressor gene that controls the initiation of the cell cycle, regulates cell division, inhibits cell growth, regulates transcription, and induces apoptosis (8).

DWI determines the cell density of tissues by detecting the diffusion of water molecules and quantifying it by using the ADC value (9). It has been shown (10) that the ADC mean value is related to the expression level of Ki-67 and p53 in esophageal squamous cell carcinoma, which can be used as a noninvasive biological indicator to predict the proliferation of esophageal squamous cell carcinoma cells and to determine the prognosis of patients. TA is a method to quantitatively measure the distribution and (or) relationship between histogram, pixel intensity or grey level of an image in the region of interest (11, 12). TA highlights subtle patterns of tissue distribution (texture features) that cannot be recognized by human eyes and extends the intrinsic value of images. Thus, in recent years, TA has been used with various cross-sectional imaging modalities and shows clinical applicability in the detection, diagnosis, prognosis, characterization and response evaluation of different cancers (1315). In the present study, five models were developed to preoperatively predict the expression levels of Ki-67 and p53 in EC. The purpose of this study was to noninvasively evaluate the expression levels of Ki-67 and p53 in EC and to provide imaging markers for the clinical diagnosis and treatment of EC.

Materials and Methods

Patients

The Institutional Review Board approved the present retrospective study, and informed consent was waived. Between January 2015 and December 2020, 172 patients who underwent conventional 1.5T MRI before surgery were enrolled by searching our pathology and radiology database in the study. The patients were screened through the medical record system of our hospital. The inclusion criteria were as follows: (a) histologic confirmation of primary EC according to the World Health Organization criteria; (b) no history of preoperative treatment; (c) preoperative available T2WI, DWI and CE-T1WI images; (d) available Ki-67 and p53 expression based on immunohistochemical detection; (e) lesion that could be measured and segmented on MRI. Patients were excluded for the following reasons: (a) incomplete medical records or did not receive treatment in our hospital (n=18); (b) other malignant tumors (n=2); (c) postoperative pathology of non-endometrial carcinoma (n=15); (d) inadequate histopathological reports (n=16); (e) without obvious lesions or maximal tumor diameter of less than 1 cm on MRI (n=11). Finally, 110 patients were included in the present study. The flowchart of the exclusion criteria is shown in Figure 1.

Figure 1

Patients were divided into a training cohort and a validation cohort at a ratio of 7:3. This division was based on the chronological order in which patients were examined to ensure the randomness of the pathological results in the training and validation cohorts. The prediction models of Ki-67 and p53 were built with a training cohort and evaluated with a validation cohort. In the p53 expression of EC, only 38 patients with negative p53 expression and 72 patients with positive p53 expression were enrolled. The ratio of negative p53 expression patients to positive p53 expression patients was about 1:1.89, revealing a sample imbalance. Therefore, the synthetic minority over-sampling technique (SMOTE) algorithm was used to balance the minority class, so that the two classes of EC patients were 1:1(72 negative p53 expression and 72 positive p53 expression).

MR Imaging

All MR imaging studies were performed with 1.5T MR imaging units (United Imaging Healthcare, uMR560, Shanghai) using a 6-channel body array coil. All patients were asked to fast for at least 4 hours, and intramuscular injection of scopolamine butyl bromide was given half an hour before the MRI examination to reduce bowel peristalsis. Detailed scanning parameters are listed as follows: (1) Axial fast spin echo (FSE) T1-weighted images (T1WI): repetition time (TR)/echo time (TE): 391 ms/9.3 ms, slice thickness: 6 mm, inter-slice gap: 2 mm, a field of view (FOV): 38 cm × 28 cm, and matrix size: 320 × 224. (2) Axial FSE T2-weighted images (T2WI): TR/TE: 4500 ms/93 ms, slice thickness: 6 mm, inter-slice gap: 2 mm, FOV: 28 cm × 25 cm, and matrix size: 320 × 224. (3) Axial DWI with b-values of 0 and 1000 s/mm2: TR/TE: 3201 ms/75.4 ms, slice thickness: 6 mm, inter-slice gap: 2 mm, FOV: 38 cm× 28 cm, and matrix size: 128 × 128. (4) Three-dimensional volumetric interpolated quick gradient echo contrast-enhanced imaging on axial: TR/TE: 8 ms/3 ms, slice thickness: 5 mm, inter-slice gap: 0 mm, FOV: 38 cm×28 cm, and matrix size 256×192. The image was exported in DICOM format.

Image Analysis

For the measurement of ADC values, two subspecialty radiologists (reader 1 and reader 2) with 8 and 10 years of experience in imaging diagnosis independently examined the images to obtain ADC values. The original ADC data were directly analyzed in the postprocessing workstation (United Imaging Healthcare). Under the guidance of T2WI, the region of interest (ROI) was delineated in the solid part of the largest tumor layer on the DWI image with b=1000 s/mm2, avoiding hemorrhage and necrotic areas as much as possible. The averages were obtained after the repetition of the measurement three times. The ADC value extraction process is shown in Figure 2.

Figure 2

For the assessment of texture features, the conventional axial T2WI image, DWI with a b value of 1000 s/mm2, and axial CE-T1WI image were segmented by reader 2. The conventional axial T2WI image, DWI with a b value of 1000 s/mm2, and axial CE-T1WI image were imported into ITK-SNAP software (version 3.6.0, http://www.itksnap.org). The three-dimensional region of interest (3D ROI) of the whole tumor was manually delineated at each level of sequence, including areas of hemorrhage and necrosis site. To reduce registration errors, a DWI map with a b value of 1000 s/mm2 was used to assist segmentation. The 3D ROI file was then imported into AK (Analysis Kit, Kinetics Version 2.1, GE Health-care) software to extract texture features. Finally, a total of 828 texture features of the whole tumor were automatically extracted from each of the three MR scanning sequences. As an example, the texture feature extraction process on T2WI is shown in Figure 3. After one week, 30 patients were randomly selected, and the images of 30 patients were re-segmented by reader 1 and reader 2. The intraclass correlation coefficients (ICCs) were used to evaluate the stability of texture features between inter- and intrareader segmentations. The texture features with ICC > 0.75 were preserved.

Figure 3

Pathological Examination

Postoperative Ki-67 and p53 testing were performed by two professional pathologists with more than 8 years of pathological diagnosis experience. The diagnostic criteria of immuno-histochemical staining for p53 included non-staining (0), faint staining (1+), moderate staining (2+) and strong staining (3+). 0 was viewed as negative expression, while 1+ ~ 3+ was viewed as positive expression. The expression of Ki-67 was localized in the nucleus of tumor cells, and five fields were randomly selected under the high-power microscope. The tumor cells were determined to be positive if there were brown particles in the cytoplasm of the tumor cells, and the staining intensity was higher than the background nonspecific staining. According to Li et al. (16), Ki-67 < 50% is the low expression group, and Ki-67 ≥ 50% is the high expression group.

Statistical Analysis

The data were analyzed by SPSS 26.0 (IBM Corporation, Armonk NY, USA), IPMS (Version 2.4.0, GE Health-care) software, and MedCalc Statistical Software version 15.2.2 (MedCalc Software bvba, Ostend, Belgium). Data that had a normal distribution are expressed as the mean ± standard deviation, while abnormally distributed data are expressed as the median. Independent sample t-tests or Mann-Whitney U-tests were used to compare the ADC values and the texture features of each group. P < 0.05 was considered statistically significant. Univariate and multivariate logistic regression analyses and cross-validation were performed on the texture features with statistical significance between each subgroup to select the optimal texture features. The ADC model, TA model, and combined model were built in the training cohorts and validation cohorts to predict the expression levels of Ki-67 and p53. The ROC curves were plotted to assess the performance of the five models in both cohorts. The AUCs were compared using DeLong’s test of equality. P < 0.05 was considered as an independent predictor. A calibration curve was plotted to evaluate the agreement between the prediction result and gold standard.

Results

Clinical and Pathological Findings

Among the 110 EC patients, there were 50 cases of low Ki-67 expression and 60 cases of high Ki-67 expression as well as 38 cases of negative p53 expression and 72 cases of positive p53 expression. The age, FIGO stage, and histological type are presented in Table 1. Regarding Ki-67 expression, there was no significant difference between the two groups in age (P=0.319), but FIGO stage (P=0.020) and histological type (P<0.001) were significant differences. Regarding p53 expression, age (P=0.634) and FIGO stage (P=0.063) were not significantly different between the two groups, but histological type (P=0.030) was statistically significant.

Table 1

Ki-67 (n=110)p53 (n=110)
Patients (n)Low expression (<50%)High expression (≥ 50%)P-valuePatients (n)Negative expression (0)Positive expression (1+ ~ 3+)P-value
Total11050 (45.5%)60 (54.5%)11038 (34.5%)72 (65.5%)
Age53.6 ± 8.3a55.3 ± 9.10.319b54.0 ± 8.654.8 ± 8.90.634
FIGOc stage
I-II9748 (96.0%)49 (81.7%)0.020d9737 (97.4%)60 (83.3%)0.063
III-IV132 (4.0%)11 (18.3%)131 (2.6%)12 (16.7%)
Histologic type
eEMC8948 (96.0%)41 (68.3%)<0.0018935 (92.1%)54 (75.0%)0.030
fNEMC212 (4.0%)19 (31.7%)213 (7.9%)18 (25.0%)

Patients clinical and pathological characteristics.

a

Data are mean ± standard deviation.

b

T-statistical test.

c

FIGO, International Federation of Gynecology and Obstetrics.

d

Chi-square test.

e

EMC, Endometrioid carcinoma.

f

NEMC, Non-endometrioid carcinoma.

ICC Analysis of ADC Values

ICC analysis of measurements of ADC values by the two attending physicians showed good agreement between the surveyors [ICC = 0.883, 95% CI (0.825-0.921), P<0.001]. Therefore, this study used only the results of the first radiologist for a full-text analysis. The ADC values of subgroups were compared by independent sample t-test, and the results are shown in Table 2. The ADC values in EC with high Ki-67 expression were lower than that of the low Ki-67 expression (P =0.007 and P<0.001) with AUCs of 0.698 and 0.853 in the training and validation cohorts, respectively. The ADC value in EC with positive p53 expression was significantly lower than that in EC with negative p53 expression, and the difference was statistically significant (P=0.039 and P=0.048). The AUC of the ADC values in differentiating the level of p53 expression were 0.626 and 0.702 in the training and validation cohorts, respectively.

Table 2

Training cohortValidation cohort
nADC value×10-3mm2/sP-valueAUCnADC value×10-3mm2/sP-valueAUC
low Ki-67 expression340.933 ± 0.1250.0070.698160.974 ± 0.131<0.0010.853
high Ki-67 expression430.844 ± 0.150170.772 ± 0.160
negative p53 expression280.929 ± 0.1530.0390.626100.962 ± 0.1830.0480.702
positive p53 expression490.858 ± 0.137230.830 ± 0.163

The ADC value in relation to Ki-67 and p53 and its predictive performance.

Texture Feature Difference and TA Prediction Model of the Expression Level of Ki-67 and p53 in EC

A total of 828 texture features were obtained from AK software based on T2WI, DWI and CE-T1WI. Sufficient dimension reduction was obtained with the use of univariate and multivariate logistic regression analysis and cross-validation in the training cohort. Three texture features based on T2WI, DWI, and CE-T1WI were extracted from Ki-67 expression, and the detailed information for these texture features is listed in Table 3. Four texture features based on T2WI, DWI, and CE-T1WI were extracted from p53 expression, and the detailed information for these texture features is listed in Table 4. The AUC, sensitivity, and specificity of the models that MR imaging-based TA used to predict the expression level of Ki-67 and p53 in the training and validation cohorts are shown in Table 5.

Table 3

Texture featureslow Ki-67 expressionhigh Ki-67 expressionMultivariate logistic regression analysisAUC
ORP
T2WI- texture features
T2WI-wavelet-HLL_firstorder_Skewness-0.805 ± 0.622-0.399 ± 0.44048.5970.0110.701
T2WI-wavelet-LLL_firstorder_Minimum1125.378 ± 359.7111277.809 ± 309.65724.5020.0240.648
T2WI-wavelet-HHL_glszm_SZN2.867(1.782,3.522)2.833(1.800,5.870)14.5570.0470.530
DWI- texture features
DWI-glcm_Correlation0.212 ± 0.2250.312 ± 0.20214.1340.0200.646
DWI-wavelet-HHL_glszm_HGLZE1.757 ± 0.3331.969 ± 0.36240.9080.0080.662
DWI-wavelet-LHL_firstorder_IR3.006 ± 1.4883.690 ± 1.63349.6230.0170.669
CE-T1WI- texture features
CE-T1WI-wavelet-HLL_glszm_SAE0.315 ± 0.0910.339 ± 0.08820.7630.0420.625
CE-T1WI-wavelet-HLH_firstorder_Kurtosis4.503 ± 1.2123.924 ± 0.6840.0370.0180.647
CE-T1WI-wavelet-LLH_glcm_Correlation0.251 ± 0.1770.344 ± 0.16924.6020.0180.644

Statistical results of texture features of Ki-67 high and low expression groups in EC.

OR, odds ratio; SZN, SizeZoneNonUniformity; HGLZE, HighGrayLevelZoneEmphasis; IR, InterquartileRange; SAE, SmallAreaEmphasis.

Table 4

Texture featuresNegative p53 expressionPositive p53 expressionMultivariate logistic regression analysisAUC
ORP
T2WI- texture features
T2WI -wavelet-LLH_glszm_GLNN0.021 ± 0.0080.026 ± 0.00957.7160.0020.696
T2WI-wavelet-HHH_firstorder_Skewness0.058 ± 0.257-0.054 ± 0.2510.0150.0040.601
T2WI -wavelet-HHL_glszm_GLN8.149 ± 3.7647.402 ± 3.4650.0320.0090.562
T2WI-wavelet-HHL_glcm_ClusterShade-0.004 ± 0.0130.004 ± 0.02256.5450.0230.645
DWI- texture features
DWI -wavelet-HLH_glszm_GLNN0.605 ± 0.0700.558 ± 0.0540.008<0.0010.703
DWI -wavelet-LHH_glszm_GLNN0.520 (0.501,0.571)0.556(0.520,0.625)20.3970.0040.611
DWI -wavelet-LHL_glrlm_LRLGLE14.721 ± 5.22619.409 ± 10.09263.6880.0060.645
DWI -wavelet-LHL_glszm_SZN2.088 (1.617,3.159)1.588(1.000,2.582)0.0350.0140.637
CE-T1WI- texture features
CE-T1WI-firstorder_Maximum138.960 ± 30.666116.564 ± 33.6540.0150.0060.700
CE-T1WI-firstorder_Kurtosis3.265 ± 0.8363.000 ± 0.6790.0120.0030.612
CE-T1WI-glszm_SAE0.219 ± 0.0850.241 ± 0.06738.5690.0090.620
CE-T1WI-glcm_InverseVariance0.435 ± 0.0680.396 ± 0.0530.0480.0220.690

Statistical results of texture features of p53 negative and positive expression groups in EC.

OR, odds ratio; GLNN, GrayLevelNon-UniformityNormalized; GLN, GrayLevelNon-Uniformity; LRLGLE, LongRunLowGrayLevelEmphasis; SZN, SizeZoneNonUniformity; SAE, SmallAreaEmphasis.

Table 5

ModelTraining cohortValidation cohort
AUCa (95%CIb)SENc (%)SPEd (%)AUCa (95%CIb)SENc (%)SPEd (%)
Ki-67T2WI model0.741 (0.629-0.835)53.585.30.688 (0.503-0.837)94.137.5
DWI model0.765 (0.655-0.854)90.755.90.691 (0.507-0.840)52.987.5
CE-T1WI model0.733 (0.619-0.827)72.164.70.651 (0.466-0.808)94.143.8
Combinede model0.922 (0.838-0.971)90.782.40.938 (0.795-0.992)76.5100.0
p53T2WI model0.763 (0.653-0.853)57.185.70.796 (0.620-0.915)65.2100.0
DWI model0.805 (0.699-0.887)71.485.70.713 (0.530-0.856)87.060.0
CE-T1WI model0.781 (0.673-0.868)57.192.90.657 (0.471-0.812)73.970.0
Combinede model0.901 (0.811-0.957)83.785.70.922 (0.773-0.986)100.080.0

The predictive performance of the models in training and validation cohort.

a

AUC, area under the curve;

b

CI, confidence interval;

c

SEN: sensitivity;

d

SPE, specificity;

e

Combined, ADC value combined with T2WI + DWI + CE-T1WI.

Efficacy of Five Models in Predicting the Expression Levels of Ki-67 and p53 in EC

ROC curves were used to evaluate the diagnostic efficacy of prediction models based on ADC values, T2WI, DWI, CE-T1WI, and ADC values combined with T2WI + DWI + CE-T1WI in predicting the expression levels of Ki-67 and p53 in EC. The AUC of ADC values in predicting Ki-67 and p53 expression levels in the training and validation cohorts were 0.698, 0.853 and 0.626, 0.702, respectively. The AUC of the TA model based on T2WI, DWI, CE-T1WI, and ADC value combined with T2WI + DWI + CE-T1WI in the training and validation cohorts for predicting the expression of Ki-67 were 0.741, 0.765, 0.733, 0.922 and 0.688, 0.691, 0.651, 0.938, respectively, and for predicting the expression of p53 were 0.763, 0.805, 0.781, 0.901 and 0.796, 0.713, 0.657, 0.922, respectively, with great performance. DeLong’s test was used to analyze the improvement resulting from the combined model compared to the other single models in the training and validation cohorts (Table 6). The diagnostic efficiency of the combined model was superior to that of the single prediction models based on ADC values, T2WI, DWI and CE-T1WI in the training and validation cohorts. In general, there was no statistical significance among each single model, but the ADC values model was superior to DWI and CE-T1WI in the training group for p53. The ROC curves for the training cohorts are shown in Figure 4. The calibration curves for each training group are shown in Figure 5, and the statistics of the Hosmer–Lemeshow test in each group were not significant (P > 0.05), indicating good calibration.4.

Table 6

Ki-67 (P-value)p53 (P-value)
TrainingValidationTrainingValidation
Combined vs. ADC<0.0010.069<0.0010.031
Combined vs. T2WI<0.0010.0090.0100.114
Combined vs. DWI0.0030.0070.0130.040
Combined vs. CE-T1WI<0.0010.0060.0070.054
ADC vs. T2WI0.6090.1400.0770.450
ADC vs. DWI0.3990.0630.0330.944
ADC vs. CE-T1WI0.6810.1130.0330.746
T2WI vs. DWI0.7710.9750.5850.545
T2WI vs. CE-T1WI0.9000.8020.7990.383
DWI vs. CE-T1WI0.7020.7830.7140.701

AUC of the five models was compared.

Figure 4

Figure 5

Discussion

In recent years, research on tumor biomarkers has increased, aiming to improve the survival rate and quality of life of cancer patients. At the same time, as early detection of EC is critical for treatment, it is important to identify reliable histopathology markers to improve diagnostic accuracy and prognosis. Previous studies have found that Ki-67 and p53 may contribute to the accuracy of cytodiagnosis of different EC lesions (17). However, Ki-67 and p53 are commonly detected using immunohistochemistry, which was not only invasive, but may also be influenced by subjective factors. In the present study, a novel MRI analysis method was proposed to detect the expression levels of Ki-67 and p53 in immunohistochemistry. Not only is MRI noninvasive, but it adds texture features without increasing the patients risk of side effects. Quantitative analysis of the texture of the image may also have an advantage over more targeted biopsies.

The ADC values reflect the diffusion of water molecules in the gap between tumor tissues (18), and they have been shown to be negatively correlated with the proliferation degree and cell density of tumors. The present study suggested that the ADC value in the high Ki-67 expression group was significantly lower than that in the low Ki-67 expression group in EC, which was similar to previous studies in EC (19). With the increase in the expression level of Ki-67, the proliferation activity of EC cells increases, tumors grow vigorously, cell density increases, the arrangement is closer, and the extracellular space decreases, resulting in lower ADC values.

p53 is located on chromosome 17p13-p3 and can be divided into two types, namely, wild type and mutant type, and the pathologically detectable form is mutant p53. When the expression of mutant p53 increases, it may give rise to higher cell proliferative activity, poorer differentiation, and more malignant degree in EC. The present study suggested that the ADC value in the positive p53 group was significantly lower than that in the negative p53 group in EC. Our results were expected and consistent with most previously reported data (10, 20). Although our study showed that the ADC value differed between different subgroups of Ki-67 and p53, the AUC of ADC value for Ki-67 and p53 were 0.698 and 0.626, respectively, which indicated poor performance. This may be because ADC values may not only be related to the movement of water molecules, but also affected by the microcirculation of the blood and tumor heterogeneity; therefore, they cannot truly reflect the movement of water molecules in tissues.

TA is a method that quantifies pixel intensity variations (heterogeneity). TA quantifies tumor heterogeneity by calculating the grey changes of pixels in the image (21). Some studies (22, 23) have found that magnetic resonance imaging–based texture features show the association among deep myometrial invasion, lymphovascular space invasion, and histological high-grade EC. However, only first-order statistical features were extracted and analyzed in these previous studies. In the present study, comprehensive texture features, including first-order statistics, shape-based, GLCM, GLRLM, and GLSZM were extracted based on T2WI, DWI, and CE-T1WI, and they were used to predict the expression levels of Ki-67 and p53 in patients with EC. The texture feature parameters are different in each sequence, and the area under the curve (AUC ≥ 0.7) of the “wavelet-HLL_firstorder_Skewness” of Ki-67 and the “firstorder_Maximum” and “wavelet-HLH_glszm_GLNN” of p53 are more significant. In our study, the Skewness value of the high Ki-67 expression group was higher than that of the low Ki-67 expression group, which indicated that the distribution of the tissue strength grade of the high Ki-67 expression group in EC was more disordered and heterogeneous. Skewness measures the asymmetry of the distribution of values about the mean value. A high skewness value indicates a more asymmetrical the strength grade distribution and greater heterogeneity. The maximum value represents the maximum grey level intensity within the ROI. In the present study, the maximum value of the p53-negative group was higher than that of the p53-positive group, indicating more tumor hemorrhage, secretion, and solid components in the p53-positive group of EC. GLNN measures the variability of grey-level intensity values in the image with a lower value indicating a greater similarity in intensity values. In the present study, the GLNN value of the p53-negative group was lower than that of the p53-positive group, which indicated that the ECs of the p53-positive group were more inhomogeneous and complex than those of the p53-negative group.

Different from previous studies, several studies (24, 25) on TA used one or two MRI sequences. We included T2WI, DWI, and CE-T1WI sequences of EC and established a series of models in which the ADC value and TA predicted the expression levels of Ki-67 and p53. When the ROC curve was used to evaluate the performance of the predictive model, the AUC of the TA model based on T2WI, DWI, and CE-T1WI in the training cohort for predicting the expression of Ki-67 were 0.741, 0.765, and 0.733, respectively, and for predicting the expression of p53 were 0.763, 0.805, and 0.781, respectively. The AUC performed well. The results showed that the TA of T2WI, DWI, and CE-T1WI was helpful to evaluate the expression levels of Ki-67 and p53 in EC. In three sequences of MRI, the AUC of the TA model based on DWI image texture features was higher, which was different from the results reported by Dong et al. (26) who distinguished solitary fibrous tumor/hemangiopericytoma and angiomatous meningioma base on texture feature. One explanation may be that DWI provides a better representation of the microscopic geometry of the EC tissue and the diffusion of water molecules inside and outside the cell and that the images contain more differential texture features with discriminating value. However, there was no significant difference in the TA diagnostic efficiency among the three sequences (T2WI, DWI, and CE-T1WI). It is noteworthy that the diagnostic efficiency of the combined models was better than that of the single model (P < 0.05). The combined model performed well in identifying the expression level of Ki-67 and p53 in EC with AUCs of 0.922 and 0.938 for Ki-67 expression and 0.901 and 0.922 for p53 expression in the training and validation cohorts, respectively. Thus, the present findings suggested that the constructed model of multiparameter MRI TA can capture higher-order interactions between data, reflect tumor heterogeneity from different aspects, and achieve better predictive efficiency. In addition, when using DeLong’s test to evaluate improvements to the combined model, the diagnostic efficiency of the combined models was better than that of the single ADC model, T2WI model, DWI model and CE-T1WI model (P < 0.05). The ADC value improved the diagnostic efficiency of the model, which was consistent with the results of previous studies (27) combining the ADC value with radiomics in EC.

However, there were several limitations in the present study. First, the sample size of this study was small, and it was a retrospective study, indicating potential bias in the selection of enrolled cases. Second, the ADC values were measured based on the maximum level of solid tumor components, which did not represent the overall tumor condition. Third, our model did not include EC FIGO stage, grading, or other clinical indicators. Therefore, future studies will be performed to further increase the sample size.

Conclusions

In conclusion, a combination of ADC values and TA based on three MRI sequences were developed to provide a noninvasive method for preoperatively predicting the expression levels of Ki-67 and p53 in EC. To some extent, this noninvasive imaging marker can compensate for the limitation of endometrium curettage biopsy and adverse impact of tumor heterogeneity, and it provide an objective imaging basis for clinical and accurate individualized treatment.

Funding

This work was supported by 2020 SKY Image Research Fund (No. Z-2014-07-2003-11) and National key research and development program (No. 2016YFB1000905).

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.

Statements

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

Ethics statement

The studies involving human participants were reviewed and approved by the Medical Ethics Committee of Anhui Provincial Cancer Hospital (2021-54). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

JD, XJ, HJ, CW, and CBW contributed to conception and design. XJ, HJ, and ZZ contributed to the collection and arrangement of data. XJ, HJ, CW, CBW, and JD contributed to data analysis and manuscript writing. All authors contributed to the article and approved the submitted version.

Acknowledgments

We are grateful to the technical assistance provided by United Imaging Healthcare.

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.

References

  • 1

    SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

  • 2

    HenleySJWardEMScottSMaJAndersonRNFirthAUet al. Annual Report to the Nation on the Status of Cancer, Part I: National Cancer Statistics. Cancer (2020) 126(10):2225–49. doi: 10.1002/cncr.32802

  • 3

    JiaMJiangPHuangZHuJDengYHuZ. The Combined Ratio of Estrogen, Progesterone, Ki-67, and P53 to Predict the Recurrence of Endometrial Cancer. J Surg Oncol (2020) 122(8):1808–14. doi: 10.1002/jso.26212

  • 4

    LiuYDuFYChenWFuPFYaoMYZhengSS. G15 Sensitizes Epithelial Breast Cancer Cells to Doxorubicin by Preventing Epithelial-Mesenchymal Transition Through Inhibition of GPR30. Am J Transl Res (2015) 7(5):967–75.

  • 5

    KoiCHachisugaTMurakamiMKuritaTNguyenTTShimajiriSet al. Overexpression of P53 in the Endometrial Gland in Postmenopausal Women. Menopause (2015) 22(1):104–7. doi: 10.1097/GME.0000000000000265

  • 6

    LiLTJiangGChenQZhengJN. Ki67 is a Promising Molecular Target in the Diagnosis of Cancer (Review). Mol Med Rep (2015) 11(3):1566–72. doi: 10.3892/mmr.2014.2914

  • 7

    XuQChenCLiuBLinYZhengPZhouDet al. Association of Irhom1 and Irhom2 Expression With Prognosis in Patients With Cervical Cancer and Possible Signaling Pathways. Oncol Rep (2020) 43(1):4154. doi: 10.3892/or.2019.7389

  • 8

    AwaisRSpillerDGWhiteMRParaoanL. P63 is Required Beside P53 for PERP-Mediated Apoptosis in Uveal Melanoma. Br J Cancer (2016) 115(8):983–92. doi: 10.1038/bjc.2016.269

  • 9

    ZhangQYuXLinMXieLZhangMOuyangHet al. Multi-B-Value Diffusion Weighted Imaging for Preoperative Evaluation of Risk Stratification in Early-Stage Endometrial Cancer. Eur J Radiol (2019) 119:108637. doi: 10.1016/j.ejrad.2019.08.006

  • 10

    WangYBaiGZhangXShanWXuLChenW. Correlation Analysis of Apparent Diffusion Coefficient Value and P53 and Ki-67 Expression in Esophageal Squamous Cell Carcinoma. Magn Reson Imaging (2020) 68:183–9. doi: 10.1016/j.mri.2020.01.011

  • 11

    VargheseBACenSYHwangDHDuddalwarVA. Texture Analysis of Imaging: What Radiologists Need to Know. AJR Am J Roentgenol (2019) 212(3):520–8. doi: 10.2214/AJR.18.20624

  • 12

    LubnerMGSmithADSandrasegaranKSahaniDVPickhardtPJ. CT Texture Analysis: Definitions, Applications, Biologic Correlates, and Challenges. Radiographics (2017) 37(5):1483–503. doi: 10.1148/rg.2017170056

  • 13

    YeRWengSLiYYanCChenJZhuYet al. Texture Analysis of Three-Dimensional MRI Images May Differentiate Borderline and Malignant Epithelial Ovarian Tumors. Korean J Radiol (2021) 22(1):106–17. doi: 10.3348/kjr.2020.0121

  • 14

    De CeccoCNGaneshanBCiolinaMRengoMMeinelFGMusioDet al. Texture Analysis as Imaging Biomarker of Tumoral Response to Neoadjuvant Chemoradiotherapy in Rectal Cancer Patients Studied With 3-T Magnetic Resonance. Invest Radiol (2015) 50(4):239–45. doi: 10.1097/RLI.0000000000000116

  • 15

    ChoiMHLeeYJYoonSBChoiJIJungSERhaSEet al. MRI of Pancreatic Ductal Adenocarcinoma: Texture Analysis of T2-Weighted Images for Predicting Long-Term Outcome. Abdom Radiol (NY) (2019) 44(1):122–30. doi: 10.1007/s00261-018-1681-2

  • 16

    LiCZhengMZhengXFangXDongJWangCet al. Predictive Ki-67 Proliferation Index of Cervical Squamous Cell Carcinoma Based on IVIM-DWI Combined With Texture Features. Contrast Media Mol Imaging (2021) 2021:8873065. doi: 10.1155/2021/8873065

  • 17

    ApostolouGApostolouNNikolaidouCKavantzasNPatsourisEAthanassiadouP. Cytodiagnosis of Endometrial Carcinoma and Hyperplasia on Imprint Smears With Additional Immunocytochemistry Using Ki-67 and P53 Biomarkers. Cytopathology (2014) 25(2):8694. doi: 10.1111/cyt.12095

  • 18

    TaouliBKohDM. Diffusion-Weighted MR Imaging of the Liver. Radiology (2010) 254(1):4766. doi: 10.1148/radiol.09090021

  • 19

    JiangJXZhaoJLZhangQQingJFZhangSQZhangYMet al. Endometrial Carcinoma: Diffusion-Weighted Imaging Diagnostic Accuracy and Correlation With Ki-67 Expression. Clin Radiol (2018) 73(4):413.e1–6. doi: 10.1016/j.crad.2017.11.011

  • 20

    WangFWangYZhouYLiuCLiangDXieLet al. Apparent Diffusion Coefficient Histogram Analysis for Assessing Tumor Staging and Detection of Lymph Node Metastasis in Epithelial Ovarian Cancer: Correlation With P53 and Ki-67 Expression. Mol Imaging Biol (2019) 21(4):731–9. doi: 10.1007/s11307-018-1295-7

  • 21

    GaneshanBMilesKA. Quantifying Tumour Heterogeneity With CT. Cancer Imaging (2013) 13(1):140–9. doi: 10.1102/1470-7330.2013.0015

  • 22

    UenoYForghaniBForghaniRDohanAZengXZChamming’sFet al. Endometrial Carcinoma: MR Imaging-Based Texture Model for Preoperative Risk Stratification-A Preliminary Analysis. Radiology (2017) 284(3):748–57. doi: 10.1148/radiol.2017161950

  • 23

    Bereby-KahaneMDautryRMatzner-LoberECornelisFSebbag-SfezDPlaceVet al. Prediction of Tumor Grade and Lymphovascular Space Invasion in Endometrial Adenocarcinoma With MR Imaging-Based Radiomic Analysis. Diagn Interv Imaging (2020) 101(6):401–11. doi: 10.1016/j.diii.2020.01.003

  • 24

    DangMLysackJTWuTMatthewsTWChandaranaSPBrocktonNTet al. MRI Texture Analysis Predicts P53 Status in Head and Neck Squamous Cell Carcinoma. AJNR Am J Neuroradiol (2015) 36(1):166–70. doi: 10.3174/ajnr.A4110

  • 25

    SunXPangPLouLFengQDingZZhouJ. Radiomic Prediction Models for the Level of Ki-67 and P53 in Glioma. J Int Med Res (2020) 48(5):300060520914466. doi: 10.1177/0300060520914466

  • 26

    DongJYuMMiaoYShenHSuiYLiuYet al. Differential Diagnosis of Solitary Fibrous Tumor/Hemangiopericytoma and Angiomatous Meningioma Using Three-Dimensional Magnetic Resonance Imaging Texture Feature Model. BioMed Res Int (2020) 2020:5042356. doi: 10.1155/2020/5042356

  • 27

    ZhangKZhangYFangXFangMShiBDongJet al. Nomograms of Combining Apparent Diffusion Coefficient Value and Radiomics for Preoperative Risk Evaluation in Endometrial Carcinoma. Front Oncol (2021) 11:705456. doi: 10.3389/fonc.2021.705456

Summary

Keywords

endometrial carcinoma, p53, Ki-67, apparent diffusion coefficient, texture analysis

Citation

Jiang X, Jia H, Zhang Z, Wei C, Wang C and Dong J (2022) The Feasibility of Combining ADC Value With Texture Analysis of T2WI, DWI and CE-T1WI to Preoperatively Predict the Expression Levels of Ki-67 and p53 of Endometrial Carcinoma. Front. Oncol. 11:805545. doi: 10.3389/fonc.2021.805545

Received

01 November 2021

Accepted

29 December 2021

Published

20 January 2022

Volume

11 - 2021

Edited by

Chuanming Li, Chongqing Medical University, China

Reviewed by

Caterina Giannitto, Humanitas Research Hospital, Italy; Jian Shu, Affiliated Hospital of Southwest Medical University, China

Updates

Copyright

*Correspondence: Jiangning Dong,

†These authors have contributed equally to this work

This article was submitted to Cancer Imaging and Image-directed Interventions, a section of the journal Frontiers in Oncology

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics