Skip to main content
Advertisement
  • Loading metrics

Genetics of cocaine and methamphetamine consumption and preference in Drosophila melanogaster

  • Chad A. Highfill,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft

    Current address: Center for Human Genetics and Department of Genetics and Biochemistry, Clemson University, Greenwood, SC, United States of America.

    Affiliation Department of Biological Sciences, W. M. Keck Center for Behavioral Biology, and Program in Genetics, North Carolina State University, Raleigh, NC, United States of America

  • Brandon M. Baker,

    Roles Conceptualization, Data curation, Investigation, Writing – original draft

    Current address: Center for Human Genetics and Department of Genetics and Biochemistry, Clemson University, Greenwood, SC, United States of America.

    Affiliation Department of Biological Sciences, W. M. Keck Center for Behavioral Biology, and Program in Genetics, North Carolina State University, Raleigh, NC, United States of America

  • Stephenie D. Stevens,

    Roles Data curation

    Affiliation Department of Biological Sciences, W. M. Keck Center for Behavioral Biology, and Program in Genetics, North Carolina State University, Raleigh, NC, United States of America

  • Robert R. H. Anholt,

    Roles Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    Current address: Center for Human Genetics and Department of Genetics and Biochemistry, Clemson University, Greenwood, SC, United States of America.

    Affiliation Department of Biological Sciences, W. M. Keck Center for Behavioral Biology, and Program in Genetics, North Carolina State University, Raleigh, NC, United States of America

  • Trudy F. C. Mackay

    Roles Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    tmackay@clemson.edu

    Current address: Center for Human Genetics and Department of Genetics and Biochemistry, Clemson University, Greenwood, SC, United States of America.

    Affiliation Department of Biological Sciences, W. M. Keck Center for Behavioral Biology, and Program in Genetics, North Carolina State University, Raleigh, NC, United States of America

Abstract

Illicit use of psychostimulants, such as cocaine and methamphetamine, constitutes a significant public health problem. Whereas neural mechanisms that mediate the effects of these drugs are well-characterized, genetic factors that account for individual variation in susceptibility to substance abuse and addiction remain largely unknown. Drosophila melanogaster can serve as a translational model for studies on substance abuse, since flies have a dopamine transporter that can bind cocaine and methamphetamine, and exposure to these compounds elicits effects similar to those observed in people, suggesting conserved evolutionary mechanisms underlying drug responses. Here, we used the D. melanogaster Genetic Reference Panel to investigate the genetic basis for variation in psychostimulant drug consumption, to determine whether similar or distinct genetic networks underlie variation in consumption of cocaine and methamphetamine, and to assess the extent of sexual dimorphism and effect of genetic context on variation in voluntary drug consumption. Quantification of natural genetic variation in voluntary consumption, preference, and change in consumption and preference over time for cocaine and methamphetamine uncovered significant genetic variation for all traits, including sex-, exposure- and drug-specific genetic variation. Genome wide association analyses identified both shared and drug-specific candidate genes, which could be integrated in genetic interaction networks. We assessed the effects of ubiquitous RNA interference (RNAi) on consumption behaviors for 34 candidate genes: all affected at least one behavior. Finally, we utilized RNAi knockdown in the nervous system to implicate dopaminergic neurons and the mushroom bodies as part of the neural circuitry underlying experience-dependent development of drug preference.

Author summary

Illicit use of cocaine and methamphetamine is a major public health problem. Whereas the neurological effects of these drugs are well characterized, it remains challenging to determine genetic risk factors for substance abuse in human populations. The fruit fly, Drosophila melanogaster, presents an excellent model for identifying evolutionarily conserved genes that affect drug consumption, since genetic background and exposure can be controlled precisely. We took advantage of natural variation in a panel of inbred wild derived fly lines with complete genome sequences to assess the extent of genetic variation among these lines for voluntary consumption of cocaine and methamphetamine and to explore whether some genetic backgrounds might show experience-dependent development of drug preference. The drug consumption traits were highly variable among the lines with strong sex-, drug- and exposure time-specific components. We identified candidate genes and gene networks associated with variation in consumption of cocaine and methamphetamine and development of drug preference. Using tissue-specific suppression of gene expression, we were able to functionally implicate candidate genes that affected at least one consumption trait in at least one drug and sex. In humans, the mesolimbic dopaminergic projection plays a role in drug addiction. We asked whether in Drosophila the mushroom bodies could play an analogous role, as they are integrative brain centers associated with experience-dependent learning. Indeed, our results suggest that variation in consumption and development of preference for both cocaine and methamphetamine is mediated, at least in part, through a neural network that comprises dopaminergic projections to the mushroom bodies.

Introduction

Illicit use of cocaine and methamphetamine constitutes a significant public health problem that incurs great socioeconomic costs in the United States and worldwide [13]. Cocaine and the amphetamine class of drugs are potent central nervous system stimulants that act by raising synaptic concentrations of biogenic amines. Cocaine inhibits neurotransmitter reuptake at dopaminergic, serotonergic and noradrenergic synapses [4,5]. Amphetamine increases neurotransmission by promoting the release of dopamine from presynaptic vesicles through its actions on the vesicular monoamine transporter and subsequent reverse flux of dopamine via the dopamine transporter and through the plasma membrane into the synaptic cleft [6,7].

Amphetamine and methamphetamine are used clinically to treat attention deficit hyperactivity disorder and narcolepsy. Long term use of these compounds, however, can lead to addiction, and ultimately death [8]. The addictive properties of these drugs are mediated through the dopaminergic mesolimbic reward pathway, which projects from the ventral tegmental area via the nucleus accumbens to prefrontal cortex [9]. Although most studies on psychostimulants focus on addiction, addiction represents only one facet of the diverse organismal effects that result from psychostimulant drug abuse. These drugs exert a wide range of physiological and behavioral effects, including suppression of appetite, which can result in malnutrition, and severe cardiovascular, respiratory and renal disorders. Use of cocaine and amphetamine can also cause mental disorders, including paranoia, anxiety, and psychosis [10,11].

Susceptibility to the effects of cocaine and methamphetamine is likely to vary among individuals and be determined both by environmental and genetic factors. However, there is limited information regarding the genetic basis of susceptibility to the effects of these drugs in human populations [12]. Twin and adoption studies have focused primarily on alcohol abuse and illicit drugs, such as cannabis, with heritability estimates ranging from ~30–70% [13,14]. Most studies on psychostimulant addiction to date have centered on candidate genes associated with neurotransmission in the mesolimbic projection [12], and many of these are inconclusive or contradictory. For example, some studies reported that alleles of the dopamine D2 receptor were associated with substance abuse [1518], whereas others did not replicate this finding [1924]. Similar contradictory results have been obtained for association analyses between polymorphisms in the dopamine transporter gene and cocaine-related phenotypes [2428]. These contradictory findings may be due in part to failure to account for multiple testing or population structure [29]. However, genetic studies of substance abuse and addiction in human populations are challenging due to diverse social conditions and physical environments, confounding factors with comorbid conditions such as alcoholism or psychiatric disorders, and difficulty to recruit large numbers of study subjects due to criminalization.

Drosophila melanogaster is an excellent model for identifying genes that affect drug consumption behaviors since both the genetic background and environment, including exposure to drugs, can be controlled precisely. These results have translational potential since 75% of disease-causing genes in humans have a fly ortholog [30]. High resolution X-ray crystallography has shown that the D. melanogaster dopamine transporter has a central conformationally pliable binding site that can accommodate cocaine, methamphetamine and their closely related analogues [31]. Similar to its effects in humans, methamphetamine suppresses sleep, causes arousal and suppresses food intake in flies [3234]. In addition, cocaine and amphetamine exert quantifiable locomotor effects in flies [3541]. Thus, despite profound differences between the neuroanatomical organization of the fly and vertebrate brains, it is likely that behavioral and physiological effects of methamphetamine and cocaine are mediated, at least in part, by analogous mechanisms.

Here, we used the inbred, sequenced lines of the D. melanogaster Genetic Reference Panel (DGRP [42,43]) to investigate the genetic basis for variation in psychostimulant drug consumption. We used a four-capillary Capillary Feeding (CAFE) assay [4446] to quantify voluntary consumption, preference and change of consumption and preference over time for cocaine and methamphetamine. Since cocaine and methamphetamine both target dopaminergic synaptic transmission, but through different mechanisms, we asked to what extent genetic networks that underlie variation in consumption of cocaine and methamphetamine incorporate the same or different genes. We also sought to determine the extent of sexual dimorphism for naïve and experience-dependent voluntary drug intake. In addition, we asked how much variation in voluntary drug consumption exists among different DGRP lines and what fraction of that variation is accounted for by genetic variation. We showed that there is naturally occurring genetic variation for all drug consumption traits with strong sex-, drug- and exposure time-specific components. We performed genome wide association (GWA) analyses to identify candidate genes associated with the drug consumption behaviors that could be mapped to a genetic interaction network. We tested the effects of RNAi mediated suppression of gene expression [47] on all consumption behaviors for 34 candidate genes and found that all affected at least one behavior in at least one drug and sex. Finally, we used RNAi to suppress gene expression in neurons, glia, the mushroom bodies and dopaminergic neurons in a subset of genes and showed that innate preference and the development of preference for psychostimulant drugs involves dopaminergic neurons and the mushroom bodies, neural elements associated with experience-dependent modulation of behavior.

