Genome-wide association study of antidepressant treatment resistance in a population-based cohort using health service prescription data and meta-analysis with GENDEP

Susceptibility to Treatment-Resistant Depression Within Families

Importance

Antidepressant responses and the phenotype of treatment-resistant depression (TRD) are believed to have a genetic basis. Genetic susceptibility between the TRD phenotype and other psychiatric disorders has also been established in previous genetic studies, but population-based cohort studies have not yet provided evidence to support these outcomes.

Objective

To estimate the TRD susceptibility and the susceptibility between TRD and other psychiatric disorders within families in a nationwide insurance cohort with extremely high coverage and comprehensive health care data.

Design, Setting, and Participants

This cohort study assessed data from the Taiwan national health insurance database across entire population (N = 26 554 001) between January 2003 and December 2017. Data analysis was performed from August 2021 to April 2023. TRD was defined as having experienced at least 3 distinct antidepressant treatments in the current episode, each with adequate dose and duration, based on the prescribing records. Then, we identified the first-degree relatives of individuals with TRD (n = 34 467). A 1:4 comparison group (n = 137 868) of first-degree relatives of individuals without TRD was arranged for the comparison group, matched by birth year, sex, and kinship.

Main Outcomes and Measures

Modified Poisson regression analyses were performed and adjusted relative risks (aRRs) and 95% CIs were calculated for the risk of TRD, the risk of other major psychiatric disorders, and different causes of mortality.

Results

This study included 172 335 participants (88 330 male and 84 005 female; mean [SD] age at beginning of follow-up, 22.9 [18.1] years). First-degree relatives of individuals with TRD had lower incomes, more physical comorbidities, higher suicide mortality, and increased risk of developing TRD (aRR, 9.16; 95% CI, 7.21-11.63) and higher risk of other psychiatric disorders than matched control individuals, including schizophrenia (aRR, 2.36; 95% CI, 2.10-2.65), bipolar disorder (aRR, 3.74; 95% CI, 3.39-4.13), major depressive disorder (aRR, 3.65; 95% CI, 3.44-3.87), attention-deficit/hyperactivity disorders (aRR, 2.38; 95% CI, 2.20-2.58), autism spectrum disorder (aRR, 2.26; 95% CI, 1.86-2.74), anxiety disorder (aRR, 2.71; 95% CI, 2.59-2.84), and obsessive-compulsive disorder (aRR, 3.14; 95% CI, 2.70-3.66). Sensitivity and subgroup analyses validated the robustness of the findings.

Conclusions and Relevance

To our knowledge, this study is the largest and perhaps first nationwide cohort study to demonstrate TRD phenotype transmission across families and coaggregation with other major psychiatric disorders. Patients with a family history of TRD had an increased risk of suicide mortality and tendency toward antidepressant resistance; therefore, more intensive treatments for depressive symptoms might be considered earlier, rather than antidepressant monotherapy.

Applying polygenic risk score methods to pharmacogenomics GWAS: challenges and opportunities

Polygenic risk scores (PRSs) have emerged as promising tools for the prediction of human diseases and complex traits in disease genome-wide association studies (GWAS). Applying PRSs to pharmacogenomics (PGx) studies has begun to show great potential for improving patient stratification and drug response prediction. However, there are unique challenges that arise when applying PRSs to PGx GWAS beyond those typically encountered in disease GWAS (e.g. Eurocentric or trans-ethnic bias). These challenges include: (i) the lack of knowledge about whether PGx or disease GWAS/variants should be used in the base cohort (BC); (ii) the small sample sizes in PGx GWAS with corresponding low power and (iii) the more complex PRS statistical modeling required for handling both prognostic and predictive effects simultaneously. To gain insights in this landscape about the general trends, challenges and possible solutions, we first conduct a systematic review of both PRS applications and PRS method development in PGx GWAS. To further address the challenges, we propose (i) a novel PRS application strategy by leveraging both PGx and disease GWAS summary statistics in the BC for PRS construction and (ii) a new Bayesian method (PRS-PGx-Bayesx) to reduce Eurocentric or cross-population PRS prediction bias. Extensive simulations are conducted to demonstrate their advantages over existing PRS methods applied in PGx GWAS. Our systematic review and methodology research work not only highlights current gaps and key considerations while applying PRS methods to PGx GWAS, but also provides possible solutions for better PGx PRS applications and future research.

Polygenic risk scores of lithium response and treatment resistance in major depressive disorder

Treatment response and resistance in major depressive disorder (MDD) are suggested to be heritable. Due to significant challenges in defining treatment-related phenotypes, our understanding of their genetic bases is limited. This study aimed to derive a stringent definition of treatment resistance and to investigate the genetic overlap between treatment response and resistance in MDD. Using electronic medical records on the use of antidepressants and electroconvulsive therapy (ECT) from Swedish registers, we derived the phenotype of treatment-resistant depression (TRD) and non-TRD within ~4500 individuals with MDD in three Swedish cohorts. Considering antidepressants and lithium are first-line treatment and augmentation used for MDD, respectively, we generated polygenic risk scores (PRS) of antidepressants and lithium response for individuals with MDD and evaluated their associations with treatment resistance by comparing TRD with non-TRD. Among 1778 ECT-treated MDD cases, nearly all (94%) used antidepressants before their first ECT and the vast majority had at least one (84%) or two (61%) antidepressants of adequate duration, suggesting these MDD cases receiving ECT were resistant to antidepressants. We did not observe a significant difference in the mean PRS of antidepressant response between TRD and non-TRD; however, we found that TRD cases had a significantly higher PRS of lithium response compared to non-TRD cases (OR = 1.10-1.12 under various definitions). The results support the evidence of heritable components in treatment-related phenotypes and highlight the overall genetic profile of lithium-sensitivity in TRD. This finding further provides a genetic explanation for lithium efficacy in treating TRD.

