Amy L. Cherry
Michael J. Wheeler
Karolina Mathisova3
Mathieu Di Miceli- 1Worcester Biomedical Research Group, School of Science and the Environment, University of Worcester, Worcester, United Kingdom
- 2Sustainable Environments Research Group, School of Science and the Environment University of Worcester, Worcester, United Kingdom
- 3School of Science and the Environment University of Worcester, Worcester, United Kingdom
Introduction: The endocannabinoid (eCB) system is named after the discovery that endogenous cannabinoids bind to the same receptors as the phytochemical compounds found in Cannabis. While endogenous cannabinoids include anandamide (AEA) and 2-arachidonoylglycerol (2-AG), exogenous phytocannabinoids include Δ-9 tetrahydrocannabinol (THC) and cannabidiol (CBD). These compounds finely tune neurotransmission following synapse activation, via retrograde signaling that activates cannabinoid receptor 1 (CB1R) and/or transient receptor potential cation channel subfamily V member 1 (TRPV1). Recently, the eCB system has been linked to several neurological diseases, such as neuro-ocular abnormalities, pain insensitivity, migraine, epilepsy, addiction and neurodevelopmental disorders. In the current study, we aim to: (i) highlight a potential link between the eCB system and neurological disorders, (ii) assess if THC exposure alters the expression of eCB-related genes, and (iii) identify evolutionary-conserved residues in CB1R or TRPV1 in light of their function.
Methods: To address this, we used several bioinformatic approaches, such as transcriptomic (Gene Expression Omnibus), protein–protein (STRING), phylogenic (BLASTP, MEGA) and structural (Phyre2, AutoDock, Vina, PyMol) analyzes.
Results: Using RNA sequencing datasets, we did not observe any dysregulation of eCB-related transcripts in major depressive disorders, bipolar disorder or schizophrenia in the anterior cingulate cortex, nucleus accumbens or dorsolateral striatum. Following in vivo THC exposure in adolescent mice, GPR55 was significantly upregulated in neurons from the ventral tegmental area, while other transcripts involved in the eCB system were not affected by THC exposure. Our results also suggest that THC likely induces neuroinflammation following in vitro application on mice microglia. Significant downregulation of TPRV1 occurred in the hippocampi of mice in which a model of temporal lobe epilepsy was induced, confirming previous observations. In addition, several transcriptomic dysregulations were observed in neurons of both epileptic mice and humans, which included transcripts involved in neuronal death. When scanning known interactions for transcripts involved in the eCB system (n = 12), we observed branching between the eCB system and neurophysiology, including proteins involved in the dopaminergic system. Our protein phylogenic analyzes revealed that CB1R forms a clade with CB2R, which is distinct from related paralogues such as sphingosine-1-phosphate, receptors, lysophosphatidic acid receptors and melanocortin receptors. As expected, several conserved residues were identified, which are crucial for CB1R receptor function. The anandamide-binding pocket seems to have appeared later in evolution. Similar results were observed for TRPV1, with conserved residues involved in receptor activation.
Conclusion: The current study found that GPR55 is upregulated in neurons following THC exposure, while TRPV1 is downregulated in temporal lobe epilepsy. Caution is advised when interpreting the present results, as we have employed secondary analyzes. Common ancestors for CB1R and TRPV1 diverged from jawless vertebrates during the late Ordovician, 450 million years ago. Conserved residues are identified, which mediate crucial receptor functions.
1 Introduction
The endocannabinoid (eCB) system acts as a regulator during neurotransmission. Indeed, following synaptic activation, eCBs will be released in the post-synapse, which, in turn, will decrease neuronal activation (Lu and Mackie, 2016). Historically, anandamide (AEA) and 2-arachidonoylglycerol (2-AG) were discovered first and named after the receptors they involve, mimicking the effects observed after “cannabis” exposure (Lu and Mackie, 2016). Expression of the cannabinoid receptor type 1 (CB1R) is mostly confined to the central nervous system (Matsuda et al., 1990), while CB2R expression is considered to be mainly located at the periphery (Munro et al., 1993), although potential clinical applications of targeting CB2R in neurological disorders have been highlighted recently (Kibret et al., 2022), likely due to microglial expression of CB2R (Tanaka et al., 2020). Other known eCB receptors include transient receptor potential vanilloid 1 (TRPV1) and G protein-coupled receptor 55 (GPR55) (Kano et al., 2009). Furthermore, GPR18 (McHugh et al., 2010), GPR119 (Overton et al., 2006) and some PPARs [peroxisome proliferator-activated receptors, reviewed in 2016 (O’Sullivan, 2016)] have also been identified as potential eCB receptors.
Derived from the omega-6 arachidonic acid (C20:4 n-6), both AEA and 2-AG are endogenous ligands of CB1R and induce decreased neuronal excitation (Ameri et al., 1999; Wilson and Nicoll, 2001; Straiker et al., 2023). Both molecules seem to involve a retrograde mechanism, whereby synthesis of AEA or 2-AG is located in the post-synapse, followed by diffusion in the synaptic cleft to then act on pre-synaptic vesicular fusion (Wilson and Nicoll, 2001). However, non-retrograde signaling has also been observed, as well as neuron-astrocyte communication (Castillo et al., 2012). The biosynthesis of AEA involves activation of N-acyltransferase (NAT) and N-acylphosphatidylethanolamine-hydrolyzing phospholipase D (NAPE-PLD), while 2-AG requires both phospholipase C (PLC) and diacylglycerol lipase (DAGL) (Kano et al., 2009). Degradation of AEA is performed by activation of fatty acid amide hydrolase (FAAH) and 2-AG is broken down by monoacyl-glycerol lipase (MAGL) (Kano et al., 2009). Interestingly, other degradation pathways can be triggered through two other enzymes: cyclooxygenase and lipoxygenase (Kano et al., 2009). The biological activity and regional distributions of these enzymes, at both the organ- and cell-specific levels, have been summarized previously (Basavarajappa, 2007).
Recent evidence has pinpointed that dysfunctions of the eCB system may cause neuro-ocular abnormalities (Bainbridge et al., 2022), pain insensitivity (Habib et al., 2019), migraine (Juhasz et al., 2009), epilepsy (Ludányi et al., 2008), addiction (Herman et al., 2006; Agrawal et al., 2009; Proudnikov et al., 2010), metabolic disease (de Miguel-Yanes et al., 2011), obesity (Russo et al., 2007) and neurodevelopmental disorders (Lu et al., 2008; Jung et al., 2012; Siniscalco et al., 2013; Gomis-González et al., 2016; Servadio et al., 2016; Wei et al., 2016; Barchel et al., 2018; Karhson et al., 2018; Melancia et al., 2018; Aran et al., 2019; Zamberletti et al., 2019; Wu et al., 2020). Many clinical observations have been paralleled and explained by animal studies, often linking AEA and neurodevelopmental disorders. In a recent meta-analysis, medical cannabis, largely used in the form of either tetrahydrocannabinol (THC) or cannabidiol (CBD), was found to be efficient for the treatment of several neurological diseases, such as chronic pain, spasticity, substance use disorder, epilepsy, Tourette’s syndrome and Parkinson’s Disease (Bilbao and Spanagel, 2022). Thus, activation of the eCB by natural or synthetic agonists provides efficient clinical benefits, which could be explained by the negative feedback role of the eCB system in neurotransmission.
The most probable evolutionary origin of the eCB system was explained in 2006 (McPartland et al., 2006). In fact, different evolutionary origins have been discovered, depending on which eCB is scrutinized. TRPV1 and GPR55 arose rather recently, likely in mammals (McPartland et al., 2006). The oldest member of the eCB family, FAAH, was present in early eukaryotic cells, while NAPE-PLD arose in opisthokonts, which comprise both animals and fungi. Recent evidence has corroborated previous observations, as two FAAH isoforms have been studied in the legume Medicago truncatula (Arias-Gaguancela et al., 2023), in addition to the known expression in the angiosperm Arabidopsis thaliana (Haq and Kilaru, 2020). Furthermore, a duplication event likely induced the appearance of CB1R and CB2R in vertebrates from an invertebrate deuterostome ancestor (Elphick et al., 2003), in which the pre-CBR protein was responsible for axonal regulation (Egertová and Elphick, 2007).
In this article, we aim to:
(i) Assess the potential transcriptomic alterations of 13 specific eCB-related transcripts in epilepsy, major depressive disorder, bipolar disorder or schizophrenia.
(ii) Assess the impact of THC on these 13 specific eCB-related transcripts.
(iii) Identify paralogues of CB1R or TRPV1 to infer evolutionary ancestry.
(iv) Identify conserved residues in CB1R or TRPV1 that mediate receptor function.
2 Methods
2.1 In silico analysis of publicly available RNA sequencing datasets
Previously-published RNA sequencing datasets were retrieved from the Gene Expression Omnibus (GEO) repository (Edgar et al., 2002; Clough and Barrett, 2016). Briefly, GEO datasets were scanned for adequate output using specific search terms. Table 1 describes the datasets used in the current study, with corresponding accession links given in Supplementary Table S1. We focused on 13 specific genes from the eCB system (Supplementary Table S2). Supplementary material M1 summarizes all DESeq2 results of the present study.
Table 1. Detailed description of all RNA sequencing datasets retrieved for the present study.
Differential gene expression analysis was performed in the freeware R Studio version 4.0.3 (R Core Team, 2013; RStudio Team, 2020). Following the workflow described by Sanchis et al. (2021), we first retrieved count data, using the GEOquery package. If necessary, count data were then normalized using the DESeq2 package. Differential expression analyzes were performed using the DESeq2 package. Gene annotations were performed with either the org.Hs.eg.bd (Homo sapiens), org.Mm.eg.db (Mus musculus) or AnnotationDbi packages, as appropriate. Volcano plots were drawn with the EnhancedVolcano package. Gene ontology (over-representation analysis, ORA) was performed with the GOplot package using Biological Processes. Gene set enrichment analysis (GSEA) was performed with clusterProfiler and enrichplot packages. Scatterplot graphs were drawn using the ggplot2 package. All of the above-mentioned packages are available under open access at the Comprehensive R Archive Network (CRAN) repository.1
2.2 Protein phylogenetics
Phylogenetic analyzes were performed in Mega X (Kumar et al., 2018) following alignment of FASTA sequences. Evolutional analysis was performed using the maximum likelihood method. Initial trees were obtained by applying Neighbor-Join (NJ) and BioNJ algorithms (Gascuel, 1997) to a matrix of pairwise distances estimated using the Jones-Taylor-Thornton (JTT)-based model (Jones et al., 1992). The final drawn trees are those with the highest log likelihoods. One hundred bootstraps (r = 100) were used to infer confidence in the output trees. BLASTP was used to scan different species with Homo sapiens as bait (reference sequence), available at the National Center for Biotechnology Information (NCBI).2
To search for conserved domains/residues, protein sequence (Homo sapiens) alignments was performed in Clustal Omega, available at the European Molecular Biology Laboratory – European Bioinformatics Institute (EMBL-EBI).3
2.3 Protein structural analyzes
Structural prediction of transmembrane protein topology (TMHMM, Transmembrane prediction using Hidden Markov Models) was performed from protein FASTA sequences using the platform provided by the Department of Health Technology at the Technical University of Denmark,4 which was made publicly available in 2001 (Krogh et al., 2001). Homology models were created using Phyre2 (Kelley and Sternberg, 2009; Kelley et al., 2015). For docking of anandamide into TRPV1, the Cryo-EM structure of human TRPV1 in complex with the analgesic drug SB-366791 (PDB: 8GFA) was used. The A chain of the multimer was used with ligand and phospholipids removed from the.pdb file. The anandamide ligand (CHEMBL15848) was retrieved from the ChEMBL webserver as a.csv file and converted to a.pdb file using the Online SMILES Translator and Structure File Generator and then to a.pdbqt file with merged non-polar hydrogens with AutoDockTools version 1.5.7. The ligand was docked into the receptor using AutoDock Vina 1.1.2 (Trott and Olson, 2010) with a docking grid 20 Å in length in the x, y, and z directions centered on the “tunnel” anandamide-binding pocket previously identified (Morales et al., 2022), flexible N438, Y555, Y554, Y487 and D708 side-chains and an exhaustiveness of 20. Structural analysis was performed by superposition in PyMol (Schrödinger and DeLano, 2020). A professional academic license was acquired for PyMol.
2.4 Protein–protein interaction networks
Protein–protein interaction networks (Szklarczyk et al., 2019) were predicted using the STRING (Search Tool for the Retrieval of Interacting Genes) database server,5 which contains both known and predicted protein interactions. Queries were pasted with either a single protein name or multiple names (as appropriate). When relevant, outer networks were searched for by expanding the nodes on display, using the built-in database tool, to reflect the immediate, intermediate and extended protein interactomes. For protein clustering within the networks, k mean clustering (based on centroids) was applied (Dalmaijer et al., 2022) directly in the server.
3 Results
3.1 Neurodevelopmental disorders and the eCB system: no change for 13 selected transcripts
We first examined dataset GSE80655 (Ramaker et al., 2017), which contains RNA sequencing of 281 postmortem samples extracted from patients with either major depressive disorder (MDD, n = 69), bipolar disorder (BP, n = 71), schizophrenia (SCZ, n = 71) or controls (n = 70), in an attempt to decipher the potential role of the eCB system in these pathologies (Table 1). Out of 13 transcripts of interest (CNR1, CNR2, GPR55, GPR110, TRPV1, MGLL, DAGLA, PLCB1, NAPEPLD, FAAH, ABHD6, PLA2G4A and ALOX8, Supplementary Table S2), none were considered significantly up- or down-regulated in MDD (Figure 1A), respectively defined as log2(fold change) >1 or < −1 and -log10(p) ≥ 1.3 (p = 0.05). Similar results were observed with BP (Figure 1B) and SCZ (Figure 1C). These results indicate that our 13 transcripts of interest are not dysregulated under these conditions.
Figure 1. The eCB system in postmortem tissues from patients with major depressive disorders (MDD), bipolar disorder (BD) or schizophrenia (SCZ). Volcano plots of the RNA sequencing dataset GSE80655 (Ramaker et al., 2017), performed on neurons in the anterior cingulate cortex, nucleus accumbens and dorsolateral striatum comparing control patients to patients with clinically-determined MDD (A), BD (B) or SCZ (C). Twelve transcripts of interest, involved in the eCB system, were extracted from the GSE80655 dataset. The horizontal and vertical bars represent the significant threshold, set at p = 0.05 (−log10(p) = 1.3) and log2(fold change) >1 or < −1, respectively.
3.2 Chronic THC induces significant up-regulation of GPR55 in the ventral tegmental area of mice
Next, we investigated whether we could detect transcriptional changes in our 13 transcripts of interest following acute or chronic THC exposure. To that end, we used GSE189821, GSE116813 and GSE70689 (Table 1). In GSE189821, 119 samples were analyzed from adolescent mice chronically exposed to THC (3 weeks at 10 mg/kg/day) or not (vehicle) and following a 14-day washout period, mimicking human behavior. In accordance with previous observations (Zuo et al., 2022), we found only minor transcriptomic dysregulations of our 13 eCB-related transcripts of interest following THC exposure in these mice (males and females pooled). This might be due to the long washout period that was applied. However, we noted 31 differentially expressed genes, with 13 transcripts down-regulated and 18 transcripts up-regulated, outside of our list of interest, when pooling all different brain regions together (Figure 2A). In the ventral tegmental area, GPR55 was found to be significantly upregulated after THC exposure, the only one from our list of interest and out of all brain regions examined. In GSE116813, 12 hippocampal samples were sequenced from mice exposed to chronic THC (21 days at 8 mg/kg/day) or not (Jouroukhin et al., 2019) and after 21 days of washout. After subsetting the dataset to only include wild-type mice, we observed 134 differentially expressed transcripts, 38 up-regulated and 96 down-regulated (Figure 2B). As observed with the previous dataset, none of these differentially-expressed transcripts were linked to our list of interest. The next dataset, GSE70689, was used to observe acute transcriptional changes in mice microglia exposed to in vitro THC (10 μM for 2 h) (Juknat et al., 2013). Here, we performed subsetting of microglia exposed to either THC or vehicle, excluding samples exposed to cannabidiol. Using such a method, 60 transcripts were found to be significantly differentially expressed, 2 up-regulated and 58 downregulated (Figure 2C). None of our transcripts of interest were found to be significantly dysregulated in such a dataset. A further analysis revealed that none of the previously differentially expressed genes (outside our 13 transcripts of interest) overlapped in the three datasets (Venn diagram in Figure 2D). In addition to these observations, we also ran ORA Gene Ontology (GO) analyzes for each dataset. GSE189821 and GSE116813 returned no GO enrichment, while GSE70689 revealed many GO enrichments. These included transcripts involved in immune response, such as T cell activation, lymphocyte migration, cytokine-mediated signaling and ERK (extracellular signal-regulated kinases) cascades (Figure 2D). These surprising results indicate that our 13 transcripts of interest are only minimally involved in the response to in vivo THC exposure, at least in mice. However, in vitro exposure to THC likely activates microglia, in an attempt to prevent THC-induced inflammation. Representations of the three datasets (GSE189821, GSE116813 and GSE70689) and our 13 transcripts of interest can be found on Figure 2E. Note that GPR55 is the only transcript of interest that was significantly upregulated after chronic in vivo THC exposure in mice, and only in the VTA (Figure 2E). This suggests that our 13 transcripts of interest were not affected by THC exposure, at least under these specific conditions, except for GPR55. GSEA analyzes did not reveal any biological processes of interest (not shown).
Figure 2. GPR55 is significantly upregulated in the ventral tegmental area of adolescent mice following chronic in vivo THC exposure. (A) Volcano plot comparing transcript expression in adolescent mice chronically treated (3 weeks) with vehicle to THC (10 mg/kg/day), extracted from the GSE189821 dataset (Zuo et al., 2022). (B) Volcano plot comparing transcript expression in adolescent mice chronically treated (3 weeks) with vehicle to THC (8 mg/kg/day), extracted from the GSE116813 dataset (Jouroukhin et al., 2019). (C) Volcano plot comparing transcript expression in mice microglia treated in vitro with THC (10 μM for 2 h) or control (untreated), extracted from the GSE70689 dataset (Juknat et al., 2013). Red dots indicate significantly up- or down-regulated transcripts. Horizontal bars represent the significant threshold, set at p = 0.05 (−log10(p) = 1.3). vertical bars represent the threshold for significant differences in transcript expression, with either down- (<−1) or up- (>1) regulated transcripts. (D) Histograms (left), Venn diagram (middle) and gene ontology (biological processes) plots (right) showing no overlap between the significantly differentially expressed transcripts for the three datasets. In both GSE186821 and GSE116813, no biological processes were returned. However, in GSE70689, gene ontology (GO) analysis revealed several biological processes that are affected by THC exposure. Of interest, transcripts involved in neuro-inflammation and neurophysiology. (E) Out of the 12 transcripts of interest involved in the eCB system, only GPR55 was significantly upregulated in the ventral tegmental area of adolescent mice exposed chronically to THC. Amyg, amygdala; DL Strm, dorsolateral striatum; NAcc, nucleus accumbens; PFC, prefrontal cortex; VTA, ventral tegmental area; GO, gene ontology; THC, Δ-9 tetrahydrocannabinol.
3.3 TPRV1 is significantly downregulated in the epileptic tissue of patients with temporal lobe epilepsy
Our next step was to investigate the eCB system with regard to temporal lobe epilepsy, in light of extensive evidence of the links between the two (Alger, 2004; Kwan Cheung et al., 2019; Bilbao and Spanagel, 2022; Martínez-Aguirre et al., 2022; Sugaya and Kano, 2022). Three datasets were thus investigated: GSE77578 (Srivastava et al., 2018), GSE190451 (Chen et al., 2023) and GSE74150 (Damasceno et al., 2020). In GSE77578, 56 hippocampal samples were compared from epileptic or control mice. Here, we only compared control mice to epileptic mice (temporal lobe epilepsy), excluding mice treated with an antiepileptic drug. We found 129 differentially expressed transcripts in such a model of temporal lobe epilepsy, with 43 transcripts up-regulated and 86 down-regulated (Figure 3A). In GSE190451, 6 neocortical samples were compared between patients with temporal lobe epilepsy (n = 3) and aged-matched non-epileptic controls (n = 3), an analysis that differs slightly from the published article, in which the authors analyzed a total of 8 samples (Chen et al., 2023). We found a total of 1873 differentially expressed transcripts in temporal lobe epilepsy, 833 up-regulated and 1,040 down-regulated (Figure 3B). In GSE74150, an animal model of genetic epilepsy was used (Wistar audiogenic rat). Here, we found 67 differentially expressed transcripts, 53 down-regulated and 14 up-regulated (Figure 3C). Up- and down-regulated transcripts in the three datasets presented some overlap (Figure 3D, Supplementary Table S3), showing that transcripts involved in, or responding to, epilepsy seem to be found across species and models. Several gene ontology enrichments (ORA) were detected in the first two datasets, ranging from regulation of neurophysiology to neuronal death. Interestingly, 37 biological process terms were found to be common across the two first datasets (Supplementary Table S4), such as synapse organization, neuron death, myelination and plasma membrane organization (Figure 3C). GSEA analyzes did not reveal any biological processes of interest (not shown). Out of our 13 transcripts of interest, TRPV1 was found to be significantly down-regulated in the neocortex of patients with temporal lobe epilepsy (Figure 3D). These results confirm, once again, that TRPV1 plays an important role in the pathophysiology of epilepsy, as observed previously (Jia et al., 2015; Saffarzadeh et al., 2015; Kong et al., 2019; Wang et al., 2019; Lazarini-Lopes et al., 2022).
Figure 3. TRPV1 is significantly downregulated in temporal lobe epilepsy. (A) Volcano plot comparing transcript expression in control mice to mice with temporal lobe epilepsy (pilocarpine model), extracted from the GSE77578 dataset (Srivastava et al., 2018). (B) Volcano plot comparing transcript expression in tissue samples from control patients to fresh tissue samples extracted from the epileptogenic foci from patients with temporal lobe epilepsy [GSE190451 dataset (Chen et al., 2023)]. (C) Volcano plot comparing transcript expression in control rats to a genetic model of epilepsy [Wistar audiogenic rat, Damasceno et al., 2020]. Red dots indicate significantly up- or down-regulated transcripts. Horizontal bars represent the significant threshold, set at p = 0.05 (−log10(p) = 1.3). Vertical bars represent the threshold for significant differences in transcript expression, with either down- (<−1) or up- (>1) regulated transcripts. (D) Histograms (left), Venn diagram (middle) and gene ontology plots (right) showing small overlap between the significantly up- or down-regulated transcripts in the three datasets. On gene ontology plots, several biological processes are found dysregulated in both datasets, including neurophysiology and neuronal death. TRPV1 (bottom, blue arrow) is significantly down-regulated in epilepsy (human, GSE190451). GO, gene ontology; TLE, temporal lobe epilepsy.
3.4 The eCB system interactome
Using the STRING database (Szklarczyk et al., 2019), containing both known and predicted protein interactions, we first fed into the database single protein names (from our list of 13 transcripts/proteins). Analysis of the immediate interactomes revealed great disparities in terms of protein overlap in pair-wise comparisons (Figure 4A). Interestingly, great overlap was observed with CB1R, CB2R and ADGRF1 (GPR110). We also noted great protein overlap in the interactomes of DAGLA, NAPE-PLD and ABHD6. As we expanded the interactome searches, similar conclusions were reached (Figure 4A). Surprisingly, TRPV1 and ALOX15B presented the lowest overlap in all levels, whether immediate, intermediated or extended interactomes (Figure 4A). We next scanned the database using our 13 proteins of interest by bulk feeding all 13 proteins into the search equation. When examining the immediate interactome of these 13 proteins, only the same 13 proteins were returned (Figure 4B), providing no additional information of interest. However, when the interactome search was expanded to an intermediate level, new proteins appeared, including GNAI1, GNAI2, GNAQ, and TRMP8 (Figure 4B). When the search was expanded to another level (extended interactome), new protein interactions arose, such as CALM3, CALML3, CALML5, CALML6 and TRPA1 (Figure 4B). Our final step was to apply k means clustering using a fixed number of clusters (k = 3). Very surprisingly, 4 proteins clustered outside of the main cluster containing the proteins of interest: TRPV1, PLC1B, ALOX15B and PLA2G4A (Figure 4C). We also noted that TRPM8 and GATB clustered with the majority of our proteins of interest, via FAAH, acting as a link within the extended interactome (Figure 4C).
Figure 4. Endocannabinoid protein interactomes. (A) Pair-wise matrix comparisons of immediate, intermediate or extended protein interactomes overlap, using single proteins as input into the STRING database, from our list of 13 proteins of interest (based on the 13 transcripts of interest). Note the high interactome overlap with CB1R, CB2R, GPR55, ADGFR1 (GPR110) on one side and NAPE-PLD, DAGLA and ABHD6 on the other. The interactome of ALOX15B appears to not overlap with the other 12 proteins of interest, except PLA2G4A (14.3–18.2%). (B) Protein interactomes at an immediate, intermediate or extended level. (C) K means clustering (k = 3) was applied to identify segregation groups within the extended endocannabinoid protein interactome.
3.5 Paralogues of human CBRs
Our next step was to try to determine how CB1R (the protein produced from the CNR1 gene) evolved within the human proteome. Using GeneCards (Stelzer et al., 2016) with Homo sapiens CB1R as input, 18 paralogues were returned. As expected, the closest paralogue was CNR2 (Figure 5A). CNR1 and CNR2 formed a distinct clade, with the next most closely related G-coupled protein receptors (GPCRs) being lysophosphatidic acid receptors (LPARs) and sphingosine-1-phosphate receptors (S1PRs). The furthest clade included melanocortin receptors 1–4 (MCRs) and GPR119. Inspection of the sequences revealed three conserved regions (Figure 5B), which could be divided into residues conserved within the clade containing CNR, LPAR and S1PR families (Figure 5B, blue) and residues only conserved within the CNR clade (Figure 5B, orange). Three conserved motifs (CMs) were identified: CM1, CM2 and CM3 are located within transmembrane (TM) 2, 3 and 7, respectively (Figure 5C). We then examined the location and function of these residues in the Homo sapiens CB1R structure bound to an anandamide analog (Krishna Kumar et al., 2023) and in the lysophosphatidic acid receptor structure bound to its native lysophosphatidic acid ligand (PDB: 7TD1) (Liu et al., 2022). The residues conserved within all three CNR, LPAR and S1PR families are largely found in the lower part of the receptor and include five of the six residues identified as being important in receptor activation (Figure 5D), as observed before (Hua et al., 2020). The ligand binding pocket in both receptors is formed by a large pocket of hydrophobic residues. Most of the residues conserved only within the CNR clade are located in this area and include some hydrophobic residues and H178 which forms a hydrogen bond with the terminal hydroxyl of the ligand (Figures 5D,E).
Figure 5. Paralogues of CNR1 in Homo sapiens and evolution of CB1R structure. (A) Phylogenic relationships of cannabinoid receptor 1 (CNR1) and the eighteen CNR1 paralogues listed in Genecards. The evolutionary history was inferred by using the Maximum Likelihood method and Jones-Taylor-Thornton matrix-based model (Jones et al., 1992). The tree with the highest log likelihood (−13679.67) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches (orange, from 100 bootstraps). Initial trees for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the JTT model, and then selecting the topology with superior log likelihood value. A discrete Gamma distribution was used to model evolutionary rate differences among sites [3 categories (+G, parameter = 1.3423)]. This analysis involved 19 amino acid sequences. There were a total of 504 positions in the final dataset. Blue values indicate branches length. (B) Clustal Omega alignment of Homo sapiens (Hs) CNR1 and CNR2. Runs (≥ 3 amino acids or longer runs with >80% identity) that are present in CNR1/2 but absent in other Hs GPCR paralogues are highlighted in blue with unique CNR amino acids (compared to S1PR and LPAR) in orange on CNR1 sequence. Three areas of interest are identified, conserved motif (CM) 1 (residues 151–181 in CNR1), CM2 (200–216) and CM3 (379–401). (C) Transmembrane Helices Hidden Markov Models scan of CNR1. CM1 is located in TM (transmembrane) 2, CM2 in TM3 and CM3 in TM7. (D) Homo sapiens cannabinoid receptor type 1 (PDB: 8GHV) bound to anandamide analog AMG315 (dark blue). Residues conserved in CNR, LPAR and S1PR clades are colored blue and those conserved in only the CNR clade are colored orange. Residues critical for receptor activation are shown as cyan sticks and additional, non-conserved residues in the anandamide binding pocket are colored yellow. Axes colored X red, Y green, Z blue. (E) Comparison of Homo sapiens cannabinoid receptor type 1 (left) and Homo sapiens lysophosphatidic acid receptor (right, PDB: 7TD1) binding pockets. Axes colored X red, Y green, Z blue. Color scheme as in (D) with lysophosphatidic acid show in magenta.
Out of 13 deuterostome species spanning from mammals to echinoderms (plus insects as protostomes), CB1R orthologues were found in placental mammals, marsupials, monotremes, birds, reptiles, amphibians, teleosts, elasmobranchs and Agnatha (Supplementary Table S5), but absent in cephalochordates, tunicates, echinoderms and insects. Furthermore, CB1R orthologues were also not found in plants, bacteria, Cnidaria and Porifera, except in Klebsiella (gram-negative bacteria) which presented 100% homology with the bird CB1R sequence (not shown), thus likely accounting for a database error in GenBank. When the same experiment was conducted with CB2R, we observed the presence of orthologues in placental mammals, marsupials, monotremes, birds, reptiles, amphibians and teleosts, the latter presenting two orthologues (Supplementary Table S5). Within the BLASTP database, unknown CBRs were found in cephalochordates and tunicates. These results suggest that the classical endocannabinoid receptors (CB1R and CB2R) are confined to a subgroup of deuterostomes, containing urochordates, cephalochordates and chordates, but absent from other deuterostomes, other animals, plants and prokaryotes, thus partly corroborating previous observations (McPartland et al., 2006; Elphick, 2012). Furthermore, it is suggested that a duplication event has occurred in teleosts since they diverged from other chordates. A further phylogenetic analysis revealed a clear division between CB1R and CB2R into two distinct clades (Figure 6A). We noted that most taxa presented copies of both proteins. There are two explanations for the lack of CNR2 paralogues in agnathans and elasmobranchs: (i) a duplication may have occurred since the divergence of other chordates and elasmobranchs or (ii) CNR2 paralogues were lost by both agnathans and elasmobranchs, the first scenario being by far the more likely explanation. Finally, we noted that the CBRs of Cephalochordates and Tunicates are very different, making it difficult to distinguish these two proteins from other related GPCRs. To address this, we performed multiple alignments with CB1R and CB2R protein sequences from the above-mentioned species. Our first alignment included cephalochordates and tunicates and led to great disruptions of the CMs (Figure 6B). The second alignment, which was performed without cephalochordates and tunicates, produced more defined, but yet not optimal, CMs (not shown). Our final alignment was performed without cephalochordates, tunicates and Danio rerio. This final approach produced robust identification of the three CMs (Figure 6C), thus corroborating the above results. Subsequently, we created homology models of the agnathan CB1R and cephalochordate unclassified CBR. Both the agnathan (Petromyzon marinus: ModelArchive ma-yvfgy) and cephalochordate (Branchiostoma floridae: ModelArchive ma-mi5fb) receptors contain the six residues important for receptor activation. However, the cephalochordate receptor lacks the majority of the residues involved in the AEA-binding pocket, including a histidine equivalent (H178) needed for hydrogen-bond formation with the terminal hydroxyl (Figures 6D,E).
Figure 6. Evolutionary analysis of CNR1/CB1R in different taxa and structural comparisons. (A) Maximum Likelihood analysis using a neighbor joining method. Sequences aligned with Clustal W and phylogenetic analysis carried out using JTT + G with 100 bootstraps (orange values). Sphingosine-1-phosphate receptor 1 (S1PR1) was used to root the tree (blue branch). (B) Alignment of all CNR proteins. The three motifs identified earlier are highlighted in blue. Sequence not containing the conserved motifs were removed from the figure. (C) Alignment of CNR proteins identified in chordates (with the removal of Teleost CNR2). Sequence not containing the motifs removed from the figure. Blue highlighting are residues that were identified earlier as conserved (from CM1, CM2 and CM3). (D) Homology models of Petromyzon marinus CB1R superposed onto the anandamide analog AMG315 (dark blue) (left) and Branchiostoma floridae cannabinoid receptor of unknown type (right). Residues conserved in CNR, LPAR and S1PR clades are colored blue and those conserved in only the CNR clade are colored orange. Residues critical for receptor activation are shown as cyan sticks and additional, non-conserved residues in the anandamide binding pocket are colored yellow. Axes colored X red, Y green, Z blue. (E) Petromyzon marinus CB1R superposed to the anandamide analog binding pocket. Axes colored X red, Y green, Z blue. Color scheme as in (D).
3.6 Paralogues of human TRPVs
Our final analyzes focused on TRPV1, as this is the second eCB receptor in the central nervous system (together with CB1R). Using similar approaches, we first determined the phylogeny of the six Homo sapiens paralogues, using TRPV1 as the reference sequence. We observed that TRPV1 and TRPV2 diverged recently, which arose from a duplication event (Figure 7A). Likewise, divergence of TPRV3, TRPV5 and TRPV6 arose earlier, with TPRV5 and TRPV6 being the most recent split (Figure 7A). Two small CMs were noted throughout those 6 different proteins, including Y375-G376-P377, L382-Y383-D384-L385, L547-G548-W549 and R557-G558 (Figure 7B). Using Homo sapiens TRPV1 sequence as bait, we found that TPRV5 or TRPV6 were present across the most diverse range of taxa (Supplementary Table S6). These receptors are present in Bilatera but absent from Porifera, the latter presenting weak protein homology. These observations suggest that the TRPV family is specific to the Bilatera. Furthermore, TPRV1 and TRPV2 resulted from later duplication events in the common ancestor of jawed vertebrates and amniotes, respectively (Figure 7C). In addition, our phylogenetic analysis highlighted that eCB-binding receptors (TRPV1 and TRPV2) form one clade, which includes taxa from elasmobranchs onwards, with a more recent duplication forming TRPV2, found in amniotes only (Figure 7C), as suggested above (Supplementary Table S6). Alignments of TPRV1 and TRPV2 from these species correctly identified conserved residues involved in channel gating and selectivity, as well as ligand binding [anandamide, vanilloid, double-knot toxin (DkTx), Figure 7D]. Conserved residues were often found throughout species, indicative of the importance of these residues for TRPVs function. Finally, out of these 28 functional residues in TRPV1, only 8 were conserved in TRPV5, whether fully conserved or present with conservative substitutions (Figure 7E). Taken together, these results pinpoint a divergence of eCB receptors between jawed and jawless vertebrates, around 450 million years ago during the late Ordovician. These observations are consistent with previous studies (Elphick et al., 2003; McPartland, 2004; McPartland et al., 2006).
Figure 7. Phylogeny of TRPVs and conserved residues. (A) Simple phylogeny illustrating the relationships of the six TRPVs human paralogues. (B) Clustal alignment of the six Homo sapiens TRPVs. Two small fully conserved motif highlighted in blue. (C) Minimum Evolution Tree of Chordate TRPV family. Aligned using ClustalW, and ME tree (JTT + G) with 500 bootstraps. DmInactive (proteostome TRPV protein) is used to root the tree (blue branch). (D) Alignment of TRPV1 and TRPV2 showing conservation of sites. In each case, highlighted residues are either fully conserved or with conservative substitutions. (E) Comparison of TRPV1 and TRPV5, aligned using ClustalW.
We examined the structures of TRPV1 and TRPV4 to identify residues important for ligand binding. TRPV1 can be activated by endogenous ligands such as AEA and exogenous ligands such as capsaicin, the component of chili peppers which creates a sensation of burning (Lam et al., 2005). Structural studies have identified the capsaicin binding site, known as the vanilloid-binding pocket (Nadezhdin et al., 2021; Kwon et al., 2022). There are no experimentally determined structures of the TRPV1 receptor bound to AEA. However, using a model of human TRPV1 based on the experimental rat structure, molecular dynamics and docking studies have identified two possible binding sites for anandamide (Muller et al., 2020; Morales et al., 2022), represented in Figure 8A. One site was the vanilloid-binding pocket (VBP) but the preferred site was a tunnel (AEA tunnel) involving interactions with Y554, Y555, Y487, D708, and N438. We docked AEA into the AEA tunnel using the recently experimentally determined structure of human TRPV1 (PDB: 8GFA) (Neuberger et al., 2023). One of nine docking solutions had a similar orientation to that shown in Morales et alia (Morales et al., 2022), with Y555, N438 and Y554 forming hydrogen bonds with the amide group hydrogen and the terminal hydroxyl group oxygen (Figure 8B, ModelArchive ma-dol4y). The five interacting residues are well conserved across species in TRPV1 and TRPV2. AEA itself does not directly bind to TPRV4, however this receptor is activated by an AEA metabolite 5′,6′-epoxyeicosatrienoic acid (5′,6′-EET) (Watanabe et al., 2003). Homo sapiens TRPV4 (PDB: 8TD1) (Nadezhdin et al., 2021) is missing equivalents to TRPV1 Y847 and Y555 which stabilize AEA with hydrogen bonds in the AEA tunnel. Simulations have suggested that residues K535, F549, Q550, Y591 and R594 of TRPV4 form the 5′,6′-EET binding pocket (Berna-Erro et al., 2017), as displayed in Figure 8C. Of these, Y591 is equivalent to Y554 in TRPV1 which forms a hydrogen bond with the terminal hydroxyl of AEA, but the remaining pocket is very different. We also examined differences in the vanilloid-binding pocket by superimposing capsaicin from the rat (Rattus norvegicus) TRPV2 complex structure (PDB: 8SLY) (Gochman et al., 2023) onto the Homo sapiens TRPV1 receptor (PDB: 8GFA). Whilst some important hydrophobic residues are conserved in the equivalent pocket in TRPV4, there are no equivalent residues for the hydrogen-bonding Y511 and T550 which were previously observed to be involved in propagating ligand-activated conformational changes (Figure 8D) (Kwon et al., 2022).
Figure 8. Structural comparisons of Homo sapiens TRPV1 and TRPV4. (A) Homo sapiens TRPV1 (PDB: 8GFA) showing one of four monomers in detail. Anandamide (AEA) tunnel residues shown in dark blue and vanilloid-binding pocket residues shown in pink. (B) Comparison of residues in the AEA tunnel for TRVP1 (white sticks with residue labels) and TRPV4 (green sticks, PDB: 8TD1). AEA is shown in blue sticks. (C) Comparison of residues in the AEA tunnel for TRVP1 (white sticks) and the 5′,6′-EET pocket for TRPV4 (green sticks with residue labels). AEA is shown in blue sticks. (D) Comparison of residues in the vanilloid-binding pocket for TRVP1 (white sticks with residue labels) and TRPV4 (green sticks). Capsaicin (natural phytochemical ligand) is shown in magenta.
4 Discussion
The location of CB1Rs is restricted to pre-synaptic terminals of both excitatory and inhibitory neurons (Mátyás et al., 2008). The physiological roles of CB1R are to mediate mood (Sbarski and Akirav, 2018), memory (Borgan et al., 2019a), stress (Arévalo et al., 2001; Barna et al., 2004) and locomotion (Han et al., 2023). These functions were reviewed in 2020 (Cristino et al., 2020). Activation of CB1Rs induces synaptic plasticity via long-term depression (Wu et al., 2015). Indeed, following synaptic activation, endogenous endocannabinoids (AEA and 2-AG) will be released from the post-synaptic element to then diffuse to the pre-synaptic element, where both molecules will act as agonists at CB1 receptors. Such activation of CB1 receptors will reduce neurotransmitter release, ending the cycle of synapse activation. This mechanism is often referred to as endocannabinoid-mediated negative feedback. Interestingly, CB1Rs were found to be expressed in the nucleus accumbens of vervet monkeys (Kucera et al., 2018), thus highlighting the role of these receptors in the reward circuitry (Méndez-Díaz et al., 2019; Braunscheidel et al., 2022; Ceccarini et al., 2022; Han et al., 2023). Dysfunction of CB1Rs has been linked to a wide array of neurological disorders such as Parkinson’s Disease, Alzheimer’s Disease, Huntington Disease and multiple sclerosis [reviewed in Cristino et al. (2020)].
In contrast to CB1R, TRPV1 is expressed in both presynaptic (Tominaga et al., 1998; Medvedeva et al., 2008) and postsynaptic (Bae et al., 2004; Grueter et al., 2010) elements. The physiological roles of TRPV1 intertwined those of CB1R, such as memory and locomotion, but also other functions, such as food intake (Cristino et al., 2020). Extensive evidence suggests that TRPV1 can be linked to epilepsy (Jia et al., 2015; Saffarzadeh et al., 2015; Kong et al., 2019; Wang et al., 2019; Lazarini-Lopes et al., 2022), although causal determination remains unclear. Indeed, in naïve rodents, chronic exposure to the TRPV1 agonist capsaicin induces spontaneous seizures (Jia et al., 2015). On the other hand, pre-treatment with the TRPV1 antagonist capsazepine decreased seizure severity and mortality in pentylenetetrazol-induced epilepsy (GABAA antagonist) (Jia et al., 2015). Another study has observed altered electrophysiological properties in microglial cells following capsaicin, which was hypothesized as driving hyperthermia-induced seizures (Kong et al., 2019), although direct experimental confirmation has not been performed. Increased expression of TRPV1 was also observed after the successful induction of temporal lope epilepsy through kindling (audiogenic epilepsy) in rats (Lazarini-Lopes et al., 2022). Similar results were observed in another study using a model based on pilocarpine exposure (Saffarzadeh et al., 2015). Another elegant confirmation of the link between epilepsy and TRPV1 was observed in animals lacking TRPV1 (knock-out). Indeed, compared to wild-type mice, mice lacking TRPV1 presented decreased seizure latency, duration and severity (Wang et al., 2019). In agreement with these observations, our current transcriptomic results report a downregulation of TPRV1 in temporal lobe epilepsy. Indeed, we found a significant downregulation of TRPV1 mRNA in the neocortex of patients with temporal lobe epilepsy compared to control tissues. As mentioned above, preclinical studies have shown that agonism of TRPV1, using capsaicin, can induce tonic–clonic seizures, while TPRV1 antagonism reduced mortality in animals in which seizures were induced (Jia et al., 2015). These results were also observed in knock-out animals for TRPV1 (Jia et al., 2015). However, contrasting results were observed in a 2020 study, in which capsaicin inhibited epileptiform activities in the prefrontal cortex (Pasierski and Szulczyk, 2020). The same authors also witnessed decreased neuronal excitation, thus explaining the above-mentioned effects (Pasierski and Szulczyk, 2020), which are mediated by decreased sodium currents. These were not isolated results, as other studies have also concluded that capsaicin possesses anti-convulsive properties in rats (Abdel-Salam et al., 2020) and mice (Lee et al., 2011; Pezzoli et al., 2014). Interestingly, hippocampal expression of CB1Rs was also found to be significantly decreased in patients suffering from epilepsy (Ludányi et al., 2008). Thus, in light of these previous results and those from the present study, TRPV1 could be of interest for developing new therapeutics for patients with refractory epilepsy.
The role of the eCB system in neurodevelopmental disorders has been described previously. Indeed, many studies observed a direct relationship between altered CB1R availability in schizophrenia (Dean et al., 2001; Zavitsanou et al., 2004; Newell et al., 2006; Eggan et al., 2008; Urigüen et al., 2009; Wong et al., 2010; Dalton et al., 2011; Jenko et al., 2012; Ceccarini et al., 2013; Ranganathan et al., 2016; Muguruza et al., 2019; Bloch Priel et al., 2023) or during psychosis (Borgan et al., 2019b, 2021), which is also corroborated by animal models of schizophrenia (Szűcs et al., 2016; Stark et al., 2022). Furthermore, a link between THC exposure and schizophrenia has been evidenced in humans (Hjorthøj et al., 2021) and animal models (Rodríguez et al., 2017). Compared to control patients, levels of AEA are increased in the blood (De Marchi et al., 2003; Ibarra-Lecue et al., 2022) and cerebrospinal fluid (Giuffrida et al., 2004; Leweke et al., 2007; Minichino et al., 2019) in patients with schizophrenia. In addition, exposure to THC alters AEA levels (Leweke et al., 2007; Seillier et al., 2020; Ibarra-Lecue et al., 2022). Furthermore, lower brain FAAH availability was observed in adolescent chronic cannabis users (Jacobson et al., 2021). Targeting FAAH has proven an efficient treatment to decrease cannabis use and withdrawal (D’Souza et al., 2019). During psychosis, increased levels of DAGL and NAPE are reported in mononuclear cells (peripheral blood) compared to control patients (Bioque et al., 2013), two enzymes controlling endocannabinoid production. In addition, such a study also reported increased FAAH levels during psychosis (Bioque et al., 2013). Altogether, these studies provide ample evidence for a link between the eCB system and neurodevelopmental disorders. Recent studies reviewed the relationship of the eCB in several neurological pathologies such as depression (Gallego-Landin et al., 2021), schizophrenia (Garani et al., 2021), autism (Su et al., 2021) and substance use disorder (Navarro et al., 2022). In the current study, our transcriptomic analysis found that GPR55 is upregulated in the ventral tegmental area of adolescent mice following THC exposure, even after a long washout period (2 weeks). These results can be explained by the fact that THC possesses psychoactive effects (Ameri et al., 1999) and that the ventral tegmental area is known to respond to THC (French et al., 1997; Lupica et al., 2004; Ostlund et al., 2022). However, what remains surprising is the fact that GPR55, and not CB1R, is upregulated after THC exposure, in light of two studies demonstrating loss of CB1R-dependent plasticity after chronic THC exposure (Friend et al., 2017; Ostlund et al., 2022). The answer appears to lie within the capacity of GPR55 to respond to THC (Ryberg et al., 2007; Lauckner et al., 2008), which has been the reason why GPR55 is classed as an endocannabinoid receptor. In fact, THC, AEA, methanandamide (stable metabolic analog of AEA) and JWH015 (cannabinoid agonist) all seem to elicit intracellular responses in GPR55-transfected HEK cells (Lauckner et al., 2008), thus proving a link between GPR55 and THC. Furthermore, CBD is also known to abolish GPR55-mediated synaptic events (Sylantyev et al., 2013), highlighting again a strong interaction between GPR55 and the eCB system.
The two phytocannabinoids THC and CBD are known to bind to CB1R (Jakowiecki and Filipek, 2016; Hua et al., 2017; Chung et al., 2019), since the only chemical difference between these two compounds is a bond disconnection in CBD (Shao et al., 2016). Functionally, CBD has a low affinity for the CB1R (McPartland et al., 2007) and is described as an antagonist at CB1R (Thomas et al., 2007). However, evidence suggests that CBD rather behaves as a negative allosteric modulator of THC (Laprairie et al., 2015; Chung et al., 2019). Interestingly, CBD also binds to TRPV1-4 channels, as observed recently (Costa et al., 2004; Qin et al., 2008; Iannotti et al., 2014; Muller et al., 2019; Starkus et al., 2019; Galaj et al., 2020; Etemad et al., 2022; Landucci et al., 2022). These observations were also made outside of the central nervous system (De Petrocellis et al., 2012; Nabissi et al., 2013; de la Harpe et al., 2022).
We also observed that the amino acid sequence of TRPV1 in Danio rerio is missing Y554 and Y555, which have been shown to allow AEA binding (Muller et al., 2020). Instead, TRPV1 in zebrafish possesses A554 and A555. The switch from polar amino acids (Tyr-Tyr) to non-polar amino acids (Ala-Ala) in the TRPV1 orthologue of zebrafish does not seem to impact heat detection, which is the main function of TRPV1 (Zheng and Wen, 2019). Despite these amino acids substitutions, in vitro and in vivo observations confirmed that TRPV1 is triggered above 25°C in zebrafish (Gau et al., 2013). In contrast, heat activation of TRPV1 in mammals occurs above 42°C (Caterina et al., 1997), which could be explained by slightly different structural conformations following amino acid substitutions. Interestingly, TRPV3 seems to possess opposite heat-detection functions between humans and reptiles (Xenopus tropicalis), as observed before (Saito et al., 2011), with the receptor activated by cold temperatures (≤ 16°C). Behaviorally, TRPA1 seems to be the receptor mediating avoidance of dangerous temperatures in frogs (Saito et al., 2022). Similar to our findings, Saito et alia also reported that, in vertebrates, TRPV channels form very distinct clades (Saito et al., 2011). Furthermore, TRPV1 and TRPA1 seem to have been co-expressed in an early vertebrate common ancestor (Saito et al., 2012).
The link between TRPV1 and heat can be experienced when eating hot chili peppers, which dose-dependently trigger receptor activation (Mohapatra and Nau, 2003; Voets et al., 2004; Neelands et al., 2005; Vetter et al., 2006; Vasquez et al., 2020; Matsuyama et al., 2021). THC can bind TRPV2-4, TRPA1 and TRPM8, but not TRPV1 (Muller et al., 2019). This seems very different from CBD, which has been shown to bind to TRPV1 (see above). What remains fascinating is the fact that TRPV channels are involved in sensory perception (chemical nociception, heat detection) in the peripheral nervous system, whilst TRPV1 in the central nervous system is involved in synaptic plasticity, neurite growth and signaling (Moran et al., 2004).
In terms of evolution, our results suggest that the eCB system appeared during the late Ordovician. McPartland and colleagues studied the evolution of all components of the eCB system in a study published in 2006 (McPartland et al., 2006). In their article, authors evidenced that TRPV1 and GPR55 are limited to mammals, whilst CB2R and DAGL were found amongst vertebrates (McPartland et al., 2006). Other components, such as MAGL, are found in chordates, while all eukaryotes express FAAH (McPartland et al., 2006). Thus, it can be safely hypothesized that the origins of the eCB appeared in a common eukaryotic ancestor. Such a study also suggested that ligand-metabolizing enzymes of the eCB system evolved prior to functional cannabinoid receptors (McPartland et al., 2006). In the absence of these receptors, such as CB1R, CB2R and TRPV1, endocannabinoids, acting as signaling molecules, such as AEA and 2-AG, might possess other functional properties (McPartland et al., 2006). A likely explanation is modulating physiological excitability in the nervous system (McPartland, 2004). 2-AG is expressed in insects and, apart from signaling properties, 2-AG might serve as a deterrent to predators, such as vertebrates, who do express functional receptors to 2-AG, such as TRPV1. This mechanism was likely driven by convergent evolution. Interestingly, this feeding-deterrent mechanism was also suggested for plants (McPartland et al., 2006), expressing endocannabinoid-like compounds but without eCB receptors. Intertwined evolution of receptors, enzymes and ligands of the eCB system was likely driven by convergence, divergence and parallel evolution (McPartland and Guy, 2004). To the best of our knowledge, pinpointing an exact timeframe for such mechanisms is currently impossible. However, appearance of the eCB in the late Ordovician seems to be the accepted consensus, and our current study further corroborates these previous suggestions.
Indeed, recent studies have observed brain-like structures in fossils from the Cambrian-Ordovician eras (Dong et al., 2022; Ortega-Hernández et al., 2022; Pates et al., 2022), which seem to support the idea of an earlier origin of the nervous system, likely occurring before the appearance of annelids (Parry and Caron, 2019). This pinpoints to the nervous system originating just before 500 million years ago, with the eCB system likely appearing 50 million years after the former. Further evidence has been gathered in invertebrates (annelids and insects), in which observations of eCB-like receptors were discovered (Stefano et al., 1997; Elphick, 1998, 2012; McPartland et al., 2006). Some components of the eCB system have been found in species much older than annelids. Indeed, FAAH, which is responsible for metabolism of AEA, has been identified in the moss Physcomitrella patens, with a high catabolic activity for AEA (Haq and Kilaru, 2020). Furthermore, 9 paralogues were identified in this species, with some closely related to rat FAAH, while others closely related to plant FAAH (Haq and Kilaru, 2020). It is now accepted that the nervous system appeared in bilaterians ancestors. Recovered fossils suggest that the nervous system was already present 540 million years ago (McPartland and Guy, 2004; Budd, 2013; Heger et al., 2020). Furthermore, several studies confirmed the presence of a nervous system in bilaterians (Adey, 1951; Holland et al., 2013; Yang et al., 2013, 2016; Martín-Durán et al., 2018), sponges (Leys, 2015; Musser et al., 2021; Yañez-Guerra et al., 2022) or ctenophores (Jékely et al., 2015; Moroz, 2015; Burkhardt et al., 2023). Interestingly, the oldest bilaterians fossils recovered so far range between 540 and 585 million years (Martin et al., 2000; Pecoits et al., 2012; Chen et al., 2019; Evans et al., 2020). Altogether, it can be safely assumed that the nervous system appeared before such a timeframe.
To conclude, we provide here new evidence that GPR55 is upregulated in the ventral tegmental area of mice following THC exposure and that TPRV1 is downregulated in the epileptogenic focus of patients with temporal lobe epilepsy. These results further support previous observations linking the eCB system to epilepsy, cannabis and neurodevelopmental disorders. We also report here an updated analysis of the evolution of CB1R and TRPV1, two major actors in the eCB system.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
AC: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing. MW: Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing. KM: Formal analysis, Investigation, Writing – review & editing. MDM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Software, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the University of Worcester, grants SSE2206 (to MDM) and ROFF22234MDM (to MDM).
Acknowledgments
The authors would like to thank the University of Worcester for support with funding. In addition, the authors would like to thank Pr. Eleanor Bradley, Dr. Derek McDougall, and Dr. Steven Coles.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fninf.2024.1294939/full#supplementary-material
Footnotes
1. ^ https://cran.r-project.org
2. ^ https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Proteins
3. ^ https://www.ebi.ac.uk/Tools/msa/clustalo/
References
Abdel-Salam, O. M. E., Sleem, A. A., Sayed, M. A. E. B. M., Youness, E. R., and Shaffie, N. (2020). Capsaicin exerts anti-convulsant and neuroprotective effects in Pentylenetetrazole-induced seizures. Neurochem. Res. 45, 1045–1061. doi: 10.1007/s11064-020-02979-3
Adey, W. R. (1951). The nervous system of the earthworm Megascolex. J. Comp. Neurol. 94, 57–103. doi: 10.1002/cne.900940104
Agrawal, A., Wetherill, L., Dick, D. M., Xuei, X., Hinrichs, A., Hesselbrock, V., et al. (2009). Evidence for association between polymorphisms in the cannabinoid receptor 1 (CNR1) gene and cannabis dependence. Am. J. Med. Genet. B Neuropsychiatr. Genet. 150B, 736–740. doi: 10.1002/ajmg.b.30881
Alger, B. E. (2004). Endocannabinoids and their implications for epilepsy. Epilepsy Curr. 4, 169–173. doi: 10.1111/j.1535-7597.2004.04501.x
Ameri, A., Wilhelm, A., and Simmet, T. (1999). Effects of the endogeneous cannabinoid, anandamide, on neuronal activity in rat hippocampal slices. Br. J. Pharmacol. 126, 1831–1839. doi: 10.1038/sj.bjp.0702478
Aran, A., Eylon, M., Harel, M., Polianski, L., Nemirovski, A., Tepper, S., et al. (2019). Lower circulating endocannabinoid levels in children with autism spectrum disorder. Mol. Autism. 10:2. doi: 10.1186/s13229-019-0256-6
Arévalo, C., de Miguel, R., and Hernández-Tristán, R. (2001). Cannabinoid effects on anxiety-related behaviours and hypothalamic neurotransmitters. Pharmacol. Biochem. Behav. 70, 123–131. doi: 10.1016/s0091-3057(01)00578-0
Arias-Gaguancela, O., Herrell, E., Aziz, M., and Chapman, K. D. (2023). Two legume fatty acid amide hydrolase isoforms with distinct preferences for microbial- and plant-derived acylamides. Sci. Rep. 13:7486. doi: 10.1038/s41598-023-34754-z
Bae, Y. C., Oh, J. M., Hwang, S. J., Shigenaga, Y., and Valtschanoff, J. G. (2004). Expression of vanilloid receptor TRPV1 in the rat trigeminal sensory nuclei. J. Comp. Neurol. 478, 62–71. doi: 10.1002/cne.20272
Bainbridge, M. N., Mazumder, A., Ogasawara, D., Abou Jamra, R., Bernard, G., Bertini, E., et al. (2022). Endocannabinoid dysfunction in neurological disease: neuro-ocular DAGLA-related syndrome. Brain 145, 3383–3390. doi: 10.1093/brain/awac223
Barchel, D., Stolar, O., De-Haan, T., Ziv-Baran, T., Saban, N., Fuchs, D. O., et al. (2018). Oral Cannabidiol use in children with autism Spectrum disorder to treat related symptoms and co-morbidities. Front. Pharmacol. 9:1521. doi: 10.3389/fphar.2018.01521
Barna, I., Zelena, D., Arszovszki, A. C., and Ledent, C. (2004). The role of endogenous cannabinoids in the hypothalamo-pituitary-adrenal axis regulation: in vivo and in vitro studies in CB1 receptor knockout mice. Life Sci. 75, 2959–2970. doi: 10.1016/j.lfs.2004.06.006
Basavarajappa, B. S. (2007). Critical enzymes involved in endocannabinoid metabolism. Protein Pept. Lett. 14, 237–246. doi: 10.2174/092986607780090829
Berna-Erro, A., Izquierdo-Serra, M., Sepúlveda, R. V., Rubio-Moscardo, F., Doñate-Macián, P., Serra, S. A., et al. (2017). Structural determinants of 5′,6′-epoxyeicosatrienoic acid binding to and activation of TRPV4 channel. Sci. Rep. 7:10522. doi: 10.1038/s41598-017-11274-1
Bilbao, A., and Spanagel, R. (2022). Medical cannabinoids: a pharmacology-based systematic review and meta-analysis for all relevant medical indications. BMC Med. 20:259. doi: 10.1186/s12916-022-02459-1
Bioque, M., García-Bueno, B., MacDowell, K. S., Meseguer, A., Saiz, P. A., Parellada, M., et al. (2013). Peripheral endocannabinoid system dysregulation in first-episode psychosis. Neuropsychopharmacology 38, 2568–2577. doi: 10.1038/npp.2013.165
Bloch Priel, S., Yitzhaky, A., Gurwitz, D., and Hertzberg, L. (2023). Cannabinoid receptor gene CNR1 is downregulated in subcortical brain samples and upregulated in blood samples of individuals with schizophrenia: a participant data systematic meta-analysis. Eur. J. Neurosci. 58, 3540–3554. doi: 10.1111/ejn.16122
Borgan, F., Beck, K., Butler, E., McCutcheon, R., Veronese, M., Vernon, A., et al. (2019a). The effects of cannabinoid 1 receptor compounds on memory: a meta-analysis and systematic review across species. Psychopharmacology 236, 3257–3270. doi: 10.1007/s00213-019-05283-3
Borgan, F., Laurikainen, H., Veronese, M., Marques, T. R., Haaparanta-Solin, M., Solin, O., et al. (2019b). In vivo availability of cannabinoid 1 receptor levels in patients with first-episode psychosis. JAMA Psychiatry 76, 1074–1084. doi: 10.1001/jamapsychiatry.2019.1427
Borgan, F., Veronese, M., Reis Marques, T., Lythgoe, D. J., and Howes, O. (2021). Association between cannabinoid 1 receptor availability and glutamate levels in healthy controls and drug-free patients with first episode psychosis: a multi-modal PET and 1H-MRS study. Eur. Arch. Psychiatry Clin. Neurosci. 271, 677–687. doi: 10.1007/s00406-020-01191-2
Braunscheidel, K. M., Okas, M. P., Floresco, S. B., and Woodward, J. J. (2022). Cannabinoid receptor type 1 antagonists alter aspects of risk/reward decision making independent of toluene-mediated effects. Psychopharmacology 239, 1337–1347. doi: 10.1007/s00213-021-05914-8
Budd, G. E. (2013). The Cambrian explosion: the reconstruction of animal biodiversity.— by Douglas H. Erwin and James W. Valentine. Syst. Biol. 62, 915–917. doi: 10.1093/sysbio/syt043
Burkhardt, P., Colgren, J., Medhus, A., Digel, L., Naumann, B., Soto-Angel, J. J., et al. (2023). Syncytial nerve net in a ctenophore adds insights on the evolution of nervous systems. Science 380, 293–297. doi: 10.1126/science.ade5645
Castillo, P. E., Younts, T. J., Chávez, A. E., and Hashimotodani, Y. (2012). Endocannabinoid signaling and synaptic function. Neuron 76, 70–81. doi: 10.1016/j.neuron.2012.09.020
Caterina, M. J., Schumacher, M. A., Tominaga, M., Rosen, T. A., Levine, J. D., and Julius, D. (1997). The capsaicin receptor: a heat-activated ion channel in the pain pathway. Nature 389, 816–824. doi: 10.1038/39807
Ceccarini, J., De Hert, M., Van Winkel, R., Peuskens, J., Bormans, G., Kranaster, L., et al. (2013). Increased ventral striatal CB1 receptor binding is related to negative symptoms in drug-free patients with schizophrenia. NeuroImage 79, 304–312. doi: 10.1016/j.neuroimage.2013.04.052
Ceccarini, J., Koole, M., and Van Laere, K. (2022). Cannabinoid receptor availability modulates the magnitude of dopamine release in vivo in the human reward system: a preliminary multitracer positron emission tomography study. Addict. Biol. 27:e13167. doi: 10.1111/adb.13167
Chen, Z.-P., Wang, S., Zhao, X., Fang, W., Wang, Z., Ye, H., et al. (2023). Lipid-accumulated reactive astrocytes promote disease progression in epilepsy. Nat. Neurosci. 26, 542–554. doi: 10.1038/s41593-023-01288-6
Chen, Z., Zhou, C., Yuan, X., and Xiao, S. (2019). Death march of a segmented and trilobate bilaterian elucidates early animal evolution. Nature 573, 412–415. doi: 10.1038/s41586-019-1522-7
Chung, H., Fierro, A., and Pessoa-Mahana, C. D. (2019). Cannabidiol binding and negative allosteric modulation at the cannabinoid type 1 receptor in the presence of delta-9-tetrahydrocannabinol: an in silico study. PLoS One 14:e0220025. doi: 10.1371/journal.pone.0220025
Clough, E., and Barrett, T. (2016). The gene expression omnibus database. Methods Mol. Biol. 1418, 93–110. doi: 10.1007/978-1-4939-3578-9_5
Costa, B., Giagnoni, G., Franke, C., Trovato, A. E., and Colleoni, M. (2004). Vanilloid TRPV1 receptor mediates the antihyperalgesic effect of the nonpsychoactive cannabinoid, cannabidiol, in a rat model of acute inflammation. Br. J. Pharmacol. 143, 247–250. doi: 10.1038/sj.bjp.0705920
Cristino, L., Bisogno, T., and Di Marzo, V. (2020). Cannabinoids and the expanded endocannabinoid system in neurological disorders. Nat. Rev. Neurol. 16, 9–29. doi: 10.1038/s41582-019-0284-z
D’Souza, D. C., Cortes-Briones, J., Creatura, G., Bluez, G., Thurnauer, H., Deaso, E., et al. (2019). Efficacy and safety of a fatty acid amide hydrolase inhibitor (PF-04457845) in the treatment of cannabis withdrawal and dependence in men: a double-blind, placebo-controlled, parallel group, phase 2a single-site randomised controlled trial. Lancet Psychiatry 6, 35–45. doi: 10.1016/S2215-0366(18)30427-9
Dalmaijer, E. S., Nord, C. L., and Astle, D. E. (2022). Statistical power for cluster analysis. BMC Bioinformatics 23:205. doi: 10.1186/s12859-022-04675-1
Dalton, V. S., Long, L. E., Weickert, C. S., and Zavitsanou, K. (2011). Paranoid schizophrenia is characterized by increased CB1 receptor binding in the dorsolateral prefrontal cortex. Neuropsychopharmacology 36, 1620–1630. doi: 10.1038/npp.2011.43
Damasceno, S., Gómez-Nieto, R., Garcia-Cairasco, N., Herrero-Turrión, M. J., Marín, F., and Lopéz, D. E. (2020). Top common differentially expressed genes in the epileptogenic nucleus of two strains of rodents susceptible to Audiogenic seizures: WAR and GASH/Sal. Front. Neurol. 11:33. doi: 10.3389/fneur.2020.00033
de la Harpe, A., Beukes, N., and Frost, C. L. (2022). CBD activation of TRPV1 induces oxidative signaling and subsequent ER stress in breast cancer cell lines. Biotechnol. Appl. Biochem. 69, 420–430. doi: 10.1002/bab.2119
De Marchi, N., De Petrocellis, L., Orlando, P., Daniele, F., Fezza, F., and Di Marzo, V. (2003). Endocannabinoid signalling in the blood of patients with schizophrenia. Lipids Health Dis. 2:5. doi: 10.1186/1476-511X-2-5
de Miguel-Yanes, J. M., Manning, A. K., Shrader, P., McAteer, J. B., Goel, A., Hamsten, A., et al. (2011). Variants at the endocannabinoid receptor CB1 gene (CNR1) and insulin sensitivity, type 2 diabetes, and coronary heart disease. Obesity (Silver Spring) 19, 2031–2037. doi: 10.1038/oby.2011.135
De Petrocellis, L., Orlando, P., Moriello, A. S., Aviello, G., Stott, C., Izzo, A. A., et al. (2012). Cannabinoid actions at TRPV channels: effects on TRPV3 and TRPV4 and their potential relevance to gastrointestinal inflammation. Acta Physiol (Oxf.) 204, 255–266. doi: 10.1111/j.1748-1716.2011.02338.x
Dean, B., Sundram, S., Bradbury, R., Scarr, E., and Copolov, D. (2001). Studies on [3H]CP-55940 binding in the human central nervous system: regional specific changes in density of cannabinoid-1 receptors associated with schizophrenia and cannabis use. Neuroscience 103, 9–15. doi: 10.1016/S0306-4522(00)00552-2
Dong, X., Duan, B., Liu, J., and Donoghue, P. C. J. (2022). Internal anatomy of a fossilized embryonic stage of the Cambrian-Ordovician scalidophoran Markuelia. R. Soc. Open Sci. 9:220115. doi: 10.1098/rsos.220115
Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30, 207–210. doi: 10.1093/nar/30.1.207
Egertová, M., and Elphick, M. R. (2007). Localization of CiCBR in the invertebrate chordate Ciona intestinalis: evidence of an ancient role for cannabinoid receptors as axonal regulators of neuronal signalling. J. Comp. Neurol. 502, 660–672. doi: 10.1002/cne.21331
Eggan, S. M., Hashimoto, T., and Lewis, D. A. (2008). Reduced cortical cannabinoid 1 receptor messenger RNA and protein expression in schizophrenia. Arch. Gen. Psychiatry 65, 772–784. doi: 10.1001/archpsyc.65.7.772
Elphick, M. R. (1998). An invertebrate G-protein coupled receptor is a chimeric cannabinoid/melanocortin receptor. Brain Res. 780, 168–171. doi: 10.1016/S0006-8993(97)01297-3
Elphick, M. R. (2012). The evolution and comparative neurobiology of endocannabinoid signalling. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 367, 3201–3215. doi: 10.1098/rstb.2011.0394
Elphick, M. R., Satou, Y., and Satoh, N. (2003). The invertebrate ancestry of endocannabinoid signalling: an orthologue of vertebrate cannabinoid receptors in the urochordate Ciona intestinalis. Gene 302, 95–101. doi: 10.1016/s0378-1119(02)01094-6
Etemad, L., Karimi, G., Alavi, M. S., and Roohbakhsh, A. (2022). Pharmacological effects of cannabidiol by transient receptor potential channels. Life Sci. 300:120582. doi: 10.1016/j.lfs.2022.120582
Evans, S. D., Hughes, I. V., Gehling, J. G., and Droser, M. L. (2020). Discovery of the oldest bilaterian from the Ediacaran of South Australia. Proc. Natl. Acad. Sci. 117, 7845–7850. doi: 10.1073/pnas.2001045117
French, E. D., Dillon, K., and Wu, X. (1997). Cannabinoids excite dopamine neurons in the ventral tegmentum and substantia nigra. Neuroreport 8, 649–652. doi: 10.1097/00001756-199702100-00014
Friend, L., Weed, J., Sandoval, P., Nufer, T., Ostlund, I., and Edwards, J. G. (2017). CB1-dependent Long-term depression in ventral tegmental area GABA neurons: a novel target for marijuana. J. Neurosci. 37, 10943–10954. doi: 10.1523/JNEUROSCI.0190-17.2017
Galaj, E., Bi, G.-H., Yang, H.-J., and Xi, Z.-X. (2020). Cannabidiol attenuates the rewarding effects of cocaine in rats by CB2, 5-HT1A and TRPV1 receptor mechanisms. Neuropharmacology 167:107740. doi: 10.1016/j.neuropharm.2019.107740
Gallego-Landin, I., García-Baos, A., Castro-Zavala, A., and Valverde, O. (2021). Reviewing the role of the endocannabinoid system in the pathophysiology of depression. Front. Pharmacol. 12:762738. doi: 10.3389/fphar.2021.762738
Garani, R., Watts, J. J., and Mizrahi, R. (2021). Endocannabinoid system in psychotic and mood disorders, a review of human studies. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 106:110096. doi: 10.1016/j.pnpbp.2020.110096
Gascuel, O. (1997). BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol. Biol. Evol. 14, 685–695. doi: 10.1093/oxfordjournals.molbev.a025808
Gau, P., Poon, J., Ufret-Vincenty, C., Snelson, C. D., Gordon, S. E., Raible, D. W., et al. (2013). The zebrafish Ortholog of TRPV1 is required for heat-induced locomotion. J. Neurosci. 33, 5249–5260. doi: 10.1523/JNEUROSCI.5403-12.2013
Giuffrida, A., Leweke, F. M., Gerth, C. W., Schreiber, D., Koethe, D., Faulhaber, J., et al. (2004). Cerebrospinal anandamide levels are elevated in acute schizophrenia and are inversely correlated with psychotic symptoms. Neuropsychopharmacology 29, 2108–2114. doi: 10.1038/sj.npp.1300558
Gochman, A., Tan, X.-F., Bae, C., Chen, H., Swartz, K. J., and Jara-Oseguera, A. (2023). Cannabidiol sensitizes TRPV2 channels to activation by 2-APB. elife 12:e86166. doi: 10.7554/eLife.86166
Gomis-González, M., Busquets-Garcia, A., Matute, C., Maldonado, R., Mato, S., and Ozaita, A. (2016). Possible therapeutic doses of cannabinoid type 1 receptor antagonist reverses key alterations in fragile X syndrome mouse model. Genes (Basel) 7:56. doi: 10.3390/genes7090056
Grueter, B. A., Brasnjo, G., and Malenka, R. C. (2010). Postsynaptic TRPV1 triggers cell type-specific long-term depression in the nucleus accumbens. Nat. Neurosci. 13, 1519–1525. doi: 10.1038/nn.2685
Habib, A. M., Okorokov, A. L., Hill, M. N., Bras, J. T., Lee, M.-C., Li, S., et al. (2019). Microdeletion in a FAAH pseudogene identified in a patient with high anandamide concentrations and pain insensitivity. Br. J. Anaesth. 123, e249–e253. doi: 10.1016/j.bja.2019.02.019
Han, X., Liang, Y., Hempel, B., Jordan, C. J., Shen, H., Bi, G.-H., et al. (2023). Cannabinoid CB1 receptors are expressed in a subset of dopamine neurons and underlie cannabinoid-induced aversion, Hypoactivity, and anxiolytic effects in mice. J. Neurosci. 43, 373–385. doi: 10.1523/JNEUROSCI.1493-22.2022
Haq, I., and Kilaru, A. (2020). An endocannabinoid catabolic enzyme FAAH and its paralogs in an early land plant reveal evolutionary and functional relationship with eukaryotic orthologs. Sci. Rep. 10:3115. doi: 10.1038/s41598-020-59948-7
Heger, P., Zheng, W., Rottmann, A., Panfilio, K. A., and Wiehe, T. (2020). The genetic factors of bilaterian evolution. elife 9:e45530. doi: 10.7554/eLife.45530
Herman, A. I., Kranzler, H. R., Cubells, J. F., Gelernter, J., and Covault, J. (2006). Association study of the CNR1 gene exon 3 alternative promoter region polymorphisms and substance dependence. Am. J. Med. Genet. B Neuropsychiatr. Genet. 141B, 499–503. doi: 10.1002/ajmg.b.30325
Hjorthøj, C., Posselt, C. M., and Nordentoft, M. (2021). Development over time of the population-attributable risk fraction for Cannabis use disorder in schizophrenia in Denmark. JAMA Psychiatry 78, 1013–1019. doi: 10.1001/jamapsychiatry.2021.1471
Holland, L. Z., Carvalho, J. E., Escriva, H., Laudet, V., Schubert, M., Shimeld, S. M., et al. (2013). Evolution of bilaterian central nervous systems: a single origin? EvoDevo 4:27. doi: 10.1186/2041-9139-4-27
Hua, T., Li, X., Wu, L., Iliopoulos-Tsoutsouvas, C., Wang, Y., Wu, M., et al. (2020). Activation and signaling mechanism revealed by cannabinoid receptor-Gi complex structures. Cell 180, 655–665.e18. doi: 10.1016/j.cell.2020.01.008
Hua, T., Vemuri, K., Nikas, S. P., Laprairie, R. B., Wu, Y., Qu, L., et al. (2017). Crystal structures of agonist-bound human cannabinoid receptor CB1. Nature 547, 468–471. doi: 10.1038/nature23272
Iannotti, F. A., Hill, C. L., Leo, A., Alhusaini, A., Soubrane, C., Mazzarella, E., et al. (2014). Nonpsychotropic plant cannabinoids, Cannabidivarin (CBDV) and Cannabidiol (CBD), activate and desensitize transient receptor potential Vanilloid 1 (TRPV1) channels in vitro: potential for the treatment of neuronal Hyperexcitability. ACS Chem. Neurosci. 5, 1131–1141. doi: 10.1021/cn5000524
Ibarra-Lecue, I., Unzueta-Larrinaga, P., Barrena-Barbadillo, R., Villate, A., Horrillo, I., Mendivil, B., et al. (2022). Cannabis use selectively modulates circulating biomarkers in the blood of schizophrenia patients. Addict. Biol. 27:e13233. doi: 10.1111/adb.13233
Jacobson, M. R., Watts, J. J., Da Silva, T., Tyndale, R. F., Rusjan, P. M., Houle, S., et al. (2021). Fatty acid amide hydrolase is lower in young cannabis users. Addict. Biol. 26:e12872. doi: 10.1111/adb.12872
Jakowiecki, J., and Filipek, S. (2016). Hydrophobic ligand entry and exit pathways of the CB1 cannabinoid receptor. J. Chem. Inf. Model 56, 2457–2466. doi: 10.1021/acs.jcim.6b00499
Jékely, G., Paps, J., and Nielsen, C. (2015). The phylogenetic position of ctenophores and the origin(s) of nervous systems. EvoDevo 6:1. doi: 10.1186/2041-9139-6-1
Jenko, K. J., Hirvonen, J., Henter, I. D., Anderson, K. B., Zoghbi, S. S., Hyde, T. M., et al. (2012). Binding of a tritiated inverse agonist to cannabinoid CB1 receptors is increased in patients with schizophrenia. Schizophr. Res. 141, 185–188. doi: 10.1016/j.schres.2012.07.021
Jia, Y.-F., Li, Y.-C., Tang, Y.-P., Cao, J., Wang, L.-P., Yang, Y.-X., et al. (2015). Interference of TRPV1 function altered the susceptibility of PTZ-induced seizures. Front. Cell Neurosci. 9:20. doi: 10.3389/fncel.2015.00020
Jones, D. T., Taylor, W. R., and Thornton, J. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275–282. doi: 10.1093/bioinformatics/8.3.275
Jouroukhin, Y., Zhu, X., Shevelkin, A. V., Hasegawa, Y., Abazyan, B., Saito, A., et al. (2019). Adolescent Δ9-tetrahydrocannabinol exposure and astrocyte-specific genetic vulnerability converge on nuclear factor-κB-Cyclooxygenase-2 signaling to impair memory in adulthood. Biol. Psychiatry 85, 891–903. doi: 10.1016/j.biopsych.2018.07.024
Juhasz, G., Lazary, J., Chase, D., Pegg, E., Downey, D., Toth, Z. G., et al. (2009). Variations in the cannabinoid receptor 1 gene predispose to migraine. Neurosci. Lett. 461, 116–120. doi: 10.1016/j.neulet.2009.06.021
Juknat, A., Pietr, M., Kozela, E., Rimmerman, N., Levy, R., Gao, F., et al. (2013). Microarray and pathway analysis reveal distinct mechanisms underlying cannabinoid-mediated modulation of LPS-induced activation of BV-2 microglial cells. PLoS One 8:e61462. doi: 10.1371/journal.pone.0061462
Jung, K.-M., Sepers, M., Henstridge, C. M., Lassalle, O., Neuhofer, D., Martin, H., et al. (2012). Uncoupling of the endocannabinoid signalling complex in a mouse model of fragile X syndrome. Nat. Commun. 3:1080. doi: 10.1038/ncomms2045
Kano, M., Ohno-Shosaku, T., Hashimotodani, Y., Uchigashima, M., and Watanabe, M. (2009). Endocannabinoid-mediated control of synaptic transmission. Physiol. Rev. 89, 309–380. doi: 10.1152/physrev.00019.2008
Karhson, D. S., Krasinska, K. M., Dallaire, J. A., Libove, R. A., Phillips, J. M., Chien, A. S., et al. (2018). Plasma anandamide concentrations are lower in children with autism spectrum disorder. Mol. Autism. 9:18. doi: 10.1186/s13229-018-0203-y
Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N., and Sternberg, M. J. E. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10, 845–858. doi: 10.1038/nprot.2015.053
Kelley, L. A., and Sternberg, M. J. E. (2009). Protein structure prediction on the web: a case study using the Phyre server. Nat. Protoc. 4, 363–371. doi: 10.1038/nprot.2009.2
Kibret, B. G., Ishiguro, H., Horiuchi, Y., and Onaivi, E. S. (2022). New insights and potential therapeutic targeting of CB2 cannabinoid receptors in CNS disorders. Int. J. Mol. Sci. 23:975. doi: 10.3390/ijms23020975
Kong, W., Wang, X., Yang, X., Huang, W., Han, S., Yin, J., et al. (2019). Activation of TRPV1 contributes to recurrent febrile seizures via inhibiting the microglial M2 phenotype in the immature brain. Front. Cell. Neurosci. 13:442. doi: 10.3389/fncel.2019.00442
Krishna Kumar, K., Robertson, M. J., Thadhani, E., Wang, H., Suomivuori, C.-M., Powers, A. S., et al. (2023). Structural basis for activation of CB1 by an endocannabinoid analog. Nat. Commun. 14:2672. doi: 10.1038/s41467-023-37864-4
Krogh, A., Larsson, B., von Heijne, G., and Sonnhammer, E. L. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J. Mol. Biol. 305, 567–580. doi: 10.1006/jmbi.2000.4315
Kucera, R., Bouskila, J., Elkrief, L., Fink-Jensen, A., Palmour, R., Bouchard, J.-F., et al. (2018). Expression and localization of CB1R, NAPE-PLD, and FAAH in the vervet monkey nucleus accumbens. Sci. Rep. 8:8689. doi: 10.1038/s41598-018-26826-2
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Kwan Cheung, K. A., Peiris, H., Wallace, G., Holland, O. J., and Mitchell, M. D. (2019). The interplay between the endocannabinoid system, epilepsy and cannabinoids. Int. J. Mol. Sci. 20:6079. doi: 10.3390/ijms20236079
Kwon, D. H., Zhang, F., Fedor, J. G., Suo, Y., and Lee, S.-Y. (2022). Vanilloid-dependent TRPV1 opening trajectory from cryoEM ensemble analysis. Nat. Commun. 13:2874. doi: 10.1038/s41467-022-30602-2
Lam, P. M. W., McDonald, J., and Lambert, D. G. (2005). Characterization and comparison of recombinant human and rat TRPV1 receptors: effects of exo- and endocannabinoids†. Br. J. Anaesth. 94, 649–656. doi: 10.1093/bja/aei098
Landucci, E., Pellegrini-Giampietro, D. E., Gianoncelli, A., and Ribaudo, G. (2022). Cannabidiol preferentially binds TRPV2: a novel mechanism of action. Neural Regen. Res. 17, 2693–2694. doi: 10.4103/1673-5374.335821
Laprairie, R. B., Bagher, A. M., Kelly, M. E. M., and Denovan-Wright, E. M. (2015). Cannabidiol is a negative allosteric modulator of the cannabinoid CB1 receptor. Br. J. Pharmacol. 172, 4790–4805. doi: 10.1111/bph.13250
Lauckner, J. E., Jensen, J. B., Chen, H.-Y., Lu, H.-C., Hille, B., and Mackie, K. (2008). GPR55 is a cannabinoid receptor that increases intracellular calcium and inhibits M current. Proc. Natl. Acad. Sci. USA 105, 2699–2704. doi: 10.1073/pnas.0711278105
Lazarini-Lopes, W., Silva-Cardoso, G. K., Leite-Panissi, C. R. A., and Garcia-Cairasco, N. (2022). Increased TRPV1 channels and FosB protein expression are associated with chronic epileptic seizures and Anxiogenic-like behaviors in a preclinical model of temporal lobe epilepsy. Biomedicines 10:416. doi: 10.3390/biomedicines10020416
Lee, T.-H., Lee, J.-G., Yon, J.-M., Oh, K.-W., Baek, I.-J., Nahm, S.-S., et al. (2011). Capsaicin prevents kainic acid-induced epileptogenesis in mice. Neurochem. Int. 58, 634–640. doi: 10.1016/j.neuint.2011.01.027
Leweke, F. M., Giuffrida, A., Koethe, D., Schreiber, D., Nolden, B. M., Kranaster, L., et al. (2007). Anandamide levels in cerebrospinal fluid of first-episode schizophrenic patients: impact of cannabis use. Schizophr. Res. 94, 29–36. doi: 10.1016/j.schres.2007.04.025
Leys, S. P. (2015). Elements of a “nervous system” in sponges. J. Exp. Biol. 218, 581–591. doi: 10.1242/jeb.110817
Liu, S., Paknejad, N., Zhu, L., Kihara, Y., Ray, M., Chun, J., et al. (2022). Differential activation mechanisms of lipid GPCRs by lysophosphatidic acid and sphingosine 1-phosphate. Nat. Commun. 13:731. doi: 10.1038/s41467-022-28417-2
Lu, H.-C., and Mackie, K. (2016). An introduction to the endogenous cannabinoid system. Biol. Psychiatry 79, 516–525. doi: 10.1016/j.biopsych.2015.07.028
Lu, A. T., Ogdie, M. N., Järvelin, M.-R., Moilanen, I. K., Loo, S. K., McCracken, J. T., et al. (2008). Association of the cannabinoid receptor gene (CNR1) with ADHD and post-traumatic stress disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 147B, 1488–1494. doi: 10.1002/ajmg.b.30693
Ludányi, A., Eross, L., Czirják, S., Vajda, J., Halász, P., Watanabe, M., et al. (2008). Downregulation of the CB1 cannabinoid receptor and related molecular elements of the endocannabinoid system in epileptic human hippocampus. J. Neurosci. 28, 2976–2990. doi: 10.1523/JNEUROSCI.4465-07.2008
Lupica, C. R., Riegel, A. C., and Hoffman, A. F. (2004). Marijuana and cannabinoid regulation of brain reward circuits. Br. J. Pharmacol. 143, 227–234. doi: 10.1038/sj.bjp.0705931
Martin, M. W., Grazhdankin, D. V., Bowring, S. A., Evans, D. A., Fedonkin, M. A., and Kirschvink, J. L. (2000). Age of Neoproterozoic bilatarian body and trace fossils, White Sea, Russia: implications for metazoan evolution. Science 288, 841–845. doi: 10.1126/science.288.5467.841
Martín-Durán, J. M., Pang, K., Børve, A., Lê, H. S., Furu, A., Cannon, J. T., et al. (2018). Convergent evolution of bilaterian nerve cords. Nature 553, 45–50. doi: 10.1038/nature25030
Martínez-Aguirre, C., Cinar, R., and Rocha, L. (2022). Targeting endocannabinoid system in epilepsy: for good or for bad. Neuroscience 482, 172–185. doi: 10.1016/j.neuroscience.2021.12.013
Matsuda, L. A., Lolait, S. J., Brownstein, M. J., Young, A. C., and Bonner, T. I. (1990). Structure of a cannabinoid receptor and functional expression of the cloned cDNA. Nature 346, 561–564. doi: 10.1038/346561a0
Matsuyama, M., Terada, Y., Yamazaki-Ito, T., and Ito, K. (2021). A luminescence-based human TRPV1 assay system for quantifying pungency in spicy foods. Food Secur. 10:151. doi: 10.3390/foods10010151
Mátyás, F., Urbán, G. M., Watanabe, M., Mackie, K., Zimmer, A., Freund, T. F., et al. (2008). Identification of the sites of 2-arachidonoylglycerol synthesis and action imply retrograde endocannabinoid signaling at both GABAergic and glutamatergic synapses in the ventral tegmental area. Neuropharmacology 54, 95–107. doi: 10.1016/j.neuropharm.2007.05.028
McHugh, D., Hu, S. S. J., Rimmerman, N., Juknat, A., Vogel, Z., Walker, J. M., et al. (2010). N-arachidonoyl glycine, an abundant endogenous lipid, potently drives directed cellular migration through GPR18, the putative abnormal cannabidiol receptor. BMC Neurosci. 11:44. doi: 10.1186/1471-2202-11-44
McPartland, J. M. (2004). Phylogenomic and chemotaxonomic analysis of the endocannabinoid system. Brain Res. Brain Res. Rev. 45, 18–29. doi: 10.1016/j.brainresrev.2003.11.005
McPartland, J. M., Glass, M., and Pertwee, R. G. (2007). Meta-analysis of cannabinoid ligand binding affinity and receptor distribution: interspecies differences. Br. J. Pharmacol. 152, 583–593. doi: 10.1038/sj.bjp.0707399
McPartland, J. M., and Guy, G. (2004). “The evolution of Cannabis and coevolution with the cannabinoid receptor - a hypothesis” in The medicinal use of Cannabis. eds. G. W. Guy, R. Robson, K. Strong, and B. Whittle (London: Royal Society of Pharmacists), 71–102.
McPartland, J. M., Matias, I., Di Marzo, V., and Glass, M. (2006). Evolutionary origins of the endocannabinoid system. Gene 370, 64–74. doi: 10.1016/j.gene.2005.11.004
Medvedeva, Y. V., Kim, M.-S., and Usachev, Y. M. (2008). Mechanisms of prolonged presynaptic Ca2+ signaling and glutamate release induced by TRPV1 activation in rat sensory neurons. J. Neurosci. 28, 5295–5311. doi: 10.1523/JNEUROSCI.4810-07.2008
Melancia, F., Schiavi, S., Servadio, M., Cartocci, V., Campolongo, P., Palmery, M., et al. (2018). Sex-specific autistic endophenotypes induced by prenatal exposure to valproic acid involve anandamide signalling. Br. J. Pharmacol. 175, 3699–3712. doi: 10.1111/bph.14435
Méndez-Díaz, M., Amancio-Belmont, O., Estrada-González, V., Ruiz-Contreras, A. E., and Prospéro-García, O. (2019). CB1R mediates oleamide’s reward while 5HT2cR mediates aversion in the nucleus accumbens shell of rats. Neurosci. Lett. 706, 189–193. doi: 10.1016/j.neulet.2019.05.033
Minichino, A., Senior, M., Brondino, N., Zhang, S. H., Godlewska, B. R., Burnet, P. W. J., et al. (2019). Measuring disturbance of the endocannabinoid system in psychosis: a systematic review and Meta-analysis. JAMA Psychiatry 76, 914–923. doi: 10.1001/jamapsychiatry.2019.0970
Mohapatra, D. P., and Nau, C. (2003). Desensitization of capsaicin-activated currents in the Vanilloid receptor TRPV1 is decreased by the cyclic AMP-dependent protein kinase pathway*. J. Biol. Chem. 278, 50080–50090. doi: 10.1074/jbc.M306619200
Morales, P., Muller, C., Jagerovic, N., and Reggio, P. H. (2022). Targeting CB2 and TRPV1: computational approaches for the identification of dual modulators. Front. Mol. Biosci. 9:841190. doi: 10.3389/fmolb.2022.841190
Moran, M. M., Xu, H., and Clapham, D. E. (2004). TRP ion channels in the nervous system. Curr. Opin. Neurobiol. 14, 362–369. doi: 10.1016/j.conb.2004.05.003
Moroz, L. L. (2015). Convergent evolution of neural systems in ctenophores. J. Exp. Biol. 218, 598–611. doi: 10.1242/jeb.110692
Muguruza, C., Morentin, B., Meana, J. J., Alexander, S. P., and Callado, L. F. (2019). Endocannabinoid system imbalance in the postmortem prefrontal cortex of subjects with schizophrenia. J. Psychopharmacol. 33, 1132–1140. doi: 10.1177/0269881119857205
Muller, C., Lynch, D. L., Hurst, D. P., and Reggio, P. H. (2020). A closer look at anandamide interaction with TRPV1. Front. Mol. Biosci. 7:144. doi: 10.3389/fmolb.2020.00144
Muller, C., Morales, P., and Reggio, P. H. (2019). Cannabinoid ligands targeting TRP channels. Front. Mol. Neurosci. 11:487. doi: 10.3389/fnmol.2018.00487
Munro, S., Thomas, K. L., and Abu-Shaar, M. (1993). Molecular characterization of a peripheral receptor for cannabinoids. Nature 365, 61–65. doi: 10.1038/365061a0
Musser, J. M., Schippers, K. J., Nickel, M., Mizzon, G., Kohn, A. B., Pape, C., et al. (2021). Profiling cellular diversity in sponges informs animal cell type and nervous system evolution. Science 374, 717–723. doi: 10.1126/science.abj2949
Nabissi, M., Morelli, M. B., Santoni, M., and Santoni, G. (2013). Triggering of the TRPV2 channel by cannabidiol sensitizes glioblastoma cells to cytotoxic chemotherapeutic agents. Carcinogenesis 34, 48–57. doi: 10.1093/carcin/bgs328
Nadezhdin, K. D., Neuberger, A., Nikolaev, Y. A., Murphy, L. A., Gracheva, E. O., Bagriantsev, S. N., et al. (2021). Extracellular cap domain is an essential component of the TRPV1 gating mechanism. Nat. Commun. 12:2154. doi: 10.1038/s41467-021-22507-3
Navarro, D., Gasparyan, A., Navarrete, F., Torregrosa, A. B., Rubio, G., Marín-Mayor, M., et al. (2022). Molecular alterations of the endocannabinoid system in psychiatric disorders. Int. J. Mol. Sci. 23:4764. doi: 10.3390/ijms23094764
Neelands, T. R., Jarvis, M. F., Han, P., Faltynek, C. R., and Surowy, C. S. (2005). Acidification of rat TRPV1 alters the kinetics of capsaicin responses. Mol. Pain 1:1744-8069-1-28. doi: 10.1186/1744-8069-1-28
Neuberger, A., Oda, M., Nikolaev, Y. A., Nadezhdin, K. D., Gracheva, E. O., Bagriantsev, S. N., et al. (2023). Human TRPV1 structure and inhibition by the analgesic SB-366791. Nat. Commun. 14:2451. doi: 10.1038/s41467-023-38162-9
Newell, K. A., Deng, C., and Huang, X.-F. (2006). Increased cannabinoid receptor density in the posterior cingulate cortex in schizophrenia. Exp. Brain Res. 172, 556–560. doi: 10.1007/s00221-006-0503-x
O’Sullivan, S. E. (2016). An update on PPAR activation by cannabinoids. Br. J. Pharmacol. 173, 1899–1910. doi: 10.1111/bph.13497
Ortega-Hernández, J., Lerosey-Aubril, R., Losso, S. R., and Weaver, J. C. (2022). Neuroanatomy in a middle Cambrian mollisoniid and the ancestral nervous system organization of chelicerates. Nat. Commun. 13:410. doi: 10.1038/s41467-022-28054-9
Ostlund, I., Von Gunten, M., Smith, C., and Edwards, J. G. (2022). Chronic Δ9-tetrahydrocannabinol impact on plasticity, and differential activation requirement for CB1-dependent long-term depression in ventral tegmental area GABA neurons in adult versus young mice. Front. Neurosci. 16:1067493. doi: 10.3389/fnins.2022.1067493
Overton, H. A., Babbs, A. J., Doel, S. M., Fyfe, M. C. T., Gardner, L. S., Griffin, G., et al. (2006). Deorphanization of a G protein-coupled receptor for oleoylethanolamide and its use in the discovery of small-molecule hypophagic agents. Cell Metab. 3, 167–175. doi: 10.1016/j.cmet.2006.02.004
Parry, L., and Caron, J.-B. (2019). Canadia spinosa and the early evolution of the annelid nervous system. Sci. Adv. 5:eaax5858. doi: 10.1126/sciadv.aax5858
Pasierski, M., and Szulczyk, B. (2020). Capsaicin inhibits sodium currents and epileptiform activity in prefrontal cortex pyramidal neurons. Neurochem. Int. 135:104709. doi: 10.1016/j.neuint.2020.104709
Pates, S., Botting, J. P., Muir, L. A., and Wolfe, J. M. (2022). Ordovician opabiniid-like animals and the role of the proboscis in euarthropod head evolution. Nat. Commun. 13:6969. doi: 10.1038/s41467-022-34204-w
Pecoits, E., Konhauser, K. O., Aubet, N. R., Heaman, L. M., Veroslavsky, G., Stern, R. A., et al. (2012). Bilaterian burrows and grazing behavior at >585 million years ago. Science 336, 1693–1696. doi: 10.1126/science.1216295
Pezzoli, M., Elhamdani, A., Camacho, S., Meystre, J., González, S. M., le Coutre, J., et al. (2014). Dampened neural activity and abolition of epileptic-like activity in cortical slices by active ingredients of spices. Sci. Rep. 4:6825. doi: 10.1038/srep06825
Proudnikov, D., Kroslak, T., Sipe, J. C., Randesi, M., Li, D., Hamon, S., et al. (2010). Association of polymorphisms of the cannabinoid receptor (CNR1) and fatty acid amide hydrolase (FAAH) genes with heroin addiction: impact of long repeats of CNR1. Pharmacogenomics J. 10, 232–242. doi: 10.1038/tpj.2009.59
Qin, N., Neeper, M. P., Liu, Y., Hutchinson, T. L., Lubin, M. L., and Flores, C. M. (2008). TRPV2 is activated by Cannabidiol and mediates CGRP release in cultured rat dorsal root ganglion neurons. J. Neurosci. 28, 6231–6238. doi: 10.1523/JNEUROSCI.0504-08.2008
R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. Available at: http://www.R-project.org/
Ramaker, R. C., Bowling, K. M., Lasseigne, B. N., Hagenauer, M. H., Hardigan, A. A., Davis, N. S., et al. (2017). Post-mortem molecular profiling of three psychiatric disorders. Genome Med. 9:72. doi: 10.1186/s13073-017-0458-5
Ranganathan, M., Cortes-Briones, J., Radhakrishnan, R., Thurnauer, H., Planeta, B., Skosnik, P., et al. (2016). Reduced brain cannabinoid receptor availability in schizophrenia. Biol. Psychiatry 79, 997–1005. doi: 10.1016/j.biopsych.2015.08.021
Rodríguez, G., Neugebauer, N. M., Yao, K. L., Meltzer, H. Y., Csernansky, J. G., and Dong, H. (2017). Δ9-tetrahydrocannabinol (Δ9-THC) administration after neonatal exposure to phencyclidine potentiates schizophrenia-related behavioral phenotypes in mice. Pharmacol. Biochem. Behav. 159, 6–11. doi: 10.1016/j.pbb.2017.06.010
RStudio Team (2020). RStudio: Integrated development for R. Available at: http://www.rstudio.com/
Russo, P., Strazzullo, P., Cappuccio, F. P., Tregouet, D. A., Lauria, F., Loguercio, M., et al. (2007). Genetic variations at the endocannabinoid type 1 receptor gene (CNR1) are associated with obesity phenotypes in men. J. Clin. Endocrinol. Metab. 92, 2382–2386. doi: 10.1210/jc.2006-2523
Ryberg, E., Larsson, N., Sjögren, S., Hjorth, S., Hermansson, N.-O., Leonova, J., et al. (2007). The orphan receptor GPR55 is a novel cannabinoid receptor. Br. J. Pharmacol. 152, 1092–1101. doi: 10.1038/sj.bjp.0707460
Saffarzadeh, F., Eslamizade, M. J., Ghadiri, T., Modarres Mousavi, S. M., Hadjighassem, M., and Gorji, A. (2015). Effects of TRPV1 on the hippocampal synaptic plasticity in the epileptic rat brain. Synapse 69, 375–383. doi: 10.1002/syn.21825
Saito, S., Fukuta, N., Shingai, R., and Tominaga, M. (2011). Evolution of vertebrate transient receptor potential vanilloid 3 channels: opposite temperature sensitivity between mammals and western clawed frogs. PLoS Genet. 7:e1002041. doi: 10.1371/journal.pgen.1002041
Saito, S., Nakatsuka, K., Takahashi, K., Fukuta, N., Imagawa, T., Ohta, T., et al. (2012). Analysis of transient receptor potential ankyrin 1 (TRPA1) in frogs and lizards illuminates both nociceptive heat and chemical sensitivities and coexpression with TRP vanilloid 1 (TRPV1) in ancestral vertebrates. J. Biol. Chem. 287, 30743–30754. doi: 10.1074/jbc.M112.362194
Saito, S., Saito, C. T., Igawa, T., Takeda, N., Komaki, S., Ohta, T., et al. (2022). Evolutionary tuning of transient receptor potential Ankyrin 1 underlies the variation in heat avoidance behaviors among frog species inhabiting diverse thermal niches. Mol. Biol. Evol. 39:msac180. doi: 10.1093/molbev/msac180
Sanchis, P., Lavignolle, R., Abbate, M., Lage-Vickers, S., Vazquez, E., Cotignola, J., et al. (2021). Analysis workflow of publicly available RNA-sequencing datasets. STAR Protoc. 2:100478. doi: 10.1016/j.xpro.2021.100478
Sbarski, B., and Akirav, I. (2018). Chronic exposure to cannabinoids before an emotional trauma may have negative effects on emotional function. Eur. Neuropsychopharmacol. 28, 955–969. doi: 10.1016/j.euroneuro.2018.05.008
Schrödinger, L., and DeLano, W. (2020). PyMol. Available at: http://www.pymol.org/pymol
Seillier, A., Martinez, A. A., and Giuffrida, A. (2020). Differential effects of Δ9-tetrahydrocannabinol dosing on correlates of schizophrenia in the sub-chronic PCP rat model. PLoS One 15:e0230238. doi: 10.1371/journal.pone.0230238
Servadio, M., Melancia, F., Manduca, A., di Masi, A., Schiavi, S., Cartocci, V., et al. (2016). Targeting anandamide metabolism rescues core and associated autistic-like symptoms in rats prenatally exposed to valproic acid. Transl. Psychiatry 6:e902. doi: 10.1038/tp.2016.182
Shao, Z., Yin, J., Chapman, K., Grzemska, M., Clark, L., Wang, J., et al. (2016). High-resolution crystal structure of the human CB1 cannabinoid receptor. Nature 540, 602–606. doi: 10.1038/nature20613
Siniscalco, D., Sapone, A., Giordano, C., Cirillo, A., de Magistris, L., Rossi, F., et al. (2013). Cannabinoid receptor type 2, but not type 1, is up-regulated in peripheral blood mononuclear cells of children affected by autistic disorders. J. Autism Dev. Disord. 43, 2686–2695. doi: 10.1007/s10803-013-1824-9
Srivastava, P. K., van Eyll, J., Godard, P., Mazzuferi, M., Delahaye-Duriez, A., Van Steenwinckel, J., et al. (2018). A systems-level framework for drug discovery identifies Csf1R as an anti-epileptic drug target. Nat. Commun. 9:3561. doi: 10.1038/s41467-018-06008-4
Stark, T., Iannotti, F. A., Di Martino, S., Di Bartolomeo, M., Ruda-Kucerova, J., Piscitelli, F., et al. (2022). Early blockade of CB1 receptors ameliorates schizophrenia-like alterations in the neurodevelopmental MAM model of schizophrenia. Biomol. Ther. 12:108. doi: 10.3390/biom12010108
Starkus, J., Jansen, C., Shimoda, L. M. N., Stokes, A. J., Small-Howard, A. L., and Turner, H. (2019). Diverse TRPV1 responses to cannabinoids. Channels (Austin) 13, 172–191. doi: 10.1080/19336950.2019.1619436
Stefano, G. B., Salzet, B., and Salzet, M. (1997). Identification and characterization of the leech CNS cannabinoid receptor: coupling to nitric oxide release. Brain Res. 753, 219–224. doi: 10.1016/s0006-8993(96)01484-9
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards suite: from gene data mining to disease genome sequence analyzes. Curr. Protoc. Bioinformatics 54, 1.30.1–1.30.33. doi: 10.1002/cpbi.5
Straiker, A., Dvorakova, M., Bosquez-Berger, T., Blahos, J., and Mackie, K. (2023). A collection of cannabinoid-related negative findings from autaptic hippocampal neurons. Sci. Rep. 13:9610. doi: 10.1038/s41598-023-36710-3
Su, T., Yan, Y., Li, Q., Ye, J., and Pei, L. (2021). Endocannabinoid system unlocks the puzzle of autism treatment via microglia. Front. Psych. 12:734837. doi: 10.3389/fpsyt.2021.734837
Sugaya, Y., and Kano, M. (2022). Endocannabinoid-mediated control of neural circuit excitability and epileptic seizures. Front. Neural Circuits 15:781113. doi: 10.3389/fncir.2021.781113
Sylantyev, S., Jensen, T. P., Ross, R. A., and Rusakov, D. A. (2013). Cannabinoid- and lysophosphatidylinositol-sensitive receptor GPR55 boosts neurotransmitter release at central synapses. Proc. Natl. Acad. Sci. USA 110, 5193–5198. doi: 10.1073/pnas.1211204110
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Szűcs, E., Dvorácskó, S., Tömböly, C., Büki, A., Kékesi, G., Horváth, G., et al. (2016). Decreased CB receptor binding and cannabinoid signaling in three brain regions of a rat model of schizophrenia. Neurosci. Lett. 633, 87–93. doi: 10.1016/j.neulet.2016.09.020
Tanaka, M., Sackett, S., and Zhang, Y. (2020). Endocannabinoid modulation of microglial phenotypes in neuropathology. Front. Neurol. 11:87. doi: 10.3389/fneur.2020.00087
Thomas, A., Baillie, G. L., Phillips, A. M., Razdan, R. K., Ross, R. A., and Pertwee, R. G. (2007). Cannabidiol displays unexpectedly high potency as an antagonist of CB1 and CB2 receptor agonists in vitro. Br. J. Pharmacol. 150, 613–623. doi: 10.1038/sj.bjp.0707133
Tominaga, M., Caterina, M. J., Malmberg, A. B., Rosen, T. A., Gilbert, H., Skinner, K., et al. (1998). The cloned capsaicin receptor integrates multiple pain-producing stimuli. Neuron 21, 531–543. doi: 10.1016/s0896-6273(00)80564-4
Trott, O., and Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J. Comput. Chem. 31, 455–461. doi: 10.1002/jcc.21334
Urigüen, L., García-Fuster, M. J., Callado, L. F., Morentin, B., La Harpe, R., Casadó, V., et al. (2009). Immunodensity and mRNA expression of A2A adenosine, D2 dopamine, and CB1 cannabinoid receptors in postmortem frontal cortex of subjects with schizophrenia: effect of antipsychotic treatment. Psychopharmacology 206, 313–324. doi: 10.1007/s00213-009-1608-2
Vasquez, V., Mugarra, R. C., Bell, B., Lee, J., and Cordero-Morales, J. F. (2020). Inhibition of inositol Monophosphatase enhances TRPV1 function in vivo. Biophys. J. 118:410a. doi: 10.1016/j.bpj.2019.11.2321
Vetter, I., Wyse, B. D., Monteith, G. R., Roberts-Thomson, S. J., and Cabot, P. J. (2006). The μ opioid agonist morphine modulates potentiation of capsaicin-evoked TRPV1 responses through a cyclic AMP-dependent protein kinase a pathway. Mol. Pain 2:22. doi: 10.1186/1744-8069-2-22
Voets, T., Droogmans, G., Wissenbach, U., Janssens, A., Flockerzi, V., and Nilius, B. (2004). The principle of temperature-dependent gating in cold- and heat-sensitive TRP channels. Nature 430, 748–754. doi: 10.1038/nature02732
Wang, X., Yang, X.-L., Kong, W.-L., Zeng, M.-L., Shao, L., Jiang, G.-T., et al. (2019). TRPV1 translocated to astrocytic membrane to promote migration and inflammatory infiltration thus promotes epilepsy after hypoxic ischemia in immature brain. J. Neuroinflammation 16:214. doi: 10.1186/s12974-019-1618-x
Watanabe, H., Vriens, J., Prenen, J., Droogmans, G., Voets, T., and Nilius, B. (2003). Anandamide and arachidonic acid use epoxyeicosatrienoic acids to activate TRPV4 channels. Nature 424, 434–438. doi: 10.1038/nature01807
Wei, D., Dinh, D., Lee, D., Li, D., Anguren, A., Moreno-Sanz, G., et al. (2016). Enhancement of anandamide-mediated endocannabinoid signaling corrects autism-related social impairment. Cannabis Cannabinoid Res. 1, 81–89. doi: 10.1089/can.2015.0008
Wilson, R. I., and Nicoll, R. A. (2001). Endogenous cannabinoids mediate retrograde signalling at hippocampal synapses. Nature 410, 588–592. doi: 10.1038/35069076
Wong, D. F., Kuwabara, H., Horti, A. G., Raymont, V., Brasic, J., Guevara, M., et al. (2010). Quantification of cerebral cannabinoid receptors subtype 1 (CB1) in healthy subjects and schizophrenia by the novel PET radioligand [11C]OMAR. NeuroImage 52, 1505–1513. doi: 10.1016/j.neuroimage.2010.04.034
Wu, Y.-W., Kim, J.-I., Tawfik, V. L., Lalchandani, R. R., Scherrer, G., and Ding, J. B. (2015). Input- and cell type-specific endocannabinoid-dependent LTD in the striatum. Cell Rep. 10, 75–87. doi: 10.1016/j.celrep.2014.12.005
Wu, H.-F., Lu, T.-Y., Chu, M.-C., Chen, P. S., Lee, C.-W., and Lin, H.-C. (2020). Targeting the inhibition of fatty acid amide hydrolase ameliorate the endocannabinoid-mediated synaptic dysfunction in a valproic acid-induced rat model of autism. Neuropharmacology 162:107736. doi: 10.1016/j.neuropharm.2019.107736
Yañez-Guerra, L. A., Thiel, D., and Jékely, G. (2022). Premetazoan origin of neuropeptide signaling. Mol. Biol. Evol. 39:msac051. doi: 10.1093/molbev/msac051
Yang, J., Ortega-Hernández, J., Butterfield, N. J., Liu, Y., Boyan, G. S., Hou, J., et al. (2016). Fuxianhuiid ventral nerve cord and early nervous system evolution in Panarthropoda. Proc. Natl. Acad. Sci. 113, 2988–2993. doi: 10.1073/pnas.1522434113
Yang, J., Ortega-Hernández, J., Butterfield, N. J., and Zhang, X. (2013). Specialized appendages in fuxianhuiids and the head organization of early euarthropods. Nature 494, 468–471. doi: 10.1038/nature11874
Zamberletti, E., Gabaglio, M., Woolley-Roberts, M., Bingham, S., Rubino, T., and Parolaro, D. (2019). Cannabidivarin treatment ameliorates autism-like behaviors and restores hippocampal endocannabinoid system and glia alterations induced by prenatal Valproic acid exposure in rats. Front. Cell. Neurosci. 13:367. doi: 10.3389/fncel.2019.00367
Zavitsanou, K., Garrick, T., and Huang, X. F. (2004). Selective antagonist [3H]SR141716A binding to cannabinoid CB1 receptors is increased in the anterior cingulate cortex in schizophrenia. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 28, 355–360. doi: 10.1016/j.pnpbp.2003.11.005
Zheng, W., and Wen, H. (2019). Heat activation mechanism of TRPV1: new insights from molecular dynamics simulation. Temperature (Austin) 6, 120–131. doi: 10.1080/23328940.2019.1578634
Keywords: CB1R, TRPV1, endocannabinoid system, evolution, protein structure, phylogenetic tree, transcriptomics, bioinformatics
Citation: Cherry AL, Wheeler MJ, Mathisova K and Di Miceli M (2024) In silico analyses of the involvement of GPR55, CB1R and TRPV1: response to THC, contribution to temporal lobe epilepsy, structural modeling and updated evolution. Front. Neuroinform. 18:1294939. doi: 10.3389/fninf.2024.1294939
Edited by:
Gabor Juhasz, Eötvös Loránd University, HungaryReviewed by:
José Luis Castañeda-Cabral, University of Guadalajara, MexicoBalapal Basavarajappa, New York University, United States
Rui Milton Patricio Silva Junior, University of São Paulo, Brazil
Copyright © 2024 Cherry, Wheeler, Mathisova and Di Miceli. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mathieu Di Miceli, bS5kaW1pY2VsaUB3b3JjLmFjLnVr
†These authors have contributed equally to this work and share first authorship