Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes.

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Targeted sequencing of the LRRTM gene family in suicide attempters with bipolar disorder.

Glutamatergic signaling is the primary excitatory neurotransmission pathway in the brain, and its relationship to neuropsychiatric disorders is of considerable interest. Our previous attempted suicide genome-wide association study, and numerous studies investigating gene expression, genetic variation, and DNA methylation have implicated aberrant glutamatergic signaling in suicide risk. The glutamatergic pathway gene LRRTM4 was an associated gene identified in our attempted suicide genome-wide association study, with association support seen primarily in females. Recent evidence has also shown that glutamatergic signaling is partly regulated by sex-related hormones. The LRRTM gene family encodes neuronal leucine-rich transmembrane proteins that localize to and promote glutamatergic synapse development. In this study, we sequenced the coding and regulatory regions of all four LRRTM gene members plus a large intronic region of LRRTM4 in 476 bipolar disorder suicide attempters and 473 bipolar disorder nonattempters. We identified two male-specific variants, one female- and five male-specific haplotypes significantly associated with attempted suicide in LRRTM4. Furthermore, variants within significant haplotypes may be brain expression quantitative trait loci for LRRTM4 and some of these variants overlap with predicted hormone response elements. Overall, these results provide supporting evidence for a sex-specific association of genetic variation in LRRTM4 with attempted suicide.

Differentially methylated regions in bipolar disorder and suicide.

The addition of a methyl group to, typically, a cytosine-guanine dinucleotide (CpG) creates distinct DNA methylation patterns across the genome that can regulate gene expression. Aberrant DNA methylation of CpG sites has been associated with many psychiatric disorders including bipolar disorder (BD) and suicide. Using the SureSelectXT system, Methyl-Seq, we investigated the DNA methylation status of CpG sites throughout the genome in 50 BD individuals (23 subjects who died by suicide and 27 subjects who died from other causes) and 31 nonpsychiatric controls. We identified differentially methylated regions (DMRs) from three analyses: (a) BD subjects compared to nonpsychiatric controls (BD-NC), (b) BD subjects who died by suicide compared to nonpsychiatric controls (BDS-NC), and (c) BDS subjects compared to BD subjects who died from other causes (BDS-BDNS). One DMR from the BDS-NC analysis, located in ARHGEF38, was significantly hypomethylated (23.4%) in BDS subjects. This finding remained significant after multiple testing (PBootstrapped = 9.0 × 10-3 ), was validated using pyrosequencing, and was more significant in males. A secondary analysis utilized Ingenuity Pathway Analysis to identify enrichment in nominally significant DMRs. This identified an association with several pathways including axonal guidance signaling, calcium signaling, ß-adrenergic signaling, and opioid signaling. Our comprehensive study provides further support that DNA methylation alterations influence the risk for BD and suicide. However, further investigation is required to confirm these associations and identify their functional consequences.

Whole-gene sequencing investigation of SAT1 in attempted suicide.

Suicidal behavior imposes a tremendous cost, with current US estimates reporting approximately 1.3 million suicide attempts and more than 40,000 suicide deaths each year. Several recent research efforts have identified an association between suicidal behavior and the expression level of the spermidine/spermine N1-acetyltransferase 1 (SAT1) gene. To date, several SAT1 genetic variants have been inconsistently associated with altered gene expression and/or directly with suicidal behavior. To clarify the role SAT1 genetic variation plays in suicidal behavior risk, we present a whole-gene sequencing effort of SAT1 in 476 bipolar disorder subjects with a history of suicide attempt and 473 subjects with bipolar disorder but no suicide attempts. Agilent SureSelect target enrichment was used to sequence all exons, introns, promoter regions, and putative regulatory regions identified from the ENCODE project within 10 kb of SAT1. Individual variant, haplotype, and collapsing variant tests were performed. Our results identified no variant or assessed region of SAT1 that showed a significant association with attempted suicide, nor did any assessment show evidence for replication of previously reported associations. Overall, no evidence for SAT1 sequence variation contributing to the risk for attempted suicide could be identified. It is possible that past associations of SAT1 expression with suicidal behavior arise from variation not captured in this study, or that causal variants in the region are too rare to be detected within our sample. Larger sample sizes and broader sequencing efforts will likely be required to identify the source of SAT1 expression level associations with suicidal behavior. © 2016 Wiley Periodicals, Inc.

