Dopamine perturbation of gene co-expression networks reveals differential response in schizophrenia for translational machinery.

The dopaminergic hypothesis of schizophrenia (SZ) postulates that positive symptoms of SZ, in particular psychosis, are due to disturbed neurotransmission via the dopamine (DA) receptor D2 (DRD2). However, DA is a reactive molecule that yields various oxidative species, and thus has important non-receptor-mediated effects, with empirical evidence of cellular toxicity and neurodegeneration. Here we examine non-receptor-mediated effects of DA on gene co-expression networks and its potential role in SZ pathology. Transcriptomic profiles were measured by RNA-seq in B-cell transformed lymphoblastoid cell lines from 514 SZ cases and 690 controls, both before and after exposure to DA ex vivo (100 µM). Gene co-expression modules were identified using Weighted Gene Co-expression Network Analysis for both baseline and DA-stimulated conditions, with each module characterized for biological function and tested for association with SZ status and SNPs from a genome-wide panel. We identified seven co-expression modules under baseline, of which six were preserved in DA-stimulated data. One module shows significantly increased association with SZ after DA perturbation (baseline: P = 0.023; DA-stimulated: P = 7.8 × 10-5; ΔAIC = -10.5) and is highly enriched for genes related to ribosomal proteins and translation (FDR = 4 × 10-141), mitochondrial oxidative phosphorylation, and neurodegeneration. SNP association testing revealed tentative QTLs underlying module co-expression, notably at FASTKD2 (top P = 2.8 × 10-6), a gene involved in mitochondrial translation. These results substantiate the role of translational machinery in SZ pathogenesis, providing insights into a possible dopaminergic mechanism disrupting mitochondrial function, and demonstrates the utility of disease-relevant functional perturbation in the study of complex genetic etiologies.

Transcriptomic signatures of schizophrenia revealed by dopamine perturbation in an ex vivo model.

The dopaminergic hypothesis of schizophrenia (SZ) postulates that dopaminergic over activity causes psychosis, a central feature of SZ, based on the observation that blocking dopamine (DA) improves psychotic symptoms. DA is known to have both receptor- and non-receptor-mediated effects, including oxidative mechanisms that lead to apoptosis. The role of DA-mediated oxidative processes in SZ has been little studied. Here, we have used a cell perturbation approach and measured transcriptomic profiles by RNAseq to study the effect of DA exposure on transcription in B-cell transformed lymphoblastoid cell lines (LCLs) from 514 SZ cases and 690 controls. We found that DA had widespread effects on both cell growth and gene expression in LCLs. Overall, 1455 genes showed statistically significant differential DA response in SZ cases and controls. This set of differentially expressed genes is enriched for brain expression and for functions related to immune processes and apoptosis, suggesting that DA may play a role in SZ pathogenesis through modulating those systems. Moreover, we observed a non-significant enrichment of genes near genome-wide significant SZ loci and with genes spanned by SZ-associated copy number variants (CNVs), which suggests convergent pathogenic mechanisms detected by both genetic association and gene expression. The study suggests a novel role of DA in the biological processes of immune and apoptosis that may be relevant to SZ pathogenesis. Furthermore, our results show the utility of pathophysiologically relevant perturbation experiments to investigate the biology of complex mental disorders.

Transcriptome outlier analysis implicates schizophrenia susceptibility genes and enriches putatively functional rare genetic variants.

We searched a gene expression dataset comprised of 634 schizophrenia (SZ) cases and 713 controls for expression outliers (i.e., extreme tails of the distribution of transcript expression values) with SZ cases overrepresented compared with controls. These outlier genes were enriched for brain expression and for genes known to be associated with neurodevelopmental disorders. SZ cases showed higher outlier burden (i.e., total outlier events per subject) than controls for genes within copy number variants (CNVs) associated with SZ or neurodevelopmental disorders. Outlier genes were enriched for CNVs and for rare putative regulatory variants, but this only explained a small proportion of the outlier subjects, highlighting the underlying presence of additional genetic and potentially, epigenetic mechanisms.

Common genetic variants and gene expression associated with white matter microstructure in the human brain.

Identifying genes that contribute to white matter microstructure should provide insights into the neurobiological processes that regulate white matter development, plasticity and pathology. We detected five significant SNPs using genome-wide association analysis on a global measure of fractional anisotropy in 776 individuals from large extended pedigrees. Genetic correlations and genome-wide association results indicated that the genetic signal was largely homogeneous across white matter regions. Using RNA transcripts derived from lymphocytes in the same individuals, we identified two genes (GNA13 and CCDC91) that are likely to be cis-regulated by top SNPs, and whose expression levels were also genetically correlated with fractional anisotropy. A transcript of HTR7 was phenotypically associated with FA, and was associated with an intronic genome-wide significant SNP. These results encourage further research in the mechanisms by which GNA13, HTR7 and CCDC91 influence brain structure, and emphasize a role for g-protein signaling in the development and maintenance of white matter microstructure in health and disease.

Transcriptome study of differential expression in schizophrenia.

Schizophrenia genome-wide association studies (GWAS) have identified common SNPs, rare copy number variants (CNVs) and a large polygenic contribution to illness risk, but biological mechanisms remain unclear. Bioinformatic analyses of significantly associated genetic variants point to a large role for regulatory variants. To identify gene expression abnormalities in schizophrenia, we generated whole-genome gene expression profiles using microarrays on lymphoblastoid cell lines (LCLs) from 413 cases and 446 controls. Regression analysis identified 95 transcripts differentially expressed by affection status at a genome-wide false discovery rate (FDR) of 0.05, while simultaneously controlling for confounding effects. These transcripts represented 89 genes with functions such as neurotransmission, gene regulation, cell cycle progression, differentiation, apoptosis, microRNA (miRNA) processing and immunity. This functional diversity is consistent with schizophrenia's likely significant pathophysiological heterogeneity. The overall enrichment of immune-related genes among those differentially expressed by affection status is consistent with hypothesized immune contributions to schizophrenia risk. The observed differential expression of extended major histocompatibility complex (xMHC) region histones (HIST1H2BD, HIST1H2BC, HIST1H2BH, HIST1H2BG and HIST1H4K) converges with the genetic evidence from GWAS, which find the xMHC to be the most significant susceptibility locus. Among the differentially expressed immune-related genes, B3GNT2 is implicated in autoimmune disorders previously tied to schizophrenia risk (rheumatoid arthritis and Graves' disease), and DICER1 is pivotal in miRNA processing potentially linking to miRNA alterations in schizophrenia (e.g. MIR137, the second strongest GWAS finding). Our analysis provides novel candidate genes for further study to assess their potential contribution to schizophrenia.