Differential Expression of Proteins Associated with Bipolar Disorder as Identified Using the PeptideShaker Software

The prevalence of bipolar disorder (BD) in modern society is growing rapidly, but due to the lack of paraclinical criteria, its differential diagnosis with other mental disorders is somewhat challenging. In this regard, the relevance of proteomic studies is increasing due to the development of methods for processing large data arrays; this contributes to the discovery of protein patterns of pathological processes and the creation of new methods of diagnosis and treatment. It seems promising to search for proteins involved in the pathogenesis of BD in an easily accessible material—blood serum. Sera from BD patients and healthy individuals were purified via affinity chromatography to isolate 14 major proteins and separated using 1D SDS-PAGE. After trypsinolysis, the proteins in the samples were identified via HPLC/mass spectrometry. Mass spectrometric data were processed using the OMSSA and X!Tandem search algorithms using the UniProtKB database, and the results were analyzed using PeptideShaker. Differences in proteomes were assessed via an unlabeled NSAF-based analysis using a two-tailed Bonferroni-adjusted t-test. When comparing the blood serum proteomes of BD patients and healthy individuals, 10 proteins showed significant differences in NSAF values. Of these, four proteins were predominantly present in BD patients with the maximum NSAF value: 14-3-3 protein zeta/delta; ectonucleoside triphosphate diphosphohydrolase 7; transforming growth factor-beta-induced protein ig-h3; and B-cell CLL/lymphoma 9 protein. Further exploration of the role of these proteins in BD is warranted; conducting such studies will help develop new paraclinical criteria and discover new targets for BD drug therapy.


Introduction
Mood disorders, such as bipolar disorder (BD), affect up to 3% of the world's population [1].They negatively impact the professional and social aspects of patients' lives, leading to a significant decrease in quality of life and an increase in healthcare costs [2,3].In addition, patients with BD have an increased risk of suicide (10-15%) [4].Therefore, early diagnosis of BD and specific treatment will help prevent potential suicidal behavior in patients and improve their quality of life [5].
However, diagnosing BD remains challenging.The literature provides strong evidence in favor of the concept of a broad spectrum of bipolar disorders [6].Many of the symptoms of bipolar spectrum disorders may also result from other conditions, including anxiety disorders, psychotic disorders, substance use disorders, personality disorders, etc. [7][8][9].In addition, cardiovascular disease, diabetes, and obesity are also strongly associated with BD [10,11].This association leads to excess cardiovascular mortality and suicide, with a loss of approximately 10-20 potential life years [9].Delayed diagnosis and, as a result, incorrect prescription of therapy leads to worsening of the course of the disease, cognitive impairment, and reduced social functioning [12,13].J.D. Lish et al. state that 73% of patients with BD are initially misdiagnosed, and it takes 8 years on average for patients to receive a correct diagnosis.The authors also note that in about 60% of all patients, the first episode manifests in childhood or adolescence, yet more than 50% of patients do not receive specific therapy for the next 5 years or more [14].Approximately 75% of cases of BD manifest in late adolescence or early adulthood [15], and the disease is the fourth leading cause of the global burden of disease in adolescents and young adults [16].
Differential diagnosis of unipolar depression with BD is particularly challenging, and the condition is commonly underdiagnosed as a result.A total of 5-20% of patients were initially diagnosed with recurrent depression and subsequently received a revised diagnosis of BD [17].Studies have revealed that 17% of patients observed by general practitioners and receiving prescribed long-term maintenance antidepressant therapy actually suffer from bipolar disorder.In a study of 600 patients with BD, it was found that 69% of patients initially received a different diagnosis, which changed on average about three times (the first diagnosis in 60% of patients was recurrent depression, and the second was anxiety disorder), and it took about 10 years for them to be diagnosed correctly [18].
Difficulties in diagnosing BD may be attributed to the fact that, at present, the diagnosis is made through a comprehensive clinical assessment solely based on clinical symptoms and anamnestic and follow-up data, which remain descriptive and generally do not have sufficient predictive validity [19,20].This is due to a lack of understanding of the pathogenetic mechanisms of BD and, as a result, the absence of paraclinical markers.
Recently, highly effective methods for studying mental disorders have been developed, in particular, proteomic methods [21,22].Proteomic analysis allows the detection of protein biomarkers associated with functional disorders involved in the pathophysiology of diseases without the need to put forward a hypothesis and limit its search area [23].In the future, this approach could be used to help identify biomarkers of affective disorders and risks of predisposition to developing the disease later in life [24,25].
There are few published works that explore the proteomics of schizophrenia and BD, but most studies have been conducted on schizophrenia and post-mortem brain tissues [26][27][28][29][30][31].For clinical practice, studies of biological fluids are more popular, according to which various physiological states of organs and systems, including the brain, are reflected in their composition [32,33].Therefore, it is important to study potential biomarkers of mental disorders in an easily accessible material-blood serum [34][35][36][37][38][39][40][41].However, all discovered proteins do not solve the problem of differential diagnosis of mental disorders due to their non-specificity.Thus, it is possible that research should focus not on the search for a specific protein but on identifying differences in the content of a set of proteins (panels) that reflect the main pathogenetic mechanisms.
Until recently, studies on the search for biomarkers of mental disorders have mainly used qualitative proteomics methods that reveal only qualitative differences in the composition of proteomes, that is, the presence or absence of specific proteins in the proteomes of the studied samples.For this reason, in recent years, quantitative proteomics has become the preferred approach for detecting proteome differences in disease [42].However, quantitative proteomics methods using radionuclide labels are expensive, making them inconvenient for screening studies.At this point, label-free approaches may be preferable because they are simpler, more replicable, and less expensive than label-based approaches [43].In addition, label-free semi-quantification does not limit the number of samples and conditions compared, making it suitable for longitudinal and clinical proteomics [44].
In this article, we searched for potential BD-specific biomarkers in serum by label-free semi-quantitative proteomics using the PeptideShaker version 2.2.9 software package.This study may be useful to future research on pathogenetic mechanisms, the development of new laboratory diagnostic criteria, and drug targets for BD and other affective disorders.