A targeted sequencing study of glutamatergic candidate genes in suicide attempters with bipolar disorder.

Suicidal behavior has been shown to have a heritable component that is partly driven by psychiatric disorders [Brent and Mann, 2005]. However, there is also an independent factor contributing to the heritability of suicidal behavior. We previously conducted a genome-wide association study (GWAS) of bipolar suicide attempters and bipolar non-attempters to assess this independent factor [Willour et al., 2012]. This GWAS implicated glutamatergic neurotransmission in attempted suicide. In the current study, we have conducted a targeted next-generation sequencing study of the glutamatergic N-methyl-D-aspartate (NMDA) receptor, neurexin, and neuroligin gene families in 476 bipolar suicide attempters and 473 bipolar non-attempters. The goal of this study was to gather sequence information from coding and regulatory regions of these glutamatergic genes to identify variants associated with attempted suicide. We identified 186 coding variants and 4,298 regulatory variants predicted to be functional in these genes. No individual variants were overrepresented in cases or controls to a degree that was statistically significant after correction for multiple testing. Additionally, none of the gene-level results were statistically significant following correction. While this study provides no direct support for a role of the examined glutamatergic candidate genes, further sequencing in expanded gene sets and datasets will be required to ultimately determine whether genetic variation in glutamatergic signaling influences suicidal behavior. © 2016 Wiley Periodicals, Inc.

Association study of serotonin pathway genes in attempted suicide.

Epidemiological studies, such as family, twin, and adoption studies, demonstrate the presence of a heritable component to both attempted and completed suicide. Some of this heritability is accounted for by the presence of comorbid psychiatric disorders, but the evidence also indicates that a portion of this heritability is specific to suicidality. The serotonergic system has been studied extensively in this phenotype, but findings have been inconsistent, possibly due to the presence of multiple susceptibility variants and/or gene-gene interactions. In this study, we genotyped 174 tag and coding single nucleotide polymorphisms (SNPs) from 17 genes within the serotonin pathway on 516 subjects with a major mood disorder and a history of a suicide attempt (cases) and 515 healthy controls, with the goal of capturing the common genetic variation across each of these candidate genes. We tested the 174 markers in single-SNP, haplotype, gene-based, and epistasis analyses. While these association analyses identified multiple marginally significant SNPs, haplotypes, genes, and interactions, none of them survived correction for multiple testing. Additional studies, including assessment in larger sample sets and deep resequencing to identify rare causal variants, may be required to fully understand the role that the serotonin pathway plays in suicidal behavior.

Distinct patterns of altered quantitative T1ρ and functional BOLD response associated with history of suicide attempts in bipolar disorder.

Despite the high risk for suicide, relatively few studies have explored the relationship between suicide and brain imaging measures in bipolar disorder. In addition, fewer studies have explored the possibility that altered brain metabolism may be associated with suicide attempt. To begin to fill in these gaps, we evaluated functional (task based fMRI) and metabolic (quantitative T1ρ) differences associated with suicide attempt in participants with bipolar disorder. Thirty-nine participants with bipolar disorder underwent fMRI during a flashing checkerboard task and 27 also underwent quantitative T1ρ. The relationship between neuroimaging and history of suicide attempt was tested using multiple regression while adjusting for age, sex, and current mood state. Differences between two measures of suicide attempt (binary: yes/no and continuous: number of attempts) were quantified using the corrected Akaike Information Criterion. Participants who had attempted suicide had greater fMRI task-related activation in visual areas and the cerebellum. The number of suicide attempts was associated with a difference in BOLD response in the amygdala, prefrontal cortex, and cerebellum. Increased quantitative T1ρ was associated with number of suicide attempts in limbic, basal ganglia, and prefrontal cortex regions. This study is a secondary analysis with a modest sample size. Differences between measures of suicide history may be due to differences in statistical power. History of suicide was associated with limbic, prefrontal, and cerebellar alterations. Results comparing those with and without suicide attempts differed from results using number of suicide attempts, suggesting that these variables have different neurobiological underpinnings.

