Identification of a potential bipolar risk haplotype in the gene encoding the winged-helix transcription factor RFX4

T cell reactivity to regulatory factor X4 in type 1 narcolepsy

Type 1 narcolepsy is strongly (98%) associated with human leukocyte antigen (HLA) class II DQA1*01:02/DQB1*06:02 (DQ0602) and highly associated with T cell receptor (TCR) alpha locus polymorphism as well as other immune regulatory loci. Increased incidence of narcolepsy was detected following the 2009 H1N1 pandemic and linked to Pandemrix vaccination, strongly supporting that narcolepsy is an autoimmune disorder. Although recent results suggest CD4+ T cell reactivity to neuropeptide hypocretin/orexin and cross-reactive flu peptide is involved, identification of other autoantigens has remained elusive. Here we study whether autoimmunity directed against Regulatory Factor X4 (RFX4), a protein co-localized with hypocretin, is involved in some cases of narcolepsy. Studying human serum, we found that autoantibodies against RFX4 were rare. Using RFX4 peptides bound to DQ0602 tetramers, antigen RFX4-86, -95, and -60 specific human CD4+ T cells were detected in 4/10 patients and 2 unaffected siblings, but not in others. Following culture with each cognate peptide, enriched autoreactive TCRαβ clones were isolated by single-cell sorting and TCR sequenced. Homologous clones bearing TRBV4-2 and recognizing RFX4-86 in patients and one twin control of patient were identified. These results suggest the involvement of RFX4 CD4+ T cell autoreactivity in some cases of narcolepsy, but also in healthy donors.

Effect of Ionizing Radiation on Transcriptome during Neural Differentiation of Human Embryonic Stem Cells

Human embryonic brain development is highly sensitive to ionizing radiation. However, detailed information on the mechanisms of this sensitivity is not available due to limited experimental data. In this study, differentiation of human embryonic stem cells (hESCs) to neural lineages was used as a model for early embryonic brain development to assess the effect of exposure to low (17 mGy) and high (572 mGy) doses of radiation on gene expression. Transcriptomes were assessed using RNA sequencing during neural differentiation at three time points in control and irradiated samples. The first time point was when the cells were still pluripotent (day 0), the second time point was during the stage of embryoid body formation (day 6), and the third and final time point was during the stage of neural rosette formation (day 10). Analysis of the transcriptomes revealed neurodifferentiation in both the control and irradiated cells. Low-dose irradiation did not result in changes in gene expression at any of the time points, whereas high-dose irradiation resulted in downregulation of some major neurodifferentiation markers on days 6 and 10. Gene ontology analysis showed that pathways related to nervous system development, neurogenesis and generation of neurons were among the most affected. Expression of such key regulators of neuronal development as NEUROG1, ARX, ASCL1, RFX4 and INSM1 was reduced more than twofold. In conclusion, exposure to a 17 mGy low dose of radiation was well tolerated by hESCs while exposure to 572 mGy significantly affected their genetic reprogramming into neuronal lineages.

Integrative approach for inference of gene regulatory networks using lasso-based random featuring and application to psychiatric disorders

Background

Inferring gene regulatory networks is one of the most interesting research areas in the systems biology. Many inference methods have been developed by using a variety of computational models and approaches. However, there are two issues to solve. First, depending on the structural or computational model of inference method, the results tend to be inconsistent due to innately different advantages and limitations of the methods. Therefore the combination of dissimilar approaches is demanded as an alternative way in order to overcome the limitations of standalone methods through complementary integration. Second, sparse linear regression that is penalized by the regularization parameter (lasso) and bootstrapping-based sparse linear regression methods were suggested in state of the art methods for network inference but they are not effective for a small sample size data and also a true regulator could be missed if the target gene is strongly affected by an indirect regulator with high correlation or another true regulator.

Results

We present two novel network inference methods based on the integration of three different criteria, (i) z-score to measure the variation of gene expression from knockout data, (ii) mutual information for the dependency between two genes, and (iii) linear regression-based feature selection.

