H Le-Niculescu

Convergent functional genomics of schizophrenia: from comprehensive understanding to genetic risk prediction

We have used a translational convergent functional genomics (CFG) approach to identify and prioritize genes involved in schizophrenia, by gene-level integration of genome-wide association study data with other genetic and gene expression studies in humans and animal models. Using this polyevidence scoring and pathway analyses, we identify top genes (DISC1, TCF4, MBP, MOBP, NCAM1, NRCAM, NDUFV2, RAB18, as well as ADCYAP1, BDNF, CNR1, COMT, DRD2, DTNBP1, GAD1, GRIA1, GRIN2B, HTR2A, NRG1, RELN, SNAP-25, TNIK), brain development, myelination, cell adhesion, glutamate receptor signaling, G-protein–coupled receptor signaling and cAMP-mediated signaling as key to pathophysiology and as targets for therapeutic intervention. Overall, the data are consistent with a model of disrupted connectivity in schizophrenia, resulting from the effects of neurodevelopmental environmental stress on a background of genetic vulnerability. In addition, we show how the top candidate genes identified by CFG can be used to generate a genetic risk prediction score (GRPS) to aid schizophrenia diagnostics, with predictive ability in independent cohorts. The GRPS also differentiates classic age of onset schizophrenia from early onset and late-onset disease. We also show, in three independent cohorts, two European American and one African American, increasing overlap, reproducibility and consistency of findings from single-nucleotide polymorphisms to genes, then genes prioritized by CFG, and ultimately at the level of biological pathways and mechanisms. Finally, we compared our top candidate genes for schizophrenia from this analysis with top candidate genes for bipolar disorder and anxiety disorders from previous CFG analyses conducted by us, as well as findings from the fields of autism and Alzheimer. Overall, our work maps the genomic and biological landscape for schizophrenia, providing leads towards a better understanding of illness, diagnostics and therapeutics. It also reveals the significant genetic overlap with other major psychiatric disorder domains, suggesting the need for improved nosology.

Identifying blood biomarkers for mood disorders using convergent functional genomics.

There are to date no objective clinical laboratory blood tests for mood disorders. The current reliance on patient self-report of symptom severity and on the clinicians’ impression is a rate-limiting step in effective treatment and new drug development. We propose, and provide proof of principle for, an approach to help identify blood biomarkers for mood state. We measured whole-genome gene expression differences in blood samples from subjects with bipolar disorder that had low mood vs those that had high mood at the time of the blood draw, and separately, changes in gene expression in brain and blood of a mouse pharmacogenomic model. We then integrated our human blood gene expression data with animal model gene expression data, human genetic linkage/association data and human postmortem brain data, an approach called convergent functional genomics, as a Bayesian strategy for cross-validating and prioritizing findings. Topping our list of candidate blood biomarker genes we have five genes involved in myelination (Mbp, Edg2, Mag, Pmp22 and Ugt8), and six genes involved in growth factor signaling (Fgfr1, Fzd3, Erbb3, Igfbp4, Igfbp6 and Ptprm). All of these genes have prior evidence of differential expression in human postmortem brains from mood disorder subjects. A predictive score developed based on a panel of 10 top candidate biomarkers (five for high mood and five for low mood) shows sensitivity and specificity for high mood and low mood states, in two independent cohorts. Our studies suggest that blood biomarkers may offer an unexpectedly informative window into brain functioning and disease state.

Convergent functional genomics of genome-wide association data for bipolar disorder: Comprehensive identification of candidate genes, pathways and mechanisms

Given the mounting convergent evidence implicating many more genes in complex disorders such as bipolar disorder than the small number identified unambiguously by the first‐generation Genome‐Wide Association studies (GWAS) to date, there is a strong need for improvements in methodology. One strategy is to include in the next generation GWAS larger numbers of subjects, and/or to pool independent studies into meta‐analyses. We propose and provide proof of principle for the use of a complementary approach, convergent functional genomics (CFG), as a way of mining the existing GWAS datasets for signals that are there already, but did not reach significance using a genetics‐only approach. With the CFG approach, the integration of genetics with genomics, of human and animal model data, and of multiple independent lines of evidence converging on the same genes offers a way of extracting signal from noise and prioritizing candidates. In essence our analysis is the most comprehensive integration of genetics and functional genomics to date in the field of bipolar disorder, yielding a series of novel (such as Klf12, Aldh1a1, A2bp1, Ak3l1, Rorb, Rora) and previously known (such as Bdnf, Arntl, Gsk3b, Disc1, Nrg1, Htr2a) candidate genes, blood biomarkers, as well as a comprehensive identification of pathways and mechanisms. These become prime targets for hypothesis driven follow‐up studies, new drug development and personalized medicine approaches.

