Association of prescription data with clinical manifestations and polygenic risk scores in patients with bipolar I disorder: An exploratory study.

BACKGROUND: We assessed the association of prescription data with clinical manifestations and polygenic risk scores (PRS) in patients with bipolar I disorder. METHODS: We enrolled 1471 individuals with BID and divided them into several groups according to treatment options and clinical manifestations. BD-PRS of each patient was calculated using the Psychiatric Genomics Consortium data. Data on single nucleotide polymorphisms, clinical manifestations, and prescriptions were extracted from BiGS. RESULTS: Chronicity, suicidality, substance misuse, mixed symptoms, and deterioration of life functioning were significantly more severe in the group that was not prescribed any mood stabilizers (MS). Chronicity, psychotic symptoms, suicidality, and deterioration of life functioning were significantly severe in the group that received two or more antipsychotics (APs). BD-PRS between the group with AP(s) only and that with other treatment options significantly differed. BD-PRS of the group with AP(s) only was significantly lower than that with other treatment options. Our linear regression results showed that high severity of particular clinical aspects, lower BD-PRS, and prescriptions with fewer MSs or more APs were independently associated with poor life functioning. LIMITATIONS: This study had a cross-sectional design, without differentiating the bipolar phase, which could influence our results. CONCLUSIONS: Poor life functioning in patients with BID was associated with a high severity of particular clinical aspects, BD-PRS, and prescriptions including fewer MSs or more APs. BD-PRS was significantly higher in the group receiving only AP(s) than that in the groups receiving other drugs.

Association of polygenic score and the involvement of cholinergic and glutamatergic pathways with lithium treatment response in patients with bipolar disorder.

Lithium is regarded as the first-line treatment for bipolar disorder (BD), a severe and disabling mental health disorder that affects about 1% of the population worldwide. Nevertheless, lithium is not consistently effective, with only 30% of patients showing a favorable response to treatment. To provide personalized treatment options for bipolar patients, it is essential to identify prediction biomarkers such as polygenic scores. In this study, we developed a polygenic score for lithium treatment response (Li+PGS) in patients with BD. To gain further insights into lithium's possible molecular mechanism of action, we performed a genome-wide gene-based analysis. Using polygenic score modeling, via methods incorporating Bayesian regression and continuous shrinkage priors, Li+PGS was developed in the International Consortium of Lithium Genetics cohort (ConLi+Gen: N = 2367) and replicated in the combined PsyCourse (N = 89) and BipoLife (N = 102) studies. The associations of Li+PGS and lithium treatment response - defined in a continuous ALDA scale and a categorical outcome (good response vs. poor response) were tested using regression models, each adjusted for the covariates: age, sex, and the first four genetic principal components. Statistical significance was determined at P < 0.05. Li+PGS was positively associated with lithium treatment response in the ConLi+Gen cohort, in both the categorical (P = 9.8 × 10-12, R2 = 1.9%) and continuous (P = 6.4 × 10-9, R2 = 2.6%) outcomes. Compared to bipolar patients in the 1st decile of the risk distribution, individuals in the 10th decile had 3.47-fold (95%CI: 2.22-5.47) higher odds of responding favorably to lithium. The results were replicated in the independent cohorts for the categorical treatment outcome (P = 3.9 × 10-4, R2 = 0.9%), but not for the continuous outcome (P = 0.13). Gene-based analyses revealed 36 candidate genes that are enriched in biological pathways controlled by glutamate and acetylcholine. Li+PGS may be useful in the development of pharmacogenomic testing strategies by enabling a classification of bipolar patients according to their response to treatment.

Focal adhesion is associated with lithium response in bipolar disorder: evidence from a network-based multi-omics analysis.