Results

Quantitative genetic analysis of drug consumption behaviors in the DGRP

We used a four-capillary CAFE assay [4446] to enable flies to choose to consume either sucrose or sucrose supplemented with 0.2 mg/ml cocaine (or 0.5 mg/ml methamphetamine), analogous to the two-bottle choice assay used in rodent studies [48] (Fig 1). We quantified consumption for three consecutive days for males and females from each of 46 DGRP lines that were unrelated, free of chromosomal inversions, and free of infection with the endosymbiont Wolbachia pipientis [43; S1 Table]. These data enabled us to assess whether there is naturally occurring genetic variation in this population for naïve consumption of each solution and preference, and change of consumption and preference upon repeated exposures (i.e., experience-dependent modification of behavior).

thumbnail
Fig 1. Consumption and preference assay.

(A) Cartoon illustrating the four capillary CAFE assay. Each of the three exposures consists of an 18 hour feeding trial with sucrose or drug + sucrose, followed by 6 hours recovery with standard culture medium. (B) Positions of capillaries with the two solutions (indicated by red and yellow).

https://doi.org/10.1371/journal.pgen.1007834.g001

We performed four-way mixed model analyses of variance (ANOVA) to partition variation in consumption between DGRP lines, males and females, drug vs. sucrose, and the three exposures. All main effects were significant for both drugs (Table 1), indicating genetic variation for consumption, difference between amount of sucrose and drug consumed, sexual dimorphism, and experience-dependent modulation of behavior. We are most interested in the two- and three-way interaction terms involving Line, as they indicate genetic variation in sexual dimorphism (L×X), change of consumption between exposures (L×E), preference for sucrose or drug solution (L×S), and change of preference for sucrose or drug between exposures (L×E×S). With the exception of L×S, these interaction terms were significant for both the cocaine and methamphetamine analyses (Table 1).

thumbnail
Table 1. Analyses of variance of consumption measured over three exposures.

Exposure, Sex, Solution, and their interaction are fixed effects, the rest are random. E: Exposure; X: Sex; S: Solution; L: DGRP Line; ε: residual; df: degrees of freedom; MS: Type III mean squares; F: F-ratio test; P: P-value; σ2: variance component estimate; SE: standard error; H2: Broad sense heritability. Significant P-values are highlighted in bold font.

https://doi.org/10.1371/journal.pgen.1007834.t001

We next performed reduced ANOVA models to quantify broad sense heritabilities (H2) for consumption and change in consumption traits (S2 Table). We found significant genetic variation in consumption of both drugs and sucrose alone within each sex and exposure, with H2 ranging between 0.20 and 0.38 for cocaine consumption and between 0.22 and 0.30 for methamphetamine consumption (Fig 2, S2 Table). Further, there was significant genetic variation for the change in consumption of sucrose alone or drug in both sexes between the third and first exposures, with H2 ranging between 0.14 and 0.18 for cocaine and between 0.17 and 0.22 for methamphetamine (Fig 2, S2 Table). Thus, there is genetic variation for both consumption and experience-dependent consumption of both drugs and sucrose alone in the DGRP.

thumbnail
Fig 2. Variation in drug consumption among 46 DGRP lines.

(A) Initial exposure. Lines are from lowest to highest consumption in females. (B) Third exposure. The line order is the same as in (A). (C) Change in consumption between exposures 3 and 1. Positive values indicate increased drug consumption in Exposure 3. The line order is the same as in (A). Pink denotes females, blue indicates males, and purple is overlap of both sexes. Error bars are ± 1SD.

https://doi.org/10.1371/journal.pgen.1007834.g002

Finally, we defined preference in two ways: as the difference between amount of drug and sucrose alone consumed (Preference A), and as this difference scaled by the total amount of both solutions consumed (Preference B). Variation in Preference A is the statistical interpretation of the line by solution interaction; Preference B is the metric commonly used in rodent studies [48]. Preference values of 0 indicate equal consumption of sucrose alone and sucrose containing drug; values > 0 represent preference for the drug and values < 0 indicate drug avoidance. Both preference metrics were significantly genetically variable for each sex and exposure for cocaine, with H2 ranging from 0.06–0.16; while for methamphetamine, both preference metrics were significantly genetically variable in females for all exposures (H2 from 0.05–0.18) and for males in the second and third exposures (H2 from 0.08–0.11) (S2 Table). For cocaine, the difference in Preference A between exposures 3 and 1 was significant only in females (H2 = 0.11) while the difference in Preference B was significant for females (H2 = 0.13) and males (H2 = 0.05). For methamphetamine, the difference in Preference A was significant in males (H2 = 0.04) and the difference in Preference B was significant in females (H2 = 0.04) (S2 Table). Thus, there is genetic variation for both innate drug preference and experience-dependent drug preference in the DGRP.

The heritabilities of consumption traits are low, as is typical for behavioral traits, indicating that environmental factors, including previous experience, predominantly contribute to the observed phenotypic variation. The advantage of performing multiple replicate measurements of each DGRP line is that the broad sense heritabilities of line means (S3 Table) used in the GWA analyses (see below) are much greater than heritabilities based on individual vials (S2 Table).

We computed the genetic and phenotypic correlations between males and females for the consumption behaviors, between exposures for consumption and preference, and between solutions (S4 Table). Cross-sex genetic correlations for consumption tended to decrease with the number of exposures for both cocaine and methamphetamine, suggesting that the experience-dependent modification of consumption is sex-specific. Consumption of drugs and sucrose is highly correlated across the three exposures (albeit significantly different from unity), while the correlations of drug preference across exposures are low to moderate for both cocaine and methamphetamine in both sexes. Although the consumption of drugs and sucrose for cocaine and methamphetamine are genetically and phenotypically correlated in both sexes, preference for the two drugs is not significantly correlated. Finally, Preference A and Preference B within each exposure are nearly perfectly correlated, as expected since the difference in consumption is in both metrics.

In summary, we found that there is extensive genetic variation in consumption and preference as well as change in consumption and preference with repeated exposures for both cocaine and methamphetamine across different genetic backgrounds, and that genetic variation for these traits has significant sex- and drug-specific components.

Genome wide association analyses of drug consumption in the DGRP

Our quantitative genetic analyses of consumption in the DGRP indicate that there is genetic variation for all traits assessed, and that the traits have a complex correlation structure indicating partially common and partially distinct genetic bases. Therefore, we performed GWA analyses for 12 traits (drug and sucrose consumption exposure 1, drug and sucrose consumption exposure 3, change in drug and sucrose consumption, Preference A exposure 1, Preference A exposure 3, Preference B exposure 1, Preference B exposure 3, change in Preference A, and change in Preference B) for cocaine and methamphetamine, separately for males and females. We performed association tests for 1,891,456 DNA sequence variants present in the 46 DGRP lines with minor allele frequencies greater than 0.05 [43].

At a lenient significance threshold of P < 5 x 10−5, we identified 1,441 polymorphisms in or near (within 1 kb of the start and end of the gene body) 725 genes for all consumption behaviors related to cocaine, and 1,413 polymorphisms in or near 774 genes for methamphetamine exposure (S5 Table). The majority of these variants had sex-specific effects. A total of 40 variants and 141 genes overlapped between cocaine and methamphetamine. The variants in or near genes implicate candidate genes affecting consumption behaviors, while the intergenic variants could potentially contain regulatory motifs for transcription factor-binding sites or chromatin structure regulating these traits. Only two variants are formally significant following a Bonferroni correction for multiple tests (0.05/1,891,456 = P < 2.64 × 10−8). 2L_10179155_SNP is located within an intronic region in CG44153 and affects experience-dependent development of methamphetamine preference in both sexes. Its human homolog ADGRB3 encodes a G-protein coupled receptor, which contributes to the formation and maintenance of excitatory synapses [49] and has been implicated in GWA studies on human addiction [50]. 3R_27215016_SNP is a synonymous SNP in the coding sequence of CG1607 and affects naïve consumption of sucrose. CG1607 encodes an amino acid transmembrane transporter. One of its human orthologs, SLC7A5, is an amino acid transporter, mutations in which are associated with autism spectrum disorder and defects in motor coordination [51].