Identification of blood biomarkers for psychosis using convergent functional genomics.

There are to date no objective clinical laboratory blood tests for psychotic disease states. We provide proof of principle for a convergent functional genomics (CFG) approach to help identify and prioritize blood biomarkers for two key psychotic symptoms, one sensory (hallucinations) and one cognitive (delusions). We used gene expression profiling in whole blood samples from patients with schizophrenia and related disorders, with phenotypic information collected at the time of blood draw, then cross-matched the data with other human and animal model lines of evidence. Topping our list of candidate blood biomarkers for hallucinations, we have four genes decreased in expression in high hallucinations states (Fn1, Rhobtb3, Aldh1l1, Mpp3), and three genes increased in high hallucinations states (Arhgef9, Phlda1, S100a6). All of these genes have prior evidence of differential expression in schizophrenia patients. At the top of our list of candidate blood biomarkers for delusions, we have 15 genes decreased in expression in high delusions states (such as Drd2, Apoe, Scamp1, Fn1, Idh1, Aldh1l1), and 16 genes increased in high delusions states (such as Nrg1, Egr1, Pvalb, Dctn1, Nmt1, Tob2). Twenty-five of these genes have prior evidence of differential expression in schizophrenia patients. Predictive scores, based on panels of top candidate biomarkers, show good sensitivity and negative predictive value for detecting high psychosis states in the original cohort as well as in three additional cohorts. These results have implications for the development of objective laboratory tests to measure illness severity and response to treatment in devastating disorders such as schizophrenia.

Towards understanding the schizophrenia code : An expanded convergent functional genomics approach

Identifying genes for schizophrenia through classical genetic approaches has proven arduous. Here, we present a comprehensive convergent analysis that translationally integrates brain gene expression data from a relevant pharmacogenomic mouse model (involving treatments with a psychomimetic agent—phencyclidine (PCP), and an anti‐psychotic—clozapine), with human genetic linkage data and human postmortem brain data, as a Bayesian strategy of cross validating findings. Topping the list of candidate genes, we have three genes involved in GABA neurotransmission (GABRA1, GABBR1, and GAD2), one gene involved in glutamate neurotransmission (GRIA2), one gene involved in neuropeptide signaling (TAC1), two genes involved in synaptic function (SYN2 and KCNJ4), six genes involved in myelin/glial function (CNP, MAL, MBP, PLP1, MOBP and GFAP), and one gene involved in lipid metabolism (LPL). These data suggest that schizophrenia is primarily a disorder of brain functional and structural connectivity, with GABA neurotransmission playing a prominent role. These findings may explain the EEG gamma band abnormalities detected in schizophrenia. The analysis also revealed other high probability candidates genes (neurotransmitter signaling, other structural proteins, ion channels, signal transduction, regulatory enzymes, neuronal migration/neurite outgrowth, clock genes, transcription factors, RNA regulatory genes), pathways and mechanisms of likely importance in pathophysiology. Some of the pathways identified suggest possible avenues for augmentation pharmacotherapy of schizophrenia with other existing agents, such as benzodiazepines, anticonvulsants and lipid modulating agents. Other pathways are new potential targets for drug development. Lastly, a comparison with our earlier work on bipolar disorder illuminates the significant molecular overlap between schizophrenia and bipolar disorder.

Discovery and validation of blood biomarkers for suicidality.

