Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data

Genome-wide expression profiling is a powerful tool for implicating novel gene ensembles in cellular mechanisms of health and disease. The most popular platform for genome-wide expression profiling is the Affymetrix GeneChip. However, its selection of probes relied on earlier genome and transcriptome annotation which is significantly different from current knowledge. The resultant informatics problems have a profound impact on analysis and interpretation the data. Here, we address these critical issues and offer a solution. We identified several classes of problems at the individual probe level in the existing annotation, under the assumption that current genome and transcriptome databases are more accurate than those used for GeneChip design. We then reorganized probes on more than a dozen popular GeneChips into gene-, transcript- and exon-specific probe sets in light of up-to-date genome, cDNA/EST clustering and single nucleotide polymorphism information. Comparing analysis results between the original and the redefined probe sets reveals ∼30–50% discrepancy in the genes previously identified as differentially expressed, regardless of analysis method. Our results demonstrate that the original Affymetrix probe set definitions are inaccurate, and many conclusions derived from past GeneChip analyses may be significantly flawed. It will be beneficial to re-analyze existing GeneChip data with updated probe set definitions

CRY2 Is Associated with Depression

Abnormalities in the circadian clockwork often characterize patients with major depressive and bipolar disorders. Circadian clock genes are targets of interest in these patients. CRY2 is a circadian gene that participates in regulation of the evening oscillator. This is of interest in mood disorders where a lack of switch from evening to morning oscillators has been postulated.We observed a marked diurnal variation in human CRY2 mRNA levels from peripheral blood mononuclear cells and a significant up-regulation (P = 0.020) following one-night total sleep deprivation, a known antidepressant. In depressed bipolar patients, levels of CRY2 mRNA were decreased (P = 0.029) and a complete lack of increase was observed following sleep deprivation. To investigate a possible genetic contribution, we undertook SNP genotyping of the CRY2 gene in two independent population-based samples from Sweden (118 cases and 1011 controls) and Finland (86 cases and 1096 controls). The CRY2 gene was significantly associated with winter depression in both samples (haplotype analysis in Swedish and Finnish samples: OR = 1.8, P = 0.0059 and OR = 1.8, P = 0.00044, respectively).We propose that a CRY2 locus is associated with vulnerability for depression, and that mechanisms of action involve dysregulation of CRY2 expression

Mitochondrial Mutations and Polymorphisms in Psychiatric Disorders

Mitochondrial deficiencies with unknown causes have been observed in schizophrenia (SZ) and bipolar disorder (BD) in imaging and postmortem studies. Polymorphisms and somatic mutations in mitochondrial DNA (mtDNA) were investigated as potential causes with next generation sequencing of mtDNA (mtDNA-Seq) and genotyping arrays in subjects with SZ, BD, major depressive disorder (MDD), and controls. The common deletion of 4,977 bp in mtDNA was compared between SZ and controls in 11 different vulnerable brain regions and in blood samples, and in dorsolateral prefrontal cortex (DLPFC) of BD, SZ, and controls. In a separate analysis, association of mitochondria SNPs (mtSNPs) with SZ and BD in European ancestry individuals (n = 6,040) was tested using Genetic Association Information Network (GAIN) and Wellcome Trust Case Control Consortium 2 (WTCCC2) datasets. The common deletion levels were highly variable across brain regions, with a 40-fold increase in some regions (nucleus accumbens, caudate nucleus and amygdala), increased with age, and showed little change in blood samples from the same subjects. The common deletion levels were increased in the DLPFC for BD compared to controls, but not in SZ. Full mtDNA genome resequencing of 23 subjects, showed seven novel homoplasmic mutations, five were novel synonymous coding mutations. By logistic regression analysis there were no significant mtSNPs associated with BD or SZ after genome wide correction. However, nominal association of mtSNPs (p < 0.05) to SZ and BD were found in the hypervariable region of mtDNA to T195C and T16519C. The results confirm prior reports that certain brain regions accumulate somatic mutations at higher levels than blood. The study in mtDNA of common polymorphisms, somatic mutations, and rare mutations in larger populations may lead to a better understanding of the pathophysiology of psychiatric disorders

Sample matching by inferred agonal stress in gene expression analyses of the brain

