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.

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.

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.