Exome-wide association study of treatment-resistant depression suggests novel treatment targets

Treatment-resistant depression (TRD) is a severe form of major depressive disorder (MDD) with substantial public health impact and poor treatment outcome. Treatment outcome in MDD is significantly heritable, but genome-wide association studies have failed to identify replicable common marker alleles, suggesting a potential role for uncommon variants. Here we investigated the hypothesis that uncommon, putatively functional genetic variants are associated with TRD. Whole-exome sequencing data was obtained from 182 TRD cases and 2021 psychiatrically healthy controls. After quality control, the remaining 149 TRD cases and 1976 controls were analyzed with tests designed to detect excess burdens of uncommon variants. At the gene level, 5 genes, ZNF248, PRKRA, PYHIN1, SLC7A8, and STK19 each carried exome-wide significant excess burdens of variants in TRD cases (q < 0.05). Analysis of 41 pre-selected gene sets suggested an excess of uncommon, functional variants among genes involved in lithium response. Among the genes identified in previous TRD studies, ZDHHC3 was also significant in this sample after multiple test correction. ZNF248 and STK19 are involved in transcriptional regulation, PHYIN1 and PRKRA are involved in immune response, SLC7A8 is associated with thyroid hormone transporter activity, and ZDHHC3 regulates synaptic clustering of GABA and glutamate receptors. These results implicate uncommon, functional alleles in TRD and suggest promising novel targets for future research.

The prediction of Alzheimer's disease through multi-trait genetic modeling

To better capture the polygenic architecture of Alzheimer's disease (AD), we developed a joint genetic score, MetaGRS. We incorporated genetic variants for AD and 24 other traits from two independent cohorts, NACC (n = 3,174, training set) and UPitt (n = 2,053, validation set). One standard deviation increase in the MetaGRS is associated with about 57% increase in the AD risk [hazard ratio (HR) = 1.577, p = 7.17 E-56], showing little difference from the HR for AD GRS alone (HR = 1.579, p = 1.20E-56), suggesting similar utility of both models. We also conducted APOE-stratified analyses to assess the role of the e4 allele on risk prediction. Similar to that of the combined model, our stratified results did not show a considerable improvement of the MetaGRS. Our study showed that the prediction power of the MetaGRS significantly outperformed that of the reference model without any genetic information, but was effectively equivalent to the prediction power of the AD GRS.

A cultural-ecosocial systems view for psychiatry

While contemporary psychiatry seeks the mechanisms of mental disorders in neurobiology, mental health problems clearly depend on developmental processes of learning and adaptation through ongoing interactions with the social environment. Symptoms or disorders emerge in specific social contexts and involve predicaments that cannot be fully characterized in terms of brain function but require a larger social-ecological view. Causal processes that result in mental health problems can begin anywhere within the extended system of body-person-environment. In particular, individuals' narrative self-construal, culturally mediated interpretations of symptoms and coping strategies as well as the responses of others in the social world contribute to the mechanisms of mental disorders, illness experience, and recovery. In this paper, we outline the conceptual basis and practical implications of a hierarchical ecosocial systems view for an integrative approach to psychiatric theory and practice. The cultural-ecosocial systems view we propose understands mind, brain and person as situated in the social world and as constituted by cultural and self-reflexive processes. This view can be incorporated into a pragmatic approach to clinical assessment and case formulation that characterizes mechanisms of pathology and identifies targets for intervention.

Investigating genetic overlap between antidepressant and lithium response and treatment resistance in major depressive disorder

Treatment response and resistance in major depressive disorder (MDD) are suggested to be heritable. Due to significant challenges in defining treatment-related phenotypes, our understanding of their genetic bases is limited. This study aimed to derive a stringent definition of treatment resistance and to investigate genetic overlap between treatment response and resistance in MDD. Using electronic medical records on the use of antidepressants and electroconvulsive therapy (ECT) from Swedish registers, we derived the phenotype of treatment-resistant depression (TRD) within ~ 4 500 individuals with MDD in three Swedish cohorts. Considering antidepressants and lithium are first-line treatment and augmentation used for MDD, respectively, we generated polygenic risk scores of antidepressant and lithium response for individuals with MDD, and evaluated their associations with treatment resistance by comparing TRD with non-TRD. Among 1 778 ECT-treated MDD cases, nearly all (94%) used antidepressants before first ECT, and the vast majority had at least one (84%) or two (61%) antidepressants of adequate duration, suggesting these MDD cases receiving ECT were resistant to antidepressants. We found that TRD cases tend to have lower genetic load of antidepressant response than non-TRD, although the difference was not significant; furthermore, TRD cases had significantly higher genetic load of lithium response (OR = 1.10-1.12 under different definitions). The results support evidence of heritable components in treatment-related phenotypes and highlight the overall genetic profile of lithium-sensitivity in TRD. This finding further provides a genetic explanation for lithium efficacy in treating TRD.

Precision Medicine in Antidepressants Treatment

Precision medicine uses innovative approaches to improve disease prevention and treatment outcomes by taking into account people's genetic backgrounds, environments, and lifestyles. Treatment of depression is particularly challenging, given that 30-50% of patients do not respond adequately to antidepressants, while those who respond may experience unpleasant adverse drug reactions (ADRs) that decrease their quality of life and compliance. This chapter aims to present the available scientific data that focus on the impact of genetic variants on the efficacy and toxicity of antidepressants. We compiled data from candidate gene and genome-wide association studies that investigated associations between pharmacodynamic and pharmacokinetic genes and response to antidepressants regarding symptom improvement and ADRs. We also summarized the existing pharmacogenetic-based treatment guidelines for antidepressants, used to guide the selection of the right antidepressant and its dose based on the patient's genetic profile, aiming to achieve maximum efficacy and minimum toxicity. Finally, we reviewed the clinical implementation of pharmacogenomics studies focusing on patients on antidepressants. The available data demonstrate that precision medicine can increase the efficacy of antidepressants and reduce the occurrence of ADRs and ultimately improve patients' quality of life.

