Clinical and genetic contributions to medical comorbidity in bipolar disorder: a study using electronic health records-linked biobank data.

Bipolar disorder is a chronic and complex polygenic disease with high rates of comorbidity. However, the independent contribution of either diagnosis or genetic risk of bipolar disorder to the medical comorbidity profile of individuals with the disease remains unresolved. Here, we conducted a multi-step phenome-wide association study (PheWAS) of bipolar disorder using phenomes derived from the electronic health records of participants enrolled in the Mayo Clinic Biobank and the Mayo Clinic Bipolar Disorder Biobank. First, we explored the conditions associated with a diagnosis of bipolar disorder by conducting a phenotype-based PheWAS followed by LASSO-penalized regression to account for correlations within the phenome. Then, we explored the conditions associated with bipolar disorder polygenic risk score (BD-PRS) using a PRS-based PheWAS with a sequential exclusion approach to account for the possibility that diagnosis, instead of genetic risk, may drive such associations. 53,386 participants (58.7% women) with a mean age at analysis of 67.8 years (SD = 15.6) were included. A bipolar disorder diagnosis (n = 1479) was associated with higher rates of psychiatric conditions, injuries and poisonings, endocrine/metabolic and neurological conditions, viral hepatitis C, and asthma. BD-PRS was associated with psychiatric comorbidities but, in contrast, had no positive associations with general medical conditions. While our findings warrant confirmation with longitudinal-prospective studies, the limited associations between bipolar disorder genetics and medical conditions suggest that shared environmental effects or environmental consequences of diagnosis may have a greater impact on the general medical comorbidity profile of individuals with bipolar disorder than its genetic risk.

IL17RB genetic variants are associated with acamprosate treatment response in patients with alcohol use disorder: A proteomics-informed genomics study.

Acamprosate is a Food and Drug Administration (FDA) approved medication for the treatment of alcohol use disorder (AUD). However, only a subset of patients achieves optimal treatment outcomes. Currently, no biological measures are utilized to predict response to acamprosate treatment. We applied our established pharmaco-omics informed genomics strategy to identify potential biomarkers associated with acamprosate treatment response. Specifically, our previous open-label acamprosate clinical trial recruited 442 patients with AUD who were treated with acamprosate for three months. We first performed proteomics using baseline plasma samples to identify potential biomarkers associated with acamprosate treatment outcomes. Next, we applied our established "proteomics-informed genome-wide association study (GWAS)" research strategy, and identified 12 proteins, including interleukin-17 receptor B (IL17RB), associated with acamprosate treatment response.​ A GWAS for IL17RB concentrations identified several genome-wide significant signals. Specifically, the top hit single nucleotide polymorphism (SNP) rs6801605 with a minor allele frequency of 38% in the European American population mapped 4 kilobase (Kb) upstream of IL17RB, and intron 1 of the choline dehydrogenase (CHDH) gene on chromosome 3 (p: 4.8E-20). The variant genotype (AA) for the SNP rs6801605 was associated with lower IL17RB protein expression. In addition, we identified a series of genetic variants in IL17RB that were associated with acamprosate treatment outcomes. Furthermore, the variantgenotypes for all of those IL17RB SNPs were protective for alcohol relapse. Finally, we demonstrated that the basal level of mRNA expression of IL17RB was inversely correlated with those of nuclear factor-κB (NF-κB) subunits, and a significantly higher expression of NF-κB subunits was observed in AUD patients who relapsed to alcohol use. In summary, this study illustrates that IL17RB genetic variants might contribute to acamprosate treatment outcomes. This series of studies represents an important step toward generating functional hypotheses that could be tested to gain insight into mechanisms underlying acamprosate treatment response phenotypes. (The ClinicalTrials.gov Identifier: NCT00662571).

Antidepressants that increase mitochondrial energetics may elevate risk of treatment-emergent mania.

