Lithium-induced gene expression alterations in two peripheral cell models of bipolar disorder

Lithium response in bipolar disorder: Epigenome-wide DNA methylation signatures and epigenetic aging

Lithium (Li) is the first-line treatment for bipolar disorder (BD) even though only 30 % of BD patients are considered excellent responders. The mechanisms by which Li exerts its action are not clearly understood, but it has been suggested that specific epigenetic mechanisms, such as methylation processes, may play a role. In this regard, DNA methylation patterns can be used to estimate epigenetic age (EpiAge), which is accelerated in BD patients and reversed by Li treatment. Our first aim was to compare the DNA methylation profile in peripheral blood between BD patients categorized as excellent responders to Li (Ex-Rp) and non-responders (N-Rp). Secondly, EpiAge was estimated to detect differential age acceleration between the two groups. A total of 130 differentially methylated positions (DMPs) and 16 differentially methylated regions (DMRs) between Ex-Rp (n = 26) and N-Rp (n = 37) were identified (FDR adjusted p-value < 0.05). We found 122 genes mapping the DMPs and DMRs, nine of which (HOXB6, HOXB3, HOXB-AS3, TENM2, CACNA1B, ANK3, EEF2K, CYP1A1, and SORCS2) had previously been linked to Li response. We found genes related to the GSK3β pathway to be highly represented. Using FUMA, we found enrichment in Gene Ontology Cell Component for the synapse. Gene network analysis highlighted functions related to the cell cycle, nervous system development and function, and gene expression. No significant differences in age acceleration were found between Ex-Rp and N-Rp for any of the epigenetic clocks analysed. Our findings indicate that a specific methylation pattern could determine the response to Li in BD patients. We also found that a significant portion of the differentially methylated genes are closely associated with the GSK3β pathway, reinforcing the role of this system in Li response. Future longitudinal studies with larger samples will help to elucidate the epigenetic mechanisms underlying Li response.

RNA Biomarkers in Bipolar Disorder and Response to Mood Stabilizers

Bipolar disorder (BD) is a severe chronic disorder that represents one of the main causes of disability among young people. To date, no reliable biomarkers are available to inform the diagnosis of BD or clinical response to pharmacological treatment. Studies focused on coding and noncoding transcripts may provide information complementary to genome-wide association studies, allowing to correlate the dynamic evolution of different types of RNAs based on specific cell types and developmental stage with disease development or clinical course. In this narrative review, we summarize findings from human studies that evaluated the potential utility of messenger RNAs and noncoding transcripts, such as microRNAs, circular RNAs and long noncoding RNAs, as peripheral markers of BD and/or response to lithium and other mood stabilizers. The majority of available studies investigated specific targets or pathways, with large heterogeneity in the included type of cells or biofluids. However, a growing number of studies are using hypothesis-free designs, with some studies also integrating data on coding and noncoding RNAs measured in the same participants. Finally, studies conducted in neurons derived from induced-pluripotent stem cells or in brain organoids provide promising preliminary findings supporting the power and utility of these cellular models to investigate the molecular determinants of BD and clinical response.

Gene expression of circadian genes and CIART in bipolar disorder: A preliminary case-control study

Based on the observed circadian rhythms disruptions and sleep abnormalities in bipolar disorders (BD), a chronobiological model has been proposed suggesting that core clock genes play a central role in the vulnerability to the disorder. In this context, the analysis of circadian genes expression levels is particularly relevant, however studies focused on the whole set of core clock genes are scarce. We compared the levels of expression of 19 circadian genes (including the recently described circadian repressor (CIART)) in 37 euthymic individuals with BD and 20 healthy controls (HC), using data obtained by RNA sequencing of lymphoblastoid cell lines and validated the results using RT-qPCR. RNA sequencing data showed that CIART gene expression was correlated with those of ARNTL, ARNTL2, DBP, PER2 and TIMELESS. Data from RNA sequencing showed that the level of expression of four circadian genes (ARNTL, ARNTL2, BHLHE41 and CIART) discriminated individuals with BD from HC. We replicated this result using RT-qPCR for ARNTL and CIART. This study suggests that an imbalance between activation/repression of the transcription within the circadian system in individuals with BD as compared to HC and as such opens avenues for further research in larger independent samples combining both expression and epigenetic analyses.

Divergent Directionality of Immune Cell-Specific Protein Expression between Bipolar Lithium Responders and Non-Responders Revealed by Enhanced Flow Cytometry