Results
In this work, mass spectrometric data of the blood serum proteome of patients with bipolar disorder (BD) and healthy individuals were analyzed in a comparative aspect.Us-ing the SearchGUI version 3.3.12software package and the OMSSA and X!Tandem search algorithms, protein lists were obtained for each examined patient and healthy volunteer.Tables of reliably identified proteins were formed from these lists using the PeptideShaker program.Next, the intensity values of the NSAF of reliably identified proteins in all the samples were analyzed.After comparing the proteomes of the examined individuals using a two-tailed unpaired Student's t-test (FDR 0.05 and S0 = 2), proteins whose peak intensities were statistically significantly different in patients with BD and healthy individuals were identified.A list of these proteins is presented in Table 1, indicating the statistical differences according to Student's t-tests, along with their NSAF values for the BD group and control patients.The identified proteins are involved in many biological processes, including those that play an important role in the development and functioning of the central nervous system.Most of these proteins are involved in the regulation of cell adhesion, organization of the extracellular matrix, growth and migration of endothelial cells, the regulation of angiogenesis in response to a stimulus of vascular endothelial growth factor, and endothelial cell apoptosis.Proteins are also involved in the regulation of the radial mobility of neurons in the cerebral cortex controlled by glia and the growth of neuronal processes, in the regulation of cytoskeleton assembly, binding of actin filaments, and in the transmission of cell signals in glutamatergic synapses.These proteins are also involved in the processes of ubiquitination and protein catabolism.They regulate protein phosphorylation and dephosphorylation by mediating the activity of G-protein-coupled receptors; regulate the activity of GTPase, deoxycytidine triphosphatase, MAP kinase and GDP phosphatase, the ERK1 and ERK2 cascade, and other protein kinase phosphorylation enzymes; and regulate the activity of nucleoside diphosphate phosphatase and ribonucleoside triphosphate phosphatase.They participate in the regulation of cell morphogenesis, transcription, proliferation and cell differentiation, and regulation of the mitotic cycle and regulate the morphogenesis of neurons, the formation of layers in the cerebral cortex, the development of dendrites, the formation and maturation of synapses, axonogenesis, and synaptic plasticity.The identified proteins are involved in intracellular transport and the regulation of protein and vesicle transport, autophagosome assembly, transmembrane transport, the intracellular position of the Golgi complex, its assembly, and the maintenance of its structural integrity.They are also involved in the directed movement of vesicles inside the Golgi apparatus and intracellular protein transport and transport of organelles through microtubules.In addition to these functions, proteins are involved in the regulation of the immune response and differentiation of T-helpers, the regulation of cytokine production, and the apoptotic process.
Figure 1 shows a functional diagram of the relationship between the genes of the identified proteins obtained using the GeneMANIA tool [45].
transport of organelles through microtubules.In addition to these functions, proteins are involved in the regulation of the immune response and differentiation of T-helpers, the regulation of cytokine production, and the apoptotic process.
Figure 1 shows a functional diagram of the relationship between the genes of the identified proteins obtained using the GeneMANIA tool [45].
Figure 1.Functional network of protein gene connections that showed statistically significant differences between groups of patients with bipolar disorder and healthy individuals.Line colors indicate the following: lilac-co-expression, pink-physical interaction, orange-predicted, green-genetic interaction.The thickness of the line reflects the strength of the interaction.The tool GeneMANIA (genemania.org) was used.
The figure shows the functional gene communication network of significantly expressed proteins.The network is built based on data from studies of gene and protein expression profiles, as well as their molecular interaction pathways, from various sources such as GEO, BioGRID, Ensembl, NCBI, MGI, Pathway Commons, and others.For a complete list of GeneMANIA sources, see http://genemania.org/,accessed on 13 July 2023.The association of genes in a data set is represented by the Pearson correlation coefficient, which reflects the strength of the interaction or the reliability of the observation that they interact, and it is color-coded for the different measures.
Co-expression reflects gene expression data.Two genes are functionally related if their expression levels change equally when the conditions of the gene expression study Functional network of protein gene connections that showed statistically significant differences between groups of patients with bipolar disorder and healthy individuals.Line colors indicate the following: lilac-co-expression, pink-physical interaction, orange-predicted, green-genetic interaction.The thickness of the line reflects the strength of the interaction.The tool GeneMANIA (genemania.org) was used.
The figure shows the functional gene communication network of significantly expressed proteins.The network is built based on data from studies of gene and protein expression profiles, as well as their molecular interaction pathways, from various sources such as GEO, BioGRID, Ensembl, NCBI, MGI, Pathway Commons, and others.For a complete list of GeneMANIA sources, see http://genemania.org/,accessed on 13 July 2023.The association of genes in a data set is represented by the Pearson correlation coefficient, which reflects the strength of the interaction or the reliability of the observation that they interact, and it is color-coded for the different measures.
Co-expression reflects gene expression data.Two genes are functionally related if their expression levels change equally when the conditions of the gene expression study are changed.The expression of such genes is probably regulated by the same mechanisms or by similar or interrelated pathways.For our proteins, the level of co-expression is 54.4%.Physical interaction reflects data on the protein-protein interactions of the products of related genes.For the proteins presented in the figure, the physical interaction is 26.1%.Predicted interaction, a similar indicator to physical interaction, reflects functional connections between genes, mainly protein-protein interactions, based on similar functional connections in another organism.For the proteins we identified, this indicator is 12.67%.Genetic interaction indicates a functional relationship between two genes if the level of expression of one gene entails a change in the expression of another gene.For the proteins shown in the figure, genetic interaction is quite low and amounts to 2.5%.

