Impact of cannabis use on brain metabolism using 31P and 1H magnetic resonance spectroscopy
Abstract
Purpose:
This prospective cross-sectional study investigated the influence of regular cannabis use on brain metabolism in young cannabis users by using combined proton and phosphorus magnetic resonance spectroscopy.
Methods:
The study was performed in 45 young cannabis users aged 18–30, who had been using cannabis on a regular basis over a period of at least 2 years and in 47 age-matched controls. We acquired 31P MRS data in different brain regions at 3T with a double-resonant 1H/31P head coil, anatomic images, and 1H MRS data with a standard 20-channel 1H head coil. Absolute concentration values of proton metabolites were obtained via calibration from tissue water as an internal reference, whereas a standard solution of 75 mmol/l KH2PO4 was used as an external reference for the calibration of phosphorus signals.
Results:
We found an overall but not statistically significant lower concentration level of several proton and phosphorus metabolites in cannabis users compared to non-users. In particular, energy-related phosphates such as adenosine triphosphate (ATP) and inorganic phosphate (Pi) were reduced in all regions under investigation. Phosphocreatine (PCr) showed lowered values mainly in the left basal ganglia and the left frontal white matter.
Conclusion:
The results suggest that the increased risk of functional brain disorders observed in long-term cannabis users could be caused by an impairment of the energy metabolism of the brain, but this needs to be verified in future studies.
Supplementary Information:
The online version contains supplementary material available at 10.1007/s00234-023-03220-y.
Article type: Research Article
Keywords: Cannabis, Marijuana, 1H MRS, 31P MRS, Brain metabolites
Affiliations: https://ror.org/01jdpyv68grid.11749.3a0000 0001 2167 7588Institute of Neuroradiology, Saarland University, 66421 Homburg, Germany; grid.4567.00000 0004 0483 2525Helmholtz Zentrum Munich, German Research Center for Environmental Health Institute of Biological and Medical Imaging, 85748 Munich, Germany; https://ror.org/01jdpyv68grid.11749.3a0000 0001 2167 7588Department of Psychiatry and Psychotherapy, Saarland University, 66421 Homburg, Germany; https://ror.org/038t36y30grid.7700.00000 0001 2190 4373Department of General Psychiatry at the Center for Psychosocial Medicine, Heidelberg University, 69115 Heidelberg, Germany; Department of Obstetrics and Gynecology, RKH Clinic Ludwigsburg, 71640 Ludwigsburg, Germany
License: © The Author(s) 2023 CC BY 4.0 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Article links: DOI: 10.1007/s00234-023-03220-y | PubMed: 37735222 | PMC: PMC10567915
Relevance: Moderate: mentioned 3+ times in text
Full text: PDF (1.7 MB)
Introduction
Cannabis is one of the most widely used recreational drugs in the world [ref. 1]. Even though there has been a concern over decades about the use of cannabis as a cause of psychiatric illness, cannabis-related disorders have been rising among the past years [ref. 2]. Partial legalization can be associated with the increasing usage and the reduction in the perception of harm [ref. 3]. Due to this development, more scientific evidence is needed to determine the degree of harmfulness to the human body, especially with respect to brain metabolism and the whole nervous system.
Delta-9-tetrahydrocannabinol (∆9-THC) is the main psychoactive component of cannabis, acting on cannabinoid (CB1) receptors which can be densely found within brain networks critical for learning, attention, memory, cognitive processing, and motor control [ref. 4]. Moderate to high concentrations of CB1-binding sites have been detected in the thalamus, cerebellum, amygdala, basal ganglia, occipito-temporal gyrus, inferior temporal gyrus, frontal cortex, and hippocampus [ref. 4–ref. 6].
Several studies have shown that long-term cannabis use negatively affects memory, motor skills, executive function, emotional processing, and attention in adolescents [ref. 7–ref. 9] and adults [ref. 10–ref. 14]. In neuroimaging studies, long-term cannabis users exhibited abnormal brain activation during performance of functional tasks, including decision-making, verbal list learning, visual attention, and response inhibition [ref. 15–ref. 18].
Proton MRS is a non-invasive technique that has been widely applied to detect and quantify important neurometabolites [ref. 19]. Using single-voxel or multi-voxel acquisition schemes, cerebral metabolites including NAA (N-acetyl-aspartate), Cr (creatine), and cytosolic choline (Cho) can be assessed. NAA plays a role as a biomarker indicating neuronal viability [ref. 20]. Total Cr (tCr, creatine plus phosphocreatine) is involved in the energy metabolism, acting as an energy buffer by distributing energy within the brain and by maintaining constant brain adenosine triphosphate (ATP) levels through the creatine kinase reaction [ref. 21, ref. 22]. The Cho signal is associated with cellular membrane synthesis and degradation.
Phosphorus MRS in addition allows in vivo evaluation of compounds directly related to the energy metabolism and the composition of cell membranes. Adenosine triphosphate (ATP), phosphocreatine (PCr), and inorganic phosphate (Pi) are linked to brain bioenergetics through biochemical energy production (i.e., ATP synthesis) and energy use (i.e., ATP utilization). The phosphomonoesters (PME) play an important role in the synthesis of membrane lipids such as phosphatidylcholine and phosphatidylethanolamine. The main PME constituents, phosphoethanolamine (PE) and phosphocholine (PC), are precursors of the corresponding phospholipids. Membrane breakdown, in turn, is indicated by the phosphodiester (PDE) and catabolic products of phospholipid metabolism, glycerol-phosphoethanolamine (GPE), and glycerol-phosphocholine (GPC). Decreased membrane turnover has been associated with elevated PDE levels [ref. 23]. PME reduction refers to altered membrane turnover rates. In bipolar depression, e.g., studies have shown significantly altered frontal lobe PME [ref. 24]. Furthermore, 31P MRS can detect nicotinamide adenine dinucleotide phosphate (NADP), which is involved in oxidative chains and in membrane phospholipid metabolism [ref. 25]. Finally, it is possible to obtain the value of intracellular pH as well as the concentration of magnesium (Mg2+) from the spectrum [ref. 26, ref. 27].
To the best of our knowledge, no literature is available on brain metabolic changes using 31P-MRS related to cannabis use. In this study, we used both single-voxel-1H MRS and multi-voxel-31P MRS at a 3 T scanner to determine absolute metabolite concentration values from five brain areas that are suspected to be affected by cannabis [ref. 28] including frontal gray (FGM) and frontal white matter (FWM), thalamus (TH), basal ganglia (BG), and temporal lobe (TL). With 31P MRS, all regions except FGM were evaluated in both hemispheres separately. 1H MRS voxels other than FGM were restricted to the right hemisphere. Comparing concentration data between long-term cannabis users and non-users, we detected considerable though not statistically significant differences which might help to better understand the impact of cannabis use on brain metabolism. In addition, sex-related differences in non-users were found.
Methods
Study subjects
The subjects were recruited through local drug counseling centers. None of the participants received treatment for substance-use disorder. Control subjects were recruited through advertisement (poster, flyer) at the hospital. Recruitment of both subject groups took place simultaneously. All subjects were interviewed by an experienced psychologists or psychiatrist to assess extent and history of their cannabis use and underwent a complex psychometric assessment (supplementary_Demographics: suppl_table 1) to ensure inclusion criteria as seen below [ref. 29].
Only right-handed study subjects and controls between 18 and 30 years without neurological, psychiatric, and systemic diseases and without further drug addictions were included. This restriction was meant to exclude the effects of handedness and medical conditions on brain metabolism.
We investigated 21 female non-users (fN) (age 23 ± 2) and 26 male non-users (mN) (age 25 ± 4), who had never been using cannabis before or less than 10 times in total (= lifetime consume).
In the consumer group, 5 female cannabis-users (fC) (age 24 ± 4) and 40 male cannabis-users (mC) (age 24 ± 3) were examined. All users had been using cannabis on a regular basis at least 1 day per month in the last 24 months.
The fC group was excluded from further evaluation because not sufficient subjects could be found during the study.
Before MRI scan, study participants had to remove all metal objects. The subjects were instructed to move as little as possible during the MR examination which lasted about 1.5 h. Smoking was prohibited on examination day.
In this study, we used the STROBE cross-sectional reporting guidelines [ref. 30].
Data acquisition
Data acquisition of the brain was performed on a 3T whole body system (Magnetom Skyra, Siemens Healthcare, Erlangen, Germany). For anatomic images and 1H MRS, the standard 20-channel 1H (receive-only) head coil was used due to quality reasons. 31P spectra were acquired using a double-resonant 1H/31P (transmit/receive) head coil (RAPID Biomedical GmbH, Rimpar, Germany).
Anatomical data included three orthogonal T2-weighted localizers and a sagittal 3D T1-weighted data set (resolution 0.9 mm × 0.9 mm × 0.9 mm) of the whole brain (MPRAGE) which allowed segmentation of the brain tissue to obtain compartment maps of gray matter, white matter, and CSF. Segmentation was obtained using the SPM software (SPM 8, statistical parametric mapping, The Wellcome Trust Centre for Neuroimaging, University College London). For all spectroscopic volumes of interest, volume fractions of the three compartments were calculated from the maps.
Single-voxel proton spectra were obtained from 4 different brain regions: frontal gray matter (FGM), right frontal white matter (r_FWM), right thalamus (r_TH), and right temporal region (r_TL). Because of time restrictions, no spectra were acquired from the left hemisphere. Mean voxel size was 15 ml, 12 ml, 10 ml, and 8 ml respectively. Depending on brain size, voxel size was individually slightly adjusted to ensure accurate coverage of the anatomical target region. We used a PRESS sequence with TR = 1500 ms, TE = 135 ms, 80 acquisitions, bandwidth = 1200 Hz, and vector size = 1024. Shim adjustment was corrected manually to achieve minimal line width. As tissue water was used as an internal reference for absolute quantification of metabolites, additional spectra without water suppression were acquired from each voxel.
After coil change and repositioning of the patient phosphorus spectra were recorded using a 3D-chemical-shift-imaging (CSI) free-induction-decay (FID) sequence (TR = 1200 ms, TE = 2.3 ms, 15 acquisitions, bandwidth = 2000 Hz, vector size = 1024). Elliptical phase encoding with a weighted acquisition scheme was employed. Matrix size was 8 × 8 × 8, FOV = 200 × 200 mm2 resulting in 25 × 25 × 25 mm3 voxels. Optimized signal intensity was achieved by applying proton decoupling using the WALTZ-4 scheme and by a reduced flip angle of 60°. Careful manual shimming of the 3D volume was applied yielding line widths lower than 30 Hz. Acquisition time was 8:42 min. For absolute quantification of phosphorus metabolites, a phosphorus phantom with 75 mmol/l KH2PO4 was used as an external reference. The phantom was placed in the headcoil close to the left fronto-parietal part of the head.
Further information concerning data quality of 31P MRS and 1H MRS spectra can be found in supplementary_P_DataQuality and supplementary_H_DataQuality, respectively.
Data processing
Evaluation of proton spectra was done using the commercial software tool LCModel [ref. 31] (http://s-provencher.com/lcmodel.shtml). The signal-to-noise ratio (SNR) and the value of Cramer-Rao lower bound (%SD) were used to discard low quality data. Only spectra with SNR higher than 3 and %SD lower than 20% both for Cr and Cho were included for further analysis. To obtain absolute metabolite concentration values, the LCModel output data were corrected for longitudinal and transversal relaxation of both metabolites and brain tissue water taking into account the fractions of GM, WM, and CSF determined separately for each voxel from the segmentation maps. Relative tissue water content of 78%, 65%, and 97% was assumed for GM, WM, and CSF, respectively [ref. 32]. Relaxation factors (RH) were calculated according to the following equation for double-spin-echo sequences [ref. 33]:
{S}_H={S}_{H0}\cdot {\mathrm{R}}_H
\]
with
{R}_H=\exp \left(-\frac{TE}{T2}\right)\left\{1-\exp \left(-\frac{TR}{T1}\right)+2\exp \left(\left[\left(\frac{TE}{2}\right)- TR\right]/T1\right)-2\exp \left(\left[\left(\frac{3\ TE}{2}\right)- TR\right]/T1\right)\right\}
\]
SH and SH0 represent the measured proton peak intensity and the peak intensity corrected for T1 and T2 relaxation, respectively. T1 and T2 values were chosen as mean values (supplementary_H_Results: suppl_table 6) from the literature [ref. 34–ref. 47]. Finally, concentrations in units of milli-mole (mMol) per kg of brain tissue were calculated by correcting all metabolite values for the CSF fraction of each spectroscopic voxel determined from the compartment maps.
Phosphorus CSI data were transferred to a Leonardo workstation (Siemens Healthcare GmbH, Erlangen, Germany) and interpolated to a 32 × 32 × 8 grid resulting in a stack of 8 axial slices with 25 mm thickness (voxel size 6.3 × 6.3 × 25 mm3 ≈ 1 ml) which were superimposed on axial T2-weighted slices (Fig. 1). Anatomical volumes of interest (VOI) for spectral evaluation were identified by manually (M.F., M.B., and S.B.) selecting appropriate voxels in the grid. Grid shift in-plane as well as in head-feet direction was applied to optimally enclose the respective anatomical region of interest. Nine different VOIs were delineated for each subject in FGM, l_FWM, r_FWM, l_TH, r_TH, l_BG, r_BG, l_TL, and r_TL. VOI size ranged from 7 to 18 ml.

