We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further ...distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites.
cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant
cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In
trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both
cis- and
trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.
The expression microarray is a frequently used approach to study gene expression on a genome-wide scale. However, the data produced by the thousands of microarray studies published annually are ...confounded by "batch effects," the systematic error introduced when samples are processed in multiple batches. Although batch effects can be reduced by careful experimental design, they cannot be eliminated unless the whole study is done in a single batch. A number of programs are now available to adjust microarray data for batch effects prior to analysis. We systematically evaluated six of these programs using multiple measures of precision, accuracy and overall performance. ComBat, an Empirical Bayes method, outperformed the other five programs by most metrics. We also showed that it is essential to standardize expression data at the probe level when testing for correlation of expression profiles, due to a sizeable probe effect in microarray data that can inflate the correlation among replicates and unrelated samples.
Schizophrenia and bipolar disorder are complex psychiatric diseases with risks contributed by multiple genes. Dysregulation of gene expression has been implicated in these disorders, but little is ...known about such dysregulation in the human brain. We analyzed three transcriptome datasets from 394 postmortem brain tissue samples from patients with schizophrenia or bipolar disorder or from healthy control individuals without a known history of psychiatric disease. We built genome-wide coexpression networks that included microRNAs (miRNAs). We identified a coexpression network module that was differentially expressed in the brain tissue from patients compared to healthy control individuals. This module contained genes that were principally involved in glial and neural cell genesis and glial cell differentiation, and included schizophrenia risk genes carrying rare variants. This module included five miRNAs and 545 mRNAs, with six transcription factors serving as hub genes in this module. We found that the most connected transcription factor gene
, also identified on a genome-wide association study for bipolar disorder, could regulate the miRNA
and other putative target mRNAs. These regulatory relationships were replicated using PsychENCODE/BrainGVEX datasets and validated by knockdown and overexpression experiments in SH-SY5Y cells and human neural progenitor cells in vitro. Thus, we identified a brain gene expression module that was enriched for rare coding variants in genes associated with schizophrenia and that contained the putative bipolar disorder risk gene
The transcription factor
may be a key regulator of gene expression in this disease-associated gene coexpression module.
Caffeine produces mild psychostimulant and sometimes anxiogenic effects by antagonizing adenosine at A(1) and A(2A) receptors, and perhaps through interactions with other transmitter systems. ...Adenosine receptors are colocalized and functionally interact with dopamine receptors in the brain. Thus, functional polymorphisms in the genes for either adenosine or dopamine receptors may affect responses to caffeine. In this study, we examined associations between self-reported anxiogenic effects of caffeine and variation in the genes for A(2A) (ADORA2A) and DRD(2) (DRD2) receptors. Healthy male and female individuals (n=102), who consumed less than 300 mg caffeine per week, ingested capsules containing 0, 50, 150, and 450 mg caffeine under double-blind conditions in four separate experimental sessions. Subjective anxiety was measured before and at repeated times after capsules were consumed. At the 150 mg dose of caffeine, we found a significant association between caffeine-induced anxiety (Visual Analog Scales, VAS) and ADORA2A rs5751876 (1976C/T), rs2298383 (intron 1a) and rs4822492 (3'-flank), and DRD2 rs1110976 (intron 6). Caffeine-induced anxiety (VAS) was also associated with two-loci interactions of selected ADORA2A and DRD2 polymorphisms. The lowest dose of caffeine did not increase ratings of anxiety while the highest dose increased anxiety in the majority of subjects. These findings provide support for an association between an ADORA2A polymorphism and self-reported anxiety after a moderate dose of caffeine. It is likely that other ADORA2A and DRD2 polymorphisms also contribute to responses to caffeine.
Autism is a childhood neurodevelopmental disorder with a strong genetic component, yet the identification of autism susceptibility loci remains elusive. We investigated 180 autism probands and 372 ...control subjects by array comparative genomic hybridization (aCGH) using a 19K whole-genome tiling path bacterial artificial chromosome microarray to identify submicroscopic chromosomal rearrangements specific to autism. We discovered a recurrent 16p11.2 microdeletion in two probands with autism and none in controls. The deletion spans ∼500-kb and is flanked by ∼147-kb segmental duplications (SDs) that are >99% identical, a common characteristic of genomic disorders. We assessed the frequency of this new autism genomic disorder by screening an additional 532 probands and 465 controls by quantitative PCR and identified two more patients but no controls with the microdeletion, indicating a combined frequency of 0.6% (4/712 autism versus 0/837 controls; Fisher exact test P = 0.044). We confirmed all 16p11.2 deletions using fluorescence in situ hybridization, microsatellite analyses and aCGH, and mapped the approximate deletion breakpoints to the edges of the flanking SDs using a custom-designed high-density oligonucleotide microarray. Bioinformatic analysis localized 12 of the 25 genes within the microdeletion to nodes in one interaction network. We performed phenotype analyses and found no striking features that distinguish patients with the 16p11.2 microdeletion as a distinct autism subtype. Our work reports the first frequency, breakpoint, bioinformatic and phenotypic analyses of a de novo 16p11.2 microdeletion that represents one of the most common recurrent genomic disorders associated with autism to date.
Background One genetic mechanism known to be associated with autism spectrum disorders (ASD) is chromosomal abnormalities. The identification of copy number variants (CNV), i.e., microdeletions and ...microduplications that are undetectable at the level of traditional cytogenetic analysis, allows the potential association of submicroscopic chromosomal imbalances and human disease. Methods We performed array comparative genomic hybridization (aCGH) utilizing a 19K whole genome tiling path bacterial artificial chromosome (BAC) microarray on 397 unrelated subjects with autism spectrum disorder. Common CNV were excluded using a control group comprised of 372 individuals from the National Institute of Mental Health (NIMH) Genetics Initiative Control samples. Confirmation studies were performed on all remaining CNV using fluorescence in situ hybridization (FISH), microsatellite analysis, and/or quantitative polymerase chain reaction (PCR) analysis. Results A total of 51 CNV were confirmed in 46 ASD subjects. Three maternal interstitial duplications of 15q11-q13 known to be associated with ASD were identified. The other 48 CNV ranged in size from 189 kilobase (kb) to 5.5 megabase (Mb) and contained from 0 to ∼40 National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) genes. Seven CNV were de novo and 44 were inherited. Conclusions Fifty-one autism-specific CNV were identified in 46 of 397 ASD patients using a 19K BAC microarray for an overall rate of 11.6%. These microdeletions and microduplications cause gene dosage imbalance in 272 genes, many of which could be considered as candidate genes for autism.
AbstractSmooth pursuit eye tracking deficits are a promising intermediate phenotype for schizophrenia and possibly for psychotic disorders more broadly. The Bipolar-Schizophrenia Network on ...Intermediate Phenotypes (B-SNIP) consortium investigated the severity and familiality of different pursuit parameters across psychotic disorders. Probands with schizophrenia (N = 265), schizoaffective disorder (N = 178), psychotic bipolar disorder (N = 231), their first-degree relatives (N = 306, N = 217, N = 273, respectively) and healthy controls (N = 305) performed pursuit tracking tasks designed to evaluate sensorimotor and cognitive/predictive aspects of pursuit. Probands from all diagnostic groups were impaired on all pursuit measures of interest compared to controls (p < 0.001). Schizophrenia probands were more impaired than other proband groups on both early pursuit gain and predictive gain. Relatives with and without enhanced psychosis spectrum personality traits were impaired on initial eye acceleration, the most direct sensorimotor pursuit measure, but not on pursuit gain measures. This suggests that alterations in early sensorimotor function may track susceptibility to psychosis even in the absence of psychosis related personality traits. There were no differences in pursuit measures between relatives of the three proband groups. Familiality estimates of pursuit deficits indicate that early pursuit gain was more familial than predictive gain, which has been the most widely used measure in previous family studies of psychotic disorders. Thus, while disease-related factors may induce significant impairments of pursuit gain, especially in schizophrenia, the pattern of deficits in relatives and their familiality estimates suggest that alterations in sensorimotor function at pursuit onset may indicate increased susceptibility across psychotic disorders.
Autism is a complex childhood neurodevelopmental disorder with a strong genetic basis. Microdeletion or duplication of a approximately 500-700-kb genomic rearrangement on 16p11.2 that contains 24 ...genes represents the second most frequent chromosomal disorder associated with autism. The role of common and rare 16p11.2 sequence variants in autism etiology is unknown.
To identify common 16p11.2 variants with a potential role in autism, we performed association studies using existing data generated from three microarray platforms: Affymetrix 5.0 (777 families), Illumina 550 K (943 families), and Affymetrix 500 K (60 families). No common variants were identified that were significantly associated with autism. To look for rare variants, we performed resequencing of coding and promoter regions for eight candidate genes selected based on their known expression patterns and functions. In total, we identified 26 novel variants in autism: 13 exonic (nine non-synonymous, three synonymous, and one untranslated region) and 13 promoter variants. We found a significant association between autism and a coding variant in the seizure-related gene SEZ6L2 (12/1106 autism vs. 3/1161 controls; p = 0.018). Sez6l2 expression in mouse embryos was restricted to the spinal cord and brain. SEZ6L2 expression in human fetal brain was highest in post-mitotic cortical layers, hippocampus, amygdala, and thalamus. Association analysis of SEZ6L2 in an independent sample set failed to replicate our initial findings.
We have identified sequence variation in at least one candidate gene in 16p11.2 that may represent a novel genetic risk factor for autism. However, further studies are required to substantiate these preliminary findings.
Bipolar disorder is an often-severe mental health condition characterized by alternation between extreme mood states of mania and depression. Despite strong heritability and the recent identification ...of 64 common variant risk loci of small effect, pathophysiological mechanisms remain unknown. Here, we analyzed genome sequences from 41 multiply-affected pedigrees and identified variants in 741 genes with nominally significant linkage or association with bipolar disorder. These 741 genes overlapped known risk genes for neurodevelopmental disorders and clustered within gene networks enriched for synaptic and nuclear functions. The top variant in this analysis - prioritized by statistical association, predicted deleteriousness, and network centrality - was a missense variant in the gene encoding D-amino acid oxidase (DAO
). Heterologous expression of DAO
in human cells resulted in decreased DAO protein abundance and enzymatic activity. In a knock-in mouse model of DAO
, Dao
, we similarly found decreased DAO protein abundance in hindbrain regions, as well as enhanced stress susceptibility and blunted behavioral responses to pharmacological inhibition of N-methyl-D-aspartate receptors (NMDARs). RNA sequencing of cerebellar tissue revealed that Dao
resulted in decreased expression of two gene networks that are enriched for synaptic functions and for genes expressed, respectively, in Purkinje neurons or granule neurons. These gene networks were also down-regulated in the cerebellum of patients with bipolar disorder compared to healthy controls and were enriched for additional rare variants associated with bipolar disorder risk. These findings implicate dysregulation of NMDAR signaling and of gene expression in cerebellar neurons in bipolar disorder pathophysiology and provide insight into its genetic architecture.