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

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

Brain-derived neurotrophic factor (BDNF) gene: no major impact on antidepressant treatment response

The brain-derived neurotrophic factor (BDNF) has been suggested to play a pivotal role in the aetiology of affective disorders. In order to further clarify the impact of BDNF gene variation on major depression as well as antidepressant treatment response, association of three BDNF polymorphisms [rs7103411, Val66Met (rs6265) and rs7124442] with major depression and antidepressant treatment response was investigated in an overall sample of 268 German patients with major depression and 424 healthy controls. False discovery rate (FDR) was applied to control for multiple testing. Additionally, ten markers in BDNF were tested for association with citalopram outcome in the STAR*D sample. While BDNF was not associated with major depression as a categorical diagnosis, the BDNF rs7124442 TT genotype was significantly related to worse treatment outcome over 6 wk in major depression (p=0.01) particularly in anxious depression (p=0.003) in the German sample. However, BDNF rs7103411 and rs6265 similarly predicted worse treatment response over 6 wk in clinical subtypes of depression such as melancholic depression only (rs7103411: TT<CC, p=0.003; rs6265: GG<AA, p=0.001). All SNPs had main effects on antidepressant treatment response in ANOVA models when the remaining SNPs were considered as covariates. The STAR*D analyses did not yield significant results at any of the ten BDNF markers. Our results do not support an association between genetic variation in BDNF and antidepressant treatment response or remission. Post-hoc analyses provide some preliminary support for a potential minor role of genetic variation in BDNF and antidepressant treatment outcome in the context of melancholic depression

Genetic variation in htr2a influences serotonin transporter binding potential as measured using pet and [11c]dasb

In a previous study we showed that genetic variation in HTR2A, which encodes the serotonin 2A receptor, influenced outcome of citalopram treatment in patients with major depressive disorder. Since chronic administration of citalopram, which selectively and potently inhibits the serotonin transporter (5-HTT), putatively enhances serotonergic transmission, it is conceivable that genetic variation within HTR2A also influences pretreatment 5-HTT function or serotonergic transmission. The present study used positron emission tomography (PET) and the selective 5-HTT ligand, [C-11]DASB, to investigate whether the HTR2A marker alleles that predict treatment outcome also predict differences in 5-HTT binding. Brain levels of 5-HTT were assessed in vivo using PET measures of the non-displaceable component of the [C-11]DASB binding potential (BPND). DNA from 43 patients and healthy volunteers, all unmedicated, was genotyped with 14 single nucleotide polymorphisms located within or around HTR2A. Allelic association with BPND was assessed in eight brain regions, with covariates to control for race and ethnicity. We detected allelic association between [(11)]DASB BPND in thalamus and three markers in a region spanning the 3\u27 untranslated region and second intron of HTR2A (rs7333412, p=0.000045; rs7997012, p=0.000086; rs977003, p=0.000069). The association signal at rs7333412 remained significant (p&amp;lt;0.05) after applying corrections for multiple testing via permutation. Genetic variation in HTR2A that was previously associated with citalopram treatment outcome was also associated with thalamic 5-HTT binding. While further work is needed to identify the actual functional genetic variants involved, these results suggest that a relationship exists between genetic variation in HTR2A and either 5-HTT expression or central serotonergic transmission that influences the therapeutic response to 5-HTT inhibition in major depression

Citalopram-induced pathways regulation and tentative treatment-outcome-predicting biomarkers in lymphoblastoid cell lines from depression patients

Antidepressant therapy is still associated with delays in symptomatic improvement and low response rates. Incomplete understanding of molecular mechanisms underlying antidepressant effects hampered the identification of objective biomarkers for antidepressant response. In this work, we studied transcriptome-wide expression followed by pathway analysis in lymphoblastoid cell lines (LCLs) derived from 17 patients documented for response to SSRI antidepressants from the Munich Antidepressant Response Signatures (MARS) study upon short-term incubation (24 and 48 h) with citalopram. Candidate transcripts were further validated with qPCR in MARS LCLs from responders (n = 33) vs. non-responders (n = 36) and afterward in an independent cohort of treatment-resistant patients (n = 20) vs. first-line responders (n = 24) from the STAR*D study. In MARS cohort we observed significant associations of GAD1 (glutamate decarboxylase 1; p = 0.045), TBC1D9 (TBC1 Domain Family Member 9; p = 0.014–0.021) and NFIB (nuclear factor I B; p = 0.015–0.025) expression with response status, remission status and improvement in depression scale, respectively. Pathway analysis of citalopram-altered gene expression indicated response-status-dependent transcriptional reactions. Whereas in clinical responders neural function pathways were primarily up- or downregulated after incubation with citalopram, deregulated pathways in non-responders LCLs mainly involved cell adhesion and immune response. Results from the STAR*D study showed a marginal association of treatment-resistant depression with NFIB (p = 0.068) but not with GAD1 (p = 0.23) and TBC1D9 (p = 0.27). Our results propose the existence of distinct pathway regulation mechanisms in responders vs. non-responders and suggest GAD1, TBC1D9, and NFIB as tentative predictors for clinical response, full remission, and improvement in depression scale, respectively, with only a weak overlap in predictors of different therapy outcome phenotypes