Associations of Genetic Polymorphisms and Neuroimmune Markers With Some Parameters of Frontal Lobe Dysfunction in Schizophrenia

We investigated the associations of DRD3 rs6280, HTR1A rs6295, BDNF rs6265, SCL6A4 rs16965628, and 5HT2A rs7322347 with schizophrenia in a case–control study, and associations of these genetic variants with several clinical features. We also investigated markers of inflammatory response (C-reactive protein, IL-2, IL-6, IL-10), the activity of leukocytic elastase (LE) and α1-proteinase inhibitor (a1-PI), antibodies to S100B and myelin basic protein (MBP) in schizophrenia. Clinical symptoms were assessed on three scales: Positive and Negative Syndrome Scale, The Bush – Francis Catatonia Rating Scale and Frontal Assessment Battery. All SNPs were typed using predesigned TaqMan SNP genotyping assays. The biomarkers related to the immune system were routinely tested using ELISA kits. The association with schizophrenia was found for DRD3 rs6280 (p = 0.05) and HTR2A rs7322347 (p = 0.0013). We found differences between groups by parameters of LE and a1-PI and LE/a1-PI (p < 0.001). And IL-6 was evaluated in the schizophrenia group (p < 0.001). We showed that patients with the TT allele (BDNF rs6265) had more severe impairments in frontal lobe function. a1-PI can serve as a marker for assessing the severity of frontal lobe damage in patients with frontal dementia. We found some biological parameters reflecting the severity of frontal dysfunction in schizophrenia.