Preclinical evidence suggests that antidepressants (ADs) may differentially influence mitochondrial energetics. This study was conducted to investigate the relationship between mitochondrial function and illness vulnerability in bipolar disorder (BD), specifically risk of treatment-emergent mania (TEM). Participants with BD already clinically phenotyped as TEM+ (n = 176) or TEM- (n = 516) were further classified whether the TEM associated AD, based on preclinical studies, increased (Mito+, n = 600) or decreased (Mito-, n = 289) mitochondrial electron transport chain (ETC) activity. Comparison of TEM+ rates between Mito+ and Mito- ADs was performed using generalized estimating equations to account for participants exposed to multiple ADs while adjusting for sex, age at time of enrollment into the biobank and BD type (BD-I/schizoaffective vs. BD-II). A total of 692 subjects (62.7% female, 91.4% White, mean age 43.0 ± 14.0 years) including 176 cases (25.3%) of TEM+ and 516 cases (74.7%) of TEM- with previous exposure to Mito+ and/or Mito- antidepressants were identified. Adjusting for age, sex and BD subtype, TEM+ was more frequent with antidepressants that increased (24.7%), versus decreased (13.5%) mitochondrial energetics (OR = 2.21; p = 0.000009). Our preliminary retrospective data suggests there may be merit in reconceptualizing AD classification, not solely based on monoaminergic conventional drug mechanism of action, but additionally based on mitochondrial energetics. Future prospective clinical studies on specific antidepressants and mitochondrial activity are encouraged. Recognizing pharmacogenomic investigation of drug response may extend or overlap to genomics of disease risk, future studies should investigate potential interactions between mitochondrial mechanisms of disease risk and drug response.

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.

Patterns of pharmacotherapy for bipolar disorder: A GBC survey.

OBJECTIVES: To understand treatment practices for bipolar disorders (BD), this study leveraged the Global Bipolar Cohort collaborative network to investigate pharmacotherapeutic treatment patterns in multiple cohorts of well-characterized individuals with BD in North America, Europe, and Australia. METHODS: Data on pharmacotherapy, demographics, diagnostic subtypes, and comorbidities were provided from each participating cohort. Individual site and regional pooled proportional meta-analyses with generalized linear mixed methods were conducted to identify prescription patterns. RESULTS: This study included 10,351 individuals from North America (n = 3985), Europe (n = 3822), and Australia (n = 2544). Overall, participants were predominantly female (60%) with BD-I (60%; vs. BD-II = 33%). Cross-sectionally, mood-stabilizing anticonvulsants (44%), second-generation antipsychotics (42%), and antidepressants (38%) were the most prescribed medications. Lithium was prescribed in 29% of patients, primarily in the Australian (31%) and European (36%) cohorts. First-generation antipsychotics were prescribed in 24% of the European versus 1% in the North American cohort. Antidepressant prescription rates were higher in BD-II (47%) compared to BD-I (35%). Major limitations were significant differences among cohorts based on inclusion/exclusion criteria, data source, and time/year of enrollment into cohort. CONCLUSIONS: Mood-stabilizing anticonvulsants, second-generation antipsychotics, and antidepressants were the most prescribed medications suggesting prescription patterns that are not necessarily guideline concordant. Significant differences exist in the prescription practices across different geographic regions, especially the underutilization of lithium in the North American cohorts and the higher utilization of first-generation antipsychotics in the European cohorts. There is a need to conduct future longitudinal studies to further explore these differences and their impact on outcomes, and to inform and implement evidence-based guidelines to help improve treatment practices in BD.

Phenotypic distinctions in depression and anxiety: a comparative analysis of comorbid and isolated cases.