Background and Objectives: There is no biomarker to predict lithium response. This study used CellPrint™ enhanced flow cytometry to study 28 proteins representing a spectrum of cellular pathways in monocytes and CD4⁺ lymphocytes before and after lithium treatment in patients with bipolar disorder (BD). Materials and Methods: Symptomatic patients with BD type I or II received lithium (serum level ≥ 0.6 mEq/L) for 16 weeks. Patients were assessed with standard rating scales and divided into two groups, responders (≥50% improvement from baseline) and non-responders. Twenty-eight intracellular proteins in CD4⁺ lymphocytes and monocytes were analyzed with CellPrint™, an enhanced flow cytometry procedure. Data were analyzed for differences in protein expression levels. Results: The intent-to-treat sample included 13 lithium-responders (12 blood samples before treatment and 9 after treatment) and 11 lithium-non-responders (11 blood samples before treatment and 4 after treatment). No significant differences in expression between the groups was observed prior to lithium treatment. After treatment, the majority of analytes increased expression in responders and decreased expression in non-responders. Significant increases were seen for PDEB4 and NR3C1 in responders. A significant decrease was seen for NR3C1 in non-responders. Conclusions: Lithium induced divergent directionality of protein expression depending on the whether the patient was a responder or non-responder, elucidating molecular characteristics of lithium responsiveness. A subsequent study with a larger sample size is warranted.

Transcriptomics and miRNomics data integration in lymphoblastoid cells highlights the key role of immune-related functions in lithium treatment response in Bipolar disorder

BACKGROUND

Bipolar Disorder (BD) is a complex mental disease characterized by recurrent episodes of mania and depression. Lithium (Li) represents the mainstay of BD pharmacotherapy, despite the narrow therapeutic index and the high variability in treatment response. However, although several studies have been conducted, the molecular mechanisms underlying Li therapeutic effects remain unclear.

METHODS

In order to identify molecular signatures and biological pathways associated with Li treatment response, we conducted transcriptome and miRNome microarray analyses on lymphoblastoid cell lines (LCLs) from 20 patients diagnosed with BD classified as Li responders (n = 11) or non-responders (n = 9).

RESULTS

We found 335 mRNAs and 77 microRNAs (miRNAs) significantly modulated in BD responders versus non-responders. Interestingly, pathway and network analyses on these differentially expressed molecules suggested a modulatory effect of Li on several immune-related functions. Indeed, among the functional molecular nodes, we found NF-κB and TNF. Moreover, networks related to these molecules resulted overall inhibited in BD responder patients, suggesting anti-inflammatory properties of Li. From the integrative analysis between transcriptomics and miRNomics data carried out using miRComb R package on the same samples from patients diagnosed with BD, we found 97 significantly and negatively correlated mRNA-miRNA pairs, mainly involved in inflammatory/immune response.

CONCLUSIONS

Our results highlight that Li exerts modulatory effects on immune-related functions and that epigenetic mechanisms, especially miRNAs, can influence the modulation of different genes and pathways involved in Li response. Moreover, our data suggest the potentiality to integrate data coming from different high-throughput approaches as a tool to prioritize genes and pathways.

ADHD-associated PARK2 copy number variants: A pilot study on gene expression and effects of supplementary deprivation in patient-derived cell lines

Recent studies show an association of Parkin RBR E3 ubiquitin protein ligase (PARK2) copy number variations (CNVs) with attention deficit hyperactivity disorder (ADHD). The aim of our pilot study to investigate gene expression associated with PARK2 CNVs in human-derived cellular models. We investigated gene expression in fibroblasts, hiPSC and dopaminergic neurons (DNs) of ADHD PARK2 deletion and duplication carriers by qRT PCR compared with healthy and ADHD cell lines without PARK2 CNVs. The selected 10 genes of interest were associated with oxidative stress response (TP53, NQO1, and NFE2L2), ubiquitin pathway (UBE3A, UBB, UBC, and ATXN3) and with a function in mitochondrial quality control (PINK1, MFN2, and ATG5). Additionally, an exploratory RNA bulk sequencing analysis in DNs was conducted. Nutrient deprivation as a supplementary deprivation stress paradigm was used to enhance potential genotype effects. At baseline, in fibroblasts, hiPSC, and DNs, there was no significant difference in gene expression after correction for multiple testing. After nutrient deprivation in fibroblasts NAD(P)H-quinone-dehydrogenase 1 (NQO1) expression was significantly increased in PARK2 CNV carriers. In a multivariate analysis, ubiquitin C (UBC) was significantly upregulated in fibroblasts of PARK2 CNV carriers. RNA sequencing analysis of DNs showed the strongest significant differential regulation in Neurontin (NNAT) at baseline and after nutrient deprivation. Our preliminary results suggest differential gene expression in pathways associated with oxidative stress, ubiquitine-proteasome, immunity, inflammation, cell growth, and differentiation, excitation/inhibition modulation, and energy metabolism in PARK2 CNV carriers compared to wildtype healthy controls and ADHD patients.

