Skip to main content

A radiomics based approach using adrenal gland and periadrenal fat CT images to allocate COVID-19 health care resources fairly

Abstract

Background

The value of radiomics features from the adrenal gland and periadrenal fat CT images for predicting disease progression in patients with COVID-19 has not been studied extensively. We assess the value of radiomics features from the adrenal gland and periadrenal fat CT images in predicting COVID-19 disease exacerbation.

Methods

A total of 1,245 patients (685 moderate and 560 severe patients) were enrolled in a retrospective study. We proposed a 3D V-net to segment adrenal glands in onset CT images automatically, and periadrenal fat was obtained using inflation operation around the adrenal gland. Next, we built a clinical model (CM), three radiomics models (adrenal gland model [AM], periadrenal fat model [PM], and fusion of adrenal gland and periadrenal fat model [FM]), and radiomics nomogram (RN) after radiomics features extracted.

Results

The auto-segmentation framework yielded a dice value 0.79 in the training set. CM, AM, PM, FM, and RN obtained AUCs of 0.717, 0.716, 0.736, 0.760, and 0.833 in the validation set. FM and RN had better predictive efficacy than CM (P < 0.0001) in the training set. RN showed that there was no significant difference in the validation set (mean absolute error [MAE] = 0.04) and test set (MAE = 0.075) between predictive and actual results. Decision curve analysis showed that if the threshold probability was between 0.4 and 0.8 in the validation set or between 0.3 and 0.7 in the test set, it could gain more net benefits using RN than FM and CM.

Conclusions

Radiomics features extracted from the adrenal gland and periadrenal fat CT images are related to disease exacerbation in patients with COVID-19.

Peer Review reports

Background

The coronavirus disease 2019 (COVID-19) has caused serious public health problems. Although a diagnosis and treatment protocol suitable for the national situation has been enacted in each country, the consensus is to separate moderate and severe patients who need special care. Predicting early disease progression and accurately identifying which patients are at risk of developing severe or critical disease is not only conducive to early disease control. However, it will also effectively allocate healthcare resources globally.

Because it was unclear that the polymerase chain reaction (PCR) technique [1], clinical symptoms, and laboratory tests were correlated with progression in patients with COVID-19 [2], researchers started exploring predicting the disease prognosis of COVID-19 combining artificial intelligence (AI) technology and CT images. Zhang et al. [3] constructed a comprehensive system and proved that the combination of chest CT and AI technology is a potential tool to predict disease progression in patients with COVID-19. Most studies have focused on the relationship between features from chest CT images and disease progression and ignored the endocrine system, which plays an essential role in disease progression [4,5,6]. However, studies showing whether the prognosis of patients with COVID-19 can be affected by direct damage or indirect change of adrenal glands or periadrenal fat are lacking.

The entire system and the organs are closely related to the immune-inflammatory state (78). Moreover, the adrenal gland can be changed not only by receiving a direct attack from SARS-CoV-2 through angiotensin-converting enzyme-2 (ACE2) [9] but also by adapting to physiological needs and pathological conditions because it has an astonishing regenerative capacity to adjust to unique micro-environments [10]. Therefore, we hypothesized that adrenal gland changes could occur due to its adaptive mechanism and SARS-CoV-2 invasion. The periadrenal fat was bound to be affected as the closest tissue around the adrenal gland. Fundamental changes can be detected by extracting radiomics features from chest CT images using AI technology [11]. Manual segmentation is the first and critical step in successfully constructing a radiomics model, but it is time-consuming and error-prone because of repeated labor. Combining the auto-segmentation framework based on AI technology with manual revision improvesefficiency and quality (1213).

We, therefore, constructed an auto-segmentation framework and developed a radiomics nomogram (RN) to predict progression in patients with COVID-19 by integrating radiomics features from adrenal glands and periadrenal fat CT images with clinical indicators and to determine that radiomics features extracted from the adrenal gland and periadrenal fat CT images were related to possibility of deterioration in patients with COVID-19.

Methods

Data resources and image grouping

