Multimodal neuroimaging features provide opportunities for accurate classification and personalized treatment options in the psychiatric domain. This study aimed to investigate whether brain features predict responses to the overall treatment of schizophrenia at the end of the first or a single hospitalization. Structural and functional magnetic resonance imaging (MRI) data from two independent samples (N = 85 and 63, separately) of schizophrenia patients at baseline were included. After treatment, patients were classified as responders and non-responders. Radiomics features of gray matter morphology and functional connectivity were extracted using Least Absolute Shrinkage and Selection Operator. Support vector machine was used to explore the predictive performance. Prediction models were based on structural features (cortical thickness, surface area, gray matter regional volume, mean curvature, metric distortion, and sulcal depth), functional features (functional connectivity), and combined features. There were 12 features after dimensionality reduction. The structural features involved the right precuneus, cuneus, and inferior parietal lobule. The functional features predominately included inter-hemispheric connectivity. We observed a prediction accuracy of 80.38% (sensitivity: 87.28%; specificity 82.47%) for the model using functional features, and 69.68% (sensitivity: 83.96%; specificity: 72.41%) for the one using structural features. Our model combining both structural and functional features achieved a higher accuracy of 85.03%, with 92.04% responder and 80.23% non-responders to the overall treatment to be correctly predicted. These results highlight the power of structural and functional MRI-derived radiomics features to predict early response to treatment in schizophrenia. Prediction models of the very early treatment response in schizophrenia could augment effective therapeutic strategies. 1962 | CUI et al.