<p>Abstract</p> <p>Background</p> <p>Gene expression patterns in the brain are strongly influenced by the severity and duration of physiological stress at the time of death. This agonal effect, if not well controlled, can lead to spurious findings and diminished statistical power in case-control comparisons. While some recent studies match samples by tissue pH and clinically recorded agonal conditions, we found that these indicators were sometimes at odds with observed stress-related gene expression patterns, and that matching by these criteria still sometimes results in identifying case-control differences that are primarily driven by residual agonal effects. This problem is analogous to the one encountered in genetic association studies, where self-reported race and ethnicity are often imprecise proxies for an individual's actual genetic ancestry.</p> <p>Results</p> <p>We developed an Agonal Stress Rating (ASR) system that evaluates each sample's degree of stress based on gene expression data, and used ASRs in <it>post hoc </it>sample matching or covariate analysis. While gene expression patterns are generally correlated across different brain regions, we found strong region-region differences in empirical ASRs in many subjects that likely reflect inter-individual variabilities in local structure or function, resulting in region-specific vulnerability to agonal stress.</p> <p>Conclusion</p> <p>Variation of agonal stress across different brain regions differs between individuals, revealing a new level of complexity for gene expression studies of brain tissues. The Agonal Stress Ratings quantitatively assess each sample's extent of regulatory response to agonal stress, and allow a strong control of this important confounder.</p

Conserved chromosome 2q31 conformations are associated with transcriptional regulation of GAD1 GABA synthesis enzyme and altered in prefrontal cortex of subjects with schizophrenia

Little is known about chromosomal loopings involving proximal promoter and distal enhancer elements regulating GABAergic gene expression, including changes in schizophrenia and other psychiatric conditions linked to altered inhibition. Here, we map in human chromosome 2q31 the 3D configuration of 200 kb of linear sequence encompassing the GAD1 GABA synthesis enzyme gene locus, and we describe a loop formation involving the GAD1 transcription start site and intergenic noncoding DNA elements facilitating reporter gene expression. The GAD1-TSS(-50kbLoop) was enriched with nucleosomes epigenetically decorated with the transcriptional mark, histone H3 trimethylated at lysine 4, and was weak or absent in skin fibroblasts and pluripotent stem cells compared with neuronal cultures differentiated from them. In the prefrontal cortex of subjects with schizophrenia, GAD1-TSS(-50kbLoop) was decreased compared with controls, in conjunction with downregulated GAD1 expression. We generated transgenic mice expressing Gad2 promoter-driven green fluorescent protein-conjugated histone H2B and confirmed that Gad1-TSS(-55kbLoop), the murine homolog to GAD1-TSS(-50kbLoop), is a chromosomal conformation specific for GABAergic neurons. In primary neuronal culture, Gad1-TSS(-55kbLoop) and Gad1 expression became upregulated when neuronal activity was increased. We conclude that 3D genome architectures, including chromosomal loopings for promoter-enhancer interactions involved in the regulation of GABAergic gene expression, are conserved between the rodent and primate brain, and subject to developmental and activity-dependent regulation, and disordered in some cases with schizophrenia. More broadly, the findings presented here draw a connection between noncoding DNA, spatial genome architecture, and neuronal plasticity in development and disease

Gene Expression Changes in the Prefrontal Cortex, Anterior Cingulate Cortex and Nucleus Accumbens of Mood Disorders Subjects That Committed Suicide

Suicidal behaviors are frequent in mood disorders patients but only a subset of them ever complete suicide. Understanding predisposing factors for suicidal behaviors in high risk populations is of major importance for the prevention and treatment of suicidal behaviors. The objective of this project was to investigate gene expression changes associated with suicide in brains of mood disorder patients by microarrays (Affymetrix HG-U133 Plus2.0) in the dorsolateral prefrontal cortex (DLPFC: 6 Non-suicides, 15 suicides), the anterior cingulate cortex (ACC: 6NS, 9S) and the nucleus accumbens (NAcc: 8NS, 13S). ANCOVA was used to control for age, gender, pH and RNA degradation, with P≤0.01 and fold change±1.25 as criteria for significance. Pathway analysis revealed serotonergic signaling alterations in the DLPFC and glucocorticoid signaling alterations in the ACC and NAcc. The gene with the lowest p-value in the DLPFC was the 5-HT2A gene, previously associated both with suicide and mood disorders. In the ACC 6 metallothionein genes were down-regulated in suicide (MT1E, MT1F, MT1G, MT1H, MT1X, MT2A) and three were down-regulated in the NAcc (MT1F, MT1G, MT1H). Differential expression of selected genes was confirmed by qPCR, we confirmed the 5-HT2A alterations and the global down-regulation of members of the metallothionein subfamilies MT 1 and 2 in suicide completers. MTs 1 and 2 are neuro-protective following stress and glucocorticoid stimulations, suggesting that in suicide victims neuroprotective response to stress and cortisol may be diminished. Our results thus suggest that suicide-specific expression changes in mood disorders involve both glucocorticoids regulated metallothioneins and serotonergic signaling in different regions of the brain