Quantitative analysis of the 31P spectra was performed with the jMRUI software tool (version 5.1) employing the AMARES algorithm [ref. 48]. The model function was composed of 14 resonances including PE, PC, Pi, GPE, GPC, PCr, ATP, and one macromolecular component to account for the broad signal baseline (Fig. 1). ATP was represented by a total of 7 peaks: a doublet γ-ATP, a doublet α-ATP, and a triplet β-ATP. Constraints for frequency, damping, coupling constants, and amplitude ratios (prior knowledge) were defined for the compounds to be estimated by the algorithm. The resulting amplitude values are proportional to the corresponding metabolite concentration. The concentration of ATP was calculated from the γ-ATP resonance. Only spectra with SNR higher than 3 were included for further analysis.
The AMARES algorithm provides Cramer-Rao lower bound (sd.amp.) values as an error estimate for all peaks in each spectrum. Whereas PCr signals always had relative error values lower than 20%, weak signals, e.g., NAD and PC, suffer from low intensities and high errors. Peaks with relative error values > 1 were excluded from further analysis.
Several postprocessing steps are required to obtain absolute quantification of metabolites: first, the signal amplitudes were corrected for the reduced flip angle and for T1 relaxation. Correction factors (RP) were calculated using the following equation:
{S}_P={S}_{P0}\cdot {R}_P
\]
with
{R}_P=\frac{\sin (x)\left\{\ 1-\exp \left(-\frac{TR}{T1}\right)\ \right\}}{\left\{\ 1-\cos (x)\exp \left(-\frac{TR}{T1}\right)\right\}},\kern0.5em x=60{}^{\circ},\mathrm{TR}=1200\ \mathrm{ms}.
\]
SP and SP0 represent the measured phosphorus peak intensity and the peak intensity corrected for T1 relaxation and flip angle, respectively (supplementary_P_Results, suppl_table 4).
Varying coil loading due to different head sizes of subjects was taken into account based on the radiofrequency transmitter amplitude required for a 90° pulse. Calibration of signal intensities was done with the phantom replacement method [ref. 49]. Finally, the calculated metabolic concentrations were corrected for partial CSF volume of each VOI to obtain concentration values in units of mMol per kg of brain tissue.
Calculated parameters
Intracellular pH was calculated from the chemical shift difference δ between the peak of inorganic phosphate (Pi) and the PCr peak [ref. 50–ref. 52] according to the equation:
\mathrm{pH}=6.75+{\log}_{10}\left[\frac{3.27-\updelta}{\updelta -5.63}\right]
\]
Free cytosolic Mg2+ was estimated from the chemical shift difference δβ between the peak of β-ATP and the PCr peak according to the formula:
\mathrm{pMg}=4.24-{\log}_{10}\left[\frac{{\left({\delta}_{\upbeta}+18.58\right)}^{0.42}}{{\left(-15.74-{\delta}_{\upbeta}\right)}^{0.84}}\right]
\]
The relation between the concentration of Mg2+ (in mol/l) and the value of pMg is given by: [Mg2+] = − log10(pMg).
Concentration ratios of PCr and Cr were estimated for those ROIs, where both phosphorus and proton spectra were acquired: FWM, r_TH, r_TL, and r_FWM. Cr values were calculated as [Cr] = [tCr] − [PCr].
Statistical methods
All statistical evaluations were performed using IBM SPSS Statistics (Version 27). Mean values for each metabolite concentration as well as for pH and Mg were calculated for every VOI separately. Concentration differences between groups were determined as relative values in percent according to:
\Delta \mathrm{mf}=\frac{\left(\mathrm{mN}-\mathrm{fN}\right)}{\mathrm{fN}}, \Delta \mathrm{CN}=\frac{\left(\mathrm{mC}-\mathrm{mN}\right)}{\mathrm{mN}}
\]
The statistical analysis was based on the General Linear Model using multivariate analysis of variance (MANOVA).
For P MRS, the metabolite values of PME, Pi, PDE, PCr, ATP, pH, and Mg were set as dependent variables, while membership to one of the three groups (fN, mN, and mC) was set as a fixed factor. NAD was excluded from the analysis because of too many low-quality data. To investigate the overall effect of the groups on all seven metabolite values, a multivariate Wilks-Lambda test was used. Paired comparisons were performed by post hoc Scheffé test. The level of significance was corrected for multiple tests using the Bonferroni approach. We analyzed nine regions simultaneously, so a p < 0.0056 was chosen as the criterion for significance.
Side related differences in metabolite values were calculated as relative values in percent according to:
\Delta \mathrm{rl}=\frac{\left(\mathrm{right}-\mathrm{left}\right)}{\mathrm{left}}.
\]
For statistical evaluation, multivariate analysis of variance with repeated measurements (RM MANOVA) was used in four regions: TH, BG, TL, and FWM. To investigate overall hemispheric effects, the metabolite values from the left and right hemisphere were set as within-subject factors, membership to the groups was set as a between-subjects factor. As we made four bilateral comparisons, the level of significance was chosen as p < 0.0125 according to the Bonferroni approach. In order to compare side related effects between the groups, additional RM ANOVAs were performed for each group separately. Hemispheric differences for individual metabolite values were evaluated with paired t test.
For H MRS, the metabolic values of tNAA, tCr, and tCho were set as dependent variables in the MANOVA. To investigate the overall effect of the groups on all three metabolite values, a multivariate Wilks-Lambda test was used. Paired comparisons were performed by post hoc Scheffé test. The level of significance was corrected for multiple tests using the Bonferroni approach. We measured four regions, so a p < 0.0125 was chosen as the criterion for significance.
Results
31P MRS results are shown in Tables 1 and 2 and Figs. 2, 3, and 4; more detailed data can be found in supplementary_P_results (suppl_table 5 and suppl_figs. 9a–d). 1H MRS results are shown in Table 3 and Fig. 5; more detailed data are given in supplementary_H_results (suppl_table 7, suppl_table 8 and suppl_figs.10a–b).
Table 1: Results of 31P MRS, comparison between groups
| FGM | MANOVA: Wilks-Lambda: p = 0.175 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.70 | 2.61 | 2.62 | 0% | 0.997 |
| Pi | 0.68 | 0.59 | 0.54 | − 9% | 0.599 |
| PDE | 2.89 | 3.04 | 2.83 | − 7% | 0.259 |
| PCr | 3.65 | 3.58 | 3.60 | 1% | 0.984 |
| ATP | 2.83 | 2.70 | 2.51 | − 7% | 0.400 |
| pH | 6.98 | 6.99 | 6.98 | − 0.1% | 0.343 |
| Mg | 0.11 | 0.10 | 0.10 | 1% | 0.923 |
| r_TH | MANOVA: Wilks-Lambda: p = 0.014 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.30 | 1.96 | 2.05 | 4% | 0.680 |
| Pi | 0.79 | 0.75 | 0.73 | − 3% | 0.892 |
| PDE | 2.61 | 2.56 | 2.63 | 3% | 0.790 |
| PCr | 3.70 | 3.23 | 3.38 | 5% | 0.558 |
| ATP | 2.43 | 2.16 | 2.03 | − 6% | 0.586 |
| pH | 6.99 | 6.99 | 6.99 | 0.0% | 0.827 |
| Mg | 0.10 | 0.11 | 0.10 | − 5% | 0.469 |
| r_BG | MANOVA: Wilks-Lambda: p = 0.358 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.17 | 2.16 | 2.11 | − 2% | 0.874 |
| Pi | 0.61 | 0.63 | 0.53 | − 15% | 0.117 |
| PDE | 2.66 | 2.70 | 2.63 | − 2% | 0.846 |
| PCr | 3.48 | 3.22 | 3.29 | 2% | 0.865 |
| ATP | 2.49 | 2.36 | 2.15 | − 9% | 0.233 |
| pH | 6.99 | 6.99 | 6.99 | 0.0% | 0.998 |
| Mg | 0.11 | 0.11 | 0.11 | 0% | 0.991 |
| r_TL | MANOVA: Wilks-Lambda: p = 0.034 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.03 | 2.18 | 2.25 | 4% | 0.764 |
| Pi | 0.59 | 0.58 | 0.53 | − 9% | 0.453 |
| PDE | 2.03 | 2.36 | 2.20 | − 7% | 0.487 |
| PCr | 3.50 | 3.78 | 3.79 | 0% | 1.000 |
| ATP | 2.30 | 2.23 | 2.07 | − 7% | 0.367 |
| pH | 7.00 | 7.00 | 6.99 | − 0.1% | 0420 |
| Mg | 0.12 | 0.11 | 0.11 | − 2% | 0.854 |
| r_FWM | MANOVA: Wilks-Lambda: p = 0.216 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.27 | 2.35 | 2.30 | − 2% | 0.894 |
| Pi | 0.61 | 0.63 | 0.54 | − 15% | 0.145 |
| PDE | 2.56 | 2.69 | 2.44 | − 9% | 0.118 |
| PCr | 3.40 | 3.46 | 3.42 | − 1% | 0.941 |
| ATP | 2.60 | 2.47 | 2.27 | − 8% | 0.199 |
| pH | 6.98 | 6.99 | 6.98 | − 0.1% | 0.738 |
| Mg | 0.11 | 0.11 | 0.11 | 1% | 0.970 |
| l_TH | MANOVA: Wilks-Lambda: p = 0.359 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.32 | 2.10 | 2.09 | − 1% | 0.986 |
| Pi | 0.84 | 0.80 | 0.78 | − 2% | 0.976 |
| PDE | 2.63 | 2.60 | 2.70 | 4% | 0.943 |
| PCr | 3.75 | 3.50 | 3.48 | − 1% | 0.996 |
| ATP | 2.44 | 2.21 | 2.12 | − 4% | 0.773 |
| pH | 6.99 | 6.99 | 6.99 | 0.0% | 0.992 |
| Mg | 0.10 | 0.11 | 0.10 | − 6% | 0.146 |
| l_BG | MANOVA: Wilks-Lambda: p = 0.122 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.41 | 2.30 | 2.29 | − 1% | 0.985 |
| Pi | 0.72 | 0.69 | 0.63 | − 9% | 0.488 |
| PDE | 2.42 | 2.54 | 2.46 | − 3% | 0.698 |
| PCr | 3.68 | 3.88 | 3.50 | − 10% | 0.015 |
| ATP | 2.53 | 2.32 | 2.23 | − 4% | 0.766 |
| pH | 7.00 | 7.00 | 6.99 | − 0.1% | 0.710 |
| Mg | 0.11 | 0.11 | 0.11 | 0% | 0.995 |
| l_TL | MANOVA: Wilks-Lambda: p = 0.579 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.13 | 2.13 | 2.23 | 4% | 0.631 |
| Pi | 0.63 | 0.67 | 0.59 | − 12% | 0.366 |
| PDE | 1.90 | 1.98 | 1.86 | − 6% | 0.518 |
| PCr | 3.62 | 3.91 | 3.78 | − 3% | 0.792 |
| ATP | 2.27 | 2.15 | 2.10 | − 3% | 0.870 |
| pH | 7.00 | 7.00 | 7.00 | 0.0% | 0.987 |
| Mg | 0.11 | 0.11 | 0.11 | − 1% | 0.976 |
| l_FWM | MANOVA: Wilks-Lambda: p = 0.086 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| PME | 2.22 | 2.23 | 2.19 | − 2% | 0.934 |
| Pi | 0.66 | 0.60 | 0.56 | − 8% | 0.666 |
| PDE | 2.14 | 2.27 | 2.30 | 1% | 0.958 |
| PCr | 3.33 | 3.73 | 3.47 | − 7% | 0.112 |
| ATP | 2.35 | 2.30 | 2.08 | − 10% | 0.158 |
| pH | 7.00 | 7.00 | 7.00 | − 0.1% | 0.839 |
| Mg | 0.11 | 0.11 | 0.11 | − 2% | 0.909 |
Absolute mean concentration values (in mmol/kg) and pH in left and right-sided voxels for fN, mN, and mC. Left half and right half of the FGM voxel were not evaluated separately
ΔCN indicates the relative difference of metabolite values between mN and mC: \(\Delta \mathrm{CN}=\frac{\left(\mathrm{mC}-\mathrm{mN}\right)}{\mathrm{mN}}\)
The Wilks-Lambda test reflects the overall effect of the three groups on all seven metabolite values included in the MANOVA. NAD was excluded because of too many low-quality data. Post hoc Scheffé test was used for paired comparison of groups
p values < 0.05 are marked in bold. *Values that remained significant after multiple comparisons correction
Table 2: Results of 31P MRS, comparison between hemispheres
| TH | RM MANOVA: p = 0.000* | |||||
| Paired comparison between hemispheres: | ||||||
| fN | mN | mC | ||||
| Δrl: p = 0.255 | Δrl: p = 0.000* | Δrl: p = 0.008* | ||||
| rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | |
| PME | − 1% | 0.844 | − 7% | 0.013* | − 2% | 0.382 |
| Pi | − 5% | 0.213 | − 5% | 0.168 | − 7% | 0.012* |
| PDE | − 1% | 0.771 | − 1% | 0.567 | − 2% | 0.189 |
| PCr | − 1% | 0.373 | − 8% | 0.000* | − 3% | 0.002* |
| ATP | 0% | 0.785 | − 2% | 0.178 | − 4% | 0.045 |
| pH | 0.0% | 0.722 | 0.0% | 0.838 | 0.0% | 0.782 |
| Mg | − 5% | 0.081 | 0% | 0.891 | 1% | 0.603 |
| BG | RM MANOVA: p = 0.000* | |||||
| Paired comparison between hemispheres: | ||||||
| fN | mN | mC | ||||
| Δrl: p = 0.153 | Δrl: p = 0.030 | Δrl: p = 0.116 | ||||
| rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | |
| PME | − 10% | 0.049 | − 6% | 0.134 | − 8% | 0.013 |
| Pi | − 15% | 0.063 | − 10% | 0.212 | − 16% | 0.004* |
| PDE | 10% | 0.017 | 6% | 0.304 | 7% | 0.037 |
| PCr | − 5% | 0.017 | − 17% | 0.000* | − 6% | 0.008 |
| ATP | − 2% | 0.534 | 2% | 0.441 | − 4% | 0.020 |
| pH | − 0.2% | 0.224 | − 0.1% | 0.234 | 0.0% | 0.739 |
| Mg | 0% | 0.956 | − 2% | 0.680 | − 2% | 0.391 |
| TL | RM MANOVA: p = 0.163 | |||||
| Paired comparison between hemispheres: | ||||||
| fN | mN | mC | ||||
| Δrl: p = 0.917 | Δrl: p = 0.192 | Δrl: p = 0.238 | ||||
| rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | |
| PME | − 4% | 0.422 | 2% | 0.677 | 1% | 0.673 |
| Pi | − 5% | 0.627 | − 12% | 0.156 | − 10% | 0.080 |
| PDE | 7% | 0.177 | 19% | 0.007* | 18% | 0.000* |
| PCr | − 3% | 0.331 | − 3% | 0.409 | 0% | 0.920 |
| ATP | 1% | 0.647 | 4% | 0.300 | − 1% | 0.481 |
| pH | 0.1% | 0.456 | 0.0% | 0.906 | − 0.1% | 0.272 |
| Mg | 9% | 0.093 | 3% | 0.646 | 1% | 0.890 |
| FWM | RM MANOVA: p = 0.011* | |||||
| Paired comparison between hemispheres: | ||||||
| fN | mN | mC | ||||
| Δrl: p = 0.149 | Δrl: p = 0.058 | Δrl: p = 0.220 | ||||
| rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | rel. diff. | p value (t-test) | |
| PME | 2% | 0.540 | 5% | 0.462 | 5% | 0.172 |
| Pi | − 7% | 0.628 | 5% | 0.656 | − 3% | 0.739 |
| PDE | 20% | 0.006* | 19% | 0.001* | 6% | 0.223 |
| PCr | 2% | 0.655 | − 7% | 0.007* | − 2% | 0.123 |
| ATP | 11% | 0.023 | 7% | 0.030 | 9% | 0.008* |
| pH | − 0.2% | 0.581 | − 0.2% | 0.156 | − 0.2% | 0.054 |
| Mg | − 3% | 0.428 | − 3% | 0.380 | − 1% | 0.903 |
Difference of mean concentration values and pH between right and left hemisphere for fN, mN, and mC. Relative differences Δrl were calculated as \(\Delta \mathrm{rl}=\frac{\left(\mathrm{right}-\mathrm{left}\right)}{\mathrm{left}}\). For the FGM voxel, no side difference could be determined because the left half and right half of the voxel were not evaluated separately. The RM MANOVA reflects the overall hemispheric effect of the three groups on all seven metabolite values included in the MANOVA. NAD was excluded because of too many low-quality data. Paired RM MANOVA was used to investigate hemispheric effects for each group separately. Hemispheric differences for single metabolites were evaluated by Student’s paired t-test. p values < 0.05 are marked in bold. *Values that remained significant after multiple comparison correction