Discussion
In our work, we used the PeptideShaker software package, which uses several search engines on one data set at once using the target-decoy search strategy [45] for probabilistic evaluation of data processing errors.The use of this strategy to unify the peptide and spectrum matching lists (PSM) of various search engines makes it possible to increase the reliability and sensitivity of peptide detection compared with data processing methods using just one search engine [46,47].PeptideShaker also delivers false discovery rates (FDRs) of peptides and proteins at the PSM level and allows relative quantification of peptides in different samples via label-free quantitation (LFQ) based on peak intensity or area under the curve of the detected peptide [48,49].
After using the PeptideShaker software and the statistical processing of the results, 10 differentially expressed proteins that showed significant differences between the groups of patients with bipolar disorder and healthy individuals were identified.Previously, most of these proteins were not associated with the pathogenesis of bipolar disorder.Therefore, below, we provide a description of the already known physiological and pathophysiological functions of these proteins so that, on the basis of this, we can make an assumption about their possible participation in the mechanisms of the pathogenesis of affective disorders.
B-cell CLL/lymphoma 9 protein (BCL9) was the most represented by the NSAF value in the BD patients in our study.This protein is involved in signal transmission through the Wnt pathway, one of the intracellular signaling pathways that regulates embryogenesis and cell differentiation [50].Wnt signals via β-catenin and lymphoid enhancer factor 1/T-cellspecific transcription factors (LEF1/TCFs).There is ample evidence linking Wnt/β-catenin signaling to mood disorders, but the mechanism of this association is still unknown.The BCL9 gene encoding the protein we found was associated with major depressive disorder.Additionally, the classical mood stabilizer lithium can activate β-catenin by inhibiting GSK3α/β 6, which is a key enzyme in the canonical Wnt 7 pathway [51].The Wnt pathway is a major regulator of adult hippocampal neurogenesis.Upregulation of Wnt3 enhances neurogenesis in the hippocampus in vitro and in vivo, and blockade of Wnt signaling reduces neurogenesis [52].The genes WNT2B and TCF7L2 have also been linked with BD [53,54].In addition, astrocyte-produced Wnt is required to maintain the structural and functional integrity of BBB endothelial cells [55,56].Thus, the protein we identified can presumably also play a role in the regulation of BBB in mental pathology.This protein can apparently also be involved in the pathogenesis of bipolar disorder, perhaps by analogy with the same mechanisms that are involved in the pathogenesis of major depressive disorder.This fact will not allow its use for the differential diagnosis of affective disorders but will undoubtedly contribute to the discovery of the mechanisms of the pathogenesis of bipolar disorder.
The second-largest NSAF in bipolar patients in our study was 14-3-3 protein zeta/delta (YWHAZ).The 14-3-3 proteins are a family of highly conserved, multifunctional proteins that are primarily expressed in the brain and are closely associated with various brain disorders, although their exact neurophysiological function is not yet fully understood [57,58].The main function of 14-3-3 proteins is their ability to close phosphorylated regions of target proteins, thus blocking the action of phosphatases and preventing their dephosphorylation [59].Proteins of the 14-3-3 family can also change the conformation and subcellular localization of their target proteins [60][61][62][63].
In neurons, proteins of the 14-3-3 family are present in the cytoplasmic reticulum, intracellular organelles, and the plasma membrane, where they play a functional role in cellular processes such as the differentiation, migration, and survival of neurons; neurite outgrowth; and regulation of ion channels [64][65][66].Some of the 14-3-3 isoforms are especially abundant in synapses, where they regulate signaling and neuroplasticity [67][68][69].The isoform found in our study, 14-3-3 ζ (14-3-3 protein zeta/delta), which is maximally represented in patients with bipolar disorder, is an activator of tyrosine and tryptophan hydroxylases, enzymes that limit the rate of synthesis of dopamine and other neurotransmitters [70][71][72].The isoforms 14-3-3 ε and ζ are important for neuronal development [73].The underlying mechanism involves the Ndel1/LIS1/14-3-3 protein complex, which is critical to proper neuronal migration and axon growth.Disturbances in the axonal transport of this protein complex lead to abnormal migration and development of neurons, which is associated with diseases of the nervous system, such as schizophrenia [74,75].It was also found that gene knockout 14-3-3 ζ in mice leads to morphological changes in the brain; these mice have enlarged lateral ventricles, aberrant connectivity of mossy fibers, and reduced density of synapses and dendritic spines of the hippocampus [76,77].These studies provide evidence of involvement 14-3-3 ζ in the formation of dendritic spines.Other genetic and post-mortem mRNA studies have identified genes for the 14-3-3 isoforms-including the gene YWHAZ, the encoding of which was discovered by us in the 14-3-3 ζ protein-as a potential susceptibility gene for schizophrenia [78][79][80][81].A significant decrease in mRNA expression was also identified for 14-3-3 (isoforms β, η, ε, σ, θ, ζ) in the prefrontal cortex and for 14-3-3η in the cerebellum of patients with schizophrenia [82,83].
In addition, chromosomal regions have been identified that share common candidate genes for the risk of BD and schizophrenia.One such associated area is 22q12-13, which contains the gene 14-3-3η (YWHAH) [84][85][86].For example, one meta-analysis found that YWHAH has a statistically significant relationship with BD [87].In addition to 14-3-3η, other 14-3-3 isoforms are involved in BD, so reduced levels of mRNA expression have been reported in 14-3-3ε, -σ, and -ζ in brain samples of patients with BD [88].
Thus, all of the above works were performed on postmortem genetic material, theoretically proving the role of 14-3-3ζ protein in the pathogenesis of mental disorders.The potential detection of an increase in the amount of 14-3-3 protein zeta/delta in the blood serum of patients with bipolar disorder makes it a promising biomarker for BD.
For transforming growth factor-beta-induced protein ig-h3 (TGFBI), the NSAF was three times higher in patients with BD compared to that in healthy individuals.It is an extracellular matrix (ECM) protein that is involved in many physiological processes, including morphogenesis, cell adhesion and migration, angiogenesis, and inflammation [89].The TGFBI protein connects various molecules of the extracellular matrix with each other and promotes the interaction of cells with collagen, fibronectin [90], various integrins [91], and proteoglycans such as biglycan and decorin [92].TGFBI also triggers phosphorylation and activates several intracellular pathways, including protein kinase AKT1 (a key enzyme in the PI3K/AKT signaling pathway), focal adhesion kinase, and paxillin [93].TGFBI is distributed in the extracellular matrices of a wide range of tissues and plays a role in the adhesion and migration of a wide range of cells, including endothelial cells [94].Additionally, TGFBI expression is induced in reactive astrocytes of the rat cerebral cortex at sites of injury [95].It appears that the detection of this protein indicates the activation of protective mechanisms in response to damage to the vascular endothelium, which may, in turn, indicate the modulation of the permeability of the BBB in BD.Until now, this protein has not been associated with psychiatric disorders.However, the activation of neuroinflammation and impaired BBB permeability in affective mental disorders has recently attracted more attention from researchers.
A significant difference in the NSAF value in patients with bipolar disorder and healthy individuals was also shown for ectonucleoside triphosphate diphosphohydrolase 7 (ENTPD7).This protein is an enzyme that catalyzes the hydrolysis of nucleoside triphosphates and diphosphates in the presence of calcium or magnesium and predominantly hydrolyzes nucleoside-5 -triphosphates.It also hydrolyzes ATP and nucleoside diphosphates to a lesser extent [96].The functional role of this protein in the human body has not been studied practically, but it can be argued with a high degree of probability that this protein hydrolyzes nucleotides in the body and participates in important regulatory processes in purine metabolism, which plays a key role in oxidative stress, DNA damage, and cellular aging [97,98].
Until now, this protein has not been associated with mental disorders, but according to the literature, oxidative stress plays a significant role in the pathogenesis of bipolar disorder.Therefore, based on the results presented, it is very likely that the protein produced by the ENTPD7 gene also contributes to the development of BD.
The remaining six proteins, which showed a statistically significant difference in the NSAF value in patients with bipolar disorder and healthy individuals, are more frequently represented in healthy individuals and in patients with minimal peak values.These are the following proteins.
Coiled-coil domain-containing protein 80 (CCDC80) and TGFBI can be involved in cell adhesion and extracellular matrix assembly; however, its mechanism of action and the role of this protein in cell adhesion have not yet been studied [99].However, its decrease in patients may also lead to an increase in BBB permeability.
Disabled homolog 2-interacting protein (DAB2IP) is a member of the Ras-GTPaseactivating protein family and is a scaffold protein involved in the regulation of a wide range of both general and specialized signaling pathways.It is involved in processes such as innate immune response, inflammation and cell growth inhibition, apoptosis, cell survival, angiogenesis, cell migration, and maturation, and also plays a role in cell cycle regulation [100].
This protein is normally highly expressed in the vascular endothelium and functions as an inhibitor of angiogenesis and endothelial cell migration by blocking VEGFR-2 activity, so DAB2IP deficiency leads to increased VEGFR-2 mediated angiogenesis [101].However, DAB2IP gene knockout reduces corneal and retinal angiogenesis, which is associated with positive regulation of DAB2IP endocytosis and VEGFR-3 expression [102].Under stress of the endoplasmic reticulum caused by the accumulation of unfolded (misfolded) proteins, DAB2IP, with the participation of TNF, mediates the activation of a signaling pathway involving the enzyme apoptosis signal-regulating kinase 1 (ASK1) [103].Activation of this mechanism under endoplasmic reticulum stress has also been described for endothelial cells [104].TNF, together with RIP1, induces phosphorylation of DAB2IP, which in turn causes dissociation of ASK1 with its inhibitory proteins 14-3-3; this activates the signal transduction cascade for endothelial cell apoptosis [105][106][107].
DAB2IP is also highly expressed in the brain and interacts with the disabled-1 (DAB1) protein, a key mediator of the reelin pathway that controls the migration and position of neurons during development.DAB1 plays a role in neuronal migration and the growth of neurites and their processes in the developing neocortex [108] and regulates dendritic development and synapse formation in the developing cerebellum [109].In addition, regulation of the DAB2IP gene via methylation is key to the differentiation and maturation of the central nervous system [110].A decrease in the expression of this protein may be the key to the development of psychiatric pathology.
One other protein, adhesion G-protein-coupled receptor B1 (ADGRB1), is part of the G-protein-coupled family of membrane receptor protein (GPCR), which allows the nervous system to accurately respond to neurotransmitters, nucleotides, amines, peptides, cytokines and hormones, and transmit specific signals across the cell membrane [111,112].The inherent ligand selectivity of neuronal GPCRs ensures proper integration between signal transduction pathways.Thus, about 35% of the drugs included in the FDA list act through the GPCR [113,114].The role of GPCR in psychiatric illnesses, including schizophrenia, BD, and depression, has been discussed [115].
The ADGRB1 we discovered has a large extracellular domain (ECD) that is responsible for adhesive function [116].This class of proteins is actively involved in the early development of the nervous system and brain [117].The protein in the GPCR receptor allows nerve cells to transmit signals between themselves and their microenvironment and migrate to destinations to perform specific functions.Thus, it has been demonstrated that in mouse Purkinje neurons, GPCR adhesion is necessary for the formation of complex dendritic structures for synaptic connections [118].AGRB1 regulates synaptogenesis by controlling the recruitment of the Par3/Tiam1 polarity complex to synaptic sites [119].In the adult, ADGRB1 regulates synaptic plasticity in learning and memory in the hippocampus [120].ADGRB1 is also involved in the formation of postsynaptic receptors that control the development of excitatory synapses [121,122].Despite the fact that the connection of this neurospecific protein with mental disorders has already been established, its role in the pathogenesis of BD has yet to be elucidated.However, the study of its pathophysiological mechanisms is deemed to be extremely promising for the development of new drugs and methods of treatment and diagnosis of BD.
The next protein, Sec1 family domain-containing protein 1 (SCFD1), is a member of the Senc1/Munk 18 protein family (SM).These are vesicle-traversing proteins that interact with the SNARE integral membrane proteins that mediate the fusion of vesicles with membrane-bound compartments [123].
SCFD1 is mainly involved in transport from the endoplasmic reticulum to the Golgi apparatus by promoting membrane fusion and allowing vesicles to pass from one compartment to another [124,125].Moreover, it plays a significant protective role in response to oxidative stress, participating in the mechanisms of protection of the endoplasmic reticulum from morphological changes associated with oxidative stress and preventing cell death.SCFD1 has been proven to be closely associated with the pathogenesis of Parkinson's disease and amyotrophic lateral sclerosis (ALS) [126].This is how the association of the rs10139154 polymorphism of the SCFD1 gene with the risk of developing ALS is described [127].The mechanisms underlying possible neurodegeneration have not yet been fully elucidated, although it has been demonstrated that ALS patients have impaired endosomal transport function and dysfunction of the Golgi network [128].These effects are exacerbated under conditions of stress, leading to neurodegeneration [129,130].In addition, in Alzheimer's disease and ALS, impairment of endocytic migration was revealed [131,132].In these neurodegenerative diseases, pathogenetic mechanisms involving SCFD1 have also been identified, leading to a disruption in the mechanism of axonal transport of proteins to the synaptic ending, which is regulated by neurofilaments [133][134][135].Thus, until now, this protein has been associated with neurodegenerative diseases, and its connection with mental disorders has been identified for the first time in this study.Based on the mechanisms described above, a decrease in its amount in patients with bipolar disorder causes a decrease in protection against OS, which in turn triggers a vicious circle (pathway): oxidative stress-inflammation-oxidative stress.This presumably leads to both the development of chronic affective disorders and the development of neurodegenerative diseases.
Ras GTPase-activating-like protein IQGAP1 (IQGAP1) is a member of the Ras superfamily of small GTPases [136], a large family of monomeric G proteins with structural and functional homology to the Gα subunits of heterotrimeric G proteins.These proteins play a crucial role in regulating the dynamics and assembly of the actin cytoskeleton in the leading processes of migrating neurons, directing and coordinating changes in both actin filaments and microtubules.IQGAP1 binds to and alters the function of several proteins, including actin, E-cadherin, beta-catenin, Cdc42, phospholipase C epsilon 1 (PLCe1), Lissencephaly-1 homolog (Lis1), and Ras-related C3 botulinum toxin substrate 1 (Rac1).
Lis1 regulates the activity of motor protein dynein by mediating microtubule glide during nuclear migration during cell division [137,138].Lis1 is known to be involved in the localization of GABAergic synaptic vesicles in the ventral nerve trunk [139] and the differentiation of neurons [140].Defects in the Lis1 gene disrupt neuronal migration, causing severe lissencephaly, a complex disorder of brain function caused by anomalies in neurogenesis.Although the role of Lis1 in this process is not fully understood, Lis1 deficiency is known to cause dysregulation of the Ras GTPases Cdc42, Rac1, and RhoA family proteins and subsequent cytoskeletal actin defects.Here, we note that IQGAP1 may be actively involved in the reorganization of the actin cytoskeleton.IQGAP1 is a substrate for the PLCe1 protein kinase, which catalyzes the phosphorylation of Ser-1443 to promote neurite outgrowth [141,142].Thus, IQGAP1 may play a role in cell cycle regulation by promoting neurite outgrowth, neurogenesis, neuronal migration, and dendritic outgrowth.A decrease in the expression of this protein may be associated with the development of bipolar disorder.
Coatomer subunit gamma-1 (COPG1) is a cytosolic protein that, as part of the coatomer complex (COPI), binds to dilysin proteins and Golgi vesicles that are not covered with clathrin.COPI is required for the budding of vesicles from the Golgi membranes and the transport of dilysin-labeled proteins from the endoplasmic reticulum and also affects the structural integrity of the Golgi apparatus and recycling of LDL receptors on the plasma membrane [143,144].
The role of the COPG1 protein in the development of mental illness has not yet been studied.It has been reported that the COPI complex may play a role in Alzheimer's disease, as it is involved in the localization of the amyloid precursor protein, and downregulation of the COPI delta subunit (δ-COP) expression leads to a decrease in amyloid plaques and improved memory in mice [145].It has also been shown in a cell model on mouse neuronal progenitors that Copg1, an analog of human COPG1, is required for the proper differentiation and maturation of neurons.Impairment of this gene leads to morphological changes in the Golgi apparatus and disruption of the formation of neuronal processes [146].Reduced expression of this protein significantly reduces neuroplasticity and, under certain conditions, likely leads to the development of bipolar disorder.
Eight proteins presented in our study were shown to be associated with BD for the first time and were not previously associated with the pathogenesis of BD.Two of the previously identified proteins with significant differences had already been shown to be associated with BD.All ten differentially expressed proteins were found to be neurospecific to one degree or another.Figure 1 shows a functional diagram of the relationships between the genes of the identified proteins, showing a large number of different connections.The proteins we identified showed a high level of co-expression of 54.4%, which indicates a close functional connection between the proteins and, most likely, mutual influence in the pathway of bipolar disorder.
Thus, the use of the PeptideShaker software package in this work made it possible to conduct a semi-quantitative analysis of the protein spectrum without using a labeled standard and to approximately estimate the protein levels of comparison groups without using expensive quantitative proteomics methods.From this work, it was concluded that PeptideShaker is an easy-to-use tool for analyzing mass spectrometry data at the initial exploratory stage of a study.The tool is suitable for both the analysis and interpretation of primary data and for the re-analysis and comparison of data from different experiments obtained over different periods of time.The most important reason we chose this tool is the ability to work with the combined results of several identification algorithms on a single set of input data.In this case, the program uses spectrum files in the general Mascot format (MGF) (which can be obtained from any MS spectrum file) and allows one to load custom databases of protein sequences in the FASTA file format and analyze experiments performed over a long time under the same conditions using one protein database.Also, to increase the reliability of search results, PeptideShaker uses a "target-bait" approach, in which the protein database is supplemented with false sequences; in our case, the reverse sequences of all the proteins were used, which reduced the rate of false discoveries.