Understanding treatment-resistant depression using "omics" techniques: A systematic review

BACKGROUND

Treatment-resistant depression (TRD) results in huge healthcare costs and poor patient clinical outcomes. Most studies have adopted a "candidate mechanism" approach to investigate TRD pathogenesis, however this is made more challenging due to the complex and heterogeneous nature of this condition. High-throughput "omics" technologies can provide a more holistic view and further insight into the underlying mechanisms involved in TRD development, expanding knowledge beyond already-identified mechanisms. This systematic review assessed the information from studies that examined TRD using hypothesis-free omics techniques.

METHODS

PubMed, MEDLINE, Embase, APA PsycInfo, Scopus and Web of Science databases were searched on July 2022. 37 human studies met the eligibility criteria, totalling 17,518 TRD patients, 571,402 healthy controls and 62,279 non-TRD depressed patients (including antidepressant responders and untreated MDD patients).

RESULTS

Significant findings were reported that implicate the role in TRD of various molecules, including polymorphisms, genes, mRNAs and microRNAs. The pathways most commonly reported by the identified studies were involved in immune system and inflammation, neuroplasticity, calcium signalling and neurotransmitters.

LIMITATIONS

Small sample sizes, variability in defining TRD, and heterogeneity in study design and methodology.

CONCLUSIONS

These findings provide insight into TRD pathophysiology, proposing future research directions for novel drug targets and potential biomarkers for clinical staging and response to antidepressants (citalopram/escitalopram in particular) and electroconvulsive therapy (ECT). Further validation is warranted in large prospective studies using standardised TRD criteria. A multi-omics and systems biology strategy with a collaborative effort will likely deliver robust findings for translation into the clinic.

Polygenic risk scores for predicting outcomes and treatment response in psychiatry: hope or hype?

Over the last years, the decreased costs and enhanced accessibility to large genome-wide association studies datasets have laid the foundations for the development of polygenic risk scores (PRSs). A PRS is calculated on the weighted sum of single nucleotide polymorphisms and measures the individual genetic predisposition to develop a certain phenotype. An increasing number of studies have attempted to utilize the PRSs for risk stratification and prognostic evaluation. The present narrative review aims to discuss the potential clinical utility of PRSs in predicting outcomes and treatment response in psychiatry. After summarizing the evidence on major mental disorders, we have discussed the advantages and limitations of currently available PRSs. Although PRSs represent stable trait features with a normal distribution in the general population and can be relatively easily calculated in terms of time and costs, their real-world applicability is reduced by several limitations, such as low predictive power and lack of population diversity. Even with the rapid expansion of the psychiatric genetic knowledge base, pure genetic prediction in clinical psychiatry appears to be out of reach in the near future. Therefore, combining genomic and exposomic vulnerabilities for mental disorders with a detailed clinical characterization is needed to personalize care.

Polygenic heterogeneity in antidepressant treatment and placebo response

The genetic architecture of antidepressant response is poorly understood. Polygenic risk scores (PRS), exploration of placebo response and the use of sub-scales might provide insights. Here, we investigate the association between PRSs for relevant complex traits and response to vortioxetine treatment and placebo using clinical scales, including sub-scales and self-reported assessments. We collected a clinical test sample of Major Depressive Disorder (MDD) patients treated with vortioxetine (N = 907) or placebo (N = 455) from seven randomized, double-blind, clinical trials. In parallel, we obtained data from an observational web-based study of vortioxetine-treated patients (N = 642) with self-reported response. PRSs for antidepressant response, psychiatric disorders, and symptom traits were derived using summary statistics from well-powered genome-wide association studies (GWAS). Association tests were performed between the PRSs and treatment response in each of the two test samples and empirical p-values were evaluated. In the clinical test sample, no PRSs were significantly associated with response to vortioxetine treatment or placebo following Bonferroni correction. However, clinically assessed treatment response PRS was nominally associated with vortioxetine treatment and placebo response given by several secondary outcome scales (improvement on HAM-A, HAM-A Psychic Anxiety sub-scale, CPFQ &amp; PDQ), (P ≤ 0.026). Further, higher subjective well-being PRS (P ≤ 0.033) and lower depression PRS (P = 0.01) were nominally associated with higher placebo response. In the self-reported test sample, higher schizophrenia PRS was significantly associated with poorer self-reported response (P = 0.0001). The identified PRSs explain a low proportion of the variance (1.2-5.3%) in placebo and treatment response. Although the results were limited, we believe that PRS associations bear unredeemed potential as a predictor for treatment response, as more well-powered and phenotypically similar GWAS bases become available.

Potential Anti-Depressive Effects and Mechanisms of Zhi-Zi Hou-Po Decoction Using Behavioral Despair Tests Combined With in Vitro Approaches