Table 3: Results of 1H MRS, comparison between groups
| FGM | MANOVA: Wilks-Lambda: p = 0.047 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| tNAA | 22.6 | 22.8 | 22.1 | − 3% | 0.759 |
| tCr | 15.5 | 16.0 | 16.2 | 2% | 0.984 |
| tCho | 4.1 | 4.5 | 4.6 | 1% | 0.973 |
| r_TL | MANOVA: Wilks-Lambda: p = 0.480 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| tNAA | 16.1 | 15.9 | 15.2 | − 4% | 0.846 |
| tCr | 10.3 | 10.3 | 10.2 | − 2% | 0.950 |
| tCho | 3.0 | 3.3 | 3.1 | − 5% | 0.978 |
| r_TH | MANOVA: Wilks-Lambda: p = 0.455 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| tNAA | 14.0 | 14.8 | 13.4 | − 10% | 0.500 |
| tCr | 8.5 | 9.2 | 7.9 | − 14% | 0.712 |
| tCho | 2.4 | 2.7 | 2.5 | − 6% | 0.610 |
| r_FWM | MANOVA: Wilks-Lambda: p = 0.593 | ||||
| Post hoc Scheffé test ΔCN | |||||
| fN | mN | mC | rel. diff. | p value | |
| tNAA | 14.7 | 15.7 | 14.7 | − 6% | 0.272 |
| tCr | 8.2 | 8.8 | 8.6 | − 2% | 0.997 |
| tCho | 2.7 | 2.9 | 2.8 | − 2% | 0.962 |
Absolute mean concentration values (in mmol/kg) for fN, mN, and mC
ΔCN indicates the relative group difference of metabolite values between male non-consumers (mN) and male cannabis-consumers (mC):\(\Delta \mathrm{CN}=\frac{\left(\mathrm{mC}-\mathrm{mN}\right)}{\mathrm{mN}}\)
The Wilks-Lambda test reflects the overall effect of the three groups on all three metabolites included in the MANOVA. Post hoc Scheffé test was used for paired comparison of groups. p values < 0.05 are marked in bold