Based on these criterion, we propose a lasso-based random feature selection algorithm (LARF) to achieve better performance overcoming the limitations of bootstrapping as mentioned above.

Conclusions

In this work, there are three main contributions. First, our z score-based method to measure gene expression variations from knockout data is more effective than similar criteria of related works. Second, we confirmed that the true regulator selection can be effectively improved by LARF. Lastly, we verified that an integrative approach can clearly outperform a single method when two different methods are effectively jointed. In the experiments, our methods were validated by outperforming the state of the art methods on DREAM challenge data, and then LARF was applied to inferences of gene regulatory network associated with psychiatric disorders.

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.

RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders

Genome-wide expression analysis using next generation sequencing (RNA-Seq) provides an opportunity for in-depth molecular profiling of fundamental biological processes, such as cellular differentiation and malignant transformation. Differentiating human neurons derived from induced pluripotent stem cells (iPSCs) provide an ideal system for RNA-Seq since defective neurogenesis caused by abnormalities in transcription factors, DNA methylation, and chromatin modifiers lie at the heart of some neuropsychiatric disorders. As a preliminary step towards applying next generation sequencing using neurons derived from patient-specific iPSCs, we have carried out an RNA-Seq analysis on control human neurons. Dramatic changes in the expression of coding genes, long non-coding RNAs (lncRNAs), pseudogenes, and splice isoforms were seen during the transition from pluripotent stem cells to early differentiating neurons. A number of genes that undergo radical changes in expression during this transition include candidates for schizophrenia (SZ), bipolar disorder (BD) and autism spectrum disorders (ASD) that function as transcription factors and chromatin modifiers, such as POU3F2 and ZNF804A, and genes coding for cell adhesion proteins implicated in these conditions including NRXN1 and NLGN1. In addition, a number of novel lncRNAs were found to undergo dramatic changes in expression, one of which is HOTAIRM1, a regulator of several HOXA genes during myelopoiesis. The increase we observed in differentiating neurons suggests a role in neurogenesis as well. Finally, several lncRNAs that map near SNPs associated with SZ in genome wide association studies also increase during neuronal differentiation, suggesting that these novel transcripts may be abnormally regulated in a subgroup of patients.

Polymorphisms in ABLIM1 are associated with personality traits and alcohol dependence

Personality traits like novelty seeking (NS), harm avoidance (HA), and reward dependence (RD) are known to be moderately heritable (30-60%). These personality traits and their comorbidities, such as alcohol dependence (AD), may share genetic components. We examined 11,120 single nucleotide polymorphisms (SNPs) genotyped in 292 nuclear families from the Genetic Analysis Workshop 14, a subset from the Collaborative Study on the Genetics of Alcoholism (COGA). A family-based association analysis was performed using the FBAT program. NS, HA, and RD were treated as quantitative traits and AD as a binary trait. Based on a multivariate association test of three quantitative traits in FBAT, we observed 20 SNPs with p < 10(-3). Interestingly, several genes (TESK2, TIPARP, THEMIS, ABLIM1, RFX4, STON2 and LILRA1) are associated with three personality traits with p < 10(-3) using single trait analysis and AD. Especially, SNP rs727532 within ABLIM1 gene at 10q25 showed the most significant association (p = 6.4 × 10(-5)) in the multivariate test and strong associations with NS, HA, RD, and AD (p = 4.48 × 10(-4), 1.2 × 10(-5), 5.6 × 10(-5), 3.12 × 10(-4), respectively) in the COGA sample. In addition, the association of rs727532 with AD was confirmed in a replication study. This study reports some newly recognized associations between several genetic loci and both AD and three personality traits.

Genetic association and sequencing of the insulin-like growth factor 1 gene in bipolar affective disorder

