Gene alteration in zebrafish exposed to a mixture of substances of abuse.

A recent surge in the use and abuse of diverse prescribed psychotic and illicit drugs necessitates the surveillance of drug residues in source water and the associated ecological impacts of chronic exposure to the aquatic organism. Thirty-six psychotic and illicit drug residues were determined in discharged wastewater from two centralized municipal wastewater treatment facilities and two wastewater receiving creeks for seven consecutive days in Kentucky. Zebrafish (Danio rerio) larvae were exposed to the environmental relevant mixtures of all drug residues, all illicit drugs, and all prescribed psychotic drugs. The extracted RNA from fish homogenates was sequenced, and differentially expressed sequences were analyzed for known or predicted nervous system expression, and screened annotated protein-coding genes to the true environmental cocktail mixture. Illicit stimulant (cocaine and one metabolite), opioids (methadone, methadone metabolite, and oxycodone), hallucinogen (MDA), benzodiazepine (oxazepam and temazepam), carbamazepine, and all target selective serotonin reuptake inhibitors including sertraline, fluoxetine, venlafaxine, and citalopram were quantified in 100% of collected samples from both creeks. The high dose cocktail mixture exposure group revealed the largest group of differentially expressed genes: 100 upregulated and 77 downregulated (p ≤ 0.05; q ≤ 0.05). The top 20 differentially expressed sequences in each exposure group comprise 82 unique transcripts corresponding to 74% annotated genes, 7% non-coding sequences, and 19% uncharacterized sequences. Among 61 differentially expressed sequences that corresponded to annotated protein-coding genes, 23 (38%) genes or their homologs are known to be expressed in the nervous system of fish or other organisms. Several of the differentially expressed sequences are associated primarily with the immune system, including several major histocompatibility complex class I and interferon-induced proteins. Interleukin-1 beta (downregulated in this study) abnormalities are considered a risk factor for psychosis. This is the first study to assess the contributions of multiple classes of psychotic and illicit drugs in combination with developmental gene expression.

Genetic vulnerability to DUSP22 promoter hypermethylation is involved in the relation between in utero famine exposure and schizophrenia.

Epigenetic changes may account for the doubled risk to develop schizophrenia in individuals exposed to famine in utero. We therefore investigated DNA methylation in a unique sample of patients and healthy individuals conceived during the great famine in China. Subsequently, we examined two case-control samples without famine exposure in whole blood and brain tissue. To shed light on the causality of the relation between famine exposure and DNA methylation, we exposed human fibroblasts to nutritional deprivation. In the famine-exposed schizophrenia patients, we found significant hypermethylation of the dual specificity phosphatase 22 (DUSP22) gene promoter (Chr6:291687-293285) (N = 153, p = 0.01). In this sample, DUSP22 methylation was also significantly higher in patients independent of famine exposure (p = 0.025), suggesting that hypermethylation of DUSP22 is also more generally involved in schizophrenia risk. Similarly, DUSP22 methylation was also higher in two separate case-control samples not exposed to famine using DNA from whole blood (N = 64, p = 0.03) and postmortem brains (N = 214, p = 0.007). DUSP22 methylation showed strong genetic regulation across chromosomes by a region on chromosome 16 which was consistent with new 3D genome interaction data. The presence of a direct link between famine and DUSP22 transcription was supported by data from cultured human fibroblasts that showed increased methylation (p = 0.048) and expression (p = 0.019) in response to nutritional deprivation (N = 10). These results highlight an epigenetic locus that is genetically regulated across chromosomes and that is involved in the response to early-life exposure to famine and that is relevant for a major psychiatric disorder.

GAD1 alternative transcripts and DNA methylation in human prefrontal cortex and hippocampus in brain development, schizophrenia.

