RNA-sequencing of the brain transcriptome implicates dysregulation of neuroplasticity, circadian rhythms and GTPase binding in bipolar disorder.

RNA-sequencing (RNA-seq) is a powerful technique to investigate the complexity of gene expression in the human brain. We used RNA-seq to survey the brain transcriptome in high-quality postmortem dorsolateral prefrontal cortex from 11 individuals diagnosed with bipolar disorder (BD) and from 11 age- and gender-matched controls. Deep sequencing was performed, with over 350 million reads per specimen. At a false discovery rate of <5%, we detected five differentially expressed (DE) genes and 12 DE transcripts, most of which have not been previously implicated in BD. Among these, Prominin 1/CD133 and ATP-binding cassette-sub-family G-member2 (ABCG2) have important roles in neuroplasticity. We also show for the first time differential expression of long noncoding RNAs (lncRNAs) in BD. DE transcripts include those of serine/arginine-rich splicing factor 5 (SRSF5) and regulatory factor X4 (RFX4), which along with lncRNAs have a role in mammalian circadian rhythms. The DE genes were significantly enriched for several Gene Ontology categories. Of these, genes involved with GTPase binding were also enriched for BD-associated SNPs from previous genome-wide association studies, suggesting that differential expression of these genes is not simply a consequence of BD or its treatment. Many of these findings were replicated by microarray in an independent sample of 60 cases and controls. These results highlight common pathways for inherited and non-inherited influences on disease risk that may constitute good targets for novel therapies.

Copy number variants and therapeutic response to antidepressant medication in major depressive disorder.

It would be beneficial to find genetic predictors of antidepressant response to help personalise treatment of major depressive disorder (MDD). Rare copy number variants (CNVs) have been implicated in several psychiatric disorders, including MDD, but their role in antidepressant response has yet to be investigated. CNV data were available for 1565 individuals with MDD from the NEWMEDS (Novel Methods leading to New Medications in Depression and Schizophrenia) consortium with prospective data on treatment outcome with either a serotonergic or noradrenergic antidepressant. No association was seen between the presence of CNV (rare or common), the overall number of CNVs or genomic CNV 'burden' and antidepressant response. Specific CNVs were nominally associated with antidepressant response, including 15q13.3 duplications and exonic NRXN1 deletions. These were associated with poor response to antidepressants. Overall burden of CNVs is unlikely to contribute to personalising antidepressant treatment. Specific CNVs associated with antidepressant treatment require replication and further study to confirm their role in the therapeutic action of antidepressant.

Meta-analysis of association between obsessive-compulsive disorder and the 3' region of neuronal glutamate transporter gene SLC1A1.

The neuronal glutamate transporter gene SLC1A1 is a candidate gene for obsessive-compulsive disorder (OCD) based on linkage studies and convergent evidence implicating glutamate in OCD etiology. The 3' end of SLC1A1 is the only genomic region with consistently demonstrated OCD association, especially when analyzing male-only probands. However, specific allele associations have not been consistently replicated, and recent OCD genome-wide association and meta-analysis studies have not incorporated all previously associated SLC1A1 SNPs. To clarify the nature of association between SLC1A1 and OCD, pooled analysis was performed on all available relevant raw study data, comprising a final sample of 815 trios, 306 cases and 634 controls. This revealed weak association between OCD and one of nine tested SLC1A1 polymorphisms (rs301443; uncorrected P = 0.046; non-significant corrected P). Secondary analyses of male-affecteds only (N = 358 trios and 133 cases) demonstrated modest association between OCD and a different SNP (rs12682807; uncorrected P = 0.012; non-significant corrected P). Findings of this meta-analysis are consistent with the trend of previous candidate gene studies in psychiatry and do not clarify the putative role of SLC1A1 in OCD pathophysiology. Nonetheless, it may be important to further examine the potential associations demonstrated in this amalgamated sample, especially since the SNPs with modest associations were not included in the more highly powered recent GWAS or in a past meta-analysis including five SLC1A1 polymorphisms. This study underscores the need for much larger sample sizes in future genetic association studies and suggests that next-generation sequencing may be beneficial in examining the potential role of rare variants in OCD.