Patient Characteristics
In the present work, patients with bipolar disorder treated in the Affective States Department of Mental Health Research Institute of Tomsk National Research Medical Center were examined.This study used a cross-sectional design.The group of patients with BD included 8 individuals.The patients' age (Me [Q25; Q75]) was 44 [38; 48] years, and the disease duration was 8 [4; 12] years.The age of disease manifestation was 33 [22; 45] years; the average number of experienced episodes was 4 [2; 5].Most of the patients were treated for a current depressive episode of varying severity with diagnoses according to ICD-10: F31.3-4 individuals, F31.5-2 individuals.Additionally, this group included patients with a current mixed affective episode (F31.6-2individuals).
The criteria for inclusion in the group of examined persons were the presence of an established diagnosis of bipolar disorder (F31) based on the ICD-10 criteria, age from 18 to 55 years, the absence of acute somatic pathology, and the presence of a signed informed consent form for participation in this study.
The control group consisted of 7 mentally healthy individuals aged 42 [37; 45] years without signs of acute somatic diseases at the time of the examination.For the control group, the inclusion criteria for this study were age from 18 to 55 years, the absence of somatic and mental diseases, and the presence of written informed consent to participate in this study.
Exclusion criteria were age over 55 years, the presence of somatic pathology in the acute stage, epilepsy, alcoholism, drug addiction, or any other mental disorder not included in this study.

