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.

Using polygenic scores and clinical data for bipolar disorder patient stratification and lithium response prediction: machine learning approach.

BACKGROUND: Response to lithium in patients with bipolar disorder is associated with clinical and transdiagnostic genetic factors. The predictive combination of these variables might help clinicians better predict which patients will respond to lithium treatment. AIMS: To use a combination of transdiagnostic genetic and clinical factors to predict lithium response in patients with bipolar disorder. METHOD: This study utilised genetic and clinical data (n = 1034) collected as part of the International Consortium on Lithium Genetics (ConLi+Gen) project. Polygenic risk scores (PRS) were computed for schizophrenia and major depressive disorder, and then combined with clinical variables using a cross-validated machine-learning regression approach. Unimodal, multimodal and genetically stratified models were trained and validated using ridge, elastic net and random forest regression on 692 patients with bipolar disorder from ten study sites using leave-site-out cross-validation. All models were then tested on an independent test set of 342 patients. The best performing models were then tested in a classification framework. RESULTS: The best performing linear model explained 5.1% (P = 0.0001) of variance in lithium response and was composed of clinical variables, PRS variables and interaction terms between them. The best performing non-linear model used only clinical variables and explained 8.1% (P = 0.0001) of variance in lithium response. A priori genomic stratification improved non-linear model performance to 13.7% (P = 0.0001) and improved the binary classification of lithium response. This model stratified patients based on their meta-polygenic loadings for major depressive disorder and schizophrenia and was then trained using clinical data. CONCLUSIONS: Using PRS to first stratify patients genetically and then train machine-learning models with clinical predictors led to large improvements in lithium response prediction. When used with other PRS and biological markers in the future this approach may help inform which patients are most likely to respond to lithium treatment.

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.

Genome-wide association study of 40,000 individuals identifies two novel loci associated with bipolar disorder.

Bipolar disorder (BD) is a genetically complex mental illness characterized by severe oscillations of mood and behaviour. Genome-wide association studies (GWAS) have identified several risk loci that together account for a small portion of the heritability. To identify additional risk loci, we performed a two-stage meta-analysis of >9 million genetic variants in 9,784 bipolar disorder patients and 30,471 controls, the largest GWAS of BD to date. In this study, to increase power we used ∼2,000 lithium-treated cases with a long-term diagnosis of BD from the Consortium on Lithium Genetics, excess controls, and analytic methods optimized for markers on the X-chromosome. In addition to four known loci, results revealed genome-wide significant associations at two novel loci: an intergenic region on 9p21.3 (rs12553324, P = 5.87 × 10 - 9; odds ratio (OR) = 1.12) and markers within ERBB2 (rs2517959, P = 4.53 × 10 - 9; OR = 1.13). No significant X-chromosome associations were detected and X-linked markers explained very little BD heritability. The results add to a growing list of common autosomal variants involved in BD and illustrate the power of comparing well-characterized cases to an excess of controls in GWAS.

Genetic variants associated with response to lithium treatment in bipolar disorder: a genome-wide association study.

BACKGROUND: Lithium is a first-line treatment in bipolar disorder, but individual response is variable. Previous studies have suggested that lithium response is a heritable trait. However, no genetic markers of treatment response have been reproducibly identified. METHODS: Here, we report the results of a genome-wide association study of lithium response in 2563 patients collected by 22 participating sites from the International Consortium on Lithium Genetics (ConLiGen). Data from common single nucleotide polymorphisms (SNPs) were tested for association with categorical and continuous ratings of lithium response. Lithium response was measured using a well established scale (Alda scale). Genotyped SNPs were used to generate data at more than 6 million sites, using standard genomic imputation methods. Traits were regressed against genotype dosage. Results were combined across two batches by meta-analysis. FINDINGS: A single locus of four linked SNPs on chromosome 21 met genome-wide significance criteria for association with lithium response (rs79663003, p=1·37 × 10(-8); rs78015114, p=1·31 × 10(-8); rs74795342, p=3·31 × 10(-9); and rs75222709, p=3·50 × 10(-9)). In an independent, prospective study of 73 patients treated with lithium monotherapy for a period of up to 2 years, carriers of the response-associated alleles had a significantly lower rate of relapse than carriers of the alternate alleles (p=0·03268, hazard ratio 3·8, 95% CI 1·1-13·0). INTERPRETATION: The response-associated region contains two genes for long, non-coding RNAs (lncRNAs), AL157359.3 and AL157359.4. LncRNAs are increasingly appreciated as important regulators of gene expression, particularly in the CNS. Confirmed biomarkers of lithium response would constitute an important step forward in the clinical management of bipolar disorder. Further studies are needed to establish the biological context and potential clinical utility of these findings. FUNDING: Deutsche Forschungsgemeinschaft, National Institute of Mental Health Intramural Research Program.

