BackgroundDue to individual differences and lack of objective biomarkers, only 30-40% patients with major depressive disorder (MDD) achieve remission after initial antidepressant medication (ADM). We aimed to employ radiomics analysis after ComBat harmonization to predict early improvement to ADM in adolescents with MDD by using brain multiscale structural MRI (sMRI) and identify the radiomics features with high prediction power for selection of selective serotonin reuptake inhibitors (SSRIs) and serotonin norepinephrine reuptake inhibitors (SNRIs).Methods121 MDD patients were recruited for brain sMRI, including three-dimensional T1 weighted imaging (3D-T1WI)and diffusion tensor imaging (DTI). After receiving SSRIs or SNRIs for 2 weeks, the subjects were divided into ADM improvers (SSRIs improvers and SNRIs improvers) and non-improvers according to reduction rate of the Hamilton Depression Rating Scale, 17 item (HAM-D17) score. Then, sMRI data were preprocessed, and conventional imaging indicators and radiomics features of gray matter (GM) based on surface-based morphology (SBM) and voxel-based morphology (VBM) and diffusion properties of white matter (WM) were extracted and harmonized with ComBat harmonization. Two-level reduction strategy with analysis of variance (ANOVA) and recursive feature elimination (RFE) was utilized sequentially to decrease high-dimensional features. Support vector machine with radial basis function kernel (RBF-SVM) was used to integrate multiscale sMRI features to construct models for early improvement prediction. Area under the curve (AUC), accuracy, sensitivity, and specificity based on the leave-one-out cross-validation (LOO-CV) and receiver operating characteristic (ROC) curve analysis were calculated to evaluate the model performance. Permutation tests were used for assessing the generalization rate.ResultsAfter 2-week ADM, 121 patients were divided into 67 ADM improvers (31 SSRIs improvers and 36 SNRIs improvers) and 54 ADM non-improvers. After two-level dimensionality reduction, 8 conventional indicators (2 VBM-based features and 6 diffusion features) and 49 radiomics features (16 VBM-based features and 33 diffusion features) were selected. The overall accuracy of RBF-SVM models based on conventional indicators and radiomics features was 74.80% and 88.19%. The radiomics model achieved the AUC, sensitivity, specificity, and accuracy of 0.889, 91.2%, 80.1% and 85.1%, 0.954, 89.2%, 87.4% and 88.5%, 0.942, 91.9%, 82.5% and 86.8% for predicting ADM improvers, SSRIs improvers and SNRIs improvers, respectively.Pvalue of permutation tests were less than 0.001. The radiomics features predicting ADM improver were mainly located in the hippocampus, medial orbitofrontal gyrus, anterior cingulate gyrus, cerebellum (lobule vii-b), body of corpus callosum, etc. The radiomics features predicting SSRIs improver were primarily distributed in hippocampus, amygdala, inferior temporal gyrus, thalamus, cerebellum (lobule vi), fornix, cerebellar peduncle, etc. The radiomics features predicting SNRIs improver were primarily located in the medial orbitofrontal cortex, anterior cingulate gyrus, ventral striatum, corpus callosum, etc.ConclusionsThese findings suggest the radiomics analysis based on brain multiscale sMRI after ComBat harmonization could effectively predict the early improvement of ADM in adolescent MDD patients with a high accuracy, which was superior to the model based on the conventional indicators. The radiomics features with high prediction power may help for the individual selection of SSRIs and SNRIs.Supplementary InformationThe online version contains supplementary material available at 10.1186/s12888-023-04966-8.