Monoamine oxidase A regulates antisocial personality in whites with no history of physical abuse.

OBJECTIVE: Preclinical and human family studies clearly link monoamine oxidase A (MAOA) to aggression and antisocial personality (ASP). The 30-base pair variable number tandem repeat in the MAOA promoter regulates MAOA levels, but its effects on ASP in humans are unclear. METHODS: We evaluated the association of the variable number tandem repeat of the MAOA promoter with Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, ASP disorder (ASPD) traits in a community sample of 435 participants from the Hopkins Epidemiology of Personality Disorders Study. RESULTS: We did not find an association between the activity of the MAOA allele and ASPD traits; however, among whites, when subjects with a history of childhood physical abuse were excluded, the remaining subjects with low-activity alleles had ASPD trait counts that were 41% greater than those with high-activity alleles (P < .05). CONCLUSION: The high-activity MAOA allele is protective against ASP among whites with no history of physical abuse, lending support to a link between MAOA expression and antisocial behavior.

Assessment of suicide attempt and death in bipolar affective disorder: a combined clinical and genetic approach.

Bipolar disorder (BP) suicide death rates are 10-30 times greater than the general population, likely arising from environmental and genetic risk factors. Though suicidal behavior in BP has been investigated, studies have not addressed combined clinical and genetic factors specific to suicide death. To address this gap, a large, harmonized BP cohort was assessed to identify clinical risk factors for suicide death and attempt which then directed testing of underlying polygenic risks. 5901 individuals of European ancestry were assessed: 353 individuals with BP and 2498 without BP who died from suicide (BPS and NBPS, respectively) from a population-derived sample along with a volunteer-derived sample of 799 individuals with BP and a history of suicide attempt (BPSA), 824 individuals with BP and no prior attempts (BPNSA), and 1427 individuals without several common psychiatric illnesses per self-report (C). Clinical and subsequent directed genetic analyses utilized multivariable logistic models accounting for critical covariates and multiple testing. There was overrepresentation of diagnosis of PTSD (OR = 4.9, 95%CI: 3.1-7.6) in BPS versus BPSA, driven by female subjects. PRS assessments showed elevations in BPS including PTSD (OR = 1.3, 95%CI:1.1-1.5, versus C), female-derived ADHD (OR = 1.2, 95%CI:1.1-1.4, versus C), and male insomnia (OR = 1.4, 95%CI: 1.1-1.7, versus BPSA). The results provide support from genetic and clinical standpoints for dysregulated traumatic response particularly increasing risk of suicide death among individuals with BP of Northern European ancestry. Such findings may direct more aggressive treatment and prevention of trauma sequelae within at-risk bipolar individuals.

Case-control association study of TGOLN2 in attempted suicide.

Family, twin, and adoption studies provide convincing evidence for a genetic contribution to suicidal behavior. The heritability for suicidal behavior depends in part on the transmission of psychiatric disorders, such as mood disorders and substance use disorders, but is also partly independent of them. Three linkage studies using the attempted suicide phenotype in pedigrees with bipolar disorder, major depression, or alcoholism have provided consistent evidence that 2p11-12 harbors a susceptibility gene for attempted suicide. A microarray expression study using postmortem brain samples has implicated a gene from the 2p11-12 candidate region, the trans-Golgi network protein 2 (TGOLN2) gene, as being consistently up-regulated in suicide cases as compared to controls. Here, we present a TGOLN2 case-control association study using nine single nucleotide polymorphisms (SNPs). These nine SNPs, which include seven tag SNPs and two coding SNPs, have been genotyped in 517 mood disorder subjects with a history of attempted suicide and 515 normal controls. Allelic and genotypic analyses of the case-control sample did not provide evidence for association with the attempted suicide phenotype. Eight of the nine SNPs provided supportive evidence for association (P-values ranging from 0.008 to 0.03) when we compared the attempted suicide cases with a history of alcoholism to the attempted suicide cases without a history of alcoholism. However, this association finding was not replicated in an independent sample. Taken together, these analyses do not provide support for the hypothesis that common genetic variation in TGOLN2 contributes significantly to the risk for attempted suicide in subjects with major mood disorders.

