ABSTRACT
PURPOSE
This study aimed to determine the accuracy of texture analysis in differentiating adrenal lesions on unenhanced computed tomography (CT) images.
METHODS
In this single-center retrospective study, 166 adrenal lesions in 140 patients (64 women, 76 men; mean age 56.58 ± 13.65 years) were evaluated between January 2015 and December 2019. The lesions consisted of 54 lipid-rich adrenal adenomas, 37 lipid-poor adrenal adenomas (LPAs), 56 adrenal metastases (ADM), and 19 adrenal pheochromocytomas (APhs). Each adrenal lesion was segmented by manually contouring the borders of the lesion on unenhanced CT images. A texture analysis of the CT images was performed using Local Image Feature Extraction software. First-order and second-order texture parameters were assessed, and 45 features were extracted from each lesion. One-Way analysis of variance with Bonferroni correction and the Mann–Whitney U test was performed to determine the relationships between the texture features and adrenal lesions. Receiver operating characteristic curves were performed for lesion discrimination based on the texture features. Logistic regression analysis was used to generate logistic models, including only the texture parameters with a high-class separation capacity (i.e., P < 0.050). SPSS software was used for all statistical analyses.
RESULTS
First-order and second-order texture parameters were identified as significant factors capable of differentiating among the four lesion types (P < 0.050). The logistic models were evaluated to ascertain the relationships between LPA and ADM, LPA and APh, and ADM and APh. The sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy of the first model (LPA vs. ADM) were 85.7%, 70.3%, 81.3%, 76.4%, and 79.5%, respectively. The sensitivity, specificity, PPV, NPV, and accuracy of the second model (LPA vs. APh) were all 100%. The sensitivity, specificity, PPV, NPV, and accuracy of the third model (ADM vs. APh) were 87.5%, 82%, 36.8%, 98.2%, and 82.7%, respectively.
CONCLUSION
Texture features may help in the characterization of adrenal lesions on unenhanced CT images.
Main points
• Adrenal masses are commonly encountered in daily practice, and it can be difficult to distinguish between lipid-poor adenoma and masses that may require surgical treatment based on non-contrast computed tomography (CT).
• Texture analysis reveals minor differences in gray levels, which can provide guidance in the differential diagnosis of visceral neoplasms.
• This study suggests that lipid-poor adenoma, adrenal metastasis, and pheochromocytoma can be distinguished through an analysis of textural features in non-contrast CT images.
Incidental adrenal masses are common lesions, with a rate of detection of 4%–6% using abdominal computed tomography (CT) in patients undergoing abdominal imaging.1 Adrenal adenomas are the most common adrenal lesions, with a prevalence up to 9% in the general population.2 The rate of adrenal adenomas increases with age.3 Nearly 25% of adrenal adenomas have insufficient intracytoplasmic lipid content for accurate detection using conventional CT, rendering these lesions difficult to differentiate from malignant adrenal lesions and adrenal pheochromocytomas (APhs).3 This can ultimately lead to unnecessary surgical resection of these benign adrenal lesions. Clinical findings and radiological features in unenhanced CT images have a limited ability to differentiate lipid-poor adrenal adenomas (LPAs) from malignant adrenal lesions and APhs because they are in the same density spectrum range. CT washout imaging and magnetic resonance imaging (MRI) are non-invasive imaging techniques used for the differential diagnosis of adrenal lesions.4,5,6
Lipid-rich adrenal adenomas (LRAs) can be easily recognized in unenhanced CT, with homogeneous and relatively low attenuation values [<10 Hounsfield units (HU)], due to their high fat content.7 LRAs can also be diagnosed through MRI with high sensitivity and specificity.1 Most LRAs show a loss of signal intensity (>20%) in out-of-phase images from conventional dual-phase T1-weighted imaging and using the Dixon technique due to intratumoral microscopic fat content.4 LPAs, however, may have a diagnostic overlap with other lipid-poor adrenal lesions, such as adrenal metastasis (ADM) and APhs, in unenhanced CT or MRI images.8 Washout CT of the adrenal gland is the most reliable imaging method for differentiating adrenal adenomas from other lesions. Most adrenal adenomas show a rapid washout pattern after administration of the contrast agent.4 However, hypervascular ADM, especially from hepatocellular carcinoma or renal cell carcinoma, may also show a rapid washout pattern.9 Therefore, the definitive preoperative diagnosis of adrenal lesions is not always possible.
Texture analysis (TA), performed using several complex mathematical processes, is an objective assessment of the structure, gray-level intensity, and pixel position of a lesion.10 The TA parameters of the radiological images can be divided into first-order and second-order statistics. First-order TA parameters consist of the uniformity (in gray-level dispersion), skewness (histogram asymmetry), entropy (irregularity in gray-level dispersion), kurtosis (histogram flatness), energy (uniformity of the distribution), mean/minimum/maximum density, and standard deviation of the gray-level histogram dispersion.11 Second-order TA, used to assess the spatial relationships among pixels, includes the following parameters: gray-level co-occurrence matrix (GLCM), neighborhood gray-level difference matrix (NGLDM), gray-level run-length matrix (GLRLM), gray-level zone length matrix (GLZLM), and their subgroup parameters.12 First-order and second-order TA parameters in adrenal lesions may help physicians to differentiate subgroup lesions based on the underlying histopathological composition. This study investigated the role of TA in differentiating adrenal lesions in unenhanced CT images.
Materials and Methods
Patient selection
The protocol of this retrospective study was approved by our institutional review board (approval number: 2020/20-31; date: 2020-08-31), and written informed consent was obtained from each patient. The endocrinology department database was reviewed to identify patients with LRA, LPA, ADM, or APh between January 2015 and December 2019. The inclusion criteria were age ≥19 years, histopathologically proven APh, diagnosis of LRA, LPA, or ADM based on radiological features, positron emission tomography (PET)–CT features, long-term follow-up, histopathological diagnosis, or biopsy, and unenhanced CT performed before resection or during follow-up screening. Eligible patients had an adrenal lesion determined by unenhanced CT or PET–CT to be ADM (for those with known primary tumors with standardized uptake value (SUV)max ratio >2.5 [the accepted value for distinguishing between benign and cancerous lesions], and SUVmax values of ADM were higher than those of normal liver parenchyma), LRA (≤10 HU in unenhanced CT images), or LPA (diagnosis made by histopathological or radiological follow-up at least 36 months following initial CT).13
Histopathological diagnosis of adrenal lesions is made by surgical intervention, and adrenal lesion biopsy should only be considered for making a diagnosis of metastatic disease in patients with known or suspected non-adrenal cancer.14 The exclusion criteria were adrenal lesion <1 cm (n = 18), lack of unenhanced CT images (n = 82), CT images with severe motion or other artifacts (n = 7), and CT imaging performed at another institution (n = 23). Both lesions of patients with bilateral adrenal lesions were included in the study. We were unable to confirm a diagnosis of adenoma versus ADM histopathologically in three patients, and these patients were also excluded. After these criteria had been applied, a total of 140 patients with 166 adrenal lesions were included (Figure 1). All patients were evaluated through non-contrast abdominal CT using the CT devices available in our department (Supplementary Table 1).
Endocrinological examinations were performed in all cases preoperatively. Metanephrine and normetanephrine were measured in 24-h urine for a differential diagnosis of pheochromocytoma in all patients. Basal plasma levels of adrenocorticotropic hormone and serum levels of cortisol, the dexamethasone suppression test, plasma renin activity, and the aldosterone/renin ratio were used to evaluate adrenal adenoma.
Texture analysis
The unenhanced CT images of each patient were reviewed by three radiologists (HAÖ, IBA, and CA) with 4, 15, and 20 years of experience in abdominal radiology, respectively. The selected images of all adrenal lesions were decided by consensus. For each adrenal lesion, a selected axial image was obtained through CT for use in TA. The selected images for the 166 adrenal masses were anonymized and exported from Sectra IDS7 PACS (Sectra AB, Linköping, Sweden). The image selection was made in the largest diameter of the lesion. Reprocessing of the images was performed using standardization: (gray value-average value)/standard deviation. The images were examined and segmented using Local Image Feature Extraction (LIFEx) software (version 5.1; http://www.lifexsoft.org) independently by two radiologists (IBA and CA) with 15 and 20 years of experience in abdominal radiology, respectively. Regions of interest (ROIs) were drawn manually on each CT image for the measurements (Figures 2, 3). To provide uniformity, the number of gray levels was adjusted to 128 (7 bits). Furthermore, to homogenize the voxel sizes in the XYZ directions, they were adjusted as X, 0.5 mm; Y, 0.5 mm; and Z, 2.5 mm after calculating their mean ± 3 standard deviations. First-order and second-order TA features were calculated in each adrenal lesion.
Statistical Analysis
IBM SPSS software for Windows, version 24.0 (IBM Corporation, Armonk, NY), was used for all statistical analyses. We recorded demographic data, including the patient’s age, sex, diagnosis, and TA measurements. Interobserver variability between the two radiologists was evaluated using intraclass correlation coefficients (0–0.4, poor agreement; 0.41–0.6, moderate agreement; 0.61–0.8, good agreement; 0.81–1, excellent agreement). We performed One-Way analysis of variance (ANOVA) with Bonferroni correction to compare the differences in the TA parameters of LRA, LPA, ADM, and APh. The non-parametric Mann–Whitney U test was performed to detect significant associations between the texture parameters and each adrenal subgroup. LPAs were excluded from the statistical evaluation because they are easily recognized in non-contrast CT examinations and do not cause diagnostic difficulties. Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic performance of the parameters in differentiating among the adrenal lesion subgroups by calculating the area under the ROC curve (AUC). A P value <0.050 was considered statistically significant.
Next, we developed three models to identify texture features predicting LPA, ADM, and APh, with P < 0.050 in the univariate analysis. A logistic regression model was developed based on the qualitative features. In the final logistic regression model, statistically significant texture parameters were used for the three comparisons (12 for LPA vs. ADM, 15 for LPA vs. APh, and 4 for ADM vs. APh). The model performance was assessed according to sensitivity, specificity, negative predictive value (NPV), positive predictive value (PPV), and accuracy.
Results
The final study population consisted of 140 patients (64 men, 76 women; mean age, 56.58 ± 13.65 years) with 166 adrenal lesions. Fifty-four adrenal lesions were diagnosed with LRA, 37 were diagnosed with LPA, 56 were diagnosed with ADM, and 19 were diagnosed with APh (Table 1). There were no significant differences in age or sex among the subgroups. In our series of 91 adrenal adenomas, 76% were non-functioning. The total interobserver agreement between two radiologists (IBA and CA) was κ = 0.73 in the TA calculation. Due to the high level of interobserver agreement and the large amount of numerical data, the data of observer 1 are presented in the tables. Twenty-seven texture parameters, consisting of 11 histogram features, 7 GLCM features, 3 GLRLM features, 3 NGLDM features, and 3 GLZLM features, showed statistically significant differences. The One-Way ANOVA identified 16 second-order texture parameters that were significantly different among all adrenal lesion subgroups (Table 2). The Mann–Whitney U test identified 12 texture parameters that could differentiate between LPA and ADM (P = 0.000–0.032), 15 that could differentiate between LPA and APh (P = 0.000–0.039), and 7 that could differentiate between ADM and APh (P = 0.002–0.042) (Table 3). Texture parameters that were significantly different between LPA and ADM, between LPA and APh, and between ADM and APh were identified by ROC analysis. The AUCs for all independent factors distinguishing the lesions in these three comparisons were between 0.25 and 0.87. The AUCs, 95% confidence intervals, optimal cut-off values, sensitivities, and specificities for each of the three comparisons are listed in Tables 3, 4, 5, respectively. Logistic regression analyses produced three logistic models. The logistic models were evaluated to ascertain the relationships between LPA and ADM, LPA and APh, and ADM and APh (Tables 4, 5, 6). The sensitivity, specificity, PPV, NPV, and accuracy of the first model (LPA vs. ADM) were 85.7%, 70.3%, 81.3%, 76.4%, and 79.5%, respectively (P = 0.050). The sensitivity, specificity, PPV, NPV, and accuracy of the second model (LPA vs. APh) were all 100% (P < 0.001). The sensitivity, specificity, PPV, NPV, and accuracy of the third model (ADM vs. APh) were 87.5%, 82%, 36.8%, 98.2%, and 82.7%, respectively (P = 0.018). We present all values of our logistic regression models in Table 7. The first logistic model found that HUmean, GLRLM-long run low gray-level emphasis (GLRLM-LGE), and NGLDM-Busyness were significant factors differentiating LPA from ADM, with odds ratios (ORs) of 3.87, 2.64, and 5.07, respectively. In the second model, GLRLM-gray-level non-uniformity (GLRLM-GLNU) and HUmean were identified as factors differentiating LPA from APh, with ORs of 119.09 and 4.98, respectively. An HU peak sphere of 0.5 mL was identified as a factor differentiating ADM from APh, with an OR of 4.2.
Discussion
We evaluated the performance of CT-based TA for differentiating adrenal lesions in unenhanced CT images. The results indicate that TA in unenhanced CT images may help to differentiate adrenal lesion subgroups. Numerous TA parameters can be useful for the differential diagnosis of adrenal lesions. GLCM-homogeneity, GLCM-contrast, GLRLM-LGE, and NGLDM-contrast are unique identifiers for discrimination, demonstrating the highest sensitivity and specificity rates (Table 1).
Unenhanced CT, CT washout, and adrenal MRI have been used to differentiate LRAs from other adrenal masses, with high sensitivities and specificities.3,4 However, there are occasionally false-positive diagnoses of LPAs that mimic ADM or APh.15 Adrenal LRAs are easily identified through contrast-enhanced CT using a cut-off of 10 HU.3 Washout CT is widely used to differentiate LRA and LPA from other adrenal lesions. Johnson et al.16 reported a sensitivity of 88% and specificity of 96% for diagnosing adrenal adenomas using absolute percentage washout in CT. However, Patel et al.17 revealed that LPA and APh have a diagnostic overlap with CT attenuation and washout criteria. Seo et al.18 compared the ability of MRI and washout CT to diagnose LPA and reported sensitivities of 75.7% and 100% and specificities of 60% and 80%, respectively. Hypervascular metastases in the adrenal gland show a rapid washout pattern in CT and cannot be differentiated from adrenal adenomas.19 Another diagnostic issue is that hepatocellular carcinomas and renal cell carcinomas are rich in intracytoplasmic fat, resulting in difficult differentiation from LRAs based on chemical-shift MRI.19 Therefore, the diagnostic capabilities of unenhanced CT, washout CT, and adrenal MRI are limited in terms of definitive differentiation among LPA, ADM, and APh.
TA enables the evaluation of a variety of image pixels used to describe the relationships among their gray-level intensity positions within an image.20 To evaluate the internal architecture of adrenal lesions, TA can be used to quantify intratumoral heterogeneity based on the distribution of gray-level values and the spatial arrangement of the pixels within a given region of interest.20 Previous studies have demonstrated that TA helps distinguish among adrenal lesions. Shi et al.1 performed TA using CT to distinguish ADM from benign adrenal masses; they used support vector machine modeling and reported an AUC of 0.85 ± 0.03 and accuracy of 77% for distinguishing metastatic from benign adrenal masses. Yi et al.21 used a texture software program (MaZda, version 4.6) to analyze CT images to differentiate LPA from APh. Their logistic regression model based on four texture parameters differentiated LPA from APh, with a sensitivity, specificity, and accuracy of 86.2%, 97.5%, and 94.4%, respectively. Sensitivity and specificity were higher in the present study; the difference may be related to the differences in the software used and the higher number of second-order parameters in our study. Elmohr et al.22 used TA in 54 cases of adrenal adenomas and carcinomas. Using Boruta random forest modeling, they achieved a validated accuracy of 82% for differentiating benign from malignant adrenal lesions. Ho et al.23 performed a multivariate logistic analysis of 21 texture features, which were combined for each modality, to show that contrast-enhanced CT and chemical-shift MRI could identify malignant lesions, whereas unenhanced CT could not, in 23 patients (malignant adrenal lesions in 8 and LPAs in 15). We included more patients in our study and a significant difference was found between adenomas and malignant lesions. Torresan et al.24 used a cutoff mean densitometry value of 22.5 HU for contrast-enhanced CT TA with a sensitivity and specificity of 95% and 100%, respectively.
To differentiate LRAs from other adrenal lesions (ADMs and APhs) using TA, 10 first-order and 16 second-order texture parameters were identified as statistically significant in our study. HU-Skewness, GLCM-Homogeneity, and NGLDM-Contrast were the most significantly different texture parameters between LRA and ADM. LRAs are easily differentiated from other lipid-poor adrenal lesions using conventional adrenal imaging techniques.
It is often difficult to differentiate LPA from ADM or APh in daily radiological practice. LPA is a benign condition that does not require surgical treatment. In contrast to previous studies, we focused on comparisons of first-order and second-order texture parameters to distinguish among LPA, ADM, and APh.
In our study, LPA was significantly different from ADM in appearance, based on 12 texture parameters in unenhanced CT images, with high sensitivity and specificity of 85.7% and 70.3%, respectively. The first-order texture parameters GLRLM-LGE, GLRLM-GLNU, NGLDM-Busyness, GLZLM low gray-level zone emphasis, and GLZLM-GLNU calculated on enhanced CT images had positive correlations with ADM. Moreover, HUmean (OR: 3.87), GLRLM-LGE (OR: 2.64), and NGLDM-Busyness (OR: 5.07) were the texture parameters that best differentiated LPA from ADM. Our study included a sufficient number of cases to evaluate TA using 45 parameters and to construct a logistic model based on these texture parameters, allowing us to predict malignancy using TA.
We also evaluated the relationship between LPA and APh. A comparison of our results with those in recent studies that evaluated CT for the same clinical questions revealed higher sensitivity and specificity values for our TA results than those identified by Yi et al.21. In this study, sensitivity and specificity for distinction between LPA and APh were 86.2% and 97.5%, respectively. In addition, our model predicts APh with high accuracy using conventional assessment methods in unenhanced CT. Moreover, the GLRLM-GLNU value in the second-order TA was a statistically significant factor differentiating LPA from APh, with a high OR of 119.09.
Finally, we also evaluated the role of TA in differentiating ADM from APh in unenhanced CT images. Our logistic regression analysis identified an HU peak sphere of 0.5 mL (OR: 4.2) as a more robust factor differentiating these two lesion types compared with the other texture parameters. Our study has several limitations. First, the retrospective study design was potentially subject to bias. Previous studies have indicated that the reproducibility and variability of TA are affected by the CT device and technique used. To overcome these limitations, we included 140 patients who underwent unenhanced CT using the same acquisition protocol performed by two different CT scanners. Additional studies involving a larger number of cases and standardized CT acquisition protocols will provide more reliable results. Second, we did not include any primary malignant adrenal tumors. Although malignant adrenal tumors are rare, ADM being the most common, the inclusion of malignant adrenal tumors of different histopathological types, especially adrenocortical cancer, is needed to improve the reliability of our results. Another limitation is that we used two-dimensional analysis, which may have resulted in less available information during the data transfer compared with three-dimensional analysis.25 Finally, our machine learning program used logistic regression analysis and SPSS software, whereas recent studies have used additional texture parameters and specialized machine learning programs (i.e., random forest, GBM), which can improve the results.
In conclusion, TA parameters in unenhanced CT images can be used to differentiate among the four most frequently encountered adrenal tumors, LRA, LPA, ADM, and APh. Differences in texture parameters may be related to differences in tumor cellularity, lipid content, and biological behavior. The addition of TA parameters to CT findings may improve the ability of radiologists to distinguish LPA from ADM or APh. Unenhanced CT-based TA has potential for the diagnosis of incidental adrenal lesions, leading to accurate diagnoses and preventing unnecessary surgical treatment. Further prospective trials in larger populations are needed to verify the role and performance of TA in patients with adrenal masses.