BACKGROUND: Anxiety and depression are frequently comorbid yet phenotypically distinct. This study identifies differences in the clinically observable phenome across a wide variety of physical and mental disorders comparing patients with diagnoses of depression without anxiety, anxiety without depression, or both depression and anxiety. METHODS: Using electronic health records for 14 994 participants with depression and/or anxiety in the Mayo Clinic Biobank, a phenotype-based phenome-wide association study (Phe2WAS) was performed to test for differences between these groups across a broad range of clinical diagnoses observed in the electronic health record. Additional analyses were performed to determine the temporal sequencing of diagnoses. RESULTS: Compared to patients diagnosed only with anxiety, those diagnosed only with depression were more likely to have diagnoses of obesity (OR 1.75; p = 1 × 10-27), sleep apnea (OR 1.71; p = 1 × 10-22), and type II diabetes (OR 1.74; p = 9 × 10-18). Compared to those diagnosed only with depression, those diagnosed only with anxiety were more likely to have diagnoses of palpitations (OR 1.91; p = 2 × 10-25), benign skin neoplasms (OR 1.61; p = 2 × 10-17), and cardiac dysrhythmias (OR 1.45; p = 2 × 10-12). Patients with comorbid depression and anxiety were more likely to have diagnoses of other mental health disorders, substance use disorders, sleep problems, and gastroesophageal reflux relative to isolated depression. CONCLUSIONS: While depression and anxiety are closely related, this study suggests that phenotypic distinctions exist between depression and anxiety. Improving phenotypic characterization within the broad categories of depression and anxiety could improve the clinical assessment of depression and anxiety.

The importance of social activity to risk of major depression in older adults.

BACKGROUND: Several social determinants of health (SDoH) have been associated with the onset of major depressive disorder (MDD). However, prior studies largely focused on individual SDoH and thus less is known about the relative importance (RI) of SDoH variables, especially in older adults. Given that risk factors for MDD may differ across the lifespan, we aimed to identify the SDoH that was most strongly related to newly diagnosed MDD in a cohort of older adults. METHODS: We used self-reported health-related survey data from 41 174 older adults (50-89 years, median age = 67 years) who participated in the Mayo Clinic Biobank, and linked ICD codes for MDD in the participants' electronic health records. Participants with a history of clinically documented or self-reported MDD prior to survey completion were excluded from analysis (N = 10 938, 27%). We used Cox proportional hazards models with a gradient boosting machine approach to quantify the RI of 30 pre-selected SDoH variables on the risk of future MDD diagnosis. RESULTS: Following biobank enrollment, 2073 older participants were diagnosed with MDD during the follow-up period (median duration = 6.7 years). The most influential SDoH was perceived level of social activity (RI = 0.17). Lower level of social activity was associated with a higher risk of MDD [hazard ratio = 2.27 (95% CI 2.00-2.50) for highest v. lowest level]. CONCLUSION: Across a range of SDoH variables, perceived level of social activity is most strongly related to MDD in older adults. Monitoring changes in the level of social activity may help identify older adults at an increased risk of MDD.

The genetic contribution to the comorbidity of depression and anxiety: a multi-site electronic health records study of almost 178 000 people.

BACKGROUND: Depression and anxiety are common and highly comorbid, and their comorbidity is associated with poorer outcomes posing clinical and public health concerns. We evaluated the polygenic contribution to comorbid depression and anxiety, and to each in isolation. METHODS: Diagnostic codes were extracted from electronic health records for four biobanks [N = 177 865 including 138 632 European (77.9%), 25 612 African (14.4%), and 13 621 Hispanic (7.7%) ancestry participants]. The outcome was a four-level variable representing the depression/anxiety diagnosis group: neither, depression-only, anxiety-only, and comorbid. Multinomial regression was used to test for association of depression and anxiety polygenic risk scores (PRSs) with the outcome while adjusting for principal components of ancestry. RESULTS: In total, 132 960 patients had neither diagnosis (74.8%), 16 092 depression-only (9.0%), 13 098 anxiety-only (7.4%), and 16 584 comorbid (9.3%). In the European meta-analysis across biobanks, both PRSs were higher in each diagnosis group compared to controls. Notably, depression-PRS (OR 1.20 per s.d. increase in PRS; 95% CI 1.18-1.23) and anxiety-PRS (OR 1.07; 95% CI 1.05-1.09) had the largest effect when the comorbid group was compared with controls. Furthermore, the depression-PRS was significantly higher in the comorbid group than the depression-only group (OR 1.09; 95% CI 1.06-1.12) and the anxiety-only group (OR 1.15; 95% CI 1.11-1.19) and was significantly higher in the depression-only group than the anxiety-only group (OR 1.06; 95% CI 1.02-1.09), showing a genetic risk gradient across the conditions and the comorbidity. CONCLUSIONS: This study suggests that depression and anxiety have partially independent genetic liabilities and the genetic vulnerabilities to depression and anxiety make distinct contributions to comorbid depression and anxiety.