Male-specific association of the 2p25 region with suicide attempt in bipolar disorder.

We previously conducted a genome-wide association study (GWAS) of attempted suicide within bipolar disorder, which implicated common variation in the 2p25 region primarily in males. The top association signal from our GWAS occurred in an intergenic region of 2p25 (p = 5.07 × 10-8) and was supported by two independent studies. In the current study, to better characterize the association of the 2p25 region with attempted suicide, we sequenced the entire 350kb 2p25 region in 476 bipolar suicide attempters and 473 bipolar non-attempters using targeted next-generation sequencing. This fine-mapping project identified 4,681 variants in the 2p25 region. We performed both gene-level and individual-variant tests on our sequencing results and identified 375 variants which were nominally significant (p < 0.05) and three common variants that were significantly associated with attempted suicide in males (corrected p = 0.035, odds ratio (OR) = 2.13). These three variants are in strong linkage disequilibrium with the top variant from our GWAS. Our top five variants are also predicted expression quantitative trait loci (eQTL) for three genes in the 2p25 region based on publicly available brain expression databases. Our sequencing and eQTL data implicate these three genes - SH3YL1, ACP1, and FAM150B - and three additional pathways - androgen receptor, Wnt signaling, and glutamatergic/GABAergic signaling - in the association of the 2p25 region with suicide. The current study provides additional support for an association of the 2p25 region with attempted suicide in males and identifies several candidate genes and pathways that warrant further investigation to understand their role in suicidal behavior.

Family-based association study of Neuregulin 1 with psychotic bipolar disorder.

The Neuregulin 1 gene (NRG1) has been associated with schizophrenia, and, to a lesser extent, with bipolar disorder (BP). We investigated the association of NRG1 with BP in a large family sample, and then performed analyses according to the presence of psychotic features or mood-incongruent psychotic features. We genotyped 116 tagSNPs and four Icelandic "core" SNPs in 1,199 subjects from 314 nuclear families. Of 515 BP offspring, 341 had psychotic features, and 103 had mood-incongruent psychotic features. In single-marker and sliding window haplotype analyses using FBAT, there was little association using the standard BP or mood-incongruent psychotic BP phenotypes, but stronger signals were seen in the psychotic BP phenotype. The most significant associations with psychotic BP were in haplotypes within the 5' "core" region. The strongest global P-value was across three SNPs: NRG241930-NRG243177-rs7819063 (P = 0.0016), with an undertransmitted haplotype showing an individual P = 0.0007. The most significant individual haplotype was an undertransmitted two-allele subset of the above (NRG243177-rs7819063, P = 0.0004). Additional associations with psychotic BP were found across six SNPs in a 270 kb central region of the gene. The most 3' of these, rs7005606 (P = 0.0029), is located approximately 4 kb from the type I NRG1 isoform promoter. In sum, our study suggests that NRG1 may be specifically associated with the psychotic subset of BP; however, our results should be interpreted cautiously since they do not meet correction for multiple testing and await independent replication.

Common and rare variants of DAOA in bipolar disorder.