CT images and clinical data were retrospectively and consecutively collected between January 1 and April 30, 2020, from two hospitals: Huoshenshan Hospital (HSH) (n = 1,209) and Maternal and Child Health Hospital Optical Valley Branch Hospital of Hubei Province (MCH) (n = 36), China. Diagnosis and clinical classification of patients with COVID-19 were confirmed according to Diagnosis and Treatment Protocol for Novel Coronavirus Pneumonia (Trial Version 8). 30 clinical indicators were collected including age, sex, lnterleukin 6 (IL-6) levels, white blood cell count (WBC), lymphocyte count (L), neutrophil count (N), hemoglobin (HB ), red blood cell (RBC), blood platelet count (PLT), C-reactive protein (CRP), procalcitonin (PCT), prothrombin time (PT), thrombin time (TT), fibrinogen (FIB), D dimer (DD), blood sugar (BS), alanine transaminase (ALT), glutamic oxaloacetic transaminase (AST), total protein (TP), albumin (ALB), total bilirubin (TBil), direct bilirubin (D-Bil), blood urea nitrogen (BUN), creatinine (Cre), B-type natriuretic peptide (BNP), myoglobin (Mb), troponin (Tn), lactate dehydrogenase (LDH), creatine kinase (CK) and creatine kinase-MB (CK-MB). Moreover, the missing value, more than 20% of the total data, were excluded. If there were incomplete clinical or blood examination data, we used the median of the variable distribution to fix it. Patients of age ≤ 14 years (n = 2), low quality or unreadable chest CT images, and lesion or changes in adrenal glands CT images ( verified by two experienced radiologists) (n = 209) were excluded. All images were non-enhanced chest CT images and reconstructed at a slice thickness of 1.00 mm. Details of CT characteristics are listed in Supplementary Table 1. Patients with adrenal lesions were excluded after being evaluated by two radiologists with more than ten years of experience. We chose chest CT images scanned within four days of the first diagnosis as the onset image. If the CT scan was done more than once, we chose the one closer to the admission date. Supplementary Fig. 1 demonstrates the inclusion and exclusion criteria. Figure 1 shows the workflow of our study.

Fig. 1
figure 1

The workflow of the study

CT images of 1,209 patients from Huoshenshan Hospital were used to construct the prediction models divided into the training set (967/1209, 80%) and the validation set (242/1209, 20%) randomly. We constructed five models: the adrenal gland radiomics model (AM), periadrenal fat radiomics model (PM), a fusion of the adrenal gland and periadrenal fat radiomics model (FM), clinical model (CM) and radiomics nomogram (RN). A total of 36 Patients (13 moderate and 23 severe patients from Maternal and Child Health Hospital Optical Valley Branch Hospital) were included as an independent test set.

Auto-segmentation framework and ROIs delineation

We proposed a cascaded V-net network framework to automatically segment the bilateral adrenal glands. The cascaded V-net network framework was constructed using two V-nets based on the coarse-to-fine principle, including a coarse segmentation network for rapidly locating the target area, namely the adrenal glands, and an exemplary segmentation network for optimization and precisely delineating the adrenal glands. Coarse segmentation adopts 3 × 3 × 3 resolution to quickly find the target region (global sampling) from the global, and fine segmentation adopts 1 × 1 × 1 to carry out more detailed segmentation (mask sampling) in the target region. Dice loss is used as a loss function. ROI of periadrenal fat was obtained using the inflation algorithm based on the adrenal gland. An experienced radiologist (Y.F.) manually delineated the ROIs of bilateral adrenal glands, which werethe ground truth labels of adrenal, and he was asked to delineate the adrenals according to the image delineating principles of BraTS 2018. They compared the automated segmentation with the ground truth; any discrepancies or errors identified during the review are manually corrected or adjusted (1415). We trained a coarse localization model as Model 1 for the auto-segmentation framework, which can perform coarse segmentation to locate the adrenal gland area. The second V-net model was used to define segmentation as Model 2 and further divided into the left and right adrenal glands. The segmentation frame is shown in Fig. 2. ROIs, including bilateral adrenal glands and periadrenal fat from all CT images without annotation, were segmented using the auto-segmentation framework. Dice value was used to evaluate the cascaded V-net network framework.

Fig. 2
figure 2

The framework of the used V-net

Radiomics feature selection and models building

Radiomic features were extracted from ROIs on the CT images using a Python package (PyRadiomics V3.0). B-spline interpolation resampling was used to normalize the voxel size, and the anisotropic voxels were resampled to form isotropic voxels of 1.0 mm × 1.0 mm × 1.0 mm in the feature extraction. A total of 2,264 radiomics features of each ROI were extracted. Next, the feature in the training set was preprocessed by standardization. The mean and variance in the training set were applied to the validation and test sets.

