Journal of Personalized Medicine Article Cannabis Use Induces Distinctive Proteomic Alterations in Olfactory Neuroepithelial Cells of Schizophrenia Patients Marta Barrera-Conde 1,2 , Karina Ausin 3 , Mercedes Lachén-Montes 3 , Joaquín Fernández-Irigoyen 3 , Liliana Galindo 4,5 , Aida Cuenca-Royo 1 , Cristina Fernández-Avilés 2 , Víctor Pérez 5,6 , Rafael de la Torre 1,2,7 , Enrique Santamaría 3 and Patricia Robledo 1,2, * 1 Integrative Pharmacology and Systems Neuroscience, Neuroscience Research Program, IMIM-Hospital del Mar Research Institute, 08003 Barcelona, Spain; mbarrera1@imim.es (M.B.-C.); acuenca@imim.es (A.C.-R.); RTorre@imim.es (R.d.l.T.) 2 Department of Experimental and Health Sciences, University Pompeu Fabra, 08003 Barcelona, Spain; cristina.fernandeza@upf.edu 3 Clinical Neuroproteomics Unit, Proteomics Platform, Navarrabiomed, Complejo Hospitalario de Navarra (CHN), Universidad Pública de Navarra (UPNA), IdisNA, Proteored-ISCIII, 31006 Pamplona, Spain; karina.ausin.perez@navarra.es (K.A.); mercedes.lachen.montes@navarra.es (M.L.-M.); jokfer@gmail.com (J.F.-I.); enrique.santamaria.martinez@navarra.es (E.S.) 4 Department of Psychiatry, University of Cambridge, Cambridgeshire and Peterborough NHS Foundation Trust, Cambridge CB2 1TN, UK; lg532@cam.ac.uk  5 Neuropsychiatry and Addictions Institute (INAD) of Parc de Salut Mar, 08003 Barcelona and CIBER de Salud  Mental, Spain; 61155@parcdesalutmar.cat 6 Department of Psychiatry and Legal Medicine, Autonomous University of Barcelona, 08193 Barcelona, Spain Citation: Barrera-Conde, M.; 7 Centro de Investigación Biomédica en Red de Fisiopatología de la Obesidad y Nutrición, Instituto de Salud Ausin, K.; Lachén-Montes, M.; Carlos III, 28029 Madrid, Spain Fernández-Irigoyen, J.; Galindo, L.; * Correspondence: probledo@imim.es; Tel.: +34-93-316-0455 Cuenca-Royo, A.; Fernández-Avilés, C.; Pérez, V.; Abstract: A close epidemiological link has been reported between cannabis use and schizophrenia de la Torre, R.; Santamaría, E.; et al. (SCZ). However, biochemical markers in living humans related to the impact of cannabis in this Cannabis Use Induces Distinctive disease are still missing. Olfactory neuroepithelium (ON) cells express neural features and offer Proteomic Alterations in Olfactory Neuroepithelial Cells of a unique advantage to study biomarkers of psychiatric diseases. The aim of our study was to Schizophrenia Patients. J. Pers. Med. find exclusively deregulated proteins in ON cells of SCZ patients with and without a history of 2021, 11, 160. https://doi.org/ cannabis use. Thus, we compared the proteomic profiles of SCZ non-cannabis users (SCZ/nc) and 10.3390/jpm11030160 SCZ cannabis users (SCZ/c) with control subjects non-cannabis users (C/nc) and control cannabis users (C/c). The results revealed that the main cascades affected in SCZ/nc were cell cycle, DNA Academic Editor: Tomiki Sumiyoshi replication, signal transduction and protein localization. Conversely, cannabis use in SCZ patients induced specific alterations in metabolism of RNA and metabolism of proteins. The levels of targeted Received: 14 January 2021 proteins in each population were then correlated with cognitive performance and clinical scores. In Accepted: 20 February 2021 SCZ/c, the expression levels of 2 proteins involved in the metabolism of RNA (MTREX and ZNF326) Published: 25 February 2021 correlated with several cognitive markers and clinical signs. Moreover, use duration of cannabis negatively correlated with ZNF326 expression. These findings indicate that RNA-related proteins Publisher’s Note: MDPI stays neutral might be relevant to understand the influence of cannabis use on SCZ. with regard to jurisdictional claims in published maps and institutional affil- Keywords: cannabis; schizophrenia; proteomics; olfactory neuroepithelium; metabolism; RNA; iations. ZNF326; MTREX Copyright: © 2021 by the authors. 1. Introduction Licensee MDPI, Basel, Switzerland. This article is an open access article The complexity of schizophrenia (SCZ), involving intricate interactions between envi- distributed under the terms and ronmental factors and genetics, hinders the identification of the molecular mechanisms conditions of the Creative Commons underlying its development [1]. State of the art techniques are being used to find biomark- Attribution (CC BY) license (https:// ers and ease clinical practice. In this sense, proteomic tools, allowing the study of a large creativecommons.org/licenses/by/ number of proteins at a time, can provide an integrated picture of the biological dysfunc- 4.0/). tions associated with SCZ [2]. Serum, plasma and cerebrospinal fluid studies reveal an J. Pers. Med. 2021, 11, 160. https://doi.org/10.3390/jpm11030160 https://www.mdpi.com/journal/jpm J. Pers. Med. 2021, 11, 160 2 of 15 effect on proteins related to the innate immune system and a highly reported downreg- ulation of Apolipoprotein A-I1 (APOA1) [3–5], which might mediate the impact of SCZ on the peripheral immune response [6,7]. On the other hand, postmortem brain studies have shown deregulations in metabolic pathways, calcium signalling and cytoskeleton assembly in several brain areas in SCZ [8–10]. However, other neuropsychiatric condi- tions present similar affectations [11]. In fact, the two most referenced biomarkers of SCZ, namely Fructose-bisphosphate aldolase C (ALDOC) and Glial fibrillary acidic protein (GFAP), are also altered in major depression and bipolar disorder [12]. Thus, postmortem brains from individuals diagnosed with different mental conditions appear to have great similarities that hamper the identification of exclusive biomarkers for SCZ. Additionally, the impossibility of correlating these biochemical findings with cognitive dysfunctions and environmental cues specific for each condition highlights the necessity to use other substrates. In this context, the olfactory neuroepithelium (ON) has emerged as a useful tool. This specialized epithelial tissue can be easily obtained from living subjects and it contains multipotent progenitors that express neural markers [13–16]. Therefore, the main advantage of ON cells is that ongoing environmental factors and molecular alterations can be evaluated during the progression of the disease. Years of epidemiological research in psychotic disorders suggest that cannabis use is a specific environmental factor for SCZ. The incidence of cannabis use is two to four times higher among SCZ patients, and its consumption doubles the risk of psychosis in a dose- and time-dependent manner [17,18]. Moreover, an endocannabinoid system imbalance, including alterations in the expression of the cannabinoid receptor 1 (CB1 ) in postmortem brains and a diminished production of the enzymes involved in the synthesis of endocannabinoids, has been observed in postmortem and peripheral samples of SCZ patients [19,20]. Additionally, cannabis and SCZ have a shared effect on some cellular func- tions. Indeed, protein metabolism via AKT/mTOR, cytoskeleton organization, and calcium signalling are affected in human samples from SCZ patients [2], and after cannabinoid ex- posure in cultured cells [21,22] and mice models [23,24]. Moreover, plasmatic samples from heavy cannabis users and SCZ patients share an effect on oxidative stress response [7,25]. However, whether the coexistence of both factors has equal or distinct consequences is less clear. In fact, in a synchrotron-based infrared spectroscopy study, significant differences in the protein spectra were reported in ON cells from SCZ cannabis users compared to non-users [26]. Additionally, we recently demonstrated that the functional signature of the CB1 and serotonin 2A receptor (CB1 -5HT2A R) heteromer in ON cells is differentially regulated in SCZ patients depending on cannabis use [27]. Therefore, the aim of this study was to find exclusively deregulated proteins in ON cells of SCZ patients with and without a history of chronic cannabis use. We applied mass spectrometry-based quantitative proteomics followed by a functional analysis to characterize the proteomic changes specifically induced by cannabis use in ON cells from SCZ patients, as compared to controls that use cannabis or not. The primary outcome of the study, i.e., the expression levels of the target proteins, were correlated with cannabis consumption, and with secondary cognitive and clinical outcomes associated with SCZ. 2. Materials and Methods 2.1. Study Design SCZ patients without a history of chronic cannabis use (SCZ/nc, n = 5), SCZ pa- tients with a reported history of chronic cannabis use (SCZ/c, n = 5), chronic cannabis users without any psychiatric diagnosis (C/c, n = 5), and control subjects non-cannabis users (C/nc, n = 5) between 18 and 45 years old (both males and females), were recruited to conduct a cross-sectional study. Every subject gave written informed consent after a complete description of the study and the procedures involved. Cannabis users had to consume more than 5 cannabis cigarettes per week during at least 6 months. On the day of testing, subjects were told to refrain from cannabis use for at least 12 h before testing to avoid the potential confounding of acute cannabis intoxication in both neuropsychological J. Pers. Med. 2021, 11, 160 3 of 15 and proteomic studies. The exclusion criteria were: (i) meeting criteria for any severe mental disorder according to the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5); (ii) history of severe mental illness among first degree relatives; (iii) history of severe congenital, medical or neurological illness; (iv) show medical condi- tions with nasal repercussions (rhinitis or bleeding); and (v) use of other drugs of abuse. SCZ-diagnosed participants were receiving antipsychotic medication at the time of the study. All subjects underwent a physical, psychiatric and neuropsychological evaluation, as previously described [16]. This study was approved by the local Institutional ethics committee (CEIC-PSMAR). 2.2. Clinical and Neuropsychological Evaluation Neuropsychological assessments were carried out in all subjects. Executive functions were tested with the semantic verbal fluency test. Attention performance was evaluated using the direct spatial span (SSP) and digit direct series. Working memory was addressed using the inverse spatial span and digit inverse series using the Cambridge neuropsycho- logical test automated battery (CANTAB 2017), and the digit span of the Wechsler Adult Intelligence Scale (WAIS-III). Emotional recognition was evaluated with the CANTAB test. Clinical outcomes included the global assessment of functioning (GAF), and the neurological soft signs (NSS) scores. 2.3. Quantification of Cannabis Metabolites To estimate the amount of cannabis consumed by cannabis users, the plasma con- centrations of its main non-psychoactive metabolite, namely, 11-nor-9-carboxy-∆-9-THC (THC-COOH), were calculated. Briefly, 1 mL of plasma was transferred into a glass tube and spiked with d3-∆-9-THC as the internal standard. A protein precipitation with 2 mL of 0.1% formic acid in acetonitrile was performed prior to a solid phase extraction with Oasis Prime HLB 3 cm3 , 60 mg column (Waters Co., Milford, MA, USA). The supernatant was diluted with MilliQ water and loaded. Subsequently, 2 mL of 25% of methanol was added twice. Elution was carried out twice with 2 mL of 90:10 acetonitrile:methanol (ACN:MeOH). The organic phase was evaporated to dryness under a nitrogen stream at < 39 ◦ C and <15 psi pressure. Analytes were reconstituted in 50 µL of 90:10 ACN:MeOH and 50 µL of MilliQ water. Quantification of THC-COOH in plasma was performed using an Agilent 1200 series HPLC system (Agilent Technologies) coupled to a 6410 Triple Quadrupole LC-MS (Agilent Technologies) mass spectrometer with an electrospray interface. 2.4. Nasal Exfoliation and Cell Culture ON samples were obtained by nasal brushing. Samples from the middle and upper turbinates were maintained in 250 µL of cold Dulbecco’s Modified Eagle Medium/Ham F-12 (DMEM/F12) enriched with 10% FBS, 2% glutamine and 1% streptomycin–penicillin (GibcoBRL), as previously described [27]. At 80% confluence, cells were expanded using 0.25% trypsin (GibcoBRL) and replated in 75 cm2 flasks. Then, cells were expanded until passage four was reached and stored in liquid nitrogen with 20% FBS and 10% dimethyl sulfoxide (Sigma-Aldrich, Madrid, Spain). The aliquots from all subjects (n = 20) were unfrozen at the same time and harvested in 75 cm2 flasks under standard conditions using enriched DMEM/F12 medium. 2.5. Quantitative Proteomics in ON Cells Physical protein extraction was performed using sterile scrappers and PBS 1X. After centrifugation, the pellet was kept frozen at −80 ◦ C until homogenization using lysis buffer (7 M urea, 2 M thiourea, and 50 mM DTT). The homogenates were then subjected to in-solution digestion, peptide purification and reconstitution prior to mass spectro- metric analysis (MS/MS). Data acquisition was performed as previously described [28] using ProteinPilot v5.0 (Sciex) as a search engine, ParagonTM Algorithm (v.4.0.0.0) [29] for database searching and a non-lineal fitting method to calculate the false discovery J. Pers. Med. 2021, 11, 160 4 of 15 rate (FDR) (1% Global FDR or better). Data analysis was performed by sequential win- dow acquisition of all theoretical mass spectra–mass spectrometry (SWATH-MS), and the TripleTOF 5600+ mass-spectrometer was configured as previously described [30]. The library generation-associated ProteinPilot group file was loaded into PeakView® 2.1 (Sciex), and peaks from SWATH runs were extracted with a peptide confidence threshold of 99% confidence (Unused Score ≥1.3) and FDR lower than 1%. ProteinPilot was used to ex- tract the MS/MS spectra of the assigned peptides, and only proteins quantified with at least two unique peptides were considered. For more detailed information about the SWATH-MS library generation, please refer to the extended experimental procedures in the Supplementary Information. 2.6. Pathway Analysis, Statistics and Bioinformatics The peer-reviewed pathway database REACTOME [31] was used to functionally characterize the proteomic alterations detected in every comparison. This database hi- erarchically classifies cellular functions in 27 cascades. We established a 0.01 p-value threshold, and the functional proteome results obtained were ordered based on their FDR. All calculations were performed using SPSS (SPSS Inc., Chicago, IL, USA). Normality and homoscedasticity for continuous variables were tested using Shapiro–Wilk W and Levene tests. Demographic categorical variables were evaluated using Chi-squared tests and continuous normally distributed variables were evaluated using one-way ANOVAs followed by the least significance difference (LSD) post-hoc test. Neuropsychological calcu- lations were corrected for tobacco use duration. Not normally distributed variables were compared using Kruskal–Wallis and Dunn’s post-hoc test. Vulcano plot representations and Spearman correlation plots were depicted using GraphPad prism 8.0 Software (La Jolla California, CA, USA). Venn diagrams for proteomic comparisons were designed using BioVenn platform [32]. 3. Results 3.1. Demographics and Neuropsychological Outcomes The demographic data are shown in Table 1. The groups did not differ in age, sex or dosage of tobacco use. However, SCZ/c showed a significantly higher duration of tobacco use (years) when compared to all other groups (p < 0.05). The analysis of cannabis use patterns between C/c and SCZ/c revealed no significant differences. The neuropsychologi- cal evaluation (Table 2) revealed significantly lower verbal fluency scores for SCZ/nc and SCZ/c as compared to C/nc (p < 0.05), but not when compared to C/c. GAF scores were significantly lower in SCZ/nc and SCZ/c when compared to C/nc and C/c (p < 0.001). Finally, significantly more NSS were present in SCZ patients regardless of cannabis use when compared to C/nc (p < 0.01) and C/c (p < 0.05). 3.2. Plasmatic Concentrations of THC-COOH The plasmatic concentrations of THC-COOH were not significantly different between the groups (C/c: 34.92 ± 17.55 ng/mL; SCZ/c: 11.22 ± 15.91 ng/mL). 3.3. Proteomic Analyses Our proteomic workflow was fundamentally based on a triangular approach [33], whereby a relatively small number of cases and controls are analyzed by hypothesis-free discovery proteomics in great depth, leading to the quantification of thousands of proteins. Using this workflow, we observed specific changes in protein abundance in C/nc vs. SCZ patients and in SCZ/nc vs. SCZ/c, underpinning the importance of analyzing independent factors (disease/cannabis), instead of merely contrasting SCZ versus control cases. Firstly, C/c, SCZ/nc and SCZ/c were compared to C/nc subtracting exclusively deregulated proteins for each population. Then, each group of SCZ patients was compared to C/c to subtract the effect of cannabis use on control subjects. In the end, each group of SCZ (SCZ/nc and SCZ/c) was compared to obtain proteomic biomarkers that were only affected J. Pers. Med. 2021, 11, 160 5 of 15 by SCZ or by SCZ plus cannabis use. For the summary followed in the proteomic data analysis, see Supplementary Figure S1. Table 1. Demographic characteristics of the different groups included in the study. Schizophrenia Control Subjects Cannabis Users Schizophrenia Non-Cannabis (C/nc) (C/c) Cannabis (SCZ/c) (SCZ/nc) Age (years) 31.4 ± 5.5 29.2 ± 5.6 37 ± 10.1 41.4 ± 6.2 Gender (M-F) 3/2 4/1 2/3 5/0 Tobacco use Users—n (%) 5 (100%) 4 (80%) 4 (80%) 5 (100%) Units per week µ 53.8 ± 50.2 36.4 ± 40.6 91.2 ± 91.0 147 ± 62.6 Use duration (years) µ 11.4 ± 4.9 10.4 ± 10.0 5.8 ± 6.8 24.2 ± 9.6 *,+,# Cannabis use Age first use (years) µ - 15.4 ± 2.5 - 17.2 ± 3 Units per week µ - 16.2 ± 10.2 - 19.4 ± 14.6 Use duration (years) µ - 12.2 ± 7.2 - 17.6 ± 11.5 Antipsychotic treatment Clozapine/Olanzapine - - 1/5 3/5 Aripiprazole - - 2/5 1/5 Paliperidone - - 1/5 1/5 Risperidone - - 1/5 - Data are shown as mean ± SD. * p < 0.05 vs. C/nc; +p < 0.05 vs. C/c; # p < 0.05 vs. SCZ/nc. (µ) Indicates continuous non-normally distributed variables. Table 2. Neuropsychological and clinical data of the different groups. Schizophrenia Control Subjects Cannabis Users Schizophrenia Non-Cannabis (C/nc) (C/c) Cannabis (SCZ/c) (SCZ/nc) Direct series score 8.8 ± 0.4 8.2 ± 1.8 7.4 ± 1.5 8.8 ± 2.2 Inverse series score µ 6.3 ± 1.3 8.2 ± 5.5 4.4 ± 1.5 8.8 ± 2.1 Verbal fluency score 26 ± 1.2 22.8 ± 4.5 18.8 ± 6.9 17.2 ± 5.3 *,+ Direct spatial span 7.4 ± 0.6 5.6 ± 1.9 5.2 ± 1.5 5.6 ± 1.7 Inverse spatial span 7.6 ± 2.9 5.8 ± 1.5 5.6 ± 3 5.6 ± 1.7 Emotional recognition 65.9 ± 8.6 73.8 ± 7.4 58.2 ± 22.5 58.9 ± 10.2 Global assessment of 99 ± 5 96 ± 5.5 67 ± 8.4 ***,+,++ 63 ± 2.7 ***,+,++ functioning µ Neurological soft signs µ 2.8 ± 1.1 5.2 ± 3.1 14.8 ± 8.6 **,+ 13.2 ± 2.3 **,+ Data are shown as mean ± SD of the mean. * p < 0.05, ** p < 0.01, *** p < 0.001 versus C/nc; + p < 0.05; ++ p < 0.01 versus C/c. (µ) Indicates continuous non-normally distributed variables. 3.3.1. Quantitative and Functional Proteomic Profile of ON Cells from SCZ Patients as Compared to C/nc The proteomic analysis of ON cells from C/c, SCZ/nc and SCZ/c as compared to C/nc showed 1185, 1209 and 1584 deregulated proteins, respectively. The Venn diagrams (Figure 1a) showed that in C/c, 550 proteins were upregulated (554 in SCZ/nc and 758 in SCZ/c). On the other hand, 635 proteins were downregulated in C/c (655 in SCZ/nc and 826 in SCZ/c); 371 proteins were commonly upregulated in C/c, SCZ/nc and SCZ/c; 211 were exclusively upregulated in SCZ/c and 44 in SCZ/nc. In addition, 474 proteins were commonly downregulated in C/c, SCZ/nc and SCZ/c; 169 were exclusively down- regulated in SCZ/c and 37 in SCZ/nc. In total, 845 proteins were commonly deregulated in C/c, SCZ/nc and SCZ/c as compared to C/nc. The pathway analysis of the 845 com- monly deregulated proteins (Figure 1b) revealed programmed cell death (37 proteins; log(10) FDR = 4.5), metabolism of proteins (225 proteins; log(10) FDR = 4.2), extracellular matrix organization (47 proteins; log(10) FDR = 3.5), vesicle-mediated transport (87 proteins; J. Pers. Med. 2021, 11, 160 6 of 15 log(10) FDR = 2.7), immune system (241 proteins; log(10) FDR = 2.0), protein localization J. Pers. Med. 2021, 11, x FOR PEER REVIEW 7 of 16 (22 proteins; log(10) FDR = 1.4) and cellular responses to external stimuli (67 proteins; log(10) FDR = 1.2). Figure Figure 1. 1. Characterization Characterization ofof thethe proteome proteome of of ON ON cells cells from from control control subjects subjects non-cannabis non-cannabis users users (C/nc), (C/nc), control control cannabis cannabis users (C/c), schizophrenia patients non-cannabis users (SCZ/nc) and schizophrenia patients cannabis users (SCZ/c). (a) users (C/c), schizophrenia patients non-cannabis users (SCZ/nc) and schizophrenia patients cannabis users (SCZ/c). (a) Proteomic alterations in C/c, SCZ/nc and SCZ/c as compared to C/nc. The upregulated and downregulated proteins in Proteomic alterations in C/c, SCZ/nc and SCZ/c as compared to C/nc. The upregulated and downregulated proteins in C/c-C/nc, SCZ/nc-C/nc and SCZ/c-C/nc are represented in two different Venn diagrams. (b) Enrichment pathway analysis C/c-C/nc, of the 845 SCZ/nc-C/nc and SCZ/c-C/nc proteins commonly are represented deregulated in twoand in C/c, SCZ/nc different SCZ/cVenn diagrams. (b) as compared Enrichment to C/nc. pathway analysis (c) Enrichment pathway ofanalysis the 845 proteins commonly deregulated in C/c, SCZ/nc and SCZ/c as compared to C/nc. (c) Enrichment of the exclusively deregulated proteins in C/c, (d) SCZ/nc (e) and SCZ/c versus C/nc. (f) Proteomic pathway analysisin alterations ofSCZ/nc the exclusively and SCZ/cderegulated as comparedproteins in C/c, to C/c. The (d) SCZ/nc (e)and upregulated anddownregulated SCZ/c versus C/nc. (f) Proteomic proteins alterations when SCZ/nc-C/c in SCZ/c-C/c and SCZ/nc and areSCZ/c as compared compared to C/c. in are represented The upregulated two and downregulated Venn diagrams. (g) Functional proteins when SCZ/nc-C/c characterization and SCZ/c-C/c of the 43 proteins arein deregulated SCZ/nc and compared the 154 proteins are represented in twoderegulated in SCZ/c Venn diagrams. as compared (g) Functional to C/c. (h) Functional characterization of the 43characterization of thein proteins deregulated 115 proteins SCZ/nc deregulated in SCZ/c as compared to SCZ/nc. Data from the pathway analysis are represented as minus the logarithm (10) of the FDR. J. Pers. Med. 2021, 11, 160 7 of 15 and the 154 proteins deregulated in SCZ/c as compared to C/c. (h) Functional characterization of the 115 proteins deregulated in SCZ/c as compared to SCZ/nc. Data from the pathway analysis are represented as minus the logarithm (10) of the FDR. The exclusive impact of cannabis on C/c was represented by 102 proteins (Figure 1c), wherein the main pathways affected included the immune system (79 proteins; log(10) FDR = 10.9), metabolism of RNA (21 proteins; log(10) FDR = 1.7), cellular responses to external stimuli (19 proteins; log(10) FDR = 1.6) and protein localization (7 proteins; log(10) FDR = 1.4). When comparing the two groups of SCZ vs. C/nc, we detected a substantial differ- ence (Figure 1d). SCZ/nc only showed 81 exclusively deregulated proteins, whereas 380 were deregulated in SCZ/c. The functional characterization of the 81 proteins deregu- lated in SCZ/nc indicated that most of them were involved in cell cycle (16 proteins; log(10) FDR = 2.0), DNA replication (6 proteins; log(10) FDR = 1.7) and cellular responses to external stimuli (11 proteins; log(10) FDR = 1.2). On the other hand, the 380 proteins deregulated in SCZ/c (Figure 1e) were associated with metabolism of RNA (61 proteins; log(10) FDR = 8.3), cellular responses to external stimuli (54 proteins; log(10) FDR = 7.1) and metabolism of proteins (104 proteins; log(10) FDR = 2.8). 3.3.2. Quantitative and Functional Proteomic Profile of ON Cells from SCZ Patients as Compared to C/nc To identify disease-exclusive protein markers, and tease out the effect of cannabis use, the exclusively deregulated proteins in SCZ/nc and SCZ/c were compared to C/c (Figure 1f). In total, 53 proteins were deregulated in SCZ/nc, whereas 164 proteins were differentially expressed in SCZ/c. In SCZ/nc, 30 proteins were upregulated and 101 in SCZ/c. On the other hand, 23 proteins were downregulated in SCZ/nc and 63 in SCZ/c. The Venn diagram revealed that only 10 proteins were commonly upregulated in SCZ/nc and SCZ/c as compared to C/c, wherein none were commonly downregulated. The functional characterization of the 43 exclusively deregulated proteins in SCZ/nc (Figure 1g) indicated that signal transduction (2 proteins; log(10) FDR = 0.6) and protein localization (2 proteins; log(10) FDR = 0.6) were the mainly enriched functions. The 154 exclusively deregulated proteins in SCZ/c unveiled metabolism of RNA (21 proteins; log(10) FDR = 1.4) as the only deregulated cascade. Thus, the profile of deregulated protein pathways in SCZ/c patients as compared to C/c is quantitatively larger and functionally different to the profile observed in SCZ/nc. To further evaluate the specific proteins involved SCZ/nc vs. SCZ/c, we assessed their functional differences. The results showed that 115 proteins were differentially expressed (Figure 1h). These proteins were associated with metabolism of RNA (28 proteins; log(10) FDR = 7.6), cellular responses to external stimuli (27 proteins; log(10) FDR = 7.2), metabolism of proteins (40 proteins; log(10) FDR = 2.7) and developmental biology (25 proteins; log(10) FDR = 2.4). A detailed list of the proteins involved in each comparison can be found in the Supplement Material. 3.3.3. Specific Proteins Markers of SCZ Depending on Cannabis Use To identify specific protein markers in ON cells of SCZ/nc, we designed a Venn diagram including the 81 proteins exclusively deregulated in SCZ/nc as compared to C/nc; the 115 proteins differently expressed in SCZ/nc as compared to SCZ/c; and the 43 proteins exclusively different between SCZ/nc and C/c (Figure 2a). In total, 3 proteins were commonly deregulated in these comparisons: (i) Cyclin-dependent-like kinase 5 (CDK5), (ii) LSM2 Homolog, U6 Small Nuclear RNA and mRNA Degradation-Associated (LSM2), and (iii) Tyrosine-tRNA ligase, mitochondrial (YARS2). The expression of CDK5 was significantly reduced in SCZ/nc (p < 0.01) as compared to C/nc. In addition, a significant increase in the expression of this protein was observed in SCZ/c when compared to SCZ/nc (p < 0.01) (Figure 2b). LSM2 was downregulated in SCZ/nc as compared to C/nc (p < 0.01) and vs. C/c (p < 0.01), and presented a different expression in SCZ/c vs. SCZ/nc J. Pers. Med. 2021, 11, 160 8 of 15 J. Pers. Med. 2021, 11, x FOR PEER REVIEW 9 of 16 (p < 0.01) (Figure 2c). YARS2 was significantly downregulated in SCZ/nc as compared to C/c (p < 0.01), and vs. SCZ/c (p < 0.01) (Figure 2d). Figure 2. Specific proteomic markers of schizophrenia patients non-cannabis users (SCZ/nc). (a) Venn diagram comparing SCZ/nc2.to Figure controlproteomic Specific subjects non-cannabis users (C/nc), patients markers of schizophrenia control cannabis users users non-cannabis (C/c)(SCZ/nc). and schizophrenia patients cannabis (a) Venn diagram compar- users ing (SCZ/c). SCZ/nc Graphical to control representation subjects of users non-cannabis CDK5(C/nc), (b), LSM2 (c) cannabis control and YARS2 (d),(C/c) users which were and commonlypatients schizophrenia deregulated in canna- bis users SCZ/c as(SCZ/c). compared Graphical to C/nc,representation C/c and SCZ/nc.of CDK5 ** p < (b), 0.01;LSM2 (c)C/nc; versus and YARS2 (d), which ++ p < 0.01 versus were commonly C/c; ## deregulated p < 0.01 versus in SCZ/nc. SCZ/c as compared to C/nc, C/c and SCZ/nc. ** p < 0.01; versus C/nc; ++ p < 0.01 versus C/c; ## p < 0.01 versus SCZ/nc. To identify specific protein markers in the ON cells of SCZ/c, the 380 proteins exclu- sively 3.4. deregulated Correlations in SCZ/c between SpecificasProtein compared to C/nc, Markers, the 115 Cannabis Useproteins differently and Cognitive expressed Performance in SCZ/c as compared to SCZ/nc, and the 154 proteins exclusively A Spearman correlation plot was designed to evaluate the association between differing between the SCZ/c and C/c, were compared (Figure 3a). Seven proteins were commonly expression levels of target proteins identified in SCZ/nc and in SCZ/c with cognitive deregulated in these three cannabis performance, comparisons; (i) CDK5 and use patterns, regulatory clinicalsubunit-associated signs in the entireprotein 3 (CDKRAP3); population. For pro- (ii) Haloacid teins related to dehalogenase-like SCZ/nc, the results hydrolase showeddomain-containing 5 (HDHD5); that higher LSM2 expression (iii)correlated was Exosome RNA helicase MTR4 (MTREX); (iv) Unconventional myosin-Ib (MYO1B); (v) 40S ribosomal with better GAF (r = 0.50, p < 0.05) and less NSS (r = −0.54, p < 0.05), whereas higher protein S20 (RPS20); (vi) Nucleoporin SEH1 (SEH1L), and (vii) DBIRD complex subunit YARS2 levels were associated with better attentional performance in the direct series ZNF326 (ZNF326). CDK5RAP3 was significantly upregulated in SCZ/c as compared to score test (r = 0.53, p < 0.05) (Data not shown). For proteins related to SCZ/c, MYO1B C/nc, C/c and SCZ/nc (p < 0.001, p < 0.05, p < 0.05, respectively) (Figure 3b). HDHD5 negatively correlated with plasmatic THC-COOH concentration (r = −0.80, p < 0.05) (Fig- was not significantly altered when the four groups were compared together. MTREX ure 4a), and ZNF326 positively correlated with cannabis use duration (r = 0.81, p < 0.01) was downregulated in SCZ/c (p < 0.001) and SCZ/nc (p < 0.05) as compared to C/nc. (Figure 4b). MTREX showed a positive correlation with attentional performance in the SCZ/c showed a significant downregulation in its expression as compared to C/c (p < 0.05) direct spatial span test (r = 0.52, p < 0.05) (Figure 4c), and with better verbal fluency scores (Figure 3c). MYO1B was significantly upregulated in SCZ/c as compared to C/nc (p < 0.01) (r = 0.50, p < 0.05) (Figure 4d). Moreover, MTREX positively correlated with GAF (r = 0.73, and C/c and SCZ/nc (p < 0.05) (Figure 3d). RPS20 was downregulated in SCZ/c only, as p < 0.001) (Figure 4e), and correlated negatively with NSS (r = −0.65, p < 0.01) (Figure 4f). compared to non-cannabis users (p < 0.001 vs. C/nc and p < 0.05 vs. SCZ/nc) (Figure 3e). In addition, ZNF326 negatively correlated with attention in the direct spatial span test (r SEH1L was significantly upregulated in SCZ/c as compared to C/nc (p < 0.01), C/c and =SCZ/nc −0.57, p(p< <0.01) 0.05)(Figure (Figure4j); 3f).with working Finally, ZNF326memory in the inverse was significantly spatial span upregulated test (r as in SCZ/c = −0.51, p < 0.05) (Figure 4g); and with emotional recognition compared to C/nc (p < 0.001) and C/c and SCZ/nc (p < 0.05) (Figure 3g). scores (r = −0.55, p < 0.05) (Figure 4h). Finally, higher ZNF326 correlated with lower GAF (r = −0.49, p < 0.05) (Figure 4i) and higher NSS (r = 0.52, p < 0.05) (Figure 4j). CDK5RAP3 correlated with GAF (r = −0.5, p < 0.05) and NSS (r = 0.47, p < 0.05) (Supplement Figure 2a,b). MYO1B levels were negatively correlated with GAF scores (r = −0.56, p < 0.05) and positively with NSS (r = 0.48, p < 0.05) (Supplementary Figure 2c,d). Finally, RPS20 was positively associated with J. Pers. Med. 2021, 11, 160 GAF (r = −0.50, p < 0.05) (Supplementary Figure 2e), while SEH1L showed a significant 9 of 15 negative correlation with this parameter (r = 0.48, p < 0.05) (Supplementary Figure 2f). Figure 3. Specific proteome markers of schizophrenia patients cannabis users (SCZ/c). (a) Venn diagram comparing SCZ/c to control Figure subjects non-cannabis 3. Specific users (C/nc), proteome markers schizophrenia of schizophrenia patients patients non-cannabis cannabis users (SCZ/nc), users (SCZ/c). (a) Venn and control diagram cannabis comparing SCZ/c to control users (C/c). subjects Graphical non-cannabisofusers representation (C/nc), (b), CDK5RAP3 schizophrenia MTREX (c),patients MYO1Bnon-cannabis (d), RPS20 (e),users SEH1L(SCZ/nc), (f), and and control ZNF326 (g), cannabis which were users (C/c). Graphical commonly representation deregulated in SCZ/c asof CDK5RAP3 compared (b), MTREX to C/nc, C/c and (c), MYO1B SCZ/c. Data(d), are RPS20 (e), SEH1L represented (f), and as the mean ± ZNF326 (g), which were commonly deregulated in SCZ/c as compared to C/nc, C/c and SCZ/c. Data are represented minimum to maximum peak intensity. * p < 0.05; ** p < 0.01; *** p < 0001 versus C/nc; + p < 0.05 versus C/c; # p < 0.05; as the mean versus ± minimum to maximum peak intensity. * p < 0.05; ** p < 0.01; *** p < 0001 versus C/nc; + p < 0.05 versus C/c; # p < SCZ/nc. 0.05; versus SCZ/nc. 3.4. Correlations between Specific Protein Markers, Cannabis Use and Cognitive Performance A Spearman correlation plot was designed to evaluate the association between the expression levels of target proteins identified in SCZ/nc and in SCZ/c with cognitive performance, cannabis use patterns, and clinical signs in the entire population. For proteins related to SCZ/nc, the results showed that higher LSM2 expression was correlated with better GAF (r = 0.50, p < 0.05) and less NSS (r = −0.54, p < 0.05), whereas higher YARS2 levels were associated with better attentional performance in the direct series score test (r = 0.53, p < 0.05) (Data not shown). For proteins related to SCZ/c, MYO1B negatively J. Pers. Med. 2021, 11, 160 10 of 15 correlated with plasmatic THC-COOH concentration (r = −0.80, p < 0.05) (Figure 4a), and ZNF326 positively correlated with cannabis use duration (r = 0.81, p < 0.01) (Figure 4b). MTREX showed a positive correlation with attentional performance in the direct spatial span test (r = 0.52, p < 0.05) (Figure 4c), and with better verbal fluency scores (r = 0.50, p < 0.05) (Figure 4d). Moreover, MTREX positively correlated with GAF (r = 0.73, p < 0.001) (Figure 4e), and correlated negatively with NSS (r = −0.65, p < 0.01) (Figure 4f). In addition, ZNF326 negatively correlated with attention in the direct spatial span test (r = −0.57, p < 0.01) (Figure 4j); with working memory in the inverse spatial span test (r = −0.51, p < 0.05) (Figure 4g); and with emotional recognition scores (r = −0.55, p < 0.05) (Figure 4h). Finally, higher ZNF326 correlated with lower GAF (r = −0.49, p < 0.05) (Figure 4i) and higher NSS (r = 0.52, p < 0.05) (Figure 4j). CDK5RAP3 correlated with GAF (r = −0.5, p < 0.05) and NSS (r = 0.47, p < 0.05) (Supplement Figure S2a,b). MYO1B levels were negatively correlated with GAF scores (r = −0.56, p < 0.05) and positively with NSS (r = 0.48, p < 0.05) (Supplementary Figure S2c,d). Finally, RPS20 was positively associated J. Pers. Med. 2021, 11, x FOR PEER REVIEW with GAF (r = −0.50, p < 0.05) (Supplementary Figure S2e), while SEH1L showed 11 ofa16 significant negative correlation with this parameter (r = 0.48, p < 0.05) (Supplementary Figure S2f). Figure 4. Association studies of protein markers, cannabis use and neurocognitive performance. (a) THC-COOH (ng/mL Figure and plasma) 4. Association studies of (b) MYO1B correlation. protein markers, and Use duration cannabis usespearman ZNF326 and neurocognitive correlation.performance. (c) MTREX and(a) THC-COOH direct spatial(ng/mL span plasma) and MYO1B correlation. (b) Use duration and ZNF326 spearman correlation. (c) MTREX and direct spatial span correlation. (d) MTREX and verbal fluency score correlation. (e) MTREX and GAF correlation. (f) MTREX and NSS correlation. (d) MTREX and verbal fluency score correlation. (e) MTREX and GAF correlation. (f) MTREX and NSS cor- correlation. (g) ZNF326 and direct spatial span correlation. (h) ZNF326 and inverse spatial span correlation. (i) ZNF326 and relation. (g) ZNF326 and direct spatial span correlation. (h) ZNF326 and inverse spatial span correlation. (i) ZNF326 and emotional emotionalrecognition recognitioncorrelation. correlation.(j)(j)ZNF326 ZNF326and andGAF GAFcorrelation. correlation.(k) (k)ZNF326 ZNF326and andNSS NSScorrelation. correlation.InInthe thecorrelation correlation matrixes, only the significant correlations are colored. Darker squares indicate lower p-values. matrixes, only the significant correlations are colored. Darker squares indicate lower p-values. 4. Discussion 4. Discussion In this study we revealed a quantitative shared impact of cannabis use and SCZ. These In this study we revealed a quantitative shared impact of cannabis use and SCZ. results are consistent with recent GWAS studies showing a significant common genetic risk These results are consistent with recent GWAS studies showing a significant common of SCZ and cannabis use [29]. Interestingly, when assessing the exclusively deregulated genetic risk of SCZ and cannabis use [29]. Interestingly, when assessing the exclusively deregulated proteins, we observed a functional selectivity among the different groups. Thus, the most significantly deregulated pathway in C/c was the immune system cas- cade, while in SCZ/c metabolism of RNA, cell responses to external stimuli and metabo- lism of proteins showed strong and significant alterations. In SCZ/nc, smaller but signif- J. Pers. Med. 2021, 11, 160 11 of 15 proteins, we observed a functional selectivity among the different groups. Thus, the most significantly deregulated pathway in C/c was the immune system cascade, while in SCZ/c metabolism of RNA, cell responses to external stimuli and metabolism of proteins showed strong and significant alterations. In SCZ/nc, smaller but significant deregulations were observed in cell cycle, DNA replication and cell responses to external stimuli cascades. These findings indicate that despite the common proteomic changes induced by cannabis use and SCZ, there might be specific changes prompted in the proteome depending on whether these two factors are present together or separately. In fact, when we subtracted the effect of cannabis comparing both groups of SCZ patients with C/c, we found a more similar profile between SCZ/nc and C/c, with small deregulations in protein localization and signal transduction, than in SCZ/c and C/c, which showed larger changes only in metabolism of RNA. Moreover, when the exclusively deregulated proteins were compared in both groups of SCZ patients, differences were obtained for metabolism of RNA and cellular responses to external stimuli. Again, these comparisons revealed functional differ- ences in the proteomic alterations induced by the presence of SCZ depending on whether patients use cannabis or not. Therefore, to identify specific proteins that could be proteomic markers of SCZ with or without the concomitant use of cannabis, we assessed separately the commonly deregulated proteins in SCZ/nc and in SCZ/c with respect to all the other groups. In SCZ/nc, we found three protein markers: CDK5, LSM2 and YARS2. These proteins were significantly downregulated in SCZ/nc, but not in SCZ/c. CDK5 is a cyclin-dependent kinase that controls the development of the central nervous system [34]. In the mature brain, it is associated with cognitive processes [35], and in vitro, it modulates oxidative stress responses [36]. According to the REACTOME database, it participates in developmental biology, which functionally differentiates SCZ/nc from SCZ/c. The observed reduction in CDK5 expression is consistent with previous studies showing a significant decline in CDK5 expression in postmortem brain samples of antipsychotic-treated SCZ patients, but not in drug-naïve individuals [37]. Markedly, our data indicate that cannabis use may counteract the effects of SCZ on CDK5. Further studies may shed light on whether cannabis use also opposes antipsychotic efficiency through this mechanism. Secondly, the downregulation of LSM2, which encodes for a key protein of the spliceosome [38], is in agreement with previously reported global changes in alternative splicing in SCZ patients [39]. LSM2 participates in metabolism of RNA, which differentiated SCZ/nc from SCZ/c, and presented a positive correlation with GAF and a negative correlation with NSS, indicating that it may be relevant for the clinical alterations observed in SCZ. Thirdly, we show a decrease in YARS2, a tyrosyl-tRNA synthetase located in the mitochondria, involved in metabolism of proteins, which differentiated SCZ/nc from SCZ/c. Its downregulation may lead to respiratory chain dysfunctions causing mitochondrial oxidative stress [40], which has been observed in postmortem brain samples from SCZ patients [41,42]. Moreover, lower levels of YARS2 could be correlated with worse attention performance, consistent with previous studies reporting that an oxidative imbalance increases negative symptoms’ severity in SCZ patients [43]. Once more, cannabis use seems to counteract the molecular changes in YARS2 in ON cells of SCZ/c, hypothetically via CB1 receptors expressed in the mitochondrial membranes [44], which have been proven to modulate respiration in this organelle [45]. In SCZ/c, we identified expression level changes in seven proteins: CDK5RAP3, HDHD5, MYO1B, MTREX, RPS20, SEH1L, and ZNF326. CDK5RAP3, a regulator of CDK5, was significantly upregulated in SCZ/c with respect to all other groups. It is genetically associated with SCZ [46], and participates in cell cycle. It modulates hippocampal aging because its transcriptional upregulation lowers neurogenesis [47]. Higher levels of this pro- tein were associated with less GAF and more NSS in the entire population. MYO1B, SEH1L and RPS20 were also deregulated, specifically in SCZ/c, and were associated with clinical outcomes (GAF and NSS). They are involved in the metabolism of proteins cascade, which includes the AKT-mTOR pathway that is thought to play a role in the interaction between J. Pers. Med. 2021, 11, 160 12 of 15 SCZ and cannabinoids exposure in humans [48] and in mice models [24]. In fact, MYO1B, SEH1L and RPS20 directly or indirectly interact with AKT-mTOR [49–51]. Interestingly, we found a negative correlation between MYO1B and plasmatic concentrations of THC-COOH in SCZ/c and C/c. THC-COOH is the last molecule arising from ∆9 -tetrahidrocannabinol degradation [52], which might point to the role of this unconventional myosin in cannabis metabolism. Finally, MTREX and ZNF326, which are involved in the metabolism of RNA cascade, were found to be deregulated in SCZ/c as compared to the rest of the groups. In fact, metabolism of RNA was consistently deregulated in SCZ/c in every functional com- parison. Although prior studies have found deregulations in ribosomal and RNA-related proteins in human samples from SCZ patients [53,54], our data revealed functional differ- ences in metabolism of RNA exclusively linked to the coexistence of cannabis consumption and SCZ. According to REACTOME, MTREX and ZNF326 are involved in the same step of pre-mRNA processing prior to protein translation. However, we found that MTREX was significantly downregulated in SCZ/nc and SCZ/c, while ZNF326 was upregulated in SCZ/c, but both of these alterations correlated with worse cognitive performance and more clinical signs of SCZ. Thus, a lower peak intensity of MTREX was associated with worse spatial attention, verbal fluency, and GAF scores, but with more NSS. On the other hand, higher peak intensity levels of ZNF326 were related with worse attention, working mem- ory, emotional recognition, and clinical signs. Moreover, cannabis use duration positively correlated with ZNF326 expression. These data indicate that either a decrease in MTREX or an increase in ZNF326 expression may have a negative effect on RNA metabolism, which could impact neurocognitive functioning. The results of our study need to be interpreted considering its strengths and limita- tions. The main strength of our data relies on the translational power of ON cell models to investigate biomarkers of neuropsychiatric disorders, such as SCZ. In addition, the evaluation of cannabis’ effects on separate cohorts of SCZ patients is unique, since there is an absence of studies comparing these populations. The main limitation of our study was the small number of subjects per group, thus, additional follow-up studies on a larger scale, including longitudinal and epidemiological studies, will be needed to further corroborate these findings. In summary, we revealed a quantitative shared effect of SCZ and cannabis use in the proteomic profile of ON cells, consistent with the genetic similarities previously de- scribed [29]. We found that cannabis use in controls has a significant impact on the immune system, while it alters metabolism of RNA and metabolism of proteins in SCZ patients. Additionally, SCZ/nc show small but significant differences in cell cycle, DNA replica- tion, proteins localization and signal transduction. In this group, we identified significant changes in two functionally relevant proteins (YARS2 and LSM2) associated to attentional processes, GAF and NSS scores. On the other hand, in ON cells from SCZ/c we found significant changes in six proteins relevant to clinical scores (MYO1B, MTREX, ZNF326, RPS20, CDK5RAP3, SEH1L) and two for cannabis use (MYO1B, ZNF326). Moreover, in SCZ/c, we identified consistent deregulation in metabolism of RNA in every comparison, and expression levels of MTREX and ZNF326, which take part in this pathway, correlated with several aspects of cognitive performance. Supplementary Materials: The following are available online at https://www.mdpi.com/2075-442 6/11/3/160/s1, Supplementary Figure S1, Supplementary Figure S2, Supplementary Material. Author Contributions: Experimental design, M.B.-C., P.R.; methodology, M.B.-C., K.A., M.L.-M., J.F.-I., E.S., P.R.; software, M.B.-C., K.A., M.L.-M., J.F.-I., E.S., P.R.; investigation, M.B.-C., C.F.-A., L.G., A.C.-R., P.R.; writing—original draft preparation, M.B.-C., P.R.; writing—review and editing, M.B.-C. and P.R.; supervision, R.d.l.T., V.P., E.S., P.R.; project administration, R.d.l.T., V.P., P.R.; funding acquisition, E.S., J.F.-I., P.R. All authors have read and agreed to the published version of the manuscript. J. Pers. Med. 2021, 11, 160 13 of 15 Funding: This work was supported by grants from DIUE de la Generalitat de Catalunya (2017-SGR- 1497), Instituto de Salud Carlos III, (PI18/00053 to P.R.) FIS-FEDER FUNDS, the Spanish Ministry of Science Innovation and Universities (PID2019-110356RB-I00 to ES and JFI), and Department of Economic and Business Development from Government of Navarra (ref. 0011-1411-2020-000028) to ES. The Proteomics Platform of Navarrabiomed is a member of Proteored (PRB3-ISCIII) and is supported by grant PT17/0019/009 to J.F.-I., of the PE I + D + I 2013-2016 funded by ISCIII and FEDER. We would like to acknowledge the Clinical Research Unit at the IMIM-Hospital del Mar Research Institute for their expertise in nasal brushing, especially Marta Pérez Otero and Dr. Ana M. Aldea Perona. The authors thank all the PRIDE Team for helping with the mass spectrometric data deposit in ProteomeXChange/PRIDE. The Clinical Neuroproteomics Unit of Navarrabiomed is a member of the Spanish Olfactory Network (ROE) (supported by grant RED2018-102662-T funded by Spanish Ministry of Science and Innovation). Institutional Review Board Statement: The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of the Parc de Salut Mar Clinical Research (CEIC; protocol code IMIMFTCL/SC/1, 13-3-2015). Informed Consent Statement: Informed consent was obtained from all subjects involved in the study. Data Availability Statement: In All MS raw data and search results files have been deposited into the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifiers PXD020739 (Reviewer account details: Username: reviewer15955@ebi.ac.uk; Password: UDhwtjVV). Conflicts of Interest: The authors declare no conflict of interest. References 1. Van Os, J.; Rutten, B.P.F.; Poulton, R. Gene-environment interactions in schizophrenia: Review of epidemiological findings and future directions. Schizophr. Bull. 2008, 34, 1066–1082. [CrossRef] 2. Nascimento, J.M.; Martins-De-Souza, D. The proteome of schizophrenia. NPJ Schizophr. 2015, 1, 14003. [CrossRef] 3. Jaros, J.A.J.; Martins-de-Souza, D.; Rahmoune, H.; Rothermundt, M.; Leweke, F.M.; Guest, P.C.; Bahn, S. Protein phosphorylation patterns in serum from schizophrenia patients and healthy controls. J. Proteom. 2012, 76, 43–55. [CrossRef] [PubMed] 4. Yang, Y.; Wan, C.; Li, H.; Zhu, H.; La, Y.; Xi, Z.; Chen, Y.; Jiang, L.; Feng, G.; He, L. Altered levels of acute phase proteins in the plasma of patients with schizophrenia. Anal. Chem. 2006, 78, 3571–3576. [CrossRef] [PubMed] 5. Huang, J.T.J.; Wang, L.; Prabakaran, S.; Wengenroth, M.; Lockstone, H.E.; Koethe, D.; Gerth, C.W.; Gross, S.; Schreiber, D.; Lilley, K.; et al. Independent protein-profiling studies show a decrease in apolipoprotein A1 levels in schizophrenia CSF, brain and peripheral tissues. Mol. Psychiatry 2008, 13, 1118–1128. [CrossRef] [PubMed] 6. Khandaker, G.M.; Dantzer, R. Is there a role for immune-to-brain communication in schizophrenia? Psychopharmacology 2016, 233, 1559–1573. [CrossRef] 7. Boiko, A.S.; Mednova, I.A.; Kornetova, E.G.; Semke, A.V.; Bokhan, N.A.; Loonen, A.J.M.; Ivanova, S.A. Apolipoprotein serum levels related to metabolic syndrome in patients with schizophrenia. Heliyon 2019, 5, e02033. [CrossRef] 8. Föcking, M.; Dicker, P.; English, J.A.; Schubert, K.O.; Dunn, M.J.; Cotter, D.R. Common proteomic changes in the hippocampus in schizophrenia and bipolar disorder and particular evidence for involvement of cornu ammonis regions 2 and 3. Arch. Gen. Psychiatry 2011, 68, 477–488. [CrossRef] [PubMed] 9. Martins-de-Souza, D.; Maccarrone, G.; Wobrock, T.; Zerr, I.; Gormanns, P.; Reckow, S.; Falkai, P.; Schmitt, A.; Turck, C.W. Proteome analysis of the thalamus and cerebrospinal fluid reveals glycolysis dysfunction and potential biomarkers candidates for schizophrenia. J. Psychiatr. Res. 2010, 44, 1176–1189. [CrossRef] 10. Wesseling, H.; Chan, M.K.; Tsang, T.M.; Ernst, A.; Peters, F.; Guest, P.C.; Holmes, E.; Bahn, S. A combined metabonomic and proteomic approach identifies frontal cortex changes in a chronic phencyclidine rat model in relation to human schizophrenia brain pathology. Neuropsychopharmacology 2013, 38, 2532–2544. [CrossRef] 11. Beasley, C.L.; Pennington, K.; Behan, A.; Wait, R.; Dunn, M.J.; Cotter, D. Proteomic analysis of the anterior cingulate cortex in the major psychiatric disorders: Evidence for disease-associated changes. Proteomics 2006, 6, 3414–3425. [CrossRef] 12. Saia-Cereda, V.M.; Cassoli, J.S.; Martins-de-Souza, D.; Nascimento, J.M. Psychiatric disorders biochemical pathways unraveled by human brain proteomics. Eur. Arch. Psychiatry Clin. Neurosci. 2017, 267, 3–17. [CrossRef] [PubMed] 13. Willard, S.L.; Sinclair, D.; Mirza, N.; Turetsky, B.; Berretta, S.; Hahn, C. Translational potential of olfactory mucosa for the study of neuropsychiatric illness. Transl. Psychiatry 2015, 5, e527-12. 14. Mackay-Sim, A. Concise review: Patient-derived olfactory stem cells: New models for brain diseases. Stem Cells 2012, 30, 2361–2365. [CrossRef] 15. Benítez-King, G.; Riquelme, A.; Ortíz-López, L.; Berlanga, C.; Rodríguez-Verdugo, M.S.; Romo, F.; Calixto, E.; Solís-Chagoyán, H.; Jímenez, M.; Montaño, L.M.; et al. A non-invasive method to isolate the neuronal linage from the nasal epithelium from schizophrenic and bipolar diseases. J. Neurosci. Methods 2011, 201, 35–45. [CrossRef] J. Pers. Med. 2021, 11, 160 14 of 15 16. Galindo, L.; Moreno, E.; Lopez-Armenta, F.; Guinart, D.; Cuenca-Royo, A.; Izquierdo-Serra, M.; Xicota, L.; Fernandez, C.; Menoyo, E.; Fernandez-Fernandez, J.M.; et al. Cannabis Users Show Enhanced Expression of CB1-5HT2A Receptor Heteromers in Olfactory Neuroepithelium Cells. Mol. Neurobiol. 2018, 55, 6347–6361. [CrossRef] 17. Ortiz-Medina, M.B.; Perea, M.; Torales, J.; Ventriglio, A.; Vitrani, G.; Aguilar, L.; Roncero, C. Cannabis consumption and psychosis or schizophrenia development. Int. J. Soc. Psychiatry 2018, 64, 690–704. [CrossRef] 18. Di Forti, M.; Quattrone, D.; Freeman, T.P.; Tripoli, G.; Gayer-Anderson, C.; Quigley, H.; Rodriguez, V.; Jongsma, H.E.; Ferraro, L.; La Cascia, C.; et al. The contribution of cannabis use to variation in the incidence of psychotic disorder across Europe (EU-GEI): A multicentre case-control study. Lancet Psychiatry 2019, 6, 427–436. [CrossRef] 19. Muguruza, C.; Morentin, B.; Meana, J.J.; Alexander, S.P.H.; Callado, L.F. Endocannabinoid system imbalance in the postmortem prefrontal cortex of subjects with schizophrenia. J. Psychopharmacol. 2019, 33, 1132–1140. [CrossRef] 20. Bioque, M.; García-Bueno, B.; MacDowell, K.S.; Meseguer, A.; Saiz, P.A.; Parellada, M.; Gonzalez-Pinto, A.; Rodriguez-Jimenez, R.; Lobo, A.; Leza, J.C.; et al. Peripheral endocannabinoid system dysregulation in first-episode psychosis. Neuropsychopharmacology 2013, 38, 2568–2577. [CrossRef] [PubMed] 21. Twitchell, W.; Brown, S.; Mackie, K. Cannabinoids inhibit n- and p/q-type calcium channels in cultured rat hippocampal neurons. J. Neurophysiol. 1997, 78, 43–50. [CrossRef] 22. Fisyunov, A.; Tsintsadze, V.; Min, R.; Burnashev, N.; Lozovaya, N. Cannabinoids modulate the P-type high-voltage-activated calcium currents in Purkinje neurons. J. Neurophysiol. 2006, 96, 1267–1277. [CrossRef] 23. Miller, M.L.; Chadwick, B.; Dickstein, D.L.; Purushothaman, I.; Egervari, G.; Rahman, T.; Tessereau, C.; Hof, P.R.; Roussos, P.; Shen, L.; et al. Adolescent exposure to ∆ 9 -tetrahydrocannabinol alters the transcriptional trajectory and dendritic architecture of prefrontal pyramidal neurons. Mol. Psychiatry 2019, 24, 588–600. [CrossRef] [PubMed] 24. Ibarra-Lecue, I.; Mollinedo-Gajate, I.; Meana, J.J.; Callado, L.F.; Diez-Alarcia, R.; Urigüen, L. Chronic cannabis promotes pro- hallucinogenic signaling of 5-HT2A receptors through Akt/mTOR pathway. Neuropsychopharmacology 2018, 43, 2028–2035. [CrossRef] [PubMed] 25. Jayanthi, S.; Buie, S.; Moore, S.; Herning, R.I.; Better, W.; Wilson, N.M.; Contoreggi, C.; Cadet, J.L. Heavy marijuana users show increased serum apolipoprotein C-III levels: Evidence from proteomic analyses. Mol. Psychiatry 2010, 15, 101–112. [CrossRef] 26. Saladrigas-Manjón, S.; Dučić, T.; Galindo, L.; Fernández-Avilés, C.; Pérez, V.; de la Torre, R.; Robledo, P. Effects of cannabis use on the protein and lipid profile of olfactory neuroepithelium cells from schizophrenia patients studied by synchrotron-based FTIR spectroscopy. Biomolecules 2020, 10, 329. [CrossRef] 27. Guinart, D.; Moreno, E.; Galindo, L.; Cuenca-Royo, A.; Barrera-Conde, M.; Pérez, E.J.; Fernández-Avilés, C.; Correll, C.U.; Canela, E.I.; Casadó, V.; et al. Altered Signaling in CB1R-5-HT2AR Heteromers in Olfactory Neuroepithelium Cells of Schizophre- nia Patients is Modulated by Cannabis Use. Schizophr. Bull. 2020, 46, 1547–1557. [CrossRef] 28. Lachén-Montes, M.; Mendizuri, N.; Ausín, K.; Pérez-Mediavilla, A.; Azkargorta, M.; Iloro, I.; Elortza, F.; Kondo, H.; Ohigashi, I.; Ferrer, I.; et al. Smelling the Dark Proteome: Functional Characterization of PITH Domain-Containing Protein 1 (C1orf128) in Olfactory Metabolism. J. Proteome Res. 2020, 19, 4826–4843. [CrossRef] 29. Pasman, J.A.; Verweij, K.J.H.; Gerring, Z.; Stringer, S.; Sanchez-roige, S.; Treur, J.L.; Abdellaoui, A.; Nivard, M.G.; Baselmans, B.M.L.; Ong, J.; et al. A Causal Influence of Schizophrenia. Nat. Neurosci. 2018, 21, 1161–1170. [CrossRef] 30. Gillet, L.C.; Navarro, P.; Tate, S.; Röst, H.; Selevsek, N.; Reiter, L.; Bonner, R.; Aebersold, R. Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: A new concept for consistent and accurate proteome analysis. Mol. Cell. Proteom. 2012, 11, 1–17. [CrossRef] 31. Jassal, B.; Matthews, L.; Viteri, G.; Gong, C.; Lorente, P.; Fabregat, A.; Sidiropoulos, K.; Cook, J.; Gillespie, M.; Haw, R.; et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020, 48, D498–D503. [CrossRef] 32. Hulsen, T.; de Vlieg, J.; Alkema, W. BioVenn—A web application for the comparison and visualization of biological lists using area-proportional Venn diagrams. BMC Genom. 2008, 9, 488. [CrossRef] 33. Geyer, P.E.; Holdt, L.M.; Teupser, D.; Mann, M. Revisiting biomarker discovery by plasma proteomics. Mol. Syst. Biol. 2017, 13, 942. [CrossRef] 34. Cruz, J.C.; Tseng, H.C.; Goldman, J.A.; Shih, H.; Tsai, L.H. Aberrant Cdk5 activation by p25 triggers pathological events leading to neurodegeneration and neurofibrillary tangles. Neuron 2003, 40, 471–483. [CrossRef] 35. Nishi, A.; Bibb, J.A.; Matsuyama, S.; Hamada, M.; Higashi, H.; Nairn, A.C.; Greengard, P. Regulation of DARPP-32 dephospho- rylation at PKA- and Cdk5-sites by NMDA and AMPA receptors: Distinct roles of calcineurin and protein phosphatase-2A. J. Neurochem. 2002, 832–841. [CrossRef] [PubMed] 36. Jimenez-Blasco, D.; Santofimia-Castanõ, P.; Gonzalez, A.; Almeida, A.; Bolanõs, J.P. Astrocyte NMDA receptors’ activity sustains neuronal survival through a Cdk5-Nrf2 pathway. Cell Death Differ. 2015, 22, 1877–1889. [CrossRef] [PubMed] 37. Ramos-Miguel, A.; Javier Meana, J.; García-Sevilla, J.A. Cyclin-dependent kinase-5 and p35/p25 activators in schizophrenia and major depression prefrontal cortex: Basal contents and effects of psychotropic medications. Int. J. Neuropsychopharmacol. 2013, 16, 683–689. [CrossRef] 38. Bertram, K.; Agafonov, D.E.; Dybkov, O.; Haselbach, D.; Leelaram, M.N.; Will, C.L.; Urlaub, H.; Kastner, B.; Lührmann, R.; Stark, H. Cryo-EM Structure of a Pre-catalytic Human Spliceosome Primed for Activation. Cell 2017, 170, 701–713.e11. [CrossRef] 39. Reble, E.; Dineen, A.; Barr, C.L. The contribution of alternative splicing to genetic risk for psychiatric disorders. Genes Brain Behav. 2018, 17, 1–12. [CrossRef] [PubMed] J. Pers. Med. 2021, 11, 160 15 of 15 40. Riley, L.G.; Cooper, S.; Hickey, P.; Rudinger-Thirion, J.; McKenzie, M.; Compton, A.; Lim, S.C.; Thorburn, D.; Ryan, M.T.; Giegé, R.; et al. Mutation of the mitochondrial tyrosyl-tRNA synthetase gene, YARS2, causes myopathy, lactic acidosis, and sideroblastic anemia—MLASA syndrome. Am. J. Hum. Genet. 2010, 87, 52–59. [CrossRef] [PubMed] 41. Rezin, G.T.; Amboni, G.; Zugno, A.I.; Quevedo, J.; Streck, E.L. Mitochondrial dysfunction and psychiatric disorders. Neurochem. Res. 2009, 34, 1021–1029. [CrossRef] [PubMed] 42. Maurer, I.; Zierz, S.; Möller, H.J. Evidence for a mitochondrial oxidative phosphorylation defect in brains from patients with schizophrenia. Schizophr. Res. 2001, 48, 125–136. [CrossRef] 43. Gunes, M.; Altindag, A.; Bulut, M.; Demir, S.; Ibiloglu, A.O.; Kaya, M.C.; Atli, A.; Aksoy, N. Oxidative metabolism may be associated with negative symptoms in schizophrenia. Psychiatry Clin. Psychopharmacol. 2017, 27, 54–61. [CrossRef] 44. Melser, S.; Zottola, A.C.P.; Serrat, R.; Puente, N.; Grandes, P.; Marsicano, G.; Hebert-Chatelain, E. Functional Analysis of Mitochondrial CB1 Cannabinoid Receptors (mtCB1) in the Brain, 1st ed.; Elsevier Inc.: Amsterdam, The Netherlands, 2017; Volume 593. 45. Bénard, G.; Massa, F.; Puente, N.; Lourenço, J.; Bellocchio, L.; Soria-Gómez, E.; Matias, I.; Delamarre, A.; Metna-Laurent, M.; Cannich, A.; et al. Mitochondrial CB 1 receptors regulate neuronal energy metabolism. Nat. Neurosci. 2012, 15, 558–564. [CrossRef] [PubMed] 46. Camargo, L.M.; Collura, V.; Rain, J.C.; Mizuguchi, K.; Hermjakob, H.; Kerrien, S.; Bonnert, T.P.; Whiting, P.J.; Brandon, N.J. Disrupted in Schizophrenia 1 interactome: Evidence for the close connectivity of risk genes and a potential synaptic basis for schizophrenia. Mol. Psychiatry 2007, 12, 74–86. [CrossRef] 47. Shetty, G.A.; Hattiangady, B.; Shetty, A.K. Neural stem cell- and neurogenesis-related gene expression profiles in the young and aged dentate gyrus. Age 2013, 35, 2165–2176. [CrossRef] [PubMed] 48. Di Forti, M.; Iyegbe, C.; Sallis, H.; Kolliakou, A.; Falcone, M.A.; Paparelli, A.; Sirianni, M.; La Cascia, C.; Stilo, S.A.; Marques, T.R.; et al. Confirmation that the AKT1 (rs2494732) genotype influences the risk of psychosis in cannabis users. Biol. Psychiatry 2012, 72, 811–816. [CrossRef] 49. Platani, M.; Samejima, I.; Samejima, K.; Kanemaki, M.T.; Earnshaw, W.C. Seh1 targets GATOR2 and Nup153 to mitotic chromo- somes. J. Cell Sci. 2018, 131, jcs213140. [CrossRef] 50. Salas-Cortes, L.; Ye, F.; Tenza, D.; Wilhelm, C.; Theos, A.; Louvard, D.; Raposo, G.; Coudrier, E. Myosin lb modulates the morphology and the protein transport within multi-vesicular sorting endosomes. J. Cell Sci. 2005, 118, 4823–4832. [CrossRef] 51. Yu, Y.; Xiong, Y.; Montani, J.P.; Yang, Z.; Ming, X.F. Arginase-II activates mTORC1 through myosin-1b in vascular cell senescence and apoptosis. Cell Death Dis. 2018, 9, 313. [CrossRef] 52. Sharma, P.; Murthy, P.; Bharath, M.M.S. Chemistry, metabolism, and toxicology of cannabis: Clinical implications. Iran. J. Psychiatry 2012, 7, 149–156. [PubMed] 53. English, J.A.; Fan, Y.; Föcking, M.; Lopez, L.M.; Hryniewiecka, M.; Wynne, K.; Dicker, P.; Matigian, N.; Cagney, G. Reduced protein synthesis in schizophrenia patient-derived olfactory cells. Transl. Psychiatry 2015, 5, e663. [CrossRef] [PubMed] 54. Ibarra-Lecue, I.; Diez-Alarcia, R.; Morentin, B.; Meana, J.J.; Callado, L.F.; Urigüen, L. Ribosomal Protein S6 Hypofunction in Postmortem Human Brain Links mTORC1-Dependent Signaling and Schizophrenia. Front. Pharmacol. 2020, 11, 344. [CrossRef] [PubMed]