Lithium (Li) is one of the most effective drugs for treating bipolar disorder (BD), however, there is presently no way to predict response to guide treatment. The aim of this study is to identify functional genes and pathways that distinguish BD Li responders (LR) from BD Li non-responders (NR). An initial Pharmacogenomics of Bipolar Disorder study (PGBD) GWAS of lithium response did not provide any significant results. As a result, we then employed network-based integrative analysis of transcriptomic and genomic data. In transcriptomic study of iPSC-derived neurons, 41 significantly differentially expressed (DE) genes were identified in LR vs NR regardless of lithium exposure. In the PGBD, post-GWAS gene prioritization using the GWA-boosting (GWAB) approach identified 1119 candidate genes. Following DE-derived network propagation, there was a highly significant overlap of genes between the top 500- and top 2000-proximal gene networks and the GWAB gene list (Phypergeometric = 1.28E-09 and 4.10E-18, respectively). Functional enrichment analyses of the top 500 proximal network genes identified focal adhesion and the extracellular matrix (ECM) as the most significant functions. Our findings suggest that the difference between LR and NR was a much greater effect than that of lithium. The direct impact of dysregulation of focal adhesion on axon guidance and neuronal circuits could underpin mechanisms of response to lithium, as well as underlying BD. It also highlights the power of integrative multi-omics analysis of transcriptomic and genomic profiling to gain molecular insights into lithium response in BD.

Association of Polygenic Score and the involvement of Cholinergic and Glutamatergic Pathways with Lithium Treatment Response in Patients with Bipolar Disorder.

Lithium is regarded as the first-line treatment for bipolar disorder (BD), a severe and disabling mental disorder that affects about 1% of the population worldwide. Nevertheless, lithium is not consistently effective, with only 30% of patients showing a favorable response to treatment. To provide personalized treatment options for bipolar patients, it is essential to identify prediction biomarkers such as polygenic scores. In this study, we developed a polygenic score for lithium treatment response (Li+PGS) in patients with BD. To gain further insights into lithium's possible molecular mechanism of action, we performed a genome-wide gene-based analysis. Using polygenic score modeling, via methods incorporating Bayesian regression and continuous shrinkage priors, Li+PGS was developed in the International Consortium of Lithium Genetics cohort (ConLi+Gen: N=2,367) and replicated in the combined PsyCourse (N=89) and BipoLife (N=102) studies. The associations of Li+PGS and lithium treatment response - defined in a continuous ALDA scale and a categorical outcome (good response vs. poor response) were tested using regression models, each adjusted for the covariates: age, sex, and the first four genetic principal components. Statistical significance was determined at P<����������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������.

Rare variants implicate NMDA receptor signaling and cerebellar gene networks in risk for bipolar disorder.

Bipolar disorder is an often-severe mental health condition characterized by alternation between extreme mood states of mania and depression. Despite strong heritability and the recent identification of 64 common variant risk loci of small effect, pathophysiological mechanisms remain unknown. Here, we analyzed genome sequences from 41 multiply-affected pedigrees and identified variants in 741 genes with nominally significant linkage or association with bipolar disorder. These 741 genes overlapped known risk genes for neurodevelopmental disorders and clustered within gene networks enriched for synaptic and nuclear functions. The top variant in this analysis - prioritized by statistical association, predicted deleteriousness, and network centrality - was a missense variant in the gene encoding D-amino acid oxidase (DAOG131V). Heterologous expression of DAOG131V in human cells resulted in decreased DAO protein abundance and enzymatic activity. In a knock-in mouse model of DAOG131, DaoG130V/+, we similarly found decreased DAO protein abundance in hindbrain regions, as well as enhanced stress susceptibility and blunted behavioral responses to pharmacological inhibition of N-methyl-D-aspartate receptors (NMDARs). RNA sequencing of cerebellar tissue revealed that DaoG130V resulted in decreased expression of two gene networks that are enriched for synaptic functions and for genes expressed, respectively, in Purkinje neurons or granule neurons. These gene networks were also down-regulated in the cerebellum of patients with bipolar disorder compared to healthy controls and were enriched for additional rare variants associated with bipolar disorder risk. These findings implicate dysregulation of NMDAR signaling and of gene expression in cerebellar neurons in bipolar disorder pathophysiology and provide insight into its genetic architecture.

HLA-DRB1 and HLA-DQB1 genetic diversity modulates response to lithium in bipolar affective disorders.

Bipolar affective disorder (BD) is a severe psychiatric illness, for which lithium (Li) is the gold standard for acute and maintenance therapies. The therapeutic response to Li in BD is heterogeneous and reliable biomarkers allowing patients stratification are still needed. A GWAS performed by the International Consortium on Lithium Genetics (ConLiGen) has recently identified genetic markers associated with treatment responses to Li in the human leukocyte antigens (HLA) region. To better understand the molecular mechanisms underlying this association, we have genetically imputed the classical alleles of the HLA region in the European patients of the ConLiGen cohort. We found our best signal for amino-acid variants belonging to the HLA-DRB1*11:01 classical allele, associated with a better response to Li (p < 1 × 10-3; FDR < 0.09 in the recessive model). Alanine or Leucine at position 74 of the HLA-DRB1 heavy chain was associated with a good response while Arginine or Glutamic acid with a poor response. As these variants have been implicated in common inflammatory/autoimmune processes, our findings strongly suggest that HLA-mediated low inflammatory background may contribute to the efficient response to Li in BD patients, while an inflammatory status overriding Li anti-inflammatory properties would favor a weak response.