Zhi-Zi Hou-Po Decoction (ZHD) has been widely used in the treatment of depression for centuries. This study aimed to investigate the antidepressant effects of the water extract of ZHD (ZHD-WE) and ethanol extract of ZHD (ZHD-EE) using behavioral despair tests in mice, and to further explore the neuroprotective effects in a PC12 cell injury model induced by corticosterone (CORT). Mice were divided into a control group (normal saline), ZHD-WE groups (4, 8, and 16 g kg^-1), ZHD-EE groups (4, 8, and 16 g kg^-1) and the fluoxetine group (20 mg kg^-1). The forced swimming test (FST) and tail suspension test (TST) were used to screen the antidepressant effects of ZHD-WE and ZHD-EE after oral administration for seven consecutive days. The level of brain-derived neurotrophic factor (BDNF) in the hippocampus was determined by ELISA. The MTT, lactate dehydrogenase (LDH) and flow cytometry analysis were performed to elucidate the neuroprotective effect of ZHD-EE on a PC12 cell injury model. Additionally, the mRNA and proteins expression of apoptotic molecules Bax, Bcl-2 and BDNF were detected by RT-PCR and western blot assay. It showed that ZHD-EE at concentrations of 8 and 16 g kg^-1 significantly decreased the immobility time in the TST and FST, and increased the BDNF levels in the hippocampus. While ZHD-WE at concentrations of 4, 8, and 16 g kg^-1 had no significant effect on the immobility time in the TST, and only the 16 g kg^-1 of extract group significantly decreased the immobility time in the FST. In vitro, the obtained results showed that PC12 cells pre-incubated with ZHD-EE at concentrations of 100 and 400 μg ml^-1 improved cell viability, decreased LDH release, and reduced apoptosis rate of PC12 cells. Moreover, ZHD-EE significantly increased the mRNA and proteins expression of Bcl-2 and BDNF, while decreased the mRNA and protein expression of Bax. ZHD-EE significantly improved despair-like behavior in mice, and its mechanism may be related to BDNF upregulation in the hippocampus. This study also showed that ZHD-EE had a protective effect on CORT-induced injury in PC12 cells by upregulating the expression of BDNF and restoring Bcl-2/Bax balance.

The pharmacogenomics of selective serotonin reuptake inhibitors

Antidepressant medications are frequently used as the first line of treatment for depression. However, their effectiveness is highly variable and influenced by genetic factors. Recently, pharmacogenetic studies, including candidate-gene, genome-wide association studies or polygenic risk scores, have attempted to uncover the genetic architecture of antidepressant response. Genetic variants in at least 27 genes are linked to antidepressant treatment response in both coding and non-coding genomic regions, but evidence is largely inconclusive due to the high polygenicity of the trait and limited cohort sizes in published studies. Future studies should increase the number and diversity of participants to yield sufficient statistical power to characterize the genetic underpinnings and biological mechanisms of treatment response, improve results generalizability and reduce racial health-related inequities.

The Potential of Polygenic Risk Scores to Predict Antidepressant Treatment Response in Major Depression: A Systematic Review

BACKGROUND

Understanding the genetic underpinnings of antidepressant treatment response in unipolar major depressive disorder (MDD) can be useful in identifying patients at risk for poor treatment response or treatment resistant depression. A polygenic risk score (PRS) is a useful tool to explore genetic liability of a complex trait such as antidepressant treatment response. Here, we review studies that use PRSs to examine genetic overlap between any trait and antidepressant treatment response in unipolar MDD.

METHODS

A systematic search of literature was conducted in PubMed, Embase, and PsycINFO. Our search included studies examining associations between PRSs of psychiatric as well as non-psychiatric traits and antidepressant treatment response in patients with unipolar MDD. A quality assessment of the included studies was performed.

RESULTS

In total, eleven articles were included which contained PRSs for 30 traits. Studies varied in sample size and endpoints used for antidepressant treatment response. Overall, PRSs for attention-deficit hyperactivity disorder, the personality trait openness, coronary artery disease, obesity, and stroke have been associated with antidepressant treatment response in patients with unipolar MDD.

LIMITATIONS

The endpoints used by included studies differed significantly, therefore it was not possible to perform a meta-analysis.

CONCLUSIONS

Associations between a PRS and antidepressant treatment response have been reported for a number of traits in patients with unipolar MDD. PRSs could be informative to predict antidepressant treatment response in this population, given advances in the field. Most importantly, there is a need for larger study cohorts and the use of standardized outcome measures.

Identifying the Common Genetic Basis of Antidepressant Response

Background

Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction.

Methods

Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA.

Results

Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response.

Conclusions

This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

A meta-analysis of polygenic risk scores for mood disorders, neuroticism, and schizophrenia in antidepressant response

About two-thirds of patients with major depressive disorder (MDD) fail to achieve symptom remission after the initial antidepressant treatment. Despite a role of genetic factors was proven, the specific underpinnings are not fully understood yet. Polygenic risk scores (PRSs), which summarise the additive effect of multiple risk variants across the genome, might provide insights into the underlying genetics. This study aims to investigate the possible association of PRSs for bipolar disorder, MDD, neuroticism, and schizophrenia (SCZ) with antidepressant non-response or non-remission in patients with MDD. PRSs were calculated at eight genome-wide P-thresholds based on publicly available summary statistics of the largest genome-wide association studies. Logistic regressions were performed between PRSs and non-response or non-remission in six European clinical samples, adjusting for age, sex, baseline symptom severity, recruitment sites, and population stratification. Results were meta-analysed across samples, including up to 3,637 individuals. Bonferroni correction was applied. In the meta-analysis, no result was significant after Bonferroni correction. The top result was found for MDD-PRS and non-remission (p = 0.004), with patients in the highest vs. lowest PRS quintile being more likely not to achieve remission (OR=1.5, 95% CI=1.11-1.98, p = 0.007). Nominal associations were also found between MDD-PRS and non-response (p = 0.013), as well as between SCZ-PRS and non-remission (p = 0.035). Although PRSs are still not able to predict non-response or non-remission, our results are in line with previous works; methodological improvements in PRSs calculation may improve their predictive performance and have a meaningful role in precision psychiatry.

Understanding genetic risk factors for common side effects of antidepressant medications

Background

Major depression is one of the most disabling health conditions internationally. In recent years, new generation antidepressant medicines have become very widely prescribed. While these medicines are efficacious, side effects are common and frequently result in discontinuation of treatment. Compared with specific pharmacological properties of the different medications, the relevance of individual vulnerability is understudied.

Methods

We used data from the Australian Genetics of Depression Study to gain insights into the aetiology and genetic risk factors to antidepressant side effects. To this end, we employed structural equation modelling, polygenic risk scoring and regressions.

Results