First, a univariate analysis named K-best was employed to reduce dimensionality and select significant features. It selected features according to K’s highest scores computed through ANOVA F-value and P-value between label and features. Features with a significant difference (P < 0.05) were selected. Next, the LASSO feature-selection algorithm extracted the most informative radiomics features to prevent the “curse of dimensionality.“ The LASSO algorithm assigns a penalty (lambda) to each feature based on its contribution to the model’s performance. The optimal regularization parameter (lambda) for the LASSO algorithm is typically determined using 10-fold cross-validation. By varying the value of lambda, only the most essential features that contribute significantly to predicting the outcome are included in the final radiomics model. After feature extraction and selection, logistic regression (LR) algorithms were trained to construct three radiomics models (AM, PM, and FM) for predicting the disease progression of COVID-19. We used five-fold cross-validation to find the optimal hyperparameters of the logistic regression model, then used the optimal hyperparameters to construct models. Lbfgs solver was used in the optimization problem. In the FM model, the RadScore of the patient was calculated according to the LASSO algorithm.

RN construction and evaluation

We used univariate analysis to assess the relationship between clinical factors plus serum biomarkers and disease outcome. Features with P < 0.05 were introduced into multivariate LR analysis, and the optimal radiomics feature subset was selected by LASSO five-fold cross-validation. The two make the best combination. The best training and validation sets assignment was chosen for the subsequent analysis. Next, we applied the multivariate LR model to build CM using valuable clinical indicators and RN using the RadScore from FM with clinical indicators to predict the disease progression of COVID-19.

We diagnosed collinearity by calculating the VIF for variables in RN to detect multicollinearity among the radiomics nomogram variables. In the end, RN was verified in the validation and test sets. Calibration curves and the Hosmer–Lemeshow test assessed the relation between the predicted risks and actual results. DCA was used to evaluate the performance of the RN.

Statistics

Before model building, differences in clinical factors and serum biomarkers between moderate and severe patient sets were assessed using the Mann-Whitney U test or Student’s t-test for continuous variables and the χ2 test or Fisher’s exact test for categorical variables. We analyzed all data using SPSS for Windows version 26.0 (IBM Corp., Armonk, New York, USA). P < 0.05 was considered a statistically significant difference.

Dice value was used to assess the effectiveness of the auto-segmentation framework. The AUC of receiver-operating characteristics (ROC) with 95% confidence interval (95% CI), sensitivity, and specificity were used to evaluate the performance of AM, PM, FM, CM, and RN. Accuracy was calculated to assess the prediction performance. Differences in AUC values among different models were estimated using the DeLong test. Calculate three statistics for each pair of models being compared: AUC - ROC, variance of AUC - ROC, and covariance of AUC - ROC. The DeLong test whether there is a significant difference between the AUC - ROC values of the two models based on their variances and covariance.

Results

Patient characteristics

The data of 1,209 patients were formed as our training and validation sets, including 672 patients from moderate COVID-19 and 537 from severe COVID-19. The patient characteristics in training and validation sets are listed in Table 1. No significant differences were observed between the training and validation set in sex (P = 0.238). IL-6, WBC, L, N, CRP, PCT, PT, DD, BS, AST, ALB, D-Bil, BNP, LDH, and CK-MB differed significantly between moderate and severe pneumonia sets both in training and validation sets (P < 0.05).

Table 1 Clinical characteristics of patients in training and validation set (n = 1209)

Adrenal gland and periadrenal fat auto-segmentation framework

For adrenal gland segmentation, we manually delineated bilateral adrenal glands from the CT images of 315 patients; 265 were used for training. The remaining data from 50 patients were used to evaluate the performance. The segmentation model yielded average Dice values of 79.48% for the left and 78.55% for the right adrenal gland. The entire adrenal gland achieved an average Dice value of 79.02%. Representative auto-segmentation results are shown in Fig. 3. The segmentation algorithm was then used to segment all the remaining data automatically.

Fig. 3
figure 3

The comparison of auto-segmentation results and ground truth on representative cases from the training set. The upper line shows the ground truth, and the lower line shows the auto-segmentation results. The green color shows the right adrenal glands, whereas red denotes the left. Yellow shows right periadrenal fat, and purple shows left

Radiomics feature and clinical indicator selection

In the training set, the number of radiomics features was reduced to 23 for building AM that included 8 first-order features and 15 texture features (Gray Level Co-occurrence Matrix [GLCM] = 3, Gray Level Size Zone Matrix [GLSZM] = 8, Gray Level Run Length Matrix [GLRLM] = 2 and Gray Level Dependence Matrix [GLDM] = 2); 68 for PM that included 11 first-order features, 2 sharp feature, and 55 texture features (GLCM = 13, GLSZM = 25, GLRLM = 7, GLDM = 9 and Neighboring Gray Tone Difference Matrix [NGTDM] = 1) and 82 for FM that included 12 first-order features and 70 texture features (GLCM = 17, GLSZM = 38, GLRLM = 3, GLDM = 5 and NGTDM = 7). These features were evaluated to construct three radiomics models.

