Identification of the BRD1 interaction network and its impact on mental disorder risk

BACKGROUND: The bromodomain containing 1 (BRD1) gene has been implicated with transcriptional regulation, brain development, and susceptibility to schizophrenia and bipolar disorder. To advance the understanding of BRD1 and its role in mental disorders, we characterized the protein and chromatin interactions of the BRD1 isoforms, BRD1-S and BRD1-L. METHODS: Stable human cell lines expressing epitope tagged BRD1-S and BRD1-L were generated and used as discovery systems for identifying protein and chromatin interactions. Protein-protein interactions were identified using co-immunoprecipitation followed by mass spectrometry and chromatin interactions were identified using chromatin immunoprecipitation followed by next generation sequencing. Gene expression profiles and differentially expressed genes were identified after upregulating and downregulating BRD1 expression using microarrays. The presented functional molecular data were integrated with human genomic and transcriptomic data using available GWAS, exome-sequencing datasets as well as spatiotemporal transcriptomic datasets from the human brain. RESULTS: We present several novel protein interactions of BRD1, including isoform-specific interactions as well as proteins previously implicated with mental disorders. By BRD1-S and BRD1-L chromatin immunoprecipitation followed by next generation sequencing we identified binding to promoter regions of 1540 and 823 genes, respectively, and showed correlation between BRD1-S and BRD1-L binding and regulation of gene expression. The identified BRD1 interaction network was found to be predominantly co-expressed with BRD1 mRNA in the human brain and enriched for pathways involved in gene expression and brain function. By interrogation of large datasets from genome-wide association studies, we further demonstrate that the BRD1 interaction network is enriched for schizophrenia risk. CONCLUSION: Our results show that BRD1 interacts with chromatin remodeling proteins, e.g. PBRM1, as well as histone modifiers, e.g. MYST2 and SUV420H1. We find that BRD1 primarily binds in close proximity to transcription start sites and regulates expression of numerous genes, many of which are involved with brain development and susceptibility to mental disorders. Our findings indicate that BRD1 acts as a regulatory hub in a comprehensive schizophrenia risk network which plays a role in many brain regions throughout life, implicating e.g. striatum, hippocampus, and amygdala at mid-fetal stages. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13073-016-0308-x) contains supplementary material, which is available to authorized users

Identification of common genetic risk variants for autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD

A genetic investigation of sex bias in the prevalence of attention-deficit/hyperactivity disorder

Background Attention-deficit/hyperactivity disorder (ADHD) shows substantial heritability and is 2-7 times more common in males than females. We examined two putative genetic mechanisms underlying this sex bias: sex-specific heterogeneity and higher burden of risk in female cases. Methods We analyzed genome-wide autosomal common variants from the Psychiatric Genomics Consortium and iPSYCH Project (20,183 cases, 35,191 controls) and Swedish populationregister data (N=77,905 cases, N=1,874,637 population controls). Results Genetic correlation analyses using two methods suggested near complete sharing of common variant effects across sexes, with rg estimates close to 1. Analyses of population data, however, indicated that females with ADHD may be at especially high risk of certain comorbid developmental conditions (i.e. autism spectrum disorder and congenital malformations), potentially indicating some clinical and etiological heterogeneity. Polygenic risk score (PRS) analysis did not support a higher burden of ADHD common risk variants in female cases (OR=1.02 [0.98-1.06], p=0.28). In contrast, epidemiological sibling analyses revealed that the siblings of females with ADHD are at higher familial risk of ADHD than siblings of affected males (OR=1.14, [95% CI: 1.11-1.18], p=1.5E-15). Conclusions Overall, this study supports a greater familial burden of risk in females with ADHD and some clinical and etiological heterogeneity, based on epidemiological analyses. However, molecular genetic analyses suggest that autosomal common variants largely do not explain the sex bias in ADHD prevalence

A Bayesian framework for multiple trait colocalization from summary association statistics