Sample Preparation
Venous blood was taken in the morning on an empty stomach in sample tubes (Becton Dickinson Vacutainer, BD, Vianen, the Netherlands) containing a clot activator.Serum was separated from blood via centrifugation for 20 min at 2000× g in a Digicen 21R centrifuge (Orto Alresa, Madrid, Spain).The serum was then divided into aliquots and stored at −80 • C.
The resulting samples were concentrated via ultrafiltration through 3 kDa Microcon ® centrifuge ultrafilters (Millipore, Molsheim, France) at 14,000 rcf for 15 min at 20 • C. The protein concentration was measured based on the absorbance at 280/260 nm using a Varioskan LUX spectrophotometer (Thermo Scientific, Waltham, MA, USA) located at the core facility Medical Genomics at Tomsk National Research Center.

One-Dimensional Laemmli PAG Electrophoresis
The staining agent (2 mL of glycerol, 2 mL of 0.5 M Tris-HCl pH 6.8, 750 µL of deionized water, 10 µL of 4% SDS, 30 µL of β-mercaptoethanol, bromophenol blue) was added to the samples after purification, as described above, in a ratio of 1: 2. The samples were heated at 95 • C for 5 min and briefly centrifuged to precipitate the condensate.Samples containing 30 µg of total protein were added to the gel and subjected to electrophoresis in a polyacrylamide gel with sodium dodecyl sulfate according to the Laemmli method in 12% PAGE [147].To control the migration of proteins in the gel and further calculate the molecular weights, a set of 14 recombinant highly purified unstained PageRuler proteins 10-200 kDa (Fermentas, Waltham, MA, USA) was used.Electrophoresis was performed in a Protean II xi Cell (Bio-Rad, Hercules, CA, USA) at a voltage of 150-180 V, powered by a PowerPac™ Universal Power Supply (Bio-Rad).The gels were stained with Coomassie brilliant blue G250 (0.1% Coomassie brilliant blue G250, 40% C 2 H 5 OH, 10% CH 3 COOH), and washing was carried out in a solution of 5% CH 3 COOH and 50% C 2 H 5 OH until the stained bands were visualized.
The molecular weights corresponding to the stained protein bands in the gel were automatically calculated using the iBright Imaging Systems FL1500 gel documentation system (Thermo Fisher Scientific) at the base of the Core Facility "Medical Genomics" Tomsk NRMC with the accompanying software relative to molecular weights of protein standards (Broad range, Fermentas, Thermo Fisher Scientific).Statistically significant differences in the distribution of proteins between patients and healthy individuals were determined using Fisher's exact test with Yates' correction.Protein bands that differed statistically significantly between the groups of patients with BD and healthy individuals were subjected to further analysis.These bands were manually cut out of the gel (size ~3 mm) using a scarifier and placed in microtubes.From 8 to 24 samples from different bands from each patient were included in this study.
Next, the Coomassie G250 dye was removed by washing three times with 50 mM NH 4 HCO 3 in 50% acetonitrile for 10-15 min on a shaker until complete discoloration was achieved.Then, the samples were lyophilized for 45 min at 40 • C.

