The relationship between case-control differential gene expression from brain tissue and genetic associations in schizophrenia.

Large numbers of genetic loci have been identified that are known to contain common risk alleles for schizophrenia, but linking associated alleles to specific risk genes remains challenging. Given that most alleles that influence liability to schizophrenia are thought to do so by altered gene expression, intuitively, case-control differential gene expression studies should highlight genes with a higher probability of being associated with schizophrenia and could help identify the most likely causal genes within associated loci. Here, we test this hypothesis by comparing transcriptome analysis of the dorsolateral prefrontal cortex from 563 schizophrenia cases and 802 controls with genome-wide association study (GWAS) data from the third wave study of the Psychiatric Genomics Consortium. Genes differentially expressed in schizophrenia were not enriched for common allelic association statistics compared with other brain-expressed genes, nor were they enriched for genes within associated loci previously reported to be prioritized by genetic fine-mapping. Genes prioritized by Summary-based Mendelian Randomization were underexpressed in cases compared to other genes in the same GWAS loci. However, the overall strength and direction of expression change predicted by SMR were not related to that observed in the differential expression data. Overall, this study does not support the hypothesis that genes identified as differentially expressed from RNA sequencing of bulk brain tissue are enriched for those that show evidence for genetic associations. Such data have limited utility for prioritizing genes in currently associated loci in schizophrenia.

GRK3 deficiency elicits brain immune activation and psychosis.

The G protein-coupled receptor kinase (GRK) family member protein GRK3 has been linked to the pathophysiology of schizophrenia and bipolar disorder. Expression, as well as protein levels, of GRK3 are reduced in post-mortem prefrontal cortex of schizophrenia subjects. Here, we investigate functional behavior and neurotransmission related to immune activation and psychosis using mice lacking functional Grk3 and utilizing a variety of methods, including behavioral, biochemical, electrophysiological, molecular, and imaging methods. Compared to wildtype controls, the Grk3-/- mice show a number of aberrations linked to psychosis, including elevated brain levels of IL-1ß, increased turnover of kynurenic acid (KYNA), hyper-responsiveness to D-amphetamine, elevated spontaneous firing of midbrain dopamine neurons, and disruption in prepulse inhibition. Analyzing human genetic data, we observe a link between psychotic features in bipolar disorder, decreased GRK expression, and increased concentration of CSF KYNA. Taken together, our data suggest that Grk3-/- mice show face and construct validity relating to the psychosis phenotype with glial activation and would be suitable for translational studies of novel immunomodulatory agents in psychotic disorders.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry.

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

Splice-Break: exploiting an RNA-seq splice junction algorithm to discover mitochondrial DNA deletion breakpoints and analyses of psychiatric disorders.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

Melanin Concentrating Hormone Signaling Deficits in Schizophrenia: Association With Memory and Social Impairments and Abnormal Sensorimotor Gating.

BACKGROUND: Evidence from anatomical, pharmacological, and genetic studies supports a role for the neuropeptide melanin concentrating hormone system in modulating emotional and cognitive functions. Genome-wide association studies revealed a potential association between the melanin concentrating hormone receptor (MCHR1) gene locus and schizophrenia, and the largest genome-wide association study conducted to date shows a credible genome-wide association. METHODS: We analyzed MCHR1 and pro-melanin concentrating hormone RNA-Seq expression in the prefrontal cortex in schizophrenia patients and healthy controls. Disruptions in the melanin concentrating hormone system were modeled in the mouse brain by germline deletion of MCHR1 and by conditional ablation of melanin concentrating hormone expressing neurons using a Cre-inducible diphtheria toxin system. RESULTS: MCHR1 expression is decreased in the prefrontal cortex of schizophrenia samples (false discovery rate (FDR) P < .05, CommonMind and PsychEncode combined datasets, n = 901) while pro-melanin concentrating hormone is below the detection threshold. MCHR1 expression decreased with aging (P = 6.6E-57) in human dorsolateral prefrontal cortex. The deletion of MCHR1 was found to lead to behavioral abnormalities mimicking schizophrenia-like phenotypes: hyperactivity, increased stereotypic and repetitive behavior, social impairment, impaired sensorimotor gating, and disrupted cognitive functions. Conditional ablation of pro-melanin concentrating hormone neurons increased repetitive behavior and produced a deficit in sensorimotor gating. CONCLUSIONS: Our study indicates that early disruption of the melanin concentrating hormone system interferes with neurodevelopmental processes, which may contribute to the pathogenesis of schizophrenia. Further neurobiological research on the developmental timing and circuits that are affected by melanin concentrating hormone may lead to a therapeutic target for early prevention of schizophrenia.