A total of 30 clinical factors and serum biomarkers were analyzed in our study. Next, univariate logistic regression analysis selected 17 clinical factors and serum biomarkers. They are LDH, CRP, age, ALB, L, N, Hb, RBC, DD, BS, WBC, CK_MB, TT, BUN, AST, TP, PT. and 7 indicators, LDH, L, HB, DD, WBC, TT, and TP, were selected using multivariate logistic regression analysis. The relationship between RadScore from FM used in constructing a radiomics nomogram (RN) and 30 clinical factors plus serum biomarkers were analyzed using Pearson correlation between training, validation, and test sets (Fig. 4). The difference in RadScores with clinical factors or serum biomarkers was not significant. Then, 17 clinical factors and serum biomarkers were selected using univariate logistic regression analysis, and 7 indicators, LDH, L, HB, DD, WBC, TT, and TP, were selected using multivariate logistic regression analysis.

Fig. 4
figure 4

Correlation heatmap in training, validation, and test set. FM = fusion of adrenal gland and periadrenal fat model; CM = clinical model; RN = radiomics nomogram

Three radiomics models and clinical model building

We developed three radiomics models (AM, PM, FM) based on radiomics features and a clinical model (CM) based on the seven selected independent predictive clinical indicators. We used three evaluation indicators (area under the curve [AUC], 95% CI, sensitivity [SEN], and specificity [SPE]) to assess AM, PM, FM, and CM for predicting the progression of patients with COVID-19 in training, validation, and test sets. In general, AM achieved an AUC of 0.692, 0.714, and 0.659 in the training set, validation set, and test set, respectively; PM achieved an AUC of 0.764, 0.736, and 0.645; FM achieved an AUC of 0.791, 0.760 and 0.686; CM obtained an AUC of 0.712,0.717 and 0.692 (Fig. 5, supplementary Table 2).

Fig. 5
figure 5

ROC curves in training, validation, and test sets. AM = adrenal gland model; PM = periadrenal fat model; FM = fusion of adrenal gland and periadrenal fat model; CM = clinical model; RN = radiomics nomogram

Box plots summarizing the RadScores and coefficients of seven clinical indicators in training, validation, and test sets directly demonstrate the difference between RadScore and coefficients of seven clinical indicators between the moderate and severe patient sets (Fig. 6).

Fig. 6
figure 6

Box plots summarizing rad score and 7 clinical indicators between moderate and severe patient sets in training, validation, and test sets. In each box plot, the horizontal line crossing the box is the median, and the bottom and top are the lower and upper quartiles. FM = fusion of adrenal gland and periadrenal fat model; CM = clinical model; RN = radiomics nomogram; LDH = lactate dehydrogenase; L = lymphocyte count; Hb = hemoglobin, DD = D dimer; WBC = white blood cells count; TT = thrombin time and TP = total protein

RN construction and validation

Multivariate analysis revealed that RadScore and seven clinical indicators were significant independent factors predicting disease progression in patients with COVID-19. We conducted collinearity diagnosis by calculating the VIF for variables in RN to detect multicollinearity among the radiomics nomogram variables, and the threshold was set to 10 [16,17,18]. Finally, the VIF value for the radiomics score and seven clinical indicators in RN ranged from 1.007 to 1.191, indicating no severe collinearity in these factors. Next, we used the RadScore from FM combined with seven clinical indicators to construct the RN to assess disease progression in patients with COVID-19 (Fig. 7). The RN showed satisfactory performance in predicting and assessing progression in patients with COVID-19 with an AUC of 0.806 (95% CI, 0.780 to 0.831) in the training set, 0.833 (95% CI, 0.780 to 0.878) in the validation set, and 0.773 (95% CI, 0.603 to 0.895) in the test set (Fig. 5, supplementary Table 2).

Fig. 7
figure 7

Radiomics nomogram developed in training set with radiomics features, lactate dehydrogenase (LDH), lymphocyte count (L), hemoglobin (Hb), D dimer (DD), white blood cells count (WBC), thrombin time (TT) and total protein(TP). Points are assigned for each variable by drawing a line upward from the corresponding variable to the Points line. The sum of points plotted on the total Points line corresponds with the severity of patients with COVID-19.

DeLong’s test was used to compare the AUCs of the training set’s three radiomics models, CM and RN. The result showed that the RN and FM were significantly better than CM (P < 0.0001). The difference between FM and RN was not statistically significant (P = 0.233) in the validation and test sets.

