Next Article in Journal
The Potential of Mesenchymal Stromal Cells in Neuroblastoma Therapy for Delivery of Anti-Cancer Agents and Hematopoietic Recovery
Next Article in Special Issue
Genetic Polymorphisms of 5-HT Receptors and Antipsychotic-Induced Metabolic Dysfunction in Patients with Schizophrenia
Previous Article in Journal
Towards Personalised Contrast Injection: Artificial-Intelligence-Derived Body Composition and Liver Enhancement in Computed Tomography
Previous Article in Special Issue
Increased Heschl’s Gyrus Duplication in Schizophrenia Spectrum Disorders: A Cross-Sectional MRI Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cannabis Use Induces Distinctive Proteomic Alterations in Olfactory Neuroepithelial Cells of Schizophrenia Patients

1
Integrative Pharmacology and Systems Neuroscience, Neuroscience Research Program, IMIM-Hospital del Mar Research Institute, 08003 Barcelona, Spain
2
Department of Experimental and Health Sciences, University Pompeu Fabra, 08003 Barcelona, Spain
3
Clinical Neuroproteomics Unit, Proteomics Platform, Navarrabiomed, Complejo Hospitalario de Navarra (CHN), Universidad Pública de Navarra (UPNA), IdisNA, Proteored-ISCIII, 31006 Pamplona, Spain
4
Department of Psychiatry, University of Cambridge, Cambridgeshire and Peterborough NHS Foundation Trust, Cambridge CB2 1TN, UK
5
Neuropsychiatry and Addictions Institute (INAD) of Parc de Salut Mar, 08003 Barcelona and CIBER de Salud Mental, Spain
6
Department of Psychiatry and Legal Medicine, Autonomous University of Barcelona, 08193 Barcelona, Spain
7
Centro de Investigación Biomédica en Red de Fisiopatología de la Obesidad y Nutrición, Instituto de Salud Carlos III, 28029 Madrid, Spain
*
Author to whom correspondence should be addressed.
Submission received: 14 January 2021 / Revised: 17 February 2021 / Accepted: 20 February 2021 / Published: 25 February 2021

Abstract

:
A close epidemiological link has been reported between cannabis use and schizophrenia (SCZ). However, biochemical markers in living humans related to the impact of cannabis in this disease are still missing. Olfactory neuroepithelium (ON) cells express neural features and offer a unique advantage to study biomarkers of psychiatric diseases. The aim of our study was to find exclusively deregulated proteins in ON cells of SCZ patients with and without a history of cannabis use. Thus, we compared the proteomic profiles of SCZ non-cannabis users (SCZ/nc) and 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 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 proteins in each population were then correlated with cognitive performance and clinical scores. In SCZ/c, the expression levels of 2 proteins involved in the metabolism of RNA (MTREX and ZNF326) 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 might be relevant to understand the influence of cannabis use on SCZ.

1. Introduction