Antidepressant-Associated Treatment Emergent Mania: A Meta-Analysis to Guide Risk Modeling Pharmacogenomic Targets of Potential Clinical Value.

BACKGROUND: The purpose of this study was to review the association between the SLC6A4 5-HTTLPR polymorphism and antidepressant (AD)-associated treatment emergent mania (TEM) in bipolar disorder alongside starting a discussion on the merits of developing risk stratification models to guide when not to provide AD treatment for bipolar depression. METHODS: Studies that examined the association between clinical and genetic risk factors, specifically monoaminergic transporter genetic variation, and TEM were identified. A meta-analysis was performed using the odds ratio to estimate the effect size under the Der-Simonian and Laird model. RESULTS: Seven studies, referencing the SLC6A4 5-HTTLPR polymorphism and TEM (total N = 1578; TEM+ =594, TEM- = 984), of 142 identified articles were included. The time duration between the start of the AD to emergence of TEM ranged from 4 to 12 weeks. There was a nominally significant association between the s allele of the 5-HTTLPR polymorphism and TEM (odds ratio, 1.434; 95% confidence interval, 1.001-2.055; P = 0.0493; I2 = 52%). No studies have investigated norepinephrine or dopamine transporters. CONCLUSION: Although the serotonin transporter genetic variation is commercially available in pharmacogenomic decision support tools, greater efforts, more broadly, should focus on complete genome-wide approaches to determine genetic variants that may contribute to TEM. Moreover, these data are exemplary to the merits of developing risk stratification models, which include both clinical and biological risk factors, to guide when not to use ADs in bipolar disorder. Future studies will need to validate new risk models that best inform the development of personalized medicine best practices treating bipolar depression.

Real-World Clinical Practice Among Patients With Bipolar Disorder and Chronic Kidney Disease on Long-term Lithium Therapy.

PURPOSE: Long-term lithium therapy (LTLT) has been associated with chronic kidney disease (CKD). We investigated changes in clinical characteristics, pharmacotherapeutic treatments for medical/psychiatric disorders, and outcomes among patients with bipolar disorder (BD) and CKD on LTLT in a 2-year mirror-image study design. METHODS: Adult BD patients on LTLT for ≥1 year who enrolled in the Mayo Clinic Bipolar Disorder Biobank and developed CKD (stage 3) were included, and our study was approved by the Mayo Clinic Institutional Review Board. The primary outcome was the time to the first mood episode after CKD diagnosis among the lithium (Li) continuers and discontinuers. Cox proportional hazards models were used to estimate the time to the first mood episode. We tested for differences in other medication changes between the Li continuers and discontinuers group using Mantel-Haenszel χ2 tests (linear associations). RESULTS: Of 38 BD patients who developed CKD, 18 (47%) discontinued Li, and the remainder continued (n = 20). The median age of the cohort was 56 years (interquartile range [IQR], 48-67 years), 63.2% were female, and 97.4% were White. As compared with continuers, discontinuers had more psychotropic medication trials (6 [IQR, 4-6] vs 3 [IQR, 2-5], P = 0.02), a higher rate of 1 or more mood episodes (61% vs 10%, P = 0.002), and a higher risk of a mood episode after CKD diagnoses (Hazard Ratio, 8.38; 95% confidence interval, 1.85-38.0 [log-rank P = 0.001]]. CONCLUSIONS: Bipolar disorder patients on LTLT who discontinued Li had a higher risk for relapse and a shorter time to the first mood episode, suggesting a need for more thorough discussion before Li discontinuation after the CKD diagnosis.

Single cell transcriptomics reveals distinct transcriptional responses to oxycodone and buprenorphine by iPSC-derived brain organoids from patients with opioid use disorder.