Here we show that participants reporting a specific side effect for one antidepressant are more likely to report the same side effect for other antidepressants, suggesting the presence of shared individual or pharmacological factors. Polygenic risk scores (PRS) for depression associated with side effects that overlapped with depressive symptoms, including suicidality and anxiety. Body Mass Index PRS are strongly associated with weight gain from all medications. PRS for headaches are associated with headaches from sertraline. Insomnia PRS show some evidence of predicting insomnia from amitriptyline and escitalopram.

Conclusions

Our results suggest a set of common factors underlying the risk for antidepressant side effects. These factors seem to be partly explained by genetic liability related to depression severity and the nature of the side effect. Future studies on the genetic aetiology of side effects will enable insights into their underlying mechanisms and the possibility of risk stratification and prophylaxis strategies.

Antidepressants are commonly prescribed medications, but adverse side effects are cause for treatment discontinuation. We analysed data from a large group of adults who have taken antidepressants to understand why some people experience specific side effects. Our results suggest that a person’s genetic characteristics play a role. For example, participants genetically predisposed to a higher body mass index were more likely to report weight gain from antidepressants. These results open up the possibility of predicting adverse side effects as we increase our knowledge on the genetics of related complex traits. Future studies can focus on performing large-scale genetic studies of antidepressant side effects to gain further insights into the mechanisms underlying antidepressant side effects and to identify genetic markers of side effects that could be used in the clinic.

Campos et al. study the genetic aetiology of antidepressant side effects. Using data from the Australian Genetics of Depression study, the authors show that polygenic risk scores for traits such as BMI, insomnia and headaches have a shared genetic basis with side effects to commonly used antidepressant drugs.

Genetic contributions to anxiety disorders: where we are and where we are heading

Anxiety disorders are among the most common psychiatric disorders worldwide. They often onset early in life, with symptoms and consequences that can persist for decades. This makes anxiety disorders some of the most debilitating and costly disorders of our time. Although much is known about the synaptic and circuit mechanisms of fear and anxiety, research on the underlying genetics has lagged behind that of other psychiatric disorders. However, alongside the formation of the Psychiatric Genomic Consortium Anxiety workgroup, progress is rapidly advancing, offering opportunities for future research.Here we review current knowledge about the genetics of anxiety across the lifespan from genetically informative designs (i.e. twin studies and molecular genetics). We include studies of specific anxiety disorders (e.g. panic disorder, generalised anxiety disorder) as well as those using dimensional measures of trait anxiety. We particularly address findings from large-scale genome-wide association studies and show how such discoveries may provide opportunities for translation into improved or new therapeutics for affected individuals. Finally, we describe how discoveries in anxiety genetics open the door to numerous new research possibilities, such as the investigation of specific gene-environment interactions and the disentangling of causal associations with related traits and disorders.We discuss how the field of anxiety genetics is expected to move forward. In addition to the obvious need for larger sample sizes in genome-wide studies, we highlight the need for studies among young people, focusing on specific underlying dimensional traits or components of anxiety.

Self-reported medication use as an alternate phenotyping method for anxiety and depression in the UK Biobank

The requirement for large sample sizes for psychiatric genetic analyses necessitates novel approaches to derive cases. Anxiety and depression show substantial genetic overlap and share pharmacological treatments. Data on prescribed medication could be effective for inferring case status when other indicators of mental health are unavailable. We investigated self-reported current medication use in UK Biobank participants of European ancestry. Medication Status cases reported using antidepressant or anxiolytic medication (n = 22,218), controls did not report psychotropic medication use (n = 168,959). A subset, "Medication Only," additionally did not meet criteria for any other mental health indicator (case n = 2,643, control n = 107,029). We assessed genetic overlap between these phenotypes and two published genetic association studies of anxiety and depression, and an internalizing disorder trait derived from symptom-based questionnaires in UK Biobank. Genetic correlations between Medication Status and the three anxiety and depression phenotypes were significant (r_g  = 0.60-0.73). In the Medication Only subset, the genetic correlation with depression was significant (r_g  = 0.51). The three polygenic scores explained 0.33% - 0.80% of the variance in Medication Status and 0.07% - 0.19% of the variance in Medication Only. This study provides evidence that self-reported current medication use offers an alternate or supplementary anxiety or depression phenotype in genetic studies where diagnostic information is sparse or unavailable.

Genetic and clinical characteristics of treatment-resistant depression using primary care records in two UK cohorts

Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Primary care electronic health records provide an easily accessible approach to investigate TRD clinical and genetic characteristics. MDD defined from primary care records in UK Biobank (UKB) and EXCEED studies was compared with other measures of depression and tested for association with MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least 6 weeks. Clinical-demographic characteristics, SNP-based heritability (h²_SNP) and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases. In 230,096 and 8926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD, respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 71-88% with other MDD definitions. In UKB, TRD vs healthy controls and non-TRD vs healthy controls h²_SNP was comparable (0.25 [SE = 0.04] and 0.19 [SE = 0.02], respectively). TRD vs non-TRD was positively associated with the PRS of attention deficit hyperactivity disorder, with lower socio-economic status, obesity, higher neuroticism and other unfavourable clinical characteristics. This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.

Higher polygenic risk scores for schizophrenia may be suggestive of treatment non-response in major depressive disorder

Up to 60% of patients with major depressive disorder (MDD) do not respond to the first treatment with antidepressants. Response to antidepressants is a polygenic trait, although its underpinning genetics has not been fully clarified. This study aimed to investigate if polygenic risk scores (PRSs) for major psychiatric disorders and trait neuroticism (NEU) were associated with non-response or resistance to antidepressants in MDD. PRSs for bipolar disorder, MDD, NEU, and schizophrenia (SCZ) were computed in 1,148 patients with MDD. Summary statistics from the largest meta-analyses of genome-wide association studies were used as base data. Patients were classified as responders, non-responders to one treatment, non-responders to two or more treatments (treatment-resistant depression or TRD). Regression analyses were adjusted for population stratification and recruitment sites. PRSs did not predict either non-response vs response or TRD vs response after Bonferroni correction. However, SCZ-PRS was nominally associated with non-response (p = 0.003). Patients in the highest SCZ-PRS quintile were more likely to be non-responders than those in the lowest quintile (OR = 2.23, 95% CI = 1.21-4.10, p = 0.02). Patients in the lowest SCZ-PRS quintile showed higher response rates when they did not receive augmentation with second-generation antipsychotics (SGAs), while those in the highest SCZ-PRS quintile had a poor response independently from the treatment strategy (p = 0.009). A higher genetic liability to SCZ may reduce treatment response in MDD, and patients with low SCZ-PRSs may show higher response rates without SGA augmentation. Multivariate approaches and methodological refinements will be necessary before clinical implementations of PRSs.