Association of polygenic score for major depression with response to lithium in patients with bipolar disorder.

Lithium is a first-line medication for bipolar disorder (BD), but only one in three patients respond optimally to the drug. Since evidence shows a strong clinical and genetic overlap between depression and bipolar disorder, we investigated whether a polygenic susceptibility to major depression is associated with response to lithium treatment in patients with BD. Weighted polygenic scores (PGSs) were computed for major depression (MD) at different GWAS p value thresholds using genetic data obtained from 2586 bipolar patients who received lithium treatment and took part in the Consortium on Lithium Genetics (ConLi+Gen) study. Summary statistics from genome-wide association studies in MD (135,458 cases and 344,901 controls) from the Psychiatric Genomics Consortium (PGC) were used for PGS weighting. Response to lithium treatment was defined by continuous scores and categorical outcome (responders versus non-responders) using measurements on the Alda scale. Associations between PGSs of MD and lithium treatment response were assessed using a linear and binary logistic regression modeling for the continuous and categorical outcomes, respectively. The analysis was performed for the entire cohort, and for European and Asian sub-samples. The PGSs for MD were significantly associated with lithium treatment response in multi-ethnic, European or Asian populations, at various p value thresholds. Bipolar patients with a low polygenic load for MD were more likely to respond well to lithium, compared to those patients with high polygenic load [lowest vs highest PGS quartiles, multi-ethnic sample: OR = 1.54 (95% CI: 1.18-2.01) and European sample: OR = 1.75 (95% CI: 1.30-2.36)]. While our analysis in the Asian sample found equivalent effect size in the same direction: OR = 1.71 (95% CI: 0.61-4.90), this was not statistically significant. Using PGS decile comparison, we found a similar trend of association between a high genetic loading for MD and lower response to lithium. Our findings underscore the genetic contribution to lithium response in BD and support the emerging concept of a lithium-responsive biotype in BD.

Ntrk1 mutation co-segregating with bipolar disorder and inherited kidney disease in a multiplex family causes defects in neuronal growth and depression-like behavior in mice.

Previously, we reported a family in which bipolar disorder (BD) co-segregates with a Mendelian kidney disorder linked to 1q22. The causative renal gene was later identified as MUC1. Genome-wide linkage analysis of BD in the family yielded a peak at 1q22 that encompassed the NTRK1 and MUC1 genes. NTRK1 codes for TrkA (Tropomyosin-related kinase A) which is essential for development of the cholinergic nervous system. Whole genome sequencing of the proband identified a damaging missense mutation, E492K, in NTRK1. Induced pluripotent stem cells were generated from family members, and then differentiated to neural stem cells (NSCs). E492K NSCs had reduced neurite outgrowth. A conditional knock-in mouse line, harboring the point mutation in the brain, showed depression-like behavior in the tail suspension test following challenge by physostigmine, a cholinesterase inhibitor. These results are consistent with the cholinergic hypothesis of depression. They imply that the NTRK1 E492K mutation, impairs cholinergic neurotransmission, and may convey susceptibility to bipolar disorder.

Interaction between adverse childhood experiences and polygenic risk in patients with bipolar disorder.

The interaction between genes and environment often occurs when they depend on one another. We hypothesized that adverse childhood experiences (ACEs) would interact with genetic predispositions to bipolar disorder (BD), demonstrating earlier age at onset (AAO) and worse clinical outcomes. We aimed to clarify the effects of the interaction between ACEs and genetic susceptibility using polygenic risk score (PRS) on AAO and clinical outcomes. Single nucleotide polymorphisms and clinical data, including ACEs, were obtained from the Bipolar Genomic Study, which contains a large sample of BD participants. A total of 1615 subjects with BD I were obtained and divided into two groups according to the presence or absence of ACEs and an additional four groups based on the number of ACEs (none versus one versus two versus ≥ three types). ACEs was evaluated using the childhood life events scale (CLES). BD-PRS was obtained from the Psychiatric Genomics Consortium, which compared BD patients and healthy controls. The BD-PRS was higher in the group with ACEs than without ACEs at most p-value thresholds. In multivariate linear regression analyses, both groups with more ACEs and higher BD-PRS were independently and interactively associated with an earlier AAO of BD; however, only greater ACEs were associated with worsened clinical outcome. These findings highlight the clinical importance of evaluating ACEs and polygenic risk in research of the etiology of BD.