Insulin-like growth factor 1 (IGF1) has been shown to have an important role in brain development and function. Studies of IGF1 administration in rodents have shown that it has an anxiolytic and antidepressant effect. A genome-wide association study (GWAS) of the first University College London (UCL) cohort of 506 bipolar affective disorder subjects and 510 controls was carried out. The exons and flanking regions of IGF1 were resequenced, any new polymorphisms found were genotyped in an enlarged UCL sample of 937 cases and 941 controls. GWAS data gave good evidence of allelic and haplotypic association between multiple IGF1 SNP's and bipolar disorder (BD). New polymorphisms were found by resequencing IGF1 region. Data from GWAS and the new markers showed that twelve out of 43 SNPs showed association with BD with the four most significant SNPs having values of 3.7 × 10(-5) , 8.4 × 10(-4) , 2.6 × 10(-4) , and 2.5 × 10(-4) . A 5' promoter microsatellite polymorphism previously correlated with plasma lipoprotein concentration was also associated with BD (P = 0.013). Haplotypic association confirmed association with BD with significance values similar to the single marker SNP values. The marker rs12426318 has also been found to be associated with BD in a second sample. A test of gene wide significance with permutation testing for all markers genotyped at IGF1 was also significant. These data implicate IGF1 as a candidate gene to cause genetic susceptibility to BD.

Genome-wide association study confirms BST1 and suggests a locus on 12q24 as the risk loci for Parkinson's disease in the European population

We performed a three-stage genome-wide association study (GWAS) to identify common Parkinson's disease (PD) risk variants in the European population. The initial genome-wide scan was conducted in a French sample of 1039 cases and 1984 controls, using almost 500 000 single nucleotide polymorphisms (SNPs). Two SNPs at SNCA were found to be associated with PD at the genome-wide significance level (P < 3 × 10(-8)). An additional set of promising and new association signals was identified and submitted for immediate replication in two independent case-control studies of subjects of European descent. We first carried out an in silico replication study using GWAS data from the WTCCC2 PD study sample (1705 cases, 5200 WTCCC controls). Nominally replicated SNPs were further genotyped in a third sample of 1527 cases and 1864 controls from France and Australia. We found converging evidence of association with PD on 12q24 (rs4964469, combined P = 2.4 × 10(-7)) and confirmed the association on 4p15/BST1 (rs4698412, combined P = 1.8 × 10(-6)), previously reported in Japanese data. The 12q24 locus includes RFX4, an isoform of which, named RFX4_v3, encodes brain-specific transcription factors that regulate many genes involved in brain morphogenesis and intracellular calcium homeostasis.

Proteomics of dense core secretory vesicles reveal distinct protein categories for secretion of neuroeffectors for cell-cell communication

Regulated secretion of neurotransmitters and neurohumoral factors from dense core secretory vesicles provides essential neuroeffectors for cell-cell communication in the nervous and endocrine systems. This study provides comprehensive proteomic characterization of the categories of proteins in chromaffin dense core secretory vesicles that participate in cell-cell communication from the adrenal medulla. Proteomic studies were conducted by nano-HPLC Chip MS/MS tandem mass spectrometry. Results demonstrate that these secretory vesicles contain proteins of distinct functional categories consisting of neuropeptides and neurohumoral factors, protease systems, neurotransmitter enzymes and transporters, receptors, enzymes for biochemical processes, reduction/oxidation regulation, ATPases, protein folding, lipid biochemistry, signal transduction, exocytosis, calcium regulation, as well as structural and cell adhesion proteins. The secretory vesicle proteomic data identified 371 proteins in the soluble fraction and 384 membrane proteins, for a total of 686 distinct secretory vesicle proteins. Notably, these proteomic analyses illustrate the presence of several neurological disease-related proteins in these secretory vesicles, including huntingtin interacting protein, cystatin C, ataxin 7, and prion protein. Overall, these findings demonstrate that multiple protein categories participate in dense core secretory vesicles for production, storage, and secretion of bioactive neuroeffectors for cell-cell communication in health and disease.

Candidate region linkage analysis in twins discordant or concordant for depression symptomatology