The promise of machine learning in predicting treatment outcomes in psychiatry

For many years, psychiatrists have tried to understand factors involved in response to medications or psychotherapies, in order to personalize their treatment choices. There is now a broad and growing interest in the idea that we can develop models to personalize treatment decisions using new statistical approaches from the field of machine learning and applying them to larger volumes of data. In this pursuit, there has been a paradigm shift away from experimental studies to confirm or refute specific hypotheses towards a focus on the overall explanatory power of a predictive model when tested on new, unseen datasets. In this paper, we review key studies using machine learning to predict treatment outcomes in psychiatry, ranging from medications and psychotherapies to digital interventions and neurobiological treatments. Next, we focus on some new sources of data that are being used for the development of predictive models based on machine learning, such as electronic health records, smartphone and social media data, and on the potential utility of data from genetics, electrophysiology, neuroimaging and cognitive testing. Finally, we discuss how far the field has come towards implementing prediction tools in real-world clinical practice. Relatively few retrospective studies to-date include appropriate external validation procedures, and there are even fewer prospective studies testing the clinical feasibility and effectiveness of predictive models. Applications of machine learning in psychiatry face some of the same ethical challenges posed by these techniques in other areas of medicine or computer science, which we discuss here. In short, machine learning is a nascent but important approach to improve the effectiveness of mental health care, and several prospective clinical studies suggest that it may be working already.

Ligand-specific pharmacogenetic effects of nonsynonymous mutations

In pharmacogenomics, variable receptor phenotypes, resulting from genetic polymorphisms, are often described as a change in protein function or regulation observed upon exposure to a drug. However, in some instances, phenotypes are defined using a class of medications rather than individual drugs. This paradigm assumes that a variation associated with a drug response phenotype will retain the magnitude and direction of the effect for other drugs with the same mechanism of action. However, nonsynonymous polymorphisms may have ligand-specific effects. The purpose of this study was to investigate the potential for point mutations to asymmetrically affect the binding of different drugs to a common target. Ligand binding data from site-directed mutagenesis studies on five G-protein coupled receptors (beta-1 and -2 adrenergic, dopamine D2, angiotensin II and mu-opioid receptor) were collected and analyzed. Binding data from 81 studies for 253 ligands with 447 mutant proteins, including 10 naturally occurring human variants, were analyzed, yielding 1989 mutation-ligand pairs. Fold change in binding affinity for mutant proteins, relative to the wild-type, for different drugs was examined for ligand-specific effects, with a fold-change difference of one or more orders of magnitude between agents considered significant. Of the mutations examined, 49% were associated with ligand-specific effects. One human variant (T164I, beta-2 adrenergic receptor) showed ligand-specific effects for antiasthmatic agents. These results indicate that ligand-specific changes in binding are a possible consequence of missense mutations. This implies that caution needs to be exercised when grouping drugs together during design or interpretation of genotype-phenotype association studies.

Prospects for personalization of depression treatment with genome sequencing

The effectiveness of antidepressants in the treatment of major depressive disorder varies considerably between patients. With these interindividual differences and a number of antidepressants to choose from, the first choice of treatment often fails to produce improvement in the patient's condition. A substantial part of the variation in response to antidepressants can be explained by genetic factors. Accordingly, variants related to drug metabolism in two pharmacogenes, CYP2D6 and CYP2C19, have already been translated into guidelines for antidepressant prescriptions. The role of variants in other genes that influence antidepressant responses is not yet understood. Furthermore, rare and individual variants account for a substantial part of genetic differences in antidepressant efficacy. Recent years have brought a tremendous increase in the accessibility of genome sequencing in terms of data availability and its clinical use. In this review, we summarize recent developments and current issues in the personalization of major depressive disorder treatment through pharmacogenomics.

A distinct transcriptional signature of antidepressant response in hippocampal dentate gyrus granule cells

Major depressive disorder is the most prevalent mental illness worldwide, still its pharmacological treatment is limited by various challenges, such as the large heterogeneity in treatment response and the lack of insight into the neurobiological pathways underlying this phenomenon. To decode the molecular mechanisms shaping antidepressant response and to distinguish those from general paroxetine effects, we used a previously established approach targeting extremes (i.e., good vs poor responder mice). We focused on the dentate gyrus (DG), a subregion of major interest in the context of antidepressant mechanisms. Transcriptome profiling on micro-dissected DG granule cells was performed to (i) reveal cell-type-specific changes in paroxetine-induced gene expression (paroxetine vs vehicle) and (ii) to identify molecular signatures of treatment response within a cohort of paroxetine-treated animals. We identified 112 differentially expressed genes associated with paroxetine treatment. The extreme group comparison (good vs poor responder) yielded 211 differentially expressed genes. General paroxetine effects could be distinguished from treatment response-associated molecular signatures, with a differential gene expression overlap of only 4.6% (15 genes). Biological pathway enrichment and cluster analyses identified candidate mechanisms associated with good treatment response, e.g., neuropeptide signaling, synaptic transmission, calcium signaling, and regulation of glucocorticoid secretion. Finally, we examined glucocorticoid receptor (GR)-dependent regulation of selected response-associated genes to analyze a hypothesized interplay between GR signaling and good antidepressant treatment response. Among the most promising candidates, we suggest potential targets such as the developmental gene Otx2 or Htr2c for further investigations into antidepressant treatment response in the future.