INTRODUCTION
Schizophrenia spectrum disorders are severe and debilitating psychiatric disorders hereditary in nature, mostly chronic. They have a heterogeneous neurobiological background affecting brain development. Schizophrenia is a multifactorial disease, where both hereditary, development and environmental risks factors are involved (1).
Due to diverse clinical phenotypes and heterogeneity of the overlapping symptoms and functional deficits in schizophrenia spectrum disorders, their biological background is still not fully understood. A large number of hypotheses regarding the pathogenesis of psychotic disorders is known nowadays. For a long time, it has been known that genetic background is a crucial factor of schizophrenia. The rationale behind genomewide association study (GWAS) is the fact that schizophrenia is mainly associated with common genetic variants (SNPs). Largescale GWAS have identified about 150 risk loci. Most of these common alleles confer only relatively with very small effect (typically odds ratios < 1.2) but cumulatively they have been estimated to explain between a quarter and half of the variance in genetic liability (2). Besides, not all of them are functional (or non-synonymous) SNPs that have a deleterious effect on protein function.
The possible functional effect of the functional missense mutation in the DRD3 gene (Ser9Gly, rs6280) -polymorphism on dopamine release is consistent with previous work demonstrating that the glycine allele yields D3 autoreceptors that have a higher affinity for dopamine and display more robust intracellular signaling (8,9).
5HT1A SNP rs6295 (C-1019G) is located in the promoter region and is known to be a functional polymorphism that regulates HTR1A transcription and region-specific modification of HTR1A expression.
Previous studies identified the association between the rs6295 polymorphism and symptomatic improvement in schizophrenia during treatment with antipsychotics. In particular, the researches showed a significantly greater negative symptoms improvement than G allele (10) carriers and improved attention (11). Meta-analysis of five studies, including the three with lurasidone, demonstrated the association of rs6295 with improvement of positive and negative symptoms in patients with schizophrenia (12).
BDNF rs6265 (Val66Met), a non-synonymous substitution at codon 66: C → T/Valine → Methionine, is the most studied SNP in psychiatric genetics. The functional polymorphism rs6265 produces a profound effect on BDNF cellular biology. This SNP has been linked to BDNF subcellular trafficking.
The effects of the BDNF Val66Met polymorphism on cognition were extensively documented in schizophrenia and healthy individuals. The Met allele was associated with impairments in attentional processing, visuospatial ability and cognitive inflexibility in patients diagnosed with schizophrenia (17). This polymorphism was also associated with depressive symptoms in schizophrenia (18), severe positive symptoms, such as delusions (19) and hallucinations (20), cognitive impairments, and increased risk for schizophrenia (21). 5-HTT is encoded by the SLC6A4 (serotonin transporter) gene. Serotonin transporter protein is primarily responsible for the termination of serotonin neurotransmission. The changes in serotonin transporter protein levels alter the amounts of serotonin receptors, serotonin synthesis and metabolism (22). SNP rs16965628, located on the intron of the SLC6A4 is associated with suicides among patients with schizophrenia. Changes caused by rs1696528 C alleles may lead to alterations in the serotonergic system (23). Moreover, this allele is a candidate polymorphic region for etiology of obsessivecompulsive behavior (24).
Obviously, schizophrenia is a neurodevelopmental disorder with multiple biochemical abnormalities. The earliest hypothesis regarding the development of schizophrenia describes genecontrolled deficiencies in neurotransmitter metabolism (dopamine, serotonin, glutamate, gamma aminobutyric acid, and others) (25).
Later, there appeared a hypothesis that the neuroplasticity and neurotrophic function was a core of disturbances in some brain areas supported by cognitive deficits due to deterioration in central nervous system (CNS) and gradual progressive alterations in gray matter, occurring in neonatal period and adolescence was presupposed as possible mechanism. Neurotrophic factors were proposed as targets to reflect the degree of neurodamage (26).
Recently the role of inflammatory processes causing alterations of brain tissue and CNS functioning has also been considered as explanation to schizophrenia development, this is why immune status is an important indicator of course of schizophrenia. For psychotic disorders complex immunebrain interactions that affect neural development, survival, and functioning could play a substantial role (1,27). High concentrations of the circulating proinflammatory cytokines and other immune markers are reported to be associated with increased risk of subsequent psychosis (27).
The given evidence suggests that schizophrenia is accompanied by changes in pro-inflammatory cytokines and other inflammatory markers in the brain (neuroinflammation) and the activation of the inflammatory responses in blood (28).
Changes in cytokine profiles in schizophrenic patients imply immune and neurotransmission dysregulation (29). IL-2, IL-6, IL-10 are known markers of immune status in schizophrenia; and the alterations of their levels are connected with the severity of clinical symptoms and potential clinical implications. IL-2, IL-6 are increased and IL-10 is decreased in the first-episode psychosis vs. controls and in acutely ill patients with chronic schizophrenia vs. controls (27). Moreover, IL-6 concentrations are positively associated with both positive and negative symptoms, as assessed on the Positive and Negative Syndrome Scale (PANSS) (30).
There is also an evidence for associations between higher blood IL-6 and C-reactive protein (CRP) levels and cognition in patients with chronic and first-episode schizophrenia, including executive function, verbal and working memory, learning, and psychomotor speed, although these associations are not always replicable (31,32).
Immunological markers of leukocyte (neutrophil) elastase (LE), α1-proteinase inhibitor (α1-PI) activities, and autoantibodies' levels against S-100-beta and myelin basic protein (MBP) were chosen because of their association with patients' clinical state assessed by PANSS. The significant correlations were found between the severity of patients' psychosis and the activity of LE in the blood, as well as anti-S100B antibodies (33,34).
LE is a proteolytic enzyme of neutrophils that hydrolyzes proteins of blood plasma, including coagulation factors, fibrinolysis, kallikreinkinin system and complement. LE can be a destructive factor in relation to the blood-brain barrier permeability (35). The biological role of LE and its contribution to brain pathology was reviewed by Shimakura et al. (36). a1-proteinase inhibitor is the acute phase protein that regulates the proteolytic activity of LE and other proteinases and also provides conditions for limiting the source of inflammation and destruction (37). The ratio of the activity of LE and its inhibitor can be predictors of tissue damage in the brain determining the outcome of the inflammatory response (38). An increase in the level of antibodies to S100B and MBP characterize more severe course of the schizophrenic process (39).
It is known that cognitive dysfunctions and abnormalities in the underlying brain processes are considered, respectively, as core clinical symptoms and biological features in the pathophysiology of psychiatric disorders.
Studies during the past 20 years have found complex interactions between the immune system, systemic inflammation and the brain, leading to changes in mood, cognition, and behavior (1).
Cognitive impairments in schizophrenia are associated with neuroinflammatory markers such as CRP, S100B and neuronspecific enolase (40,41), and with lower concentrations of BDNF (40).
The aim of this study was to research the associations between symptomatology and SNP of genes related to the systems of neurotransmission (DRD3 rs6280, HTR1A rs6295, 5HT2A rs7322347) and neurotrophic factor (BDNF rs6265) in schizophrenia patients. Also we investigated associations between some immunological patterns and markers of inflammatory response (CRP, IL-2, IL-6, IL-10, the activity of LE and a1-PI, antibodies to S100B and MBP) with clinical symptoms in schizophrenia.