The Hosmer–Lemeshow test was not significant in the validation set (mean absolute error [MAE] = 0.075) or test set (MAE = 0.04), which suggests that there was no significant departure from actual values (Fig. 8). Decision curve analysis (DCA) (Fig. 9) showed that if the threshold probability was between 0.4 and 0.8 in the validation set, the RN could get more net benefits than FM and CM. If the threshold probability was between 0.3 and 0.7 in the test set, RN can still get more net benefits than FM and CM. The threshold cannot be set above 0.8; otherwise, the net benefit would become negative values.

Fig. 8
figure 8

Calibration curves of the radiomics nomogram for predicting the disease progression of COVID-19 pneumonia in validation set (A) and test set (B). The y-axis represents the probability of COVID-19 pneumonia becoming severe, and the x-axis represents the predicted risk. The dashed line was the reference line where an ideal nomogram would lie. The dotted line was the performance of the radiomics nomogram, while the solid line corrected for any bias in the radiomics nomogram. MAE is 0.04 in the validation set and 0.075 in the test set

Fig. 9
figure 9

Decision curve analysis (DCA) for the radiomic model, clinical model, and radiomics nomogram. The y-axis measures the net benefit. Using the clinical model, radiomic model, and radiomics nomogram in the study to predict COVID-19 pneumonia progress adds more benefit than the treat all patients as severity patients scheme or the treat none scheme. The net benefit of the radiomics nomogram was better than clinical and radiomic models in both sets, with several overlaps in the training set

Discussion

In our study, the auto-segmentation framework could robustly localize the adrenal glands and accurately refine their boundary. The RN using adrenal glands and periadrenal fat onset CT images performed well, reflecting that the microscopic changes in adrenal glands and periadrenal fat in patients with COVID-19 can be detected using radiomics features. We hope this model will assist radiologists and clinicians in early intervention, monitoring, and management, ultimately leading to improved clinical outcomes. It can be used in the COVID-19 epidemic, especially when there is a shortage of healthcare workers.

Our results suggested that adrenal gland and periadrenal fat changes on the onset of CT images in severe patients differed from those in moderate patients. Results from autopsies in 10 patients from COVID-19 performed by Zinserling et al. [19] showed that inflammation—small proliferations of cells with enlarged light nuclei and mononuclear infiltration, CD3 + and CD8 + in different layers of adrenal glands and their surrounding tissue, such as periadrenal fat. These changes may be related to direct damage of adrenal glands and indirect changes stimulated by systemic inflammation caused by SARS-CoV-2 or the immune state. In one respect, the changes were caused by direct damage done by SARS-CoV-2. The SARS virus has been identified in adrenal cells, suggesting a direct local replication-mediated cytopathic effect of the virus in adrenal tissue [20]. Moreover, the virus may cause hemorrhage, necrosis, or thrombosis at the adrenal level [21]. Adrenal lesions due to influenza and other viruses have been described previously [22]. Recent findings have indicated the possibility of venous thromboembolism in patients with COVID-19, which may cause an acute adrenal insufficiency that may be an indicator of making the disease worse [23]. All of these changes caused by the COVID-19 virus can change the physical properties of the adrenal tissue and the heterogeneity between tissues and cells in patients with different conditions. First-order and higher-order radiomics features can reflect the changes in physical properties. We selected 82 features, concluding with 12 first-order features and 70 texture features. In FM, the most relevant to the severity of the disease is Robust - Mean - Absolute - Deviation and there is a negative correlation. This feature reflects the amount of variation or dispersion from the mean value which means the higher the feature is, the less severe the disease is. Moreover, in the AM model, we got the same feature with the highest negative correlation with disease severity. The most relevant positive correlation feature is Zone Entropy which reflects the uncertainty or randomness in the distribution of zone sizes and gray levels. A higher value indicates more heterogeneneity in the texture patterns which means more severity of disease. In PM, The feature with the largest positive correlation is also Zone Entropy. This feature reflects the heterogeneity of adrenal glands, which is the same as that extracted feature by tumor radiomics analysis [24]. Our predicting results showed that AM’s AUC value was slightly better than PM’s, but there was no statistical significance. That may indicate that the inner changes in adrenal glands were not different from the changes in periadrenal fat. Further research is needed to confirm whether changes in periadrenal fat were caused by adrenal gland inflammatory cells’ infiltration or direct damage of SARS-CoV-2 [25].