Effect of the Type and Number of Adverse Childhood Experiences and the Timing of Adverse Experiences on Clinical Outcomes in Individuals with Bipolar Disorder.

: Studies have reported an association between adverse childhood experiences (ACEs) and the clinical outcomes of bipolar disorder (BD). However, these studies have several limitations; therefore, we aimed to clarify the effect of the type and number of ACEs and the timing of adverse experiences on clinical outcomes in patients with BD. We analyzed the data of patients with BD (N = 2675) obtained from the National Institute of Mental Health: Bipolar Disorder Genetic Association Information Network, Translational Genomic Institute-I, and Translational Genomic Institute-II. All patients had been diagnosed using the Diagnostic Interview for Genetic Studies. ACEs were evaluated using the Childhood Life Events Scale (CLES). We analyzed the relationship between childhood trauma and clinical outcome in patients with and without exposure to ACEs. We found that ACEs had a robust negative effect on clinical outcomes, including earlier age at onset, presence of psychotic episodes, suicide attempts, mixed symptoms or episodes, substance misuse comorbidity, and worse life functioning. Specifically, the number of ACEs had the most significant effect on clinical outcomes; however, specific ACEs, such as physical abuse, had a considerable influence. Moreover, post-childhood adverse experiences had a weaker effect on clinical outcomes than ACEs did. There was an association of ACEs with negative clinical outcomes in patients with BD. This indicates the importance of basic and clinical research on ACEs in patients with BD.

A functional variant in the serotonin receptor 7 gene (HTR7), rs7905446, is associated with good response to SSRIs in bipolar and unipolar depression.

Predicting antidepressant response has been a clinical challenge for mood disorder. Although several genome-wide association studies have suggested a number of genetic variants to be associated with antidepressant response, the sample sizes are small and the results are difficult to replicate. Previous animal studies have shown that knockout of the serotonin receptor 7 gene (HTR7) resulted in an antidepressant-like phenotype, suggesting it was important to antidepressant action. In this report, in the first stage, we used a cost-effective pooled-sequencing strategy to sequence the entire HTR7 gene and its regulatory regions to investigate the association of common variants in HTR7 and clinical response to four selective serotonin reuptake inhibitors (SSRIs: citalopram, paroxetine, fluoxetine and sertraline) in a retrospective cohort mainly consisting of subjects with bipolar disorder (n = 359). We found 80 single-nucleotide polymorphisms (SNPs) with false discovery rate < 0.05 associated with response to paroxetine. Among the significant SNPs, rs7905446 (T/G), which is located at the promoter region, also showed nominal significance (P < 0.05) in fluoxetine group. GG/TG genotypes for rs7905446 and female gender were associated with better response to two SSRIs (paroxetine and fluoxetine). In the second stage, we replicated this association in two independent prospective samples of SSRI-treated patients with major depressive disorder: the MARS (n = 253, P = 0.0169) and GENDEP studies (n = 432, P = 0.008). The GG/TG genotypes were consistently associated with response in all three samples. Functional study of rs7905446 showed greater activity of the G allele in regulating expression of HTR7. The G allele displayed higher luciferase activity in two neuronal-related cell lines, and estrogen treatment decreased the activity of only the G allele. Electrophoretic mobility shift assay suggested that the G allele interacted with CCAAT/enhancer-binding protein beta transcription factor (TF), while the T allele did not show any interaction with any TFs. Our results provided novel pharmacogenomic evidence to support the role of HTR7 in association with antidepressant response.

SCN11A mRNA levels in female bipolar disorder PBMCs as tentative biomarker for distinct patient sub-phenotypes.