MATERIALS AND METHODS
Patients (SZ) (n = 655, male-345, female-310) were recruited from Mental-health Clinic No. 1 named after N.A. Alekseev of Moscow Healthcare Department, Moscow (Russia). The diagnoses included in the study were (F20 -627 patients) Schizophrenia, (F22 -3 patients) Persistent delusional disorders, (F23 -8 patients) Acute and transient psychotic disorders, (F25 -17 patients) Schizoaffective disorders. The average age of patients was 33.6 ± 9.8 in women and 30.1 ± 7.9 in men. The average age of onset was 24.6 ± 8.1 in women and 22.3 ± 5.7 in men. The average total PANSS was determined to be 105 ± 24 for women and 109 ± 22 for men.
Healthy controls (HC) (n = 768) were recruited from the same geographical area as SZ group. All participants gave written informed consent for participation in the project.
Exclusion criteria for patients were age < 18 years, serious medical or surgical illness, previous episode of psychosis due to substance abuse, and psychotic symptomatology within a clearly diagnosed affective or borderline personality disorder. For healthy controls, exclusion criteria were age < 18 years, current or past psychiatric disorder, family history of any psychiatric disorder, head trauma, neurological illness, serious medical or surgical illness, or substance abuse. HC and SZ had no addictions including <80 mg in ethanol equivalent per week and <10 cigarettes per day.
SZ status was determined based on the International Classification of Disease (ICD-10). Patients were evaluated using a structured interview for the Positive and Negative Syndrome Scale (PANSS), including the PANSS positive, PANSS negative, and PANSS general psychopathology subscales, Frontal Assessment Battery (FAB), The Bush -Francis Catatonia Rating Scale (BFCRS). The Russian version of PANSS, with acceptable validity and reliability, was used to evaluate the severity of psychosis symptoms of schizophrenia patients in this study.
The protocol of this study was approved by the Interdisciplinary Ethics Committee, Moscow (22/07/2017). All data is available from the authors on request.
The antipsychotics used were those with the highest clinical response in each individual patient, medication status was described in previous article (42). We also added treatment information to Table 1 [Chlorpromazine equivalents are taken from the work Gardner et al. (43)].
We took samples from patients in the hospital who needed urgent drug treatment. They mostly receive such treatment as early as the first aid stage before admission to the hospital. Therefore, it was difficult to completely exclude antipsychotic medication in our study. For our research we always take blood within 3-4 days of the patient's admission to hospital. A total of 10 ml fasting blood was collected from each subject by venipuncture in the morning before breakfast. Five milliliter was introduced into a free-anticoagulant vacuum tube and immediately centrifuged at 2,000 × g for 10 min. Five milliliter was placed into an ethylenediaminetetraacetic acid (EDTA) vacuum tube for DNA analysis. Serum and total blood were stored at −80 • C for further analysis.

Genotyping
Genomic DNA samples were obtained from peripheral blood lymphocytes using an automatic DNA extraction (QIAGEN QIAcube) system, according to the manufacturer's recommendations. DNA concentration and sample quality were assessed spectrophotometrically (NanoVue, GE Healthcare). The obtained DNA samples were normalized in TE buffer to a final concentration of 4 ng/µl in the 384-well plate format. All SNPs were typed using predesigned TaqMan SNP genotyping assays (Applied Biosystems, Thermo Fisher, USA). Assays ID: C__32472275_10, C__11904666_10, C__11592758_10, C____949770_20, C__29706403_10. PCR and allelic discrimination assays were run using a QuantStudio 5 Real-Time PCR system (Applied Biosystems, Thermo Fisher, USA).