Genetic variations and adverse environmental events in utero or shortly after birth can lead to abnormal brain development and increased risk of schizophrenia. γ-Aminobutyric acid (GABA), the major inhibitory neurotransmitter in the mammalian brain, plays a vital role in normal brain development. GABA synthesis is controlled by enzymes derived from two glutamic acid decarboxylase (GAD) genes, GAD1 and GAD2, both of which produce transcript isoforms. While the full-length GAD1 transcript (GAD67) has been implicated in the neuropathology of schizophrenia, the transcript structure of GAD1 in the human brain has not been fully characterized. In this study, with the use of RNA sequencing and PCR technologies, we report the discovery of 10 novel transcripts of GAD1 in the human brain. Expression levels of four novel GAD1 transcripts (8A, 8B, I80 and I86) showed a lifespan trajectory expression pattern that is anticorrelated with the expression of the full-length GAD1 transcript. In addition, methylation levels of two CpG loci within the putative GAD1 promoter were significantly associated with the schizophrenia-risk SNP rs3749034 and with the expression of GAD25 in dorsolateral prefrontal cortex (DLPFC). Moreover, schizophrenia patients who had completed suicide and/or were positive for nicotine exposure had significantly higher full-length GAD1 expression in the DLPFC. Alternative splicing of GAD1 and epigenetic state appear to play roles in the developmental profile of GAD1 expression and may contribute to GABA dysfunction in the PFC and hippocampus of patients with schizophrenia.

Investigating the neuroimmunogenic architecture of schizophrenia.

The role of the immune system in schizophrenia remains controversial despite numerous studies to date. Most studies have profiled expression of select genes or proteins in peripheral blood, but none have focused on the expression of canonical pathways that mediate overall immune response. The current study used a systematic genetic approach to investigate the role of the immune system in a large sample of post-mortem brain of patients with schizophrenia: RNA sequencing was performed to assess the differential expression of 561 immune genes and 20 immune pathways in dorsolateral prefrontal cortex (DLPFC) (144 schizophrenia and 196 control subjects) and hippocampus (83 schizophrenia and 187 control subjects). The effect of RNA quality on gene expression was found to be highly correlated with the effect of diagnosis even after adjustment for observable RNA quality parameters (i.e. RNA integrity), thus this confounding relationship was statistically controlled using principal components derived from the gene expression matrix. In DLPFC, 23 immune genes were found to be differentially expressed (false discovery rate <0.05), of which seven genes replicated in both directionality and at nominal significance (P<0.05) in an independent post-mortem DLPFC data set (182 schizophrenia and 212 control subjects), although notably at least five of these genes have prominent roles in pathways other than immune function and overall the effect sizes were minimal (fold change <1.1). In the hippocampus, no individual immune genes were identified to be differentially expressed, and in both DLPFC and hippocampus none of the individual immune pathways were relatively differentially expressed. Further, genomic schizophrenia risk profiles scores were not correlated with the expression of individual immune pathways or differentially expressed genes. Overall, past reports claiming a primary pathogenic role of the immune system intrinsic to the brain in schizophrenia could not be confirmed.

Longitudinal analyses of the DNA methylome in deployed military servicemen identify susceptibility loci for post-traumatic stress disorder.

In order to determine the impact of the epigenetic response to traumatic stress on post-traumatic stress disorder (PTSD), this study examined longitudinal changes of genome-wide blood DNA methylation profiles in relation to the development of PTSD symptoms in two prospective military cohorts (one discovery and one replication data set). In the first cohort consisting of male Dutch military servicemen (n=93), the emergence of PTSD symptoms over a deployment period to a combat zone was significantly associated with alterations in DNA methylation levels at 17 genomic positions and 12 genomic regions. Evidence for mediation of the relation between combat trauma and PTSD symptoms by longitudinal changes in DNA methylation was observed at several positions and regions. Bioinformatic analyses of the reported associations identified significant enrichment in several pathways relevant for symptoms of PTSD. Targeted analyses of the significant findings from the discovery sample in an independent prospective cohort of male US marines (n=98) replicated the observed relation between decreases in DNA methylation levels and PTSD symptoms at genomic regions in ZFP57, RNF39 and HIST1H2APS2. Together, our study pinpoints three novel genomic regions where longitudinal decreases in DNA methylation across the period of exposure to combat trauma marks susceptibility for PTSD.