Insulin resistance and outcome in bipolar disorder.

BACKGROUND: Little is known about the impact of insulin resistance on bipolar disorder. AIMS: To examine the relationships between insulin resistance, type 2 diabetes and clinical course and treatment outcomes in bipolar disorder. METHOD: We measured fasting glucose and insulin in 121 adults with bipolar disorder. We diagnosed type 2 diabetes and determined insulin resistance. The National Institute of Mental Health Life Chart was used to record the course of bipolar disorder and the Alda scale to establish response to prophylactic lithium treatment. RESULTS: Patients with bipolar disorder and type 2 diabetes or insulin resistance had three times higher odds of a chronic course of bipolar disorder compared with euglycaemic patients (50% and 48.7% respectively v. 27.3%, odds ratio (OR) = 3.07, P = 0.007), three times higher odds of rapid cycling (38.5% and 39.5% respectively v. 18.2%, OR = 3.13, P = 0.012) and were more likely to be refractory to lithium treatment (36.8% and 36.7% respectively v. 3.2%, OR = 8.40, P<0.0001). All associations remained significant after controlling for antipsychotic exposure and body mass index in sensitivity analyses. CONCLUSIONS: Comorbid insulin resistance may be an important factor in resistance to treatment in bipolar disorder.

Exploring the genetics of lithium response in bipolar disorders.

BACKGROUND: Lithium (Li) remains the treatment of choice for bipolar disorders (BP). Its mood-stabilizing effects help reduce the long-term burden of mania, depression and suicide risk in patients with BP. It also has been shown to have beneficial effects on disease-associated conditions, including sleep and cardiovascular disorders. However, the individual responses to Li treatment vary within and between diagnostic subtypes of BP (e.g. BP-I and BP-II) according to the clinical presentation. Moreover, long-term Li treatment has been linked to adverse side-effects that are a cause of concern and non-adherence, including the risk of developing chronic medical conditions such as thyroid and renal disease. In recent years, studies by the Consortium on Lithium Genetics (ConLiGen) have uncovered a number of genetic factors that contribute to the variability in Li treatment response in patients with BP. Here, we leveraged the ConLiGen cohort (N = 2064) to investigate the genetic basis of Li effects in BP. For this, we studied how Li response and linked genes associate with the psychiatric symptoms and polygenic load for medical comorbidities, placing particular emphasis on identifying differences between BP-I and BP-II. RESULTS: We found that clinical response to Li treatment, measured with the Alda scale, was associated with a diminished burden of mania, depression, substance and alcohol abuse, psychosis and suicidal ideation in patients with BP-I and, in patients with BP-II, of depression only. Our genetic analyses showed that a stronger clinical response to Li was modestly related to lower polygenic load for diabetes and hypertension in BP-I but not BP-II. Moreover, our results suggested that a number of genes that have been previously linked to Li response variability in BP differentially relate to the psychiatric symptomatology, particularly to the numbers of manic and depressive episodes, and to the polygenic load for comorbid conditions, including diabetes, hypertension and hypothyroidism. CONCLUSIONS: Taken together, our findings suggest that the effects of Li on symptomatology and comorbidity in BP are partially modulated by common genetic factors, with differential effects between BP-I and BP-II.

Lithium response in bipolar disorder is associated with focal adhesion and PI3K-Akt networks: a multi-omics replication study.

Lithium is the gold standard treatment for bipolar disorder (BD). However, its mechanism of action is incompletely understood, and prediction of treatment outcomes is limited. In our previous multi-omics study of the Pharmacogenomics of Bipolar Disorder (PGBD) sample combining transcriptomic and genomic data, we found that focal adhesion, the extracellular matrix (ECM), and PI3K-Akt signaling networks were associated with response to lithium. In this study, we replicated the results of our previous study using network propagation methods in a genome-wide association study of an independent sample of 2039 patients from the International Consortium on Lithium Genetics (ConLiGen) study. We identified functional enrichment in focal adhesion and PI3K-Akt pathways, but we did not find an association with the ECM pathway. Our results suggest that deficits in the neuronal growth cone and PI3K-Akt signaling, but not in ECM proteins, may influence response to lithium in BD.

Exploring the genetics of lithium response in bipolar disorders.