Motivation:Most genetic variants implicated in complex diseases by genome-wide association studies (GWAS) are non-coding, making it challenging to understand the causative genes involved in disease. Integrating external information such as quantitative trait locus (QTL) mapping of molecular traits (e.g. expression, methylation) is a powerful approach to identify the subset of GWAS signals explained by regulatory effects. In particular, expression QTLs (eQTLs) help pinpoint the responsible gene among the GWAS regions that harbor many genes, while methylation QTLs (mQTLs) help identify the epigenetic mechanisms that impact gene expression which in turn affect disease risk. In this work, we propose multiple-trait-coloc (moloc), a Bayesian statistical framework that integrates GWAS summary data with multiple molecular QTL data to identify regulatory effects at GWAS risk loci. Results:We applied moloc to schizophrenia (SCZ) and eQTL/mQTL data derived from human brain tissue and identified 52 candidate genes that influence SCZ through methylation. Our method can be applied to any GWAS and relevant functional data to help prioritize disease associated genes. Availability and implementation: moloc is available for download as an R package (https://github.com/clagiamba/moloc). We also developed a web site to visualize the biological findings (icahn.mssm.edu/moloc). The browser allows searches by gene, methylation probe and scenario of interest. Supplementary information:Supplementary data are available at Bioinformatics online

A Bayesian framework for multiple trait colocalization from summary association statistics

Motivation:Most genetic variants implicated in complex diseases by genome-wide association studies (GWAS) are non-coding, making it challenging to understand the causative genes involved in disease. Integrating external information such as quantitative trait locus (QTL) mapping of molecular traits (e.g. expression, methylation) is a powerful approach to identify the subset of GWAS signals explained by regulatory effects. In particular, expression QTLs (eQTLs) help pinpoint the responsible gene among the GWAS regions that harbor many genes, while methylation QTLs (mQTLs) help identify the epigenetic mechanisms that impact gene expression which in turn affect disease risk. In this work, we propose multiple-trait-coloc (moloc), a Bayesian statistical framework that integrates GWAS summary data with multiple molecular QTL data to identify regulatory effects at GWAS risk loci. Results:We applied moloc to schizophrenia (SCZ) and eQTL/mQTL data derived from human brain tissue and identified 52 candidate genes that influence SCZ through methylation. Our method can be applied to any GWAS and relevant functional data to help prioritize disease associated genes. Availability and implementation: moloc is available for download as an R package (https://github.com/clagiamba/moloc). We also developed a web site to visualize the biological findings (icahn.mssm.edu/moloc). The browser allows searches by gene, methylation probe and scenario of interest. Supplementary information:Supplementary data are available at Bioinformatics online

Dysregulation of miRNA-9 in a Subset of Schizophrenia Patient-Derived Neural Progenitor Cells

Converging evidence indicates that microRNAs (miRNAs) may contribute to disease risk for schizophrenia (SZ). We show that microRNA-9 (miR-9) is abundantly expressed in control neural progenitor cells (NPCs) but also significantly downregulated in a subset of SZ NPCs. We observed a strong correlation between miR-9 expression and miR-9 regulatory activity in NPCs as well as between miR-9 levels/activity, neural migration, and diagnosis. Overexpression of miR-9 was sufficient to ameliorate a previously reported neural migration deficit in SZ NPCs, whereas knockdown partially phenocopied aberrant migration in control NPCs. Unexpectedly, proteomic- and RNA sequencing (RNA-seq)-based analysis revealed that these effects were mediated primarily by small changes in expression of indirect miR-9 targets rather than large changes in direct miR-9 targets; these indirect targets are enriched for migration-associated genes. Together, these data indicate that aberrant levels and activity of miR-9 may be one of the many factors that contribute to SZ risk, at least in a subset of patients

Cell-specific histone modification maps in the human frontal lobe link schizophrenia risk to the neuronal epigenome

Risk variants for schizophrenia affect more than 100 genomic loci, yet cell- and tissue-specific roles underlying disease liability remain poorly characterized. We have generated for two cortical areas implicated in psychosis, the dorsolateral prefrontal cortex and anterior cingulate cortex, 157 reference maps from neuronal, neuron-depleted and bulk tissue chromatin for two histone marks associated with active promoters and enhancers, H3-trimethyl-Lys4 (H3K4me3) and H3-acetyl-Lys27 (H3K27ac). Differences between neuronal and neuron-depleted chromatin states were the major axis of variation in histone modification profiles, followed by substantial variability across subjects and cortical areas. Thousands of significant histone quantitative trait loci were identified in neuronal and neuron-depleted samples. Risk variants for schizophrenia, depressive symptoms and neuroticism were significantly over-represented in neuronal H3K4me3 and H3K27ac landscapes. Our Resource, sponsored by PsychENCODE and CommonMind, highlights the critical role of cell-type-specific signatures at regulatory and disease-associated noncoding sequences in the human frontal lobe