Abstract
Only 50% of patients with depression respond to the first antidepressant drug administered. Thus, biomarkers for prediction of antidepressant responses are needed, as predicting which patients will not respond to antidepressants can optimize selection of alternative therapies. We aimed to identify biomarkers that could predict antidepressant responsiveness using a novel data-driven approach based on statistical pattern recognition. We retrospectively divided patients with major depressive disorder into antidepressant responder and non-responder groups. Comprehensive gene expression analysis was performed using peripheral blood without narrowing the genes. We designed a classifier according to our own discrete Bayes decision rule that can handle categorical data. Nineteen genes showed differential expression in the antidepressant non-responder group (n = 15) compared to the antidepressant responder group (n = 15). In the training sample of 30 individuals, eight candidate genes had significantly altered expression according to quantitative real-time polymerase chain reaction. The expression of these genes was examined in an independent test sample of antidepressant responders (n = 22) and non-responders (n = 12). Using the discrete Bayes classifier with the HERC5, IFI6, and IFI44 genes identified in the training set yielded 85% discrimination accuracy for antidepressant responsiveness in the 34 test samples. Pathway analysis of the RNA sequencing data for antidepressant responsiveness identified that hypercytokinemia- and interferon-related genes were increased in non-responders. Disease and biofunction analysis identified changes in genes related to inflammatory and infectious diseases, including coronavirus disease. These results strongly suggest an association between antidepressant responsiveness and inflammation, which may be useful for future treatment strategies for depression.