The D-amino acid oxidase activator (DAOA, previously known as G72) gene, mapped on 13q33, has been reported to be genetically associated with bipolar disorder (BP) in several populations. The consistency of associated variants is unclear and rare variants in exons of the DAOA gene have not been investigated in psychiatric diseases. We employed a conditional linkage method-STatistical Explanation for Positional Cloning (STEPC) to evaluate whether any associated single nucleotide polymorphisms (SNPs) account for the evidence of linkage in a pedigree series that previously has been linked to marker D13S779 at 13q33. We also performed an association study in a sample of 376 Caucasian BP parent-proband trios by genotyping 38 common SNPs in the gene region. Besides, we resequenced coding regions and flanking intronic sequences of DAOA in 555 Caucasian unrelated BP patients and 564 mentally healthy controls, to identify putative functional rare variants that may contribute to disease. One SNP rs1935058 could "explain" the linkage signal in the family sample set (P = 0.055) using STEPC analysis. No significant allelic association was detected in an association study by genotyping 38 common SNPs in 376 Caucasian BP trios. Resequencing identified 53 SNPs, of which 46 were novel SNPs. There was no significant excess of rare variants in cases relative to controls. Our results suggest that DAOA does not have a major effect on BP susceptibility. However, DAOA may contribute to bipolar susceptibility in some specific families as evidenced by the STEPC analysis.

Further development of YBOCS dimensions in the OCD Collaborative Genetics study: symptoms vs. categories.

Despite progress in identifying homogeneous subphenotypes of obsessive-compulsive disorder (OCD) through factor analysis of the Yale Brown Obsessive-Compulsive Scale Symptom Checklist (YBOCS-SC), prior solutions have been limited by a reliance on presupposed symptom categories rather than discrete symptoms. Furthermore, there have been few attempts to evaluate the familiality of OCD symptom dimensions. The purpose of this study was to extend prior work by this collaborative group in category-based dimensions by conducting the first-ever exploratory dichotomous factor analysis using individual OCD symptoms, comparing these results to a refined category-level solution, and testing the familiality of derived factors. Participants were 485 adults in the six-site OCD Collaborative Genetics Study, diagnosed with lifetime OCD using semi-structured interviews. YBOCS-SC data were factor analyzed at both the individual item and symptom category levels. Factor score intraclass correlations were calculated using a subsample of 145 independent affected sib pairs. The item- and category-level factor analyses yielded nearly identical 5-factor solutions. While significant sib-sib associations were found for four of the five factors, Hoarding and Taboo Thoughts were the most robustly familial (r ICC>or=0.2). This report presents considerable converging evidence for a five-factor structural model of OCD symptoms, including separate factor analyses employing individual symptoms and symptom categories, as well as sibling concordance. The results support investigation of this multidimensional model in OCD genetic linkage studies.

Clock genes may influence bipolar disorder susceptibility and dysfunctional circadian rhythm.

Several previous studies suggest that dysfunction of circadian rhythms may increase susceptibility to bipolar disorder (BP). We conducted an association study of five circadian genes (CRY2, PER1-3, and TIMELESS) in a family collection of 36 trios and 79 quads (Sample I), and 10 circadian genes (ARNTL, ARNTL2, BHLHB2, BHLHB3, CLOCK, CRY1, CSNK1D, CSNK1E, DBP, and NR1D1) in an extended family collection of 70 trios and 237 quads (Sample II), which includes the same 114 families but not necessarily the same individuals as Sample I. In Sample II, the Sibling-Transmission Disequilibrium Test (sib-tdt) analysis showed nominally significant association of BP with three SNPs within or near the CLOCK gene (rs534654, P = 0.0097; rs6850524, P = 0.012; rs4340844, P = 0.015). In addition, SNPs in the ARNTL2, CLOCK, DBP, and TIMELESS genes and haplotypes in the ARNTL, CLOCK, CSNK1E, and TIMELESS genes showed suggestive evidence of association with several circadian phenotypes identified in BP patients. However, none of these associations reached gene-wide or experiment-wide significance after correction for multiple-testing. A multi-locus interaction between rs6442925 in the 5' upstream of BHLHB2, rs1534891 in CSNK1E, and rs534654 near the 3' end of the CLOCK gene, however, is significantly associated with BP (P = 0.00000172). It remains significant after correcting for multiple testing using the False Discovery Rate method. Our results indicate an interaction between three circadian genes in susceptibility to bipolar disorder.

Neurotransmission and bipolar disorder: a systematic family-based association study.