Background: Lithium (Li) remains the treatment of choice for bipolar disorders (BP). Its mood-stabilizing effects help reduce the long-term burden of mania, depression and suicide risk in patients with BP. It also has been shown to have beneficial effects on disease-associated conditions, including sleep and cardiovascular disorders. However, the individual responses to Li treatment vary within and between diagnostic subtypes of BP (e.g. BP-I and BP-II) according to the clinical presentation. Moreover, long-term Li treatment has been linked to adverse side-effects that are a cause of concern and non-adherence, including the risk of developing chronic medical conditions such as thyroid and renal disease. In recent years, studies by the Consortium on Lithium Genetics (ConLiGen) have uncovered a number of genetic factors that contribute to the variability in Li treatment response in patients with BP. Here, we leveraged the ConLiGen cohort (N=2,064) to investigate the genetic basis of Li effects in BP. For this, we studied how Li response and linked genes associate with the psychiatric symptoms and polygenic load for medical comorbidities, placing particular emphasis on identifying differences between BP-I and BP-II. Results: We found that clinical response to Li treatment, measured with the Alda scale, was associated with a diminished burden of mania, depression, substance and alcohol abuse, psychosis and suicidal ideation in patients with BP-I and, in patients with BP-II, of depression only. Our genetic analyses showed that a stronger clinical response to Li was modestly related to lower polygenic load for diabetes and hypertension in BP-I but not BP-II. Moreover, our results suggested that a number of genes that have been previously linked to Li response variability in BP differentially relate to the psychiatric symptomatology, particularly to the numbers of manic and depressive episodes, and to the polygenic load for comorbid conditions, including diabetes, hypertension and hypothyroidism. Conclusions: Taken together, our findings suggest that the effects of Li on symptomatology and comorbidity in BP are partially modulated by common genetic factors, with differential effects between BP-I and BP-II.

Exemplar scoring identifies genetically separable phenotypes of lithium responsive bipolar disorder.

Predicting lithium response (LiR) in bipolar disorder (BD) may inform treatment planning, but phenotypic heterogeneity complicates discovery of genomic markers. We hypothesized that patients with "exemplary phenotypes"-those whose clinical features are reliably associated with LiR and non-response (LiNR)-are more genetically separable than those with less exemplary phenotypes. Using clinical data collected from people with BD (n = 1266 across 7 centers; 34.7% responders), we computed a "clinical exemplar score," which measures the degree to which a subject's clinical phenotype is reliably predictive of LiR/LiNR. For patients whose genotypes were available (n = 321), we evaluated whether a subgroup of responders/non-responders with the top 25% of clinical exemplar scores (the "best clinical exemplars") were more accurately classified based on genetic data, compared to a subgroup with the lowest 25% of clinical exemplar scores (the "poor clinical exemplars"). On average, the best clinical exemplars of LiR had a later illness onset, completely episodic clinical course, absence of rapid cycling and psychosis, and few psychiatric comorbidities. The best clinical exemplars of LiR and LiNR were genetically separable with an area under the receiver operating characteristic curve of 0.88 (IQR [0.83, 0.98]), compared to 0.66 [0.61, 0.80] (p = 0.0032) among poor clinical exemplars. Variants in the Alzheimer's amyloid-secretase pathway, along with G-protein-coupled receptor, muscarinic acetylcholine, and histamine H1R signaling pathways were informative predictors. This study must be replicated on larger samples and extended to predict response to other mood stabilizers.

Association of Attention-Deficit/Hyperactivity Disorder and Depression Polygenic Scores with Lithium Response: A Consortium for Lithium Genetics Study.

Response to lithium varies widely between individuals with bipolar disorder (BD). Polygenic risk scores (PRSs) can uncover pharmacogenomics effects and may help predict drug response. Patients (N = 2,510) with BD were assessed for long-term lithium response in the Consortium on Lithium Genetics using the Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder score. PRSs for attention-deficit/hyperactivity disorder (ADHD), major depressive disorder (MDD), and schizophrenia (SCZ) were computed using lassosum and in a model including all three PRSs and other covariates, and the PRS of ADHD (ß = -0.14; 95% confidence interval [CI]: -0.24 to -0.03; p value = 0.010) and MDD (ß = -0.16; 95% CI: -0.27 to -0.04; p value = 0.005) predicted worse quantitative lithium response. A higher SCZ PRS was associated with higher rates of medication nonadherence (OR = 1.61; 95% CI: 1.34-1.93; p value = 2e-7). This study indicates that genetic risk for ADHD and depression may influence lithium treatment response. Interestingly, a higher SCZ PRS was associated with poor adherence, which can negatively impact treatment response. Incorporating genetic risk of ADHD, depression, and SCZ in combination with clinical risk may lead to better clinical care for patients with BD.

Combining schizophrenia and depression polygenic risk scores improves the genetic prediction of lithium response in bipolar disorder patients.