Altered expression of histamine signaling genes in autism spectrum disorder.

The histaminergic system (HS) has a critical role in cognition, sleep and other behaviors. Although not well studied in autism spectrum disorder (ASD), the HS is implicated in many neurological disorders, some of which share comorbidity with ASD, including Tourette syndrome (TS). Preliminary studies suggest that antagonism of histamine receptors 1-3 reduces symptoms and specific behaviors in ASD patients and relevant animal models. In addition, the HS mediates neuroinflammation, which may be heightened in ASD. Together, this suggests that the HS may also be altered in ASD. Using RNA sequencing (RNA-seq), we investigated genome-wide expression, as well as a focused gene set analysis of key HS genes (HDC, HNMT, HRH1, HRH2, HRH3 and HRH4) in postmortem dorsolateral prefrontal cortex (DLPFC) initially in 13 subjects with ASD and 39 matched controls. At the genome level, eight transcripts were differentially expressed (false discovery rate <0.05), six of which were small nucleolar RNAs (snoRNAs). There was no significant diagnosis effect on any of the individual HS genes but expression of the gene set of HNMT, HRH1, HRH2 and HRH3 was significantly altered. Curated HS gene sets were also significantly differentially expressed. Differential expression analysis of these gene sets in an independent RNA-seq ASD data set from DLPFC of 47 additional subjects confirmed these findings. Understanding the physiological relevance of an altered HS may suggest new therapeutic options for the treatment of ASD.

DRD2 co-expression network and a related polygenic index predict imaging, behavioral and clinical phenotypes linked to schizophrenia.

Genetic risk for schizophrenia (SCZ) is determined by many genetic loci whose compound biological effects are difficult to determine. We hypothesized that co-expression pathways of SCZ risk genes are associated with system-level brain function and clinical phenotypes of SCZ. We examined genetic variants related to the dopamine D2 receptor gene DRD2 co-expression pathway and associated them with working memory (WM) behavior, the related brain activity and treatment response. Using two independent post-mortem prefrontal messenger RNA (mRNA) data sets (total N=249), we identified a DRD2 co-expression pathway enriched for SCZ risk genes. Next, we identified non-coding single-nucleotide polymorphisms (SNPs) associated with co-expression of this pathway. These SNPs were associated with regulatory genetic loci in the dorsolateral prefrontal cortex (P<0.05). We summarized their compound effect on co-expression into a Polygenic Co-expression Index (PCI), which predicted DRD2 pathway co-expression in both mRNA data sets (all P<0.05). We associated the PCI with brain activity during WM performance in two independent samples of healthy individuals (total N=368) and 29 patients with SCZ who performed the n-back task. Greater predicted DRD2 pathway prefrontal co-expression was associated with greater prefrontal activity and longer WM reaction times (all corrected P<0.05), thus indicating inefficient WM processing. Blind prediction of treatment response to antipsychotics in two independent samples of patients with SCZ suggested better clinical course of patientswith greater PCI (total N=87; P<0.05). The findings on this DRD2 co-expression pathway are a proof of concept that gene co-expression can parse SCZ risk genes into biological pathways associated with intermediate phenotypes as well as with clinically meaningful information.

Variation in the Williams syndrome GTF2I gene and anxiety proneness interactively affect prefrontal cortical response to aversive stimuli.