Transcriptional biomarkers of response to pharmacological treatments in severe mental disorders: A systematic review

Variation in the expression level and activity of genes involved in drug disposition and action in tissues of pharmacological importance have been increasingly investigated in patients treated with psychotropic drugs. Findings are promising, but reliable predictive biomarkers of response have yet to be identified. Here we conducted a PRISMA-compliant systematic search of PubMed, Scopus and PsycInfo up to 12 September 2020 for studies investigating RNA expression levels in cells or biofluids from patients with major depressive disorder, schizophrenia or bipolar disorder characterized for response to psychotropic drugs (antidepressants, antipsychotics or mood stabilizers) or adverse effects. Among 5497 retrieved studies, 123 (63 on antidepressants, 33 on antipsychotics and 27 on mood stabilizers) met inclusion criteria. Studies were either focused on mRNAs (n = 96), microRNAs (n = 19) or long non-coding RNAs (n = 1), with only a minority investigating both mRNAs and microRNAs levels (n = 7). The most replicated results include genes playing a role in inflammation (antidepressants), neurotransmission (antidepressants and antipsychotics) or mitochondrial function (mood stabilizers). Compared to those investigating response to antidepressants, studies focused on antipsychotics or mood stabilizers more often showed lower sample size and lacked replication. Strengths and limitations of available studies are presented and discussed in light of the specific designs, methodology and clinical characterization of included patients for transcriptomic compared to DNA-based studies. Finally, future directions of transcriptomics of psychopharmacological interventions in psychiatric disorders are discussed.

Altered estradiol-dependent cellular Ca2+ homeostasis and endoplasmic reticulum stress response in Premenstrual Dysphoric Disorder

Premenstrual Dysphoric Disorder (PMDD) is characterized by debilitating mood symptoms in the luteal phase of the menstrual cycle. Prior studies of affected women have implicated a differential response to ovarian steroids. However, the molecular basis of these patients' differential response to hormone remains poorly understood. We performed transcriptomic analyses of lymphoblastoid cell lines (LCLs) derived from women with PMDD and asymptomatic controls cultured under untreated (steroid-free), estradiol-treated (E2), and progesterone-treated (P4) conditions. Weighted gene correlation network analysis (WGCNA) of transcriptomes identified four gene modules with significant diagnosis x hormone interactions, including one enriched for neuronal functions. Next, in a gene-level analysis comparing transcriptional response to hormone across diagnoses, a generalized linear model identified 1522 genes differentially responsive to E2 (E2-DRGs). Among the top 10 E2-DRGs was a physically interacting network (NUCB1, DST, GCC2, GOLGB1) involved in endoplasmic reticulum (ER)-Golgi function. qRT-PCR validation reproduced a diagnosis x E2 interaction (F(1,24)=7.01, p = 0.014) for NUCB1, a regulator of cellular Ca²⁺ and ER stress. Finally, we used a thapsigargin (Tg) challenge assay to test whether E2 induces differences in Ca²⁺ homeostasis and ER stress response in PMDD. PMDD LCLs had a 1.36-fold decrease in Tg-induced XBP1 splicing response compared to controls, and a 1.62-fold decreased response (p = 0.005), with a diagnosis x treatment interaction (F(3,33)=3.51, p = 0.026) in the E2-exposed condition. Altered hormone-dependent in cellular Ca²⁺ dynamics and ER stress may contribute to the pathophysiology of PMDD.

Whole blood transcriptome analysis in bipolar disorder reveals strong lithium effect

BACKGROUND

Bipolar disorder (BD) is a highly heritable mood disorder with complex genetic architecture and poorly understood etiology. Previous transcriptomic BD studies have had inconsistent findings due to issues such as small sample sizes and difficulty in adequately accounting for confounders like medication use.