Bipolar disorder (BD) is a complex neuropsychiatric disorder characterized by recurrent mania and depression episodes and requiring lifelong treatment with mood stabilizing drugs. Several lines of evidence, including with BD patient iPSC-derived neurons, suggest that neuronal hyperexcitability may underlie the key clinical symptoms of BD. Indeed, higher mRNA levels of SCN11A, coding for the voltage-gated sodium channel NaV 1.9 implicated in nociception, were detected in iPSC-derived neurons from BD patients, and were normalized by in vitro lithium. Here we studied SCN11A expression in peripheral blood mononuclear cells (PBMCs) from well-phenotyped female BD patients and controls and evaluated their association with several clinical sub-phenotypes. We observed higher mRNA levels of SCN11A in PBMCs from female BD patients with no records of alcohol dependence (p = .0050), no records of psychosis (p = .0097), or no records of suicide attempts (p = .0409). A trend was observed for higher SCN11A expression (FD = 1.91; p = .052) in BD PBMCs compared with controls. Datamining of published postmortem gene expression datasets indicated higher SCN11A expression in dorsolateral prefrontal cortex and orbitofrontal cortex tissues from BD patients compared with controls. Higher phenotype-associated expression levels in PBMC from BD patients were also observed for ID2 (alcohol dependence, suicide attempts) and HDGFRP3 (seasonal BD pattern). Our findings suggest that higher PBMC SCN11A expression levels may be associated with certain behavioral BD sub-phenotypes, including lack of alcohol dependence and psychosis, among BD patients. The NaV 1.9 voltage-gated sodium channel thus deserves consideration as a tentative phenotype modifier in BD.

Entrainment of Circadian Rhythms to Temperature Reveals Amplitude Deficits in Fibroblasts from Patients with Bipolar Disorder and Possible Links to Calcium Channels.

Bipolar disorder (BD) is characterized by recurrent mood episodes, and circadian rhythm disturbances. Past studies have identified calcium channel genes as risk loci for BD. CACNA1C encodes an L-type calcium channel (LTCC) involved in the entrainment of circadian rhythms to light. Another calcium channel, i.e., the ryanodine receptor (RYR), is involved in -circadian phase delays. It is unknown whether variants in CACNA1C or other calcium channels contribute to the circadian phenotype in BD. We hypothesized that, by using temperature cycles, we could model circadian entrainment in fibroblasts from BD patients and controls to interrogate the circadian functions of LTCCs. Using Per2-luc, a bioluminescent reporter, we verified that cells entrain to temperature rhythms in vitro. Under constant temperature conditions, the LTCC antagonist verapamil shortened the circadian period, and the RYR antagonist dantrolene lengthened the period. However, neither drug affected temperature entrainment. Fibroblasts from BD patients and controls also entrained to temperature. In cells from BD patients, the rhythm amplitude was lower under entrained, but not constant, conditions. Temperature entrainment was otherwise similar between BD and control cells. However, the CACNA1C genotype among BD cells predicted the degree to which cells entrained. We conclude that assessment of rhythms under entrained conditions reveals additional rhythm abnormalities in BD that are not observable under constant temperature conditions.

Study of 45 candidate genes suggests CACNG2 may be associated with lithium response in bipolar disorder.

BACKGROUND: Bipolar disorder is a neuropsychiatric disorder that is characterized by fluctuations between manic and depressive phases. Lithium is the original and best mood stabilizing treatment for bipolar disorder. While its mechanism is not well understood, it is believed to have a strong genetic component, as several studies suggest that lithium responsiveness, in bipolar disorder, is heritable. In this study we aimed to identify genetic variants that are associated with lithium responsiveness in bipolar disorder. METHODS: Here we present two cohorts; a retrospective cohort in which patients were surveyed about their response to lithium, and a prospective cohort, in which patients were placed on a lithium monotherapy and monitored for their response to lithium. In both cohorts, patients were stratified into two categories in terms of lithium response; good responders and poor responders. 45 genes were selected based on previous associations with lithium pathways or bipolar disorder and 684 SNPs within these genes were selected to test for association with lithium response. RESULTS: While no single SNP was significant after correcting for multiple comparisons, there were several that were nominally significant (p < 0.05). Of these nominally significant SNPs, the most highly significant SNP in both the prospective and retrospective cohorts were found to be in CACNG2, or Stargazin. The second best association with lithium response was several SNPs in NRG1, a gene that has previously been associated with schizophrenia. CONCLUSIONS: Evidence for the association of lithium response with SNPs in CACNG2 is consistent with previous findings that have identified CACNG2 as associated with both bipolar disorder and lithium responsiveness.

