Integrative Gene and Isoform Co-expression Networks Reveal Regulatory Rewiring in Stress-related Psychiatric Disorders [PBMC]
Isoform-specific expression patterns have been implicated in stress-related psychiatric disorders like major  depressive disorder (MDD), yet the extent of their involvement and their interrelationships remain unclear. We constructed co-expression networks for individuals affected (n=210, 81% with depressive symptoms) and unaffected (n=95) by stress-related psychiatric disorders. We incorporated total gene expression (TE) and isoform ratio (IR) data and validated the inferred networks using advanced graph generation techniques. Our analysis revealed distinct network topology and structure between the two groups. Investigation of the 127 shared hubs (degree >= 10) found that these hubs exhibit co-regulatory patterns unique to each network. The affected individuals' network also contained 61 hub nodes with a minimum absolute fold increase of two in connectivity compared to the unaffected individuals' network. Notably, 49% of these hubs showed evidence for association with psychiatric disorders. Gene Ontology enrichment analysis revealed distinct biological processes associated with hubs, such as mRNA processing for affected and immune response and cell adhesion for unaffected individuals. Enrichment analysis of GWAS loci further supported these findings. Analysis of the isoform-specific nodes showed distinct protein-protein interactions compared to gene-level analysis. Overall design: The total RNA-seq of PBMCs from 229 participants with stress-related mental disorders and 70 healthy controls was determined.