Search

We characterized the heterogeneity and risk factors of cognitive decline in euthymic bipolar disorder (BD), and their magnitude of associations with subjective daily functions. In this retrospective cohort, BD type I patients ( N = 128) were followed for an average of 6.5 years. Intelligence quotient (IQ) at index date was recorded, and premorbid IQ was estimated. We used Brief Assessment of Cognition in Affective Disorders (BAC-A) to assess cognition at follow-up. We evaluated current functions with World Health Organization Disability Assessment Schedule 2.0. Clinical and sociodemographic factors were examined for their independent effects on longitudinal cognitive decline. In addition, we employed multivariate adaptive regression spline to detect inflection points for the nature of slope changes in cognitive decline among BD patients. During follow-up years, 21 BD patients (16.4%) showed longitudinal cognitive decline. In cognitive decline group, all cognitive domains of BAC-A were significantly worsened. We found that density of episodes with psychotic features was an independent risk factor for cognitive decline after adjusted for age, gender and dose of mood stabilizer. After the age of 42 years, a steeper cognitive change was observed in the cognitive decline group. The correlation pattern between cognitive domains and functional outcomes differed between patients with and without cognitive decline. The present study characterized cognitive heterogeneity longitudinally in BD patients. As density of episodes play roles for cognitive decline, our results emphasize the importance of relapse prevention. Our findings provide hints for future personalized interventions and facilitating genetic and biological studies for dissecting the heterogeneity of bipolar illness., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Chen, Huang, Lee, Chang, Lin, Chiu, Liu, Kuo, Weng, Chen and Kuo.)

Bipolar disorder (BD) is highly heritable and well known for its recurrent manic and depressive episodes. The present study focused on manic episode in BD patients and aimed to investigate state-specific transcriptome alterations between acute episode and remission, including messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and micro-RNAs (miRNAs), using microarray and RNA sequencing (RNA-Seq) platforms. BD patients were enrolled with clinical information, and peripheral blood samples collected at both acute and remission status spanning for at least 2 months were confirmed by follow-ups. Symptom severity was assessed by Young Mania Rating Scale. We enrolled six BD patients as the discovery samples and used the Affymetrix Human Transcriptome Array 2.0 to capture transcriptome data at the two time points. For replication, expression data from Gene Expression Omnibus that consisted of 11 BD patients were downloaded, and we performed a mega-analysis for microarray data of 17 patients. Moreover, we conducted RNA sequencing (RNA-Seq) in additional samples of 7 BD patients. To identify intraindividual differentially expressed genes (DEGs), we analyzed data using a linear model controlling for symptom severity. We found that noncoding genes were of majority among the top DEGs in microarray data. The expression fold change of coding genes among DEGs showed moderate to high correlations (∼0.5) across platforms. A number of lncRNAs and two miRNAs ( MIR181B1 and MIR103A1 ) exhibited high levels of gene expression in the manic state. For coding genes, we reported that the taste function-related genes, including TAS2R5 and TAS2R3 , may be mania state-specific markers. Additionally, four genes showed a nominal p -value of less than 0.05 in all our microarray data, mega-analysis, and RNA-Seq analysis. They were upregulated in the manic state and consisted of MS4A14 , PYHIN1 , UTRN , and DMXL2 , and their gene expression patterns were further validated by quantitative real-time polymerase chain reaction (PCR) (qRT-PCR). We also performed weight gene coexpression network analysis to identify gene modules for manic episode. Genes in the mania-related modules were different from the susceptible loci of BD obtained from genome-wide association studies, and biological pathways in relation to these modules were mainly related to immune function, especially cytokine-cytokine receptor interaction. Results of the present study elucidated potential molecular targets and genomic networks that are involved in manic episode. Future studies are needed to further validate these biomarkers for their roles in the etiology of bipolar illness.

Background: High sedative use in a major depressive episode may imply specific clinical features. This study aims to examine the correlation between sedative use and clinical severity indicators in the initial treatment phase of first-onset major depressive disorder. Methods: A study cohort in the first episode of major depressive disorder was used to conduct pharmacological dissection. All participants had at least a 2-year follow-up period with a complete treatment record. The defined daily dose of antidepressants and augmentation agents were calculated as the antidepressant load and augmentation load, respectively. Sedative use, which was calculated as the equivalent dosage of lorazepam, were defined as the sedative load. These psychotropic loads were measured monthly and the averaged psychotropic loads for each day were obtained. Results: A total of 106 individuals (75.5% female) were included. The mean duration of disease course in participants was 5.5 ± 3.5 years. In the multiple regression analysis, after controlling for other classes of psychotropics and comorbid anxiety disorders, the sedative load independently correlated with higher number of antidepressants used, higher number of antidepressant used with an adequate dose and duration, more psychiatric emergency and outpatient visits within 2 years of disease onset. Conclusion: High loading of sedatives correlated with several indicators of clinical severity in major depressive disorder. The sedative load may be used as a specifier to identify subgroups in patients with major depressive disorder.

In the wake of recent advances in scientific research, personalized medicine using deep learning techniques represents a new paradigm. In this work, our goal was to establish deep learning models which distinguish responders from non-responders, and also to predict possible antidepressant treatment outcomes in major depressive disorder (MDD). To uncover relationships between the responsiveness of antidepressant treatment and biomarkers, we developed a deep learning prediction approach resulting from the analysis of genetic and clinical factors such as single nucleotide polymorphisms (SNPs), age, sex, baseline Hamilton Rating Scale for Depression score, depressive episodes, marital status, and suicide attempt status of MDD patients. The cohort consisted of 455 patients who were treated with selective serotonin reuptake inhibitors (treatment-response rate = 61.0%; remission rate = 33.0%). By using the SNP dataset that was original to a genome-wide association study, we selected 10 SNPs (including ABCA13 rs4917029, BNIP3 rs9419139, CACNA1E rs704329, EXOC4 rs6978272, GRIN2B rs7954376, LHFPL3 rs4352778, NELL1 rs2139423, NUAK1 rs2956406, PREX1 rs4810894, and SLIT3 rs139863958) which were associated with antidepressant treatment response. Furthermore, we pinpointed 10 SNPs (including ARNTL rs11022778, CAMK1D rs2724812, GABRB3 rs12904459, GRM8 rs35864549, NAALADL2 rs9878985, NCALD rs483986, PLA2G4A rs12046378, PROK2 rs73103153, RBFOX1 rs17134927, and ZNF536 rs77554113) in relation to remission. Then, we employed multilayer feedforward neural networks (MFNNs) containing 1-3 hidden layers and compared MFNN models with logistic regression models. Our analysis results revealed that the MFNN model with 2 hidden layers (area under the receiver operating characteristic curve (AUC) = 0.8228 ± 0.0571; sensitivity = 0.7546 ± 0.0619; specificity = 0.6922 ± 0.0765) performed maximally among predictive models to infer the complex relationship between antidepressant treatment response and biomarkers. In addition, the MFNN model with 3 hidden layers (AUC = 0.8060 ± 0.0722; sensitivity = 0.7732 ± 0.0583; specificity = 0.6623 ± 0.0853) achieved best among predictive models to predict remission. Our study indicates that the deep MFNN framework may provide a suitable method to establish a tool for distinguishing treatment responders from non-responders prior to antidepressant therapy.

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