Previous research also reported that rather than direct damage by SARS-CoV-2 itself, adrenal gland changes were believed to be caused by immune overreaction or cytokine storm-inducing endocrinological pathway impaired through a tight binding mechanism of SARS-CoV-2 and ACE2 [26]. Endocrinological pathway impairment may be related to the following aspects. First, ACE2 is highly expressed in human adrenals, which serves as the entry receptor for SARS-CoV-2 that is bound to be significantly affected. Some work has proposed that the imbalanced ACE/ACE2 axis may mediate lung tissue repair and wound healing pathways [27]. Second, the primary substrate for ACE2 is angiotensin II. ACE2 is a negative regulator of the renin-angiotensin-aldosterone system (RAAS) by converting the active angiotensin and angiotensin II to the inactive angiotensin 1–7 [28]. Finally, the SARS virus contains several permutations of amino acid sequences with homology to the antigenic relevant residues of ACTH [29] that will have potential and significant pathophysiological effects due to molecular mimicry. Potential cross-reactivity of antibodies may cause local infiltration with immunocompetent cells in adrenal glands [19].

In summary, these aspects suggested that adrenal cells and surrounding tissue damage in patients with COVID-19 may be caused by a viral infection and further secondary inflammatory plus autoimmune processes located in the adrenal glands. Due to the pandemic of COVID-19, scholars have done much research on computer modeling of COVID-19. Many diagnostic models are machine learning models based on chest CT for detecting COVID-19 infection and predictive models for predicting the risk of death, progression of severe disease, or the length of hospital stay (30). Among COVID-19 patients, the most common predictors of severe prognosis include age, sex and radiomics features, C-reactive protein, lactate dehydrogenase, and lymphocyte counts. These results are consistent with the univariate logistic regression analysis (31). There was little evaluation of the calibration of the predictions, and many prediction models are mainly based on clinical data (32). Compared with other models, we adopted the adrenal grand CT images as the study’s data, which could be used in multi-organ fusion prediction in the future.

Conclusions

We automatically obtained the ROI from the onset of CT images through the auto-segmentation framework and inflation algorithm. We proposed and constructed an adrenal gland auto-segmentation framework based on chest CT images and AI technology. Until now, we have not found any research using adrenal gland CT parameters as indicators to evaluate the progression in patients with COVID-19, especially in the radiomics field. Our work has some limitations in the data collection process. First, we used the chest CT images as data resources, which would affect the visualization and delineation of adrenal tissue. However, considering clinical practicality and radiation to patients, there is no need to perform another CT scan using professional adrenal glands CT parameters to observe adrenal lesions better because CT examination of patients with COVID-19 is mainly to detect and observe pulmonary lesions. Second, unlike other studies, we selected the entire organ as ROIs rather than the lesion. The periadrenal fat area may not be precise and may contain some other indistinguishable tissue from the human eye.

Data Availability

The data used and/or analyzed during the current study are available from the corresponding author on reasonable request. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. The code used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

LASSO:

Least absolute shrinkage and selection operator

CM:

Clinical model AM:adrenal gland model

PM:

Periadrenal fat model

FM:

Fusion of adrenal gland and periadrenal fat model

RN:

Radiomics nomogram model

ROI:

Region of interest

MAE:

Mean absolute error

PCR:

Polymerase chain reaction

AI:

Artificial intelligence

AUC:

Under the curve

ACTH:

Adrenocorticotropic hormone

SNS:

Sympathetic nervous system

ACE2:

Angiotensin-converting enzyme-2

HSH:

Huoshenshan Hospital

MCH:

Maternal and Child Health Hospital Optical Valley Branch Hospital

TTH:

Taikang Tongji (Wuhan) Hospital

GLCM:

Gray-level co-occurrence matrix

GLRLM:

Gray-level run-length matrix

GLSZM:

Gray-level size zone matrix

GLDM:

Gray-level dependence matrix

NGTDM:

Neighboring Gray Tone Difference Matrix

LR:

Logistic regression

DCA:

Decision curve analysis

IL-6:

Interleukin 6

L:

Lymphocyte count

CRP:

C-reactive protein

PCT:

Procalcitonin

DD:

D dimer

AST:

Glutamic oxaloacetic transaminase

BNP:

B-type natriuretic peptide

LDH:

Lactate dehydrogenase

CK:

Creatine kinase

CK-MB:

Creatine kinase-MB

WBC:

White blood cell count

N:

Neutrophil count

HB:

Hemoglobin

RBC:

Red blood cell count

PLT:

Blood platelet count

PT:

Prothrombin time

TT:

Thrombin time

FIB:

Fibrinogen

BS:

Blood sugar

ALT:

Alanine transaminase

TP:

Total protein

ALB:

Albumin

TBil:

Total bilirubin

D-Bil:

Direct bilirubin

BUN:

Blood urea nitrogen

Cre:

Creatinine

Mb:

Myoglobin

Tn:

Troponin

CK:

Creatine kinase