Expired Drugs in the Remote Environment.

INTRODUCTION: The British Antarctic Survey Medical Unit works in a very remote area of the world, with several Antarctic bases receiving only a single annual resupply of consumable goods. Pharmaceuticals supplied in this manner will often be approaching or past the end of their nominal shelf life before the following year's resupply. Drugs are transported from the UK via ship; the hold is not temperature controlled, and the ship crosses through the tropics (air temperature 25-30°C for approximately 3 wk). The drugs then must be transported from the ship to the base, often in temperatures substantially below freezing. This study assessed the stability of 5 expired drugs (atropine, nifedipine, flucloxacillin, naproxen, and bendroflumethiazide) returned from Antarctic bases. METHODS: Drugs were opportunistically obtained and tested using stability-indicating assays. RESULTS: All tested drugs were stable. CONCLUSIONS: The results suggest that the studied drugs may be stable beyond expiry, even when not maintained in strictly temperature-controlled conditions.

Circadian patterns of gene expression in the human brain and disruption in major depressive disorder.

A cardinal symptom of major depressive disorder (MDD) is the disruption of circadian patterns. However, to date, there is no direct evidence of circadian clock dysregulation in the brains of patients who have MDD. Circadian rhythmicity of gene expression has been observed in animals and peripheral human tissues, but its presence and variability in the human brain were difficult to characterize. Here, we applied time-of-death analysis to gene expression data from high-quality postmortem brains, examining 24-h cyclic patterns in six cortical and limbic regions of 55 subjects with no history of psychiatric or neurological illnesses ("controls") and 34 patients with MDD. Our dataset covered ~12,000 transcripts in the dorsolateral prefrontal cortex, anterior cingulate cortex, hippocampus, amygdala, nucleus accumbens, and cerebellum. Several hundred transcripts in each region showed 24-h cyclic patterns in controls, and >100 transcripts exhibited consistent rhythmicity and phase synchrony across regions. Among the top-ranked rhythmic genes were the canonical clock genes BMAL1(ARNTL), PER1-2-3, NR1D1(REV-ERBa), DBP, BHLHE40 (DEC1), and BHLHE41(DEC2). The phasing of known circadian genes was consistent with data derived from other diurnal mammals. Cyclic patterns were much weaker in the brains of patients with MDD due to shifted peak timing and potentially disrupted phase relationships between individual circadian genes. This transcriptome-wide analysis of the human brain demonstrates a rhythmic rise and fall of gene expression in regions outside of the suprachiasmatic nucleus in control subjects. The description of its breakdown in MDD suggests potentially important molecular targets for treatment of mood disorders.

Molecular and bioenergetic differences between cells with African versus European inherited mitochondrial DNA haplogroups: implications for population susceptibility to diseases.

The geographic origins of populations can be identified by their maternally inherited mitochondrial DNA (mtDNA) haplogroups. This study compared human cybrids (cytoplasmic hybrids), which are cell lines with identical nuclei but mitochondria from different individuals with mtDNA from either the H haplogroup or L haplogroup backgrounds. The most common European haplogroup is H while individuals of maternal African origin are of the L haplogroup. Despite lower mtDNA copy numbers, L cybrids had higher expression levels for nine mtDNA-encoded respiratory complex genes, decreased ATP (adenosine triphosphate) turnover rates and lower levels of reactive oxygen species production, parameters which are consistent with more efficient oxidative phosphorylation. Surprisingly, GeneChip arrays showed that the L and H cybrids had major differences in expression of genes of the canonical complement system (5 genes), dermatan/chondroitin sulfate biosynthesis (5 genes) and CCR3 (chemokine, CC motif, receptor 3) signaling (9 genes). Quantitative nuclear gene expression studies confirmed that L cybrids had (a) lower expression levels of complement pathway and innate immunity genes and (b) increased levels of inflammation-related signaling genes, which are critical in human diseases. Our data support the hypothesis that mtDNA haplogroups representing populations from different geographic origins may play a role in differential susceptibilities to diseases.

A splice donor mutation in NAA10 results in the dysregulation of the retinoic acid signalling pathway and causes Lenz microphthalmia syndrome.