The opioid epidemic represents a national crisis. Oxycodone is one of the most prescribed opioid medications in the United States, whereas buprenorphine is currently the most prescribed medication for opioid use disorder (OUD) pharmacotherapy. Given the extensive use of prescription opioids and the global opioid epidemic, it is essential to understand how opioids modulate brain cell type function at the single-cell level. We performed single nucleus RNA-seq (snRNA-seq) using iPSC-derived forebrain organoids from three male OUD subjects in response to oxycodone, buprenorphine, or vehicle for seven days. We utilized the snRNA-seq data to identify differentially expressed genes following drug treatment using the Seurat integrative analysis pipeline. We utilized iPSC-derived forebrain organoids and single-cell sequencing technology as an unbiased tool to study cell-type-specific and drug-specific transcriptional responses. After quality control filtering, we analyzed 25787 cells and identified sixteen clusters using unsupervised clustering analysis. Our results reveal distinct transcriptional responses to oxycodone and buprenorphine by iPSC-derived brain organoids from patients with OUD. Specifically, buprenorphine displayed a significant influence on transcription regulation in glial cells. However, oxycodone induced type I interferon signaling in many cell types, including neural cells in brain organoids. Finally, we demonstrate that oxycodone, but not buprenorphine activated STAT1 and induced the type I interferon signaling in patients with OUD. These data suggest that elevation of STAT1 expression associated with OUD might play a role in transcriptional regulation in response to oxycodone. In summary, our results provide novel mechanistic insight into drug action at single-cell resolution.

Shared brain transcriptomic signature in TDP-43 type A FTLD patients with or without GRN mutations.

Frontotemporal lobar degeneration with TDP-43 inclusions (FTLD-TDP) is a complex heterogeneous neurodegenerative disorder for which mechanisms are poorly understood. To explore transcriptional changes underlying FTLD-TDP, we performed RNA-sequencing on 66 genetically unexplained FTLD-TDP patients, 24 FTLD-TDP patients with GRN mutations and 24 control participants. Using principal component analysis, hierarchical clustering, differential expression and coexpression network analyses, we showed that GRN mutation carriers and FTLD-TDP-A patients without a known mutation shared a common transcriptional signature that is independent of GRN loss-of-function. After combining both groups, differential expression as compared to the control group and coexpression analyses revealed alteration of processes related to immune response, synaptic transmission, RNA metabolism, angiogenesis and vesicle-mediated transport. Deconvolution of the data highlighted strong cellular alterations that were similar in FTLD-TDP-A and GRN mutation carriers with NSF as a potentially important player in both groups. We propose several potentially druggable pathways such as the GABAergic, GDNF and sphingolipid pathways. Our findings underline new disease mechanisms and strongly suggest that affected pathways in GRN mutation carriers extend beyond GRN and contribute to genetically unexplained forms of FTLD-TDP-A.

Polygenic prediction of bipolar disorder in a Latin American sample.

To date, bipolar disorder (BD) genetic studies and polygenic risk scores (PRSs) for BD are based primarily on populations of European descent (EUR) and lack representation from other ancestries including Latin American (LAT). Here, we describe a new LAT cohort from the Mayo Clinic Bipolar Biobank (MCBB), a multisite collaboration with recruitment sites in the United States (EUR; 1,443 cases and 777 controls) and Mexico and Chile (LAT; 211 cases and 161 controls) and use the sample to explore the performance of a BD-PRS in a LAT population. Using results from the largest genome-wide association study of BD in EUR individuals, PRSice2 and LDpred2 were used to compute BD-PRSs in the LAT and EUR samples from the MCBB. PRSs explained up to 1.4% (PRSice) and 4% (LDpred2) of the phenotypic variance on the liability scale in the LAT sample compared to 3.8% (PRSice2) and 3.4% (LDpred2) in the EUR samples. Future larger studies should further explore the differential performance of different PRS approaches across ancestries. International multisite studies, such as this one, have the potential to address diversity-related limitations of prior genomic studies and ultimately contribute to the reduction of health disparities.

Testing and estimation of X-chromosome SNP effects: Impact of model assumptions.