Suicides are a leading cause of death in psychiatric patients, and in society at large. Developing more quantitative and objective ways (biomarkers) for predicting and tracking suicidal states would have immediate practical applications and positive societal implications. We undertook such an endeavor. First, building on our previous blood biomarker work in mood disorders and psychosis, we decided to identify blood gene expression biomarkers for suicidality, looking at differential expression of genes in the blood of subjects with a major mood disorder (bipolar disorder), a high-risk population prone to suicidality. We compared no suicidal ideation (SI) states and high SI states using a powerful intrasubject design, as well as an intersubject case–case design, to generate a list of differentially expressed genes. Second, we used a comprehensive Convergent Functional Genomics (CFG) approach to identify and prioritize from the list of differentially expressed gene biomarkers of relevance to suicidality. CFG integrates multiple independent lines of evidence—genetic and functional genomic data—as a Bayesian strategy for identifying and prioritizing findings, reducing the false-positives and false-negatives inherent in each individual approach. Third, we examined whether expression levels of the blood biomarkers identified by us in the live bipolar subject cohort are actually altered in the blood in an age-matched cohort of suicide completers collected from the coroner’s office, and report that 13 out of the 41 top CFG scoring biomarkers (32%) show step-wise significant change from no SI to high SI states, and then to the suicide completers group. Six out of them (15%) remained significant after strict Bonferroni correction for multiple comparisons. Fourth, we show that the blood levels of SAT1 (spermidine/spermine N1–acetyltransferase 1), the top biomarker identified by us, at the time of testing for this study, differentiated future as well as past hospitalizations with suicidality, in a live cohort of bipolar disorder subjects, and exhibited a similar but weaker pattern in a live cohort of psychosis (schizophrenia/schizoaffective disorder) subjects. Three other (phosphatase and tensin homolog (PTEN), myristoylated alanine-rich protein kinase C substrate (MARCKS), and mitogen-activated protein kinase kinase kinase 3 (MAP3K3)) of the six biomarkers that survived Bonferroni correction showed similar but weaker effects. Taken together, the prospective and retrospective hospitalization data suggests SAT1, PTEN, MARCKS and MAP3K3 might be not only state biomarkers but trait biomarkers as well. Fifth, we show how a multi-dimensional approach using SAT1 blood expression levels and two simple visual-analog scales for anxiety and mood enhances predictions of future hospitalizations for suicidality in the bipolar cohort (receiver-operating characteristic curve with area under the curve of 0.813). Of note, this simple approach does not directly ask about SI, which some individuals may deny or choose not to share with clinicians. Lastly, we conducted bioinformatic analyses to identify biological pathways, mechanisms and medication targets. Overall, suicidality may be underlined, at least in part, by biological mechanisms related to stress, inflammation and apoptosis.

Phenomic, convergent functional genomic, and biomarker studies in a stress-reactive genetic animal model of bipolar disorder and co-morbid alcoholism.

We had previously identified the clock gene D‐box binding protein (Dbp) as a potential candidate gene for bipolar disorder and for alcoholism, using a Convergent Functional Genomics (CFG) approach. Here we report that mice with a homozygous deletion of DBP have lower locomotor activity, blunted responses to stimulants, and gain less weight over time. In response to a chronic stress paradigm, these mice exhibit a diametric switch in these phenotypes. DBP knockout mice are also activated by sleep deprivation, similar to bipolar patients, and that activation is prevented by treatment with the mood stabilizer drug valproate. Moreover, these mice show increased alcohol intake following exposure to stress. Microarray studies of brain and blood reveal a pattern of gene expression changes that may explain the observed phenotypes. CFG analysis of the gene expression changes identified a series of novel candidate genes and blood biomarkers for bipolar disorder, alcoholism, and stress reactivity.

Candidate genes, pathways and mechanisms for alcoholism: an expanded convergent functional genomics approach

We describe a comprehensive translational approach for identifying candidate genes for alcoholism. The approach relies on the cross-matching of animal model brain gene expression data with human genetic linkage data, as well as human tissue data and biological roles data, an approach termed convergent functional genomics. An analysis of three animal model paradigms, based on inbred alcohol-preferring (iP) and alcohol-non-preferring (iNP) rats, and their response to treatments with alcohol, was used. A comprehensive analysis of microarray gene expression data from five key brain regions (frontal cortex, amygdala, caudate–putamen, nucleus accumbens and hippocampus) was carried out. The Bayesian-like integration of multiple independent lines of evidence, each by itself lacking sufficient discriminatory power, led to the identification of high probability candidate genes, pathways and mechanisms for alcoholism. These data reveal that alcohol has pleiotropic effects on multiple systems, which may explain the diverse neuropsychiatric and medical pathology in alcoholism. Some of the pathways identified suggest avenues for pharmacotherapy of alcoholism with existing agents, such as angiotensin-converting enzyme (ACE) inhibitors. Experiments we carried out in alcohol-preferring rats with an ACE inhibitor show a marked modulation of alcohol intake. Other pathways are new potential targets for drug development. The emergent overall picture is that physical and physiological robustness may permit alcohol-preferring individuals to withstand the aversive effects of alcohol. In conjunction with a higher reactivity to its rewarding effects, they may able to ingest enough of this nonspecific drug for a strong hedonic and addictive effect to occur.