The complexity of schizophrenia (SCZ), involving intricate interactions between environmental factors and genetics, hinders the identification of the molecular mechanisms underlying its development [1]. State of the art techniques are being used to find biomarkers and ease clinical practice. In this sense, proteomic tools, allowing the study of a large number of proteins at a time, can provide an integrated picture of the biological dysfunctions associated with SCZ [2]. Serum, plasma and cerebrospinal fluid studies reveal an effect on proteins related to the innate immune system and a highly reported downregulation of Apolipoprotein A-I1 (APOA1) [3,4,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,9,10]. However, other neuropsychiatric conditions 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,14,15,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 functions. Indeed, protein metabolism via AKT/mTOR, cytoskeleton organization, and calcium signalling are affected in human samples from SCZ patients [2], and after cannabinoid exposure 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-5HT2AR) 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 patients 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 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 conditions 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 neuropsychological 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 concentrations 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 spectrometric 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 rate (FDR) (1% Global FDR or better). Data analysis was performed by sequential window 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 extract 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 hierarchically 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 calculations 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 neuropsychological 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 by SCZ or by SCZ plus cannabis use. For the summary followed in the proteomic data analysis, see Supplementary Figure S1.

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 downregulated 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 commonly 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; log(10)FDR = 2.7), immune system (241 proteins; log(10)FDR = 2.0), protein localization (22 proteins; log(10)FDR = 1.4) and cellular responses to external stimuli (67 proteins; log(10)FDR = 1.2).
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 difference (Figure 1d). SCZ/nc only showed 81 exclusively deregulated proteins, whereas 380 were deregulated in SCZ/c. The functional characterization of the 81 proteins deregulated 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 (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).
To identify specific protein markers in the ON cells of SCZ/c, the 380 proteins exclusively deregulated in SCZ/c as compared to C/nc, the 115 proteins differently expressed in SCZ/c as compared to SCZ/nc, and the 154 proteins exclusively differing between SCZ/c and C/c, were compared (Figure 3a). Seven proteins were commonly deregulated in these three comparisons; (i) CDK5 regulatory subunit-associated protein 3 (CDKRAP3); (ii) Haloacid dehalogenase-like hydrolase domain-containing 5 (HDHD5); (iii) Exosome RNA helicase MTR4 (MTREX); (iv) Unconventional myosin-Ib (MYO1B); (v) 40S ribosomal protein S20 (RPS20); (vi) Nucleoporin SEH1 (SEH1L), and (vii) DBIRD complex subunit ZNF326 (ZNF326). CDK5RAP3 was significantly upregulated in SCZ/c as compared to C/nc, C/c and SCZ/nc (p < 0.001, p < 0.05, p < 0.05, respectively) (Figure 3b). HDHD5 was not significantly altered when the four groups were compared together. MTREX was downregulated in SCZ/c (p < 0.001) and SCZ/nc (p < 0.05) as compared to C/nc. SCZ/c showed a significant downregulation in its expression as compared to C/c (p < 0.05) (Figure 3c). MYO1B was significantly upregulated in SCZ/c as compared to C/nc (p < 0.01) and C/c and SCZ/nc (p < 0.05) (Figure 3d). RPS20 was downregulated in SCZ/c only, as compared to non-cannabis users (p < 0.001 vs. C/nc and p < 0.05 vs. SCZ/nc) (Figure 3e). SEH1L was significantly upregulated in SCZ/c as compared to C/nc (p < 0.01), C/c and SCZ/nc (p < 0.05) (Figure 3f). Finally, ZNF326 was significantly upregulated in SCZ/c as compared to C/nc (p < 0.001) and C/c and SCZ/nc (p < 0.05) (Figure 3g).

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 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 with GAF (r = −0.50, p < 0.05) (Supplementary Figure S2e), while SEH1L showed a significant negative correlation with this parameter (r = 0.48, p < 0.05) (Supplementary Figure S2f).

4. Discussion

