- Research
- Open access
- Published:
CT-based radiomics nomogram for differentiation of adrenal hyperplasia from lipid-poor adenoma: an exploratory study
BMC Medical Imaging volumeĀ 23, ArticleĀ number:Ā 4 (2023)
Abstract
Background
To establish and verify a radiomics nomogram for differentiating isolated micronodular adrenal hyperplasia (iMAD) from lipid-poor adenoma (LPA) based on computed tomography (CT)-extracted radiomic features.
Methods
A total of 148 patients with iMAD or LPA were divided into three cohorts: a training cohort (nā=ā72; 37 iMAD and 35 LPA), a validation cohort (nā=ā36; 22 iMAD and 14 LPA), and an external validation cohort (nā=ā40; 20 iMAD and 20 LPA). Radiomics features were extracted from contrast-enhanced and non-contrast CT images. The least absolute shrinkage and selection operator (LASSO) method was applied to develop a triphasic radiomics model and unenhanced radiomics model using reproducible radiomics features. A clinical model was constructed using certain laboratory variables and CT findings. Radiomics nomogram was established by selected radiomics signature and clinical factors. Nomogram performance was assessed by calibration curve, the areas under receiver operating characteristic curves (AUC), and decision curve analysis (DCA).
Results
Eleven and eight extracted features were finally selected to construct an unenhanced radiomics model and a triphasic radiomics model, respectively. There was no significant difference in AUC between the two models in the external validation cohort (0.838 vs. 0.843, pā=ā0.949). The radiomics nomogram inclusive of the unenhanced model, maximum diameter, and aldosterone showed the AUC of 0.951, 0.938, and 0.893 for the training, validation, and external validation cohorts, respectively. The nomogram showed good calibration, and the DCA demonstrated the superiority of the nomogram compared with the clinical factors model alone in terms of clinical usefulness.
Conclusions
A radiomics nomogram based on unenhanced CT images and clinical variables showed favorable performance for distinguishing iMAD from LPA. In addition, an efficient unenhanced model can help avoid extra contrast-enhanced scanning and radiation risk.
Background
An adrenal adenoma is the most common type ofĀ adrenal gland tumor. It is a benign tumor derived from the adrenal cortical tissue [1], usually presenting as a small rounded or oval mass (often measuringā<ā3Ā cm) [2]. Due to the existence of intrinsic fat content [3], more than 70% adrenal adenomas can be evidentially identified using a threshold ofāā¤ā10 Hounsfield units (Hu) on unenhanced computed tomography (CT), while those withā>ā10 Hu are classified into lipid-poor adenoma (LPA) [4,5,6].
Adrenal hyperplasia is a genetic disorder and physiologic overgrowth of adrenocortical tissue that affects the adrenal glands. It can be nodular or diffuse. Diffuse type manifests as a diffuse and smooth increase in the size of adrenal glands [6]; the nodular forms can be further classified into the primary pigmented micronodular adrenal type and the isolated micronodular hyperplasia (iMAD) type [7]. iMAD is characterized by a hyperplasia of the internodular tissue [8]. Patients with iMAD can be easily misdiagnosed with other high-attenuation diseases, such as LPA [3, 9]. The diagnosis of this condition is mainly based on conventional clinic and radiological methods, which are not sensitive enough to differentiate iMAD from LPA.
Patients with adrenal adenoma and adrenal hyperplasia can have similar clinical symptoms, including primary hypertension, primary aldosteronism, or Cushing's syndrome [10, 11]. Thus, differentiating these conditions may be challenging. While adrenal hyperplasia is resolved using a medication, patients with adrenal adenoma usually require surgery, i.e., adrenalectomy [12]. Adrenalectomy lowers blood pressure levels [13], ameliorates hypokalemia, reduces plasma aldosterone levels [14], and drastically increases survival in patients without atrial fibrillation [15]. On the other hand, the use of adrenalectomy in the management of adrenal hyperplasia is controversial. Surgery has been associated with certain side effects in these patients, including abnormal circadian rhythm [16] and corticotrophin deficiency [17]; also, these patients have a higher risk of having anesthesia accidents and extra medical expenses. Thus, accurate preoperative diagnosis is of utmost importance for therapy decision-making [18]. However, diagnosing based on conventional clinic and radiological methods remains challenging. Hence, a high-performance diagnostic approach needs to be established and certified.
Radiomics is a relatively new quantitative approach for medical imaging that extracts many features (e.g., intensity, geometry, and texture) from images using data-characterization algorithms.Ā Over the years, this approach has been mainly applied in oncology but also in some other medical areas [19, 20]. Successful studies on adrenal tumors have been reported using radiomics to distinguish between benign adrenal tumors with malignant tumors [21] and pheochromocytoma with the other benign tumors [22, 23]. Yet, to our knowledge, no study have reported CT-based radiomics to distinguish iMAD from LPA. Thus, in this study, we developed and validated a radiomics nomogram incorporating the radiomics signature and the clinical factors for preoperative differential diagnosis of iMAD and LPA.
Methods
Patients
Data for all cases were acquired from two medical centers: center 1 and center 2, from May 2018 to May 2022. Five hundred and sixty-eight patients who accepted unilateral adrenalectomy were included in the study, and they were all diagnosed with adrenal hyperplasia or adenoma pathologically. All patients underwent contrast-enhanced CT before surgery. CT images of all patients were further analyzed to identify iMAD or LPA using the following inclusion criteria: (1) patients with a single nodular lesion; (2) CT attenuation of lesions areā>ā10 Hu on unenhanced CT images. The exclusion criteria were: (1) patients with lesionsāā¤ā10 Hu on unenhanced CT images; (2) bilateral or diffuse lesions; (3) lesions with necrosis or hemorrhage based on the pathological diagnosis; (4) patients without complete preoperative clinical and imaging data before surgery. The patient recruitment pathway is presented in Fig.Ā 1.
Finally, 148 patients were recruited: 108 patients from center 1 were divided into two cohorts (training and testing cohort) according to the ratio of 2:1 using computer-generated random numbers, and 40 patients from center 2 were incorporated into the external validation cohort.
The institutional review board of our hospital approved this retrospective study. Therefore, the requirement for obtaining informed consent was waived.
CT image acquisition
In center 1, patients underwent CT scanning using multidetector-row CT systems (Ingenuity CT, Philips; Aquilion ONE, TOSHIBA; Somatom Force, Siemens Healthcare; Lispeed 64, GE Healthcare). The scanning parameters were as follows: 120Ā kV tube voltage; 250ā400Ā mA tube current (automatic tube current modulation was applied); 0.5Ā s or 0.6Ā s rotation time; 192āĆā0.6Ā mm or 64āĆā0.625Ā mm detector collimation; a matrix of 512āĆā512, and a pitch of 0.6 or 1.0. Axial images of 1Ā mm slice thickness were reconstructed. An 80ā90Ā mL volume of iodinated contrast medium (Omnipaque 350, GE Healthcare, Shanghai, China) was injected via antecubital vein by a power injector (at a rate of 3.0Ā mL/s). Pre-enhanced CT images were acquired before two postcontrast CT scannning in the arterial phase (25ā30Ā s) and the venous phase (60ā70Ā s).
In center 2, patients were examined using the multidetector-row CT systems (Somatom Force, Siemens Healthcare; Discovery 750, GE Healthcare). The scanning parameters were as follows: 100 or 120Ā kV tube voltage; 0.5 or 0.6Ā s rotation time; 250ā400Ā mA tube current (automatic tube current modulation was applied); 192āĆā0.6Ā mm or 64āĆā0.625Ā mm detector collimation; a matrix of 512āĆā512, and a pitch of 1. Axial images of 1Ā mm slice thickness were reconstructed. An 80ā90Ā mL volume of iodinated contrast medium (Iopromide, Ultravist 300; Bayer, Germany) was injected via antecubital vein using a power injector (at a rate of 2.5Ā mL/s). Unenhanced CT images were acquired before two postcontrast CT scanning in the arterial phase (25ā30Ā s) and venous phase (55ā60Ā s).
Conventional CT feature evaluation
All of the image analysis were performed by the same radiology resident (Reader 1, HY) and a radiologist (Reader 2, CJ) with 5 and 10Ā years of abdominal imaging experience, respectively. The radiologists were blinded to the radiological reports and pathologic details. Reader 1 construed the following CT features by consensus: the maximum diameter (MD) and the CT attenuation of unenhanced phase (CTpre) and venous phase (CTV) on the axial CT image. All quantitative values were measured 3 times, and an average figure was applied.
Clinical factors selection and construction of the clinical factor model
Clinical data included specific assay indexes (cortisol, aldosterone, and renin) and data regarding age, sex, and history of hypertension. All information were obtained from medical records. Univariate logistic analysis was used to compare the differences in the clinical factors (including clinical data and CT features) in the training cohort between the two groups; a multiple logistic regression analysis was applied to build the clinical model by using the significant variables from the univariate analysis as inputs. Odds ratios (OR) with 95% confidence intervals (CI) for each independent factor were calculated as relative risk estimates.
Three-dimensional segmentation and radiomics feature extraction
Image resampling was performed before feature extraction to decrease the variability of radiomics features. Images were resampled to 1āĆā1āĆā1 mm3 voxels using the B-Spline interpolation method. A soft tissue CT window was modified with a level of 40HU and a width of 300 Hu. Three-dimensional lesion segmentation was manually delineated as a region of interest (ROI) on images of three phases at ITK-SNAP (http://www.itksnap.org/pmwiki/pmwiki.php). Images of corresponding arterial and venous phases were used as reference when the contour on the unenhanced image was delineated. An example of manual segmentation is shown in Fig.Ā 2.
Radiomics extraction was performed with Deepwise Research Portal (Deepwise Bolian, Co., Ltd.), integrated with PyRadiomics. A total of 1288 radiomics features were extracted from each phase. Twenty patients (10 iMAD and 10 LPA) were randomly chosen for further evaluation of feature stability: both the two radiologists repeated the same delineation procedure two weeks later. Inter- and intra- class correlation coefficients (ICC) were used to assess the extracted featuresā inter-observer reliability and intra-observer reproducibility. It is commonly accepted that intra- and inter-class correlation coefficient (ICC)ā<ā0.5 indicates poor reliability, 0.5ā0.75: moderate reliability, andā>ā0.75: good or excellent reliability [24]. Thus, the features with ICCā<ā0.75 were excluded.
Construction and comparison of the radiomics signatures
Z-score normalization was first used to standardize the features before the following analysis. Dimension reduction was conducted before signature construction in order to eliminate features overfitting. The stable radiomics features were assessed by one-way analysis of variance, and the features that were significantly different between the two groups (pā<ā0.05) were enrolled into the least absolute shrinkage and selection operator (LASSO) regression model in order to select the most valuable features in the training cohort. Then, the selected features were applied to build radiomics signatures. Two radiomics models were established: a triphasic CT feature based on unenhanced and enhanced CT features and an unenhanced radiomics model based on unenhanced CT features. A radiomics score (Rad-score) was calculated for each patient through a linear combination of selected features weighted by their respective LASSO coefficients.
Development of a radiomics nomogram and assessment of the performance of different models
Delongās test was first used to compare two radiomics models before the radiomics nomogram was performed. The triphasic radiomics model was chosen when establishing a nomogram if there was a significant difference between two radiomics signatures; otherwise, an unenhanced radiomics model was chosen. Then, a radiomics nomogram was establish by incorporating the chosen radiomics model and the selected clinical factors. A radiomics nomogram score (Nomo-score) was calculated for each patient. The calibration of the nomogram was verified using a calibration curve. The HosmerāLemeshow test was used to validate the goodness-of-fit of the nomogram. The diagnostic performance of the clinical model, the chosen radiomics signature, and the radiomics nomogram was assessed based on the area under the receiver operator characteristic curve (AUC) in the training, testing, and external validation sets. Decision curve analysis (DCA) was performed to evaluate the clinical usefulness of three models by calculating the net benefits for a range of threshold probabilities in the external validation cohort [21].
Statistical analysis
Statistical tests were performed using MedCalc statistical software (version 20.110, https://www.medcalc.org), R statistical software (version 3.3.3, https://www.r-project.org), and MATLAB (version 2021a). Studentās t-tests or non-parametric tests (where appropriate) were used for continuous variables. Chi-squared test or Fisherās exact test (where appropriate) were used for categorical variables. Categorical and continuous variables are shown in frequency (percentages), meanāĀ±āstandard deviation, or median (interquartile range), where appropriate. The AUC was compared by Delongās test. A p valueā<ā0.05 was considered to be statistically significant.
Results
Clinical baseline and construction of the clinical model
The patientsā demographic baseline characteristics are shown in Table 1. The rates of LPA in the training, testing, and external validation cohort were 48.6% (35 of 72), 38.9% (14 of 36), and 50.0% (20 of 40), respectively. The results of univariate and multiple logistic regression analysis are listed in Table 2. Multiple logistic regression showed that maximum diameter and aldosterone levels are independent predictors of LPA (all pā<ā0.05), thus were chosen to create a clinical model. Lesions with larger maximum diameter (OR, 1.296; 95% CI, 1.036ā1.621) or higher aldosterone (OR, 1.005; 95% CI, 1.000 -1.010) were likely to be LPA.
Feature extraction, selection, and radiomics signature establishment and comparison
A total of 710 unenhanced CT radiomics features, 863 arterial CT features and 707 venous CT features were considered robust features (ICCsā>ā0.75). Sixteen triphasic features and eleven unenhanced features were selected by LASSO logistic regression model to construct the triphasic radiomics model and the unenhanced radiomics model, respectively (Fig.Ā 3).
Feature selection and cross-validation using LASSO. For different values of lambda, different numbers of features will be included in the final model and a decision has to be made. Figures A, C illustrate the process of cross validation in order to choose the most appropriate lambda value. Lambda resulting in the smallest mean squared error (MSE) is chosen, as the green lines shows. Figures B, D plot different values of lambda on the x-axis and each line represents a feature, showing when it enters the model and its level of influence on the outcome. A Selection of the tuning parameter (lambda) in the LASSO of unenhanced model. B The unenhanced CT feature coefficients varied according to lambda. C Selection of the tuning parameter (lambda) in the LASSO of enhanced model. D The enhanced CT feature coefficients varied according to lambda
For the unenhanced model, radiomics score [median (interquartile range)] differed significantly between iMAD and LAP groups in the training cohort [āā0.11 (āā0.15,āāā0.03) vs. 0.10 (0.04, 0.15), respectively, pā<ā0.001]; this was further verified in the validation cohort [āā0.07 (āā0.12, 0,03) vs. 0.07 (0.01, 0.19), respectively, pā<ā0.001] and the external validation cohort [āā1.74 (āā4.14,āāā0.38) vs. 0.65 (0.22, 4.54) respectively, pā<ā0.001]. The unenhanced model yielded an AUC of 0.916 (95% CI 0.826, 0.968) in the training cohort, 0.860 (95% CI 0.704, 0.953) in the validation cohort and 0.838 (95% CI 0.687, 0.935) in the external validation cohort, showing favorable discriminative efficacy. The triphasic radiomics model was also identified as a good classifier, yielding an AUC of 0.964, 0.880, and 0.843 in three cohorts, respectively.
Delongās test was performed for comparison. No significant differences were found between the unenhanced and triphasic models (pā=ā0.949). The ROC analysis of two radiomics models in testing and external validation cohort is shown in Fig.Ā 4. The AUC, 95% CI of AUC, sensitivity, specificity, accuracy, and a cut-off value are shown in Table 3. The radiomics score of unenhanced radiomics model was calculated by following formula: āRad-scoreā=ā0.1324āĆāoriginal_shape_LeastAxisLengthā+ā0.1302āĆāgradient_glcm_ClusterTendencyā+ā0.1020āĆāwavelet-HLH_firstorder_Kurtosisā+ā0.0950āĆāwavelet-LHH_glrlm_LongRunLowGrayLevelEmphasisā+ā0. 0322āĆāsquare_glcm_Idmnā+ā0.0247āĆāwavelet-LLH_glrlm_RunEntropyā+ā0.0201āĆāwavelet-LHH_firstorder_Rangeā+ā0.0189āĆāwavelet-HLH_firstorder_Rangeāāā0.0219āĆāwavelet-LHH_firstorder_Minimumāāā0.0323āĆāsquareroot_gldm_SmallDependenceEmphasisāāā0.0707āĆāoriginal_firstorder_Skewnessā+ā0.4861ā.
Performance of the unenhanced model and triphasic model. Receiver operating characteristic (ROC) curves of the radiomics signature in the testing and external validation cohorts of the unenhanced model A, B and triphasic model C, D, respectively. AUC, area under the receiver operating characteristic curve
The radiomics nomogram establishment and assessment of the performance of different models
Aldosterone count, maximum diameter, and Rad-score (unenhanced radiomics model) were incorporated to develop a radiomics nomogram in the training cohort (Fig.Ā 5a). The calibration curve of the radiomics nomogram demonstrated good agreement between the predicted and expected probabilities in three cohorts (Fig.Ā 5bād). The p values of HosmerāLemeshow test in training, validation, and external validation cohorts were 0.077, 0.424, and 0.410, respectively.
Radiomics nomogram and calibration curves. A The radiomics nomogram, combining maximum diameter, aldosterone and Rad-score, was developed in the training set. BāD The nomogram calibration curves in training (B), testing (C), and external validation (D) sets. Calibration curves indicate the goodness-of-fit of the model
The diagnostic performance of each model is shown in Table 4. The ROC curves of the clinical model and radiomics nomogram are presented in Fig.Ā 6. The radiomics nomogram showed the highest discrimination in the training cohort, with an AUC of 0.951 (95% CI 0.873, 0.998), and also achieved satisfactory predictive efficacy in the validation cohort [AUC, 0.938 (95% CI 0.805, 0.980)] and in the external validation cohort [AUC, 0.893 (95% CI 0.754, 0.968)]. The observed AUC value was higher than that of the clinical model in external validation [AUC, 0.680 (95% CI 0.514, 0.818); pā=ā0.030]. The nomo-score was acquired using the following formula: āNomo- scoreā=ā27.5663āĆāRad-scoreāāā0.4115āĆāmaximum diameterā+ā0.0094āĆāaldosteroneā+ā2.0275ā.
Performance of the clinical model and the radiomics nomogram. Comparison of receiver operating characteristic (ROC) curves of the clinical model (A, B) and radiomics nomogram (C, D) for the prediction of iMAD and LPA in the testing and external validation cohorts, respectively. AUC, area under the receiver operating characteristic curve
The decision curve analyses for the clinical model, unenhanced radiomics model, and radiomics nomogram in the external validation are presented in Fig.Ā 7. Radiomics nomogram had the highest overall net benefit in a large range of threshold probability.
Decision curve analysis for different models. The y-axis shows the net benefit; the x-axis shows the threshold probability. The decision curves indicate that applying a radiomics nomogram to predict LPA adds great benefit than treating all or none of the patients across the large range of reasonable threshold probabilities
Discussion
Our study established a radiomics nomogram that includes two clinical factors and an unenhanced radiomics signature for preoperative discrimination of iMAD and LPA. The proposed nomogram had better diagnostic performance than the clinical model alone or radiomics signature in differentiating iMAD from LPA. This nomogram can further the treatment decision-making and improve clinical outcomes.
In clinical practice, an isolated adrenal nodule can be incidentally observed in a routine abdominal or chest CT scan. A nodule that has a smaller (generallyāā¤ā4Ā cm) size, higher attenuation (generallyāā„ā20Hu), or no progression during a follow-up visit can be used as an indicator of the benign lesion [25]. However, the differential diagnosis between small benign adrenal lesions can be challenging, considering that patients may present similar symptoms, including hypokalemia, hypertension, and aldosteronism. Due to the lack of fat contents, LPA presents atypical features and can be easily misdiagnosed with iMAD [3, 9]. However, accurate diagnosis is extremely important considering that patients with LPA and iMAD require different treatment approaches [12]. Surgery in patients with adrenal hyperplasia has been associated with certain side effects [16]. The benefit of surgery performed on iMAD patients is still indeterminate.
To our knowledge, only a few studies have investigated the difference between adrenal adenoma and hyperplasia. Leung et al. [26] created a clinical prediction score using 3 parameters (plasma renin activity before saline infusion in saline infusion testāā¤ā0.26Ā ng/mL/h, age at diagnosisā<ā50Ā years, and aldosterone level after saline infusion in saline infusion testāā„ā424Ā pmol/L) to predict aldosteroneāproducing adenoma from idiopathic adrenal hyperplasia. Other clinical prediction criterion also emphasized the effect of aldosterone [27]. This is consistent with our result that patients with LPA had higher aldosterone levels and aldosterone was an independent predictor. Park et al. analyzed several CT parameters, pointing out significant differences in lesion size between adrenal adenoma and adrenal hyperplasia [12]. In the present study, we found that maximum diameter was an independent predictive factor and LPA is larger than iMAD. CT attenuation of the venous phase can serve as a potential risk predictor. However, multiple logistic regression analysis suggested that CT attenuation of the venous phase is not an independent risk predictor. This is proved by other researches. Triphasic CT attenuation and percent washout could not be used for differentiation [12]. And an MRI-based study demonstrated that the quantitative evaluation of washout parameters is not sensitive enough to differentiate LPA from other non-adenomas [9]. However, Seo et al. [1] reported differences in absolute percent washout and relative percent washout between adenoma and non-adenoma, which contradicts aforementioned and our researches. This indicates that heterogeneity may exist based on traditional image factors.
Radiomics has been widely applied in the field of adrenal tumors. For example, texture analysis was used for differentiating malignant from benign adrenal tumors [21]. Furthermore, Chen et al. [28] and He et al. [29] used radiomics signature and nomogram to confirm the risk of functional adenoma. Moreover, radiomics can be an effective medical approach for predicting subclinical pheochromocytoma from lipid-poor adenoma [22]. In this study, radiomics was used to differentiate iMAD from LPA. Diffuse adrenal hyperplasia and adenoma with fat content can easily be identified by most radiologists, thus were excluded in this study. A nomogram was established and yielded a higher AUC than the clinical model did, showing a favorable clinical application value.
In this study, radiomics features concerning shape, grey level, and texture heterogeneity conduced to the radiomics signature. āoriginal_shape_LeastAxisLengthā indicates that tumor size is an important factor. This is consistent with our discovery that size is of great importance. It was also used in another research explaining difference of lesionsā shape [30]. āwavelet-HLH_firstorder_Kurtosisā is histogram-based and represented the position of peak height that indicates CT attenuation value of the maximum number of voxels [31]. A more homogeneous pixel distribution reflected a more regular nodular architecture. In addition, āoriginal_firstorder_Skewnessā represented pixel asymmetry, indicating different heterogeneity of two kinds of nodules. This was also proved in Caruso et al.ās research as malignant tumor had higher heterogeneity [32]. Briefly, features of first-order and texture (GLCM, GLSZM, etc.) features reflected tumor heterogeneity and microenvironment, which is consistent with other studies [33,34,35].
We developed two radiomics signatures: an unenhanced radiomics model deriving from unenhanced CT features and a triphasic radiomics model deriving from unenhanced and enhanced CT features. Both models showed satisfying differential ability. Delongās test indicated no significant difference between the two signatures, which implied that enhanced CT images might fail to provide more necessary texture information compared to unenhanced CT. Whether enhanced radiomics signatures are useful seems to be dubious and needs reconsidering. Yi and colleagues [22] observed the similar diagnostic efficacy of the unenhanced model and enhanced model in differentiating subclinical pheochromocytoma from adrenal lipid-poor adenoma. Sui et al. [36] found that the unenhanced CT model showed higher accuracy, sensitivity, and specificity in distinguishing high-risk anterior mediastinal lesions. He and his team [37] obtained identical conclusions when examining lung diseases. The above-mentioned studies seem to be supportive of our conclusion. The possible interpretation was that contrast enhancement was related to angiogenic properties of tumors [38] and both two diseases are rich in blood supply. So, the minor difference of angiogenic property under the microscope between iMAD and LPA was not powerful enough to provide more favorable details in regard to the difference in vascularity for the enhanced model. Moreover, the biological heterogeneity within the tumor depicted by radiomics features may be confounded by the intravenously injected contrast material, which may have trivial benefits in discrimination between iMAD and LPA. In brief, the cost-effectiveness of the enhanced-CT should be accurately assessed and the potential risks of an additional CT scan, such as radiation hazards and the potential risks associated with contrast media (allergy, potential renal damage and so on), should be considered. The risks related to contrast agent have been reported to be exacerbated in the pediatric and elderly populations.
This study has a few limitations. First, this was a retrospective study with small sample size; training and validation groups were also derived from one institute. Second, owing to the inherent deficiency of retrospective research, the CT images were from different scanners and parameters, which may influence the texture analysis. Also, resampling was arranged to minimize the effect as much as possible. Finally, 3-D manual ROI segmentation can be time-consuming and complicated to perform. Accordingly, an automatic segmentation method for adrenal lesions with favorable reliability and reproducibility needs to be developed and applied.
Conclusions
In this study, we developed a CT-based radiomics nomogram that has favorable predictive effectiveness for preoperative differentiation of iMAD from LPA. Unenhanced CT protocol is of favorable clinical use to avoid external radiation and potential contrast agent risk. As a non-invasive and quantitative method, radiomics nomogram may serve as an effective tool to supplement the conventional imaging modalities for the clinical decision-making process.
Availability of data and materials
The datasets used and analysed during the current study available from the corresponding author on reasonable request.
Abbreviations
- iMAD:
-
Isolated micronodular adrenal hyperplasia
- LPA:
-
Lipid-poor adenoma
- LASSO:
-
Least absolute shrinkage and selection operator
- ROC:
-
Receiver operator characteristic curve
- DCA:
-
Decision curve analysis
- AUC:
-
Areas under receiver operating characteristic curves
- OR:
-
Odds ratios
- CI:
-
Confidence intervals
- ICC:
-
Inter- and intra- class correlation coefficients
References
Seo JM, Park BK, Park SY, et al. Characterization of lipid-poor adrenal adenoma: chemical-shift MRI and washout CT. AJR Am J Roentgenol. 2014;202(5):1043ā50.
Johnson P, Horton K, Fishman E. Adrenal mass imaging with multidetector CT: pathologic conditions, pearls, and pitfalls. Radiogr Rev Publ Radiol Soc N Am. 2009;29(5):1333ā51.
Elbanan MG, Javadi S, Ganeshan D, et al. Adrenal cortical adenoma: current update, imaging features, atypical findings, and mimics. Abdom Radiol (NY). 2020;45(4):905ā16.
Song JH, Chaudhry FS, Mayo-Smith WW. The incidental adrenal mass on CT: prevalence of adrenal disease in 1,049 consecutive adrenal masses in patients with no known malignancy. AJR Am J Roentgenol. 2008;190(5):1163ā8.
Chaudhary V, Bano S. Anatomical and functional imaging in endocrine hypertension. Indian J Endocrinol Metab. 2012;16(5):713ā21.
Panda A, Das CJ, Dhamija E, et al. Adrenal imaging (part 1): imaging techniques and primary cortical lesions. Indian J Endocrinol Metab. 2015;19(1):8ā15.
Bourdeau I, Parisien-La Salle S, Lacroix A. Adrenocortical hyperplasia: a multifaceted disease. Best Pract Res Clin Endocrinol Metab. 2020;34(3):101386.
Chevalier B, Vantyghem MC, Espiard S. Bilateral adrenal hyperplasia: pathogenesis and treatment. Biomedicines. 2021;9(10):1397.
Romeo V, Maurea S, Guarino S, et al. The role of dynamic post-contrast T1-w MRI sequence to characterize lipid-rich and lipid-poor adrenal adenomas in comparison to non-adenoma lesions: preliminary results. Abdom Radiol (NY). 2018;43(8):2119ā29.
Maghrabi A, Yaqub A, Denning KL, et al. Challenges in the diagnostic work-up and management of patients with subclinical Cushingās syndrome and bilateral adrenal masses. Endocr Pract. 2013;19(3):515ā21.
Reincke M, Bancos I, Mulatero P, et al. Diagnosis and treatment of primary aldosteronism. Lancet Diabetes Endocrinol. 2021;9(12):876ā92.
Park SY, Park BK, Park JJ, et al. Differentiation of adrenal hyperplasia from adenoma by use of CT densitometry and percentage washout. AJR Am J Roentgenol. 2016;206(1):106ā12.
Amar L, Plouin P, Steichen O. Aldosterone-producing adenoma and other surgically correctable forms of primary aldosteronism. Orphanet J Rare Dis. 2010;5:9.
Komada H, Hirota Y, So A, et al. Insulin secretion and sensitivity before and after surgical treatment for aldosterone-producing adenoma. Diabetes Metab. 2020;46(3):236ā42.
Yang Y, Reincke M, Williams TA. Prevalence, diagnosis and outcomes of treatment for primary aldosteronism. Best Pract Res Clin Endocrinol Metab. 2020;34(2):101365.
Meloche-Dumas L, Mercier F, Lacroix A. Role of unilateral adrenalectomy in bilateral adrenal hyperplasias with Cushingās syndrome. Best Pract Res Clin Endocrinol Metab. 2021;35(2):101486.
Kumorowicz-Czoch M, Dolezal-Oltarzewska K, Roztoczynska D, et al. Causes and consequences of abandoning one-stage bilateral adrenalectomy recommended in primary pigmented nodular adrenocortical diseaseācase presentation. J Pediatr Endocrinol Metab. 2011;24(7ā8):565ā7.
Ma D, Liu X, Zeng L, et al. The role of adrenal venous sampling and computed tomography in the management of primary aldosteronism. J Hypertens. 2021;39(2):310ā7.
Gillies R, Kinahan P, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278(2):563ā77.
Lambin P, Leijenaar RTH, Deist TM, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. 2017;14(12):749ā62.
Yu H, Parakh A, Blake M, et al. Texture analysis as a radiomic marker for differentiating benign from malignant adrenal tumors. J Comput Assist Tomogr. 2020;44(5):766ā71.
Yi X, Guan X, Zhang Y, et al. Radiomics improves efficiency for differentiating subclinical pheochromocytoma from lipid-poor adenoma: a predictive, preventive and personalized medical approach in adrenal incidentalomas. EPMA J. 2018;9(4):421ā9.
Kong J, Zheng J, Wu J, et al. Development of a radiomics model to diagnose pheochromocytoma preoperatively: a multicenter study with prospective validation. J Transl Med. 2022;20(1):31.
Benchoufi M, Matzner-Lober E, Molinari N, et al. Interobserver agreement issues in radiology. Diagn Interv Imaging. 2020;101(10):639ā41.
Torresan F, Crimi F, Ceccato F, et al. Radiomics: a new tool to differentiate adrenocortical adenoma from carcinoma. BJS Open. 2021;5(1):zraa061.
Leung HT, Woo YC, Fong CHY, et al. A clinical prediction score using age at diagnosis and saline infusion test parameters can predict aldosterone-producing adenoma from idiopathic adrenal hyperplasia. J Endocrinol Invest. 2020;43(3):347ā55.
Jl P, Mm W, Jc P, et al. Predictive value of preoperative tests in discriminating bilateral adrenal hyperplasia from an aldosterone-producing adrenal adenoma. J Clin Endocrinol Metab. 2000;85(12):4526ā33.
Chen PT, Chang D, Liu KL, et al. Radiomics utilization to differentiate nonfunctional adenoma in essential hypertension and functional adenoma in primary aldosteronism. Sci Rep. 2022;12(1):8892.
He K, Zhang ZT, Wang ZH, et al. A clinical-radiomic nomogram based on unenhanced computed tomography for predicting the risk of aldosterone-producing adenoma. Front Oncol. 2021;11:634879.
Martens RM, Koopman T, Noij DP, et al. Predictive value of quantitative (18)F-FDG-PET radiomics analysis in patients with head and neck squamous cell carcinoma. EJNMMI Res. 2020;10(1):102.
Chong Y, Kim JH, Lee HY, et al. Quantitative CT variables enabling response prediction in neoadjuvant therapy with EGFR-TKIs: Are they different from those in neoadjuvant concurrent chemoradiotherapy? PLoS ONE. 2014;9(2):e88598.
Caruso D, Zerunian M, Daffina J, et al. Radiomics and functional imaging in lung cancer: the importance of radiological heterogeneity beyond FDG PET/CT and lung biopsy. Eur J Radiol. 2021;142:109874.
Sun K, Shi L, Qiu J, et al. Multi-phase contrast-enhanced magnetic resonance image-based radiomics-combined machine learning reveals microscopic ultra-early hepatocellular carcinoma lesions. Eur J Nucl Med Mol Imaging. 2022;49(8):2917ā28.
Wang W, Gu D, Wei J, et al. A radiomics-based biomarker for cytokeratin 19 status of hepatocellular carcinoma with gadoxetic acid-enhanced MRI. Eur Radiol. 2020;30(5):3004ā14.
Wu J, Liu A, Cui J, et al. Radiomics-based classification of hepatocellular carcinoma and hepatic haemangioma on precontrast magnetic resonance images. BMC Med Imaging. 2019;19(1):23.
Sui H, Liu L, Li X, et al. CT-based radiomics features analysis for predicting the risk of anterior mediastinal lesions. J Thorac Dis. 2019;11(5):1809ā18.
He L, Huang Y, Ma Z, et al. Effects of contrast-enhancement, reconstruction slice thickness and convolution kernel on the diagnostic performance of radiomics signature in solitary pulmonary nodule. Sci Rep. 2016;6:34921.
Tomaszewski MR, Gillies RJ. The biological meaning of radiomic features. Radiology. 2021;298(3):505ā16.
Acknowledgements
Not applicable.
Funding
This study has received funding by a grant from the National Natural Science Foundation of China Grant (81871354 and 81571672), Academic Promotion Program of Shandong First Medical University (2019QL023).
Author information
Authors and Affiliations
Contributions
HT-Y and CS-J contributed to the data analysis and manuscript writing. BK, KS and SN-Q were involved in performing data collection and quality control. XM-W was responsible for the study coordination. All authors have read and approved the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The ethics committee of the Shandong Provincial Hospital Affiliated to Shandong First Medical University approved this retrospective study. All methods were carried out in accordance with relevant guidelines and regulations. Informed consent (written consent) was obtained from all individual participants included in the study.
Consent for publication
Patients signed informed consent regarding publishing their clinical data and images.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yuan, H., Kang, B., Sun, K. et al. CT-based radiomics nomogram for differentiation of adrenal hyperplasia from lipid-poor adenoma: an exploratory study. BMC Med Imaging 23, 4 (2023). https://doi.org/10.1186/s12880-022-00951-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12880-022-00951-x