Genetic risk factors contribute considerably to both clinical affective disorders and subsyndromal mood level. There is moreover evidence to suggest that the genetic basis of bipolar disorder and unipolar depression overlap to some extent, and several linkage analyses have suggested evidence for a common susceptibility locus in affective disorders on chromosome 12q24. In this study we investigated the chromosome 12 candidate region for linkage to the mean level of depression symptomatology, over a 10-year follow-up, using a highly informative sample of concordant and discordant twin pairs selected from 4,731 participants of the Longitudinal Study of Ageing Danish Twins. Our results showed suggestive evidence of linkage to this region with a peak LOD score of 1.91 for marker D12S1634 located at 148 cM, and thus indicates that the previously identified disease locus at 12q24 is also a general vulnerability locus affecting the normal range of mood.

Case-control studies show that a non-conservative amino-acid change from a glutamine to arginine in the P2RX7 purinergic receptor protein is associated with both bipolar- and unipolar-affective disorders

Three linkage studies of bipolar disorder have implicated chromosome 12q24.3 with lod scores of over 3.0 and several other linkage studies have found lods between 2 and 3. Fine mapping within the original chromosomal linkage regions has identified several loci that show association with bipolar disorder. One of these is the P2RX7 gene encoding a central nervous system-expressed purinergic receptor. A non-synonymous single nucleotide polymorphism, rs2230912 (P2RX7-E13A, G allele) and a microsatellite marker NBG6 were both previously found to be associated with bipolar disorder (P=0.00071 and 0.008, respectively). rs2230912 has also been found to show association with unipolar depression. The effect of the polymorphism is non-conservative and results in a glutamine to arginine change (Gln460Arg), which is likely to affect P2RX7 dimerization and protein-protein interactions. We have confirmed the allelic associations between bipolar disorder and the markers rs2230912 (P2RX7-E13A, G allele, P=0.043) and NBG6 (P=0.010) in a London-based sample of 604 bipolar cases and 560 controls. When we combined these data with the published case-control studies of P2RX7 and mood disorder (3586 individuals) the association between rs2230912 (Gln460Arg) and affective disorders became more robust (P=0.002). The increase in Gln460Arg was confined to heterozygotes rather than homozygotes suggesting a dominant effect (odds ratio 1.302, CI=1.129-1.503). Although further research is needed to prove that the Gln460Arg change has an aetiological role, it is so far the most convincing mutation to have been found with a role for increasing susceptibility to bipolar and genetically related unipolar disorders.

The genetics of bipolar disorder: genome 'hot regions,' genes, new potential candidates and future directions