Characterizing the molecular mechanisms underlying the heritability of complex behavioral traits such as human anxiety remains a challenging endeavor for behavioral neuroscience. Copy-number variation (CNV) in the general transcription factor gene, GTF2I, located in the 7q11.23 chromosomal region that is hemideleted in Williams syndrome and duplicated in the 7q11.23 duplication syndrome (Dup7), is associated with gene-dose-dependent anxiety in mouse models and in both Williams syndrome and Dup7. Because of this recent preclinical and clinical identification of a genetic influence on anxiety, we examined whether sequence variation in GTF2I, specifically the single-nucleotide polymorphism rs2527367, interacts with trait and state anxiety to collectively impact neural response to anxiety-laden social stimuli. Two hundred and sixty healthy adults completed the Tridimensional Personality Questionnaire Harm Avoidance (HA) subscale, a trait measure of anxiety proneness, and underwent functional magnetic resonance imaging (fMRI) while matching aversive (fearful or angry) facial identity. We found an interaction between GTF2I allelic variations and HA that affects brain response: in individuals homozygous for the major allele, there was no correlation between HA and whole-brain response to aversive cues, whereas in heterozygotes and individuals homozygous for the minor allele, there was a positive correlation between HA sub-scores and a selective dorsolateral prefrontal cortex (DLPFC) responsivity during the processing of aversive stimuli. These results demonstrate that sequence variation in the GTF2I gene influences the relationship between trait anxiety and brain response to aversive social cues in healthy individuals, supporting a role for this neurogenetic mechanism in anxiety.

DEGS2 polymorphism associated with cognition in schizophrenia is associated with gene expression in brain.

A genome-wide association study of cognitive deficits in patients with schizophrenia in Japan found association with a missense genetic variant (rs7157599, Asn8Ser) in the delta(4)-desaturase, sphingolipid 2 (DEGS2) gene. A replication analysis using Caucasian samples showed a directionally consistent trend for cognitive association of a proxy single-nucleotide polymorphism (SNP), rs3783332. Although the DEGS2 gene is expressed in human brain, it is unknown how DEGS2 expression varies during human life and whether it is affected by psychiatric disorders and genetic variants. To address these questions, we examined DEGS2 messenger RNA using next-generation sequencing in postmortem dorsolateral prefrontal cortical tissue from a total of 418 Caucasian samples including patients with schizophrenia, bipolar disorder and major depressive disorder. DEGS2 is expressed at very low levels prenatally and increases gradually from birth to adolescence and consistently expressed across adulthood. Rs3783332 genotype was significantly associated with the expression across all subjects (F3,348=10.79, P=1.12 × 10(-)(3)), particularly in control subjects (F1,87=13.14, P=4.86 × 10(-4)). Similar results were found with rs715799 genotype. The carriers of the risk-associated minor allele at both loci showed significantly lower expression compared with subjects homozygous for the non-risk major allele and this was a consistent finding across all diagnostic groups. DEGS2 expression showed no association with diagnostic status after correcting for multiple testing (P>0.05). Our findings demonstrate that a SNP showing genome-wide association study significant association with cognition in schizophrenia is also associated with regulation of DEGS2 expression, implicating a molecular mechanism for the clinical association.

Genetic neuropathology of obsessive psychiatric syndromes.

Anorexia nervosa (AN), bulimia nervosa (BN) and obsessive-compulsive disorder (OCD) are complex psychiatric disorders with shared obsessive features, thought to arise from the interaction of multiple genes of small effect with environmental factors. Potential candidate genes for AN, BN and OCD have been identified through clinical association and neuroimaging studies; however, recent genome-wide association studies of eating disorders (ED) so far have failed to report significant findings. In addition, few, if any, studies have interrogated postmortem brain tissue for evidence of expression quantitative trait loci (eQTLs) associated with candidate genes, which has particular promise as an approach to elucidating molecular mechanisms of association. We therefore selected single-nucleotide polymorphisms (SNPs) based on candidate gene studies for AN, BN and OCD from the literature, and examined the association of these SNPs with gene expression across the lifespan in prefrontal cortex of a nonpsychiatric control cohort (N=268). Several risk-predisposing SNPs were significantly associated with gene expression among control subjects. We then measured gene expression in the prefrontal cortex of cases previously diagnosed with obsessive psychiatric disorders, for example, ED (N=15) and OCD/obsessive-compulsive personality disorder or tics (OCD/OCPD/Tic; N=16), and nonpsychiatric controls (N=102) and identified 6 and 286 genes that were differentially expressed between ED compared with controls and OCD cases compared with controls, respectively (false discovery rate (FDR) <5%). However, none of the clinical risk SNPs were among the eQTLs and none were significantly associated with gene expression within the broad obsessive cohort, suggesting larger sample sizes or other brain regions may be required to identify candidate molecular mechanisms of clinical association in postmortem brain data sets.