Interest in analyzing X chromosome single nucleotide polymorphisms (SNPs) is growing and several approaches have been proposed. Prior studies have compared power of different approaches, but bias and interpretation of coefficients have received less attention. We performed simulations to demonstrate the impact of X chromosome model assumptions on effect estimates. We investigated the coefficient biases of SNP and sex effects with commonly used models for X chromosome SNPs, including models with and without assumptions of X chromosome inactivation (XCI), and with and without SNP-sex interaction terms. Sex and SNP coefficient biases were observed when assumptions made about XCI and sex differences in SNP effect in the analysis model were inconsistent with the data-generating model. However, including a SNP-sex interaction term often eliminated these biases. To illustrate these findings, estimates under different genetic model assumptions are compared and interpreted in a real data example. Models to analyze X chromosome SNPs make assumptions beyond those made in autosomal variant analysis. Assumptions made about X chromosome SNP effects should be stated clearly when reporting and interpreting X chromosome associations. Fitting models with SNP × Sex interaction terms can avoid reliance on assumptions, eliminating coefficient bias even in the absence of sex differences in SNP effect.

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.

Clinical Phenotype of Tardive Dyskinesia in Bipolar Disorder.

PURPOSE: Recognizing the negative impact that antipsychotic-induced movement disorders have on the quality of life and treatment outcomes in bipolar disorder (BD), this study aimed to assess clinical correlates and antipsychotic use patterns of tardive dyskinesia (TD+) in BD. MATERIALS AND METHODS: Participants with and without TD were included. Clinical variables were compared using t-test and χ2 test. Antipsychotic use patterns in TD+, including number of trials, mean doses, and estimated cumulative exposure, were assessed in a case-only analysis. RESULTS: The prevalence rate of TD was 5.1%. In comparison to the TD- group (n = 1074), TD+ participants (n = 58) were older, more likely to be female and have type I bipolar illness. There were 60.3% of the TD+ group that continued using antipsychotics at study entry and had a mean cumulative exposure to antipsychotics of 18.2 ± 15.6 years. Average dose, in haloperidol equivalents, was 5.9 ± 3.5 mg and 77.7% of the trials were second-generation antipsychotics. CONCLUSIONS: This study confirms previously identified TD risk factors, such as age, sex, and bipolar subtype in a large BD cohort. Limitations included a cross-sectional design and the lack of tardive illness severity assessment. As atypical antipsychotics continue to be primary mood stabilization treatment, attempting to harmonize large data sets to identify additional biomarkers of tardive risk will optimize individualized care for patients with BD.

TSPAN5 influences serotonin and kynurenine: pharmacogenomic mechanisms related to alcohol use disorder and acamprosate treatment response.

We previously reported that SNPs near TSPAN5 were associated with plasma serotonin (5-HT) concentrations which were themselves associated with selective serotonin reuptake inhibitor treatment outcomes in patients with major depressive disorder (MDD). TSPAN5 SNPs were also associated with alcohol consumption and alcohol use disorder (AUD) risk. The present study was designed to explore the biological function of TSPAN5 with a focus on 5-HT and kynurenine concentrations in the tryptophan pathway. Ethanol treatment resulted in decreased 5-HT concentrations in human induced pluripotent stem cell (iPSC)-derived neuron culture media, and the downregulation of gene expression of TSPAN5, DDC, MAOA, MAOB, TPH1, and TPH2 in those cells. Strikingly, similar observations were made when the cells were treated with acamprosate-an FDA approved drug for AUD therapy. These results were replicated in iPSC-derived astrocytes. Furthermore, TSPAN5 interacted physically with proteins related to clathrin and other vesicle-related proteins, raising the possibility that TSPAN5 might play a role in vesicular function in addition to regulating expression of genes associated with 5-HT biosynthesis and metabolism. Downregulation of TSPAN5 expression by ethanol or acamprosate treatment was also associated with decreased concentrations of kynurenine, a major metabolite of tryptophan that plays a role in neuroinflammation. Knockdown of TSPAN5 also influenced the expression of genes associated with interferon signaling pathways. Finally, we determined that TSPAN5 SNPs were associated with acamprosate treatment outcomes in AUD patients. In conclusion, TSPAN5 can modulate the concentrations of 5-HT and kynurenine. Our data also highlight a potentially novel pharmacogenomic mechanism related to response to acamprosate.

TCF7L2 lncRNA: a link between bipolar disorder and body mass index through glucocorticoid signaling.