Trypsinolysis
Trypsinolysis was performed using trypsin modified for sequencing (#V511A, Promega, Madison, WI, USA) diluted with the supplied solution (50 mM CH 3 COOH) and then sequentially with 50 mM NH 4 HCO 3 , pH = 8, to a concentration of 0.01-0.025mcg/mL.After that, 20 µL of trypsin solution was added to each sample and incubated at 4 • C for 1 h to swell the gel.Then, the samples were incubated at 37 • C for 18 h to achieve trypsinolysis.After the reaction was complete, 25 mM NH 4 HCO 3 was added to each sample and vortexed, and the supernatants were placed in separate tubes.Next, the peptide mixtures were extracted from the gels with 50% acetonitrile in 5% formic acid; the procedure was repeated three times.The extracts were lyophilized and frozen.

Mass Spectrometry Analysis
Mass spectrometric analysis was performed based on the Advanced Mass Spectrometry Core Facility at the Skolkovo Institute of Science and Technology (Moscow, Russia) according to the Rusanov protocol with minimal changes [148].
The resulting peptide samples were analyzed using an Ultimate 3000 Nano LC HPLC system (Thermo Scientific, Rockwell, IL, USA) coupled to a Q Exactive HF-X hybrid quadrupole-Orbitrap (Thermo Fisher Scientific).
Peptide separation was performed on a C18 column with an inner diameter of 75 µm and a length of 150 mm (Acclaim ® PepMap™ RSLC, Thermo Fisher Scientific, Rockwell, IL, USA).One microliter of sample, equivalent to one microgram of peptides, was loaded directly onto the column and equilibrated isocratically with mobile phase C (2% acetonitrile, 0.1% formic acid).The peptides were then eluted using a linear gradient of 5 to 55% solution B (0.1% formic acid and 80% acetonitrile) at a flow rate of 0.3 µL/min for 75 min, then for another 6 min to reach 99% solution B. Before the next sample was loaded, a 10 min wash with 99% solution B and a 7 min re-equilibration of the column with solution A (0.1% formic acid) was performed.
Then, the eluting peptides were loaded into the mass spectrometer through a capillary at a temperature of 240 • C and an emitter voltage of 2.1 kV.
Samples were analyzed in triplicate in full MS mode followed by a single DDA MS2.Mass spectra were obtained in the mass range 320-1500 m/z with a resolution of 120,000 (MS).Precursor ions were fragmented in the HCD mode (Higher-Energy Collision Dissociation).Tandem mass spectra of the fragments were obtained at a resolution of 15,000 (MS/MS) in the range from 140 m/z to the m/z value determined by the charge state of the precursor, but not more than 2000 m/z.The maximum accumulation time of precursor ions was 50 ms, and that of fragment ions was 100 ms.The target fill value of the automatic gain control (AGC) for precursor and fragment ions was set to 1 × 10 6 and 2 × 10 5 , respectively.
An isolation intensity threshold of 50,000 arbitrary units was defined for precursor selection, and up to 20 of the best precursors were selected for fragmentation at a normalized collision energy (NCE) of 30.The precursor ion isolation width was 2 m/z.Precursors with a charge state of 1+ and more than 5+ were rejected, and all measured precursors were dynamically excluded from the launch of the subsequent MS/MS for 20 s.

Protein Identification and Statistical Analysis
Peak lists obtained from MS/MS spectra were identified using OMSSA version 2.1.9[149] and X!Tandem version X! Tandem Vengeance [150].The search was conducted using SearchGUI version 3.3.12[151].
Protein identification was conducted against a concatenated target/decoy version of the UniProtKB [152] (release 2023_03; 20,423 (target) sequences) database considering the species Homo sapiens).The decoy sequences were created by reversing the target sequences in SearchGUI.The identification settings were as follows: trypsin, with a maximum of 1 missed cleavage 10.0 ppm as MS1 and 0.05 Da as MS2 tolerances; variable modifications: deamidation of N (+0.984016Da), deamidation of Q (+0.984016Da), oxidation of M (+15.994915Da), and propionamide of C (+71.037114 Da).Variable modifications during refinement procedure: acetylation of protein N-term (+42.010565Da), pyrrolidone from E (−18.010565Da), and pyrolidone from Q (−17.026549Da).The rest of the parameters were set to the default values.
Peptides and proteins were inferred from the spectrum identification results using PeptideShaker version 2.2.9 [153].Peptide-Spectrum matches (PSMs), peptides, and proteins were validated at a 1.0% false discovery rate (FDR) estimated using the decoy hit distribution.Proteins were considered reliably identified if at least two of their peptides could be identified.Post-translational modification localizations were scored using the D-score [154] and the phosphoRS score [155] with a threshold of 95.0, as implemented in the compomics-utilities package [156].A phosphoRS score above this threshold was considered a confident localization.
The mass spectrometry data, along with the identification results, were submitted to the Figshare repository.These data can be accessed at https://figshare.com/articles/ journal_contribution/Analysis_of_Mass_Spectrometric_Data_of_Proteins_from_Serum_in_ Patients_with_Bipolar_Disorder_and_and_ealthy_Individuals_Using_the_PeptideShaker_ Software/24033309 (posted date 25 August 2023).
An unlabeled analysis based on the normalized spectral abundance factor (NSAF) was used to detect differences in the relative abundance of detected proteins between study groups [157].
The NSAF depends on the number of identified peptide spectra for each protein and makes it possible to compare the contents of individual proteins in several independent samples, which allows it to be widely used in semi-quantitative proteomics [158][159][160].This index showed good specificity and sensitivity for low-protein samples, which may be useful for determining minor proteins in serum [159,[161][162][163]. Statistically significant differences in the mean NSAF values for each protein in the study group were assessed using a two-tailed unpaired Student's t-test with Bonferroni correction (p = 0.05) in the Statistica 10.0 software package (StatSoft, Hamburg, Germany).

Conclusions
Following the application of the PeptideShaker software package, during the comparison of the serum proteomes of patients with bipolar disorder and healthy individuals, several proteins were identified that are likely to be involved in the pathogenetic processes of BD.When comparing the proteomes of the blood serum of patients with bipolar disorder and healthy individuals, 10 proteins showed significant differences in NSAF values.Of these, four proteins were predominantly present in BD patients with the maximum NSAF value: 14-3-3 protein zeta/delta (NSAF = 0.009346; p = 0.06); ectonucleoside triphosphate diphosphohydrolase 7 (NSAF = 0.002105; p = 0.02); transforming growth factor-beta-induced protein ig-h3 (NSAF = 0.005941; p = 0.02); and B-cell CLL/lymphoma 9 protein (NSAF = 0.018721; p = 0.01).Meanwhile, another six proteins were minimally represented in the serum of BD patients (NSAF ≈ 0.00001).In total, 8 out of 10 proteins were found to be associated with BD for the first time.As these proteins are generally considered to be neurospecific, we believe there is great potential for further studies of the role of these proteins in the pathogenesis of BD and their quantitative contents in a larger number of patients.Conducting such studies will help with the development of new paraclinical criteria for the differential diagnosis of bipolar disorder, among other mental disorders, and the discovery of new targets for BD drug therapy.This will significantly improve the quality and timeliness of BD diagnosis, as well as the choice of the optimal therapy algorithm, significantly improving the quality of life and decreasing the risk of suicidal behavior in these patients.Limitations.This study used a cross-sectional design to recruit patients.Because this was a pilot study, the cohort size was limited.In this study, we tested the PeptideShaker tool to perform a rapid semi-quantitative comparative analysis of serum proteomes involving a large number of proteins.

Figure 1 .
Figure 1.Functional network of protein gene connections that showed statistically significant differences between groups of patients with bipolar disorder and healthy individuals.Line colors indicate the following: lilac-co-expression, pink-physical interaction, orange-predicted, green-genetic interaction.The thickness of the line reflects the strength of the interaction.The tool GeneMANIA (genemania.org) was used.
This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethical Committee of Mental Health Research Institute at Tomsk National Research Medical Center of the Russian Academy of Sciences No. 163 from 12 May 2023 (No. 163/3.2023).

Table 1 .
Statistically significant differences found between the group of patients with bipolar disorder and the control group.