Contrasting changes in DRD1 and DRD2 splice variant expression in schizophrenia and affective disorders, and associations with SNPs in postmortem brain.

Dopamine 2 receptor (DRD2) is of major interest to the pathophysiology of schizophrenia (SCZ) both as a target for antipsychotic drug action as well as a SCZ-associated risk gene. The dopamine 1 receptor (DRD1) is thought to mediate some of the cognitive deficits in SCZ, including impairment of working memory that relies on normal dorsolateral prefrontal cortex (DLPFC) function. To better understand the association of dopamine receptors with SCZ, we studied the expression of three DRD2 splice variants and the DRD1 transcript in DLPFC, hippocampus and caudate nucleus in a large cohort of subjects (~700), including patients with SCZ, affective disorders and nonpsychiatric controls (from 14th gestational week to 85 years of age), and examined genotype-expression associations of 278 single-nucleotide polymorphisms (SNPs) located in or near DRD2 and DRD1 genes. Expression of D2S mRNA and D2S/D2-long (D2L) ratio were significantly increased in DLPFC of patients with SCZ relative to controls (P<0.0001 and P<0.0001, respectively), whereas D2L, D2Longer and DRD1 were decreased (P<0.0001). Patients with affective disorders showed an opposite pattern: reduced expression of D2S (major depressive disorder, P<0.0001) and increased expression of D2L and DRD1 (bipolar disorder, P<0.0001). Moreover, SCZ-associated risk alleles at rs1079727, rs1076560 and rs2283265 predicted increased D2S/D2L expression ratio (P<0.05) in control individuals. Our data suggest that altered splicing of DRD2 and expression of DRD1 may constitute a pathophysiological mechanism in risk for SCZ and affective disorders. The association between SCZ risk-associated polymorphism and the ratio of D2S/D2L is consistent with this possibility.

Revisiting DARPP-32 in postmortem human brain: changes in schizophrenia and bipolar disorder and genetic associations with t-DARPP-32 expression.

Dopamine- and cAMP-regulated phosphoprotein of molecular weight 32 kDa (DARPP-32 or PPP1R1B) has been of interest in schizophrenia owing to its critical function in integrating dopaminergic and glutaminergic signaling. In a previous study, we identified single-nucleotide polymorphisms (SNPs) and a frequent haplotype associated with cognitive and imaging phenotypes that have been linked with schizophrenia, as well as with expression of prefrontal cortical DARPP-32 messenger RNA (mRNA) in a relatively small sample of postmortem brains. In this study, we examined the association of expression of two major DARPP-32 transcripts, full-length (FL-DARPP-32) and truncated (t-DARPP-32), with genetic variants of DARPP-32 in three brain regions receiving dopaminergic input and implicated in schizophrenia (the dorsolateral prefrontal cortex (DLPFC), hippocampus and caudate) in a much larger set of postmortem samples from patients with schizophrenia, bipolar disorder, major depression and normal controls (>700 subjects). We found that the expression of t-DARPP-32 was increased in the DLPFC of patients with schizophrenia and bipolar disorder, and was strongly associated with genotypes at SNPs (rs879606, rs90974 and rs3764352), as well as the previously identified 7-SNP haplotype related to cognitive functioning. The genetic variants that predicted worse cognitive performance were associated with higher t-DARPP-32 expression. Our results suggest that variation in PPP1R1B affects the abundance of the splice variant t-DARPP-32 mRNA and may reflect potential molecular mechanisms implicated in schizophrenia and affective disorders.

