Next Article in Journal
Whole-Transcriptome Analysis Reveals Potential CeRNA Regulatory Mechanism in Takifugu rubripes against Cryptocaryon irritans Infection
Next Article in Special Issue
Integrating Molecular Perspectives: Strategies for Comprehensive Multi-Omics Integrative Data Analysis and Machine Learning Applications in Transcriptomics, Proteomics, and Metabolomics
Previous Article in Journal
Quantifying the Ecological Performance of Migratory Bird Conservation: Evidence from Poyang Lake Wetlands in China
Previous Article in Special Issue
Contemporary Update on Clinical and Experimental Prostate Cancer Biomarkers: A Multi-Omics-Focused Approach to Detection and Risk Stratification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Codes between Poles: Linking Transcriptomic Insights into the Neurobiology of Bipolar Disorder

by
Jon Patrick T. Garcia
1,2 and
Lemmuel L. Tayo
1,3,*
1
School of Chemical, Biological, and Materials Engineering and Sciences, Mapúa University, Manila 1002, Philippines
2
School of Graduate Studies, Mapúa University, Manila 1002, Philippines
3
Department of Biology, School of Health Sciences, Mapúa University, Makati 1200, Philippines
*
Author to whom correspondence should be addressed.
Biology 2024, 13(10), 787; https://doi.org/10.3390/biology13100787
Submission received: 17 August 2024 / Revised: 2 September 2024 / Accepted: 30 September 2024 / Published: 30 September 2024
(This article belongs to the Special Issue Multi-omics Data Integration in Complex Diseases)

Abstract

:

Simple Summary

Bipolar disorder is a psychiatric condition in which one prominent symptom is the regular occurrence of mood swings. Its root cause remains unclear; thus, this study was intended to provide new insights to better understand the genetics underlying such a disorder. Using samples obtained from three regions in the brain, it was found that some molecular and cellular mechanisms are involved in the disruption of neurobiological processes. The dysregulated expression of certain genes eventually leads to neurogenesis and neurotransmission impairment events in humans. Furthermore, these genes significant in the onset of bipolar disorder were identified and evaluated for the presence of variants, which may be targeted to engineer better curative treatment strategies for the disorder.

Abstract

Bipolar disorder (BPD) is a serious psychiatric condition that is characterized by the frequent shifting of mood patterns, ranging from manic to depressive episodes. Although there are already treatment strategies that aim at regulating the manifestations of this disorder, its etiology remains unclear and continues to be a question of interest within the scientific community. The development of RNA sequencing techniques has provided newer and better approaches to studying disorders at the transcriptomic level. Hence, using RNA-seq data, we employed intramodular connectivity analysis and network pharmacology assessment of disease-associated variants to elucidate the biological pathways underlying the complex nature of BPD. This study was intended to characterize the expression profiles obtained from three regions in the brain, which are the nucleus accumbens (nAcc), the anterior cingulate cortex (AnCg), and the dorsolateral prefrontal cortex (DLPFC), provide insights into the specific roles of these regions in the onset of the disorder, and present potential targets for drug design and development. The nAcc was found to be highly associated with genes responsible for the deregulated transcription of neurotransmitters, while the DLPFC was greatly correlated with genes involved in the impairment of components crucial in neurotransmission. The AnCg did show association with some of the expressions, but the relationship was not as strong as the other two regions. Furthermore, disease-associated variants or single nucleotide polymorphisms (SNPs) were identified among the significant genes in BPD, which suggests the genetic interrelatedness of such a disorder and other mental illnesses. DRD2, GFRA2, and DCBLD1 were the genes with disease-associated variants expressed in the nAcc; ST8SIA2 and ADAMTS16 were the genes with disease-associated variants expressed in the AnCg; and FOXO3, ITGA9, CUBN, PLCB4, and RORB were the genes with disease-associated variants expressed in the DLPFC. Aside from unraveling the molecular and cellular mechanisms behind the expression of BPD, this investigation was envisioned to propose a new research pipeline in studying the transcriptome of psychiatric disorders to support and improve existing studies.

1. Introduction

Manic depression, commonly known as bipolar disorder (BPD), is a serious mental health condition characterized by a recurrent shifting of mood between two emotional polarities, which are mania or hypomania (highs) and depression (lows) [1]. Globally, around 2.4% of the population is clinically diagnosed with BPD [2]. No cure has been discovered for the disorder yet, although current medications have shown promising results for regulating its symptoms. Several studies have already been performed that elucidate the genetic underpinnings of BPD, but there are only limited studies discussing the manifestations of the disorder at the transcriptomic level. Krebs et al. performed differential expression analysis (DEA) of the transcriptomic profiles obtained from whole blood samples in which the results suggested the possible role of lithium in the cellular mechanisms pertinent among BPD patients [3]. Moreover, it was observed from a study by Holmgren et al. that dysregulation of stoichiometric balance has implications for the pathology of BPD [4]. On the other hand, Ardesch et al. analyzed the normative expression profiles of risk genes linked to four neuropsychiatric disorders, including BPD, and concluded brain dysfunction as one of the shared pathways among the illnesses [5]. Most transcriptomic studies of this disorder are commonly investigated as a group, along with other brain conditions [6,7,8]. However, when studied in isolation, some of these only present the genes and pathways underlying the disorder but fail to interpret the systematic relationship of gene sets to elucidate possible molecular and cellular associations in bringing new insights and routes for inventing novel medications or improving existing treatments against the deregulated expression of BPD.
The past decades have shown a growing interest in the use of RNA sequencing as a technique in the study of transcriptomics. As a newer and more advanced technology compared to conventional sequencing methods such as microarrays, RNA-seq data utilizes next-generation sequencing (NGS), which provides more accurate and reliable results that could lead to deeper insight into the transcriptome [9,10]. Most transcriptomic studies of BPD used blood samples for sequencing; however, this raises concerns as whole blood may confer limitations in explaining the deviations in brain function among people with the disorder [11,12,13]. However, it is undeniable that blood is commonly utilized for sampling because it is more accessible and less invasive than other biological sources. To provide a more holistic understanding of the manifestations of BPD, post-mortem examination of brain samples may be favored. A study by Luykx et al. performed RNA sequencing and DEA of medial frontal gyrus brain tissues to determine the implication of protein-coding and lncRNAs among patients with BPD [14]. Furthermore, using hippocampal and orbitofrontal cortex samples, Darby et al. highlighted the potential of mRNAs and protein synthesis as therapeutic targets for treating BPD, schizophrenia (SCZ), and major depressive disorder (MDD) [15]. Ellis et al., on the other hand, utilized cortical brain tissues of patients with autism, SCZ, and BPD to quantify the extent of genetic correlation among these three [16]. Similar to what was mentioned above, most of these studies analyzed and compared BPD relative to other psychiatric disorders. Therefore, results from the current literature may not support medical scenarios in which BPD is regarded as an independent case.
Unfortunately, the use of post-mortem brain samples possesses considerable limitations since the rate of degradation of the tissue can incur a significant impact on the analysis. This method is greatly time-dependent; therefore, the post-mortem interval (PMI), which is the interval between the time of death and the time the sample will be analyzed, must be kept short. In the study of Liharska et al., they inferred that post-mortem brain samples should not be used as a proxy for living brain samples because both show vast differences in expression profiles [17]. However, such a conclusion was considered premature and unsupported by relevant data. True enough that RNA transcripts readily start degrading shortly after death, but this is not an absolute reason to completely consider post-mortem brain samples as insignificant for human brain research. In fact, several studies suggest otherwise. Fromer et al. found that expression profiles sequenced from post-mortem brain samples have a genetic correlation with the occurrence of brain disorders [18]. Furthermore, transcriptomic analyses performed in the studies of Raj et al. and Gandal et al. showed that predicted associations of genetic underpinnings among certain diseases were proved factual based on the differential expression of genes in post-mortem-controlled cases [19,20].
With the general objective of providing a deeper understanding of the molecular and cellular causes of BPD, this study aimed to perform intramodular connectivity analysis and network pharmacology assessment of diseases-associated variants using different samples from three regions in the brain: the nucleus accumbens (nAcc), the anterior cingulate cortex (AnCg), and the dorsolateral prefrontal cortex (DLPFC). We intended to use post-mortem samples that have been previously RNA-sequenced to contribute more relevant information into the transcriptome of BPD and further support existing studies with a similar goal but utilizing whole blood transcripts. The intramodular connectivity analysis was employed to elucidate what gene expression patterns are highly correlated to which region of the brain that would give insights into their specific roles and present potential therapeutic targets for treating the disorder. Aside from that, network pharmacology was carried out to determine which genes from these regions are greatly associated with BPD. These were then further investigated by obtaining their first-degree neighbors and identifying the presence of disease-associated variants or single nucleotide polymorphisms (SNPs) that would explain the genetic interrelatedness of BPD and other psychiatric disorders.

2. Materials and Methods

2.1. Nature of the Dataset