Exploring the Role of Contactins across Psychological, Psychiatric and Cardiometabolic Traits within UK Biobank

Individuals with severe mental illness have an increased risk of cardiometabolic diseases compared to the general population. Shared risk factors and medication effects explain part of this excess risk; however, there is growing evidence to suggest that shared biology (including genetic variation) is likely to contribute to comorbidity between mental and physical illness. Contactins are a family of genes involved in development of the nervous system and implicated, though genome-wide association studies, in a wide range of psychological, psychiatric and cardiometabolic conditions. Contactins are plausible candidates for shared pathology between mental and physical health. We used data from UK Biobank to systematically assess how genetic variation in contactin genes was associated with a wide range of psychological, psychiatric and cardiometabolic conditions. We also investigated whether associations for cardiometabolic and psychological traits represented the same or distinct signals and how the genetic variation might influence the measured traits. We identified: A novel genetic association between variation in CNTN1 and current smoking; two independent signals in CNTN4 for BMI; and demonstrated that associations between CNTN5 and neuroticism were distinct from those between CNTN5 and blood pressure/HbA1c. There was no evidence that the contactin genes contributed to shared aetiology between physical and mental illness.

The influence of common polygenic risk and gene sets on social skills group training response in autism spectrum disorder

Social skills group training (SSGT) is a frequently used behavioral intervention in autism spectrum disorder (ASD), but the effects are moderate and heterogeneous. Here, we analyzed the effect of polygenic risk score (PRS) and common variants in gene sets on the intervention outcome. Participants from the largest randomized clinical trial of SSGT in ASD to date were selected (N = 188, 99 from SSGT, 89 from standard care) to calculate association between the outcomes in the SSGT trial and PRSs for ASD, attention-deficit hyperactivity disorder (ADHD), and educational attainment. In addition, specific gene sets were selected to evaluate their role on intervention outcomes. Among all participants in the trial, higher PRS for ADHD was associated with significant improvement in the outcome measure, the parental-rated Social Responsiveness Scale. The significant association was due to better outcomes in the standard care group for individuals with higher PRS for ADHD (post-intervention: β = −4.747, P = 0.0129; follow-up: β = −5.309, P = 0.0083). However, when contrasting the SSGT and standard care group, an inferior outcome in the SSGT group was associated with higher ADHD PRS at follow-up (β = 6.67, P = 0.016). Five gene sets within the synaptic category showed a nominal association with reduced response to interventions. We provide preliminary evidence that genetic liability calculated from common variants could influence the intervention outcomes. In the future, larger cohorts should be used to investigate how genetic contribution affects individual response to ASD interventions.

Reviewing the genetics of heterogeneity in depression: operationalizations, manifestations and etiologies

With progress in genome-wide association studies of depression, from identifying zero hits in ~16 000 individuals in 2013 to 223 hits in more than a million individuals in 2020, understanding the genetic architecture of this debilitating condition no longer appears to be an impossible task. The pressing question now is whether recently discovered variants describe the etiology of a single disease entity. There are a myriad of ways to measure and operationalize depression severity, and major depressive disorder as defined in the Diagnostic and Statistical Manual of Mental Disorders-5 can manifest in more than 10 000 ways based on symptom profiles alone. Variations in developmental timing, comorbidity and environmental contexts across individuals and samples further add to the heterogeneity. With big data increasingly enabling genomic discovery in psychiatry, it is more timely than ever to explicitly disentangle genetic contributions to what is likely 'depressions' rather than depression. Here, we introduce three sources of heterogeneity: operationalization, manifestation and etiology. We review recent efforts to identify depression subtypes using clinical and data-driven approaches, examine differences in genetic architecture of depression across contexts, and argue that heterogeneity in operationalizations of depression is likely a considerable source of inconsistency. Finally, we offer recommendations and considerations for the field going forward.

Investigating an in silico approach for prioritizing antidepressant drug prescription based on drug-induced expression profiles and predicted gene expression

In clinical practice, an antidepressant prescription is a trial and error approach, which is time consuming and discomforting for patients. This study investigated an in silico approach for ranking antidepressants based on their hypothetical likelihood of efficacy. We predicted the transcriptomic profile of citalopram remitters by performing an in silico transcriptomic-wide association study on STAR*D GWAS data (N = 1163). The transcriptional profile of remitters was compared with 21 antidepressant-induced gene expression profiles in five human cell lines available in the connectivity-map database. Spearman correlation, Pearson correlation, and the Kolmogorov-Smirnov test were used to determine the similarity between antidepressant-induced profiles and remitter profiles, subsequently calculating the average rank of antidepressants across the three methods and a p value for each rank by using a permutation procedure. The drugs with the top ranks were those having a high positive correlation with the expression profiles of remitters and that may have higher chances of efficacy in the tested patients. In MCF7 (breast cancer cell line), escitalopram had the highest average rank, with an average rank higher than expected by chance (p = 0.0014). In A375 (human melanoma) and PC3 (prostate cancer) cell lines, escitalopram and citalopram emerged as the second-highest ranked antidepressants, respectively (p = 0.0310 and 0.0276, respectively). In HA1E (kidney) and HT29 (colon cancer) cell types, citalopram and escitalopram did not fall among top antidepressants. The correlation between citalopram remitters' and (es)citalopram-induced expression profiles in three cell lines suggests that our approach may be useful and with future improvements, it can be applicable at the individual level to tailor treatment prescription.

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.

Genetic stratification of depression in UK Biobank