Epistatic interaction between COMT and DTNBP1 modulates prefrontal function in mice and in humans.

Cognitive functions are highly heritable and the impact of complex genetic interactions, though undoubtedly important, has received little investigation. Here we show in an animal model and in a human neuroimaging experiment a consistent non-linear interaction between two genes--catechol-O-methyl transferase (COMT) and dysbindin (dys; dystrobrevin-binding protein 1 (DTNBP1))--implicated through different mechanisms in cortical dopamine signaling and prefrontal cognitive function. In mice, we found that a single genetic mutation reducing expression of either COMT or DTNBP1 alone produced working memory advantages, while, in dramatic contrast, genetic reduction of both in the same mouse produced working memory deficits. We found evidence of the same non-linear genetic interaction in prefrontal cortical function in humans. In healthy volunteers (N=176) studied with functional magnetic resonance imaging during a working memory paradigm, individuals homozygous for the COMT rs4680 Met allele that reduces COMT enzyme activity showed a relatively more efficient prefrontal engagement. In contrast, we found that the same genotype was less efficient on the background of a dys haplotype associated with decreased DTNBP1 expression. These results illustrate that epistasis can be functionally multi-directional and non-linear and that a putatively beneficial allele in one epistastic context is a relatively deleterious one in another. These data also have important implications for single-locus association analyses of complex traits.

Brain-derived neurotrophic factor (BDNF) Val(66)Met polymorphism differentially predicts hippocampal function in medication-free patients with schizophrenia.

A Val(66)Met single-nucleotide polymorphism (SNP) in the brain-derived neurotrophic factor (BDNF) gene impairs activity-dependent BDNF release in cultured hippocampal neurons and predicts impaired memory and exaggerated basal hippocampal activity in healthy humans. Several clinical genetic association studies along with multi-modal evidence for hippocampal dysfunction in schizophrenia indirectly suggest a relationship between schizophrenia and genetically determined BDNF function in the hippocampus. To directly test this hypothesized relationship, we studied 47 medication-free patients with schizophrenia or schizoaffective disorder and 74 healthy comparison individuals with genotyping for the Val(66)Met SNP and [(15)O]H(2)O positron emission tomography (PET) to measure resting and working memory-related hippocampal regional cerebral blood flow (rCBF). In patients, harboring a Met allele was associated with significantly less hippocampal rCBF. This finding was opposite to the genotype effect seen in healthy participants, resulting in a significant diagnosis-by-genotype interaction. Exploratory analyses of interregional resting rCBF covariation revealed a specific and significant diagnosis-by-genotype interaction effect on hippocampal-prefrontal coupling. A diagnosis-by-genotype interaction was also found for working memory-related hippocampal rCBF change, which was uniquely attenuated in Met allele-carrying patients. Thus, both task-independent and task-dependent hippocampal neurophysiology accommodates a Met allelic background differently in patients with schizophrenia than in control subjects. Potentially consistent with the hypothesis that cellular sequelae of the BDNF Val(66)Met SNP interface with aspects of schizophrenic hippocampal and frontotemporal dysfunction, these results warrant future investigation to understand the contributions of unique patient trait or state variables to these robust interactions.

Searching for a consensus five-factor model of the Positive and Negative Syndrome Scale for schizophrenia.

Although the developers of the Positive and Negative Syndrome Scale (PANSS) grouped items into three subscales, factor analyses indicate that a five-factor model better characterizes PANSS data. However, lack of consensus on which model to use limits the comparability of PANSS variables across studies. We counted "votes" from published factor analyses to derive consensus models. One of these combined superior fit in our Caucasian sample (n=458, CFI=.970), and in distinct Japanese sample (n=164, CFI=.964), relative to the original three-subscale model, with a sorting of items into factors that was highly consistent across the studies reviewed.