The RNA-seq data utilized in this study were obtained from the National Center for Biotechnology Information–Gene Expression Omnibus (NCBI–GEO) [https://www.ncbi.nlm.nih.gov/geo/ (accessed on 15 June 2024)]. GEO from NCBI is an online repository of genomics data from various research and medical institutes where sequencing data of a vast range of organisms are stored for public use. GSE80655 was the chosen dataset derived from a previous study by Ramaker et al. [21]. Only data pertaining to BPD were included, and the rest were excluded. Twenty-four post-mortem brain tissues of human (Homo sapiens) patients diagnosed with BPD, with each sample obtained from the nAcc, AnCg, and DLPFC regions, were sequenced using Illumina HiSeq 2000, equating to a total of 72 samples used for this study. A statistical description of the samples to explain in detail the clinical data controls and demographic profiles of their sources is presented in Table S1.

2.2. Pre-Processing of Dataset

The raw dataset was loaded into the R program v4.2.3 [https://cran.r-project.org/bin/windows/base/ (downloaded on 21 June 2024)], where the initial computation of the dataset was performed. Using Bioconductor v3.19 [https://www.bioconductor.org/ (accessed on 21 June 2024)], which is a free, open-source software for facilitating comprehensive analyses of biological assays, we filtered and modified the dataset for quality control before it was subjected to intramodular analysis. The “goodSamplesGenes” function was employed to check and remove genes with zero variance and missing values or entries. Moreover, hierarchal clustering using the “hclust” function was performed to detect any outlying samples. Principal component analysis was then followed to provide a more quantitative basis for excluding the outliers. After which, DESeq2 was utilized to examine the expression patterns among the genes from the RNA-seq data, and all those genes with counts of less than 15 in more than 75% of the samples were removed.

2.3. Intramodular Connectivity Analysis

Variance stabilization and normalization were performed as the final steps before proceeding with WGCNA. Using the “pickSoftThreshold” function, the network topology was constructed to determine the soft-thresholding power (β) to which the co-expression similarity values of the genes will be raised. The mean connectivity and correlation coefficient values were extracted and plotted to visualize the power of choice based on a scale-free topology model fit limit of 80% to ensure that the adjacency matrices that were calculated in estimating the expression modules were highly accurate. Furthermore, the “blockwiseModules” function was employed using a block size of 20,000 to obtain the modules that were constructed from the correlation of genes with respect to the calculated module eigengenes. The “plotDendroAndColors” function was then used to plot the distribution of the modules in the dendrogram and describe how these were identified from the hierarchal clustering of the genes. After which, the respective module membership and p-value of each gene were calculated, and the categorical variables, the nAcc, AnCg, and DLPFC, were binarized to compute the module–gene associations that would elucidate what genetic expressions have high intramodular connectivity. Finally, using “CorLevelPlot”, a heatmap was generated to determine only the significant modules that would accurately characterize any of the three brain regions.

2.4. Functional Annotation and Pathway Enrichment of Modules

From the heatmap, those modules labeled in red with three asterisks were recorded and subjected to functional annotation and pathway enrichment. The higher the p-value, the darker the color of the module, which suggests a greater positive correlation among the genes grouped in that module. Meanwhile, the higher the number of asterisks, the more likely that that region will characterize the expression patterns in that module. Each module was then imported to Cytoscape v.3.10.2 [https://cytoscape.org/index.html (downloaded on 15 July 2024)] to construct their PPI network under a confidence level of 0.90 with the aid of stringApp v2.1.1 [https://apps.cytoscape.org/apps/stringapp (downloaded on 10 July 2024)]. From the generated large network of interactions, the top 50 hub genes were determined using the degree topological algorithm from cytoHubba v0.1 [https://apps.cytoscape.org/apps/cytohubba (downloaded on 15 July 2024)]. Subsequently, the hub genes of each module were sent to ShinyGO v0.80 [http://bioinformatics.sdstate.edu/go/ (accessed on 17 July 2024)] for functional annotation and pathway enrichment [22]. These gene sets were characterized based on three gene ontology classifications: biological process (BP), cellular component (CC), and molecular function (MF). Aside from that, KEGG pathway analysis was employed to further identify the expression patterns in the modules [23,24]. An FDR cutoff of 0.05 was used, and any redundant category among the lists was removed. Only those clusters with the highest fold enrichment score from each annotation were sent to the SRplot web server [http://www.bioinformatics.com.cn/srplot (accessed on 16 July 2024)] for enrichment visualization.

2.5. Network Pharmacology Assessment of Disease-Associated Variants

The gene significance of the expression profiles obtained from the three brain regions was calculated and ranked amongst each other. The top 100 genes from the nAcc, AnCg, and DLPFC were recorded and sent to DisGeNET [https://disgenet.com (accessed on 17 July 2024)] for disease-association analysis [25]. Consequently, those genes that are found positive for disease-associated variants or SNPs were further evaluated by determining their first-degree neighbors using a whole-genome PPI network constructed from the expression profiles of BPD patients under an interaction score of 0.7.

3. Results

3.1. Intramodular Connectivity Analysis

After pre-processing the raw RNA-seq data, the network topology and mean connectivity were calculated to determine the soft-thresholding power that was used to identify the modules. Figure S1 shows the scale-free topology model fit and mean connectivity plots of the dataset from which the soft-thresholding power was chosen. Since β = 12 obtained the highest mean connectivity among the values that reached the 80% threshold, it was utilized to calculate the gene co-expression networks. From there, the module eigengenes were analyzed, and each gene was quantitatively classified based on their connectivity relative to the module eigengene to construct the modules that significantly characterize the expression profiles prominent among the three brain regions. Figure 1 displays the dendrogram, which describes how the modules were identified based on the height of the branches that correlate to the strength of connectivity among sets of genes.
A total of 25 modules, including the grey module, were obtained from the WGCNA. Furthermore, the selection was filtered down by calculating the respective module membership and p-value of the genes within each module. Seven out of the twenty-five modules showed high intramodular connectivity. The pink, dark green, grey60, red, and purple modules were highly correlated to the nAcc, while the turquoise and blue modules were highly correlated to the DLPFC. The AnCg did not show a high correlation to any of the modules, as observed by the absence of modules with three asterisks. Figure 2 presents the module–gene connectivity heatmap of the samples wherein positive correlations are labeled in pink while negative correlations are labeled in blue. The darker the color of the module, the greater the connectivity among its genes. Moreover, the number of asterisks indicates the likelihood that the module is highly expressed in that region. A greater number of asterisks in that module suggests greater intramodular connectivity, or in other words, a higher probability of that certain expression occurring in that region.

3.2. Functional Annotation and Pathway Enrichment

The seven modules obtained from the previous analysis were sent to Cytoscape to obtain their top 50 hub genes using the degree algorithm for functional annotation and pathway enrichment. The hub genes were then sent to ShinyGO to characterize the modules based on BP, CC, and MF annotations. These were also subjected to KEGG pathway analysis to determine the possible routes to which these sets of genes are similar. Figure 3 presents the annotation bar plots of the cluster with the highest fold enrichment from each module. The top five highest clusters are also listed in Table S2 with their respective enrichment FDR score, ID, pathway description, and genes.
In the BP annotation, the blue module corresponded to synaptic vesicle exocytosis; the dark green module to epithelial cilium movement involved in extracellular fluid movement; the grey60 module to myeloid leukocyte activation; the pink module to epinephrine metabolic process; the purple module to Golgi lumen acidification; the red module to carboxylic acid catabolic process; and the turquoise module to mitochondrial electron transport ubiquinol to cytochrome c. In the CC annotation, the blue module corresponded to synaptobrevin 2-SNAP-25-syntaxin-1a-complexin I complex; the dark green module to outer dynein arm; the grey60 to tertiary granule membrane; the pink module to beta-catenin-TCF complex; the purple module to proteosome regulatory particle lid subcomplex; the red module to heterotrimetric G-protein complex; and the turquoise module to mitochondrial respiratory chain complex III. In the MF annotation, the blue module corresponded to phosphatidylinositol-3,4-bisphosphate 5-kinase activity; the dark green module to minus-end-directed microtubule motor activity; the grey60 module to ICAM-3 receptor activity; the pink module to opioid peptide activity; the red module to glutamate dehydrogenase (NAD+) activity; and the turquoise module to succinate dehydrogenase (ubiquinone) activity.
On the other hand, Figure 4 shows the bubble plot of the cluster with the highest enrichment FDR from each module obtained from the KEGG pathway analysis.
The top three clusters associated with the blue module were insulin secretion, synaptic vesicle cycle, and inflammatory mediator regulation of TRP channels; the dark green module were Huntington's disease, amyotrophic lateral sclerosis, and pathways of neurodegeneration-multiple diseases; the grey60 module were staphylococcus aureus infection, complement and coagulation cascades, and leishmaniasis; the pink module were tyrosine metabolism, cocaine addiction, and serotonergic synapse; the purple module were vibrio cholerae infection, homologous recombination, and proteasome; the red module were histidine metabolism, beta-alanine metabolism, and nitrogen metabolism; and the turquoise module were longevity regulating pathway-multiple species, EGFR tyrosine kinase inhibitor resistance, and ErbB signaling pathway.

3.3. Network Pharmacology Assessment of Disease-Associated Variants

The genes with the highest significance value among the nAcc, AnCg, and DLPFC were sent to DisGeNET for disease-association analysis. In the nAcc, the genes found to be associated with BPD were DRD2, PPP1R1B, CHRNA2, GFRA2, DCBLD1, PENK, NOS1AP, PART1, and NRN1; in the AnCg, the genes found to be associated with BPD were GRIK4, ST8SIA2, FABP7, TSHZ1, ADRA1A, ADAMTS16, HTR7, PCDH17, and PNPLA3; and in the DLPFC, the genes found to be associated with BPD were IGFBP2, SLC17A6, FOXO3, ITGA9, CARTPT, PREP, NEFM, CUBN, PLCB4, NOS1, NOS2, FGF9, RORB, SGCG, and ZNF365. Moreover, these genes were further evaluated for the presence of SNPs. Among the reported genes in the nAcc, DRD2, GFRA2, and DCBLD1 presented potential SNP interference; among the reported genes in the AnCg, ST8SIA2 and ADAMTS16 presented potential SNP interference; and among the reported genes in the DLPFC, FOXO3, ITGA9, CUBN, PLCB4, and RORB presented potential SNP interference. Table 1 shows the summary of disease-associated variants found among the three brain regions. Meanwhile, the complete list of significant genes, along with their associated SNPs, is listed in Table S3.
Furthermore, those genes with disease-associated variants were sent again to Cytoscape to determine their first-degree neighbors using the PPI network constructed from the expression profiles of BPD patients. Figure 5 displays the network interactions of the first-degree neighbors of DRD2, which was the gene with the highest GDA score among the genes with SNPs found in the three brain regions. The PPI network of the other genes is shown in Figure S2.

4. Discussion

4.1. Genetic Evidence of Bipolar Disorder Occurrence

Contrary to popular belief, genetics does play a major role in the increased risk of BPD in children who happen to have parents with the disorder. Although environmental stimulation indeed aggravates the manifestation of BPD, disparities in the transcriptome among positive patients relative to healthy individuals present strong evidence of its genetic cause. Recent discoveries from genomic studies of BPD have further proven its intricate and complex nature. A genome-wide association study (GWAS) by Stahl et al. found 30 loci that were genome-wide significant, of which 20 were newly reported [26]. Moreover, a much earlier GWAS conducted by the Psychiatric Genomics Consortium, involving 41,917 BPD patients and 371,549 negative controls, revealed 64 genomic locations that are associated with an increased risk of BPD, 33 of which were never found in previous literature [27]. These results greatly improved the initial findings of the former study, bringing better insights into determining more risk locations in the DNA. Furthermore, a study by the Stanley Center for Psychiatric Research compared the exomes of 14,000 BPD patients against 14,000 negative controls with the hope of identifying the sequences that cause dysregulation in the synthesis of certain proteins in BPD [28]. The analysis showed that the AKAP11 gene is a strong genetic factor that increases the risk of having the disorder.
These are only some of the several studies that prove the genetic underpinnings of BPD. Newer technologies can bring more novel findings that could further the scientific understanding of such disorders, especially since researchers now are not limited to using genomic approaches. Analyzing the transcriptome could also shed light on discovering molecular and cellular mechanisms that are evident in its occurrence. Moreover, delving into certain samples, such as the brain, where psychiatric disorders are most likely linked, may aid in identifying specific regions that would be best targeted for treating BPD.

4.2. Deregulated Transcription of Neurotransmitters in the Nucleus Accumbens

Adjacent to the septal nuclei, the nAcc is a region that is part of the ventral striatum, located in the rostrobasal forebrain [29,30]. It acts as the neural interface between motivation and action and is considered a significant part of the reward system of the brain [31]. Dysfunction of the nAcc has shown associations with various psychiatric disorders (i.e., obsessive–compulsive disorder, addiction, and MDD) [32,33,34]. The disruption of reward processing in mood disorders, such as that in MDD, can lead to anhedonia, which is the inability to perceive pleasure, and it is caused by an abnormal function of the nAcc [35]. Atypical stimulation of stressors in this region may also prime anxiety disorders (ADs), ultimately influencing emotional and behavioral responses. Moreover, the fluctuating activity of the nAcc in the brain can affect the instability of mood, a common manifestation in BPD, in which higher activation may exacerbate manic periods while lower activation may lead to depressive episodes [36]. Thus, in a variety of psychiatric diseases, dysregulation of the nAcc can take a toll on the intricate molecular and cellular interactions that drive motivation, emotional regulation, and reward processing.
It has long been known how important neurotransmitters are in neural homeostasis. Some of these neurotransmitters that play a significant function in the manifestation of mood swings in BPD are dopamine, serotonin, and norepinephrine. The manic polarity of BPD is associated with elevated levels of dopamine and increased norepinephrine activation, which result in impulsivity and the feeling of euphoria [37]. Meanwhile, the depressive polarity of BPD is linked to low levels of serotonin and increased secretion and conversion of norepinephrine, which result in irritability and the feeling of melancholy [38]. The complex interplay of these neurotransmitters, along with glutamate, an excitatory neurotransmitter, and gamma-aminobutyric acid (GABA), an inhibitory neurotransmitter, further underscores the multifaceted nature of the disorder [39,40]. Unlike other neurotransmitters, epinephrine is rarely mentioned in the neurobiology of BPD. Dopamine and norepinephrine, which are the precursors of epinephrine, are more involved in the central nervous system (CNS), which explains their direct effect on the brain. The blood–brain barrier (BBB) is impermeable against epinephrine, making it a more essential neurotransmitter in the peripheral nervous system (PNS) [41]. Dopamine and norepinephrine are significantly synthesized in the brain through the amino acid tyrosine, and their relative concentrations greatly affect the shifting of moods in BPD. Although epinephrine could be synthesized from dopamine and norepinephrine inside the brain, it is only in small amounts and its effects are less pronounced than its precursors. The genes found to be associated with the epinephrine metabolic process in the BP annotation simply suggest that the dysregulation of the nAcc could be influencing the elevated levels of dopamine and norepinephrine in the brain. This notion implies that the nAcc might be one of those overstimulated regions that contribute to the increased concentrations of these neurotransmitters, causing manic episodes in BPD. Otherwise, epinephrine metabolism would not occur if there was an absence of dopamine and norepinephrine.
The β-catenin/TCF complex is a transcriptional regulator that activates the signaling of Wnt target gene transcription in the nucleus [42]. This complex is essential in synthesizing the proteins necessary for the normal functioning of the body, thus, the disruption in its regulation may lead to serious disorders. Impairment in the Wnt pathway has been associated with various psychiatric disorders in humans, such as WNT1, WNT2, and WNT7A in autism spectrum disorder [43], WNT7B and LRP5 in SCZ [44], and LRP5 and LRP6 in attention-deficit/hyperactivity disorder [45]. In BPD, WNT2B and WNT7A are reported to be involved [46]. The β-catenin/TCF complex has major implications on neurotransmission and neurogenesis in the brain. On the one hand, the interaction of β-catenin with TCF transcription factors influences the expression of genes involved in neurotransmission and synaptic function [47]. Deviations occurring in these interactions may result in deregulated secretion and activation of neurotransmitters that eventually manifest in BPD. On the other hand, signaling pathways induced by the β-catenin/TCF complex are important in promoting the differentiation of neural progenitor cells during neurogenesis [48]. Impairment of these signaling pathways could lead to structural brain changes; this is characterized by smaller hippocampal volumes reported in patients with BPD [49]. The genes found to be associated with the beta-catenin-TCF complex in the CC annotation simply suggest that problems in transcription regulation could be occurring in the nAcc. This further explains that the dysregulation of Wnt pathways in this brain region might be a cause of the unrestricted production and release of neurotransmitters in BPD, that is, signifying its relationship with the epinephrine metabolic process in the BP annotation.
Opioid peptides, including endorphins, enkephalins, and dynorphins, are also neurotransmitters that interact with various opioid receptors in the brain and nervous system to modulate a wide range of neural processes. In psychiatric disorders, opioid peptides are implicated in the deregulation of emotional and behavioral responses [50]. Just like the role of dopamine, serotonin, and norepinephrine in the occurrence of mood swings in BPD, concentration levels of opioids in the brain also characterize the complex shifting patterns between mania and depression [51]. Moreover, opioid peptide–receptor interactions regulate the activation of the hypothalamus–pituitary–adrenal (HPA) gland, and this influences the production and secretion of crucial neurotransmitters in BPD [52]. Although there are still limited studies that back up the role of the opioid system in the symptoms of BPD, a few pieces of literature have paved solid and promising evidence of its genetic foundation. In the study of Qi et al., they concluded that OPRM1, a Mu-opioid receptor (MOR) gene, is one of the candidate biomarkers for the increased risk of BPD [53]. Another study by Lindtröm et al. and Genazzani et al. reported that a fluctuation of endorphin levels in the cerebrospinal fluid (CSF) is manifested among manic patients [54,55]. The genes found to be associated with opioid peptide activity in the MF annotation simply suggest the important role of these neurotransmitters in the molecular etiology of BPD. Along with the result from the BP annotation, this greatly enlightens how the dysregulation of the transcription process in the nAcc affects the mechanistic interplay of nervous and endocrine ligand–receptor complexes in the expression of BPD.
In the KEGG pathway analysis, the pink module presented genes from the nAcc responsible for tyrosine metabolism. This implies that there may be disruption in the metabolic pathway of tyrosine in the body, which could further explain the deregulated conditions of neurotransmitter secretion manifested in BPD. Tyrosine is an essential element in the production of dopamine and norepinephrine. Interference in the metabolic activity of tyrosine can negatively impact the synthesis of dopamine, which could eventually affect norepinephrine concentrations in the brain [56,57].

4.3. Neurotransmission Impairment in the Dorsolateral Prefrontal Cortex

The DLPFC is a region located in the frontal lobe, in the prefrontal cortex specifically, that is responsible for the higher cognitive functions of the brain. It is adjacent to the limbic system, which is involved in the regulation of behavior and emotion [58]. The DLPFC often exhibits decreased activity during manic periods, which may impede the capacity of the brain to manage impulses. It is also associated with a reduction in inhibitory control that can exacerbate risk-taking behaviors and impair judgment [59]. Meanwhile, during depressive conditions, the DLPFC becomes inefficient in interacting with other brain regions responsible for mood regulation. This leads to hypoactivity, which, in turn, results in difficulties in concentration and decision-making [60]. Aside from modulating emotional responses, the DLPFC has become an interest in the study of BPD because of its implications in the cognitive aspect of the brain, and the dysregulation of expression patterns in this region may have significant effects on certain manifestations of the disorder.
Neurotransmitters, the chemical messengers involved in neurotransmission, are crucial for regulating various brain functions and processes. Imbalance or dysfunction in neurotransmitter systems can lead to the symptoms expressed among psychiatric disorders. Synaptic vesicle exocytosis is a fundamental element in the transmission process of nerve impulses to send the necessary signals to and from neurons to respond to external stimuli [61]. Signal neurotransmission occurs when neurotransmitters enclosed in a vesicle diffuse across the synaptic cleft and bind to receptors in the postsynaptic membrane to be delivered from one neuron to another. This interaction results in the opening of ion channels which generate electrical signals to carry out cellular responses. Dysregulation of such process could lead to neural miscommunication that may result in unrestricted activation, secretion, production, or differentiation of various neurotransmitters, in which increased or decreased concentrations of these chemicals could alter the normal functioning of the brain.
Moreover, impairment in synaptic exocytosis could also be a result of dysfunctional transmembrane proteins responsible for the specific attachment of exosomes. The SNARE complex is an important cellular component in exocytotic activities during neurotransmission, and alterations in its function may contribute to the neurochemical imbalances that underlie the characteristic mood swings and emotional instability in BPD. Three key proteins make up the SNARE complex: synaptobrevin, syntaxin, and SNAP-25. Synaptobrevin is found on the synaptic vesicle, while syntaxin and SNAP-25 are found on the presynaptic membrane [62]. In BPD, unregulated neurotransmission may cause an abnormal relay of signals between neurons, and this might be potentially correlated to an overactive or dysfunctional SNARE complex [63]. As a consequence, neurotransmitters become unregulated in the brain and manifest into various mood symptoms.
The genes found to be associated with synaptic vesicle exocytosis in the BP annotation, synaptobrevin 2-SNAP-25-syntaxin-1a-complexin I complex in the CC annotation, and phosphatidylinositol-3,4-bisphosphate 5-kinase activity in the MF annotation simply suggest that in BPD, the DLPFC could be characterized by the expression of genes responsible for the disruption of neurotransmission processes prominent in the brain. This notion may be a reason why, during either manic or depressive episodes, the DLPFC exhibits molecular and cellular disjunction with other brain regions that lead to poor cognitive functions. This could also be a cause of dysregulated phosphorylation of proteins, resulting in the eventual disruption of the said process. Such involvement of the DLPFC proves its role in aggravating the complex mood patterns observed in BPD.
Furthermore, the KEGG pathway analysis showed that the DLPFC may express genes associated with insulin secretion. The excitatory and inhibitory impulses in the brain are stabilized by insulin receptors on presynaptic neurons, which regulate neurotransmitter release. Insulin facilitates the translocation of glucose transporters to the cell surface, which is essential for glucose uptake in neurons and glial cells [64]. This ensures that neurons have enough energy for generating action potentials, recycling synaptic vesicles, and synthesizing neurotransmitters, especially glutamate and GABA, which greatly rely on glucose metabolism [65]. These genes determined from the DPLFC suggest the presence of irregularities in the secretion pathway of insulin in the body, and this problem may interfere with the generation of action potentials between neurons, which could be one component in the manifestation of neurotransmission impairment in BPD.

4.4. Risk Genes and Polymorphisms in Bipolar Disorder

Scientists have now been targeting gene candidates and biomarkers that present the risk of BPD at the earlier stages. From the analysis, among the genes with the highest significance value, we found that the nAcc has the greatest association with DRD2, the AnCg with GRIK4, and the DLPFC with IGFBP2.
DRD2 (dopamine receptor D2) is an important component in the dopamine system of the body [66]. Alterations in this gene have associated it with various psychiatric disorders such as SCZ and drug addiction disorders. In a study by Huang et al., they reported DRD2 as one of the candidate genes in the etiology of BPD [67]. DRD2 interacts with ANKK1 to initialize a certain dopaminergic pathway that causes the differentiation of BPD-1 and BPD-2. Similar results were observed in the study of Lee et al., wherein they concluded that disparities in the interaction between DRD2/ANKK1 and ALDH2 genes characterize the subtypes of BPD [68].
Moreover, GRIK4 (glutamate receptor, ionotropic kainite 4) is a type of ion channel that facilitates excitatory neurotransmission processes in the brain by binding to glutamate, a neurotransmitter [69]. One of the prominent effects of a dysfunctional GRIK4 is the disruption of glutamatergic signaling, which is essential for efficient neurotransmission and synaptic plasticity. Knight et al. reported that the reduced expression of KA1, as a result of the absence of a deletion variant in GRIK4, was associated with an increased risk of BPD [70]. KA1 is essential in the increase of GRIK4 mRNA abundance which mitigates symptoms of BPD. A 3′ UTR deletion variant that exhibited protection against this disorder was also identified by Pickard et al. [71]. The results of their study suggested that an increase in kainate receptor activity caused by drugs, which can mimic the expression-modifying effect of such a deletion variant, may offer a treatment target for BPD and other psychiatric disorders.
Meanwhile, IGFBP2 (insulin-like growth factor binding protein 2) is found to be involved in certain cellular mechanisms in the brain that are vital in neurogenesis and neuroprotection. Alterations in this gene have shown associations with neuroinflammation and neurodegeneration, resulting in a loss of emotional and cognitive function among multifarious psychiatric disorders [72]. In the study conducted by Bezchlibnyk et al., they observed reduced expressions of IGFBP2 in patients with BPD relative to controls, suggesting the role of this gene in the pathophysiology of the disorder [73]. They also reported a greater reduction of expression profiles in patients who were not treated with lithium, which further supports their claim regarding the function of IGFBP2 in the etiology of BPD. Additionally, the same conclusions were drawn from the study of Fernandez and Torres-Alemán, in which they stated that the unregulated expression of IGFBP2 has implications for the proper development of the brain, thereby affecting the survival, proliferation, and differentiation of neurons [74].
The presence of SNPs can impact how the manifestations of BPD are phenotypically expressed, such that these can affect the onset and intensity of its symptoms. These genetic variations may be linked to an increased risk of rapid emotional cycling or mood episodes, more severe manic or depressive periods, or earlier inception of the condition. Due to the vast range of clinical symptoms caused by SNPs, BPD diagnosis and treatment can become more challenging.
One study presented the association of BPD with multiple SNPs of the IL1B and IL6R genes [75]. The researchers proved that the impact of SNPs in BPD has a potential influence on the relationship between the immune and nervous systems, which can characterize molecular and cellular mechanisms involved in psychiatric disorders. The presence of SNPs can also explain why BPD, in any way, is correlated to other mood disorders such as SCZ and MDD. These genetic variations can modify biological pathways in the brain that may result in similar mechanisms found in other disorders. As previously mentioned, in the study of Mikhalitskaya et al., they found rs2228145, SNP of the ILFR6 gene, to be linked with BPD, SCZ, and MDD [75]. de Marco et al., on the other hand, recorded 66 SNPs among 29 genes that have associations with both BPD and substance use disorder [76]. Some of these were rs11600996 of the ARNTL gene, rs6971 of the TSPO gene, and rs945032 of the BDKRB2 gene. Moreover, from the genome assessment analysis conducted by Wang et al., they reported the presence of rs11789399, SNP of the ASTN2 gene, in both BPD and SCZ cases [77].
SNPs contribute to the high rates of comorbidity and clinical overlap between BPD and other psychiatric conditions [78]. For instance, individuals with BPD may exhibit symptoms of anxiety aside from manic and depressive episodes. Such SNPs responsible for the disruption of genes that control mood, stress response, and emotional regulation could also increase the likelihood of developing AD or MDD. Thus, studying the transcriptome of BPD is one way of evaluating the genetic overlap among psychiatric disorders to assess biomarkers that could be targeted for therapy.
Most available drugs in the market are designed to target only specific proteins: for instance, remoxipride, which is an antipsychotic drug selective for dopamine D2 receptors [79], ketamine, which is a dissociative anesthetic that targets ionotropic glutamate receptors [80,81], and mecasermin, which is used as a treatment for growth failure in patients with IGF-1 deficiency [81]. These drugs are commonly utilized for maintaining proper regulation of disrupted proteins and not as curative strategies. Current medications for BPD are long-term, depending on the severity of the disorder and the patient’s response to the treatment. However, emerging drug therapy approaches for BPD are aiming to achieve longer periods of remission [82]. The risk genes and variants identified in this study are useful in leveraging pharmacogenomic research for drug design optimization that will target the cause of BPD at the molecular level, which is contrary to most drugs that only deal with the cellular ramifications of the disorder. Targeting the location of the genes dysregulated in BPD may provide better and longer curative effects because the root cause is directed and not the translated proteins themselves. The CRISPR technology is also one promising technique for gene-editing and -regulation, and its advancement serves as a breakthrough in the field of pharmacogenomics [83]. The risk genes and variants found in this study may serve as a preliminary step in determining the guide RNA sequences that could be used for editing the genome of patients with BPD. Additionally, gene replacement, gene silencing, or gene knock-out may be employed to edit these BPD-causing genes and variants [84]. However, research on gene therapy is still ongoing, but despite the current limitations of this method, its purpose still holds great promise in the future of pharmacogenomic techniques for engineering curative strategies for BPD.
Using a computational and pharmacological assessment of the transcriptome of BPD, we have found that the nAcc region is involved in the expression of genes responsible for the deregulated transcription of neurotransmitters, while the DLPFC region is involved in the expression of genes responsible for neurotransmission impairment. These results primarily underscore the neurobiological aspect of BPD, which could further enlighten the molecular and cellular factors that lead to the phenotypic manifestations of the disorder. Aside from that, among the risk genes identified, DRD2 was determined to be greatly associated with the nAcc, GRIK4 with the AnCg, and IGFBP2 with the DLPFC, and most of them showed the presence of SNPs.

5. Conclusions

The complex nature of BPD is one reason why the study of its etiology is challenging, just like other psychiatric conditions. However, advances in laboratory technology, such as NGS, paved the way for researchers to better understand the transcriptome of this disorder, enlightening the biological mechanisms that drive the severity of its manifestations. In this study, we performed intramodular connectivity analysis and network pharmacology assessment of the expression profiles of patients diagnosed with BPD to elucidate the molecular and cellular causes of the disorder. The results suggest that the nAcc was highly associated with the expression of genes responsible for the deregulated transcription of neurotransmitters. Meanwhile, the DLPFC was greatly correlated to genes involved in the impairment of components in neurotransmission. Relative to the two brain regions, the AnCg was not found to be highly correlated to any of the identified expressions, though it showed considerable association with other mechanisms evident in BPD. Disease-associated variants were also found present among the significant genes of the disorder and may be targeted for innovating better treatment strategies.
Future objectives may investigate other brain regions that would better elucidate the neurobiological manifestations of BPD. The SNPs reported in this study may also be further evaluated to expound on the genetic connectivity of BPD with other psychiatric disorders. This study only covered the genetic factors affecting BPD; thus, it would greatly benefit the current literature if the environmental and epigenetic conditions that influence the disorder were tackled as well. Finally, DEA may be employed to identify which genes may be upregulated or downregulated, which could provide more insights into future transcriptomic studies of BPD. We also emphasize that post-mortem brain samples were used in this study, and that may incur biological limitations posed by certain inevitable factors such as PMI, patient history, or ethnicity. The data were also initially analyzed through in vitro methods and then assessed using in silico computations, so the pharmacological findings in this study may require further in vivo experimentations or be subjected to advanced clinical research. In the future, as more transcriptomic data become available, the complex interactions among genes, RNAs, and proteins will become clearer and eventually lead to a better understanding of life.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology13100787/s1, Figure S1: Scale-free topology model fit and mean connectivity plots of the pre-processed dataset. Based on the figure, β = 12 was chosen as the soft-thresholding power in identifying the modules among the three brain regions because it attained the highest mean connectivity among the values that reached the 80% limit; Figure S2: Modified PPI network of the genes that were identified to have disease-associated variants among the three brain regions: (a) CUBN; (b) ST8SIA1; (c) GFRA2; (d) ITGA9; (e) RORB; (f) PLCB4; (g) FOXO3; Table S1: Summarized statistic table of the clinical data controls and demographic profiles of the samples. All samples were of Caucasian descent; Table S2: The top 5 clusters from each annotation of the modules identified to have strong expression of gene sets among the three brain regions; Table S3: List of genes found to have a significant association with BPD. Those which showed the presence of SNPs were also recorded. Moreover, degree corresponds to the number of first-degree neighbors adjacent to the gene, which was obtained from the PPI network constructed from the expression profiles of BPD patients.

Author Contributions

Conceptualization, J.P.T.G.; methodology, J.P.T.G.; software, J.P.T.G.; validation, J.P.T.G. and L.L.T.; formal analysis, J.P.T.G.; investigation, J.P.T.G.; resources, J.P.T.G.; data curation, J.P.T.G.; writing—original draft preparation, J.P.T.G.; writing—review and editing, J.P.T.G. and L.L.T.; visualization, J.P.T.G.; supervision, L.L.T.; project administration, L.L.T.; funding acquisition, L.L.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original data presented in the study are openly available in GEO-NCBI at 10.1186/s13073-017-0458-5 or GSE80655.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Siegel-Ramsay, J.E.; Sharp, S.J.; Ulack, C.J.; Chiang, K.S.; di Scalea, T.L.; O’Hara, S.; Carberry, K.; Strakowski, S.M.; Suarez, J.; Teisberg, E.; et al. Experiences that matter in bipolar disorder: A qualitative study using the capability, comfort and calm framework. Int. J. Bipolar Disord. 2023, 11, 13. [Google Scholar] [CrossRef] [PubMed]
  2. Zhong, Y.; Chen, Y.; Su, X.; Wang, M.; Li, Q.; Shao, Z.; Sun, L. Global, regional and national burdens of bipolar disorders in adolescents and young adults: A trend analysis from 1990 to 2019. Gen. Psychiatry 2024, 37, e101255. [Google Scholar] [CrossRef]
  3. Krebs, C.E.; Ori, A.P.S.; Vreeker, A.; Wu, T.; Cantor, R.M.; Boks, M.P.M.; Kahn, R.S.; Loohuis, L.M.O.; Ophoff, R.A. Whole blood transcriptome analysis in bipolar disorder reveals strong lithium effect. Psychol. Med. 2020, 50, 2575–2586. [Google Scholar] [CrossRef]
  4. Holmgren, A.; Akkouh, I.; O’Connell, K.S.; Osete, J.R.; Bjørnstad, P.M.; Djurovic, S.; Hughes, T. Bipolar patients display stoichiometric imbalance of gene expression in post-mortem brain samples. Mol. Psychiatry 2024, 29, 1128–1138. [Google Scholar] [CrossRef]
  5. Ardesch, D.J.; Libedinsky, I.; Scholtens, L.H.; Wei, Y.; van den Heuvel, M.P. Convergence of Brain Transcriptomic and Neuroimaging Patterns in Schizophrenia, Bipolar Disorder, Autism Spectrum Disorder, and Major Depressive Disorder. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2023, 8, 630–639. [Google Scholar] [CrossRef] [PubMed]
  6. Wongkoblap, A.; Vadillo, M.A.; Curcin, V. Researching Mental Health Disorders in the Era of Social Media: Systematic Review. J. Med. Internet Res. 2017, 19, E228. [Google Scholar] [CrossRef] [PubMed]
  7. Malla, A.; Joober, R.; Garcia, A. “Mental illness is like any other medical illness”: A critical examination of the statement and its impact on patient care and society. J. Psychiatry Neurosci. 2015, 40, 147. [Google Scholar] [CrossRef] [PubMed]
  8. Lyu, N.; Wang, H.; Zhao, Q.; Fu, B.; Li, J.; Yue, Z.; Huang, J.; Yang, F.; Liu, H.; Zhang, L.; et al. Peripheral biomarkers to differentiate bipolar depression from major depressive disorder: A real-world retrospective study. BMC Psychiatry 2024, 24, 1–11. [Google Scholar] [CrossRef] [PubMed]
  9. Stark, R.; Grzelak, M.; Hadfield, J. RNA sequencing: The teenage years. Nat. Rev. Genet. 2019, 20, 631–656. [Google Scholar] [CrossRef] [PubMed]
  10. Garcia, J.P.T.; Tayo, L.L. Theoretical Studies of DNA Microarray Present Potential Molecular and Cellular Interconnectivity of Signaling Pathways in Immune System Dysregulation. Genes 2024, 15, 393. [Google Scholar] [CrossRef] [PubMed]
  11. Sakrajda, K.; Bilska, K.; Czerski, P.M.; Narożna, B.; Dmitrzak-Węglarz, M.; Heilmann-Heimbach, S.; Brockschmidt, F.F.; Herms, S.; Nöthen, M.M.; Cichon, S.; et al. Abelson Helper Integration Site 1 haplotypes and peripheral blood expression associates with lithium response and immunomodulation in bipolar patients. Psychopharmacology 2024, 241, 727–738. [Google Scholar] [CrossRef] [PubMed]
  12. Petty, F.; Kramer, G.L.; Fulton, M.; Moeller, F.G.; Rush, A.J. Low Plasma GABA Is a Trait-Like Marker for Bipolar Illness. Neuropsychopharmacology 1993, 9, 125–132. [Google Scholar] [CrossRef] [PubMed]
  13. Bella, F.; Muscatello, M.R.A.; D’Ascola, A.; Campo, S. Gene Expression Analysis of nc-RNAs in Bipolar and Panic Disorders: A Pilot Study. Genes 2023, 14, 1778. [Google Scholar] [CrossRef] [PubMed]
  14. Luykx, J.; Giuliani, F.; Veldink, J.; Kahn, R. RNA Sequencing in Bipolar Disorder: From Long Non-coding to Circular rnas. Eur. Psychiatry 2017, 41, S56. [Google Scholar] [CrossRef]
  15. Darby, M.M.; Yolken, R.H.; Sabunciyan, S. Consistently altered expression of gene sets in postmortem brains of individuals with major psychiatric disorders. Transl. Psychiatry 2016, 6, e890. [Google Scholar] [CrossRef] [PubMed]
  16. Ellis, S.E.; Panitch, R.; West, A.B.; Arking, D.E. Transcriptome analysis of cortical tissue reveals shared sets of downregulated genes in autism and schizophrenia. Transl. Psychiatry 2016, 6, e817. [Google Scholar] [CrossRef] [PubMed]
  17. Liharska, L.E.; Park, Y.J.; Ziafat, K.; Wilkins, L.; Silk, H.; Linares, L.M.; Thompson, R.C.; Vornholt, E.; Sullivan, B.; Cohen, V.; et al. A study of gene expression in the living human brain. MedRxiv 2023, 2, 23288916. [Google Scholar] [CrossRef]
  18. Fromer, M.; Roussos, P.; Sieberts, S.K.; Johnson, J.S.; Kavanagh, D.H.; Perumal, T.M.; Ruderfer, D.M.; Oh, E.C.; Topol, A.; Shah, H.R.; et al. Gene Expression Elucidates Functional Impact of Polygenic Risk for Schizophrenia. Nat. Neurosci. 2016, 19, 1442. [Google Scholar] [CrossRef] [PubMed]
  19. Raj, T.; Li, Y.I.; Wong, G.; Humphrey, J.; Wang, M.; Ramdhani, S.; Wang, Y.C.; Ng, B.; Gupta, I.; Haroutunian, V.; et al. Integrative transcriptome analyses of the aging brain implicate altered splicing in Alzheimer’s disease susceptibility. Nat. Genet. 2018, 50, 1584–1592. [Google Scholar] [CrossRef]
  20. Gandal, M.J.; Zhang, P.; Hadjimichael, E.; Walker, R.L.; Chen, C.; Liu, S.; Won, H.; Van Bakel, H.; Varghese, M.; Wang, Y.; et al. Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science 2018, 362, eaat8127. [Google Scholar] [CrossRef]
  21. Ramaker, R.C.; Bowling, K.M.; Lasseigne, B.N.; Hagenauer, M.H.; Hardigan, A.A.; Davis, N.S.; Gertz, J.; Cartagena, P.M.; Walsh, D.M.; Vawter, M.P.; et al. Post-mortem molecular profiling of three psychiatric disorders. Genome Med. 2017, 9, 72. [Google Scholar] [CrossRef]
  22. Ge, S.X.; Jung, D.; Jung, D.; Yao, R. ShinyGO: A graphical gene-set enrichment tool for animals and plants. Bioinformatics 2020, 36, 2628–2629. [Google Scholar] [CrossRef] [PubMed]
  23. Kanehisa, M.; Furumichi, M.; Sato, Y.; Ishiguro-Watanabe, M.; Tanabe, M. KEGG: Integrating viruses and cellular organisms. Nucleic Acids Res. 2021, 49, D545–D551. [Google Scholar] [CrossRef] [PubMed]
  24. Luo, W.; Brouwer, C. Pathview: An R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics 2013, 29, 1830–1831. [Google Scholar] [CrossRef] [PubMed]
  25. Piñero, J.; Ramírez-Anguita, J.M.; Saüch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef]
  26. Stahl, E.A.; Breen, G.; Forstner, A.J.; McQuillin, A.; Ripke, S.; Trubetskoy, V.; Mattheisen, M.; Wang, Y.; Coleman, J.R.I.; Gaspar, H.A.; et al. Genome-wide association study identifies 30 Loci Associated with Bipolar Disorder. Nat. Genet. 2019, 51, 793. [Google Scholar] [CrossRef] [PubMed]
  27. Mullins, N.; Forstner, A.J.; O’Connell, K.S.; Coombes, B.; Coleman, J.R.I.; Qiao, Z.; Als, T.D.; Bigdeli, T.B.; Børte, S.; Bryois, J.; et al. Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology. Nat. Genet. 2021, 53, 817–829. [Google Scholar] [CrossRef] [PubMed]
  28. Palmer, D.S.; Howrigan, D.P.; Chapman, S.B.; Adolfsson, R.; Bass, N.; Blackwood, D.; Boks, M.P.M.; Chen, C.Y.; Churchhouse, C.; Corvin, A.P.; et al. Exome sequencing in bipolar disorder identifies AKAP11 as a risk gene shared with schizophrenia. Nat. Genet. 2022, 54, 541–547. [Google Scholar] [CrossRef] [PubMed]
  29. Scofield, M.D.; Heinsbroek, J.A.; Gipson, C.D.; Kupchik, Y.M.; Spencer, S.; Smith, A.C.W.; Roberts-Wolfe, D.; Kalivas, P.W. The Nucleus Accumbens: Mechanisms of Addiction across Drug Classes Reflect the Importance of Glutamate Homeostasis. Pharmacol. Rev. 2016, 68, 816. [Google Scholar] [CrossRef]
  30. Shivacharan, R.S.; Rolle, C.E.; Barbosa, D.A.N.; Cunningham, T.N.; Feng, A.; Johnson, N.D.; Safer, D.L.; Bohon, C.; Keller, C.; Buch, V.P.; et al. Pilot study of responsive nucleus accumbens deep brain stimulation for loss-of-control eating. Nat. Med. 2022, 28, 1791–1796. [Google Scholar] [CrossRef]
  31. Harris, H.; Peng, Y. Evidence and explanation for the involvement of the nucleus accumbens in pain processing. Neural Regen. Res. 2020, 15, 597. [Google Scholar] [CrossRef] [PubMed]
  32. Shi, Y.; Wang, M.; Xiao, L.; Gui, L.; Zheng, W.; Bai, L.; Su, B.; Li, B.; Xu, Y.; Pan, W.; et al. Potential therapeutic mechanism of deep brain stimulation of the nucleus accumbens in obsessive-compulsive disorder. Front. Cell. Neurosci. 2023, 16, 1057887. [Google Scholar] [CrossRef] [PubMed]
  33. Shirayama, Y.; Chaki, S. Neurochemistry of the Nucleus Accumbens and its Relevance to Depression and Antidepressant Action in Rodents. Curr. Neuropharmacol. 2006, 4, 277. [Google Scholar] [CrossRef] [PubMed]
  34. Koob, G.F.; Volkow, N.D. Neurocircuitry of Addiction. Neuropsychopharmacology 2010, 35, 217–238. [Google Scholar] [CrossRef]
  35. Liu, R.; Wang, Y.; Chen, X.; Zhang, Z.; Xiao, L.; Zhou, Y. Anhedonia correlates with functional connectivity of the nucleus accumbens subregions in patients with major depressive disorder. Neuroimage Clin. 2021, 30, 102599. [Google Scholar] [CrossRef]
  36. Misaki, M.; Suzuki, H.; Savitz, J.; Drevets, W.C.; Bodurka, J. Individual Variations in Nucleus Accumbens Responses Associated with Major Depressive Disorder Symptoms. Sci. Rep. 2016, 6, 21227. [Google Scholar] [CrossRef] [PubMed]
  37. Wolfenberger, T.; Diaz, A.P.; Bockmann, T.; Selvaraj, S.; Sanches, M.; Soares, J.C. Predominant polarity and associated post-traumatic stress disorder in patients with comorbid bipolar disorder and borderline personality disorder: A cross-sectional study. Braz. J. Psychiatry 2022, 44, 557. [Google Scholar] [CrossRef] [PubMed]
  38. Azorin, J.M.; Kaladjian, A.; Adida, M.; Fakra, E.; Belzeaux, R.; Hantouche, E.; Lancrenon, S. Factors associated with borderline personality disorder in major depressive patients and their relationship to bipolarity. Eur. Psychiatry 2013, 28, 463–468. [Google Scholar] [CrossRef]
  39. Shen, J.; Tomar, J.S. Elevated Brain Glutamate Levels in Bipolar Disorder and Pyruvate Carboxylase-Mediated Anaplerosis. Front. Psychiatry 2021, 12, 640977. [Google Scholar] [CrossRef]
  40. Brady, R.O.; Mccarthy, J.M.; Prescot, A.P.; Jensen, J.E.; Cooper, A.J.; Cohen, B.M.; Renshaw, P.F.; Ongür, D. Brain gamma-aminobutyric acid (GABA) abnormalities in bipolar disorder. Bipolar Disord. 2013, 15, 434. [Google Scholar] [CrossRef]
  41. Martinho, R.; Oliveira, A.; Correia, G.; Marques, M.; Seixas, R.; Serrão, P.; Moreira-Rodrigues, M. Epinephrine May Contribute to the Persistence of Traumatic Memories in a Post-traumatic Stress Disorder Animal Model. Front. Mol. Neurosci. 2020, 13, 588802. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, H.J.; Hsu, L.S.; Shia, Y.T.; Lin, M.W.; Lin, C.M. The β-catenin/TCF complex as a novel target of resveratrol in the Wnt/β-catenin signaling pathway. Biochem. Pharmacol. 2012, 84, 1143–1153. [Google Scholar] [CrossRef] [PubMed]
  43. Martin, P.M.; Yang, X.; Robin, N.; Lam, E.; Rabinowitz, J.S.; Erdman, C.A.; Quinn, J.; Weiss, L.A.; Hamilton, S.P.; Kwok, P.Y.; et al. A rare WNT1 missense variant overrepresented in ASD leads to increased Wnt signal pathway activation. Transl. Psychiatry 2013, 3, e301. [Google Scholar] [CrossRef] [PubMed]
  44. Levchenko, A.; Davtian, S.; Freylichman, O.; Zagrivnaya, M.; Kostareva, A.; Malashichev, Y. Beta-catenin in schizophrenia: Possibly deleterious novel mutation. Psychiatry Res. 2015, 228, 843–848. [Google Scholar] [CrossRef] [PubMed]
  45. Grünblatt, E.; Nemoda, Z.; Werling, A.M.; Roth, A.; Angyal, N.; Tarnok, Z.; Thomsen, H.; Peters, T.; Hinney, A.; Hebebrand, J.; et al. The involvement of the canonical Wnt-signaling receptor LRP5 and LRP6 gene variants with ADHD and sexual dimorphism: Association study and meta-analysis. Am. J. Med. Genet. 2019, 180, 365. [Google Scholar] [CrossRef] [PubMed]
  46. Zhao, H.; Nyholt, D.R. Gene-based analyses reveal novel genetic overlap and allelic heterogeneity across five major psychiatric disorders. Hum. Genet. 2017, 136, 263–274. [Google Scholar] [CrossRef] [PubMed]
  47. Wisniewska, M.B. Physiological Role of β-Catenin/TCF Signaling in Neurons of the Adult Brain. Neurochem. Res. 2013, 38, 1144–1155. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, J.; Xiao, Q.; Xiao, J.; Niu, C.; Li, Y.; Zhang, X.; Zhou, Z.; Shu, G.; Yin, G. Wnt/β-catenin signalling: Function, biological mechanisms, and therapeutic opportunities. Signal Transduct. Target. Ther. 2022, 7, 3. [Google Scholar] [CrossRef] [PubMed]
  49. Hajek, T.; Kopecek, M.; Höschl, C.; Alda, M. Smaller hippocampal volumes in patients with bipolar disorder are masked by exposure to lithium: A meta-analysis. J. Psychiatry Neurosci. 2012, 37, 333. [Google Scholar] [CrossRef]
  50. Emrich, H.M. Possible Role of Opioids in Mental Disorders. In Psychoneuroendocrine Dysfunction; Springer, Boston, MA, USA, 1984; pp. 293–307. [CrossRef]
  51. Nakamoto, K.; Tokuyama, S. Stress-Induced Changes in the Endogenous Opioid System Cause Dysfunction of Pain and Emotion Regulation. Int. J. Mol. Sci. 2023, 24, 11713. [Google Scholar] [CrossRef]
  52. Lutz, P.E.; Kieffer, B.L. Opioid receptors: Distinct roles in mood disorders. Trends Neurosci. 2013, 36, 195. [Google Scholar] [CrossRef]
  53. Qi, X.; Wen, Y.; Li, P.; Liang, C.; Cheng, B.; Ma, M.; Cheng, S.; Zhang, L.; Liu, L.; Kafle, O.P.; et al. An integrative analysis of genome-wide association study and regulatory SNP annotation datasets identified candidate genes for bipolar disorder. Int. J. Bipolar Disord. 2020, 8, 6. [Google Scholar] [CrossRef] [PubMed]
  54. Lindtröm, L.H.; Winderlöv, E.; Gunne, L.-M.; Wahlström, A.; Terenius, L. Endorphins in human cerebrospinal fluid: Clinical correlations to some psychotic states. Acta Psychiatr. Scand. 1978, 57, 153–164. [Google Scholar] [CrossRef] [PubMed]
  55. Genazzani, A.R.; Petraglia, F.; Facchinetti, F.; Monittola, C.; Scarone, S.; Brambilla, F. Opioid Plasma Levels in Primary Affective Disorders: Effect of Desimipramine Therapy. Neuropsychobiology 1984, 12, 78–85. [Google Scholar] [CrossRef] [PubMed]
  56. Fernstrom, J.D.; Fernstrom, M.H. Tyrosine, phenylalanine, and catecholamine synthesis and function in the brain. J. Nutr. 2007, 137, 1539S–1547S. [Google Scholar] [CrossRef]
  57. Daubner, S.C.; Le, T.; Wang, S. Tyrosine Hydroxylase and Regulation of Dopamine Synthesis. Arch. Biochem. Biophys. 2011, 508, 1. [Google Scholar] [CrossRef] [PubMed]
  58. Gamo, N.J.; Arnsten, A.F.T. Molecular Modulation of Prefrontal Cortex: Rational Development of Treatments for Psychiatric Disorders. Behav. Neurosci. 2011, 125, 282. [Google Scholar] [CrossRef]
  59. Klumpp, H.; Deldin, P. Review of brain functioning in depression for semantic processing and verbal fluency. Int. J. Psychophysiol. 2010, 75, 77–85. [Google Scholar] [CrossRef] [PubMed]
  60. Pomarol-Clotet, E.; Alonso-Lana, S.; Moro, N.; Sarró, S.; Bonnin, M.C.; Goikolea, J.M.; Fernández-Corcuera, P.; Amann, B.L.; Romaguera, A.; Vieta, E.; et al. Brain functional changes across the different phases of bipolar disorder. Br. J. Psychiatry 2015, 206, 136–144. [Google Scholar] [CrossRef] [PubMed]
  61. Lin, R.C.; Scheller, R.H. Mechanisms of Synaptic Vesicle Exocytosis. Annu. Rev. Cell Dev. Biol. 2000, 16, 19–49. [Google Scholar] [CrossRef]
  62. Ramakrishnan, N.A.; Drescher, M.J.; Drescher, D.G. The SNARE complex in neuronal and sensory cells. Mol. Cell Neurosci. 2012, 50, 58. [Google Scholar] [CrossRef] [PubMed]
  63. Cupertino, R.B.; Kappel, D.B.; Bandeira, C.E.; Schuch, J.B.; da Silva, B.S.; Müller, D.; Bau, C.H.D.; Mota, N.R. SNARE complex in developmental psychiatry: Neurotransmitter exocytosis and beyond. J. Neural Transm. 2016, 123, 867–883. [Google Scholar] [CrossRef] [PubMed]
  64. de Bartolomeis, A.; De Simone, G.; De Prisco, M.; Barone, A.; Napoli, R.; Beguinot, F.; Billeci, M.; Fornaro, M. Insulin effects on core neurotransmitter pathways involved in schizophrenia neurobiology: A meta-analysis of preclinical studies. Implic. Treat. Mol. Psychiatry 2023, 28, 2811. [Google Scholar] [CrossRef] [PubMed]
  65. Ferrario, C.R.; Finnell, J.E. Beyond the hypothalamus: Roles for insulin as a regulator of neurotransmission, motivation, and feeding. Neuropsychopharmacology 2023, 48, 232–233. [Google Scholar] [CrossRef]
  66. Blum, K.; Bowirrat, A.; Elman, I.; Baron, D.; Thanos, P.K.; Gold, M.S.; Hanna, C.; Makale, M.T.; Sunder, K.; Jafari, N.; et al. Evidence for the DRD2 Gene as a Determinant of Reward Deficiency Syndrome (RDS). Clin. Exp. Psychol. 2023, 9, 8. [Google Scholar] [PubMed]
  67. Huang, C.C.; Chang, Y.H.; Lee, S.Y.; Chen, S.L.; Chen, S.H.; Chu, C.H.; Huang, S.Y.; Tzeng, N.S.; Lee, I.H.; Yeh, T.L.; et al. The interaction between BDNF and DRD2 in Bipolar II disorder but not in bipolar i disorder. Am. J. Med. Genet. Part B Neuropsychiatr. Genet. 2012, 159B, 501–507. [Google Scholar] [CrossRef] [PubMed]
  68. Lee, S.Y.; Chen, S.L.; Chang, Y.H.; Chen, S.H.; Chu, C.H.; Huang, S.Y.; Tzeng, N.S.; Wang, C.L.; Lee, I.H.; Yeh, T.L.; et al. The ALDH2 and DRD2/ANKK1 genes interacted in bipolar II but not bipolar I disorder. Pharmacoge. Genom. 2010, 20, 500–506. [Google Scholar] [CrossRef]
  69. Sun, Q.; Yuan, F.; Yuan, R.; Ren, D.; Zhu, Y.; Bi, Y.; Hu, J.; Guo, Z.; Xu, F.; Niu, W.; et al. GRIK4 and GRM7 gene may be potential indicator of venlafaxine treatment reponses in Chinese of Han ethnicity. Medicine 2019, 98, e15456. [Google Scholar] [CrossRef] [PubMed]
  70. Knight, H.M.; Walker, R.; James, R.; Porteous, D.J.; Muir, W.J.; Blackwood, D.H.R.; Pickard, B.S. GRIK4/KA1 protein expression in human brain and correlation with bipolar disorder risk variant status. Am. J. Med. Genet. Part B Neuropsychiatr. Genet. 2012, 159B, 21–29. [Google Scholar] [CrossRef] [PubMed]
  71. Pickard, B.S.; Knight, H.M.; Hamilton, R.S.; Soares, D.C.; Walker, R.; Boyd, J.K.F.; Machell, J.; Maclean, A.; McGhee, K.A.; Condie, A.; et al. A common variant in the 3′UTR of the GRIK4 glutamate receptor gene affects transcript abundance and protects against bipolar disorder. Proc. Natl. Acad. Sci. USA 2008, 105, 14940–14945. [Google Scholar] [CrossRef]
  72. Li, T.; Forbes, M.E.; Fuller, G.N.; Li, J.; Yang, X.; Zhang, W. IGFBP2: Integrative hub of developmental and oncogenic signaling network. Oncogene 2020, 39, 2243. [Google Scholar] [CrossRef] [PubMed]
  73. Bezchlibnyk, Y.B.; Xu, L.; Wang, J.F.; Young, L.T. Decreased expression of insulin-like growth factor binding protein 2 in the prefrontal cortex of subjects with bipolar disorder and its regulation by lithium treatment. Brain Res. 2007, 1147, 213–217. [Google Scholar] [CrossRef] [PubMed]
  74. Fernandez, A.M.; Torres-Alemán, I. The many faces of insulin-like peptide signalling in the brain. Nat. Rev. Neurosci. 2012, 13, 225–239. [Google Scholar] [CrossRef] [PubMed]
  75. Mikhalitskaya, E.V.; Vyalova, N.M.; Ermakov, E.A.; Levchuk, L.A.; Simutkin, G.G.; Bokhan, N.A.; Ivanova, S.A. Association of Single Nucleotide Polymorphisms of Cytokine Genes with Depression, Schizophrenia and Bipolar Disorder. Genes 2023, 14, 1460. [Google Scholar] [CrossRef]
  76. de Marco, A.; Scozia, G.; Manfredi, L.; Conversi, D. A Systematic Review of Genetic Polymorphisms Associated with Bipolar Disorder Comorbid to Substance Abuse. Genes 2022, 13, 1303. [Google Scholar] [CrossRef] [PubMed]
  77. Wang, K.S.; Liu, X.F.; Aragam, N. A genome-wide meta-analysis identifies novel loci associated with schizophrenia and bipolar disorder. Schizophr. Res. 2010, 124, 192–199. [Google Scholar] [CrossRef] [PubMed]
  78. Kohshour, M.O.; Papiol, S.; Ching, C.R.K.; Schulze, T.G. Genomic and neuroimaging approaches to bipolar disorder. BJPsych Open 2022, 8, e36. [Google Scholar] [CrossRef] [PubMed]
  79. Nadal, R. Pharmacology of the Atypical Antipsychotic Remoxipride, a Dopamine D2 Receptor Antagonist. CNS Drug Rev. 2001, 7, 265. [Google Scholar] [CrossRef]
  80. Simmons, B.; Kuo, A. Analgesics, Tranquilizers, and Sedatives, Cardiac Intensive Care; Elsevier: Amsterdam, The Netherlands, 2019; pp. 421–431.e5. [Google Scholar] [CrossRef]
  81. Mecasermin for primary insulin-like growth factor-1 deficiency. Aust. Prescr. 2022, 45, 215. [CrossRef]
  82. Scheller, E.L.; Krebsbach, P.H. Gene Therapy: Design and Prospects for Craniofacial Regeneration. J. Dent. Res. 2009, 88, 585. [Google Scholar] [CrossRef]
  83. Chanchal, D.K.; Chaudhary, J.S.; Kumar, P.; Agnihotri, N.; Porwal, P. CRISPR-Based Therapies: Revolutionizing Drug Development and Precision Medicine. Curr. Gene Ther. 2024, 24, 193–207. [Google Scholar] [CrossRef] [PubMed]
  84. Madigan, V.; Zhang, F.; Dahlman, J.E. Drug delivery systems for CRISPR-based genome editors. Nat. Rev. Drug Discov. 2023, 22, 875–894. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Dendrogram based on the hierarchal clustering of genes from which the modules that display high gene interconnectivity were identified.
Figure 1. Dendrogram based on the hierarchal clustering of genes from which the modules that display high gene interconnectivity were identified.
Biology 13 00787 g001
Figure 2. Network heatmap of the identified modules. Darker colors indicate greater connectivity; positive correlation means strong association among the genes in that module; otherwise, negative correlation suggests weak association. Moreover, the number of asterisks indicates the likelihood that the expression patterns identified in that module greatly characterize that brain region.
Figure 2. Network heatmap of the identified modules. Darker colors indicate greater connectivity; positive correlation means strong association among the genes in that module; otherwise, negative correlation suggests weak association. Moreover, the number of asterisks indicates the likelihood that the expression patterns identified in that module greatly characterize that brain region.
Biology 13 00787 g002
Figure 3. Bar plots of the top clusters from each annotation of the 7 modules identified to have significant expression among the three brain regions.
Figure 3. Bar plots of the top clusters from each annotation of the 7 modules identified to have significant expression among the three brain regions.
Biology 13 00787 g003
Figure 4. Bubble plot visualization of the top cluster in the KEGG pathway analysis based on their respective fold enrichment score.
Figure 4. Bubble plot visualization of the top cluster in the KEGG pathway analysis based on their respective fold enrichment score.
Biology 13 00787 g004
Figure 5. Modified PPI network of DRD2 constructed from the expression profiles of BPD patients, showing the interactions of the gene with its first-degree neighbors.
Figure 5. Modified PPI network of DRD2 constructed from the expression profiles of BPD patients, showing the interactions of the gene with its first-degree neighbors.
Biology 13 00787 g005
Table 1. Identified genes among the three brain regions which showed the presence of disease-associated variants.
Table 1. Identified genes among the three brain regions which showed the presence of disease-associated variants.
SampleGeneVariantType
nAccDRD2rs1801028missense variant
GFRA2rs7833426intron variant
DCBLD1rs62433108intron variant
AnCgST8SIA2rs4777989intron variant
ADAMTS16rs16875288intron variant
DLPFCFOXO3rs1536057intron variant
rs1935952
rs2802292
ITGA9rs166508intron variant
CUBNrs7904579intron variant
PLCB4rs2299682intron variant
RORBrs1327836intron variant
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Garcia, J.P.T.; Tayo, L.L. Codes between Poles: Linking Transcriptomic Insights into the Neurobiology of Bipolar Disorder. Biology 2024, 13, 787. https://doi.org/10.3390/biology13100787

AMA Style

Garcia JPT, Tayo LL. Codes between Poles: Linking Transcriptomic Insights into the Neurobiology of Bipolar Disorder. Biology. 2024; 13(10):787. https://doi.org/10.3390/biology13100787

Chicago/Turabian Style

Garcia, Jon Patrick T., and Lemmuel L. Tayo. 2024. "Codes between Poles: Linking Transcriptomic Insights into the Neurobiology of Bipolar Disorder" Biology 13, no. 10: 787. https://doi.org/10.3390/biology13100787

APA Style

Garcia, J. P. T., & Tayo, L. L. (2024). Codes between Poles: Linking Transcriptomic Insights into the Neurobiology of Bipolar Disorder. Biology, 13(10), 787. https://doi.org/10.3390/biology13100787

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop