Residual Effects of THC via Novel Measures of Brain Perfusion and Metabolism in a Large Group of Chronic Cannabis Users.

Given the known vascular effects of cannabis, this study examined the neurophysiological factors that may affect studies of brain activity in cannabis users. We conducted a systematic evaluation in 72 h abstinent, chronic cannabis users (N=74) and nonusing controls (N=101) to determine the association between prolonged cannabis use and the following neurophysiological indicators: (1) global and regional resting cerebral blood flow (CBF), (2) oxygen extraction fraction (OEF), and (3) cerebral metabolic rate of oxygen (CMRO2). We found that cannabis users had greater global OEF and CMRO2 compared with nonusers. Regionally, we found higher CBF in the right pallidum/putamen of the cannabis users compared with nonusers. Global resting CBF and regional CBF of right superior frontal cortex correlated positively with creatinine-normalized Δ9-tetrahydrocannabinol (THC) levels. These findings demonstrate residual effects of cannabis use whereby global and regional brain metabolism are altered in those with prolonged cannabis exposure. These neurophysiological alterations should be considered in both research and clinical applications.

Introduction

Cannabis is the most commonly used illicit psychoactive drug, and therefore understanding its effects on behavior and brain is paramount. The primary psychoactive ingredient in cannabis, Δ-9-tetrahydrocannabinol (THC), relaxes arterial walls resulting in lower blood pressure and increased blood flow to tissues. In the brain, THC binds to ubiquitous cannabinoid receptors (CB1) that are present in arterial tissue (Bilfinger ) and regulate the microvascular environment via dose-dependent dilation of cerebral arterioles (Ellis ). These cerebrovascular effects are variable across brain regions (Bloom ; O'Leary ), but have been most reported in frontal regions (Mathew and Wilson, 1993). These findings indicate that neurophysiological factors may contribute to cannabis’ effects on brain activity, particularly as measured by the blood-oxygen-level-dependent (BOLD) signal.

To date, the effects of cannabis use on neurophysiological factors have primarily been measured using cerebral blood flow (CBF). In general, findings have been highly variable, perhaps because of the diverse and complex processes by which THC affects CBF including relaxation of vascular smooth muscle, alteration of partial pressure of carbon dioxide, and sympathetic stimulation (Zwillich ). THC and CBF changes also have discordant time courses (Mathew ). Neuroimaging studies have found differential effects of acute vs chronic THC exposure. Specifically, acute exposure to cannabis has been associated with increased resting CBF in various brain regions including prefrontal areas and basal ganglia (Volkow ), whereas chronic use has been associated with lower resting CBF in prefrontal and cerebellar metabolism as compared with nonusers (Volkow ; Mathew and Wilson, 1991). CBF alterations have also been shown to be THC dose dependent (Mathew ) and to normalize following abstinence (Jacobus ). Altogether, these findings suggest that prolonged exposure to THC alters CBF that affects brain function and structure. At this time, however, the literature remains sparse with regard to how prolonged THC use may affect the brain in terms of blood oxygenation extraction and metabolism beyond CBF studies. This is of high relevance given the growing number of cannabis use recreationally and medically.

Given the known cerebrovascular effects of cannabis, this study examined putative residual effects of THC on neurophysiological parameters during rest. We measured global/regional CBF, oxygen extraction fraction (OEF), and cerebral metabolic rate of oxygen (CMRO2) using several novel MRI techniques including pseudo-continuous arterial spin labeling (pCASL) MRI, T2-relaxation-under-spin-tagging (TRUST) MRI, and combined TRUST/phase contrast (PC) MRI, respectively. PC MRI estimates the total amount of blood flow delivered to the brain via bilateral internal carotid (ICA) and vertebral arteries (VA). PC MRI, however, does not provide regional CBF information. Thus, we used pCASL MRI to evaluate regional CBF. Next, we used TRUST MRI to first measure total brain’s blood oxygenation (Yv) and subsequently estimate oxygen extraction fraction and CMRO2. We investigated how the homeostasis among the aforementioned neurophysiological measures differed in 72 h abstinent, long-term cannabis users compared with nonusers. Furthermore, given the known CBF differences between males and females (Lu ) as well as sex-dependent physiological and behavioral effects of THC (Ketcherside ; Fattore and Fratta, 2010), we also explored potential group × sex interactions. Finally, we examined the association between the neurophysiological measures and patterns of cannabis use as well as cognitive functioning.

Materials and methods

Participants

A total of 175 participants, 74 cannabis users and 101 nonusers, were recruited for this study (see Table 1). Cannabis users were recruited based on self-reported history of regular cannabis use with a minimum of 5000 lifetime occasions of use in addition to daily use in the preceding 60 days before study enrollment. Verification of cannabis use was conducted via quantification of THC in urine via gas chromatography/mass spectrometry (GC/MS) normalized by urine creatinine (Cr) concentration as per the recommendation by Huestis and Cone (1998). GC/MS was conducted by Quest Diagnostics (https://www.questdiagnostics.com/). The nonusing controls were recruited based on the absence of lifetime daily cannabis use. Exclusionary criteria for all participants were: history of neurological or psychiatric conditions, MR scanning contraindications, drug abuse, regular tobacco use (>1 pack of cigarettes per month), left handedness, and non-native English speakers. As part of a larger study aimed to determine cue-elicited craving (see Filbey ), cannabis users were scanned following a verified 72 h abstinence from cannabis. Abstinence was verified by self-report and decrease in THC/Cr from baseline levels after the 72 h abstinence period. Those who were not able to comply with the abstinence requirement were excluded from this study.

Table 1

Participant Characteristics (Mean±SD)

	Nonusers	Users	P-value
N	101	74	—
Age	30.3±10.3	31.3±7.9	0.45
Gender (M/F)	50/51	46/28	0.10
Full IQ	109.7±13.8	106.7±12.0	0.14
Duration of use (years)	—	10.6±7.3	—
No. of lifetime occurrence of cannabis use	—	14 173.8±10 866.0	—
Baseline THC/Cr	—	6.0±6.5	—
3-day abstinent THC/Cr (at time of scan)	—	1.9±1.4	—
No. of hours since last cannabis use	—	75.8±4.7	—
Number of participants per scan: N (M/F)
PC MRI	98 (49/49)	67 (40/27)
TRUST MRI	93 (47/46)	66 (39/27)
pCASL MRI	94 (47/47)	64 (40/24)

Abbreviations: M/F, male/female; IQ, Intelligence Quotient via Wechsler Abbreviated Scale of Intelligence as a measure of general intellect (Wechsler, 1999); THC, Δ9-tetrahydrocannabinol; Cr, creatinine; PC MRI, phase contrast MRI to measure global cerebral blood flow; TRUST MRI, T2-relaxation-under-spin-tagging MRI to measure global venous oxygenation; pCASL MRI, pseudo-continuous arterial spin labeling MRI to measure regional cerebral blood flow.

Written informed consent was obtained from all participants in accordance with the institutional review board (IRB) of The University of Texas at Dallas and The University of Texas Southwestern Medical Center.

Measures of Cannabis Use and IQ

Duration of cannabis use, total number of lifetime occurrences, and THC/Cr levels measured at scan day (with 72 h abstinence) were collected as measures of cannabis use. We also collected the two-subtest version of the Wechsler Abbreviated Scale of Intelligence as a measure of general ability or IQ (Wechsler, 1999).

MRI Acquisition

MRI scans were performed on a 3 Tesla MR system (Philips Medical System, Best, The Netherlands). A body coil was used for radiofrequency (RF) transmission and an 8-channel head coil with parallel imaging capability was used for signal reception. We used several MRI techniques to investigate brain physiology at rest: a time-of-flight (TOF) angiogram to visualize the ICA and VA, PC MRI to measure whole-brain blood flow by positioning the imaging slab perpendicular to internal and vertebral arteries (Aslan ), TRUST MRI to assess blood oxygenation of the brain via superior sagittal sinus (Lu and Ge, 2008), and pCASL MRI sequence to measure regional blood flow of the brain (Aslan ). Finally, a high-resolution T1-weighted image was acquired as an anatomical reference.

The details of imaging parameters for each technique are as follow: the TOF angiogram imaging parameters were: TR/TE/flip angle=23 ms/3.45 ms/18°, field of view (FOV)=160 × 160 × 70 mm3, voxel size 1.0 × 1.0 × 1.5 mm3, number of slices=47, one saturation slab of 60 mm positioned above the imaging slab to suppress the venous vessels, duration 1 min and 26 s. The slice of the PC velocity MRI was oriented perpendicular to the ICA and VA and the parameters were: single slice, voxel size=0.50 × 0.50 mm2, FOV=200 × 200 mm2, TR/TE=20/7 msec, flip angle=15°, slice thickness=5 mm, maximum velocity encoding=80 cm/s, and scan duration=30 s. Global venous oxygenation, Yv, was noninvasively assessed from the superior sagittal sinus using a validated approach TRUST MRI (Lu and Ge, 2008; Lu ; Xu ). The TRUST MRI imaging parameters were: voxel size 3.44 × 3.44 × 5 mm3, TR=3000 ms, TI=1022 ms, four effective TEs: 0, 40, 80, and 160 ms, labeling thickness 100 mm, gap 22.5 mm, and scan duration 1.2 min. Imaging parameters for pCASL experiments were: single-shot gradient-echo EPI, FOV=240 × 240, matrix=80 × 80, voxel size=3 × 3 mm2, 29 slices acquired in ascending order, slice thickness=5 mm, no gap between slices, labeling duration=1650 ms, time interval between consecutive slice acquisitions=35.5 ms, TR/TE=4211/14 ms, SENSE factor 2.5, number of controls/labels=40 pairs, RF duration=0.5 ms, pause between RF pulses=0.5 ms, labeling pulse flip angle=18°, bandwidth=2.7 kHz, echo train length=35, and scan duration 5.7 min.

MRI Analyses

Global measures

Global CBF was calculated using PC MRI. Briefly, PC MRI provides a quantitative measurement of blood flow velocity in internal carotid artery and vertebral artery. The blood velocity was converted to flow rate by integrating over the cross-section of the vessels. Then, we divided the flow rate by the parenchyma intracranial mass (calculated from T1-weighted image and by assuming brain density of ρ=1.06 g/ml) and obtained the whole-brain averaged blood flow in units of ml blood/min/100 g brain (Aslan ).

A series of steps were carried out to determine total brain oxygen extraction and cerebral oxygen metabolism. First, TRUST MRI was utilized to noninvasively obtain global venous oxygenation, Yv, from the superior sagittal sinus (Liu ; Lu and Ge, 2008; Xu ). Briefly, TRUST data were processed by subtracting control and label images to yield pure venous blood signal. Then, the venous blood signals were fitted to a monoexponential function to obtain T2 that was in turn converted to Yv (ie global venous oxygenation) via a calibration plot (Lu ). Finally, total CMRO2 was then estimated from global cerebral blood flow and Yv based on Fick’s principle in units of μmol O2/100 g/min (Kety and Schmidt, 1948; Peng ; Xu ):

where CMRO and tCBF are total CMRO2 and cerebral blood flow, respectively; Ya and Yv are oxygen saturation percentage in arterial and venous blood, respectively; and Ch is a constant representing the capacity of blood to carry O2 and is well established in physiology literature (Guyton and Hall, 2005). Both Ya and Ch were estimated based on prior findings and were adjusted for age and sex (Peng ). In addition, we estimated the OEF by calculating the difference between arterial blood oxygenation and global venous oxygenation in sagittal sinus: (Ya−Yv)/Ya.

Regional measures

In terms of regional brain blood flow, we analyzed pCASL MRI data using routine processing (Aslan ). The pCASL image series were realigned to the first volume for motion correction (SPM5’s realign function, University College London, UK). ASL images with motion of >3 mm and >3° were excluded from further processing. Next, an in-house MATLAB (Mathworks, Natick, MA) program was used to calculate the difference between averaged control and label images. Then, the difference image was corrected for imaging slice delay time to yield CBF-weight image that was normalized to the Montreal Neurological Institute (MNI) brain template. Finally, the absolute CBF (aCBF) was estimated in the units of ml blood/min/100 g of brain tissue. The voxel-wise analyses were performed on the relative CBF (rCBF) maps that included dividing the aCBF spatial maps by the whole-brain absolute CBF value. In a prior publication, we showed that such technique improves the sensitivity of regional differences by reducing physiological variations (Aslan and Lu, 2010a). During voxel-based analyses (VBAs), the individual rCBF maps were spatially smoothed (with full-width half-maximum (FWHM) of 4 mm) to account for small differences in sulci/gyri location across subjects. To define a threshold for the VBA results, we used 3dClustsim (with -acf option) in AFNI (NIMH Scientific and Statistical Computing Core, Bethesda, MD) that controls for false positive activation clusters over the set of all activation clusters throughout the whole-brain volume given voxel size, whole-brain volume, and effective smoothness (inherent smoothness plus additional smoothness applied). We refer to this procedure in the Results section as cluster-level family-wise error correction (FWE corrected). The error correction was conditional on three criteria: smoothness of the voxel map, a magnitude threshold for defining clusters of contiguous voxels, and a minimum volume for each cluster. For the CBF VBA, we used p(FWE corrected) value of <0.05 that is based on an effective smoothness of 8.7 mm FWHM, a cluster-defining magnitude corresponding to the 99.5th percentile of t-statistic distribution at p<0.0002, and a minimum volume of 34 voxels (272 mm3).

Statistical Analyses

A general linear model (GLM) was applied to assess the effects of cannabis on the neurophysiological brain measures (ie, CBF, OEF, and CMRO2). For each of these dependent variables, the model included group status (user, nonuser), sex (male, female), and group × sex interaction. Finally, voxel-wise correlation analyses were conducted in the cannabis users to determine the relationship between neurophysiological indexes and duration of cannabis use, total number of lifetime occurrences, scan day THC/Cr, and IQ. A statistical threshold of p(FWE corrected)<0.05 and cluster size of >272 mm3 was applied for these analyses.

Results

Participant Characteristics

Only cannabis users who (1) tested positive for THC at baseline and (2) had a decrease in THC at the time of scan were included in this study. Of the 74 cannabis users, one had an incomplete scan and was excluded from further analyses. For each imaging scan, there were varying number of participants with high motion and were subsequently excluded from respective analyses. Table 1 summarizes the number of participants included per analysis. There were no significant differences in age, gender, and IQ (p>0.05). Of the nonusers, 73 reported alcohol use and 2 reported regular tobacco use. Of the nonusers, 63 reported alcohol use, and 21 of the users reported regular tobacco use.

Global Neurophysiological Markers

Table 2 summarizes the mean global neurophysiological measures and Table 3 lists the statistical results for the factorial analysis.

Table 2

Global Brain Measurements (Mean±SEM)

	Nonusers	Users
CBF (ml/min/100 g)
All	59.9±9.8	57.8±9.2
Males	56.8±9.5	55.2±7.7
Females	62.9±9.8	61.7±10.0
OEF (%)
All	32.5±7.6	35.0±7.3
Males	31.1±7.2	35.1±7.6
Females	33.9±7.8	34.8±7.0
CMRO2 (μmol/100 g/min)
All	154.9±27.2	162.2±25.7
Males	145.0±27.5	158.3±25.7
Females	165.1±23.0	167.8±25.3

Table 3

Statistical Test Results for Global Neurophysiological Measures Are Shown

	F-value	P-value
CBF (ml/min/100 g)
Group	F(1, 161)=0.20	0.66
Sex	F(1, 161)=18.80	<0.001
Group × sex	F(1, 161)=0.01	0.91
OEF (%)
Group	F(1, 155)=4.69	0.03
Sex	F(1, 155)=0.99	0.32
Group × sex	F(1, 155)=1.73	0.19
CMRO2 (μmol/100 g/min)
Group	F(1, 155)=4.20	0.04
Sex	F(1, 155)=12.84	<0.001
Group × sex	F(1, 155)=1.66	0.20

Abbreviations: CBF, cerebral blood flow; OEF, oxygenation extraction fraction; CMRO2, cerebral metabolic rate of oxygen.

CBF

There was no significant global CBF effect for group (F(1, 161)=0.20, p=0.66) or group × sex interaction (F(1, 161)=0.01, p=0.91). As expected, there was an effect of sex indicating higher CBF in females compared with males (F(1, 161)=18.80, p<0.001) (Lu ).

OEF

We found a significant group effect in total OEF. Users had higher OEF compared with nonusers (F(1, 155)=4.69, p=0.03). There was no significant effect of sex (F(1, 155)=0.99, p=0.32) or group × sex interaction (F(1, 155)=1.73, p=0.19) in OEF.

CMRO2

We found a significant group effect in CMRO2. Users had higher CMRO2 compared with controls (F(1, 155)=4.20, p=0.04). We also found that females have higher CMRO2 compared with males (F(1, 155)=12.84, p<0.001). There was no group × sex interaction effect in CMRO2 (F(1, 155)=1.66, p=0.20).

CBF vs OEF

We conducted a confirmatory analysis to verify the coupling between global CBF and OEF. Our results paralleled the findings of Peng of an inverse correlation in CBF and OEF in both the user (r=−0.58 and p<0.001) and nonuser group (r=−0.58 and p<0.001).

Correlation analyses

As shown in Figure 1, a significantly positive association emerged in the users between resting global CBF and THC/Cr (r=0.47 and p=0.001). No significant correlations were found between global neurophysiological markers and duration of cannabis use, total number of lifetime occurrences, or IQ.

Figure 1

Scatterplot of global cerebral blood flow (CBF) and THC/Cr in cannabis users. A significant correlation was found between CBF and THC/Cr level (r=0.47 and p=0.001).

Regional Neurophysiological Markers

Figure 2 and Table 4 summarize the results of voxel-based analyses. First, users had significantly higher blood flow in the right pallidum and putamen compared with nonusers: peak MNI (+20 −4 −2), F(1, 154)=27.05, cluster size=384 mm3 at p(FWE-corrected)<0.05, and k⩾272 mm3 (Figure 2). Second, we found a significant effect of sex: males had a higher regional CBF in the right insula compared with females whereas females had a higher CBF in the left posterior cingulate and bilateral precuneus compared with males. No significant interaction between group and sex was found.

Figure 2

Results of cerebral blood flow (CBF) voxel-based comparisons superimposed on average CBF map of all participants. (a) The cannabis users showed higher CBF in right pallidum compared with nonusers. (b) Males showed higher CBF in the right insula, whereas (c) females showed higher CBF in the right posterior cingulate cortex and bilateral precuneus (FWE cluster level corrected at p<0.05).

Table 4

Regional CBF Differences at Rest between Users and Nonusers

			MNI
Brain regions	BA	Cluster size (mm3)	X	Y	Z	F-score
Main effect of group
Users >Nonusers
Right pallidum	—	384	20	−4	−2	27.05
Main effect of sex
Males >females
Right insula	13	720	48	2	4	23.05
Females >males
Left posterior cingulate cortex	23	360	−4	−56	16	26.50
Left/right precuneus	7/19	1288	−10	−78	48	25.89
Group × sex interaction
No significant difference

FWE cluster level, corrected p<0.05, voxel threshold at p=0.0002 and K=272 mm3.

Voxel-wise correlation analysis showed a positive correlation between scan day THC/Cr level and the CBF of right superior frontal cortex at peak MNI coordinates: (+20 +48 +18), t-score of 4.45, cluster size=328 mm3, p(FWE corrected)<0.05, and cluster size ⩾272 mm3 (Figure 3). No significant correlations were found between regional neurophysiological markers and duration of cannabis use, total number of lifetime occurrences, or IQ.

Figure 3

The cannabis users showed a significantly positive association between right superior frontal cortex cerebral blood flow and THC/Cr level (p<0.05, FWE corrected).

Discussion

In this study, we evaluated neurophysiological indicators via novel measures of brain perfusion and metabolism in a large and well-characterized group of 72 h abstinent cannabis users relative to nonusers. We found greater regional CBF in the right pallidum/putamen in users compared with nonusers, in line with reported functional and structural changes in this critical region in long-term cannabis users (Filbey ). These findings of regional but not global CBF differences between users and nonusers are concordant with region-specific effects described in studies of both chronic (Bloom ; Mathew and Wilson, 1993) and acute effects of cannabis use (O'Leary ). Increases in resting CBF in basal ganglia structures have been linked to elevated salience and motivation processes that are mechanisms that contribute to compulsive use of cannabis (Volkow ). The pallidum plays a critical role as the interface between the ventral striatum and the thalamus with the ventral pallidum, specifically, mediating reward and motivation processes (Smith ). The putamen is dense with cannabinoid receptors (Miller and Walker, 1996) and has been shown to have functional and structural alterations that are associated with cannabis use (Filbey ; Wetherill ; Yip ). Because THC relaxes arterial walls and increases blood flow to tissues, our finding adds to the literature by demonstrating that the pathomechanism underlying the reported functional and structural alterations in the pallidum/putamen may be due to altered perfusion that is closely coupled with neural activity.

Notably, other studies have reported decreases rather than increases in CBF in chronic cannabis users (Martin-Santos ; Mathew and Wilson, 1991; Lundqvist ). The discrepancy with existing reports may be because of heterogeneity both within and between study samples. For example, Mathew reported dissociable CBF effects where experienced users had greater CBF whereas inexperienced users had decreased CBF following THC exposure. Length of abstinence can also affect CBF effects. Tunving ) noted changes during the course of abstinence whereby rCBF level in cannabis users was significantly (11%) lower than that of age- and sex-matched healthy controls, although (12%) increases after 9–60 days of detoxification of the CBF level was noted at follow-up. Given that our sample of cannabis users were heavy, long-term users and that data were collected following 3 days of required abstinence, this observation may be an ‘overshoot of recovery’ of neurophysiological function during the early stages of abstinence in chronic cannabis users. In other words, these findings reflect the residual effects of both long-term cannabis exposure and adaptive changes that occur after a short period of abstinence. Future studies should measure the trajectory of these effects to fully determine whether this phenomenon is indeed a recovery of typical function.

It is also important to note that in addition to chronicity and cessation from cannabis use that complicate the understanding of THC effects on CBF, individual variation in response to THC also plays a role (Volkow ). One important individual factor to consider is sex-specific effects of cannabis (see review by Ketcherside ). Although we found concordant findings of a main effect of sex on CBF and CMRO2 where females had greater CBF and CMRO2 (Lu ), we did not find a group × sex interaction in any of the brain measures. This could be due, in part, to the smaller number of females in the cannabis using group. Another possibility for the absence of a group × sex effect may be that any increased neuronal needs reflected by greater OEF/CMRO2 and CBF in the pallidum/putamen in the cannabis users relative to nonusers may be compensated for by the typically elevated global CBF and CMRO2 in females. Greater global CBF and CMRO2 in females have been suggested to reflect hormonal, muscle, and tissue content differences between sexes (Fattore and Fratta, 2010).

We found positive associations between global and regional CBF and THC metabolite levels that further support the relationship between changes in blood flow and exposure to cannabis (Mathew ; Volkow ). Alterations in the superior frontal gyrus have previously been reported in cannabis users in terms of reduced thickness (Lopez-Larson ) and enhanced activation during memory and inhibitory tasks (Smith ; Tapert ). This finding highlights the differential effects of acute and prolonged THC exposure that brings neural adaptations in long-term cannabis users to bear.

These findings have several implications. First, because CBF and brain function are closely coupled, it is likely that widely observed changes in brain response (eg, BOLD fMRI) may be partly explained by these latent cerebrovascular changes. There is emergent literature that has begun to shed light on the interplay between neural and vascular components in substance use disorders. Similar investigations have linked changes in EEG patterns with cerebral perfusion deficits in long-term cannabis users during abstinence that were interpreted as under-arousal of these systems (Herning ). Using CMRO2, long-term cocaine users have also been shown to have reduced CMRO2 relative to nonusing controls that was associated with lifetime frequency of cocaine use (Liu ).

It is worth noting that at our selected threshold, no group effects emerged in cortical areas. However, a less conservative threshold of p<0.01 and K>200 revealed that nonusers had higher regional CBF in cortical areas including right middle temporal, left middle frontal, and right superior frontal cortices compared with users. Thus, although differences were observed in other regions, our findings suggest that neurophysiological changes as a result of cannabis use are most robust in the pallidum/putamen area. Indeed, structural and functional alterations in the putamen have been reported to be a risk factor for cannabis use disorder (Wetherill ), predict development of problems related to cannabis use (Vingerhoets ), and are related to treatment success (Yip ).

Limitations and conclusions

Given the unique effects of CBF and THC on sex, it was surprising that we did not find a group × sex interaction. It is possible that our unequal sample sizes, which was because of smaller number of females in the users group, contributed to the failure to find an effect. Future studies should focus on obtaining equal sample sizes across group and sex. Although we found that CBF was associated with THC/Cr, no other association between measures of cannabis use or IQ emerged. Future studies should consider evaluating other measures of cognitive function related to cannabis use, such as inhibitory control or reward sensitivity. Finally, interpretation of these findings is limited by the cross-sectional nature of this study. Thus, we cannot infer causality in these effects. To conclude, cerebrovascular markers are promising for the assessment of the long-term effects of cannabis use on the brain. Understanding these neurophysiological mechanisms in cannabis users and their interplay with other mechanisms related to cannabis use will inform the widely reported alterations in brain structure and function in this population.

Funding and disclosure

This research was funded by NIH R01 DA030344 (to FMF) and NIH R01 MH084021 (to HL). The authors declare no conflict of interest.