While not formally significant, we identified genes previously associated with cocaine-related behaviors (Bx [Lmo], loco, Tao) and ethanol-related behaviors (Bx, DopR, Egfr, hppy, Tao, Tbh) [52] in D. melanogaster. In addition, the genes implicated by the GWA analyses are enriched for multiple gene ontology (GO) categories and pathways [53,54] at a false discovery rate < 0.05 (S5 Table). GO terms involved in nervous system development and function were among the most highly enriched, consistent with the known neurobiological mechanisms of action of these drugs. Finally, we note that ~ 70% of the candidate genes from the GWA analyses have human orthologs, and many of these genes have previously been associated with cocaine or methamphetamine abuse in humans or with behaviors associated with intake and response to various psychoactive substances (alcohol, cannabis, nicotine, opioids) in humans as well as zebrafish, mouse and rat models (S6 Table). This suggests that cocaine and methamphetamine exert their effects in flies and humans through evolutionarily conserved neural mechanisms.

These results suggest a highly polygenic architecture for variation in consumption and drug preference, and that the genetic underpinnings for variation in consumption or preference are both shared and distinct for cocaine and methamphetamine, consistent with the quantitative genetic analyses.

A genetic interaction network for consumption behaviors

We next asked whether the genes we identified in the GWA analyses belonged to a known genetic interaction network. Since the consumption behaviors are highly inter-correlated, we queried whether all 1,358 candidate genes from the GWA analyses for both cocaine and methamphetamine combined could be clustered into significant sub-networks based on curated genetic interactions in Drosophila. If we do not allow any missing genes, we find a significant (P = 9.99 × 10−4) network of 81 candidate genes (Fig 3, S7 Table), most of which (88.9%) are predicted to have human orthologs [55].

thumbnail
Fig 3. Significant genetic interaction network of genes identified in the GWA analyses for all cocaine and methamphetamine related traits combined.

Borders indicate the strength of the evidence for a human ortholog. Black: DIOPT score < 3; Blue: DIOPT score 3–6; Green: DIOPT score 7–9; Orange: DIOPT score 10–12; Red: DIOPT score 13–15. See S7 Table for the complete list of human orthologs and their DIOPT scores. Grey boxes have effects on at least one drug-seeking behavior from RNAi knockdown of gene expression.

https://doi.org/10.1371/journal.pgen.1007834.g003

We performed enrichment analyses [53,54] to gain insight in the biological context for genes in the network using a false discovery rate < 0.05. Surprisingly, many canonical signaling pathways are highly enriched, including the Wingless (Wnt), Cadherin, Cholecystokinin Receptor (CCKR), Transforming Growth factor beta (TGF), and Fibroblast Growth Factor (FGF) signaling pathways. Concomitantly, we find high enrichment of molecular function GO terms associated with regulation of transcription and DNA and protein binding, and biological function GO terms associated with development (including the development of the nervous system; S7 Table). These results suggest that naturally occurring genetic variation in nervous system development is associated with variation in propensity to consume psychostimulant drugs. Furthermore, our results indicate that natural variants in key genes regulating all aspects of fly development and function can be associated with variation in drug consumption behaviors.

Functional evaluation of candidate genes

We used RNA interference (RNAi) to functionally test whether reduced expression of candidate genes implicated by the GWA analyses affect consumption phenotypes. We selected 34 candidate genes for RNAi mediated suppression of gene expression from the set of GWA analyses (S5 Table). A total of nine of the candidate genes were in the network; the others were chosen based on gene expression in the nervous system and their known role in nervous system function, as well as belonging to enriched pathways and gene ontology categories. We measured consumption of cocaine and sucrose (S8 Table) and methamphetamine and sucrose (S9 Table) for three consecutive days, separately for males and females, for each of the RNAi and control genotypes, exactly as described for the DGRP lines.

We performed three-way fixed effect ANOVAs for each UAS-RNAi and control genotype, separately for males and females (S10 and S11 Tables). The main effects in these models are genotype (L, RNAi and control), solution (S, sucrose and drug) and exposure (E, first and third). A significant L effect denotes a difference in overall consumption between the RNAi and control genotypes; a significant S effect indicates a difference in preference between sucrose alone and sucrose with drug; and a significant E effect indicates a difference in consumption between exposures 1 and 3. Significant L×S and L×E interaction terms denote, respectively, a difference in preference between the RNAi and control genotypes, and a difference in consumption between exposures 1 and 3 between the two genotypes. A significant L×S×E interaction indicates a change in preference with repeated exposure between the RNAi and control genotypes. We are most interested in the main effect of genotype and interactions with genotype; i.e., consumption, preference, change of consumption and change of preference.

First, we used a weak ubiquitous GAL4 driver crossed to all 34 UAS-RNAi genotypes and their respective controls. All candidate genes had a significant (P < 0.05) effect on at least one of the consumption traits in at least one drug or sex combination. A total of 22 (25) genes affected consumption of cocaine (methamphetamine), 21 (23) affected a change of consumption with exposure to cocaine (methamphetamine), 16 (10) affected cocaine (methamphetamine) preference, and 11 (11) affected a change in cocaine (methamphetamine) preference with exposure in males and/or females (S10 and S11 Tables, S1S3 Figs). There were pronounced sex- and drug-specific effects for all drug-related traits. The majority of RNAi genotypes showed reduced consumption of cocaine and/or methamphetamine compared to their controls, dependent on exposure and sex. If consumption is positively associated with gene expression, this suggests that the products of these genes contribute to drug consumption. On the other hand, several RNAi constructs caused increased drug consumption, suggesting that naturally occurring variants that decrease expression of these genes could predispose to drug preference. Finally, several RNAi-targeted genes exhibit a relative increase or decrease in drug consumption compared to the control at the third exposure, indicating experience-dependent change in preference.

To extend and refine our RNAi analysis, we next selected 10 genes (Dop1R1, Ect4, ed, mld, msi, Oct-TyrR, olf413, Snoo, Vha100-1, wmd) from among those that showed phenotypic effects when targeted by RNAi under the ubiquitous driver and which have known effects on the nervous system. We assessed functional effects of these genes on consumption traits when their corresponding RNAi constructs were expressed under the control of the neuronal-specific elav driver or glial-specific repo driver. All of these genes had a significant (P < 0.05) effect on at least one of the consumption traits in at least one drug or sex combination under the elav driver, and all but Snoo had significant effects on at least one of the consumption traits in at least one drug or sex combination under the repo driver. With neuronal-specific suppression of gene expression, 9 (10) genes affected consumption of cocaine (methamphetamine), 6 (7) affected a change in consumption with exposure to cocaine (methamphetamine), 2 (7) affected cocaine (methamphetamine) preference, and 3 (6) affected a change in cocaine (methamphetamine) preference with exposure in males and/or females (S10 and S11 Tables, Fig 4, S4 Fig). With glia-specific suppression of gene expression, 4 (7) genes affected consumption of cocaine (methamphetamine), 7 (6) affected a change in consumption with exposure to cocaine (methamphetamine), 3 (0) affected cocaine (methamphetamine) preference, and 2 (3) affected a change in cocaine (methamphetamine) preference with exposure in males and/or females (S10 and S11 Tables, Figs 4 and 5, S4 Fig). These effects were largely sex-, drug- and driver-specific. We infer from these results that variation in gene expression in both neurons and glia contributes to phenotypic variation in drug intake behaviors.

thumbnail
Fig 4. Differences in cocaine preference and change in cocaine preference between the third and first exposures between RNAi and control genotypes.

(A) Female Preference A. (B) Male Preference A. (C) Female change of Preference A. (D) Male change of Preference A. Red, black, blue, and green bars denote elav-GAL4, repo-GAL4, 201Y-GAL4 and TH-GAL4 drivers, respectively. Asterisks represent significant L×S terms (A, B) or significant L×S×E terms from the full ANOVA models. Exact P-values are given in S10 Table.

https://doi.org/10.1371/journal.pgen.1007834.g004

thumbnail
Fig 5. Differences in methamphetamine preference and change in methamphetamine preference between the third and first exposures between RNAi and control genotypes.

(A) Female Preference A. (B) Male Preference A. (C) Female change of Preference A. (D) Male change of Preference A. Red, black, blue, and green bars denote elav-GAL4, repo-GAL4, 201Y-GAL4 and TH-GAL4 drivers, respectively. Asterisks represent significant L×S terms (A, B) or significant L×S×E terms from the full ANOVA models. Exact P-values are given in S11 Table.