Efficient region-based test strategy uncovers genetic risk factors for functional outcome in bipolar disorder.

Genome-wide association studies of case-control status have advanced the understanding of the genetic basis of psychiatric disorders. Further progress may be gained by increasing sample size but also by new analysis strategies that advance the exploitation of existing data, especially for clinically important quantitative phenotypes. The functionally-informed efficient region-based test strategy (FIERS) introduced herein uses prior knowledge on biological function and dependence of genotypes within a powerful statistical framework with improved sensitivity and specificity for detecting consistent genetic effects across studies. As proof of concept, FIERS was used for the first genome-wide single nucleotide polymorphism (SNP)-based investigation on bipolar disorder (BD) that focuses on an important aspect of disease course, the functional outcome. FIERS identified a significantly associated locus on chromosome 15 (hg38: chr15:48965004 - 49464789 bp) with consistent effect strength between two independent studies (GAIN/TGen: European Americans, BOMA: Germans; n = 1592 BD patients in total). Protective and risk haplotypes were found on the most strongly associated SNPs. They contain a CTCF binding site (rs586758); CTCF sites are known to regulate sets of genes within a chromatin domain. The rs586758 - rs2086256 - rs1904317 haplotype is located in the promoter flanking region of the COPS2 gene, close to microRNA4716, and the EID1, SHC4, DTWD1 genes as plausible biological candidates. While implication with BD is novel, COPS2, EID1, and SHC4 are known to be relevant for neuronal differentiation and function and DTWD1 for psychopharmacological side effects. The test strategy FIERS that enabled this discovery is equally applicable for tag SNPs and sequence data.

RNA sequencing of bipolar disorder lymphoblastoid cell lines implicates the neurotrophic factor HRP-3 in lithium's clinical efficacy.

Objectives: Lithium remains the oldest and most effective treatment for mood stabilisation in bipolar disorder (BD), even though at least half of patients are only partially responsive or do not respond. This study aimed to identify biomarkers associated with lithium response in BD, based on comparing RNA sequencing information derived from lymphoblastoid cell lines (LCLs) of lithium-responsive (LR) versus lithium non-responsive (LNR) BD patients, to assess gene expression variations that might bear on treatment outcome. Methods: RNA sequencing was carried out on 24 LCLs from female BD patients (12 LR and 12 LNR) followed by qPCR validation in two additional independent cohorts (41 and 17 BD patients, respectively). Results: Fifty-six genes showed nominal differential expression comparing LR and LNR (FC ≥ |1.3|, P ≤ 0.01). The differential expression of HDGFRP3 and ID2 was validated by qPCR in the independent cohorts. Conclusions: We observed higher expression levels of HDGFRP3 and ID2 in BD patients who favourably respond to lithium. Both of these genes are involved in neurogenesis, and HDGFRP3 has been suggested to be a neurotrophic factor. Additional studies in larger BD cohorts are needed to confirm the potential of HDGFRP3 and ID2 expression levels in blood cells as tentative favourable lithium response biomarkers.

Detecting significant genotype-phenotype association rules in bipolar disorder: market research meets complex genetics.

BACKGROUND: Disentangling the etiology of common, complex diseases is a major challenge in genetic research. For bipolar disorder (BD), several genome-wide association studies (GWAS) have been performed. Similar to other complex disorders, major breakthroughs in explaining the high heritability of BD through GWAS have remained elusive. To overcome this dilemma, genetic research into BD, has embraced a variety of strategies such as the formation of large consortia to increase sample size and sequencing approaches. Here we advocate a complementary approach making use of already existing GWAS data: a novel data mining procedure to identify yet undetected genotype-phenotype relationships. We adapted association rule mining, a data mining technique traditionally used in retail market research, to identify frequent and characteristic genotype patterns showing strong associations to phenotype clusters. We applied this strategy to three independent GWAS datasets from 2835 phenotypically characterized patients with BD. In a discovery step, 20,882 candidate association rules were extracted. RESULTS: Two of these rules-one associated with eating disorder and the other with anxiety-remained significant in an independent dataset after robust correction for multiple testing. Both showed considerable effect sizes (odds ratio ~ 3.4 and 3.0, respectively) and support previously reported molecular biological findings. CONCLUSION: Our approach detected novel specific genotype-phenotype relationships in BD that were missed by standard analyses like GWAS. While we developed and applied our method within the context of BD gene discovery, it may facilitate identifying highly specific genotype-phenotype relationships in subsets of genome-wide data sets of other complex phenotype with similar epidemiological properties and challenges to gene discovery efforts.