In this study we revealed a quantitative shared impact of cannabis use and SCZ. These results are consistent with recent GWAS studies showing a significant common 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 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 differences 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 protein 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 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,50,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 comparison. Although prior studies have found deregulations in ribosomal and RNA-related proteins in human samples from SCZ patients [53,54], our data revealed functional differences 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 memory, 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 limitations. 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 described [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 replication, 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://0-www-mdpi-com.brum.beds.ac.uk/2075-4426/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.

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: [email protected]; 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. [Google Scholar] [CrossRef] [Green Version]
  2. Nascimento, J.M.; Martins-De-Souza, D. The proteome of schizophrenia. NPJ Schizophr. 2015, 1, 14003. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed]
  6. Khandaker, G.M.; Dantzer, R. Is there a role for immune-to-brain communication in schizophrenia? Psychopharmacology 2016, 233, 1559–1573. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar]
  14. Mackay-Sim, A. Concise review: Patient-derived olfactory stem cells: New models for brain diseases. Stem Cells 2012, 30, 2361–2365. [Google Scholar] [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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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 Schizophrenia Patients is Modulated by Cannabis Use. Schizophr. Bull. 2020, 46, 1547–1557. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  33. Geyer, P.E.; Holdt, L.M.; Teupser, D.; Mann, M. Revisiting biomarker discovery by plasma proteomics. Mol. Syst. Biol. 2017, 13, 942. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  35. Nishi, A.; Bibb, J.A.; Matsuyama, S.; Hamada, M.; Higashi, H.; Nairn, A.C.; Greengard, P. Regulation of DARPP-32 dephosphorylation at PKA- and Cdk5-sites by NMDA and AMPA receptors: Distinct roles of calcineurin and protein phosphatase-2A. J. Neurochem. 2002, 832–841. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar]
  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. [Google Scholar] [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. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef]
  49. Platani, M.; Samejima, I.; Samejima, K.; Kanemaki, M.T.; Earnshaw, W.C. Seh1 targets GATOR2 and Nup153 to mitotic chromosomes. J. Cell Sci. 2018, 131, jcs213140. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  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. [Google Scholar] [CrossRef] [Green Version]
  52. Sharma, P.; Murthy, P.; Bharath, M.M.S. Chemistry, metabolism, and toxicology of cannabis: Clinical implications. Iran. J. Psychiatry 2012, 7, 149–156. [Google Scholar] [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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Characterization of the proteome of ON cells from control subjects non-cannabis users (C/nc), control cannabis 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 C/c-C/nc, SCZ/nc-C/nc and SCZ/c-C/nc are represented in two different Venn diagrams. (b) Enrichment pathway analysis of the 845 proteins commonly deregulated in C/c, SCZ/nc and SCZ/c as compared to C/nc. (c) Enrichment pathway analysis of the exclusively deregulated proteins in C/c, (d) SCZ/nc (e) and SCZ/c versus C/nc. (f) Proteomic alterations in SCZ/nc and SCZ/c as compared to C/c. The upregulated and downregulated proteins when SCZ/nc-C/c and SCZ/c-C/c are compared are represented in two Venn diagrams. (g) Functional characterization of the 43 proteins deregulated in SCZ/nc 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.
Figure 1. Characterization of the proteome of ON cells from control subjects non-cannabis users (C/nc), control cannabis 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 C/c-C/nc, SCZ/nc-C/nc and SCZ/c-C/nc are represented in two different Venn diagrams. (b) Enrichment pathway analysis of the 845 proteins commonly deregulated in C/c, SCZ/nc and SCZ/c as compared to C/nc. (c) Enrichment pathway analysis of the exclusively deregulated proteins in C/c, (d) SCZ/nc (e) and SCZ/c versus C/nc. (f) Proteomic alterations in SCZ/nc and SCZ/c as compared to C/c. The upregulated and downregulated proteins when SCZ/nc-C/c and SCZ/c-C/c are compared are represented in two Venn diagrams. (g) Functional characterization of the 43 proteins deregulated in SCZ/nc 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.
Jpm 11 00160 g001
Figure 2. Specific proteomic markers of schizophrenia patients non-cannabis users (SCZ/nc). (a) Venn diagram comparing SCZ/nc to control subjects non-cannabis users (C/nc), control cannabis users (C/c) and schizophrenia patients cannabis users (SCZ/c). Graphical representation of CDK5 (b), LSM2 (c) and YARS2 (d), which were commonly deregulated in 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.
Figure 2. Specific proteomic markers of schizophrenia patients non-cannabis users (SCZ/nc). (a) Venn diagram comparing SCZ/nc to control subjects non-cannabis users (C/nc), control cannabis users (C/c) and schizophrenia patients cannabis users (SCZ/c). Graphical representation of CDK5 (b), LSM2 (c) and YARS2 (d), which were commonly deregulated in 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.
Jpm 11 00160 g002
Figure 3. Specific proteome markers of schizophrenia patients cannabis users (SCZ/c). (a) Venn diagram comparing SCZ/c to control subjects non-cannabis users (C/nc), schizophrenia patients non-cannabis users (SCZ/nc), and control cannabis users (C/c). Graphical representation of CDK5RAP3 (b), MTREX (c), MYO1B (d), RPS20 (e), SEH1L (f), and ZNF326 (g), which were commonly deregulated in SCZ/c as compared to C/nc, C/c and SCZ/c. Data are represented as the mean ± 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; versus SCZ/nc.
Figure 3. Specific proteome markers of schizophrenia patients cannabis users (SCZ/c). (a) Venn diagram comparing SCZ/c to control subjects non-cannabis users (C/nc), schizophrenia patients non-cannabis users (SCZ/nc), and control cannabis users (C/c). Graphical representation of CDK5RAP3 (b), MTREX (c), MYO1B (d), RPS20 (e), SEH1L (f), and ZNF326 (g), which were commonly deregulated in SCZ/c as compared to C/nc, C/c and SCZ/c. Data are represented as the mean ± 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; versus SCZ/nc.
Jpm 11 00160 g003
Figure 4. Association studies of protein markers, cannabis use and neurocognitive performance. (a) THC-COOH (ng/mL 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. (g) ZNF326 and direct spatial span correlation. (h) ZNF326 and inverse spatial span correlation. (i) ZNF326 and emotional recognition correlation. (j) ZNF326 and GAF correlation. (k) ZNF326 and NSS correlation. In the correlation matrixes, only the significant correlations are colored. Darker squares indicate lower p-values.
Figure 4. Association studies of protein markers, cannabis use and neurocognitive performance. (a) THC-COOH (ng/mL 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. (g) ZNF326 and direct spatial span correlation. (h) ZNF326 and inverse spatial span correlation. (i) ZNF326 and emotional recognition correlation. (j) ZNF326 and GAF correlation. (k) ZNF326 and NSS correlation. In the correlation matrixes, only the significant correlations are colored. Darker squares indicate lower p-values.
Jpm 11 00160 g004
Table 1. Demographic characteristics of the different groups included in the study.
Table 1. Demographic characteristics of the different groups included in the study.
Control Subjects (C/nc)Cannabis Users
(C/c)
Schizophrenia Non-Cannabis (SCZ/nc)Schizophrenia Cannabis (SCZ/c)
Age (years)31.4 ± 5.529.2 ± 5.637 ± 10.141.4 ± 6.2
Gender (M-F)3/24/12/35/0
Tobacco use
Users—n (%)5 (100%)4 (80%)4 (80%)5 (100%)
Units per week µ53.8 ± 50.236.4 ± 40.691.2 ± 91.0147 ± 62.6
Use duration (years) µ11.4 ± 4.910.4 ± 10.05.8 ± 6.824.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/53/5
Aripiprazole--2/51/5
Paliperidone--1/51/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.
Table 2. Neuropsychological and clinical data of the different groups.
Control Subjects
(C/nc)
Cannabis Users
(C/c)
Schizophrenia Non-Cannabis (SCZ/nc)Schizophrenia Cannabis (SCZ/c)
Direct series score8.8 ± 0.48.2 ± 1.87.4 ± 1.58.8 ± 2.2
Inverse series score µ6.3 ± 1.38.2 ± 5.54.4 ± 1.58.8 ± 2.1
Verbal fluency score26 ± 1.222.8 ± 4.518.8 ± 6.917.2 ± 5.3 *,+
Direct spatial span7.4 ± 0.65.6 ± 1.95.2 ± 1.55.6 ± 1.7
Inverse spatial span7.6 ± 2.95.8 ± 1.55.6 ± 35.6 ± 1.7
Emotional recognition65.9 ± 8.673.8 ± 7.458.2 ± 22.558.9 ± 10.2
Global assessment of functioning µ99 ± 596 ± 5.567 ± 8.4 ***,+,++63 ± 2.7 ***,+,++
Neurological soft signs µ2.8 ± 1.15.2 ± 3.114.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.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barrera-Conde, M.; Ausin, K.; Lachén-Montes, M.; Fernández-Irigoyen, J.; Galindo, L.; Cuenca-Royo, A.; Fernández-Avilés, C.; Pérez, V.; de la Torre, R.; Santamaría, E.; et al. Cannabis Use Induces Distinctive Proteomic Alterations in Olfactory Neuroepithelial Cells of Schizophrenia Patients. J. Pers. Med. 2021, 11, 160. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030160

AMA Style

Barrera-Conde M, Ausin K, Lachén-Montes M, Fernández-Irigoyen J, Galindo L, Cuenca-Royo A, Fernández-Avilés C, Pérez V, de la Torre R, Santamaría E, et al. Cannabis Use Induces Distinctive Proteomic Alterations in Olfactory Neuroepithelial Cells of Schizophrenia Patients. Journal of Personalized Medicine. 2021; 11(3):160. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030160

Chicago/Turabian Style

Barrera-Conde, Marta, Karina Ausin, Mercedes Lachén-Montes, Joaquín Fernández-Irigoyen, Liliana Galindo, Aida Cuenca-Royo, Cristina Fernández-Avilés, Víctor Pérez, Rafael de la Torre, Enrique Santamaría, and et al. 2021. "Cannabis Use Induces Distinctive Proteomic Alterations in Olfactory Neuroepithelial Cells of Schizophrenia Patients" Journal of Personalized Medicine 11, no. 3: 160. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030160

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

Article Metrics

Back to TopTop