INTRODUCTION: Lenz microphthalmia syndrome (LMS) is a genetically heterogeneous X-linked disorder characterised by microphthalmia/anophthalmia, skeletal abnormalities, genitourinary malformations, and anomalies of the digits, ears, and teeth. Intellectual disability and seizure disorders are seen in about 60% of affected males. To date, no gene has been identified for LMS in the microphthalmia syndrome 1 locus (MCOPS1). In this study, we aim to find the disease-causing gene for this condition. METHODS AND RESULTS: Using exome sequencing in a family with three affected brothers, we identified a mutation in the intron 7 splice donor site (c.471+2T→A) of the N-acetyltransferase NAA10 gene. NAA10 has been previously shown to be mutated in patients with Ogden syndrome, which is clinically distinct from LMS. Linkage studies for this family mapped the disease locus to Xq27-Xq28, which was consistent with the locus of NAA10. The mutation co-segregated with the phenotype and cDNA analysis showed aberrant transcripts. Patient fibroblasts lacked expression of full length NAA10 protein and displayed cell proliferation defects. Expression array studies showed significant dysregulation of genes associated with genetic forms of anophthalmia such as BMP4, STRA6, and downstream targets of BCOR and the canonical WNT pathway. In particular, STRA6 is a retinol binding protein receptor that mediates cellular uptake of retinol/vitamin A and plays a major role in regulating the retinoic acid signalling pathway. A retinol uptake assay showed that retinol uptake was decreased in patient cells. CONCLUSIONS: We conclude that the NAA10 mutation is the cause of LMS in this family, likely through the dysregulation of the retinoic acid signalling pathway.

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.

Loss of nuclear factor E2-related factor 1 in the brain leads to dysregulation of proteasome gene expression and neurodegeneration.

The ubiquitin-proteasome pathway plays an important role in the pathogenesis of neurodegeneration, but mechanisms controlling expression of components in this pathway remain poorly understood. Nuclear factor E2-related factor 1 (Nrf1) transcription factor has been shown to regulate expression of antioxidant and cytoprotective genes. To determine the function of Nrf1 in the brain, mice with a late-stage deletion of Nrf1 in neuronal cells were generated. Loss of Nrf1 leads to impaired proteasome function and neurodegeneration. Gene expression profiling and RT-PCR analysis revealed a coordinate down-regulation of various proteasomal genes including PsmB6, which encodes a catalytic subunit of the proteasome. Transcriptional analysis and chromatin immunoprecipitation experiments demonstrated that PsmB6 is an Nrf1 target gene. These findings reveal Nrf1 as a key transcriptional regulator required for the expression of proteasomal genes in neurons and suggest that perturbations of Nrf1 function may contribute to the pathogenesis of neurodegenerative diseases.

Increased CNV-region deletions in mild cognitive impairment (MCI) and Alzheimer's disease (AD) subjects in the ADNI sample.

We investigated the genome-wide distribution of CNVs in the Alzheimer's disease (AD) Neuroimaging Initiative (ADNI) sample (146 with AD, 313 with Mild Cognitive Impairment (MCI), and 181 controls). Comparison of single CNVs between cases (MCI and AD) and controls shows overrepresentation of large heterozygous deletions in cases (p-value<0.0001). The analysis of CNV-Regions identifies 44 copy number variable loci of heterozygous deletions, with more CNV-Regions among affected than controls (p=0.005). Seven of the 44 CNV-Regions are nominally significant for association with cognitive impairment. We validated and confirmed our main findings with genome re-sequencing of selected patients and controls. The functional pathway analysis of the genes putatively affected by deletions of CNV-Regions reveals enrichment of genes implicated in axonal guidance, cell-cell adhesion, neuronal morphogenesis and differentiation. Our findings support the role of CNVs in AD, and suggest an association between large deletions and the development of cognitive impairment.

Altered transcriptomes, cell type proportions, and dendritic spine morphology in hippocampus of suicide decedents.

BACKGROUND: Suicide is a manner of death resulting from complex environmental and genetic risks that affect millions of people globally. Both structural and functional studies identified the hippocampus as one of the vulnerable brain regions contributing to suicide risk. METHODS: We have identified the hippocampal tissue transcriptomes, gene ontology, cell type proportions, and dendritic spine morphology in controls (n = 28) and suicide decedents (n = 22). In addition, the transcriptomic signature in iPSC-derived neuronal precursor cells (NPCs) and neurons were also investigated in controls (n = 2) and suicide decedents (n = 2). RESULTS: The hippocampal tissue transcriptomic data revealed that NPAS4 gene expression was downregulated while ALDH1A2, NAAA, and MLXIPL gene expressions were upregulated in hippocampal tissue of suicide decedents. The gene ontology identified 29 significant pathways including NPAS4-associated gene ontology terms "excitatory post-synaptic potential", "regulation of postsynaptic membrane potential" and "long-term memory" indicating alteration of glutamatergic synapses in the hippocampus of suicide decedents. The cell type deconvolution identified decreased excitatory neuron proportion and an increased inhibitory neuron proportion providing evidence of excitation/inhibition imbalance in the hippocampus of suicide decedents. In addition, suicide decedents had increased dendric spine density in the hippocampus, due to an increase of thin (relatively unstable) dendritic spines, compared to controls. The transcriptomes of iPSC-derived hippocampal-like NPCs and neurons revealed 31 and 33 differentially expressed genes in NPC and neurons, respectively, of suicide decedents. CONCLUSIONS: Our findings will provide new insights into the hippocampal neuropathology of suicide.

