Abstract Treatment resistant (TR) psychosis is considered to be a significant cause of disability and functional impairment. Numerous efforts have been made to identify the clinical predictors of TR. However, the exploration of molecular and biological markers is still at an early stage. To understand the TR condition and identify potential molecular and biological markers, we analyzed demographic information, clinical data, structural brain imaging data, and molecular brain imaging data in 7 Tesla magnetic resonance spectroscopy, from a first episode psychosis cohort that includes 138 patients. Age, gender, race, smoking status, duration of illness, and antipsychotic dosages were controlled in the analyses. We found that TR patients had a younger age at onset, more hospitalizations, more severe negative symptoms, a significant reduction in the volumes of the hippocampus (HP) and superior frontal gyrus (SFG), and a significant reduction in glutathione (GSH) levels in the anterior cingulate cortex (ACC), when compared to non-TR patients. The combination of multiple markers provided a better classification between TR and non-TR patients compared to any individual marker. Our study shows that ACC GSH, HP and SFG volumes, and age at onset could potentially be trait biomarkers for TR diagnosis, while hospitalization and negative symptoms could be used to evaluate the progression of the disease. Multimodal cohorts are essential in obtaining a comprehensive understanding of brain disorders.