Evidence for genetic association of RORB with bipolar disorder

BackgroundBipolar disorder, particularly in children, is characterized by rapid cycling and switching, making circadian clock genes plausible molecular underpinnings for bipolar disorder. We previously reported work establishing mice lacking the clock gene D-box binding protein (DBP) as a stress-reactive genetic animal model of bipolar disorder. Microarray studies revealed that expression of two closely related clock genes, RAR-related orphan receptors alpha (RORA) and beta (RORB), was altered in these mice. These retinoid-related receptors are involved in a number of pathways including neurogenesis, stress response, and modulation of circadian rhythms. Here we report association studies between bipolar disorder and single-nucleotide polymorphisms (SNPs) in RORA and RORB.MethodsWe genotyped 355 RORA and RORB SNPs in a pediatric cohort consisting of a family-based sample of 153 trios and an independent, non-overlapping case-control sample of 152 cases and 140 controls. Bipolar disorder in children and adolescents is characterized by increased stress reactivity and frequent episodes of shorter duration; thus our cohort provides a potentially enriched sample for identifying genes involved in cycling and switching.ResultsWe report that four intronic RORB SNPs showed positive associations with the pediatric bipolar phenotype that survived Bonferroni correction for multiple comparisons in the case-control sample. Three RORB haplotype blocks implicating an additional 11 SNPs were also associated with the disease in the case-control sample. However, these significant associations were not replicated in the sample of trios. There was no evidence for association between pediatric bipolar disorder and any RORA SNPs or haplotype blocks after multiple-test correction. In addition, we found no strong evidence for association between the age-at-onset of bipolar disorder with any RORA or RORB SNPs.ConclusionOur findings suggest that clock genes in general and RORB in particular may be important candidates for further investigation in the search for the molecular basis of bipolar disorder.

Convergent functional genomics of anxiety disorders: translational identification of genes, biomarkers, pathways and mechanisms

Anxiety disorders are prevalent and disabling yet understudied from a genetic standpoint, compared with other major psychiatric disorders such as bipolar disorder and schizophrenia. The fact that they are more common, diverse and perceived as embedded in normal life may explain this relative oversight. In addition, as for other psychiatric disorders, there are technical challenges related to the identification and validation of candidate genes and peripheral biomarkers. Human studies, particularly genetic ones, are susceptible to the issue of being underpowered, because of genetic heterogeneity, the effect of variable environmental exposure on gene expression, and difficulty of accrual of large, well phenotyped cohorts. Animal model gene expression studies, in a genetically homogeneous and experimentally tractable setting, can avoid artifacts and provide sensitivity of detection. Subsequent translational integration of the animal model datasets with human genetic and gene expression datasets can ensure cross-validatory power and specificity for illness. We have used a pharmacogenomic mouse model (involving treatments with an anxiogenic drug—yohimbine, and an anti-anxiety drug—diazepam) as a discovery engine for identification of anxiety candidate genes as well as potential blood biomarkers. Gene expression changes in key brain regions for anxiety (prefrontal cortex, amygdala and hippocampus) and blood were analyzed using a convergent functional genomics (CFG) approach, which integrates our new data with published human and animal model data, as a translational strategy of cross-matching and prioritizing findings. Our work identifies top candidate genes (such as FOS, GABBR1, NR4A2, DRD1, ADORA2A, QKI, RGS2, PTGDS, HSPA1B, DYNLL2, CCKBR and DBP), brain–blood biomarkers (such as FOS, QKI and HSPA1B), pathways (such as cAMP signaling) and mechanisms for anxiety disorders—notably signal transduction and reactivity to environment, with a prominent role for the hippocampus. Overall, this work complements our previous similar work (on bipolar mood disorders and schizophrenia) conducted over the last decade. It concludes our programmatic first pass mapping of the genomic landscape of the triad of major psychiatric disorder domains using CFG, and permitted us to uncover the significant genetic overlap between anxiety and these other major psychiatric disorders, notably the under-appreciated overlap with schizophrenia. PDE10A, TAC1 and other genes uncovered by our work provide a molecular basis for the frequently observed clinical co-morbidity and interdependence between anxiety and other major psychiatric disorders, and suggest schizo-anxiety as a possible new nosological domain.Translational Psychiatry (2011) 1, e9; doi:10.1038/tp.2011.9; published online 24 May 2011