https://doi.org/10.1371/journal.pgen.1007834.g005

In humans, the mesolimbic dopaminergic projection plays a role in drug addiction. In Drosophila, the mushroom bodies could play an analogous role, as they are integrative centers in the fly brain associated with experience-dependent learning [56,57], dependent on dopaminergic input. To test whether the mushroom bodies and dopaminergic projection neurons could serve as neural substrates that contribute to variation in drug consumption or preference, we focused on four genes (Dop1R1, ed, msi, Snoo,) that showed robust phenotypic effects when targeted with a corresponding elav-driven RNAi. Knockdown of all four genes with a mushroom body specific driver resulted in significant effects on consumption of cocaine and/or methamphetamine for at least one drug and sex combination (S10 and S11 Tables, Figs 4 and 5, S5 Fig). Expression of RNAi in mushroom bodies affected change in consumption of cocaine and methamphetamine for Dop1R1; cocaine preference and change of methamphetamine preference for ed; change in consumption of cocaine for msi; and cocaine and methamphetamine preference, cocaine preference, change of cocaine preference and change of consumption of methamphetamine for Snoo. Expression of RNAi in dopaminergic neurons affected change of consumption of cocaine and change in methamphetamine preference for Dop1R1; consumption for cocaine and methamphetamine, change of consumption of methamphetamine and cocaine preference for ed; consumption of cocaine and methamphetamine, change of consumption of cocaine, and cocaine preference for msi; and all four traits for Snoo (S10 and S11 Tables, Figs 4 and 5, S5 Fig). These effects are largely sex-, drug- and driver-specific.

These results suggest that, despite differences in the genetic underpinnings of susceptibility to cocaine and methamphetamine, phenotypic manifestation of genetic variation in consumption and development of preference for both drugs is channeled in part through a neural network that comprises dopaminergic projections to the mushroom bodies.

Discussion

Although studies using mice [58,59], rats [60,61], primates [62] and humans [63] provide important information about the cellular, developmental, physiological, and behavioral effects of psychostimulants, these systems are less suited to dissecting the relationship between naturally occurring genetic variation and phenotypic variation in individual susceptibility to drug consumption and/or preference. Here, we show that D. melanogaster harbors substantial naturally occurring variation for all consumption-related behaviors, including experience-dependent change in consumption, innate drug preference and experience-dependent change in preference, under conditions where we can obtain replicated measurements of consumption for each genotype in a choice assay performed over three successive days under controlled environmental conditions. We show that genetic variation for consumption and preference metrics is both shared between males and females and the different exposures, but is also sex-, exposure- and drug-specific. Sex differences in drug self-administration and addiction have also been shown in humans and mammalian animal models [6472].

The Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-V) defines 11 criteria for substance use disorder in humans, all related to continuing to use of the substance despite adverse social and physiological effects and the development of tolerance with repeated exposure. The DSM-V also recognizes that there is individual variability of unknown etiology for the propensity both to experiment with psychostimulants and to develop symptoms of substance abuse following initial exposure. Previous studies of effects of cocaine [35,3739,7376] and methamphetamine [77] in Drosophila examined mutations and pharmacological interventions using locomotor-based assays, clearly demonstrating an adverse effect of these substances. However, previous Drosophila studies have not assessed naturally occurring variation in drug self-administration and change in this behavior on repeated exposure, which may better model the genetic basis of individual susceptibility–or resistance–to substance abuse and the development of tolerance (increased drug preference over time).

To begin to understand the nature of the genetic basis for variation in drug consumption and preference, we performed GWA analyses for all consumption traits, separately for cocaine and methamphetamine, using 1,891,456 DNA sequence variants present in the 46 DGRP lines with minor allele frequencies greater than 0.05 [43]. We identified 1,358 unique candidate genes using a lenient significance threshold of P < 5 x 10−5. We hypothesized that these candidate genes would be enriched for true positive associations despite the low power of the GWA analyses and that choosing genes for functional evaluation from this list would be more productive than choosing genes at random. Observations supporting this hypothesis are that mutations in several candidate genes have previously been shown to affect cocaine or ethanol-related phenotypes in Drosophila [52], that the candidate genes are highly enriched for GO terms involved in the development and function of the nervous system, and that 81 candidate genes can be assembled into a known genetic interaction network (Fig 3), which is highly unlikely (P = 9.9 x 10−3) to occur by chance. The candidate genes in the significant genetic interaction network are enriched for several canonical signaling pathways as well as all aspects of development, including nervous system development. These observations suggest that subtle genetic variation in nervous system development is associated with variation in propensity for consumption of psychostimulant drugs. Nearly 90% of the genes in the network have human orthologs and are candidates for future translational studies.

We selected nine candidate genes in the significant genetic network and 25 additional candidate genes to assess whether RNAi reduction using a weak ubiquitous GAL4 driver affected consumption traits, using the same experimental design as for the DGRP lines. All of these genes affected at least one consumption trait/sex/drug. However, there is considerable variation in the effects of different drivers on consumption, preference and change in preference for cocaine and methamphetamine, which likely reflects variation in the effects of RNA interference on different neural elements of a complex integrated neural circuitry. Indeed, several candidate genes, functionally implicated by RNAi, are associated with neural development and represent several early developmental signaling pathways. Snoo has been identified as a negative regulator of the decapentaplegic signaling pathway [78,79] and has been implicated in dendritic patterning [80]. Echinoid, the gene product of ed, is an immunoglobulin domain containing membrane protein of adherens junctions that interacts with multiple developmental signaling pathways, including Egfr, Notch and Hippo signaling [8183]. Musashi, encoded by msi, is a neural RNA binding protein that interacts with Notch signaling to determine cell fate [84]. RNAi targeting of expression of these genes under MB-GAL4 or TH-GAL4 drivers show different effects on consumption, change in consumption, preference and change in preference for the two drugs (S5 Fig).

Among the functionally validated candidate genes, Oct-TyrR and Dop1R1 are of special interest. Oct-TyrR encodes an octopamine-tyramine receptor expressed in mushroom bodies [85], and Dop1R1, which encodes a dopamine receptor enriched in the mushroom bodies, has previously been implicated in aversive and appetitive conditioning [86], innate courtship behavior [87] and sleep-wake arousal [88]. Loss-of-function mutations of Dop1R1 increase sleep and these effects are reversed by administration of cocaine [88]. Octopamine and tyramine act on astrocytes via the Oct-Tyr1 receptor and this activation of astrocytes can in turn modulate dopaminergic neurons [89]. Thus, we can hypothesize that combinations of octopaminergic and dopaminergic signaling in the mushroom bodies can modulate drug consumption and/or experience-dependent changes in consumption or preference following repeated exposure to cocaine or methamphetamine.

Finally, genes which were functionally validated with RNAi represent evolutionarily conserved processes. Future studies can assess whether their human counterparts play a role in variation in susceptibility to psychostimulant drug use in human populations.

Materials and methods

Drosophila stocks

The DGRP, UAS-RNAi and GAL4 driver lines used are listed in S12 Table. The DGRP lines are maintained in the Mackay laboratory. RNAi lines [47] were obtained from the Vienna Drosophila Resource Center and the GAL4 driver lines from the Bloomington, Indiana Drosophila stock center. All lines were maintained on standard cornmeal/yeast/molasses medium at 25°C on a 12 hour light/dark cycle with constant humidity of 50%.

Consumption assay

We used a four-capillary Capillary Feeder (CAFE) assay [4446] to measure drug consumption. Briefly, five 3–5 day old flies per genotype/sex were anesthetized using CO2 and placed on cornmeal/yeast/molasses/agar medium one day prior to the assay. Flies were transferred without anesthesia 45 minutes prior to the assay to vials containing 4-5ml of 1.5% agar (Sigma Aldrich). Two capillaries (VWR International: 12.7 cm long, 5 μl total volume) containing 4% sucrose (Sigma Aldrich) + 1% yeast (Fisher Scientific) or 4% sucrose + 1% yeast + drug, with a mineral oil (Sigma Aldrich) overlay (to minimize evaporation), were inserted in the top of each vial. Cocaine and methamphetamine were obtained from the National Institute on Drug Abuse under Drug Enforcement Administration license RA0443159. Flies were allowed to feed for 16–18 hours with the vials placed in an enclosed plastic chamber wrapped in a plastic bag under a 12 hour light/dark cycle with constant humidity of 50%. For each experiment, an identical set of vials without flies was included in each chamber to determine evaporation loss. The capillaries were then removed and the volume of food consumed (1 mm = 0.067 μl) in each calculated as described previously [90]. The capillaries were replaced with a Drosophila activity monitor (DAM) tube (TriKinetics, Inc. Waltham, MA) containing standard cornmeal/yeast/molasses medium for a recovery period of 4–6 hours. The assay was performed on three consecutive days for each vial of flies. The order in which the assays were performed was preserved throughout the three exposures: the first vial replicate to enter the consumption assay was the first to be allowed to recover for each of the three exposures; and the last vial replicate to enter the consumption assay was the last to be allowed to recover for each of the three exposures. In this way flies in all vials were exposed to the same duration of feeding and recovery. A total of 10 replicate vials were tested for each genotype and sex.