SEN:

Sensitivity

SPE:

Specificity

VIF:

Variance inflation factor

RAAS:

Renin-angiotensin-aldosterone system

References

  1. Tian M, Long Y, Hong Y, Zhang X, Zha Y. The treatment and follow-up of ‘recurrence’ with discharged COVID ‐ 19 patients: data from Guizhou, China. Environ Microbiol. 2020;22(8):3588–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Giorgia Guglielmi. Fast coronavirus tests: what they can and can’t do. Nature. 2020;585:496–8.

    Article  PubMed  Google Scholar 

  3. Zhang K, Clinically Applicable AI. System for Accurate diagnosis, quantitative measurements and prognosis of COVID-19 Pneumonia using computed tomography. Cell. 2020;181(6):1423–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Ting DSW, Carin L, Dzau V, Wong TY. Digital technology and COVID-19. Nat Med. 2020;26:459–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Agarwal S, Agarwal SK. Endocrine changes in SARS-CoV-2 patients and lessons from SARS-CoV. Postgrad Med J. 2020;96(1137):412–6.

    Article  CAS  PubMed  Google Scholar 

  6. Bellastella G, Maiorino MI, Esposito K. Endocrine Complications of COVID-19: what happens to the thyroid and adrenal glands? J Endocrinol Investig. 2020;43(8):1169–7110.

    Article  CAS  Google Scholar 

  7. Waldemar, Kanczkowski, Mariko Sue SRB. 2. The adrenal gland microenvironment in health. Disease and during Regeneration. 2017;16(3):251–65.

    Google Scholar 

  8. Pal R, Banerjee M. COVID-19 and the endocrine system: exploring the unexplored. J Endocrinol Investig. 2020;43(7):1027–31.

    Article  CAS  Google Scholar 

  9. Sever P, Johnston SL. The renin-angiotensin system and SARS-CoV-2 infection: a role for the ACE2 receptor? J Renin Angiotensin Aldosterone Syst. 2020;21(2):1470320320926911.

  10. Berger I, Werdermann M, Bornstein SR, Steenblock C. The adrenal gland in stress – adaptation on a cellular level. J Steroid Biochem Mol Biol. 2019;190:198–206.

    Article  CAS  PubMed  Google Scholar 

  11. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures. They Are Data Radiology. 2016;278(2):563–77.

    PubMed  Google Scholar 

  12. Han M, Yao G, Zhang W, Mu G, Zhan Y, Zhou X, Gao Y. Segmentation of CT Thoracic Organs by Multi-resolution VB-nets. SegTHOR@ISBI. 2019. Corpus ID: 139100153.

  13. Hua R, Huo Q, Gao Y, Sui H, Zhang B, Sun Y, et al. Segmenting brain tumor using cascaded V-Nets in multimodal MR images. Front Comput Neurosci. 2020;14:9.

  14. Rezaeijo SM, Jafarpoor Nesheli S, Fatan Serj M, Tahmasebi Birgani MJ. Segmentation of the prostate, its zones, anterior fibromuscular stroma, and urethra on the MRIs and multimodality image fusion using U-Net model. Quant Imaging Med Surg. 2022; 4786–804.

  15. Salmanpour MR, Hosseinzadeh M, Rezaeijo SM, Rahmim A. Fusion-based Tensor radiomics using reproducible features: application to survival prediction in Head and Neck cancer. Comput Methods Programs Biomed. 2023;240:107714.

    Article  PubMed  Google Scholar 

  16. Salmanpour MR, Rezaeijo SM, Hosseinzadeh M, Rahmim A. Deep versus handcrafted Tensor Radiomics features: prediction of Survival in Head and Neck Cancer using Machine Learning and Fusion techniques. Diagnostics (Basel). 2023;13(10):1696.

    Article  PubMed  Google Scholar 

  17. Jahangirimehr A, Abdolahi Shahvali E, Rezaeijo SM, Khalighi A, Honarmandpour A, Honarmandpour F, Labibzadeh M, Bahmanyari N, Heydarheydari S. Machine learning approach for automated predicting of COVID-19 severity based on clinical and paraclinical characteristics: serum levels of zinc, calcium, and vitamin D. Clin Nutr ESPEN. 2022;51:404–11.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Rezaeijo SM, Abedi-Firouzjah R, Ghorvei M, Sarnameh S. Screening of COVID-19 based on the extracted radiomics features from chest CT images. J Xray Sci Technol. 2021;29(2):229–43.

    CAS  PubMed  Google Scholar 

  19. Zinserling VA, Semenova NY, Markov AG, et al. Inflammatory cell infiltration of adrenals in COVID-19. Horm Metab Res. 2020;52(9):639–41.

    Article  CAS  PubMed  Google Scholar 

  20. Gu J. Pathology and Pathogenesis of severe Acute Respiratory Syndrome. Am J Pathol. 2007;170(4):1136–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Henry BM, Vikse J, Benoit S, et al. Hyperinflammation and derangement of renin-angiotensin-aldosterone system in COVID-19: a novel hypothesis for clinically suspected hypercoagulopathy and microvascular immunothrombosis. Clin Chim Acta. 2020;507:167–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Zinserling AV, Aksenov OA, Melnikova VF, Zinserling VA. Extrapulmonary lesions in Influenza. Tohoku J Exp Med. 1983;140(3):259–72.

    Article  CAS  PubMed  Google Scholar 

  23. Tseng YH, Yang RC, Lu TS. Two hits to the renin-angiotensin system may play a key role in severe COVID-19. Kaohsiung J Med Sci. 2020;36(6):389–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ardakani AA, Bureau NJ, Ciaccio EJ, Acharya UR. Interpretation of radiomics features-A pictorial review. Comput Methods Programs Biomed. 2022;215:106609.

    Article  PubMed  Google Scholar 

  25. Battisti S, Pedone C, Napoli N, Russo E, Agnoletti V, Nigra SG, et al. Computed tomography highlights increased visceral adiposity Associated with critical Illness in COVID-19. Diabetes Care. 2020;43(10):e129–30.

    Article  CAS  PubMed  Google Scholar 

  26. Datta PK, Liu F, Fischer T, et al. SARS-CoV-2 pandemic and research gaps: understanding SARS-CoV-2 interaction with the ACE2 receptor and implications for therapy. Theranostics. 2020;10(16):7448–64.

    Article  CAS  PubMed  Google Scholar 

  27. Delpino MV, Quarleri J. SARS-CoV-2 Pathogenesis: Imbalance in the renin-angiotensin system favors lung fibrosis. Front Cell Infect Microbiol. 2020;10:340.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Cheng H, Wang Y, Wang GQ. Organ - protective effect of angiotensin ‐ converting enzyme 2 and its effect on the prognosis of COVID ‐ 19. J Med Virol. 2020;92(7):726–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Pal R. COVID-19, hypothalamo – pituitary - adrenal axis and clinical implications. Endocrine. 2020;68(2):251–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to express their heartfelt gratitude to their colleague, Mashaal, who kindly took the time to polish the manuscript about the English writing. The study was supported by the Guiyang Science and Technology Project (ZKHT [2022]-4-1-3);  Guizhou Senior Innovative Talent Project (QKHPTRC-GCC[2022]041-1); Beijing Medical and Health Foundation (YWJKJJHKYJJ-B20261CS) and Chongqing Science and Health Joint Medical Research Project(2021MSXM052).