Immunological status
Immunological parameters were assessed in serum from blood collected from 363 patients and 150 healthy controls. Samples were stored at −20 • C for 1-2 weeks before assays. LE activity (a measure of neutrophil degranulation activity) was assayed by an enzymatic method using the chromogenic substrate specific N-tert-butoxy-carbonyl-alanine-β-nitrophenyl ester (BOC-Ala-ONp) (Sigma and was expressed in nmol/min × ml. The LE-α1-PI complex was dissociated by adding acetonitrile to each serum sample (44). α1-PI functional activity was assayed spectrophotometrically and was determined based on its ability to inhibit esterase activity of trypsin, using Nα-Benzoyl-L-arginine ethyl ester (BAEE) as a substrate. The measurements were performed using an Ultrospec 1100 spectrophotometer (Amercham) and expressed in IU/ml (45). Antibodies to S-100-beta protein and myelin basic protein were assayed using an immunoenzyme assay (46). S-100B protein from bovine brain (product No. S6677, Sigma-Aldrich, USA) and myelin basic protein bovine (product No. M1891-1MG) were used for ELISA. The proteins were dissolved in deionized water at 1 mg/ml. Antibody titers were expressed in optical density units (OD) and each sample was assayed in duplicate. C-reactive protein (CRP) (product No. A-9002 CRP-EIA-BEST highly sensitive) was assayed using a solidphase immunoenzyme assay using a reagent kit (Vector-Best, Novosibirsk, Russia) and was expressed in IU/l. IL-2, IL-6, and IL-10 were determined using the diagnostic kits for

0.065
The were expressed in pg/ml. These kits are not evaluated the reference proteins.

Statistical Analysis
All parameters were checked for normal distribution using a Shapiro-Wilk test. The association between genotype or immunological markers and clinical characteristics was evaluated using ANOVA Fisher test. To compare the genetic variations in the case-control study, the chi-square test was used. Descriptive statistics of neuroinflammatory biomarkers in the group of patients and healthy control was presented median with low and upper quartile (Me, LQ, UQ) because of non-normal distribution. For comparison with immunological marker controls, the Kruskal-Wallis test was used. Dispersion analysis of association of immunological markers with parameters of frontal dysfunction on the FAB, in the group of patients with signs of frontal dementia was evaluated using ANOVA Fisher test. The classification of objects (the values of immunological markers) based on the distances between them was carried out by two-step cluster analysis with first standardization of variables. Bonferroni and FDR-corrections were applied to correct the pvalues to control inflation of the type I error rate. The differences were considered significant at p < 0.05. Statistical analyses were performed with the R program ( © 2016 The R Foundation) and IBM SPSS Statistics 26. Table 1 depicts the demographic profile of the patients and healthy subject at baseline. No differences were found between patients and controls according to age, gender, and ethnicity. Comparison of allele distribution frequencies of the investigated SNPs revealed significant differences for DRD3 rs6280 (p = 0.05) and HTR2A rs7322347 (p = 0.0013) genes ( Table 2). Comparison of allele frequencies in the case-control study: for SCL6A4 rs16965628 chi-squared = 0.3, p-value = 0.86; for HTR1A rs6295 chi -squared = 3.5, p-value = 0.18; for BDNF rs6265 chi -squared = 0.9, p-value = 0.64; for DRD3 rs6280 chi -squared = 5.92, p-value = 0.05; for HTR2A rs7322347 chi -squared = 13.3, p-value = 0.0013.

RESULTS
During the analysis of SNPs relationship with PANSS we did not find significant associations. Nevertheless, we detected the associations with parameters of BFCRS: Waxy flexibility (F = 4.83 p = 0.0085), Negativism (F = 4.28 p = 0.01), Stupor (F = 5.9 p = 0.003), Catalepsy (F = 5.63 p = 0.004). Patients with genotype T/T (Met66Met polymorphism) demonstrated higher scores, which reflect more severe symptoms. Additionally, we observed that these symptoms were not related to the age of the patient.
We also found associations between BDNF rs6265 and FAB. In particular, patients with genotype T/T (Met66Met) showed lower FAB scores (p = 0.011 for female and 0.0018 for male).  Go (F = 2.14 p = 0.12); 6. Prehension Behavior (did not pass Bonferroni correction, F = 4.11 p = 0.017). We did not observe any differences between SZ and HC in terms of parameters IL-10, CRP, antibodies to S100B and MBP ( Table 3). We detected the statistical differences between groups by parameters of the activity of LE and a1-PI, and LE/a1-PI (p < 0.001). IL-6 was evaluated in both male and female in comparison to HC groups (p < 0.001) ( Table 3).
The correlation between the activity of LE and a1-PI, a-S100B, and a-MBP and the total score of the PANSS scale (including the subscale of positive, negative symptoms and subscale of general psychopathology), NSA, BFCRS did not reveal statistically significant differences.
For the activity of LE this dependence was demonstrated with F = 6.018 p = 0.0148, but this result did not pass for multiple comparisons.
The considerable variability for each of the analyzed markers in the total group of patients served as the basis for cluster analysis. Three immunological subgroups differing in the activity of LE, LE/a1-PI, and the level of autoimmune reactions in comparison to the HC group were revealed. A significant increase in the activity of LE in all groups compared to controls was shown (p < 0.001). The differences in profiles of immunological markers in subgroups of schizophrenia patients reflect the different types of the inflammatory response to pathological process ( Table 4).
Analysis of associations of markers in immunological subgroups with the severity of clinical symptoms revealed no differences in the total FAB score and in the total BFCRS score.
We showed that a subgroup of patients with signs of frontal dementia (11 or lower points on the FAB scale) (38 male, 39 female) significantly differs by immunological parameters from others. Moreover, these differences did not depend on the age of the patients. For them the dispersion analysis of the connection of these inflammation markers with the total FAB score and with each parameter of this scale (see Supplementary Materials) was carried out.
No statistically significant associations between genotype and immune status parameters were found (data on request).

DISCUSSION
In the first part of our study, we demonstrated some associations between genetic polymorphism with cognitive and motor disturbances in patients with schizophrenia spectrum disorder and confirmed the involvement of certain genetic factors in symptomatology of the disease. In the second part of the study we verified the involvement of neuroimmunological biomarkers in schizophrenia spectrum disorders.
We have shown the associations of HTR2A rs7322347 and DRD3 rs6280 with schizophrenia spectrum disorder in the case-control study. This data confirms our earlier findings (16). Despite the fact that the associations of BDNF rs6265,   HTR1A rs6295, SCL6A4 rs16965628 with schizophrenia in case-control study was not found, we discovered interesting correlations between BDNF rs6265 and some features of the disease. Met66Met polymorphism was associated with high scores of catalepsy, waxy flexibility, stupor and negativismserious symptoms of catatonia measured by BFRS. These signs of motor abnormalities can indicate the severity of disease (47). Besides, it is known that the catatonia symptoms can reflect the extreme severity of affective deficit (48). According to other studies patients with schizophrenia accompanied by catatonic symptoms are found to have lower serum BDNF levels (49,50), which proves the link between BDNF and catatonia symptoms (51). According to findings of other researcher, motor functions significantly correlate with cognitive function in patients with schizophrenia (52). In addition to the association of Met66Met with motor disturbance, we also found a correlation with frontal lobe dysfunction. Patients with Met66Met had low scores of the FAB. Other studies also showed causality between Met allele of BDNF and cognitive disability (53)(54)(55)(56). Kim et al. (57) demonstrated the association between BDNF rs6265 and higher scores in PANSS scale and lower in cognitive tests in female patients with schizophrenia, but no performance difference was observed in male patients. In our study frontal lobe cognitive impairment was associated with Met allele in both male and female, regardless of their age.
In the second part of our research we have analyzed the serum CRP, IL-2, IL-6, IL-10 levels, the activity of LE and a1-PI, antibodies to S100B and MBP and their associations with clinical symptoms. We showed significant increase of IL-6 level in patients with schizophrenia spectrum disorder compared to healthy controls. Moreover, we demonstrated that the increase of serum activity of LE and a1-PI may be related to the diagnosis of schizophrenia, which confirms the data we obtained earlier (33) so as the other authors (58,59). It has previously been demonstrated that patients with schizophrenia have lower serum BDNF, MBP, and GFAP levels, but higher serum IL-6 and S100B (58). The Ota 2015 study showed an increase of MBP gene expression in the peripheral blood of antipsychotic-naïve patients with first-episode psychosis (59).
The ratio of the activity of LE and its inhibitor characterizes the state of the leukocyte-inhibitory system and is considered as an unfavorable factor in the maintenance of tissue damage (38).
The differences of immunological profiles are associated with the different state of the leukocyte-inhibitory system of inflammation in schizophrenia.
As a result of cluster analysis, three immunological groups (clusters) associated with different state of leukocyte-inhibitory system in the response to the pathological process in brain were revealed.
The results are consistent with previously obtained data and indicate the involvement of the immune system (inflammation) in pathogenesis of schizophrenia. The state of leukocyteinhibitory system of inflammation assessed by leukocyteinhibitory index determines the course and the outcome of the inflammatory reaction. Pro-inflammatory potential associated with the increase in the activity of both LE and α1-PI compared to control group (p < 0.001 and p < 0.001) in patients of cluster 1 indicates the activation of systemic inflammatory reactions in the response to impaired homeostasis.
Two and 3 subgroups are characterized by an insufficient functional activity of neutrophils (by the activity of LE) and a significant increase of a1-PI activity (p < 0.05, p < 0.001, respectively). These profiles reflect a different degree of imbalance in the inflammatory response to impaired brain homeostasis.
Many studies indicate the assosiation between neural inflammation markers and cognitive functions (31,32). In our study, cognitive functions were evaluated using the FAB, which indicates the degree of frontal dysfunction.
Bao demonstrated that serum levels of a1-PI in patients with chronic disorder of consciousness may have prognostic potential to disease severity (60). The prefrontal cortex has extensive mutual relations with wake-promoting centers in the brainstem and diencephalon, and therefore is in a unique position to modulate the level of consciousness. Thus, cholinergic mechanisms in the prefrontal cortex can regulate the level of consciousness (61). It has been shown that loss of consciousness after some brain injuries is associated with the reduced functional connectivity of the default mode network (DMN), fronto-parietal network, and thalamo-cortical network (62). Cognitive functions have been shown to be related to consciousness, suggesting that the neural correlates of memory formation and retrieval coincide with the neural correlates of conscious perception (63).
Presumably, the imbalance of leukocyte-inhibitory system associated with insufficient activity of LE (p < 0.05) and high activity of its inhibitor α1-PI (p < 0.001) in patients of cluster 3 may be due to functional depletion of neutrophils. This immune profile is associated with the severity of clinical symptoms in schizophrenia patients with frontal dysfunctions. Obtained results indicate that the imbalance of leukocyte-inhibitory system can be considered as a marker of the severe disease course and as an unfavorable prognostic factor. A similar immunological feature was previously identified in patients with catatonic disorders in schizophrenia (64).
Thus, the fact that the frontal lobe dysfunction altered consciousness, and data on a significant increase of this indicator activity accompanied by an insufficient increase of LE activity, as well as the correlation of a1-PI in blood serum with the severity of the consciousness disorder suggest that functional activity of this marker may serve as a potential prognostic factor for frontal function disorders in schizophrenia in our patients. According to the literature, patients with scores 11 and below have signs of frontal dementia. The group of patients with frontal dementia showed a decrease in the activity of a1-PI while increasing the score for each of the scale parameters. However, no such associations were found during the analysis of all patients. It was shown for the first time that activity of a1-PI can serve as a potential biomarker of cognitive deficit severity in the group of patients who demonstrates the signs of frontal dementia, but not in patients with normal frontal lobes function.
Despite the absence of differences in MBP serum concentration with healthy controls, in the group of patients with frontal dementia the inverse dependence of serum concentration on parameters such as 3. Motor Series ("Luria's Test"), 4. Conflicting Instructions, 5. Go-No Go (inhibitory control), 6. Prehension Behavior and FAB total score were found. We demonstrated that the serum activity of LE and a1-PI, and antibodies to the protein S100B and MBP may be related to the diagnosis of schizophrenia, which confirms the data obtained earlier by us (33).
The prefrontal cortex has extensive mutual relations with wake-promoting centers in the brainstem and diencephalon, and therefore is in a unique position to modulate the level of consciousness. Thus, cholinergic mechanisms in the prefrontal cortex can regulate the level of consciousness (61). It has been shown that loss of consciousness after some brain injuries is associated with the reduced functional connectivity of the default mode network (DMN), fronto-parietal network, and thalamocortical network (62).

CONCLUSIONS
The feature of this study is the combined investigation of genotype and immune markers in the patients with schizophrenia spectrum disorder. As statistically significant associations between Met of BDNF gene and a1-PI and the degree of frontal dysfunction were demonstrated, these parameters can serve as markers for assessing the severity of frontal lobe damage in schizophrenia patients. Observed association between Met of BDNF gene and some severe signs of the motor dysfunction (catatonia signs) requires further study and verification.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are publicly available. This data can be found here: https://www.ebi.ac.uk/ ena/browser/view/PRJEB42835.

ETHICS STATEMENT
The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Interdisciplinary Ethics Committee, Moscow (22/07/2017).
The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
AM, YZ, and SZ: writing-original draft preparation. OA and VU: writing-review and editing. KP and OP: methodology-genotyping. AVM, AR, and GK: methodology-patients. SZ, IO, and ZS: methodology-immune analysis. YZ and SZ: data analysis. GK and VC: supervision. AM and GK: project administration. GK: funding acquisition. All authors contributed to the article and approved the submitted version.