Common mitochondrial deletions in RNA-Seq: evaluation of bulk, single-cell, and spatial transcriptomic datasets.

Common mitochondrial DNA (mtDNA) deletions are large structural variants in the mitochondrial genome that accumulate in metabolically active tissues with age and have been investigated in various diseases. We applied the Splice-Break2 pipeline (designed for high-throughput quantification of mtDNA deletions) to human RNA-Seq datasets and describe the methodological considerations for evaluating common deletions in bulk, single-cell, and spatial transcriptomics datasets. A robust evaluation of 1570 samples from 14 RNA-Seq studies showed: (i) the abundance of some common deletions detected in PCR-amplified mtDNA correlates with levels observed in RNA-Seq data; (ii) RNA-Seq library preparation method has a strong effect on deletion detection; (iii) deletions had a significant, positive correlation with age in brain and muscle; (iv) deletions were enriched in cortical grey matter, specifically in layers 3 and 5; and (v) brain regions with dopaminergic neurons (i.e., substantia nigra, ventral tegmental area, and caudate nucleus) had remarkable enrichment of common mtDNA deletions.

Antipsychotic drug use complicates assessment of gene expression changes associated with schizophrenia.

Recent postmortem transcriptomic studies of schizophrenia (SCZ) have shown hundreds of differentially expressed genes. However, the extent to which these gene expression changes reflect antipsychotic drug (APD) exposure remains uncertain. We compared differential gene expression in the prefrontal cortex of SCZ patients who tested positive for APDs at the time of death with SCZ patients who did not. APD exposure was associated with numerous changes in the brain transcriptome, especially among SCZ patients on atypical APDs. Brain transcriptome data from macaques chronically treated with APDs showed that APDs affect the expression of many functionally relevant genes, some of which show expression changes in the same directions as those observed in SCZ. Co-expression modules enriched for synaptic function showed convergent patterns between SCZ and some of the APD effects, while those associated with inflammation and glucose metabolism exhibited predominantly divergent patterns between SCZ and APD effects. In contrast, major cell-type shifts inferred in SCZ were primarily unaffected by APD use. These results show that APDs may confound SCZ-associated gene expression changes in postmortem brain tissue. Disentangling these effects will help identify causal genes and improve our neurobiological understanding of SCZ.

Altered transcriptomes, cell type proportions, and dendritic spine morphology in hippocampus of suicide deaths.

Suicide is a condition resulting from complex environmental and genetic risks that affect millions of people globally. Both structural and functional studies identified the hippocampus as one of the vulnerable brain regions contributing to suicide risk. Here, we have identified the hippocampal transcriptomes, gene ontology, cell type proportions, dendritic spine morphology, and transcriptomic signature in iPSC-derived neuronal precursor cells (NPCs) and neurons in postmortem brain tissue from suicide deaths. The hippocampal tissue transcriptomic data revealed that NPAS4 gene expression was downregulated while ALDH1A2, NAAA, and MLXIPL gene expressions were upregulated in tissue from suicide deaths. The gene ontology identified 29 significant pathways including NPAS4-associated gene ontology terms "excitatory post-synaptic potential", "regulation of postsynaptic membrane potential" and "long-term memory" indicating alteration of glutamatergic synapses in the hippocampus of suicide deaths. The cell type deconvolution identified decreased excitatory neuron proportion and an increased inhibitory neuron proportion providing evidence of excitation/inhibition imbalance in the hippocampus of suicide deaths. In addition, suicide deaths had increased dendric spine density, due to an increase of thin (relatively unstable) dendritic spines, compared to controls. The transcriptomes of iPSC-derived hippocampal-like NPCs and neurons revealed 31 and 33 differentially expressed genes in NPC and neurons, respectively, of suicide deaths. The suicide-associated differentially expressed genes in NPCs were RELN, CRH, EMX2, OXTR, PARM1 and IFITM2 which overlapped with previously published results. The previously-known suicide-associated differentially expressed genes in differentiated neurons were COL1A1, THBS1, IFITM2, AQP1, and NLRP2. Together, these findings would help better understand the hippocampal neurobiology of suicide for identifying therapeutic targets to prevent suicide.