METHODS

We performed a differential expression analysis in a well-characterized BD case-control sample (Nsubjects = 480) by RNA sequencing of whole blood. We further performed co-expression network analysis, functional enrichment, and cell type decomposition, and integrated differentially expressed genes with genetic risk.

RESULTS

While we observed widespread differential gene expression patterns between affected and unaffected individuals, these effects were largely linked to lithium treatment at the time of blood draw (FDR < 0.05, Ngenes = 976) rather than BD diagnosis itself (FDR < 0.05, Ngenes = 6). These lithium-associated genes were enriched for cell signaling and immune response functional annotations, among others, and were associated with neutrophil cell-type proportions, which were elevated in lithium users. Neither genes with altered expression in cases nor in lithium users were enriched for BD, schizophrenia, and depression genetic risk based on information from genome-wide association studies, nor was gene expression associated with polygenic risk scores for BD.

CONCLUSIONS

These findings suggest that BD is associated with minimal changes in whole blood gene expression independent of medication use but emphasize the importance of accounting for medication use and cell type heterogeneity in psychiatric transcriptomic studies. The results of this study add to mounting evidence of lithium's cell signaling and immune-related mechanisms.

Exploring lithium's transcriptional mechanisms of action in bipolar disorder: a multi-step study

Lithium has been the first-line treatment for bipolar disorder (BD) for more than six decades. Although the molecular effects of lithium have been studied extensively and gene expression changes are generally believed to be involved, the specific mechanisms of action that mediate mood regulation are still not known. In this study, a multi-step approach was used to explore the transcriptional changes that may underlie lithium's therapeutic efficacy. First, we identified genes that are associated both with lithium exposure and with BD, and second, we performed differential expression analysis of these genes in brain tissue samples from BD patients (n = 42) and healthy controls (n = 42). To identify genes that are regulated by lithium exposure, we used high-sensitivity RNA-sequencing of corpus callosum (CC) tissue samples from lithium-treated (n = 8) and non-treated (n = 9) rats. We found that lithium exposure significantly affected 1108 genes (FDR < 0.05), 702 up-regulated and 406 down-regulated. These genes were mostly enriched for molecular functions related to signal transduction, including well-established lithium-related pathways such as mTOR and Wnt signaling. To identify genes with differential expression in BD, we performed expression quantitative trait loci (eQTL) analysis on BD-associated genetic variants from the most recent genome-wide association study (GWAS) using three different gene expression databases. We found 307 unique eQTL genes regulated by BD-associated variants, of which 12 were also significantly modulated by lithium treatment in rats. Two of these showed differential expression in the CC of BD cases: RPS23 was significantly down-regulated (p = 0.0036, fc = 0.80), while GRIN2A showed suggestive evidence of down-regulation in BD (p = 0.056, fc = 0.65). Crucially, GRIN2A was also significantly up-regulated by lithium in the rat brains (p = 2.2e-5, fc = 1.6), which suggests that modulation of GRIN2A expression may be a part of the therapeutic effect of the drug. These results indicate that the recent upsurge in research on this central component of the glutamatergic system, as a target of novel therapeutic agents for affective disorders, is warranted and should be intensified.

The identification of biomarkers predicting acute and maintenance lithium treatment response in bipolar disorder: A plea for further research attention

The prediction of acute and maintenance lithium treatment response carries major clinical and neurobiological implications, warranting systematic review. A Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) compliant review searched major electronic databases from inception until December 2017 for studies documenting a clinical diagnosis of bipolar disorder (BD) made according to the mainstream diagnostic manuals and confirmed by a structured interview. Eligible studies allowed a quantitative comparison of endpoint vs baseline mean values of a given biomarker, regardless of the mood phase of patients with BD, and the disorder was assessed for severity using validated rating tool(s). Owing to the purposely applied stringent selection criteria, 16 acute and 12 maintenance studies could be included. The anticipated publication bias limited the chances of reportable generalizable findings, hindering a side-by-side comparison of different records across varying biomarkers and subsequent meta-analyses. The PRISMA approach was nonetheless preferred; it aimed at enhancing the homogeneity of the included results and minimizing the chances of "apples and oranges" with respect to the present research theme. The present critical review confirms the need for future research to specifically assess either pretreatment and/or posttreatment putative biomarkers of patients with BD and treated with lithium.