Bipolar disorder (BD) and obesity are highly comorbid. We previously performed a genome-wide association study (GWAS) for BD risk accounting for the effect of body mass index (BMI), which identified a genome-wide significant single-nucleotide polymorphism (SNP) in the gene encoding the transcription factor 7 like 2 (TCF7L2). However, the molecular function of TCF7L2 in the central nervous system (CNS) and its possible role in the BD and BMI interaction remained unclear. In the present study, we demonstrated by studying human induced pluripotent stem cell (hiPSC)-derived astrocytes, cells that highly express TCF7L2 in the CNS, that the BD-BMI GWAS risk SNP is associated with glucocorticoid-dependent repression of the expression of a previously uncharacterized TCF7L2 transcript variant. That transcript is a long non-coding RNA (lncRNA-TCF7L2) that is highly expressed in the CNS but not in peripheral tissues such as the liver and pancreas that are involved in metabolism. In astrocytes, knockdown of the lncRNA-TCF7L2 resulted in decreased expression of the parent gene, TCF7L2, as well as alterations in the expression of a series of genes involved in insulin signaling and diabetes. We also studied the function of TCF7L2 in hiPSC-derived astrocytes by integrating RNA sequencing data after TCF7L2 knockdown with TCF7L2 chromatin-immunoprecipitation sequencing (ChIP-seq) data. Those studies showed that TCF7L2 directly regulated a series of BD risk genes. In summary, these results support the existence of a CNS-based mechanism underlying BD-BMI genetic risk, a mechanism based on a glucocorticoid-dependent expression quantitative trait locus that regulates the expression of a novel TCF7L2 non-coding transcript.

Genotype-dependent epigenetic regulation of DLGAP2 in alcohol use and dependence.

Alcohol misuse is a major public health problem originating from genetic and environmental risk factors. Alterations in the brain epigenome may orchestrate changes in gene expression that lead to alcohol misuse and dependence. Through epigenome-wide association analysis of DNA methylation from human brain tissues, we identified a differentially methylated region, DMR-DLGAP2, associated with alcohol dependence. Methylation within DMR-DLGAP2 was found to be genotype-dependent, allele-specific and associated with reward processing in brain. Methylation at the DMR-DLGAP2 regulated expression of DLGAP2 in vitro, and Dlgap2-deficient mice showed reduced alcohol consumption compared with wild-type controls. These results suggest that DLGAP2 may be an interface for genetic and epigenetic factors controlling alcohol use and dependence.

ERICH3: vesicular association and antidepressant treatment response.

Selective serotonin reuptake inhibitors (SSRIs) are standard of care for major depressive disorder (MDD) pharmacotherapy, but only approximately half of these patients remit on SSRI therapy. Our previous genome-wide association study identified a single-nucleotide polymorphism (SNP) signal across the glutamate-rich 3 (ERICH3) gene that was nearly genome-wide significantly associated with plasma serotonin (5-HT) concentrations, which were themselves associated with SSRI response for MDD patients enrolled in the Mayo Clinic PGRN-AMPS SSRI trial. In this study, we performed a meta-analysis which demonstrated that those SNPs were significantly associated with SSRI treatment outcomes in four independent MDD trials. However, the function of ERICH3 and molecular mechanism(s) by which it might be associated with plasma 5-HT concentrations and SSRI clinical response remained unclear. Therefore, we characterized the human ERICH3 gene functionally and identified ERICH3 mRNA transcripts and protein isoforms that are highly expressed in central nervous system cells. Coimmunoprecipitation identified a series of ERICH3 interacting proteins including clathrin heavy chain which are known to play a role in vesicular function. Immunofluorescence showed ERICH3 colocalization with 5-HT in vesicle-like structures, and ERICH3 knock-out dramatically decreased 5-HT staining in SK-N-SH cells as well as 5-HT concentrations in the culture media and cell lysates without changing the expression of 5-HT synthesizing or metabolizing enzymes. Finally, immunofluorescence also showed ERICH3 colocalization with dopamine in human iPSC-derived neurons. These results suggest that ERICH3 may play a significant role in vesicular function in serotonergic and other neuronal cell types, which might help explain its association with antidepressant treatment response.