Large Common Mitochondrial DNA Deletions Are Associated with a Mitochondrial SNP T14798C Near the 3' Breakpoints.

Introduction: Large somatic deletions of mitochondrial DNA (mtDNA) accumulate with aging in metabolically active tissues such as the brain. We have cataloged the breakpoints and frequencies of large mtDNA deletions in the human brain. Methods: We quantified 112 high-frequency mtDNA somatic deletions across four human brain regions with the Splice-Break2 pipeline. In addition, we utilized PLINK/Seq to test the association of mitochondrial genotypes with the abundance of these high-frequency mtDNA deletions. A conservative p value threshold of 5E-08 was used to find the significant loci. Results: One mtDNA SNP (T14798C) was significantly associated with mtDNA deletions in two brain regions, the dorsolateral prefrontal cortex (DLPFC) and the superior temporal gyrus. Since the DLPFC showed the most robust association between T14798C and two deletion breakpoints (7816-14807 and 5462-14807), this association was tested in the DLPFC of a replication sample and validated the first results. Incorporating the C allele at 14,798 bp increased the perfect/imperfect length of the repeat at the 3' breakpoint of the two associated deletions. Conclusion: This is the first study to identify the association of mtDNA SNP with large mtDNA deletions in the human brain. The T14798C allele located in the MT-CYB gene is a common polymorphism that occurs in several mitochondrial haplogroups. We hypothesize that the T14798C association with two deletions occurs by extending the repeat length around the 3' deletion breakpoints. This simple mechanism suggests that mtDNA SNPs can affect the mitochondrial genome structure, especially in brain where high levels of reactive oxygen species lead to deletion accumulation with aging.

Neurotransmission-related gene expression in the frontal pole is altered in subjects with bipolar disorder and schizophrenia.

The frontal pole (Brodmann area 10, BA10) is the largest cytoarchitectonic region of the human cortex, performing complex integrative functions. BA10 undergoes intensive adolescent grey matter pruning prior to the age of onset for bipolar disorder (BP) and schizophrenia (SCHIZ), and its dysfunction is likely to underly aspects of their shared symptomology. In this study, we investigated the role of BA10 neurotransmission-related gene expression in BP and SCHIZ. We performed qPCR to measure the expression of 115 neurotransmission-related targets in control, BP, and SCHIZ postmortem samples (n = 72). We chose this method for its high sensitivity to detect low-level expression. We then strengthened our findings by performing a meta-analysis of publicly released BA10 microarray data (n = 101) and identified sources of convergence with our qPCR results. To improve interpretation, we leveraged the unusually large database of clinical metadata accompanying our samples to explore the relationship between BA10 gene expression, therapeutics, substances of abuse, and symptom profiles, and validated these findings with publicly available datasets. Using these convergent sources of evidence, we identified 20 neurotransmission-related genes that were differentially expressed in BP and SCHIZ in BA10. These results included a large diagnosis-related decrease in two important therapeutic targets with low levels of expression, HTR2B and DRD4, as well as other findings related to dopaminergic, GABAergic and astrocytic function. We also observed that therapeutics may produce a differential expression that opposes diagnosis effects. In contrast, substances of abuse showed similar effects on BA10 gene expression as BP and SCHIZ, potentially amplifying diagnosis-related dysregulation.

The first decade and beyond of transcriptional profiling in schizophrenia.

Gene expression changes in brains of individuals with schizophrenia (SZ) have been hypothesized to reflect possible pathways related to pathophysiology and/or medication. Other factors having robust effects on gene expression profiling in brain and possibly influence the schizophrenia transcriptome such as age and pH are examined. Pathways of curated gene expression or gene correlation networks reported in SZ (white matter, apoptosis, neurogenesis, synaptic plasticity, glutamatergic and GABAergic neurotransmission, immune and stress-response, mitochondrial, and neurodevelopment) are not unique to SZ and have been associated with other psychiatric disorders. Suggestions going forward to improve the next decade of profiling: consider multiple brain regions that are carefully dissected, release large datasets from multiple brain regions in controls to better understand neurocircuitry, integrate genetics and gene expression, measure expression variants on genome wide level, peripheral biomarker studies, and analyze the transcriptome across a developmental series of brains. Gene expression, while an important feature of the genomic landscape, requires further systems biology to advance from control brains to a more precise definition of the schizophrenia interactome.