Bipolar disorder (BP) is a complex disorder caused by a number of liability genes interacting with the environment. In recent years, a large number of linkage and association studies have been conducted producing an extremely large number of findings often not replicated or partially replicated. Further, results from linkage and association studies are not always easily comparable. Unfortunately, at present a comprehensive coverage of available evidence is still lacking. In the present paper, we summarized results obtained from both linkage and association studies in BP. Further, we indicated new potential interesting genes, located in genome 'hot regions' for BP and being expressed in the brain. We reviewed published studies on the subject till December 2007. We precisely localized regions where positive linkage has been found, by the NCBI Map viewer (http://www.ncbi.nlm.nih.gov/mapview/); further, we identified genes located in interesting areas and expressed in the brain, by the Entrez gene, Unigene databases (http://www.ncbi.nlm.nih.gov/entrez/) and Human Protein Reference Database (http://www.hprd.org); these genes could be of interest in future investigations. The review of association studies gave interesting results, as a number of genes seem to be definitively involved in BP, such as SLC6A4, TPH2, DRD4, SLC6A3, DAOA, DTNBP1, NRG1, DISC1 and BDNF. A number of promising genes, which received independent confirmations, and genes that have to be further investigated in BP, have been also systematically listed. In conclusion, the combination of linkage and association approaches provided a number of liability genes. Nevertheless, other approaches are required to disentangle conflicting findings, such as gene interaction analyses, interaction with psychosocial and environmental factors and, finally, endophenotype investigations.

G-protein pathway suppressor 2 (GPS2) interacts with the regulatory factor X4 variant 3 (RFX4_v3) and functions as a transcriptional co-activator

RFX4_v3 (regulatory factor X4 variant 3) is a brain-specific isoform of the transcription factor RFX4. Insertional mutagenesis in mice demonstrates that Rfx4_v3 is crucial for normal brain development. Many genes involved in critical processes during brain morphogenesis are dysregulated in Rfx4_v3 mutant brains. For example, Cx3cl1 is a CX3C-type chemokine that is abundant in brain and is a direct transcriptional target of RFX4_v3 through a specific promoter X-box (X-box 1), the responsive element for RFX proteins. To identify potential interacting partners for RFX4_v3, we performed yeast two-hybrid analysis. Nine candidate interactors were identified, including GPS2 (G-protein pathway suppressor 2). Indirect immunofluorescence demonstrated that GPS2 and RFX4_v3 co-localized to the nucleus. Both GPS2 and RFX4_v3 mRNAs were also present in most portions of the adult mouse brain as well as in brains at different ages, suggesting that the two proteins could bind to each other. Co-immunoprecipitation assays indicated that physical interactions between GPS2 and RFX4_v3 did indeed occur. Furthermore, GPS2 was recruited to the Cx3cl1 promoter by RFX4_v3 and potentiated RFX4_v3 transactivation on this promoter through X-box 1, suggesting that the protein-protein interaction was functionally relevant. GPS2 bound to both the carboxyl-terminal region (amino acids 575-735) and the middle region (amino acids 250-574) of the RFX4_v3 protein. RFX4_v3 amino acids 1-574 stimulated the Cx3cl1 promoter to a similar extent as the full-length RFX4_v3 protein; however, deletion of the carboxyl-terminal region of RFX4_v3 impaired the co-activating abilities of GPS2. Based on these data, we conclude that GPS2 interacts with RFX4_v3 to modulate transactivation of genes involved in brain morphogenesis, including Cx3Cl1.

Regulatory factor X4 variant 3: a transcription factor involved in brain development and disease

Regulatory factor X4 variant 3 (RFX4_v3) is a recently identified transcription factor specifically expressed in the brain. Gene disruption in mice demonstrated that interruption of a single allele (heterozygous, +/-) prevented formation of the subcommissural organ (SCO), resulting in congenital hydrocephalus, whereas interruption of two alleles (homozygous, -/-) caused fatal failure of dorsal midline brain structure formation. These mutagenesis studies implicated RFX4_v3 in early brain development as well as the genesis of the SCO. Rfx4_v3 deficiency presumably causes abnormalities in brain by altering the expression levels of many genes that are crucial for brain morphogenesis, such as the signaling components in the Wnt, bone morphogenetic protein, and retinoic acid pathways. RFX4_v3 might affect these critical signaling pathways in brain development. Cx3cl1, a chemokine gene highly expressed in brain, was identified as a direct target for RFX4_v3, indicating that RFX4_v3 possesses trans-acting activity to stimulate gene expression. Rfx4_v3 is highly expressed in the suprachiasmatic nucleus and might be involved in regulating the circadian clock. One haplotype in RFX4_v3 gene is linked to a higher risk of bipolar disorder, suggesting that this protein might contribute to the pathogenesis of the disease. This Mini-Review describes our current knowledge about RFX4_v3, an important protein that appears to be involved in many aspects of brain development and disease.

Genetic correlates of brain aging on MRI and cognitive test measures: a genome-wide association and linkage analysis in the Framingham Study

BACKGROUND

Brain magnetic resonance imaging (MRI) and cognitive tests can identify heritable endophenotypes associated with an increased risk of developing stroke, dementia and Alzheimer's disease (AD). We conducted a genome-wide association (GWA) and linkage analysis exploring the genetic basis of these endophenotypes in a community-based sample.

METHODS

A total of 705 stroke- and dementia-free Framingham participants (age 62 +9 yrs, 50% male) who underwent volumetric brain MRI and cognitive testing (1999-2002) were genotyped. We used linear models adjusting for first degree relationships via generalized estimating equations (GEE) and family based association tests (FBAT) in additive models to relate qualifying single nucleotide polymorphisms (SNPs, 70,987 autosomal on Affymetrix 100K Human Gene Chip with minor allele frequency > or = 0.10, genotypic call rate > or = 0.80, and Hardy-Weinberg equilibrium p-value > or = 0.001) to multivariable-adjusted residuals of 9 MRI measures including total cerebral brain (TCBV), lobar, ventricular and white matter hyperintensity (WMH) volumes, and 6 cognitive factors/tests assessing verbal and visuospatial memory, visual scanning and motor speed, reading, abstract reasoning and naming. We determined multipoint identity-by-descent utilizing 10,592 informative SNPs and 613 short tandem repeats and used variance component analyses to compute LOD scores.

RESULTS

The strongest gene-phenotype association in FBAT analyses was between SORL1 (rs1131497; p = 3.2 x 10(-6)) and abstract reasoning, and in GEE analyses between CDH4 (rs1970546; p = 3.7 x 10(-8)) and TCBV. SORL1 plays a role in amyloid precursor protein processing and has been associated with the risk of AD. Among the 50 strongest associations (25 each by GEE and FBAT) were other biologically interesting genes. Polymorphisms within 28 of 163 candidate genes for stroke, AD and memory impairment were associated with the endophenotypes studied at p < 0.001. We confirmed our previously reported linkage of WMH on chromosome 4 and describe linkage of reading performance to a marker on chromosome 18 (GATA11A06), previously linked to dyslexia (LOD scores = 2.2 and 5.1).

CONCLUSION

Our results suggest that genes associated with clinical neurological disease also have detectable effects on subclinical phenotypes. These hypothesis generating data illustrate the use of an unbiased approach to discover novel pathways that may be involved in brain aging, and could be used to replicate observations made in other studies.

Analysis of single nucleotide polymorphisms in genes in the chromosome 12Q24.31 region points to P2RX7 as a susceptibility gene to bipolar affective disorder

Previous results from our genetic analyses using pedigrees from a French Canadian population suggested that the interval delimited by markers on chromosome 12, D12S86 and D12S378, was the most probable genomic region to contain a susceptibility gene for affective disorders. Association studies with microsatellite markers using a case/control sample from the same population (n = 427) revealed significant allelic associations between the bipolar phenotype and marker NBG6. Since this marker is located in intron 9 of the P2RX7 gene, we analyzed the surrounding genomic region for the presence of polymorphisms in regulatory, coding and intron/exon junction sequences. Twenty four (24) SNPs were genotyped in a case/control sample and 12 SNPs in all pedigrees used for linkage analysis. Allelic, genotypic or family-based association studies suggest the presence of two susceptibility loci, the P2RX7 and CaMKK2 genes. The strongest association was observed in bipolar families at the non-synonymous SNP P2RX7-E13A (rs2230912, P-value = 0.000708), which results from an over-transmission of the mutant G-allele to affected offspring. This Gln460Arg polymorphism occurs at an amino acid that is conserved between humans and rodents and is located in the C-terminal domain of the P2X7 receptor, known to be essential for normal P2RX7 function.

Genetic Linkage and Association Analysis for Loneliness in Dutch Twin and Sibling Pairs Points to a Region on Chromosome 12q23–24

We obtained evidence from a large study in Dutch twins (N=8,387) and siblings (N=2,295) that variation in loneliness has a genetic component. The heritability estimate for loneliness, which was assessed as an ordinal trait, was 40% and did not differ between males and females. There were 682 sibling pairs with genotypic (around 400 microsatellite markers) data. We combined phenotypic and genotypic data to carry out a genome scan to localize QTLs for loneliness. One region on chromosome 12q23.3-24.3, showed near suggestive linkage. Genetic association tests within this region revealed significant association (p-value 0.009) with one of the alleles of marker D12S79 and with one of the alleles of neighbouring marker D12S395 (p-value 0.043). We review evidence for linkage in this region for psychiatric disorders and discuss our findings within this context.