Phosphorus MRS results
Statistical analysis with MANOVA showed overall significant group differences of metabolic values only in the r_TH (p=0.014) and r_TL (p = 0.034) (Table 1). These results were no longer significant after Bonferroni correction. For some individual metabolites, concentration differences were found by post hoc tests comparing mN with fN and mN with mC, respectively.
Differences between male non-users (mN) and female non-users (fN): Δmf
ATP levels in male non-users tended to be lower than in females in nearly all regions, most noticeable in the thalamus (Δmf = − 11% in r_TH and − 10% in l_TH) (Fig. 2) where Δmf was negative for all other metabolites, too. MANOVA revealed significant differences for PME (p = 0.014) and PCr (p = 0.018) in the right thalamus, although this result did not survive Bonferroni correction.
Males tended to have lower PME mainly in the thalamus but higher PDE values in most regions. Pi was lower in the frontal lobe; pH tended to be slightly higher in males, particularly in the right thalamus and the right FWM. The higher PCr value for males (Δmf = + 12%) in the l_FWM (p = 0.018) was no longer significant after Bonferroni correction.
Differences between male cannabis users (mC) and male non-users (mN): ΔCN
Cannabis users had consistently lower ATP and Pi levels in all regions (Fig. 3). PDE values tended to be lower in mC except for the thalamus. pH levels were slightly reduced in most regions. None of these differences reached statistical significance. In the mC group, PCr was lower in the left BG (ΔCN = − 10 %, p = 0.015) and in the left FWM (− 7%), but this result was no longer significant after Bonferroni correction.
Hemispheric differences
Statistical analysis with RM MANOVA showed overall hemispheric asymmetries of metabolic values in TH (p = 0.000), BG (p = 0.000), and the FWM (p = 0.011) (Table 2) remaining significant even after Bonferroni correction. Comparing hemispheres in each group separately, significant overall side effects were found only in the thalamus for mN (p = 0.000) and mC (p = 0.008), where Δrl was negative for all metabolites (Fig. 4).
Looking at single metabolites, PCr and Pi values appeared overall lower on the right side compared to the left. ATP concentration was significantly higher in the right FWM than in the left. PME tended to be lower on the right side in TH and BG, whereas PDE was higher in the right hemisphere except for the thalamus.
In total, hemispheric differences of comparable extent could be detected in all three groups. For PCr, however, the side effect was strikingly larger in the mN group than in fN and mC, most notably in the basal ganglia (Δrl = − 17% for mN).
Proton MRS results
Statistical analysis with MANOVA revealed overall significant group differences of proton metabolic values only in FGM (p = 0.047) (Table 3), but this result remained no longer significant after Bonferroni correction.
Differences between male non-users (mN) and female non-users (fN): Δmf
Among non-users, males tended to have overall higher levels of proton metabolites than females (Fig. 5a), most pronounced for Cho (Δmf = + 11%) in the FGM.
Differences between male cannabis users (mC) and male non-users (mN): ΔCN
Compared to mN, proton metabolite data of mC showed a tendency of overall slightly lower concentrations (Fig. 5b) without reaching significance.
Overview of the results
In summary, we could not find statistically significant differences of metabolite concentrations in the brain of male cannabis-users compared to male non-users, although the data showed some tendencies.
mC showed a reduction in ATP (− 3 to − 12%) and Pi (− 2 to − 15%) in all evaluated regions. PCr concentrations were reduced only in the left BG, the left TH, and in the left FWM. For proton metabolites (tNAA, tCr, and tCho), mC tended to have slightly lower values than mN.
Some differences of metabolite values could also be detected between male and female non-users. An overall lower ATP concentration (− 2 to − 11%) was observed in mN compared to fN. In the TH region, males had generally lower metabolite concentrations than females. Proton metabolites tended to be higher in males.
Hemispheric comparison revealed statistically significant asymmetries of phosphorus metabolite values between right and left. PCr and Pi had a generally lower level in the right hemisphere, most strikingly in the BG (up to − 17%). We could not find consistent discrepancies of lateralization between the groups except for PCr which exhibited a much larger asymmetry in male non-users than in fN and mC.
Discussion
A review of the MRS literature concerning cannabis abuse clearly shows the paucity of data in this field [ref. 53–ref. 58]. To date, few 1H MRS studies characterizing proton neurometabolite concentrations in cannabis users have been published, but 31P MRS data are completely lacking until now. This study represents the first attempt to combine 1H MRS and 31P MRS in order to evaluate and compare neurometabolism in young cannabis users by performing an absolute quantification of several metabolites in different anatomic regions of the brain.
Proton spectroscopy studies dealing with cannabis consumption that have been published so far have focused on regions such as the frontal lobe, basal ganglia, hippocampus, and temporal lobe. Reduced NAA is the most frequently observed finding in cannabis users [ref. 59–ref. 62]. Particularly in the youngest subjects, reduced NAA levels were detected in frontal lobe regions, including the dorsolateral prefrontal cortex, anterior cingulate gyrus, inferior frontal gyrus, and midfrontal gray matter. Greater amount of cannabis use was associated with lower NAA and lower Cho. These results are confirmed by our study which found slightly (but not significantly) reduced NAA and Cho levels in all examined regions in the mC group.
N-acetylaspartate (NAA) is the second-most-concentrated molecule in the brain after the amino acid glutamate; its physiological function though still remains not absolutely clear [ref. 63]. NAA is detectable not only in neurons in the adult brain [ref. 64] but also in oligodendrocytes and myelin [ref. 65].
As a contributor to energy production from the amino acid glutamate, NAA correlates with the integrity of neuronal mitochondrial function [ref. 66]. Reduction of NAA concentration in the brain of cannabis users as observed both in previous and in our study might reflect neurotoxic effects of cannabis compromising neural viability.
Choline has many functions within humans and other organisms with the key feature of serving as a synthetic precursor for phospholipids that form cell membranes, the neurotransmitter acetylcholine, and trimethylglycine. Lower Cho refers to a reduced membrane turnover or increased cellular/neuronal senescence. Subsequently, lowered acetylcholine concentrations interfere with neuronal integrity, metabolism, cognition, consciousness [ref. 67] and are a predisposing factor in neurodegenerative illnesses, e.g., Alzheimer’s disease [ref. 68–ref. 70].
In our study, we did not find a significant impact of sex on proton metabolite concentrations, but the data showed some tendencies in several regions of the brain. Males (mN) tended to have higher levels than females (fN), especially for creatine and choline. These results could be interpreted in the context of well-known sex differences in brain function and structure [ref. 71]. In addition, metabolic effects of menstrual cycle have to be considered [ref. 72, ref. 73].
The most interesting result of our 31P MRS measurements is the consistent, but not statistically significant, trend to a reduction of ATP and Pi levels in mC compared to mN. PDE values were decreased mainly in the frontal and in the temporal lobe. Lowering of PCr was observed in the left part of BG, TL, and FWM. ATP, which is provided by oxidative chain reactions on the inner mitochondrial membrane, is essential for the cellular energy supply, especially for brain neurons. PCr serves as a cellular energy reservoir which can quickly provide ATP through hydrolysis. Depletion of ATP and Pi was observed in mC compared to mN which could probably point to an energy shortage in neurons, axons, and the neuroglial cells. As PDE mainly represents phospholipid breakdown products [ref. 74], reduced PDE levels as found in the frontal and in the temporal lobe could indicate lower membrane turnover, probably as a result of disturbed phospholipid generation rather than accelerated phospholipid degradation [ref. 75, ref. 76]. Cannabis-induced metabolic changes in the TL are of particular interest regarding auditory perception and language processing. The temporal lobe includes many important functions, such as the primary auditory cortex and Wernicke area which represent an integrated part of the speech recognition and speech production; the concrete function is still seen controversial [ref. 77]. The reduction of metabolite concentrations (except PME) observed in the TL is consistent with the fMRI study of Winton-Brown et al. 2011 [ref. 78] which found an attenuation of temporal auditory activation after administration of THC. Whereas their study showed an increase in psychotic symptoms associated with the attenuation of temporal activation, there were no signs of psychosis in our subject group.
Furthermore, a slight tendency of reduced pH values could be determined in the FGM, FWM, TL, and BG similar to decreased pH values reported in the frontal lobe of patients with bipolar disorder [ref. 79]. Reduction of pH could be the result of the above-mentioned energy shortage that leads to an increased anaerobic glycolysis with elevated lactate levels and reduced pH value. As we did not detect lactate in our proton MRS measurements, these presumably increased lactate levels are still very low.
In summary, our 31P MRS results can be interpreted—with all due caution—as an indication of reduced energy supply and decreased membrane turnover particularly in the frontal lobe and in the BG of cannabis users. As shown by several studies, the frontal lobe is an important part of the neuronal network responsible for social function [ref. 80], cognitive skills [ref. 81, ref. 82], and general intelligence [ref. 83]. The BG are associated with the control of movements but also a variety of cognitive and affective functions [ref. 84–ref. 86]. In conclusion, our findings might help to understand the negative impact of cannabis use on a variety of brain functions observed in long-term cannabis users. In general, the obtained 31P MRS results correlate with the results of the FDG PET findings [ref. 28] which showed a decreased glucose uptake in several brain regions of young cannabis users.
Structural T2-weighted images did not reveal any visible correlate to the metabolic changes we found in cannabis users, but several structural brain changes on cellular levels have been found in other studies, for example, the impact of cannabis use on white matter integrity [ref. 87], corpus callosum [ref. 88], gray matter density [ref. 89], and brain tissue composition [ref. 16]. Microscopy studies analyzing brain tissue of cannabis users are not available up to now. Thus, more MRI studies are needed to determine whether brain lesions might occur in the long term.
The analysis of sex influence on phosphorus metabolites yielded an inhomogeneous pattern. Major effects were found in the thalamus where mN exhibited overall lower concentrations than fN. ATP values were lower in males than in females in all examined regions. As explained above in the case of proton metabolites, sex-related differences may be partially related to hormonal conditions.
In all of our three subject groups, we found significant asymmetries of phosphorus metabolite concentrations between right and left hemisphere. In general, concentrations of PCr, Pi, and PME are higher on the left side, while PDE levels are lower, indicating an intensified energy metabolism and an elevated rate of membrane synthesis in the left hemisphere compared to the right. These side differences can be explained in the framework of functional and structural lateralization of the brain. As only right-handers were included in the study, we can assume that their left hemisphere is dominant. Complex functions like the control of behavioral structures, movement, language, and cognition are primarily located in the dominant hemisphere [ref. 90–ref. 92] potentially resulting in an asymmetric distribution of energy and membrane metabolism between the hemispheres, in accordance with our findings.
It is worth noting that the extent of metabolic asymmetry is about the same in all three analyzed subject groups except for PCr. Side differences of PCr concentration are fairly low in fN as well as in mC. In mN, however, the level of PCr was detected to be much lower on the right side than on the (dominant) left, especially in the BG. Behavioral studies have shown reduced left-hemispheric language dominance in schizophrenia as well as in healthy schizotypal subjects [ref. 93]. As cannabis use is considered as a risk factor for the development of psychosis, it may also influence the extent of lateralization for specific metabolites. Thus, the reduced asymmetry of PCr values in mC compared to mN might be interpreted in this context.
In conclusion, combined 1H/31P-MRS showed a trend towards decreased concentrations of Pi, ATP, and PCr in the frontal lobe region, as well as the right and left basal ganglia in young cannabis users compared to non-cannabis users. The results suggest that functional brain disorders observed in long-term cannabis users might be caused by an impairment of the energy metabolism of the brain, interfering neuronal integrity and viability, cognition, motoric, and sensual perception. Some of the results indicate that this impact on brain metabolism might accelerate neuronal senescence and subsequently could be a predisposing factor for neurodegenerative diseases. The extent of the observed metabolite concentration differences between the groups did not reach the level of significance; only the hemispheric asymmetric effects were statistically significant. Thus, more 31P-MRS long-term studies are required for verification.
Limitations
The results presented should be interpreted with caution. Due to imperfect magnetic field homogeneity and relatively small size of ROIs, the SNR of spectroscopic signals is fairly low in some regions of the brain, particularly for those metabolites that always produce relatively small peaks such as PME, PDE, and Pi. Moreover, there was a lack of homogeneity of the mC subject group caused by a wide range of cannabis use. This might be one of the reasons that metabolic differences between cannabis users and non-users did not reach the level of significance. In contrast, metabolic differences between the right and left hemisphere of the brain could be established in all three groups with high significance, because paired comparison within groups is less susceptible to inter-subject variations.
Another limitation concerning our 1H MRS data is caused by the use of a long echo time TE, as we originally focused on the main metabolites NAA, Cho, and Cr. Choosing a value of TE = 135 ms impeded the detection of other relevant metabolites such as glutamine and glutamate which seem to be important in cannabis use, according to recent publication [ref. 55].
One additional factor that needs to be considered in our study is the possibly confounding role of smoking. Nicotine consumption (package years) was not balanced between the groups but strongly associated with cannabis use (supplementary_Demographics: suppl_table 1). So, we cannot exclude that nicotine consumption could have contributed to the difference of metabolite values between mN and mC found in this study. On the other hand, to our knowledge, there is no comparable 31P MRS study investigating the influence of nicotine on brain metabolism in young adults but only 1H MRS studies focusing mainly on the anterior cingulate cortex with reported inconsistencies in the findings [ref. 94]. Furthermore, the half-life of nicotine in the brain is approximately 1 to 2 h [ref. 95–ref. 97] temporally restricting the effect on cerebral blood flow and metabolism. In contrast, THC is expected to have a long-term impact on brain metabolism due to its very long half-life of 5–13 days [ref. 98]. Moreover, nicotine consumption in the mC group was fairly low with a median of 1.2 py (supplementary_Demographics: suppl_table 1). In order to reduce the potential nicotine effect, subjects were instructed to abstain from smoking on examination day; noncompliance led to exclusion from the study.
As well as moving to enhanced techniques of investigation there is also the need for standardization of the populations being studied and better metabolite quantification. Future studies using 1H MRS should definitely be based on short echo times in order to extend the range of detectable metabolites. Finally, the constantly growing use of cannabis and worries about school performance in young users demand further research in this field.
Supplementary Materials
References
- 1.(UNODC) UNOoDaCWorld drugs report2018ViennaUnited Nations Publications
- DS Hasin, TD Saha, BT Kerridge, RB Goldstein, SP Chou, H Zhang. Prevalence of marijuana use disorders in the United States between 2001-2002 and 2012-2013. JAMA Psychiatry, 2015. [DOI | PubMed]
- M Cerdá, M Wall, T Feng, KM Keyes, A Sarvet, J Schulenberg. Association of state recreational marijuana laws with adolescent marijuana use. JAMA Pediatr, 2017. [DOI | PubMed]
- M Glass, M Dragunow, RL Faull. Cannabinoid receptors in the human brain: a detailed anatomical and quantitative autoradiographic study in the fetal, neonatal and adult human brain. Neuroscience, 1997. [DOI | PubMed]
- M Herkenham, AB Lynn, MR Johnson, LS Melvin, BR de Costa, KC Rice. Characterization and localization of cannabinoid receptors in rat brain: a quantitative in vitro autoradiographic study. J Neurosci, 1991. [DOI | PubMed]
- GE Terry, J Hirvonen, JS Liow, SS Zoghbi, R Gladding, JT Tauscher. Imaging and quantitation of cannabinoid CB1 receptors in human and monkey brains using (18)F-labeled inverse agonist radioligands. J Nucl Med, 2010. [DOI | PubMed]
- MH Meier, A Caspi, A Ambler, H Harrington, R Houts, RS Keefe. Persistent cannabis users show neuropsychological decline from childhood to midlife. Proc Natl Acad Sci U S A, 2012. [DOI | PubMed]
- R Gonzalez, JM Swanson. Long-term effects of adolescent-onset and persistent use of cannabis. Proc Natl Acad Sci U S A, 2012. [DOI | PubMed]
- N Solowij, KA Jones, ME Rozman, SM Davis, J Ciarrochi, PC Heaven. Verbal learning and memory in adolescent cannabis users, alcohol users and non-users. Psychopharmacology, 2011. [DOI | PubMed]
- CT Whitlow, A Liguori, LB Livengood, SL Hart, BJ Mussat-Whitlow, CM Lamborn. Long-term heavy marijuana users make costly decisions on a gambling task. Drug Alcohol Depend, 2004. [DOI | PubMed]
- HG Pope, D Yurgelun-Todd. The residual cognitive effects of heavy marijuana use in college students. Jama, 1996. [DOI | PubMed]
- RA Battisti, S Roodenrys, SJ Johnstone, C Respondek, DF Hermens, N Solowij. Chronic use of cannabis and poor neural efficiency in verbal memory ability. Psychopharmacology, 2010. [DOI | PubMed]
- RD Crean, NA Crane, BJ Mason. An evidence based review of acute and long-term effects of cannabis use on executive cognitive functions. J Addict Med, 2011. [DOI | PubMed]
- HG Pope, AJ Gruber, JI Hudson, MA Huestis, D Yurgelun-Todd. Neuropsychological performance in long-term cannabis users. Arch Gen Psychiatry, 2001. [DOI | PubMed]
- RI Block, DS O’Leary, RD Hichwa, JC Augustinack, LL Boles Ponto, MM Ghoneim. Effects of frequent marijuana use on memory-related regional cerebral blood flow. Pharmacol Biochem Behav, 2002. [DOI | PubMed]
- JA Matochik, DA Eldreth, JL Cadet, KI Bolla. Altered brain tissue composition in heavy marijuana users. Drug Alcohol Depend, 2005. [DOI | PubMed]
- L Chang, R Yakupov, C Cloak, T Ernst. Marijuana use is associated with a reorganized visual-attention network and cerebellar hypoactivation. Brain, 2006. [DOI | PubMed]
- SA Gruber, DA Yurgelun-Todd. Neuroimaging of marijuana smokers during inhibitory processing: a pilot investigation. Brain Res Cogn Brain Res, 2005. [DOI | PubMed]
- M Bittšanský, D Výbohová, D Dobrota. Proton magnetic resonance spectroscopy and its diagnostically important metabolites in the brain. Gen Physiol Biophys, 2012. [DOI | PubMed]
- JR Moffett, B Ross, P Arun, CN Madhavarao, AM Namboodiri. N-Acetylaspartate in the CNS: from neurodiagnostics to neurobiology. Prog Neurobiol, 2007. [DOI | PubMed]
- M Wyss, R Kaddurah-Daouk. Creatine and creatinine metabolism. Physiol Rev, 2000. [DOI | PubMed]
- RH Andres, AD Ducray, U Schlattner, T Wallimann, HR Widmer. Functions and effects of creatine in the central nervous system. Brain Res Bull, 2008. [DOI | PubMed]
- JS Wasser, L Vogel, SS Guthrie, N Stolowich, M Chari. 31P-NMR determinations of cytosolic phosphodiesters in turtle hearts. Comp Biochem Physiol A Physiol, 1997. [DOI | PubMed]
- T Kato, S Takahashi, T Shioiri, T Inubushi. Brain phosphorous metabolism in depressive disorders detected by phosphorus-31 magnetic resonance spectroscopy. J Affect Disord, 1992. [DOI | PubMed]
- KE Conley, AS Ali, B Flores, SA Jubrias, EG Shankland. Mitochondrial NAD(P)H in vivo: identifying natural indicators of oxidative phosphorylation in the (31)P magnetic resonance spectrum. Front Physiol, 2016. [DOI | PubMed]
- H Hamakawa, J Murashita, N Yamada, T Inubushi, N Kato, T Kato. Reduced intracellular pH in the basal ganglia and whole brain measured by 31P-MRS in bipolar disorder. Psychiatry Clin Neurosci, 2004. [DOI | PubMed]
- S Iotti, E Malucelli. In vivo assessment of Mg2+ in human brain and skeletal muscle by 31P-MRS. Magnes Res, 2008. [PubMed]
- SR Parkar, S Ramanathan, N Nair, SA Batra, SA Adarkar, AG Pandit. Cannabis dependence: effects of cannabis consumption on inter-regional cerebral metabolic relationships in an Indian population. Indian J Psychiatry, 2010. [DOI | PubMed]
- 29.Wittemann M, Brielmaier J, Rubly M, Kennel J, Werler F, Schmitgen MM, Kubera KM, Hirjak D, Wolf ND, Reith W, Wolf RC (2021) Cognition and cortical thickness in heavy cannabis users. Eur Addict Res 27(2):115–122. 10.1159/000509987
- E Von Elm, DG Altman, M Egger, SJ Pocock, PC Gøtzsche, JP Vandenbroucke. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Ann Intern Med, 2007. [DOI | PubMed]
- SW Provencher. Estimation of metabolite concentrations from localized in vivo proton NMR spectra. Magn Reson Med, 1993. [DOI | PubMed]
- T Ernst, R Kreis, B Ross. Absolute quantitation of water and metabolites in the human brain. I. Compartments and water. J Magn Reson Ser B, 1993. [DOI]
- T Mills, D Ortendahl, N Hylton, L Crooks, J Carlson, L Kaufman. Partial flip angle MR imaging. Radiology, 1987. [DOI | PubMed]
- JP Wansapura, SK Holland, RS Dunn, WS Ball. NMR relaxation times in the human brain at 3.0 tesla. Journal of Magnetic Resonance Imaging: An Official Journal of the International Society for. Magn Reson Med, 1999
- V Mlynárik, S Gruber, E Moser. Proton T 1 and T 2 relaxation times of human brain metabolites at 3 tesla. NMR Biomed, 2001. [DOI | PubMed]
- T Ethofer, I Mader, U Seeger, G Helms, M Erb, W Grodd. Comparison of longitudinal metabolite relaxation times in different regions of the human brain at 1.5 and 3 tesla. Magn Reson Med, 2003. [DOI | PubMed]
- F Träber, W Block, R Lamerichs, J Gieseke, HH Schild. 1H metabolite relaxation times at 3.0 tesla: measurements of T1 and T2 values in normal brain and determination of regional differences in transverse relaxation. J Magn Reson Imaging, 2004. [DOI | PubMed]
- H Lu, LM Nagae-Poetscher, X Golay, D Lin, M Pomper, PC Van Zijl. Routine clinical brain MRI sequences for use at 3.0 tesla. Journal of Magnetic Resonance Imaging: An Official Journal of the International Society for. Magn Reson Med, 2005
- GJ Stanisz, EE Odrobina, J Pun, M Escaravage, SJ Graham, MJ Bronskill. T1, T2 relaxation and magnetization transfer in tissue at 3T. Magn Reson Med, 2005. [DOI | PubMed]
- C Choi, NJ Coupland, PP Bhardwaj, S Kalra, CA Casault, K Reid. T2 measurement and quantification of glutamate in human brain in vivo. Magn Reson Med, 2006. [DOI | PubMed]
- SY Tsai, S Posse, YR Lin, CW Ko, R Otazo, HW Chung. Fast mapping of the T2 relaxation time of cerebral metabolites using proton echo-planar spectroscopic imaging (PEPSI). Magn Reson Med, 2007. [DOI | PubMed]
- W Zaaraoui, L Fleysher, R Fleysher, S Liu, BJ Soher, O Gonen. Human brain-structure resolved T2 relaxation times of proton metabolites at 3 tesla. Magn Reson Med, 2007. [DOI | PubMed]
- Y Li, R Srinivasan, H Ratiney, Y Lu, SM Chang, SJ Nelson. Comparison of T1 and T2 metabolite relaxation times in glioma and normal brain at 3T. J Magn Reson Imaging, 2008. [DOI | PubMed]
- II Kirov, S Liu, R Fleysher, L Fleysher, JS Babb, J Herbert. Brain metabolite proton T2 mapping at 3.0 T in relapsing-remitting multiple sclerosis. Radiology., 2010. [DOI | PubMed]
- SK Ganji, A Banerjee, AM Patel, YD Zhao, IE Dimitrov, JD Browning. T2 measurement of J-coupled metabolites in the human brain at 3T. NMR Biomed, 2012. [DOI | PubMed]
- JZ Bojorquez, S Bricq, C Acquitter, F Brunotte, PM Walker, A Lalande. What are normal relaxation times of tissues at 3 T?. Magn Reson Imaging, 2017. [DOI | PubMed]
- PO Wyss, C Bianchini, M Scheidegger, IA Giapitzakis, A Hock, A Fuchs. In vivo estimation of transverse relaxation time constant (T2) of 17 human brain metabolites at 3T. Magn Reson Med, 2018. [DOI | PubMed]
- L Vanhamme, A van den Boogaart, S Van Huffel. Improved method for accurate and efficient quantification of MRS data with use of prior knowledge. J Magn Reson, 1997. [DOI | PubMed]
- T Michaelis, K Merboldt, H Bruhn, W Hänicke, J Frahm. Absolute concentrations of metabolites in the adult human brain in vivo: quantification of localized proton MR spectra. Radiology, 1993. [DOI | PubMed]
- S Smith, P Martin, J Davies, R Edwards, A Stevens. The assessment of treatment response in non-Hodgkin’s lymphoma by image guided 31 P magnetic resonance spectroscopy. Br J Cancer, 1990. [DOI | PubMed]
- M Rata, SL Giles, NM deSouza, MO Leach, GS Payne. Comparison of three reference methods for the measurement of intracellular pH using 31P MRS in healthy volunteers and patients with lymphoma. NMR Biomed, 2014. [DOI | PubMed]
- IR Lanza, S Bhagra, KS Nair, JD Port. Measurement of human skeletal muscle oxidative capacity by 31P-MR spectroscopy: a cross-validation with in vitro measurements. J Magn Reson Imaging, 2011. [DOI | PubMed]
- 53.Sneider JT, Mashhoon Y, Silveri MM (2013) A review of magnetic resonance spectroscopy studies in marijuana using adolescents and adults. J Addict Res Ther Suppl 4. 10.4172/2155-6105.S4-010
- SD Newman, H Cheng, A Schnakenberg Martin, U Dydak, S Dharmadhikari, W Hetrick. An investigation of neurochemical changes in chronic cannabis users. Front Hum Neurosci, 2019. [DOI | PubMed]
- JJ Watts, R Garani, T Da Silva, N Lalang, S Chavez, R Mizrahi. Evidence that cannabis exposure, abuse, and dependence are related to glutamate metabolism and glial function in the anterior cingulate cortex: a (1)H-magnetic resonance spectroscopy study. Front Psychiatry, 2020. [DOI | PubMed]
- L Cupo, E Plitman, E Guma, MM Chakravarty. A systematic review of neuroimaging and acute cannabis exposure in age-of-risk for psychosis. Transl Psychiatry, 2021. [DOI | PubMed]
- CM Pretzsch, J Freyberg, B Voinescu, D Lythgoe, J Horder, MA Mendez. Effects of cannabidiol on brain excitation and inhibition systems; a randomised placebo-controlled single dose trial during magnetic resonance spectroscopy in adults with and without autism spectrum disorder. Neuropsychopharmacology., 2019. [DOI | PubMed]
- G Zoccatelli, F Alessandrini, C Rimondo, A Beltramello, G Serpelloni, EF Ciceri. Magnetic resonance spectroscopy in adolescent cannabis users: metabolites in the anterior cingulate cortex reflects individual differences in personality traits and can affect rehabilitation compliance. Neurol India, 2020. [DOI | PubMed]
- L Chang, C Cloak, R Yakupov, T Ernst. Combined and independent effects of chronic marijuana use and HIV on brain metabolites. J NeuroImmune Pharmacol, 2006. [DOI | PubMed]
- D Hermann, A Sartorius, H Welzel, S Walter, G Skopp, G Ende. Dorsolateral prefrontal cortex N-acetylaspartate/total creatine (NAA/tCr) loss in male recreational cannabis users. Biol Psychiatry, 2007. [DOI | PubMed]
- AP Prescot, AE Locatelli, PF Renshaw, DA Yurgelun-Todd. Neurochemical alterations in adolescent chronic marijuana smokers: a proton MRS study. Neuroimage., 2011. [DOI | PubMed]
- YH Sung, PD Carey, DJ Stein, HL Ferrett, BS Spottiswoode, PF Renshaw. Decreased frontal N-acetylaspartate levels in adolescents concurrently using both methamphetamine and marijuana. Behav Brain Res, 2013. [DOI | PubMed]
- G von Jonquieres, ZHT Spencer, BD Rowlands, CB Klugmann, A Bongers, AE Harasta. Uncoupling N-acetylaspartate from brain pathology: implications for Canavan disease gene therapy. Acta Neuropathol, 2018. [DOI | PubMed]
- ML Simmons, CG Frondoza, JT Coyle. Immunocytochemical localization of N-acetyl-aspartate with monoclonal antibodies. Neuroscience., 1991. [DOI | PubMed]
- K Nordengen, C Heuser, JE Rinholm, R Matalon, V Gundersen. Localisation of N-acetylaspartate in oligodendrocytes/myelin. Brain Struct Funct, 2015. [DOI | PubMed]
- G Paslakis, F Träber, J Roberz, W Block, F Jessen. N-acetyl-aspartate (NAA) as a correlate of pharmacological treatment in psychiatric disorders: a systematic review. Eur Neuropsychopharmacol, 2014. [DOI | PubMed]
- NJ Woolf. Acetylcholine, cognition, and consciousness. J Mol Neurosci, 2006. [DOI | PubMed]
- H Soreq. Checks and balances on cholinergic signaling in brain and body function. Trends Neurosci, 2015. [DOI | PubMed]
- EJ Cohen, E Quarta, G Fulgenzi, D Minciacchi. Acetylcholine, GABA and neuronal networks: a working hypothesis for compensations in the dystrophic brain. Brain Res Bull, 2015. [DOI | PubMed]
- BM Cohen, PF Renshaw, AL Stoll, RJ Wurtman, D Yurgelun-Todd, SM Babb. Decreased brain choline uptake in older adults. An in vivo proton magnetic resonance spectroscopy study. Jama., 1995. [DOI | PubMed]
- KP Cosgrove, CM Mazure, JK Staley. Evolving knowledge of sex differences in brain structure, function, and chemistry. Biol Psychiatry, 2007. [DOI | PubMed]
- LR Baxter, JC Mazziotta, ME Phelps, CE Selin, BH Guze, L Fairbanks. Cerebral glucose metabolic rates in normal human females versus normal males. Psychiatry Res, 1987. [DOI | PubMed]
- TC Hsieh, WY Lin, HJ Ding, SS Sun, YC Wu, KY Yen. Sex-and age-related differences in brain FDG metabolism of healthy adults: an SPM analysis. J Neuroimaging, 2012. [DOI | PubMed]
- JW Pettegrew, SJ Kopp, NJ Minshew, T Glonek, JM Feliksik, JP Tow. 31P nuclear magnetic resonance studies of phosphoglyceride metabolism in developing and degenerating brain: preliminary observations. J Neuropathol Exp Neurol, 1987. [DOI | PubMed]
- MS van der Knaap, J van der Grond, PC van Rijen, JA Faber, J Valk, K Willemse. Age-dependent changes in localized proton and phosphorus MR spectroscopy of the brain. Radiology, 1990. [DOI | PubMed]
- 76.Lee J-H, Komoroski R, Chu W-J, Dudley J (2012) Methods and applications of phosphorus NMR spectroscopy In Vivo. 10.1016/B978-0-12-397018-3.00003-X
- JR Binder. Current controversies on Wernicke’s area and its role in language. Curr Neurol Neurosci Rep, 2017. [DOI | PubMed]
- TT Winton-Brown, P Allen, S Bhattacharyya, SJ Borgwardt, P Fusar-Poli, JA Crippa. Modulation of auditory and visual processing by delta-9-tetrahydrocannabinol and cannabidiol: an FMRI study. Neuropsychopharmacology, 2011. [DOI | PubMed]
- T Kato, J Murashita, A Kamiya, T Shioiri, N Kato, T Inubushi. Decreased brain intracellular pH measured by 31P-MRS in bipolar disorder: a confirmation in drug-free patients and correlation with white matter hyperintensity. Eur Arch Psychiatry Clin Neurosci, 1998. [DOI | PubMed]
- D Koshiyama, M Fukunaga, N Okada, K Morita, K Nemoto, F Yamashita. Role of frontal white matter and corpus callosum on social function in schizophrenia. Schizophr Res, 2018. [DOI | PubMed]
- M Ishibashi, N Kimura, Y Aso, E Matsubara. Effects of white matter lesions on brain perfusion in patients with mild cognitive impairment. Clin Neurol Neurosurg, 2018. [DOI | PubMed]
- M Daianu, MF Mendez, VG Baboyan, Y Jin, RJ Melrose, EE Jimenez. An advanced white matter tract analysis in frontotemporal dementia and early-onset Alzheimer’s disease. Brain Imaging Behav, 2016. [DOI | PubMed]
- T Ohtani, PG Nestor, S Bouix, Y Saito, T Hosokawa, M Kubicki. Medial frontal white and gray matter contributions to general intelligence. PLoS One, 2014. [DOI | PubMed]
- AC Bostan, RP Dum, PL Strick. Functional anatomy of basal ganglia circuits with the cerebral cortex and the cerebellum. Prog Neurol Surg, 2018. [DOI | PubMed]
- RM Cotterill. Cooperation of the basal ganglia, cerebellum, sensory cerebrum and hippocampus: possible implications for cognition, consciousness, intelligence and creativity. Prog Neurobiol, 2001. [DOI | PubMed]
- TC Ford, AC Hayley, LA Downey, AC Parrott. Cannabis: an overview of its adverse acute and chronic effects and its implications. Curr Drug Abuse Rev, 2017. [DOI | PubMed]
- JM Orr, CJ Paschall, MT Banich. Recreational marijuana use impacts white matter integrity and subcortical (but not cortical) morphometry. Neuroimage Clin, 2016. [DOI | PubMed]
- S Rigucci, TR Marques, M Di Forti, H Taylor, F Dell’Acqua, V Mondelli. Effect of high-potency cannabis on corpus callosum microstructure. Psychol Med, 2016. [DOI | PubMed]
- T Demirakca, A Sartorius, G Ende, N Meyer, H Welzel, G Skopp. Diminished gray matter in the hippocampus of cannabis users: possible protective effects of cannabidiol. Drug Alcohol Depend, 2011. [DOI | PubMed]
- C Helmstaedter, M Kurthen, DB Linke, CE Elger. Right hemisphere restitution of language and memory functions in right hemisphere language-dominant patients with left temporal lobe epilepsy. Brain, 1994. [DOI | PubMed]
- MC Ridding, SC Flavel. Induction of plasticity in the dominant and non-dominant motor cortices of humans. Exp Brain Res, 2006. [DOI | PubMed]
- L Zago, PY Hervé, R Genuer, A Laurent, B Mazoyer, N Tzourio-Mazoyer. Predicting hemispheric dominance for language production in healthy individuals using support vector machine. Hum Brain Mapp, 2017. [DOI | PubMed]
- DA Herzig, S Sullivan, G Lewis, R Corcoran, R Drake, J Evans. Hemispheric language asymmetry in first episode psychosis and schizotypy: the role of cannabis consumption and cognitive disorganization. Schizophr Bull, 2015. [DOI | PubMed]
- T Hellem, X Shi, G Latendresse, PF Renshaw. The utility of magnetic resonance spectroscopy for understanding substance use disorders: a systematic review of the literature. J Am Psychiatr Nurses Assoc, 2015. [DOI | PubMed]
- JE Rose, FM Behm, EC Westman, RJ Mathew, ED London, TC Hawk. PET studies of the influences of nicotine on neural systems in cigarette smokers. Am J Psychiatry, 2003. [DOI | PubMed]
- OB Paulson, I Vigdis. Cigarette smoking and cerebral blood flow in a cohort of middle-aged adults. J Cereb Blood Flow Metab, 2020. [DOI | PubMed]
- SR Dager, SD Friedman. Brain imaging and the effects of caffeine and nicotine. Ann Med, 2000. [DOI | PubMed]
- P Sharma, P Murthy, MM Bharath. Chemistry, metabolism, and toxicology of cannabis: clinical implications. Iran J Psychiatry, 2012. [PubMed]
