Abstract
Recent trends in cannabis legalization have increased the necessity to better understand the effects of cannabis use. Animal models involving traditional cannabinoid self-administration approaches have been notoriously difficult to establish and differences in the drug used and its route of administration have limited the translational value of preclinical studies. To address this challenge in the field, we have developed a novel method of cannabis self-administration using response-contingent delivery of vaporized Δ9-tetrahydrocannabinol-rich (CANTHC) or cannabidiol-rich (CANCBD) whole-plant cannabis extracts. Male Sprague-Dawley rats were trained to nose-poke for discrete puffs of CANTHC, CANCBD, or vehicle (VEH) in daily 1 h sessions. Cannabis vapor reinforcement resulted in strong discrimination between active and inactive operanda. CANTHC maintained higher response rates under fixed ratio schedules and higher break points under progressive ratio schedules compared with CANCBD or VEH, and the number of vapor deliveries positively correlated with plasma THC concentrations. Moreover, metabolic phenotyping studies revealed alterations in locomotor activity, energy expenditure, and daily food intake that are consistent with effects in human cannabis users. Furthermore, both cannabis regimens produced ecologically relevant brain concentrations of THC and CBD and CANTHC administration decreased hippocampal CB1 receptor binding. Removal of CANTHC reinforcement (but not CANCBD) resulted in a robust extinction burst and an increase in cue-induced cannabis-seeking behavior relative to VEH. These data indicate that volitional exposure to THC-rich cannabis vapor has bona fide reinforcing properties and collectively support the utility of the vapor self-administration model for the preclinical assessment of volitional cannabis intake and cannabis-seeking behaviors.
SIGNIFICANCE STATEMENT The evolving legal landscape concerning recreational cannabis use has increased urgency to better understand its effects on the brain and behavior. Animal models are advantageous in this respect; however, current approaches typically used forced injections of synthetic cannabinoids or isolated cannabis constituents that may not capture the complex effects of volitional cannabis consumption. We have developed a novel model of cannabis self-administration using response-contingent delivery of vaporized cannabis extracts containing high concentrations of Δ9 tetrahydrocannabinol (THC) or cannabidiol. Our data indicate that THC-rich cannabis vapor has reinforcing properties that support stable rates of responding and conditioned drug-seeking behavior. This approach will be valuable for interrogating effects of cannabis and delineating neural mechanisms that give rise to aberrant cannabis-seeking behavior.
Introduction
With several states recently passing legislation allowing for the use of cannabis for recreational purposes, there is a pressing need to better understand the effects of cannabis use/misuse. Animal models are valuable in this regard because they afford precise control over extraneous variables that complicate interpretation of cross-sectional human data. However, current approaches have limited translational value (McLaughlin, 2018). Synthetic cannabinoid receptor 1 (CB1R) agonists or isolated cannabis constituents [e.g., Δ9-tetrahydrocannabinol(THC), cannabidiol (CBD)] have become the drug of choice in rodent models of cannabis exposure; even though the pharmacological properties of these compounds differ greatly from those of inhaled cannabis. THC and synthetic CB1R agonists have different pharmacological profiles (partial vs full agonists) and recruit different intracellular signaling pathways (Laprairie et al., 2014). Thus, synthetic CB1R agonists may fail to recapitulate the effects of THC, let alone the effects of cannabis.
Animal models involving access to THC alone have drawbacks as well. Over 120 unique phytocannabinoids are present in Cannabis sativa in addition to THC. Proportions of phytocannabinoids can vary considerably across cannabis chemovars and their distinct pharmacological profiles and potential to interact may give rise to their unique psychotropic effects (Fisar, 2009; Russo, 2011; Hill et al., 2012; Lewis et al., 2018). For example, CBD is a CB1R negative allosteric modulator and inhibits THC-dependent intracellular signaling and β-arrestin-2-mediated CB1R internalization (Laprairie et al., 2015, 2016). CBD attenuates THC-induced paranoia and memory impairments in humans (Englund et al., 2013) but can also increase THC serum concentrations when co-administered with THC (Greene et al., 2018). Therefore, THC administration alone may not produce effects that are representative of the effects of cannabis exposure in humans. Accordingly, THC self-administration has been difficult to demonstrate at the preclinical level (Tanda, 2016). Rats acquire stable rates of intravenous THC–CBD (10:1) self-administration only after extended passive pre-exposure to vaporized THC–CBD (Spencer et al., 2018). Thus, studying the effects of THC in the presence of other phytocannabinoids is important for modeling cannabis use.
Current animal models typically use the intravenous route for drug delivery, even though the most common route of cannabis use is inhalation (Sexton et al., 2016), and the pharmacokinetics of cannabinoids vary considerably depending on administration route (Grotenhermen, 2003; Huestis, 2007; Hložek et al., 2017). In humans, intravenous THC administration produces adverse effects often because of high dosing and fast infusion rates (Carbuto et al., 2012). Similarly, vaporized THC administration produces conditioned place preference, whereas intraperitoneal THC administration produces conditioned place avoidance in rodents (Manwell et al., 2014). Thus, both dose and route of administration may fundamentally influence the degree to which cannabis can support self-administration.
A more translationally relevant approach is needed that uses cannabis and mimics the most common route of administration in human users. With this in mind, we have developed a novel, ecologically valid model of cannabis vapor self-administration that uses “e-cigarette” technology to deliver discrete puffs of vaporized cannabis extracts to rodents in a response-contingent manner. We used this approach to examine whether vaporized cannabis extracts have reinforcing properties that support stable drug-taking behavior. We used whole-plant cannabis extracts predominantly containing THC or CBD, but also containing trace amounts of other phytocannabinoids often found in natural cannabis products. We characterized the metabolic phenotype of cannabis-trained rats and because human cannabis users and rodents treated with THC exhibit reduced CB1R availability (Ceccarini et al., 2015) and decreased CB1R expression/binding, particularly within the hippocampus (Romero et al., 1997; Silva et al., 2015; Farquhar et al., 2019; Kruse et al., 2019), we also measured hippocampal CB1 receptor binding 24 h after vapor self-administration. Given that acute abstinence can unmask withdrawal-related affective symptoms for other drugs (Hasin et al., 2008), we examined whether acute forced abstinence from cannabis vapor increases anxiety-like behavior. Finally, we tested whether vaporized cannabis maintains drug seeking in the drug-predictive context (i.e., under extinction conditions) or upon response-contingent presentation of a cannabis-paired light stimulus (i.e., cue-induced reinstatement).
Materials and Methods
Animals.
Male Sprague-Dawley rats (Simonsen Laboratories; 350–400 g) were pair-housed in a humidity-controlled vivarium on a 12 h reverse light/dark schedule (lights off at 07:00). Food and water were available ad libitum. All procedures followed the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Washington State University Institutional Animal Care and Use Committee. Table 1 describes the assignment of animals to various experiments, sample sizes, and relevant experimental parameters.
Drugs.
Cannabis extracts containing high concentrations of THC (CANTHC) or CBD (CANCBD) were obtained from the National Institute on Drug Abuse (NIDA) Drug Supply Program. According to the certificate of analyses provided upon shipment, the CANTHC extract contained 28.4% THC, 1.38% CBD, and 1.8% cannabinol (CBN). The CANCBD extract contained 1.16% THC, 59.34% CBD, 2.1% cannabichromene (CBC), 1.1% cannabigerol (CBG), and <0.01% tetrahydrocannabivarin and CBN. Extracts were heated to 60°C under constant stirring and dissolved in 80% propylene glycol/20% vegetable glycerol vehicle (VEH) at a concentration of 200 or 400 mg/ml based on previous studies (Nguyen et al., 2016; Javadi-Paydar et al., 2018). The 400 mg/ml concentration was used for all studies, except for studies involving plasma cannabinoid quantification where both 200 and 400 mg/ml concentrations were used. The final estimated concentrations of phytocannabinoids in 400 mg/ml CANTHC were as follows (in mg/ml): 113.6 THC, 5.5 CBD, and 7.2 CBN. The final concentrations in 200 mg/ml CANTHC were as follows (in mg/ml): 56.8 THC, 2.75 CBD, and 3.6 CBN. The final estimated concentrations in 400 mg/ml CANCBD were as follows (in mg/ml): 3.9 THC, 237.4 CBD, 8.4 CBC, and 4.4 CBG. The final concentrations of THC and CBD in 200 mg/ml CANCBD were as follows (in mg/ml): 1.85 THC, 118.7 CBD, 4.2 CBC, and 2.2 CBG. The CB1R antagonist, AM251 (Cayman Chemical), was dissolved in dimethyl sulfoxide:Tween-80:saline (1:1:18) vehicle and administered at a dose of 0, 1, or 3 mg/kg (1 ml/kg, i.p.).
Vapor chamber apparatus.
Eight 13.5 × 9.0 × 8.25 inch (L × W × H) 16.4 L vapor self-administration chambers (La Jolla Alcohol Research) were programmed using Med Associates IV software to deliver response-contingent puffs of vapor. An uninterrupted unidirectional flow of air entered through a port in the front of the chamber, and air was removed by vacuum in the rear of the chamber lid. The air intake port pulled air through tubing connected to an air flow meter and tubing connected to a commercial e-cigarette cartridge (first generation: Protank 3 Dual Coil; 2.2 Ω atomizer; Kanger Tech; second generation: SMOK Tank Baby Beast TFV8 with 0.2Ω M2 atomizer, 40–60 W range) filled with CANTHC, CANCBD, or VEH. Experiments involving radio telemetry, extinction, cue-induced reinstatement, CB1R antagonism, and quantification of plasma cannabinoid concentrations were conducted using first-generation vaporizers that delivered 10 s puffs of vapor. All other experiments were conducted using second-generation vaporizers that delivered 3 s puffs. These puff durations were chosen to approximate the volume of vapor produced by each vaporizer based on (1) the change in tank mass observed before versus after puff delivery, and (2) visual determination of full vapor clearance in 60 s. Two nose-poke operanda and associated cue lights were located on the rear wall of the chamber. The cartridge was connected to a vaporizer box (La Jolla Alcohol Research), and vapor puffs were delivered through the air intake port. Chamber air was evacuated through an activated charcoal filter (Carbatrol; first generation) or in-line Whatman HEPA-Cap filters (Millipore-Sigma; second generation).
Self-administration training.
Rats were trained to nose-poke for CANTHC, CANCBD, or VEH vapor puffs under a fixed ratio (FR)-1 reinforcement schedule during daily 1 h sessions on 11 consecutive days. Rats then progressed to an FR-2 schedule (Days 12–16) and then to an FR-4 schedule (Days 17–21). Nose-poke responses made on one (active) operandum resulted in a 3 s activation of the vaporizer and illumination of a cue light. The cue light remained illuminated during a 60 s timeout period, during which responses were not reinforced. Nose-poke responses made on the other (inactive) operandum had no programmed consequences. On the final self-administration day (Day 22), responding was reinforced under a progressive ratio (PR) schedule for 180 min. Schedule demand increased after each vapor delivery according to the following schedule: 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 7, 7, 9, 9, 11, 11, 13, 13, 15, 15, 18, 18, 21, 21, 24, 24, etc. (Walker and Koob, 2007). The breakpoint was defined as the total number of vapor deliveries obtained until responding ceased for a minimum of 15 min (Movie 1).
Cannabinoid quantification.
A different cohort of rats (N = 24) was trained to self-administer CANTHC or CANCBD (200 or 400 mg/ml) for 12–16 consecutive days on an FR-1 reinforcement schedule, and immediately following a single 1 h self-administration session, blood samples (∼100 μl) were collected to evaluate the relationship between the total number of vapor deliveries earned and circulating THC and CBD concentrations. Blood was collected vial the tail vein in sterile tubes containing 10 ml of ethylenediaminetetraacetic acid, centrifuged at 4°C at 4000 × g for 15 min, and stored at −20°C. THC and CBD concentrations in plasma were quantified immediately after the session as described previously (Britch et al., 2017; Greene et al., 2018).
A subset of rats from the FR/PR experiments (N = 21) were killed 24 h after their final self-administration session. Quantification of THC and CBD in brain tissue was performed as previously described (Baglot et al., under review). Briefly, frozen brain tissue was weighed and then manually homogenized (with a glass rod) in borosilicate glass culture tubes containing 2 ml of acetonitrile with 1 ng of THC-d3 and CBD-d3. Samples were then sonicated for 30 min in an ice bath and incubated overnight at −20°C to precipitate proteins. The following day samples were centrifuged at 1500 × g to remove particulates. The supernatant from each sample was transferred to a new glass tube and evaporated under nitrogen, the tube was then washed once with 350 μl of acetonitrile (to recapture any lipids adhering to the glass wall), and the acetonitrile was dried under nitrogen gas again. After completely drying, the samples were resuspended in 200 μl of 1:1 methanol–water and stored at −80°C until analysis by liquid chromatography mass spectrometry. Quantification of these molecules using mass spectrometry was performed using an Eksigent Micro LC200 coupled with an AB Sciex QTRAP 5500 mass spectrometry as previously described (Baglot SL, Petrie GN, Aukema RJ, Hume C, Zhuo R, Cole M, Parker LA, Rho J, Borgland SL, McLaughlin RJ, Brechenmacher L, Hill MN, unpublished observations). The data were acquired in positive electrospray ionization and multiple reaction monitoring mode and amount of each molecule was normalized to frozen tissue weight.
CB1R radioligand binding assay.
Frozen whole brain tissue was collected from some rats (N = 12) 24 h after the final PR challenge to assess hippocampal CB1R binding. Tissue was homogenized in TME buffer and centrifuged to generate the crude membrane fraction. Protein concentrations were determined using the Bradford method (Bio-Rad). Membranes (10 μg protein/sample) were incubated in TME buffer with [3H]CP55,940 (0.25, 0.5, 1.25, or 2.5 nm) in the absence or presence of AM251 (10 μm) to assess total and nonspecific binding, respectively. Bmax (maximal binding site density) and Kd (binding affinity) were calculated by nonlinear curve fitting to the single site Michaelis–Menten equation using GraphPad Prism as described previously (Lee and Hill, 2013; Berger et al., 2018).
Radio telemetry recordings.
Rats (N = 8) were anesthetized with an isoflurane/oxygen vapor mixture (isoflurane 5% induction, 1–3% maintenance) and a 1 cm midline vertical incision was made inferior to the xyphoid space. Sterile radio telemetry transmitters (Starr Life Sciences, PTD 4000) were inserted under the muscular layer and sutured in place using absorbable 4-0 silk sutures. The muscle layer was closed with 5-0 vicryl suture, and the skin was closed with non-absorbable wax-coated 4-0 silk suture. Rats received meloxicam (2 mg/kg, s.c.) for 3 d for postoperative pain management. After at least 5 d of recovery, rats were randomly assigned to treatment groups and trained to nose-poke for CANTHC, CANCBD, or VEH vapor on an FR-1 schedule of reinforcement. Radio telemetry transmissions indexing locomotor activity and body temperature were collected daily during 1 h self-administration sessions over the final 10 d of self-administration by Respironics ER-4000 receiver plates placed under the chambers.
Metabolic phenotyping.
A subset of rats (N = 17) was housed individually in metabolic cages (dimensions = L 18 × W 9.5 × H 8.1 inches; Promethion, Sable Systems) over the final 10 d of FR/PR self-administration training, during which feeding behavior, water intake, energy expenditure, physical activity and respiratory quotient (RQ) were monitored. Cages included a ceiling-mounted food hopper (3 mg resolution) and a water spigot connected to load cells (MM-1, Sable Systems) for food and water intake monitoring, respectively. Ad libitum access to the food hopper and water were allowed throughout the study. X- and y-axes (horizontal plane) photoelectric beam motion detectors were positioned around each cage to assess ambulatory activity. In addition to total distance traveled, bouts of inactivity lasting >15 s but <60 s were tabulated and averaged for each hour. Respiratory gases were measured with an integrated fuel cell oxygen analyzer, spectrophotometric CO2 analyzer, and capacitive water vapor partial pressure analyzer (GA3m1, Sable Systems). The Promethion system uses a pull-mode, negative pressure system. The multichannel mass flow generator measures and controls air flow (FR8-1, Sable Systems). The incurrent flow rate was set at 2000 ml/min. Water vapor was continuously measured and its dilution effect on O2 and CO2 was compensated for mathematically in the analysis (Lighton and Turner, 2008). RQ was calculated as the ratio of CO2 production to O2 consumption. Energy expenditure (i.e., the number of calories burned) was calculated using the Weir equation: Kcal/h = 60 × (0.003941 × VO2 + 0.001106 × VCO2; Weir, 1949). Data were acquired using MetaScreen v2.2.8, the raw data obtained were processed using ExpeData v1.8.2 (Sable Systems), and Macros 10 and 13 were used for data organization and transformation. Each day, rats were weighed, transferred to individual holding cages, and transported to the vapor self-administration system between 11:00 and 14:00. Daily values are averaged for active and inactive phases excluding the period encompassing transportation and self-administration training.
CB1R antagonism.
Rats (N = 16) were trained to self-administer vaporized CANTHC or CANCBD (400 mg/ml) under an FR-1 schedule over 26 daily sessions. Rats received mock intraperitoneal injections on the 2 d before testing the effects of systemic CB1R antagonism on training Days 16, 21, and 26. One hour before testing, rats received AM251 (1 or 3 mg/kg, i.p.) or VEH using a counterbalanced within-subjects design. Additional training sessions were conducted between test sessions to re-establish stable responding. All data were converted to a percentage change score relative to the previous mock injection day.
Elevated plus maze test.
The elevated plus maze (EPM) apparatus consisted of a raised Plexiglas platform (28.5 inches high) with two open exposed arms and two darker enclosed arms of equal length (21.5 inches/arm; Med Associates). The floors were made of clear Plexiglas and the walls of both closed arms were black. Twenty-four hours after their final FR-1 self-administration session, rats (N = 37) were individually placed in the center of the maze and allowed to freely explore the maze for 5 min. All tests were run in dim lighting (∼10 lux), and behaviors were recorded with Noldus EthoVision XT behavioral tracking software. The number of entries and percentage time spent in the open and closed arms, and distance traveled in the open arms of the EPM were compared across groups. The frequency of risk assessment behaviors (i.e., head-dips, stretch-attend postures, and rearing) was also scored manually by trained research assistants blinded to treatment conditions as described previously (Henricks et al., 2017; Berger et al., 2018).
Extinction training and cue-induced reinstatement.
Rats (N = 35) were trained to self-administer CANTHC, CANCBD, or VEH under an FR-1 schedule over 19 daily sessions. Rats then underwent daily 1 h extinction training sessions during which nose-poke responses had no programmed consequences. Extinction training continued until rats reached the extinction criterion (i.e., ≥7 extinction training sessions with ≥50% decrease in active nose-poke responses during the final 2 sessions). Rats were then tested for cue-induced reinstatement of extinguished drug-seeking behavior. During the 1 h test session, active nose-poke responses resulted in cue light presentations without vapor delivery.
Experimental design and statistical analyses.
Studies described herein were conducted as separate experiments using different cohorts of rats. Data pertaining to schedules of reinforcement, CB1R binding, tissue cannabinoid concentration, and metabolic phenotyping were collected in the same cohort using second-generation vaporizers. Data pertaining to extinction, cue-induced reinstatement, radio telemetry, and anxiety-like behavior were collected in the same cohort using first generation vaporizers. Data pertaining to plasma cannabinoid concentration and CB1R antagonism were each collected in separate cohorts using first generation vaporizers. Vapor deliveries, active and inactive responses, trials to extinction, nose-poke discrimination, body weight, within-session locomotor activity, food intake, water intake, lounge time, locomotor activity, energy expenditure, RQ, and EPM data were compared across groups using one-way or mixed factorial ANOVA with treatment group (CANTHC, CANCBD, VEH) as the between-subjects factor and time (day, reinforcement schedule, 15 min bin, light/dark phase) as the within-subjects factor. Nose-poke discrimination was calculated based on a formula used by Spencer et al. (2018): nose-poke discrimination index = (active nose-pokes − inactive nose-pokes)/(active nose-pokes + inactive nose-pokes), where 0 indicated no discrimination between active and inactive operanda, 1 indicated complete preference for the active operandum, and −1 indicated complete preference for inactive operandum. A discrimination index of ≥0.33 is equivalent to 2:1 active–inactive nose-pokes. Effects of AM251 were analyzed using separate repeated-measures one-way ANOVAs. Significant effects were further probed using Tukey or Bonferroni post hoc tests. Vapor delivery and plasma cannabinoid concentrations were correlated using Pearson's correlation coefficient (r). Alpha was set at 0.05. Effect sizes are reported as ηP2.
Results
THC-rich cannabis vapor supports stable rates of self-administration
A mixed-factorial ANOVA conducted with reinforcement schedule as the within-subjects factor (averaged over all FR days for each schedule) revealed main effects of treatment (F(2,27) = 22.20, p < 0.0001, ηP2 = 0.62) and schedule (F(2,54) = 48.37, p < 0.0001, ηP2 = 0.62) on active responding, as well as a significant treatment × schedule interaction (F(4,54) = 9.36, p < 0.0001, ηP2 = 0.41; Fig. 1C, colored symbols). Tukey's post hoc tests indicated that CANTHC elicited more active responses than CANCBD on FR-1 (p = 0.02), FR-2 (p < 0.0001), and FR-4 (p < 0.0001) schedules, and more active responses than VEH on FR-2 (p = 0.0008) and FR-4 (p < 0.0001) schedules. A mixed-factorial ANOVA also revealed a main effect of treatment on inactive responding (F(2,27) = 6.49, p = 0.005, ηP2 = 0.32), but no effect of schedule (F(2,54) = 2.79, p = 0.07, ηP2 = 0.09) or treatment × schedule interaction (F(4,54) = 1.39, p = 0.25, ηP2 = 0.09; Fig. 1C, open symbols). Tukey's post hoc tests indicated that CANCBD elicited fewer inactive responses than CANTHC (p = 0.02) and VEH (p = 0.01), regardless of reinforcement schedule.
The majority of active responding occurred during the first 15 min of the session under each reinforcement schedule (main effect of bin: F(3,81) = 48.1–86.6, p values < 0.0001, ηP2 = 0.64–0.76). Significant treatment × bin interactions were observed for each reinforcement schedule (F(6,81) = 5.82–9.05, p values < 0.0001, ηP2 = 0.30–0.40), and Tukey's post hoc tests indicated that CANTHC elicited more active responses than CANCBD during the first, second, and fourth bins (all p values ≤ 0.02) and elicited more active responses than VEH and CANCBD during each bin under FR-2 and FR-4 schedules (all p values ≤ 0.03). Conversely, CANCBD elicited fewer active responses than VEH during the first 15 min bin under each reinforcement schedule (all p values ≤ 0.04). There were also main effects of treatment (F(2,27) = 4.84–5.58, p values ≤ 0.02, ηP2 = 0.26–0.29) and bin (F(3,81) = 17.6–26.4, p values < 0.0001, ηP2 = 0.39–0.49) on the number of inactive responses made, as well as significant treatment × bin interactions for each reinforcement schedule (F(6,81) = 6.26–9.47, p values ≤ 0.0001, ηP2 = 0.32–0.41). Tukey's post hoc tests revealed significant differences in inactive responding during the first 15 min bin, with VEH > CANTHC > CANCBD (all p values ≤ 0.02). During the second 15 min bin under the FR-4 schedule, CANTHC elicited significantly more inactive responses than CANCBD (p = 0.01) and VEH (p = 0.04), but no differences were observed between CANCBD and VEH for this bin or any other bin.
Stable rates of vapor delivery were achieved under each reinforcement schedule, as indicated by a lack of effect of day over the final 3 d under each schedule (Fig. 1D). A mixed-factorial ANOVA revealed main effects of treatment (F(2,27) = 18.99, p < 0.0001, ηP2 = 0.58) and schedule (F(2,54) = 19.01, p < 0.0001, ηP2 = 0.41) on the number of vapor deliveries earned, as well as a significant treatment × schedule interaction (F(4,54) = 4.83, p = 0.002, ηP2 = 0.26; Fig. 1D, colored symbols). CANTHC elicited more vapor deliveries than CANCBD on FR-1, FR-2, and FR-4 schedules, and more vapor deliveries than VEH on FR-2 and FR-4 schedules (Tukey: all p values < 0.0001). There were also main effects of treatment (F(2,27) = 8.72–22.0, p values < 0.001, ηP2 = 0.39–0.62) and bin (F(3,81) = 53.6–224.7, p values < 0.0001, ηP2 = 0.67–0.89), and significant treatment × bin interactions (F(6,81) = 6.79–12.3, p values < 0.0001, ηP2 = 0.67–0.89) at each reinforcement schedule (Fig. 1E–G). Tukey's post hoc tests indicated that CANTHC elicited more vapor deliveries than CANCBD during each bin under all reinforcement schedules (all p values ≤ 0.03), and elicited more vapor deliveries than VEH during each bin under FR-2 and FR-4 schedules (all p values ≤ 0.001). CANCBD elicited fewer vapor deliveries than VEH during for first 15 min bin under the FR-1 schedule (p < 0.0001).
One-sample t tests revealed that the average nose-poke discrimination index was significantly >0.33 for CANTHC (t(20) = 4.39, p < 0.001) and CANCBD (t(20) = 8.77, p < 0.001), but not for VEH. Thus, only the cannabis vapor self-administering groups obtained active to inactive responding ratios that were >2:1. A mixed-factorial ANOVA revealed main effects of treatment (F(2,27) = 8.33, p = 0.002, ηP2 = 0.38) and reinforcement schedule (F(2,54) = 20.9, p < 0.0001, ηP2 = 0.44) on response discrimination, but no treatment × schedule interaction (F(4,54) = 0.15, p = 0.96; Fig. 1H). Tukey's post hoc tests indicated that discrimination was significantly better under FR-2 and FR-4 schedules relative to FR-1 (p values < 0.05) and that CANCBD elicited better discrimination for the active operandum than VEH (p = 0.001; Fig. 1H).
THC-rich cannabis vapor exhibits motivational properties
Cumulative responding under the PR schedule was analyzed using a mixed-factorial ANOVA that revealed main effects of treatment (F(2,27) = 4.94, p = 0.02, ηP2 = 0.27) and time (F(11,297) = 32.0, p < 0.0001, ηP2 = 0.54), as well as a significant treatment × time interaction (F(22,297) = 4.58, p < 0.0001, ηP2 = 0.25; Fig. 2A). Bonferroni post hoc tests correcting for multiple comparisons indicated that CANTHC produced higher cumulative responding compared with VEH from 105 to 180 min (all p values ≤ .03) and compared with CANCBD from 135 to 180 min (all p values ≤ 0.02; Fig. 2A). A one-way ANOVA indicated that CANTHC also elicited higher break points than VEH (F(2,27) = 7.17, p = 0.001, ηP2 = 0.34; Tukey post hoc, p = 0.001) and a trend for higher break points than CANCBD (Tukey post hoc, p = 0.058; Fig. 2B). Moreover, a one-way ANOVA indicated that CANTHC elicited shorter latencies to re-initiate responding following vapor deliveries compared with VEH and CANCBD (F(2,27) = 5.01, p = 0.01, ηP2 = 0.27; Tukey, p values ≤ 0.04; Fig. 2C). Importantly, responding for CANCBD was not significantly different from VEH on any endpoint measured during the PR challenge.
Cannabis vapor self-administration produces biologically relevant increases in plasma THC and CBD concentrations
The number of CANTHC vapor deliveries (200 or 400 mg/ml) measured during the 1 h self-administration session in a separate cohort of rats positively correlated with the concentration of THC measured in plasma immediately post-session (CANTHC-200: r = 0.51, p = 0.03; CANTHC-400: r = 0.86, p < 0.001; Fig. 3A). A similar relationship was observed with the number of CANCBD vapor deliveries and the plasma concentration of CBD immediately post-session (CANCBD-200: r = 0.58, p = 0.01; CANCBD-400: r = 0.51, p = 0.18; Fig. 3B). Brain tissue concentrations of THC and CBD were measured 24 h after the PR test session from Experiment 1. Independent samples t tests revealed that brain THC concentrations were not significantly different between CANTHC and CANCBD groups (t(20) = 0.77, p = 0.45, η2 = 0.03; Fig. 3C). However, brain tissue concentrations of CBD were higher in the CANCBD group compared with the CANTHC group (t(20) = 2.24, p = 0.04, η2 = 0.20; Fig. 3D).
CANTHC self-administration alters hippocampal CB1 receptor binding
CB1R radioligand binding assays were conducted on hippocampal tissue taken from the same brains from Experiment 1 (i.e., 24 h after the final PR session). Planned t test comparisons between CANTHC and VEH groups revealed that CANTHC self-administration significantly reduced CB1R maximal binding site density (Bmax; t(6) = 2.90, p = 0.03, η2 = 0.58; Fig. 3E) without altering CB1R binding affinity (KD; t(6) = 0.53, p = 0.61, η2 = 0.05; Fig. 3F) compared with VEH, whereas CANCBD did not significantly alter CB1R binding site density (t(6) = 1.31, p = 0.24, η2 = 0.22; Fig. 3E) or binding affinity (t(6) = 1.12, p = 0.31, η2 = 0.17; Fig. 3F) when measured 24 h after the final vapor self-administration session.
CANTHC vapor self-administration alters physical activity and daily food intake
Radio telemetry recordings were taken over the final 10 d of FR-1 self-administration (D10–19) using a separate cohort of rats. A mixed-factorial ANOVA revealed a main effect of treatment (F(2,5) = 5.97, p = 0.047, ηP2 = 0.70), with CANTHC significantly reducing within-session locomotor activity relative to VEH (Tukey, p < 0.05; Fig. 4A). There was no effect of day (F(9,45) = 1.64, p = 0.13, ηP2 = 0.25) or treatment × day interaction (F(18,45) = 1.09, p = 0.39, ηP2 = 0.30). The number of vapor deliveries negatively correlated with locomotor activity on Days 13, 16, and 18 (r values = 0.73–98; p values < 0.05).
Metabolic phenotyping of a subset of rats from Experiment 1 was conducted over the final 10 d of self-administration. When time spent inactive was measured daily during the first 3 h post-vapor exposure, a mixed-factorial ANOVA revealed main effects of treatment (F(2,14) = 10.2, p = 0.002, ηP2 = 0.59) and time (F(2,28) = 51.2, p < 0.0001, ηP2 = 0.79), but no significant interaction (F(4,28) = 1.37, p = 0.27, ηP2 = 0.16). Tukey's post hoc tests indicated that CANTHC increased inactivity relative to CANCBD and VEH (p values < 0.05; Fig. 4B). When cumulative daily inactivity time was calculated and averaged over the final 10 d, a one-way ANOVA revealed a main effect of treatment (F(2,14) = 8.0, p = 0.005, ηP2 = 0.53), with CANTHC self-administration resulting in more daily inactivity time compared with VEH (Tukey, p < 0.05; Fig. 4C). As expected, there was a main effect of light/dark phase on the average distance traveled per day (F(1,14) = 166, p < 0.0001, ηP2 = 0.92), but no effect of treatment (F(2,14) = 2.79, p = 0.10, ηP2 = 0.28) or treatment × phase interaction (F(2,14) = 1.82, p = 0.20, ηP2 = 0.21; Fig. 4D). A one-way ANOVA indicated a main effect of treatment on total daily food intake (F(2,14) = 5.73, p = 0.02, ηP2 = 0.42) with CANTHC increasing mean food intake over the last 10 d of training relative to VEH and CANCBD (Tukey, both p < 0.03; Fig. 4E). A one-way ANOVA also revealed a main effect of day on weight gain (F(13,169) = 80.8, p < 0.0001, ηP2 = 0.86), but there was no effect of treatment (F(2,13) = 1.11, p = 0.36, ηP2 = 0.15) or treatment × day interaction (F(26,169) = 0.88, p = 0.64, ηP2 = 0.12), which indicates that neither CANTHC nor CANCBD significantly altered bodyweight gain. Pearson correlation analyses examining associations between hourly/daily food intake and average vapor deliveries earned over the last 10 d did not reveal any significant relationships between these variables. O2 consumption (VO2) and CO2 consumption (VCO2) differed according to phase (F(1,14) = 98.8–140.00, p values < 0.001, ηP2 = 0.88–0.91) and treatment (F(1,14) = 4.04–4.49, p values = 0.03, ηP2 = 0.22–0.24), with CANTHC producing higher VO2 and VCO2 values compared with CANCBD when averaged over the final 10 d (Tukey: p values < 0.05; Fig. 4F,G). A mixed-factorial ANOVA assessing energy expenditure also revealed main effects of phase (F(1,14) = 179.00, p < 0.001, ηP2 = 0.93) and treatment (F(2,14) = 5.17, p = 0.021, ηP2 = 0.42), but no phase × treatment interaction (F(2,14) = 0.28, p = 0.76; Fig. 4H). Tukey's post hoc tests revealed that CANTHC increased energy expenditure compared with both VEH and CANCBD (both p < 0.05; Fig. 4H). Total food and water intake were higher during the active versus inactive phase (F(1,14) = 25.0–166, p values < 0.001, ηP2 = 0.64–0.92) but did not differ according to treatment, and there were no effects of phase or treatment on RQs.
The reinforcing effects of vaporized CANTHC require CB1 receptor stimulation
Systemic CB1R antagonism differentially impaired the reinforcing effects of vaporized CANTHC and CANCBD. Baseline active responding and vapor deliveries were not statistically different between conditions, and VEH treatment did not alter these measures relative to mock injection. Repeated-measures ANOVAs revealed a main effect of AM251 treatment on active responding (F(2,14) = 6.05, p = 0.01, ηP2 = 0.46) and vapor deliveries earned (F(2,14) = 7.44, p = 0.01, ηP2 = 0.52), with the 1 and 3 mg/kg dose of AM251 decreasing CANTHC-reinforced active responses (Tukey, p values = 0.04 and 0.02, respectively) and vapor deliveries (Tukey, p values = 0.03 and 0.01, respectively) relative to baseline (Fig. 5A,B). In contrast, there was no effect of AM251 treatment on CANCBD-reinforced active responding (F(2,14) = 1.17, p = 0.34, ηP2 = 0.14; Fig. 5C) or the number of CANCBD vapor deliveries earned (F(2,14) = 2.11, p = 0.16, ηP2 = 0.23; Fig. 5D). AM251 treatment did not significantly alter inactive responding for CANTHC (F(2,14) = 0.05, p = 0.95) or CANCBD (F(2,14) = 1.98, p = 0.18).
Acute withdrawal from cannabis vapor does not elicit anxiety-like behavior
Behaviors in the EPM were measured 24 h after the final vapor self-administration session, and one-way ANOVAs indicated that neither CANTHC nor CANCBD history altered the percentage of time spent in the open arms of the EPM (F(2,34) = 0.53, p = 0.48, ηP2 = 0.03), the number of open arm entries made (F(2,34) = 1.12, p = 0.34, ηP2 = 0.06), the number of rearing events (F(2,34) = 0.12, p = 0.89, ηP2 = 0.007), or the frequency of stretch-attend postures (F(2,34) = 2.90, p = 0.07, ηP2 = 0.15; for descriptive statistics, see Table 2).
Removal of cannabis vapor reinforcement elicits an extinction burst
Rates of vapor self-administration did not differ between groups on the final self-administration day (F(2,32) = 3.1, p = 0.06, ηP2 = 0.16; Fig. 6A). A mixed factorial ANOVA examining effects of treatment on active responding on the last day of self-administration relative to the first day of extinction revealed a main effect of day (F(1,32) = 27.1, p < 0.0001,ηP2 = 0.46) and a significant treatment × day interaction (F(2,32) = 3.47, p = 0.04, ηP2 = 0.17). Bonferroni post hoc tests indicated that removal of CANTHC or CANCBD increased active responding on the first extinction day relative to the last self-administration day (both p < 0.005), whereas removal of VEH did not alter active responses (Fig. 6A). A mixed factorial ANOVA conducted over the first seven extinction days revealed a main effect of day (F(6,192) = 7.40, p < 0.0001, ηP2 = 0.33) and a significant treatment × day interaction (F(12,192) = 1.95, p = 0.03, ηP2 = 0.11). Bonferroni post hoc tests indicated that removal of CANTHC elicited more active responses than VEH on extinction Days 1–3 (all p values ≤ 0.03). There was also a significant treatment × day interaction for inactive responding (F(12,192) = 3.13, p = 0.0004, ηP2 = 0.16), such that removal of CANCBD elicited fewer inactive responses than removal of VEH on Day 6 (Bonferroni: p = 0.044) and fewer than removal of CANTHC on Day 4 (Bonferroni, p = 0.004; Fig. 6A). Interestingly, removal of CANCBD increased the number of sessions rats needed to reach the extinction criterion compared with CANTHC (F(2,32) = 4.09, p = 0.03, ηP2 = 0.20, Tukey post hoc, p = 0.03; Fig. 6B).
CANTHC- or CANCBD-paired stimuli elicit reinstatement of cannabis vapor-seeking behavior
A mixed factorial ANOVA examining vapor cue-induced active responding revealed that response-contingent presentation of either the CANTHC- or the CANCBD-paired (but not the VEH-paired) light cue increased active responses at test relative to the last extinction training day (interaction: F(2,32) = 3.14, p = 0.05, ηP2 = 0.18; Bonferroni post hoc, p values < 0.05; Fig. 6C). Response-contingent cue presentation also altered inactive responding over the final extinction and reinstatement test days (F(2,32) = 3.34, p = 0.05, ηP2 = 0.17), with CANCBD rats making fewer inactive responses than VEH rats regardless of cue presentation (p = 0.04). Planned comparison t tests examining effects of CANTHC and CANCBD on cue-induced active responding relative to VEH indicated that the CANTHC-paired cue elicited more active responses than the VEH-paired cue (t(22) = 2.11, p = 0.05, ηP2 = 0.16; Fig. 6C). There were no differences between CANCBD and VEH for cue-induced active responding (t(20) = 1.14, p = 0.27), and no significant differences in inactive responding between CANTHC and VEH (t(22) = 0.13, p = 0.90) or CANCBD and VEH (t(20) = 1.42, p = 0.17) groups.
Discussion
Obstacles in establishing a model of cannabis use that closely mimics the human experience have limited our ability to study the neural mechanisms and impact of cannabis use (Melis et al., 2017; McLaughlin, 2018). In the current study, we present evidence supporting the feasibility of a novel preclinical model of cannabis self-administration that employs response-contingent delivery of vaporized cannabis extracts. Our findings indicate that THC-dominant cannabis vapor (CANTHC) has reinforcing properties. Despite variability in rates of self-administration, both cannabis self-administering groups exhibited stable rates of responding and demonstrated robust discrimination between active and inactive operanda. Although all rats show increased responding early in the session, responding for CANTHC was significantly higher than VEH and CANCBD during this time, which is consistent with a preserved loading dose phenomenon. Importantly, active responding was highest for CANTHC, and only the CANTHC self-administering group defended the number of vapor deliveries earned under more demanding schedules of reinforcement. Finally, CANTHC elicited higher break points under the PR reinforcement schedule and supported the acquisition of conditioned motivational effects by an initially neutral environmental stimulus. These findings support previous efforts demonstrating the feasibility of vapor self-administration for other drugs of abuse, such as opioids (Jaffe et al., 1989; Weinhold et al., 1993; Vendruscolo et al., 2018). To our knowledge however, this is the first evidence supporting vapor self-administration as a viable means to investigate cannabis-seeking behavior.
The validity of this model is further strengthened by our findings that response-contingent delivery of vaporized cannabis extracts produced behavioral and physiological changes that are consistent with studies in rodents using passive THC delivery, as well as observations from human cannabis users. Despite differences in peak intoxication times due to inherent variability in the patterns of responding, our data indicate that the number of CANTHC and CANCBD vapor deliveries obtained positively correlated with THC and CBD plasma concentrations immediately after the session. Self-administered CANTHC reliably reduced locomotor activity during the session and increased the duration of inactivity in the metabolic chambers during the first 2 h following the self-administration session, which is congruent with recent studies using passive THC vapor delivery (Nguyen et al., 2016). CANTHC also increased daily food intake and energy expenditure. These data mirror observations of increased caloric intake and energy expenditure of chronic cannabis users (Smit and Crespo, 2001; Rodondi et al., 2006; Le Strat and Le Foll, 2011; Meier et al., 2019). Moreover, CB1R binding was significantly reduced in the hippocampus 24 h after the final self-administration in CANTHC rats compared with VEH rats, which is in line with studies in human cannabis users (Ceccarini et al., 2015) and rodents repeatedly injected with THC (Romero et al., 1997; Silva et al., 2015; Farquhar et al., 2019; Kruse et al., 2019). Notably, biologically relevant THC and CBD concentrations were observed in brain tissue 24 h after the final self-administration session, which may account for the lack of abstinence-induced anxiety-like behavior observed in this study. Brain CBD concentrations were higher than brain THC concentrations, probably because of the long half-life of CBD (i.e., 24–48 h; Millar et al., 2018). THC tissue concentration was low but detectable, and similar in both CANTHC and CANCBD rats. CBD inhibits the cytochrome P450 family of liver enzymes that are responsible for THC metabolism (Zendulka et al., 2016). Thus, comparable brain THC concentrations following CANTHC and CANCBD regimens might reflect greater CBD-mediated inhibition of THC metabolism in the CANCBD group. This possibility will need to be systematically evaluated in future studies.
The vapor self-administration procedure likely facilitated volitional cannabis exposure. It has been well documented for other drugs of abuse that control over drug administration profoundly alters the subjective effects of drug intake, as well as the associated neurochemical responses (Stefański et al., 1999, 2007; Donny et al., 2000). Following passive intravenous THC administration, humans report aversive effects that are most often caused by the dose and rate of infusion (Carbuto et al., 2012). These factors likely also influence the feasibility of intravenous self-administration in rodents. Notably, we used cannabis extracts containing THC, CBD, and other naturally-occurring phytocannabinoids, which could also mitigate the aversive effects of high-dose THC exposure (Russo and Guy, 2006). Notably, CBD facilitates intravenous THC self-administration in rodents (Spencer et al., 2018; but see Wakeford et al., 2017) and offsets some of the pharmacological and behavioral effects of THC in humans (Morgan et al., 2010) and rodents (Englund et al., 2013). Furthermore, a passive THC+CBD vapor pre-exposure regimen facilitates the acquisition of intravenous THC+CBD self-administration (Spencer et al., 2018). However, it is unknown whether CBD is required for self-administration of vaporized THC. It should be noted that difficulties in establishing intravenous THC self-administration may be due more to differences in dosing, procedures, and species, rather than the presence of particular phytocannabinoids per se (for review, see Panlilio et al., 2010). Thus, future studies should test both whole cannabis extracts, as well as isolated constituents in combination to evaluate interactions between THC, CBD, and other phytocannabinoids.
THC and CBD concentrations vary widely across commercially available cannabis products. As such, it will be important to study the extent to which modifying cannabinoid constituent concentrations alters the propensity for self-administration. In the current study, we compared responding for two different cannabis extracts and observed key differences between their effects. First, CANTHC, but not CANCBD, elicited significant reinforcing effects as indicated by augmented operant responding upon increases in schedule demand. Second, CANTHC had greater motivational effects as indicated by higher breakpoints and a reduced latency to initiate responding following vapor delivery compared with CANCBD. Although both CANTHC and CANCBD self-administration regimens were sufficient to increase drug seeking upon the removal of the reinforcer (i.e., extinction) or in the presence of vapor-associated cues (i.e., reinstatement) relative to responding during the last self-administration session and relative to the absence of cues, respectively, only the CANTHC regimen elicited response rates above that of the VEH group. Thus, CANTHC vapor has greater reinforcing properties and stronger conditioned motivational effects than CANCBD and VEH vapor.
Interestingly, the CANCBD regimen produced the strongest discrimination between active and inactive operanda despite reduced rates of responding. Furthermore, CANCBD elicited drug-seeking behavior that was more resistant to extinction, as indicated by a larger number of sessions required to reach extinction criterion relative to CANTHC. Detectable brain concentrations of THC and a CBD-mediated decrease in the aversive properties of THC likely reconcile these findings with studies indicating a lack of CBD-mediated reinforcement (Haney et al., 2016; Viudez-Martínez et al., 2019). Additionally, CBD enhances the efficacy of THC to produce discriminative stimulus effects in rhesus monkeys (McMahon, 2016). Thus, the high concentration of CBD in the CANCBD extract may have augmented the discriminative effects of THC, resulting in more persistent, inflexible drug seeking. If this is the case, then high-CBD cannabis products that also contain THC may still have abuse liability concerns akin to more THC-dominant products. More rigorous interrogation into effects of vaporized CBD will undoubtedly provide insight into potential reinforcing properties of CBD in the presence and absence of THC (Javadi-Paydar et al., 2018, 2019).
The reinforcing properties of cannabis vapor self-administration require CB1R stimulation, as our data indicate that rates of CANTHC- (but not CANCBD-) reinforced responding were decreased following acute AM251 administration. This is consistent with studies in nonhuman primates indicating that CB1R blockade selectively attenuates THC reinforcement without altering cocaine- or food-reinforced responding (Tanda et al., 2000; Schindler et al., 2016). AM251 administration may have also induced a state of precipitated withdrawal that could have influenced responding. Although we did not find evidence for spontaneous withdrawal in the EPM, precipitated withdrawal has been more reliably observed in THC-treated rodents (Cook et al., 1998; Lichtman et al., 2001) and nonhuman primates (Stewart and McMahon, 2010). Thus, future studies should examine whether CB1R antagonism causes precipitated withdrawal in this model.
Although the use of whole-plant cannabis extracts more accurately models human consumption, this approach does introduce additional sources of variability. For instance, although concentrations of some primary phytocannabinoids were provided in the NIDA certificate of analysis, we were only provided with details on a subset of phytocannabinoids and did not perform additional onsite analytics. We also used two different vaporizers/vapor delivery regimens and multiple batches of extracts for these experiments. Thus, the vapor quantity administered per puff and the concentration of phytocannabinoids present in these extracts may have varied slightly between experiments. However, despite these added sources of variability, the data from these experiments are remarkably consistent in demonstrating the reinforcing effects of vaporized cannabis delivery. Nevertheless, more rigorous investigation into cannabis vapor dosing parameters and kinetics is needed to better inform future studies using this approach.
Altogether, findings from the present study strongly support the utility of a response-contingent vapor delivery protocol as a means to model cannabis use and further examine the neurobiological mechanisms underlying cannabis-seeking behavior. Ultimately, this model will permit finer interrogation of the effects of cannabis on the brain and behavior and help to identify causal factors that increase the susceptibility for developing cannabis use disorders.
Footnotes
This work was supported by NIH NIDA Grants R21 DA043722-01A1 (R.J.M.) and R01 DA025646-07 (R.A.F.), by a Foundation Grant from the Canadian Institutes of Health Research to M.N.H., and with funding for medical and biological research by the State of Washington Initiative Measure No. 171 (R.J.M.). We thank Maury Cole for technical support and the Southern Alberta Mass Spectrometry Centre, located in and supported by the Cumming School of Medicine, University of Calgary, for their services in targeted liquid chromatography tandem mass spectrometry.
The authors declare no competing financial interests.
- Correspondence should be addressed to Ryan J. McLaughlin at ryan.mclaughlin{at}wsu.edu