Many genes are linked to psychiatric disorders, but genome-wide association studies (GWAS) and differential gene expression (DGE) analyses in post-mortem brain tissue often implicate distinct gene sets. This disconnect impedes therapeutic development, which relies on integrating genetic and genomic insights. We address this issue using a novel multivariate technique that reduces DGE bias by leveraging gene co-expression networks and controlling for confounds such as drug exposure. Deep RNA sequencing was performed in bulk post-mortem sgACC from individuals with bipolar disorder (BD
N=35), major depression (MDD
N=51), schizophrenia (SCZ
N=44), and controls (N=55). Toxicology data dimensionality was reduced using multiple correspondence analysis
case-control gene expression was then analyzed using 1) traditional DGE and 2) group regularized canonical correlation analysis (GRCCA) - a multivariate regression method that accounts for feature interdependence. Gene set enrichment analyses compared results with established neuropsychiatric risk genes, gene ontology pathways, and cell type enrichments. GRCCA revealed a significant association with SCZ (