Δ 9-Tetrahydrocannabinol During Adolescence Reprograms the Nucleus Accumbens Transcriptome, Affecting Reward Processing, Impulsivity, and Specific Aspects of Cocaine Addiction-Like Behavior in a Sex-Dependent Manner
Department of Psychobiology, School of Psychology, National University for Distance Learning (UNED), Madrid, Spain
Institut de Neurosciences de la Timone, Marseille, France
Correspondence: Alejandro Higuera-Matas, PhD, Department of Psychobiology, School of Psychology, National University for Distance Learning (UNED), C/Juan del Rosal 10, Madrid, Spain (ahiguera@psi.uned.es).Abstract
Background
Cannabis exposure during adolescence is associated with emotional and motivational alterations that may entail an enhanced risk of developing psychiatric disorders. In rodent models, exposure to cannabinoids during adolescence leads to increased self-administration of opiates and cocaine, however, the psychological and neural mechanisms and the sex-specificity of this phenomenon are largely unknown.
Methods
We exposed male and female adolescent rats to Δ9-tetrahydrocannabinol (THC) and studied at adulthood the effects of such treatment on psychological processes related to reward, such as Pavlovian conditioned approach, Pavlovian to instrumental transfer, habit formation and waiting impulsivity. In the light of these data and given the involvement of the nucleus accumbens in the processes examined, we performed an RNASeq transcriptomic study and assessed cocaine addiction-like behavior.
Results
THC exposure increased goal-tracking (in males and females) and enhanced Pavlovian to instrumental transfer (especially in males) but did not affect habit formation. THC-exposed rats exhibited subtle, state-dependent changes in premature responding in the 2-CSRTT task. RNASeq data showed gene expression alterations in a marked sex-specific manner. While no effects were found on the acquisition of cocaine self-administration or punished drug-seeking, rats exposed to THC self-administered more cocaine under a progressive ratio schedule (males), had a higher rebound upon returning to continuous access to the drug (females) and showed reduced drug-seeking after 30 days of withdrawal (females).
Conclusions
Adolescent THC affects specific aspects of reward- (and cocaine-) guided behavior and the function of a key brain region mediating these effects, in a remarkable sex-specific manner.
Boxed Text
Cannabis is one of the psychoactive drugs most widely used by adolescents. Although a great deal of research has been carried out regarding its long-term effects, it is unknown if a chronic cannabinoid exposure during adolescence could affect the psychological processes governing reward-guided actions and cocaine addiction, and what could be the underlying neurobiological mechanisms. Here, we found that THC-exposed male rats showed more potentiation of reward seeking by stimuli predictive of those rewards but were less attracted to them. THC also decreased the tolerance to delays, especially in females. All this was accompanied by sex-specific alterations in the activation or deactivation of different families of genes in the nucleus accumbens, a key region of the reward circuit. We also observed a potentiation but also reduction of specific aspects of cocaine addiction, so these results do not fully support the Gateway Hypothesis of drug use.
Introduction
Adolescence is a crucial period of development characterized by profound changes in psychological and neural processes (Spear, 2000; Paus et al., 2008; Blakemore, 2012). As a result, any insult such as stressful events or drug use during this period will have several consequences at the psychological and neurobiological levels. Cannabis is the drug of abuse—other than alcohol and tobacco—most widely consumed by adolescents (EMCDA, 2019), and the exposure to this drug during adolescence has profound consequences for the developing individual (Higuera-Matas et al., 2015; Rubino and Parolaro, 2016; Ferland and Hurd, 2020; Hurd, 2020; Stringfield and Torregrossa, 2021), which may be more severe than when exposure occurs exclusively during adulthood. A potential consequence of adolescent cannabinoid exposure (ACE) is an increase in the use and/or abuse liability of other drugs later in life (termed Gateway Hypothesis, which is under intense debate (Kandel et al., 1992, 2003; Fergusson et al., 2006; Tarter et al., 2006; Vanyukov et al., 2012; Kleinig, 2015; Mayet et al., 2016; Nkansah-Amankra and Minelli, 2016; Lynskey and Agrawal, 2018)). Previous experiments by our group and others have suggested that animals with ACE show increased morphine (Biscaia et al., 2008a), heroin (Ellgren et al., 2007; Tomasiewicz et al., 2012; Lecca et al., 2020), fentanyl (but not oxycodone) (Nguyen et al., 2020), and cocaine self-administration (SA) (Higuera-Matas et al., 2008; Friedman et al., 2019). However, others have reported a delayed cocaine SA acquisition (Kononoff et al., 2018) or no changes in heroin SA (Stopponi et al., 2014). In addition, these studies were typically performed on male rats, neglecting the sex-dependent effects that are common when both sexes are included (Biscaia et al., 2008b; Higuera-Matas et al., 2008). Thus, a more detailed examination of the psychological and neurobiological processes involved in the increased SA of cocaine and opiates after ACE is warranted.
Several reward-related processes may be responsible for the increased use of the drug or may facilitate the development of addictive behavioral patterns. These processes include Pavlovian to instrumental transfer (PIT; i.e., the ability of classically conditioned cues to affect instrumental responses) (Cartoni et al., 2016), Pavlovian conditioned approach (a measure of incentive salience) (Fitzpatrick and Morrow, 2016), habit formation (Belin et al., 2013; Everitt and Robbins, 2013), and impulsivity, a core endophenotype that predicts the development of cocaine addiction (Belin et al., 2008; Jupp and Dalley, 2014). Noteworthy, the nucleus accumbens (NAcc) participates in many of these processes, and prior studies of ACE have highlighted alterations in this structure (Higuera-Matas et al., 2015; Stringfield and Torregrossa, 2021). Recent evidence suggests that exposure to a cannabinoid agonist modifies the initial responses of this region to cocaine (Scherma et al., 2020). However, this study did not examine the potential sex-specific alterations induced by ACE in the transcriptomic landscape of the NAcc, and they used WIN 55,212-2 instead of the actual phytocannabinoid, Δ 9-tetrahydrocannabinol (THC). Therefore, to gain a deeper and broader understanding of the neurochemical alterations induced by adolescent THC in the NAcc and given the crucial role of sex differences in the effects of cannabinoids (Viveros et al., 2011), we performed an RNASeq study in the NAcc of adult male and female rats that had been exposed to THC as adolescents.
In addition to exploring the potential alterations in reward-related behaviors, impulsivity, and the potential accumbal alterations involved, we decided to explore in more depth several characteristics of cocaine SA that could indicate an altered tendency to develop cocaine addiction-like behavior in THC-exposed rats. Indeed, in spite of the initial findings previously mentioned regarding the increased SA of drugs in cannabinoid pre-exposed animals, these studies have not always examined the complex full array of behaviors that are indicative of addiction, especially compulsive seeking or taking (Deroche-Gamonet et al., 2004; Everitt et al., 2018), a cardinal feature of addiction-like behavior typically evaluated using punished CSA procedures (Deroche-Gamonet et al., 2004; Belin et al., 2008). Therefore, the last goal of the present work was to experimentally examine the different features of addiction-like behavior that may be potentiated by ACE.
Our results provide extensive evidence that exposure to THC during adolescence causes sex-dependent changes in reward processing, impulsivity, and specific features of addiction-like behavior, together with gene-expression alterations in the NAcc, providing additional experimental support to the data gathered in clinical and epidemiological studies.
METHODS
Animals and THC Treatment
Subjects were the offspring of Wistar albino rats (35 males and 35 females) from Charles River S.A. (Saint-Germain-sur-l’Arbresle, France) that were mated in our laboratory 2 weeks after their arrival. Different sets of animals, belonging to different litters, were randomly assigned to each experiment, thus minimizing litter effects. The final sample size for each experiment is indicated in the sections below.
Chronic Δ 9-tetrahydrocannabinol (THCPharm, Frankfurt, Germany) treatment took place every other day from postnatal day (PND) 28 to PND 44. THC (3 mg/kg; 1 mL/kg) or its vehicle (kolliphor:ethanol:saline; 1:1:18) were administered i.p. Animals were left undisturbed until PND 90.
All procedures involving laboratory rats were conducted in accordance with the European Union legislation on the protection of animals used for scientific purposes (2010/63/EU Directive) and approved by the Ethics Board of the National University of Distance Learning. Every attempt was made to minimize the pain and discomfort of the experimental animals. See supplementary Methods for more information.
Transcriptomic Studies (Experiment 3)
RNA from the NAcc of P90 rats was extracted and sequenced in a 1 × 75 single-read sequencing run on a NextSeq500 sequencer. Sample sizes as follows: male VEH (n = 4), male THC (n = 4), female VEH (n = 4), and female THC (n = 4).
We used the Chipster analysis suite (Kallio et al., 2011) to perform data processing and analysis.
Differential gene expression analysis was performed using CUFFDIFF with replicates analysis to explore the differences in transcriptomic profiles between factor levels. Gene ontologies and pathways enrichment and overrepresentation were calculated with the online tools and databases of PANTHER Classification System (Mi et al., 2019) for every gene subset obtained in the differential analysis. See supplementary information for detailed methodological information.
Cocaine Addiction–Like Behavior (Experiment 4)
On PND 90, animals from a different batch underwent a single food-reinforced fixed ratio 1 instrumental training session limited to 10 reinforcers. After this, an i.v. polyvinylchloride tubing (0.064 mm i.d.) catheter was implanted into the right jugular vein.
The cocaine SA protocol was carried out in Coulbourn boxes. Cocaine (Alcaliber, Madrid, Spain) infusions (0.5 mg/kg in 100 µL of sterile saline solution) were administered by an electronic pump. The protocol consisted of 6 consecutive phases: (1) acquisition (12 daily 2-hour fixed ratio 1 sessions); (2) motivation for consumption (progressive ratio schedule (Sánchez-Cardoso et al., 2007); (3) rebound consumption: three 2-hour sessions under fixed ratio 1; (4) compulsive (punished) seeking: a single 1-hour session under a fixed ratio 3 schedule in which the animal randomly received an infusion or a 0.5-mA plant shock for 0.5 seconds; (5) extended access: 10 sessions of 6 hours each under fixed ratio 1; and (6) cue-induced reinstatement: 4 sessions of 1 hour each with response-contingent cues (same cues as those used during acquisition) but without drug delivery, occurring after 1, 30, 60, and 100 days of forced abstinence. See supplementary information for detailed methodological information. Initial sample sizes as follows: male VEH (n = 15), male THC (n = 18), female VEH (n = 15), and female THC (n = 15).
Statistical Analysis
In general, for the experiments involving repeated measures, we used a mixed ANOVA with 2 between-subject factors (sex and ACE) and 1 within-subject factor (session or test). For the indices without repeated measures, we used standard 2-way ANOVAs. Significant interactions were followed using simple effects analysis.
RESULTS
Transcriptome Profile in the Shell of the NAcc
After the results obtained in all these reward-related processes, we decided to examine how ACE affected the transcriptome of the NAcc, a key region regulating reward-guided behavior and impulsivity. There were 95 differentially expressed genes (27 upregulated and 68 downregulated) in THC males compared with VEH males and 84 (30 upregulated and 54 downregulated) differentially expressed genes in the females’ comparison. Only 9 of these differentially expressed genes were present in both differential analyses (see Figures 3 and 4). In the males, the categories with higher fold enrichment included biological processes such as “drug transport,” “learning and memory,” and “chemical synaptic transmission” or were restricted to cell compartments such as the axon (see Figure 3). In the females, a completely different set of categories was affected. The ontologies with higher fold enrichment were related to “hormonal activity,” the organization of cellular projections (including the “axoneme”) and the “cytoskeleton” (see Figure 4).
Cocaine Addiction–Like Behavior
Following our behavioral and transcriptomic results, we proceeded to examine cocaine addiction–like behavior. All rats acquired cocaine SA in a similar way (Figure 5B). However, during progressive ratio sessions, a between-subjects analysis of infusions showed a sex × ACE interaction (F1,25 = 5.215; P = .031; η p2 = 0.173), revealing that THC-exposed males had a higher overall cocaine intake than VEH-treated males (F1,25 = 6.197; P = .032; η p2 = 0.382) and that VEH-treated females had a higher cocaine intake than VEH-exposed males (F1,25 = 7.717; P = .018; η p2 = 0.412) during these high-effort conditions (see Figure 5C). Interestingly, this sex difference was not observed among THC-treated rats. We then returned rats to continuous drug access (fixed ratio 1) for 3 days. When we compared the relative increase during these fixed ratio 1 sessions and the average of the last 3 acquisition sessions (also under fixed-ratio 1 schedule), we observe a significant sex × ACE interaction (F1,29 = 7.507; P = .010; η p2 = 0.21) that revealed that THC-exposed females had a higher rebound than VEH-exposed controls (F1,29 = 9.497; P = .004; η p2 = 0.25). We also observed that, among VEH-exposed rats, there was a significant effect of sex (with VEH-exposed males having higher rebound than VEH-exposed females [F1,29 = 5.165; P = .015; η p2 = 0.15]) (Figure 5D).
During the punished seeking test, all rats reduced the number of infusions achieved compared with the last reacquisition session, but there were no effects due to sex or ACE (see supplementary information for additional measures and graphs). After this single session, we allowed the rats to self-administer cocaine for 6 h/d under a fixed-ratio 1 schedule of reinforcement for 10 days. All groups similarly escalated their intake (F1,20 = 4.349; P = .05; η p2 = 0.179) (see Figure S5D). We did not observe a significant effect of ACE on total cocaine intake across sessions across all 10 extended-access sessions (see Figure 5E).
We then withdrew the rats from cocaine and analyzed their (non-reinforced) seeking responses after 1, 30, 60, and 100 days of forced withdrawal. There was a progression of seeking responses increasing from withdrawal day 1 and peaking around withdrawal day 30—reproducing the incubation of seeking phenomenon—statistically evidenced by the significant effect of session (F2.05,43.05 = 6.618; P = .003; η p2 = 0.24). Noteworthy, females showed a more robust seeking behavior (significant effect of sex) (F1,21 = 11.607; P = .003; η p2 = 0.36) (see Figure 5F). We did not obtain a significant sessions × sex × ACE interaction, but the ad hoc analysis of the withdrawal day-30 session showed a sex × ACE interaction (F1,22 = 4.847; P = .038; η p2 = 0.18) with significant simple effects suggesting a significantly lower seeking behavior of THC-exposed females compared with VEH females (F1,22 = 11.924; P = .002; η p2 = 0.35) and also a significantly higher seeking VEH females compared with VEH males (F1,22 = 17.751; P < .000; η p2 = 0.45).
Discussion
We have provided evidence for a causal relationship between the exposure to THC, the main psychoactive component of cannabis, during adolescence and alterations at adulthood in a set of psychological mechanisms related to reward processing, impulsivity, and some features of cocaine addiction–like behavior.
Transcriptome Profile in the Shell of the NAcc
Our RNA-seq study provides, for the first time, an exploration of the sex-dependent differential effects of ACE on the striatal transcriptome. We will focus our discussion on reward processes, response to drugs, and substance use disorders, which is the aim of the present work. However, given the importance of some of the transcriptional signatures obtained for schizophrenia, an important comorbid condition of substance used disorders (Khokhar et al., 2018), we also provide some discussion of the relevance of our findings to this disorder in the supplementary Discussion.
In male rats, the upregulated gene with the highest fold change was Satb2 (SATB homeobox 2), involved in transcription regulation and chromatin remodeling. CB1 receptors are coupled to the regulation of the Ctip2–Satb2 transcriptional regulatory code (Diaz-Alonso et al., 2012), and in so doing, they guide corticospinal motor neuron differentiation. The alteration of the Satb2 gene in the NAcc of our animals could also have developmental consequences in the morphology or function of accumbal neurons as suggested by the ontologies affected by THC treatment. Moreover, Satb2 in the paraventricular thalamus is also sensitive to cocaine-rewarding actions (Salti et al., 2018), so it could be speculated that this upregulated gene in accumbal cells may affect the rewarding actions of cocaine under specific circumstances (such as progressive ratio schedules; see below). Another gene with potential implications for our behavioral results was Notch3 (notch receptor 3), which was downregulated. This gene belongs to the notch signaling pathway that is also involved in brain development (Androutsellis-Theotokis et al., 2006). Interestingly, Notch3 is downregulated in striatal territories in spontaneously hypertensive rats treated with methylphenidate during adolescence (and that further self-administered methylphenidate as adults) (dela Peña et al., 2014), suggesting that this gene is responsive to several pharmacological challenges during adolescence (not just cannabinoids), with dopamine acting as a potential common link (Gottlieb, 2001; Bossong et al., 2009; Wahlstrom et al., 2010), and also that its downregulation may predispose to psychostimulant consumption.
In addition, we found transcriptional and translational alterations in adult animals exposed to adolescent cannabinoids that may influence, in a sex-dependent manner, elements of the dopaminergic signaling pathway and shape drug-related behaviors. In this regard, the transcriptional factor Zinc Finger Homeobox 3, Zfhx3 (downregulated in THC males, upregulated in THC females and associated with neurogenesis and ion-binding Gene Ontologies extracted from the list of differentially expressed genes in the males), is distinctively present in a subtype of D2-expressing neurons of the adult midbrain (Poulin et al., 2014). Thus, this difference may suggest a potential modulation of this specific subtype of D2-expressing neurons in the NAcc shell of rats exposed to THC. Additionally, THC males showed an upregulation of the nuclear receptor gene Nr4a2 involved in behavior and neuron differentiation gene ontologies within the male differentially expressed genes. Nr4a2 can be modulated by neuronal firing and dopamine signaling; moreover, the loss of D2 signaling also contributes to Nr4a2 upregulation (Tseng et al., 2000).
ACE also altered several elements belonging to glutamate, GABA signaling, and other ion channels relevant for the expression of motivated behaviors and drug use. Among the glutamatergic alterations, we find it relevant to highlight the upregulation in both males and females of the solute carrier Slc17a6, which encodes the vesicular glutamate transporter 2(VGlu2) protein involved in glutamate uptake into synaptic vesicles at presynaptic nerve terminals. Noteworthy, dopamine neuronal subtypes express this protein, and the presence of VGlu2 seems to be required for psychostimulant-induced behavioral activation (Birgner et al., 2010).
In addition, there are several relevant changes in the male NAcc shell ion channel expression profile. In this regard, THC produced a protracted downregulation in the male NAcc shell of the solute carrier Slc1a2, a glial transporter that clears glutamate from the synaptic cleft. The expression of this gene is altered by many drugs of abuse (cocaine, amphetamines, nicotine, opioids, ethanol, and cannabinoids), and it has received attention as a potential target for pharmacological interventions in substance use disorders (Roberts-Wolfe and Kalivas, 2015); the voltage-gated potassium channel subunit beta-2 Kcnab2, which is similarly depleted after chronic morphine exposure (Mazei-Robison et al., 2011) and has been involved in motivated behaviors (O’Donovan et al., 2019); the ATP-sensitive inward rectifier potassium channel 10 Kcnj10, also involved in substance use disorders and ethanol preference (Zou et al., 2009); the sodium channel protein type 8 subunit alpha, Scn8a, which plays an important role in regulating excitability in the brain; and the potassium-chloride transporter member 5, Slc12a5, associated with the formation and maturation of glutamatergic and GABAergic synaptic connections (Medina et al., 2014).
In the females, ACE was associated with upregulation of the GABA A Receptor Epsilon Subunit, Gabre. In the context of substance use disorders, rats with a genetic predisposition to alcohol consumption showed a Gabre upregulation (Spence et al., 2018). The hormone activity gene ontology was enriched in the female subset of DEGs, and, noteworthy, the neuropeptides included in this subset may be determining the dopaminergic activity in the NAcc shell of THC-treated females. Thyrotropin-releasing hormone Trh, upregulated by THC in females, participates in energy metabolism and affects different hormonal functions but also enhances dopamine release in the NAcc (Puga et al., 2016). We also detected an upregulation of Agt, which encodes angiotensinogen, the precursor protein of angiotensin I, which is further converted to the peptide angiotensin II. Reductions of angiotensin II, and consequent less activation of the angiotensin II type 1 and type 2 receptors are associated with lower levels of dopamine in the ventral tegmental area and linked to lower alcohol consumption (Maul et al., 2005). Similarly, the cocaine- and amphetamine-regulated transcript seems to exert a neuromodulatory role in the NAcc attenuating dopamine release (Rakovska et al., 2017), and cocaine- and amphetamine-regulated transcript injections into NAcc inhibit the behavioral effects of cocaine (Yu et al., 2017).
We have also detected an upregulation of the Gal gene, encoding the neuropeptide galanin. Noteworthy, galanin has been involved in pathological food consumption and addiction (Gosnell et al., 1986a, 1986b; Sandi et al., 1988). In this regard, an overabundance of galanin has been shown to decrease the sensitivity to amphetamine-induced behavioral effects (Clarke et al., 1988). It is also important to mention that agonists of the galanin receptor can reduce reinstatement of cocaine-seeking (Ogbonmwan et al., 2015) and cocaine-conditioned place preference (Narasimhaiah et al., 2009), which is important considering the sex-specific effects in cocaine SA after ACE (Higuera-Matas et al., 2008). Finally, THC produced a sex-dependent change in the expression of the Cck gene (upregulated in males and downregulated in females treated with THC), which encodes the peptide hormone cholecystokinin. Cholecystokinin (Cck) signaling pathways have been related to food intake, but also with reward and anxiety and even panic (Bradwejn and Vasar, 1995; Rotzinger and Vaccarino, 2003), and have also been studied in the context of drug-related behaviors (Lu et al., 2001, 2002; Wunderlich et al., 2004). Moreover, cocaine behavioral sensitization is accompanied by increasing levels of Cck in the NAcc shell (Beinfeld et al., 2002).
Cocaine Addiction–Like Behavior
In accordance with a previous study (Kononoff et al., 2018), we found no differences in the acquisition of cocaine SA between THC- or vehicle-exposed rats. We have used an intermediate dose of cocaine (0.5 mg/kg), which may explain the divergence between our results and those of Friedman and colleagues (Friedman et al., 2019), who found potentiated cocaine SA with lower doses (0.1 mg/kg) but not with a higher dose (0.32 mg/kg).
ACE was associated with higher intake during progressive ratio sessions in males but not females and effect that may rely on the previously mentioned increase of D1 receptors in the NAcc shell after ACE specifically in the males (Bari and Pierce, 2005; Higuera-Matas et al., 2010). Previous findings showed by Friedman and Kononoff showed unaltered motivation for consumption, although the lower doses (0.1 and 0.32 mg/kg) employed by Friedman et al. (Friedman et al., 2019) and the different timing of the cannabinoid treatment or the cannabinoid agent (WIN 55,512-2) in the study of Kononoff et al. (Kononoff et al., 2018) are 2 probable sources of this divergence. After returning to continuous access to the drug, females showed an increased rebound in consumption compared with the last cocaine SA sessions. This enhanced vulnerability may be relevant for situations of difficult drug access followed by resumption of availability (such as the transition from lock-down in the COVID-19 pandemic to a normal situation).
All the rats diminished their cocaine intake in a similar proportion during compulsive drug taking, ruling out potential changes in compulsivity, although this feature may require further investigation using repeated testing sessions to reveal the effect.
As far as we know, we are the first to explore cocaine-seeking incubation after ACE and to include females in the study. Our data suggesting higher incubation in the females are consistent with previous research showing more robust incubation in females and a higher tendency to reinstate seeking by conditioned cues and drug priming (Lynch and Carroll, 1999; Kerstetter et al., 2008; Nicolas et al., 2019). Previous ACE studies with cocaine showed that adult mice with an adolescent exposition to WIN55,212-2 were less susceptible to the anxiogenic effects of cocaine abstinence (Aguilar et al., 2017), suggesting a potential mechanism for our effects that needs to be further explored, especially concerning its potential sex-specific nature.
CAVEATS AND CONCLUDING REMARKS
An important caveat to consider in this work is the fact that we have opted for a passive i.p. administration route. We chose this route of administration to ensure a homogenous and comparable exposure with the THC across subjects, something that may have been difficult to achieve using i.v. SA procedures or operant vapor SA protocols; however, this is a limitation that should be kept in mind when considering the general translatability of our results
In this work, we have shown that exposure to THC during adolescence profoundly affects the transcriptomic programs of the NAcc and, concomitantly, modulates the influence of rewards and reward-related cues on behavior with a subtle and transient impact on a specific form impulsiveness. Exposure to the main active psychoactive component of cannabis during adolescence also affects certain aspects of addiction-like behavior differentially in males and females and may protect females from the incubation of seeking. These results should be taken into consideration for the ongoing debate about the validity of the Gateway Hypothesis, for tailoring sex-specific treatment approaches for cocaine use disorder depending on previous cannabis consumption during adolescence and, in general, for the evaluation of the long-term consequences of cannabis use by adolescents, an especially vulnerable population.
Supplementary Material
Acknowledgments
We thank Alberto Marcos, Rosa Ferrado Luis Carrillo, Luis Troca, and Gonzalo Moreno for their excellent technical assistance. This manuscript was prepared during the COVID-19 pandemic and is dedicated to those who helped our society. We honor them all.
This work was supported by the Spanish Ministry of Health, Social Services and Equality (Network of Addictive Disorders – Project no: RTA-RD16/0017/0022 of the Institute of Health Carlos III to E.A. and Plan Nacional Sobre Drogas, Project no: 2016I073 to E.A. and 2017I042 to A.H.-M.); Ministry of Science (PID2019-104523RB-I00 to A.-H.M. and PID2019-111594RB-100 to E.A.); UNED (Plan for the Promotion of Research to E.A. and A.H.-M.); the European Union (Project no: JUST/2017/AG-DRUG-806996-JUSTSO) and the BBVA Foundation (2017 Leonardo Grant for Researchers and Cultural Creators to A.H.-M.). These agencies funded the study but had no further role in the study design; the collection, analysis, and interpretation of data; the writing of the report; or the decision to submit the paper for publication. J.O. received funding from Instituto de Salud Carlos III and is now an alumnus of the graduate program of the International Graduate School of UNED. D.-R.M. received a predoctoral fellowship granted by UNED and M.U. received a predoctoral fellowship awarded by the Ministry of Science and Innovation (BES-2011-043814).
Interest Statement
The authors have no conflict of interest that may affect the results or conclusions reported in this work.