Genome-wide association study meta-analysis of European and Asian-ancestry samples identifies three novel loci associated with bipolar disorder.

Meta-analyses of bipolar disorder (BD) genome-wide association studies (GWAS) have identified several genome-wide significant signals in European-ancestry samples, but so far account for little of the inherited risk. We performed a meta-analysis of ∼750,000 high-quality genetic markers on a combined sample of ∼14,000 subjects of European and Asian-ancestry (phase I). The most significant findings were further tested in an extended sample of ∼17,700 cases and controls (phase II). The results suggest novel association findings near the genes TRANK1 (LBA1), LMAN2L and PTGFR. In phase I, the most significant single nucleotide polymorphism (SNP), rs9834970 near TRANK1, was significant at the P=2.4 × 10(-11) level, with no heterogeneity. Supportive evidence for prior association findings near ANK3 and a locus on chromosome 3p21.1 was also observed. The phase II results were similar, although the heterogeneity test became significant for several SNPs. On the basis of these results and other established risk loci, we used the method developed by Park et al. to estimate the number, and the effect size distribution, of BD risk loci that could still be found by GWAS methods. We estimate that >63,000 case-control samples would be needed to identify the ∼105 BD risk loci discoverable by GWAS, and that these will together explain <6% of the inherited risk. These results support previous GWAS findings and identify three new candidate genes for BD. Further studies are needed to replicate these findings and may potentially lead to identification of functional variants. Sample size will remain a limiting factor in the discovery of common alleles associated with BD.

Genome-wide association study of obsessive-compulsive disorder.

Obsessive-compulsive disorder (OCD) is a common, debilitating neuropsychiatric illness with complex genetic etiology. The International OCD Foundation Genetics Collaborative (IOCDF-GC) is a multi-national collaboration established to discover the genetic variation predisposing to OCD. A set of individuals affected with DSM-IV OCD, a subset of their parents, and unselected controls, were genotyped with several different Illumina SNP microarrays. After extensive data cleaning, 1465 cases, 5557 ancestry-matched controls and 400 complete trios remained, with a common set of 469,410 autosomal and 9657 X-chromosome single nucleotide polymorphisms (SNPs). Ancestry-stratified case-control association analyses were conducted for three genetically-defined subpopulations and combined in two meta-analyses, with and without the trio-based analysis. In the case-control analysis, the lowest two P-values were located within DLGAP1 (P=2.49 × 10(-6) and P=3.44 × 10(-6)), a member of the neuronal postsynaptic density complex. In the trio analysis, rs6131295, near BTBD3, exceeded the genome-wide significance threshold with a P-value=3.84 × 10(-8). However, when trios were meta-analyzed with the case-control samples, the P-value for this variant was 3.62 × 10(-5), losing genome-wide significance. Although no SNPs were identified to be associated with OCD at a genome-wide significant level in the combined trio-case-control sample, a significant enrichment of methylation QTLs (P<0.001) and frontal lobe expression quantitative trait loci (eQTLs) (P=0.001) was observed within the top-ranked SNPs (P<0.01) from the trio-case-control analysis, suggesting these top signals may have a broad role in gene expression in the brain, and possibly in the etiology of OCD.

Association of SLC6A4 variants with obsessive-compulsive disorder in a large multicenter US family study.

