Abstract
Major depressive disorder (MDD) represents a major social and economic health issue and constitutes a major risk factor for suicide. The molecular pathology of suicidal depression remains poorly understood, although it has been hypothesised that regulatory genomic processes are involved in the pathology of both MDD and suicidality. In this study, genome-wide patterns of DNA methylation were assessed in depressed suicide completers (n=20) and compared with non-psychiatric, sudden-death controls (n=20) using tissue from two cortical brain regions (Brodmann Area 11 (BA11) and Brodmann Area 25 (BA25)). Analyses focused on identifying differentially methylated regions (DMRs) associated with suicidal depression and epigenetic variation were explored in the context of polygenic risk scores for major depression and suicide. Weighted gene co-methylation network analysis was used to identify modules of co-methylated loci associated with depressed suicide completers and polygenic burden for MDD and suicide attempt. We identified a DMR upstream of the PSORS1C3 gene, subsequently validated using bisulfite pyrosequencing and replicated in a second set of suicide samples, which is characterised by significant hypomethylation in both cortical brain regions in MDD suicide cases. We also identified discrete modules of co-methylated loci associated with polygenic risk burden for suicide attempt, but not major depression. Suicide-associated co-methylation modules were enriched among gene networks implicating biological processes relevant to depression and suicidality, including nervous system development and mitochondria function. Our data suggest that there are coordinated changes in DNA methylation associated with suicide that may offer novel insights into the molecular pathology associated with depressed suicide completers.