We defined four behaviors: total amount of each solution consumed, drug preference, and change in consumption and change of preference between exposures 3 and 1. Preference was quantified in two ways: as the difference between the amount of drug and sucrose consumed (Preference A), and as this difference scaled by the total amount consumed (Preference B).

Genetic variation in drug consumption behaviors in the DGRP

We performed four-way factorial mixed model analyses of variance (ANOVA) to partition variation in consumption in the DGRP: Y = μ + L + E + S + X + (L × E) + (L × S) + (L × X) + (E × S) + (E × X) + (S × X) + (L × E × S) + (L × E × X) + (L × S × X) + (E × S × X) + (L × E × S × X) + ε, where Y is consumption; μ is the overall mean; L is the random effect of line; E, S, and X are the fixed effects of exposure (day 1–3), solution (drug, sucrose), and sex (males, females); and ε is the residual variation between replicate vials. The main effect of L and all interaction terms with L are genetic factors affecting drug consumption. We also ran the same ANOVA models to compare the effects of cocaine and methamphetamine on consumption, separately for males and females. The full model for variation in change in consumption over time is Y = μ + L + S + X + (L × S) + (L × X) + (S × X) + (L × S × X) + ε. We assessed variation in the development of preference using the model Y = μ + L + E + X + (L × E) + (L × X) + (E × X) + (L × E × X) + ε. We also assessed whether there is natural variation in the change of preference over time using the model Y = μ + L + X + (L × X) + ε. We also ran reduced models for each trait. All ANOVAs were performed using the PROC GLM function in SAS. We used the R function pf to assign exact P-values.

Quantitative genetic analyses in the DGRP

We used the SAS PROC MIXED function to estimate variance components for each of the random effect terms in the full and reduced models. The R package lmer and lmerTest were utilized in combination with the pchisq function to assign P-values for the segregating genetic variation for each trait. We computed broad sense heritabilities as the sum of all genetic variance components divided by the total phenotypic variance for each model, and broad sense heritabilities of line means as the sum of all genetic variance components divided by the sum of all genetic variance components plus the environmental variance/10, where 10 is the number of replicate vials per line, sex, exposure and treatment. We computed pairwise genetic correlations as , where is the among line variance from the appropriate two-way factorial ANOVA and σL1 and σL2 are the among line standard deviations from the one-way ANOVA for each condition. We computed Pearson product-moment correlations of line means to estimate phenotypic correlations between different traits.

Genome wide association mapping in the DGRP