Neurotransmission pathways/systems have been proposed to be involved in the pathophysiology and treatment of bipolar disorder for over 40 years. In order to test the hypothesis that common variants of genes in one or more of five neurotransmission systems confer risk for bipolar disorder, we analyzed 1,005 tag single nucleotide polymorphisms in 90 genes from dopaminergic, serotonergic, noradrenergic, GABAergic, and glutamatergic neurotransmitter systems in 101 trios and 203 quads from Caucasian bipolar families. Our sample has 80% power to detect ORs >or= 1.82 and >or=1.57 for minor allele frequencies of 0.1 and 0.5, respectively. Nominally significant allelic and haplotypic associations were found for genes from each neurotransmission system, with several reaching gene-wide significance (allelic: GRIA1, GRIN2D, and QDPR; haplotypic: GRIN2C, QDPR, and SLC6A3). However, none of these associations survived correction for multiple testing in an individual system, or in all systems considered together. Significant single nucleotide polymorphism associations were not found with sub-phenotypes (alcoholism, psychosis, substance abuse, and suicide attempts) or significant gene-gene interactions. These results suggest that, within the detectable odds ratios of this study, common variants of the selected genes in the five neurotransmission systems do not play major roles in influencing the risk for bipolar disorder or comorbid sub-phenotypes.

Family-based SNP association study on 8q24 in bipolar disorder.

Previous linkage studies have identified chromosome 8q24 as a promising positional candidate region to search for bipolar disorder (BP) susceptibility genes. We, therefore, sought to identify BP susceptibility genes on chromosome 8q24 using a family-based association study of a dense panel of SNPs selected to tag the known common variation across the region of interest. A total of 1,458 SNPs across 16 Mb of 8q24 were examined in 3,512 subjects, 1,954 of whom were affected with BP, from 737 multiplex families. Single-locus tests were carried out with FBAT and Geno-PDT, and multi-locus test were carried out with HBAT and multi-locus Geno-PDT. None of the SNPs were associated with BP in the single-locus tests at a level that exceeded our threshold for study-wide significance (P < 3.00 x 10(-5)). However, there was consistent evidence at our threshold for the suggestive level (P < 7.00 x 10(-4)) from both the single locus and multi-locus tests of associations with SNPs in the genes ADCY8, ST3GAL1, and NSE2. Multi-locus analyses suggested joint effects between ADCY8 and ST3GAL1 (P = 3.00 x 10(-4)), with at least one copy of the "high risk" allele required at both genes for association with BP, consistent with a jointly dominant-dominant model of action. These findings with ADCY8 and ST3GAL1 warrant further investigation in order to confirm the observed associations and their functional significance for BP susceptibility.

Chronic social stress induces DNA methylation changes at an evolutionary conserved intergenic region in chromosome X.

Chronic stress resulting from prolonged exposure to negative life events increases the risk of mood and anxiety disorders. Although chronic stress can change gene expression relevant for behavior, molecular regulators of this change have not been fully determined. One process that could play a role is DNA methylation, an epigenetic process whereby a methyl group is added onto nucleotides, predominantly cytosine in the CpG context, and which can be induced by chronic stress. It is unknown to what extent chronic social defeat, a model of human social stress, influences DNA methylation patterns across the genome. Our study addressed this question by using a targeted-capture approach called Methyl-Seq to investigate DNA methylation patterns of the dentate gyrus at putative regulatory regions across the mouse genome from mice exposed to 14 days of social defeat. Findings were replicated in independent cohorts by bisulfite-pyrosequencing. Two differentially methylated regions (DMRs) were identified. One DMR was located at intron 9 of Drosha, and it showed reduced methylation in stressed mice. This observation replicated in one of two independent cohorts. A second DMR was identified at an intergenic region of chromosome X, and methylation in this region was increased in stressed mice. This methylation difference replicated in two independent cohorts and in Major Depressive Disorder (MDD) postmortem brains. These results highlight a region not previously known to be differentially methylated by chronic social defeat stress and which may be involved in MDD.