Funding

The study was supported by the Guiyang Science and Technology Project (ZKHT [2022]-4-1-3);  Guizhou Senior Innovative Talent Project (QKHPTRC-GCC[2022]041-1); Beijing Medical and Health Foundation (YWJKJJHKYJJ-B20261CS) and Chongqing Science and Health Joint Medical Research Project(2021MSXM052).

Author information

Authors and Affiliations

Authors

Contributions

First author, Mudan Zhang: design of the work, conducting experiments, analyzing data and writing the manuscript; Co-first author, Xuntao Yin: acquiring data, design of the work, supervising the work and review the manuscript. Wuchao Li: statistics analysis; Yan Zha: design of the work; Xianchun Zeng: Resources; Xiaoyong Zhang: Technology analysis; Jingjing Cui: statistics analysis; Xue Zhong: statistics analysis.Corresponding Author, Chen Liu: Writing - Review & Editing, Project administration; Verified underlying data. Co-corresponding Author, Rongpin Wang: Writing - Review & Editing, Project administration; Verified underlying data. All authors reviewed the manuscript.

Corresponding authors

Correspondence to Rongpin Wang or Chen Liu.

Ethics declarations

Ethics approval and consent to participate

This multicenter study was approved by the ethics committees of Guizhou Provincial People’s hospitals (2020,NO.01). Because of its retrospective nature, the need to obtain informed consent in advance was waived by Guizhou Provincial People’s Hospital. The study was performed according to the principles of the declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, M., Yin, X., Li, W. et al. A radiomics based approach using adrenal gland and periadrenal fat CT images to allocate COVID-19 health care resources fairly. BMC Med Imaging 23, 181 (2023). https://doi.org/10.1186/s12880-023-01145-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12880-023-01145-9

Keywords