Genetic association studies of SLC6A4 (SERT) and obsessive-compulsive disorder (OCD) have been equivocal. We genotyped 1241 individuals in 278 pedigrees from the OCD Collaborative Genetics Study for 13 single-nucleotide polymorphisms, for the linked polymorphic region (LPR) indel with molecular haplotypes at rs25531, for VNTR polymorphisms in introns 2 and 7 and for a 381-bp deletion 3' to the LPR. We analyzed using the Family-Based Association Test (FBAT) under additive, dominant, recessive and genotypic models, using both OCD and sex-stratified OCD as phenotypes. Two-point FBAT analysis detected association between Int2 (P = 0.0089) and Int7 (P = 0.0187) (genotypic model). Sex-stratified two-point analysis showed strong association in females with Int2 (P<0.0002), significant after correction for linkage disequilibrium, and multiple marker and model testing (P(Adj) = 0.0069). The SLC6A4 gene is composed of two haplotype blocks (our data and the HapMap); FBAT whole-marker analysis conducted using this structure was not significant. Several noteworthy nonsignificant results have emerged. Unlike Hu et al., we found no evidence for overtransmission of the LPR L(A) allele (genotype relative risk = 1.11, 95% confidence interval: 0.77-1.60); however, rare individual haplotypes containing L(A) with P<0.05 were observed. Similarly, three individuals (two with OCD/OCPD) carried the rare I425V SLC6A4 variant, but none of them passed it on to their six OCD-affected offspring, suggesting that it is unlikely to be solely responsible for the 'OCD plus syndrome', as reported by Ozaki et al. In conclusion, we found evidence of genetic association at the SLC6A4 locus with OCD. A noteworthy lack of association at the LPR, LPR-rs25531 and rare 425V variants suggests that hypotheses about OCD risk need revision to accommodate these new findings, including a possible gender effect.

Structural variation of the monoamine oxidase A gene promoter repeat polymorphism in nonhuman primates.

By conferring allele-specific transcriptional activity on the monoamine oxidase A (MAOA) gene in humans, length variation of a repetitive sequence [(variable number of tandem repeat (VNTR)] in the MAOA promoter influences a constellation of personality traits related to aggressive and antisocial behavior and increases the risk of neurodevelopmental and psychiatric disorders. Here, we have analyzed the presence and variability of this MAOA promoter repeat in several species of nonhuman primates. Sequence analysis of MAOA's transcriptional control region revealed the presence of the VNTR in chimpanzee (Pan troglodytes), bonobo (Pan paniscus), gorilla (Gorilla gorilla), orangutan (Pongo pygmaeus), rhesus macaque (Macaca mulatta) and Gelada baboon (Theropithecus gelada). The majority of P. troglodytes and P. paniscus showed a single repeat with a sequence identical to the VNTR sequence in humans. In contrast, analyses of the remaining species revealed shorter sequences similar to the first 18 bp of human VNTR. Compared with other nonhuman primates, the VNTR sequence of M. mulatta showed the highest length variability with allele frequencies of 35, 25 and 40% for the five, six and seven repeat variants, respectively. The extent of variability of the MAOA promoter repeat in both rhesus monkeys and humans supports the notion that there may be a relationship between functional MAOA expression and aggression-related traits in humans and rhesus macaque populations.

Low-threshold heat response antagonized by capsazepine in chick sensory neurons, which are capsaicin-insensitive.

The heat-transducing receptor VR1 cloned from rat sensory neurons can be activated by both noxious heat and capsaicin. As the response of sensory neurons to capsaicin is species dependent, it is conceivable that the responses to noxious heat and to capsaicin are transduced by distinct receptors across different species. Therefore, we investigated responses to noxious heat from a capsaicin-insensitive (chick) and a capsaicin-sensitive (rat) species. In chick, whole-cell patch-clamp experiments in isolated dorsal root ganglion neurons revealed two populations of neurons with different thresholds to noxious heat, activated at approximately 43 degrees C and approximately 53 degrees C. In cobalt uptake experiments, the proportion of neurons showing a heat-induced response increased with increasing heat stimuli. Application of capsaicin (1-10 microM) did not result in inward currents or cobalt uptake. Rat neurons yielded comparable results in heat experiments, but were capsaicin-sensitive. Although chick neurons are insensitive to capsaicin, the competitive capsaicin antagonist capsazepine (1-10 microM) was effective in blocking heat-induced responses, verified by patch-clamp and cobalt uptake methods. The noncompetitive capsaicin antagonist ruthenium red (10 microM) reduced to almost nil the proportion of heat-responsive neurons identified with the cobalt uptake method. These findings suggest that chick DRG neurons express a low-threshold heat-transducing receptor with a pharmacological profile distinct from the low-threshold heat receptor VR1 cloned from rat DRG neurons. The data support the idea that there might be heat receptor subtypes with differences in the capsaicin binding site.