A gene co-expression module implicating the mitochondrial electron transport chain is associated with long-term response to lithium treatment in bipolar affective disorder.

Lithium is the first-line treatment for bipolar affective disorder (BPAD) but two-thirds of patients respond only partially or not at all. The reasons for this high variability in lithium response are not well understood. Transcriptome-wide profiling, which tests the interface between genes and the environment, represents a viable means of exploring the molecular mechanisms underlying lithium response variability. Thus, in the present study we performed co-expression network analyses of whole-blood-derived RNA-seq data from n = 50 lithium-treated BPAD patients. Lithium response was assessed using the well-validated ALDA scale, which we used to define both a continuous and a dichotomous measure. We identified a nominally significant correlation between a co-expression module comprising 46 genes and lithium response represented as a continuous (i.e., scale ranging 0-10) phenotype (cor = -0.299, p = 0.035). Forty-three of these 46 genes had reduced mRNA expression levels in better lithium responders relative to poorer responders, and the central regulators of this module were all mitochondrially-encoded (MT-ND1, MT-ATP6, MT-CYB). Accordingly, enrichment analyses indicated that genes involved in mitochondrial functioning were heavily over-represented in this module, specifically highlighting the electron transport chain (ETC) and oxidative phosphorylation (OXPHOS) as affected processes. Disrupted ETC and OXPHOS activity have previously been implicated in the pathophysiology of BPAD. Our data adds to previous evidence suggesting that a normalisation of these processes could be central to lithium's mode of action, and could underlie a favourable therapeutic response.

Association of Polygenic Score for Schizophrenia and HLA Antigen and Inflammation Genes With Response to Lithium in Bipolar Affective Disorder: A Genome-Wide Association Study.

Importance: Lithium is a first-line mood stabilizer for the treatment of bipolar affective disorder (BPAD). However, the efficacy of lithium varies widely, with a nonresponse rate of up to 30%. Biological response markers are lacking. Genetic factors are thought to mediate treatment response to lithium, and there is a previously reported genetic overlap between BPAD and schizophrenia (SCZ). Objectives: To test whether a polygenic score for SCZ is associated with treatment response to lithium in BPAD and to explore the potential molecular underpinnings of this association. Design, Setting, and Participants: A total of 2586 patients with BPAD who had undergone lithium treatment were genotyped and assessed for long-term response to treatment between 2008 and 2013. Weighted SCZ polygenic scores were computed at different P value thresholds using summary statistics from an international multicenter genome-wide association study (GWAS) of 36 989 individuals with SCZ and genotype data from patients with BPAD from the Consortium on Lithium Genetics. For functional exploration, a cross-trait meta-GWAS and pathway analysis was performed, combining GWAS summary statistics on SCZ and response to treatment with lithium. Data analysis was performed from September 2016 to February 2017. Main Outcomes and Measures: Treatment response to lithium was defined on both the categorical and continuous scales using the Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder score. The effect measures include odds ratios and the proportion of variance explained. Results: Of the 2586 patients in the study (mean [SD] age, 47.2 [13.9] years), 1478 were women and 1108 were men. The polygenic score for SCZ was inversely associated with lithium treatment response in the categorical outcome, at a threshold P < 5 × 10-2. Patients with BPAD who had a low polygenic load for SCZ responded better to lithium, with odds ratios for lithium response ranging from 3.46 (95% CI, 1.42-8.41) at the first decile to 2.03 (95% CI, 0.86-4.81) at the ninth decile, compared with the patients in the 10th decile of SCZ risk. In the cross-trait meta-GWAS, 15 genetic loci that may have overlapping effects on lithium treatment response and susceptibility to SCZ were identified. Functional pathway and network analysis of these loci point to the HLA antigen complex and inflammatory cytokines. Conclusions and Relevance: This study provides evidence for a negative association between high genetic loading for SCZ and poor response to lithium in patients with BPAD. These results suggest the potential for translational research aimed at personalized prescribing of lithium.