We performed GWA analyses on line means for all consumption traits using the DGRP pipeline (http://dgrp2.gnets.ncsu.edu/) [43]. This pipeline implements single-variant tests of association for additive effects of variants with minor allele frequencies ≥ 0.05. We tested effects of 1,891,456 DNA sequence variants on each trait.

Network analysis

We annotated candidate genes identified by the GWA analyses using Flybase release 5.57 [56] and mapped gene-gene networks through the genetic interaction database downloaded from Flybase. We then constructed a subnetwork using Cytoscape 3.5.1 where candidate genes directly interact with each other. We evaluated the significance (α = 0.05) of the constructed subnetwork by a randomization test [9193].

Gene ontology analysis

We carried out gene ontology (GO) enrichment analysis with PANTHER 11.1 (http://pantherdb.org/) [53,54].

RNAi knockdown of gene expression

We used the binary GAL4-UAS system for RNAi-targeted knockdown of expression of candidate genes associated with variation in consumption of cocaine or methamphetamine with a weak ubiquitous driver (Ubi156-GAL4) and drivers specific for neurons (elav-GAL4), glia (repo-GAL4), mushroom bodies (201Y-GAL4) and dopaminergic neurons (TH-GAL4). We crossed 3 homozygous GAL4 driver males to 5–7 homozygous females harboring a unique UAS-RNAi transgene or the progenitor control to generate F1 GAL4-UAS-RNAi and GAL4 control progeny. We assessed the consumption traits exactly as described above for the DGRP lines. Differences between RNAi lines and their corresponding control lines for consumption were assessed with a fixed-effect ANOVA, separately for males and females. The full model was: Y = μ + L + E + S + (L × E) + (E × S) + (L × S) + (L × E × S) + ɛ, where Y denotes the mean consumption, E denotes the different exposures, L is the line (Control or RNAi), S denotes the different solutions (sucrose or cocaine/methamphetamine), and ɛ the error variance. Differences between RNAi lines and controls for change in consumption and preference were also assessed with fixed-effect ANOVAs. The full model for change in consumption was: Y = μ + L + S + (L × S) + ɛ, while the full model for preference was Y = μ + L + E + (L × E) + ɛ. All ANOVAs were run using R.

Supporting information

S1 Table. DGRP raw consumption data.

(A) Cocaine experiment. (B) Methamphetamine experiment. F: female; M: male.

https://doi.org/10.1371/journal.pgen.1007834.s001

(XLSX)

S2 Table. Analyses of variance of consumption, change in consumption, preference and change in preference of cocaine and methamphetamine.

Exposure, Sex, Solution, and their interaction are fixed effects, the rest are random. Mixed model three-way factorial ANOVAs are given for males and females, as well as reduced models by Exposure, Sex, and Solution. E: Exposure; X: Sex; S: Solution; L: DGRP Line; df: degrees of freedom; MS: Type III mean squares; F: F-ratio test; P: P-value; σ2: variance component estimate; SE: standard error; H2: Broad sense heritability. Significant P-values are shown in red font. (A) Cocaine experiment. (B) Methamphetamine experiment.

https://doi.org/10.1371/journal.pgen.1007834.s002

(XLSX)

S3 Table. DGRP line means for all traits.

(A) Cocaine experiment. (B) Methamphetamine experiment. Means are given in mm; 1 mm = 0.067 μl.

https://doi.org/10.1371/journal.pgen.1007834.s003

(XLSX)

S4 Table. Genetic and phenotypic correlations between traits.

(A) Cross-sex, cross-exposure and cross-solution genetic correlations. Significant P-values are indicated in red font. (B) Pair-wise phenotypic correlations. Entries in the cells are the correlation coefficients and the cell color denotes the P-value. Red: P < 0.0001; orange: P < 0.001; yellow: P < 0.01; green: P < 0.05; white: P > 0.05.

https://doi.org/10.1371/journal.pgen.1007834.s004

(XLSX)

S5 Table. Results of genome wide association (GWA) analyses for consumption behaviors.

(A) Top variants (P < 5 ×10−5) and associated genes for each trait. (B) Variants and genes for the cocaine traits, the methamphetamine traits, and variants and genes overlapping between the two experiments. (C) Pathway and gene ontology enrichment analysis for the cocaine GWA analyses. (D) Pathway and gene ontology enrichment analysis for the methamphetamine GWA analyses.

https://doi.org/10.1371/journal.pgen.1007834.s005

(XLSX)

S6 Table. DGRP candidate genes and human orthologs.

The references indicate which of the human orthologs have been associated with addictive phenotypes.

https://doi.org/10.1371/journal.pgen.1007834.s006

(XLSX)

S7 Table. A significant genetic interaction network with no missing genes.

(A) Genes in network. (B) Human orthologs. (C) Pathway and gene ontology enrichment analysis.

https://doi.org/10.1371/journal.pgen.1007834.s007

(XLSX)

S8 Table. Raw cocaine and sucrose consumption data for RNAi and control genotypes.

(A) Ubi156-GAL4. (B) elav-GAL4. (C) repo-GAL4. (D) 201Y-GAL4. (E) TH-GAL4. Data are given in mm; 1 mm = 0.067 μl.

https://doi.org/10.1371/journal.pgen.1007834.s008

(XLSX)

S9 Table. Raw methamphetamine and sucrose consumption data for RNAi and control genotypes.

(A) Ubi156-GAL4. (B) elav-GAL4. (C) repo-GAL4. (D) 201Y-GAL4. (E) TH-GAL4. Data are given in mm; 1 mm = 0.067 μl.

https://doi.org/10.1371/journal.pgen.1007834.s009

(XLSX)

S10 Table. Analyses of variance of consumption, change in consumption, preference and change in preference of cocaine and sucrose in RNAi lines and their controls.

Fixed effect three-way factorial ANOVAs are given for males and females as well as reduced models by Exposure and Solution. E: Exposure; S: Solution; L: RNAi or control genotype; df: degrees of freedom; MS: Type III mean squares; F: F-ratio test; P: P-value. Significant P-values are shown in red font. (A) Ubi156-GAL4. (B) elav-GAL4. (C) repo-GAL4. (D) 201Y-GAL4. (E) TH-GAL4.

https://doi.org/10.1371/journal.pgen.1007834.s010

(XLSX)

S11 Table. Analyses of variance of consumption, change in consumption, preference and change in preference of methamphetamine and sucrose in RNAi lines and their controls.

Fixed effect three-way factorial ANOVAs are given for males and females as well as reduced models by Exposure and Solution. E: Exposure; S: Solution; L: RNAi or control genotype; df: degrees of freedom; MS: Type III mean squares; F: F-ratio test; P: P-value. Significant P-values are shown in red font. (A) Ubi156-GAL4. (B) elav-GAL4. (C) repo-GAL4. (D) 201Y-GAL4. (E) TH-GAL4.

https://doi.org/10.1371/journal.pgen.1007834.s011

(XLSX)

S12 Table. Drosophila lines used in this study.

(A) DGRP lines. (B) RNAi lines and control genotypes. (C) GAL4 driver lines.

https://doi.org/10.1371/journal.pgen.1007834.s012

(XLSX)

S1 Fig. P-value summary from three-way ANOVA models of consumption for UAS-RNAi and control genotypes of candidate genes crossed to a weak ubiquitous GAL4 driver (Ubi156-GAL4).

Red: P < 0.0001; orange: P < 0.001; yellow: P < 0.01; green: P < 0.05; white: P > 0.05.

https://doi.org/10.1371/journal.pgen.1007834.s013

(PDF)

S2 Fig. Differences between Ubi156-GAL4 RNAi and control genotypes for 34 candidate genes.

(A) Cocaine preference, females. (B) Cocaine preference, males. (C) Change in cocaine preference between third and first exposures, females. (D) Change in cocaine preference between third and first exposures, males. Asterisks represent significant L×S terms (A, B) or significant L×S×E terms from the full ANOVA models. Exact P-values are given in S10 Table.

https://doi.org/10.1371/journal.pgen.1007834.s014

(PDF)

S3 Fig. Differences between Ubi156-GAL4 RNAi and control genotypes for 34 candidate genes.

(A) Methamphetamine preference, females. (B) Methamphetamine preference, males. (C) Change in methamphetamine preference between third and first exposures, females. (D) Change in methamphetamine preference between third and first exposures, males. Asterisks represent significant L×S terms (A, B) or significant L×S×E terms from the full ANOVA models. Exact P-values are given in S11 Table.

https://doi.org/10.1371/journal.pgen.1007834.s015

(PDF)

S4 Fig. P-value summary from three-way ANOVA models of consumption for UAS-RNAi and control genotypes of candidate genes crossed to neuronal (elav-GAL4) and glial (repo-GAL4) GAL4 drivers.

Red: P < 0.0001; orange: P < 0.001; yellow: P < 0.01; green: P < 0.05; white: P > 0.05.

https://doi.org/10.1371/journal.pgen.1007834.s016

(PDF)

S5 Fig. P-value summary from three-way ANOVA models of consumption for UAS-RNAi and control genotypes of candidate genes crossed to mushroom body (201Y-GAL4) and dopaminergic (TH-GAL4) GAL4 drivers.

Red: P < 0.0001; orange: P < 0.001; yellow: P < 0.01; green: P < 0.05; white: P > 0.05.

https://doi.org/10.1371/journal.pgen.1007834.s017

(PDF)

Acknowledgments

We thank Mathew J. Eddinger, Emily L. Davis, Kishan N. Patel, and Shaunaci Stevens for technical assistance.

References

  1. 1. Galbraith N. The methamphetamine problem. BJPsych Bull. 2015;39: 218–220. pmid:26755964
  2. 2. World Drug Report 2018 (United Nations publication, Sales No. E.18.XI.9)
  3. 3. National Institute on Drug Abuse. Trends & Statistics. 2017; https://www.drugabuse.gov/related-topics/trends-statistics
  4. 4. Horne MK, Lee J, Chen F, Lanning K, Tomas D, Lawrence AJ. Long-term administration of cocaine or serotonin reuptake inhibitors results in anatomical and neurochemical changes in noradrenergic, dopaminergic, and serotonin pathways. J Neurochem. 2008;106: 1731–1744. pmid:18624911
  5. 5. Verma V. Classic studies on the interaction of cocaine and the dopamine transporter. Clin Psychopharmacol Neurosci. 2015;13: 227–238. pmid:26598579
  6. 6. Seiden LS, Sabol KE, Ricaurte GA. Amphetamine: Effects on catecholamine systems and behavior. Annu Rev Pharmacol Toxicol. 1993;33: 639–676. pmid:8494354
  7. 7. Sulzer D, Sonders MS, Poulsen NW, Galli A. Mechanisms of neurotransmitter release by amphetamines: a review. Prog Neurobiol. 2005;75: 406–433. pmid:15955613
  8. 8. Mariani JJ, Levin FR. Treatment strategies for co-occurring ADHD and substance use disorders. Am J Addict. 16: 45–56. pmid:17453606
  9. 9. Bonci A, Bernardi G, Grillner P, Mercuri NB. The dopamine-containing neuron: maestro or simple musician in the orchestra of addiction? Trends Pharmacol Sci. 2003; 24: 172–177. pmid:12707003
  10. 10. Roncero C, Daigre C, Grau-López L, Barral C, Pérez-Pazos J, Martínez-Luna N, Casas M. An international perspective and review of cocaine-induced psychosis: a call to action. Subst Abus. 2014; 35: 321–327. pmid:24927026
  11. 11. Hellem TL. A review of methamphetamine dependence and withdrawal treatment: A focus on anxiety outcomes. J Subst Abuse Treat. 2016;71: 16–22. pmid:27776672
  12. 12. Hall FS, Drgonova J, Jain S, Uhl GR. Implications of genome wide association studies for addiction: Are our a priori assumptions all wrong? Pharmacol Ther. 2013;140: 267–279. pmid:23872493
  13. 13. Tsuang MT, Bar JL, Harley RM, Lyons MJ. The Harvard twin study of substance abuse: what we have learned. Harv Rev Psychiatry. 2001;9: 267–279. pmid:11600486
  14. 14. Agrawal A, Lynskey MT. Are there genetic influences on addiction: evidence from family, adoption and twin studies. Addiction. 2008;103: 1069–1081. pmid:18494843
  15. 15. Noble EP, Blum K, Khalsa ME, Ritchie T, Montgomery A, Wood RC, et al. Allelic association of the D2 dopamine receptor gene with cocaine dependence. Drug Alcohol Depend. 1993;33: 271–285. pmid:8261891
  16. 16. Comings DE, Muhleman D, Ahn C, Gysin R, Flanagan SD. The dopamine D2 receptor gene: a genetic risk factor in substance abuse. Drug Alcohol Depend. 1994; 34: 175–180. pmid:8033754
  17. 17. Persico AM, Bird G, Gabbay FH, Uhl GR. D2 dopamine receptor gene TaqI A1 and B1 restriction fragment length polymorphisms: enhanced frequencies in psychostimulant-preferring polysubstance abusers”. Biol Psychiatry. 1996;40: 776–784. pmid:8894071
  18. 18. Moyer RA, Wang D, Papp AC, Smith RM, Duque L, Mash DC, et al. Intronic polymorphisms affecting alternative splicing of human dopamine D2 receptor are associated with cocaine abuse. Neuropsychopharmacol. 2011;36: 753–762.
  19. 19. Gelernter J, Southwick S, Goodson S, Morgan A, Nagy L, Charney DS. No association between D2 dopamine receptor (DRD2) “A” system alleles, or DRD2 haplotypes, and posttraumatic stress disorder. Biol Psychiatry. 1999;45: 620–625. pmid:10088049
  20. 20. Serý O, Vojtová V, Zvolský P. The association study of DRD2, ACE and AGT gene polymorphisms and methamphetamine dependence. Physiol Res. 2001;50: 43–50. pmid:11300226
  21. 21. Hong C-J, Cheng C-Y, Shu L-RR, Yang C-Y, Tsai S-J. Association study of the dopamine and serotonin transporter genetic polymorphisms and methamphetamine abuse in Chinese males. J Neural Transm. 2003;110: 345–351. pmid:12658362
  22. 22. Messas G, Meira-Lima I, Turchi M, Franco O, Guindalini C, Castelo A, et al. Association study of dopamine D2 and D3 receptor gene polymorphisms with cocaine dependence. Psychiatr Genet. 2005;15: 171–174. pmid:16094250
  23. 23. Fernàndez-Castillo N, Ribasés M, Roncero C, Casas M, Gonzalvo B, Cormand B. Association study between the DAT1, DBH and DRD2 genes and cocaine dependence in a Spanish sample. Psychiatr Genet. 2010;20: 317–320. pmid:20505554
  24. 24. Lohoff FW, Bloch PJ, Hodge R, Nall AH, Ferraro TN, Kampman KM, et al. Association analysis between polymorphisms in the dopamine D2 receptor (DRD2) and dopamine transporter (DAT1) genes with cocaine dependence. Neurosci Lett. 2010;473: 87–91. pmid:20170711
  25. 25. Gelernter J, Kranzler HR, Satel SL, Rao PA. Genetic association between dopamine transporter protein alleles and cocaine-induced paranoia. Neuropsychopharmacol. 1994;11: 195–200.
  26. 26. Liu HC, Lin SK, Liu SK, Chen SL, Hu CJ, Chang JG, et al. DAT polymorphism and diverse clinical manifestations in methamphetamine abusers. Psychiatr Genet. 2004; 14:33–37. pmid:15091313
  27. 27. Guindalini C, Howard M, Haddley K, Laranjeira R, Collier D, Ammar N, et al. A dopamine transporter gene functional variant associated with cocaine abuse in a Brazilian sample. Proc Natl Acad Sci USA. 2006;103: 4552–4557. pmid:16537431
  28. 28. Martinez D, Gelernter J, Abi-Dargham A, van Dyck CH, Kegeles L, Innis RB, et al. The variable number of tandem repeats polymorphism of the dopamine transporter gene is not associated with significant change in dopamine transporter phenotype in humans. Neuropsychopharmacol. 2001;24: 553–560.
  29. 29. Patkar AA, Berrettini WH, Hoehe M, Hill KP, Sterling RC, Gottheil E, et al. Serotonin transporter (5-HTT) gene polymorphisms and susceptibility to cocaine dependence among African-American individuals. Addict Biol. 2001;6: 337–345. pmid:11900612
  30. 30. Pandey UB, Nichols CD. Human disease models in Drosophila melanogaster and the role of the fly in therapeutic drug discovery. Pharmacol Rev. 2011;63: 411–436. pmid:21415126
  31. 31. Wang KH, Penmatsa A, Gouaux E. Neurotransmitter and psychostimulant recognition by the dopamine transporter. Nature. 2015;521: 322–327. pmid:25970245
  32. 32. Andretic R, van Swinderen B, Greenspan RJ. Dopaminergic modulation of arousal in Drosophila. Curr Biol. 2005;15: 1165–1175. pmid:16005288
  33. 33. Andretic R, Kim Y-C, Jones FS, Han K-A, Greenspan RJ. Drosophila D1 dopamine receptor mediates caffeine-induced arousal. Proc Natl Acad Sci USA. 2008;105: 20392–20397. pmid:19074291
  34. 34. Walters KR, Rupassara SI, Markelz RJC, Leakey ADB, Muir WM, Pittendrigh BR. Methamphetamine causes anorexia in Drosophila melanogaster, exhausting metabolic reserves and contributing to mortality. J Toxicol Sci. 2012;37: 773–790. pmid:22863857
  35. 35. McClung C, Hirsh J. Stereotypic behavioral responses to free-base cocaine and the development of behavioral sensitization in Drosophila. Curr Biol. 1998;8: 109–112. pmid:9427649
  36. 36. Torres G, Horowitz JM. Activating properties of cocaine and cocaethylene in a behavioral preparation of Drosophila melanogaster. Synapse. 1998;29: 148–161. pmid:9593105
  37. 37. Bainton RJ, Tsai LT, Singh CM, Moore MS, Neckameyer WS, Heberlein U. Dopamine modulates acute responses to cocaine, nicotine and ethanol in Drosophila. Curr Biol. 2000;10: 187–194. pmid:10704411
  38. 38. Rothenfluh A, Heberlein U. Drugs, flies, and videotape: the effects of ethanol and cocaine on Drosophila locomotion. Curr Opin Neurobiol. 2002;12: 639–645. pmid:12490253
  39. 39. Wolf FW, Heberlein U. Invertebrate models of drug abuse. Journal of Neurobiology. 54: 161–178. pmid:12486703
  40. 40. Dimitrijevic N, Dzitoyeva S, Manev H. An automated assay of the behavioral effects of cocaine injections in adult Drosophila. J Neurosci Methods. 2004;137: 181–184. pmid:15262059
  41. 41. George R, Lease K, Burnette J, Hirsh J. A “bottom-counting” video system for measuring cocaine-induced behaviors in Drosophila. In: Young MW, editor. Methods in Enzymology. Academic Press; 2005. pp. 841–851. http://www.sciencedirect.com/science/article/pii/S0076687905930449 https://doi.org/10.1016/S0076-6879(05)93044-9 pmid:15817327
  42. 42. Mackay TFC, Richards S, Stone EA, Barbadilla A, Ayroles JF, Zhu D, et al. The Drosophila melanogaster Genetic Reference Panel. Nature. 2012;482: 173–178. pmid:22318601
  43. 43. Huang W, Massouras A, Inoue Y, Peiffer J, Ràmia M, Tarone AM, et al. Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines. Genome Res. 2014;24: 1193–1208. pmid:24714809
  44. 44. Ja WW, Carvalho GB, Mak EM, Rosa NN de la, Fang AY, Liong JC, et al. Prandiology of Drosophila and the CAFE assay. Proc Natl Acad Sci USA. 2007;104: 8253–8256. pmid:17494737
  45. 45. Devineni AV, Heberlein U. Preferential ethanol consumption in Drosophila models features of addiction. Curr Biol. 2009;19(24):2126–2132. pmid:20005106
  46. 46. Devineni AV, Heberlein U. Addiction-like behavior in Drosophila. Commun Integr Biol. 2010;3(4):357–359. pmid:20798826
  47. 47. Dietzl G, Chen D, Schnorrer F, Su K-C, Barinova Y, Fellner M, et al. A genome-wide transgenic RNAi library for conditional gene inactivation in Drosophila. Nature. 2007;448: 151–156. pmid:17625558
  48. 48. Holgate JY, Garcia H, Chatterjee S, Bartlett SE. Social and environmental enrichment has different effects on ethanol and sucrose consumption in mice. Brain Behav. 2017;7. pmid:28828224
  49. 49. Martinelli DC, Chew KS, Rohlmann A, Lum MY, Ressl S, Hattar S, et al. Expression of C1ql3 in discrete neuronal populations controls efferent synapse numbers and diverse behaviors. Neuron. 2016;91: 1034–1051. pmid:27478018
  50. 50. Liu QR, Drgon T, Johnson C, Walther D, Hess J, Uhl GR. Addiction molecular genetics: 639,401 SNP whole genome association identifies many “cell adhesion” genes. Am J Med Genet B Neuropsychiatr Genet. 2006;141B: 918–925. pmid:17099884
  51. 51. Tarlungeanu DC, Deliu E, Dotter CP, Kara M, Janiesch PC, Scalise M, et al. Impaired amino acid transport at the blood brain barrier is a cause of autism spectrum disorder. Cell. 2016;167: 1481–1494.e18. pmid:27912058
  52. 52. Kaun KR, Devineni AV, Heberlein U. Drosophila melanogaster as a model to study drug addiction. Hum Genet. 2012;131: 959–975. pmid:22350798
  53. 53. Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13: 2129–2141. pmid:12952881
  54. 54. Mi H, Muruganujan A, Casagrande JT, Thomas PD. Large-scale gene function analysis with the PANTHER classification system. Nat Protoc. 2013;8: 1551–1566. pmid:23868073
  55. 55. McQuilton P, St Pierre SE, Thurmond J. FlyBase 101 –the basics of navigating FlyBase. Nucleic Acids Res. 2012;40: D706–D714. pmid:22127867
  56. 56. Davis RL. Traces of Drosophila memory. Neuron. 2011;70: 8–19. pmid:21482352
  57. 57. Widmer YF, Bilican A, Bruggmann R, Sprecher SG. Regulators of long-term memory revealed by mushroom body-specific gene expression profiling in Drosophila melanogaster. Genetics. 2018;209:1167–1181. pmid:29925565
  58. 58. Salahpour A, Ramsey AJ, Medvedev IO, Kile B, Sotnikova TD, Holmstrand E, et al. Increased amphetamine-induced hyperactivity and reward in mice overexpressing the dopamine transporter. Proc Natl Acad Sci USA. 2008;105: 4405–4410. pmid:18347339
  59. 59. Hoops D, Reynolds LM, Restrepo-Lozano J-M, Flores C. Dopamine development in the mouse orbital prefrontal cortex Is protracted and sensitive to amphetamine in adolescence. eNeuro. 2018;5. pmid:29333488
  60. 60. Walker QD, Kuhn CM. Cocaine increases stimulated dopamine release more in periadolescent than adult rats. Neurotoxicol Teratol. 2008;30: 412–418. pmid:18508233
  61. 61. Koulchitsky S, De Backer B, Quertemont E, Charlier C, Seutin V. Differential effects of cocaine on dopamine neuron firing in awake and anesthetized rats. Neuropsychopharmacol. 2012;37: 1559–1571. pmid:22298123
  62. 62. Jedema HP, Narendran R, Bradberry CW. Amphetamine-induced release of dopamine in primate prefrontal cortex and striatum: striking differences in magnitude and time course. J Neurochem. 2014;130: 490–497. pmid:24749782
  63. 63. Ashok AH, Mizuno Y, Volkow ND, Howes OD. Association of stimulant use with dopaminergic alterations in users of cocaine, amphetamine, or methamphetamine: A systematic review and meta-analysis. JAMA Psychiatry. 2017;74: 511–519. pmid:28297025
  64. 64. Becker JB, Hu M. Sex differences in drug abuse. Front Neuroendocrinol. 2008;29: 36–47. pmid:17904621
  65. 65. Fattore L, Altea S, Fratta W. Sex differences in drug addiction: a review of animal and human studies. Womens Health (Lond). 2008;4:51–65. pmid:19072451
  66. 66. Fattore L, Fadda P, Fratta W. Sex differences in the self-administration of cannabinoids and other drugs of abuse. Psychoneuroendocrinology. 2009;34 Suppl 1:S227–S236. pmid:19744795
  67. 67. Cox BM, Young AB, See RE, Reichel CM. Sex differences in methamphetamine seeking in rats: Impact of oxytocin. Psychoneuroendocrinol. 2013;38: 2343–2353. pmid:23764194
  68. 68. Fox HC, Morgan PT, Sinha R. Sex differences in guanfacine effects on drug craving and stress arousal in cocaine-dependent individuals. Neuropsychopharmacol. 2014;39: 1527–1537. pmid:24395021
  69. 69. Kennedy AP, Epstein DH, Phillips KA, Preston KL. Sex differences in cocaine/heroin users: drug-use triggers and craving in daily life. Drug Alcohol Depend. 2013;132: 29–37. pmid:23357742
  70. 70. Becker JB, McClellan ML, Reed BG. Sex differences, gender and addiction. J Neurosci Res. 2017;95: 136–147. pmid:27870394
  71. 71. Calipari ES, Juarez B, Morel C, Walker DM, Cahill ME, Ribeiro E, et al. Dopaminergic dynamics underlying sex-specific cocaine reward. Nat Commun. 2017;8: 13877. pmid:28072417
  72. 72. Kogachi S, Chang L, Alicata D, Cunningham E, Ernst T. Sex differences in impulsivity and brain morphometry in methamphetamine users. Brain Struct Funct. 2017;222: 215–227. pmid:27095357
  73. 73. McClung C, Hirsh J. The trace amine tyramine is essential for sensitization to cocaine in Drosophila. Curr Biol. 1999;9: 853–860. pmid:10469593
  74. 74. Andretic R, Chaney S, Hirsh J. Requirement of circadian genes for cocaine sensitization in Drosophila. Science. 1999;285: 1066–1068. pmid:10446052
  75. 75. Bainton RJ, Tsai LT, Schwabe T, DeSalvo M, Gaul U, Heberlein U. moody encodes two GPCRs that regulate cocaine behaviors and blood-brain barrier permeability in Drosophila. Cell. 2005;123:145–156. pmid:16213219
  76. 76. Heberlein U, Tsai LT, Kapfhamer D, Lasek AW. Drosophila, a genetic model system to study cocaine-related behaviors: a review with focus on LIM-only proteins. Neuropharmacol. 2009;56 Suppl 1:97–106. pmid:18694769
  77. 77. Vanderwerf SM, Buck DC, Wilmarth PA, Sears LM, David LL, Morton DB, Neve KA. Role for Rab10 in methamphetamine-induced behavior. PLoS One. 2015;10: e0136167. pmid:26291453
  78. 78. Takaesu NT, Hyman-Walsh C, Ye Y, Wisotzkey RG, Stinchfield MJ, O'Connor MB, et al. dSno facilitates baboon signaling in the Drosophila brain by switching the affinity of Medea away from Mad and toward dSmad2. Genetics. 2006;174: 1299–1313. pmid:16951053
  79. 79. Barrio R, López-Varea A, Casado M, de Celis JF. Characterization of dSnoN and its relationship to Decapentaplegic signaling in Drosophila. Dev Biol. 2007;306: 66–81. pmid:17434471
  80. 80. Iyer EP, Iyer SC, Sullivan L, Wang D, Meduri R, Graybeal LL, et al. Functional genomic analyses of two morphologically distinct classes of Drosophila sensory neurons: post-mitotic roles of transcription factors in dendritic patterning. PLoS ONE. 2013;8: e72434. pmid:23977298
  81. 81. Bai J, Chiu W, Wang J, Tzeng T, Perrimon N, Hsu J. The cell adhesion molecule Echinoid defines a new pathway that antagonizes the Drosophila EGF receptor signaling pathway. Development. 2001;128: 591–601. pmid:11171342
  82. 82. Rawlins EL, Lovegrove B, Jarman AP. Echinoid facilitates Notch pathway signalling during Drosophila neurogenesis through functional interaction with Delta. Development. 2003;130: 6475–6484. pmid:14627723
  83. 83. Wei SY, Escudero LM, Yu F, Chang LH, Chen LY, Ho YH, et al. Echinoid is a component of adherens junctions that cooperates with DE-Cadherin to mediate cell adhesion. Dev Cell. 2005;8: 493–504. pmid:15809032
  84. 84. Nakamura M, Okano H, Blendy JA, Montell C. Musashi, a neural RNA-binding protein required for Drosophila adult external sensory organ development. Neuron. 1994;13: 67–81. pmid:8043282
  85. 85. Han KA, Millar NS, Davis RL. A novel octopamine receptor with preferential expression in Drosophila mushroom bodies. J Neurosci. 1998;18: 3650–3658. https://doi.org/10.1523/JNEUROSCI.18-10-03650.1998 pmid:9570796
  86. 86. Kim YC, Lee HG, Han KA. D1 dopamine receptor dDA1 is required in the mushroom body neurons for aversive and appetitive learning in Drosophila. J Neurosci. 2007;27: 7640–7647. pmid:17634358
  87. 87. Lim J, Fernandez AI, Hinojos SJ, Aranda GP, James J, Seong CS, et al. The mushroom body D1 dopamine receptor controls innate courtship drive. Genes Brain Behav. 2018;17: 158–167. pmid:28902472
  88. 88. Lebestky T, Chang JS, Dankert H, Zelnik L, Kim YC, Han KA, et al. Two different forms of arousal in Drosophila are oppositely regulated by the dopamine D1 receptor ortholog DopR via distinct neural circuits. Neuron 2009;64: 522–536. pmid:19945394
  89. 89. Ma Z, Stork T, Bergles DE, Freeman MR. Neuromodulators signal through astrocytes to alter neural circuit activity and behaviour. Nature. 2016;539: 428–432. pmid:27828941
  90. 90. Garlapow ME, Huang W, Yarboro MT, Peterson KR, Mackay TFC. Quantitative genetics of food intake in Drosophila melanogaster. PLoS One. 2015;10: e0138129. pmid:26375667
  91. 91. Antonov AV, Dietmann S, Mewes HW. KEGG spider: interpretation of genomics data in the context of the global gene metabolic network. Genome Biol. 2008;9: R179. pmid:19094223
  92. 92. Zhou S, Morozova TV, Hussain YN, Luoma SE, McCoy L, Yamamoto A, et al. The genetic basis for variation in sensitivity to lead toxicity in Drosophila melanogaster. Environ Health Perspect. 2016;124: 1062–1070. pmid:26859824
  93. 93. Morozova TV, Hussain Y, McCoy LJ, Zhirnov EV, Davis MR, Pray VA, et al. A Cyclin E centered genetic network contributes to alcohol-induced variation in Drosophila development. G3 (Bethesda). 2018;8: 2643–2653. pmid:29871898