Lithium is the gold standard therapy for Bipolar Disorder (BD) but its effectiveness differs widely between individuals. The molecular mechanisms underlying treatment response heterogeneity are not well understood, and personalized treatment in BD remains elusive. Genetic analyses of the lithium treatment response phenotype may generate novel molecular insights into lithium's therapeutic mechanisms and lead to testable hypotheses to improve BD management and outcomes. We used fixed effect meta-analysis techniques to develop meta-analytic polygenic risk scores (MET-PRS) from combinations of highly correlated psychiatric traits, namely schizophrenia (SCZ), major depression (MD) and bipolar disorder (BD). We compared the effects of cross-disorder MET-PRS and single genetic trait PRS on lithium response. For the PRS analyses, we included clinical data on lithium treatment response and genetic information for n = 2283 BD cases from the International Consortium on Lithium Genetics (ConLi+Gen; www.ConLiGen.org ). Higher SCZ and MD PRSs were associated with poorer lithium treatment response whereas BD-PRS had no association with treatment outcome. The combined MET2-PRS comprising of SCZ and MD variants (MET2-PRS) and a model using SCZ and MD-PRS sequentially improved response prediction, compared to single-disorder PRS or to a combined score using all three traits (MET3-PRS). Patients in the highest decile for MET2-PRS loading had 2.5 times higher odds of being classified as poor responders than patients with the lowest decile MET2-PRS scores. An exploratory functional pathway analysis of top MET2-PRS variants was conducted. Findings may inform the development of future testing strategies for personalized lithium prescribing in BD.

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.

Antidepressant monotherapy in pre-bipolar depression; predictive value and inherent risk.

OBJECTIVE: To identify specific treatment-emergent symptoms in response to antidepressant therapy in depression preceding bipolar disorder. METHODS: Retrospective chart review of response to antidepressants in "pre-bipolar" depression, compared to a matched unipolar sample. RESULTS: Family history of completed suicide (p=0.0003) and bipolar disorder (p=0.004) were more common in the pre-bipolar subgroup. Earlier age of onset of diagnosed depression (p=0.005) as well as even earlier episodes of untreated retrospectively diagnosed major depression (p<0.0001) were associated with a future bipolar course. The pre-bipolar group was less likely to respond to antidepressant treatment (p=0.009). Treatment-emergent "mixed" symptoms (two or more symptoms of DSM IV mania, mood lability, irritability/rage with co-existing depression) and in particular, "serious symptoms" (treatment emergent or increased agitation, rage or suicidality) occurred more commonly in the bipolar group. The two variables that best accounted for the between-group differences in logistic regression, were early age at first symptoms of depression and treatment-emergent agitation. CONCLUSIONS: Family history of completed suicide and/or bipolar disorder, early onset of depressive symptoms as well as treatment-emergent "mixed" symptoms are common in depression preceding the diagnosis of bipolar disorder.

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.

Analysis of the Influence of microRNAs in Lithium Response in Bipolar Disorder.

Bipolar disorder (BD) is a common, highly heritable neuropsychiatric disease characterized by recurrent episodes of mania and depression. Lithium is the best-established long-term treatment for BD, even though individual response is highly variable. Evidence suggests that some of this variability has a genetic basis. This is supported by the largest genome-wide association study (GWAS) of lithium response to date conducted by the International Consortium on Lithium Genetics (ConLiGen). Recently, we performed the first genome-wide analysis of the involvement of miRNAs in BD and identified nine BD-associated miRNAs. However, it is unknown whether these miRNAs are also associated with lithium response in BD. In the present study, we therefore tested whether common variants at these nine candidate miRNAs contribute to the variance in lithium response in BD. Furthermore, we systematically analyzed whether any other miRNA in the genome is implicated in the response to lithium. For this purpose, we performed gene-based tests for all known miRNA coding genes in the ConLiGen GWAS dataset (n = 2,563 patients) using a set-based testing approach adapted from the versatile gene-based test for GWAS (VEGAS2). In the candidate approach, miR-499a showed a nominally significant association with lithium response, providing some evidence for involvement in both development and treatment of BD. In the genome-wide miRNA analysis, 71 miRNAs showed nominally significant associations with the dichotomous phenotype and 106 with the continuous trait for treatment response. A total of 15 miRNAs revealed nominal significance in both phenotypes with miR-633 showing the strongest association with the continuous trait (p = 9.80E-04) and miR-607 with the dichotomous phenotype (p = 5.79E-04). No association between miRNAs and treatment response to lithium in BD in either of the tested conditions withstood multiple testing correction. Given the limited power of our study, the investigation of miRNAs in larger GWAS samples of BD and lithium response is warranted.