Depression is a common and clinically heterogeneous mental health disorder that is frequently comorbid with other diseases and conditions. Stratification of depression may align sub-diagnoses more closely with their underling aetiology and provide more tractable targets for research and effective treatment. In the current study, we investigated whether genetic data could be used to identify subgroups within people with depression using the UK Biobank. Examination of cross-locus correlations were used to test for evidence of subgroups using genetic data from seven other complex traits and disorders that were genetically correlated with depression and had sufficient power (>0.6) for detection. We found no evidence for subgroups within depression for schizophrenia, bipolar disorder, attention deficit/hyperactivity disorder, autism spectrum disorder, anorexia nervosa, inflammatory bowel disease or obesity. This suggests that for these traits, genetic correlations with depression were driven by pleiotropic genetic variants carried by everyone rather than by a specific subgroup.

Widespread Morphometric Abnormalities in Major Depression: Neuroplasticity and Potential for Biomarker Development

Major depression is common and debilitating. Identifying neurobiological subtypes that comprise the disorder and predict clinical outcome are key challenges. Genetic and environmental factors leading to major depression are expressed in neural structure and function. Volumetric decreases in gray matter have been demonstrated in corticolimbic circuits involved in emotion regulation. MR imaging observable abnormalities reflect cytoarchitectonic alterations within a local neuroendocrine milieu with systemic effects. Multivariate pattern analysis offers the potential to identify the neurobiological subtypes and predictors of clinical outcome. It is essential to characterize disease heterogeneity by incorporating data-driven inductive and symptom-based deductive approaches in an iterative process.

Pharmacogenetics and Depression: A Critical Perspective

Depression leads the higher personal and socio-economical burden within psychiatric disorders. Despite the fact that over 40 antidepressants (ADs) are available, suboptimal response still poses a major challenge and is thought to be partially a result of genetic variation. Pharmacogenetics studies the effects of genetic variants on treatment outcomes with the aim of providing tailored treatments, thereby maximizing efficacy and tolerability. After two decades of pharmacogenetic research, variants in genes coding for the cytochromes involved in ADs metabolism (CYP2D6 and CYP2C19) are now considered biomarkers with sufficient scientific support for clinical application, despite the lack of conclusive cost/effectiveness evidence. The effect of variants in genes modulating ADs mechanisms of action (pharmacodynamics) is still controversial, because of the much higher complexity of ADs pharmacodynamics compared to ADs metabolism. Considerable progress has been made since the era of candidate gene studies: the genomic revolution has made possible to assess genetic variance on an unprecedented scale, throughout the whole genome, and to analyze the cumulative effect of different variants. The results have revealed key information on the biological mechanisms mediating ADs effect and identified hypothetical new pharmacological targets. They also paved the way for future availability of polygenic pharmacogenetic panels to predict treatment outcome, which are expected to explain much higher variance in ADs response compared to CYP2D6 and CYP2C19 only. As the demand and availability of AD pharmacogenetic testing is projected to increase, it is important for clinicians to keep abreast of this evolving area to facilitate informed discussions with their patients.

Machine learning methods for developing precision treatment rules with observational data

Clinical trials have identified a variety of predictor variables for use in precision treatment protocols, ranging from clinical biomarkers and symptom profiles to self-report measures of various sorts. Although such variables are informative collectively, none has proven sufficiently powerful to guide optimal treatment selection individually. This has prompted growing interest in the development of composite precision treatment rules (PTRs) that are constructed by combining information across a range of predictors. But this work has been hampered by the generally small samples in randomized clinical trials and the use of suboptimal analysis methods to analyze the resulting data. In this paper, we propose to address the sample size problem by: working with large observational electronic medical record databases rather than controlled clinical trials to develop preliminary PTRs; validating these preliminary PTRs in subsequent pragmatic trials; and using ensemble machine learning methods rather than individual algorithms to carry out statistical analyses to develop the PTRs. The major challenges in this proposed approach are that treatment are not randomly assigned in observational databases and that these databases often lack measures of key prescriptive predictors and mental disorder treatment outcomes. We proposed a tiered case-cohort design approach that uses innovative methods for measuring and balancing baseline covariates and estimating PTRs to address these challenges.

A genome-wide association meta-analysis of prognostic outcomes following cognitive behavioural therapy in individuals with anxiety and depressive disorders

Major depressive disorder and the anxiety disorders are highly prevalent, disabling and moderately heritable. Depression and anxiety are also highly comorbid and have a strong genetic correlation (rg ≈ 1). Cognitive behavioural therapy is a leading evidence-based treatment but has variable outcomes. Currently, there are no strong predictors of outcome. Therapygenetics research aims to identify genetic predictors of prognosis following therapy. We performed genome-wide association meta-analyses of symptoms following cognitive behavioural therapy in adults with anxiety disorders (n = 972), adults with major depressive disorder (n = 832) and children with anxiety disorders (n = 920; meta-analysis n = 2724). We estimated the variance in therapy outcomes that could be explained by common genetic variants (h2SNP) and polygenic scoring was used to examine genetic associations between therapy outcomes and psychopathology, personality and learning. No single nucleotide polymorphisms were strongly associated with treatment outcomes. No significant estimate of h2SNP could be obtained, suggesting the heritability of therapy outcome is smaller than our analysis was powered to detect. Polygenic scoring failed to detect genetic overlap between therapy outcome and psychopathology, personality or learning. This study is the largest therapygenetics study to date. Results are consistent with previous, similarly powered genome-wide association studies of complex traits.

Uncovering the Genetic Architecture of Major Depression

There have been several recent studies addressing the genetic architecture of depression. This review serves to take stock of what is known now about the genetics of depression, how it has increased our knowledge and understanding of its mechanisms, and how the information and knowledge can be leveraged to improve the care of people affected. We identify four priorities for how the field of MD genetics research may move forward in future years, namely by increasing the sample sizes available for genome-wide association studies (GWASs), greater inclusion of diverse ancestries and low-income countries, the closer integration of psychiatric genetics with electronic medical records, and the development of the neuroscience toolkit for polygenic disorders.