Abstract
Individual differences in most behavioral and pharmacological responses to abused drugs are dependent on both genetic and environmental factors. The genetic influences on the complex phenotypes related to drug abuse have been difficult to study using classical genetic analyses. Quantitative trait locus (QTL) mapping is a method that has been used successfully to examine genetic contributions to some of these traits by correlating allelic variation in polymorphic genetic markers of known chromosomal location with variation in drug-response phenotypes. We evaluated several behavioral responses to multiple doses of methamphetamine (METH) in C57BL/6J (B6), DBA/2J (D2), and 25 of their recombinant inbred (BXD RI) strains. Stereotyped chewing, horizontal home cage activity, and changes in body temperature after 0, 4, 8, or 16 mg/kg METH, as well as stereotyped climbing behavior after 16 mg/kg METH, were examined. Associations (p < 0.01) between METH sensitivity and allelic status at multiple microsatellite genetic markers were subsequently determined for each response. QTLs were provisionally identified for each phenotype, some unique to a particular behavior and others that appeared to influence multiple phenotypes. Candidate genes suggested by these analyses included several that mapped near genes relevant for the neurotransmitters acetylcholine and glutamate. The locations of QTLs provisionally identified by this analysis were compared with QTLs hypothesized in other studies to influence methamphetamine- and cocaine-related phenotypes. In several instances, QTLs appeared to overlap, which is consistent with idea that common neural substrates underlie some responses to psychostimulants.
- methamphetamine
- QTL
- quantitative trait locus
- recombinant inbred strain
- mouse
- gene mapping
- stereotypy
- locomotor activity
- body temperature pharmacogenetics
Most responses to drugs of abuse reflect both genetic and environmental influences (Crabbe and Harris, 1991). Furthermore, the genetic influence on drug responses generally results from the collective influences of many genes, each of which may contribute a relatively small amount to the genotypic variance. These quantitative traits have been effectively studied in recombinant inbred (RI) strains of mice. The BXD RI strains [derived by inbreeding novel inbred strains from the F2 cross between C57BL/6J (B6) and DBA/2J (D2) mice] have been the most commonly used in drug abuse research. Under controlled laboratory conditions, phenotypic differences among RI strains largely reflect their genotypic differences: the pattern of B6 and D2 alleles captured in each RI strain is unique because of distinct recombinations of these alleles. Phenotypic strain means can be correlated with allelic status at known genetic markers: a significant relationship indicates a putative quantitative trait locus (QTL) in the same chromosome region as the marker. Each such QTL contains one or more genes that may influence the trait. This approach enables the provisional identification of candidate genes that may influence a particular behavior or phenotype in the absence of any prior hypotheses about the mechanisms by which such phenotypes are expressed. Furthermore, it is not necessary for the progenitor strains to differ phenotypically; as long as the trait is polygenic (influenced by multiple genes) and variable within the RI panel, QTL analysis can be performed successfully. This might be the case, for instance, if both progenitor strains possess genes for high and low response. Another major strength of this technique is that because inbred strains are genetically stable, all information obtained in these studies is both cumulative and comparable to earlier data, thus providing a rich resource and incentive for collaborative efforts, because genetic correlations can be determined for different measures across different laboratories. Perhaps most important, genes identified in these murine analyses can usually be readily mapped to particular human chromosomal regions because of the high degree of synteny (Copeland et al., 1993). A description of the strategy for QTL mapping has been published previously (Crabbe and Belknap, 1993; Belknap et al., 1996a,b).
Variability in the response to amphetamine in mice is known to depend in part on genetic differences (Kitahama and Valatx, 1979; Crabbe, 1986; DeWit et al., 1986) (for review, see Seale, 1991). Seale et al. (1985) were the first to attempt a systematic genetic analysis of amphetamine sensitivity. Body temperature was measured before and several times after d-amphetamine administration in B6, D2, and 10 of their RI derivatives. Later, Gora-Maslak et al. (1991)subjected the Seale et al. (1985) data to QTL analysis and found several promising associations. Although the power of this early experiment to map the responsible genes was limited because only 10 of the RI strains were studied and there were only 173 polymorphic loci available at the time, it demonstrated the feasibility of locating drug-response genes. Several studies have since used the QTL mapping strategy to locate genes influencing other abused drugs, including those affecting morphine preference (Berettini et al., 1994), alcohol preference (Phillips et al., 1994; Rodriguez et al., 1995; Melo et al., 1996), and alcohol withdrawal (K. Buck, P. Metten, J. Belknap, and J. Crabbe, unpublished observations).
The purpose of the present experiments was to investigate further the genetic substrates of amphetamine-induced responses in mice. Behavioral and physiological responses (stereotyped chewing and climbing, alterations in thermoregulation, and home cage locomotor activity) after administration of 0, 4, 8, or 16 mg/kg methamphetamine (METH) were evaluated using B6, D2, and 25 of their RI strains. Genetic correlations were calculated between these measures. QTL analysis was performed on these strain means, and QTLs suggested by these experiments were then compared with QTLs hypothesized to influence psychostimulant-related phenotypes from other studies, including other measures of locomotor activation and stereotypy to amphetamine and cocaine.
MATERIALS AND METHODS
Animals and drugs. Mice were bred at the Portland VA Veterinary Medical Unit from stock originally obtained from The Jackson Laboratory (Bar Harbor, ME). They were group housed (2–5 animals per cage) at the time of weaning and maintained with ad libitumfood and water in a 12 hr light/dark cycle (6:00 A.M. to 6:00 P.M.) at 21 ± 2°C. Adult (10–14 week old) mice were used in all studies, and testing was conducted between 8:00 A.M. and 12:00 P.M. (see below). Animals in each experiment were tested over several weeks, with strains roughly counterbalanced across the experimental period. Each mouse was tested once, but each strain was tested at multiple times. On test days, cages of mice were transported to the testing room and mice were weighed 45–60 min before the start of the test session and then randomly assigned to either the METH or the saline (SAL) group, with multiple doses represented in each cage. Animals were given either SAL or METH (4, 8, or 16 mg/kg, i.p.; 0.8, 1.2, or 1.6 mg/ml in a volume of 0.01 ml/gm body weight; Sigma, St. Louis MO) and tested as described below.
The development of the BXD RI strains has been described previously (Taylor, 1978). Briefly, inbreeding from the F2 generation of B6 × D2 crosses has resulted in 26 viable strains that show unique recombinations of B6 and D2 alleles: on average, three to four crossovers per chromosome have been preserved in each RI strain. All members of each RI strain are homozygous for either the B6 or the D2 allele at each locus. Twenty-five RI strains were available in sufficient number for testing in Experiment 1, and 24 strains for Experiment 2. All procedures were approved by the VA Institutional Animal Care and Use Committee and performed in accordance with National Institutes of Health guidelines for the care and use of animals in research.
Experiment 1
Assessment of stereotyped chewing, home cage locomotor activity, and temperature changes induced by METH. One of the primary aims of this study was to characterize the alterations in thermoregulation after METH administration in the BXD mice. Pilot studies had determined that repeated temperature assessment in B6 and D2 mice (using a rectal probe) resulted in lower temperatures than a single temperature determination at a comparable time point. Further, this effect was largely strain-dependent. To circumvent any confounds caused by repeated testing, we used only a single temperature measurement per mouse. We also chose to use behavioral measures that involved neither handling of animals nor novel environments and, thus, they were unlikely to affect body temperature determinations. These observational measures were exophthalmos (eye protrusion), chewing, and home cage horizontal activity. The analyses of these behaviors were conducted in 862 male mice of 25 RI strains and both progenitor strains. The N per strain per dose was 8 in almost all cases. The measurement of exophthalmos proved to be unreliable (reflected in a nonsignificant split-half reliability estimate ofr = 0.36) and, thus, genetic influences on this phenotype could not be evaluated. Therefore, the observation of exophthalmos in these mice was excluded from our analysis and will be omitted from the following discussion.
Each animal was taken from its home cage, weighed, injected with either SAL or METH (4, 8, or 16 mg/kg), and returned to its home cage with its normal cagemates. Only one mouse at a time in each cage was injected per experimental session, with four sessions per day. Starting 17 min after injection, activity was assessed for 1 min as the number of quadrant crossings in the home cage (each time an injected mouse entered a different quarter of the cage, 1 activity count was scored). Thirty-three minutes after injection, the injected subjects were again assessed for the number of chewing or gnawing occurrences for the next minute. Repetitive chewing involved paw-to-mouth movements most often directed at the corn cob bedding material, although frequently (∼25% of the time) it did not involve any object. Although food was present, it was never seen to be the object of chewing. Forty-eight minutes after injection, body temperature was determined with a rectal probe inserted 2 cm using an analog thermometer.
Measurement of brain METH concentration. Brain METH levels were determined by gas chromatography from samples taken immediately after the test for mice given either 4 or 8 mg/kg METH. Brains were excised and homogenized in 0.1 m potassium phosphate buffer, pH 6.0. A 550 ml aliquot was centrifuged at 13,000 rpm for 10 min, and 400 ml of the supernatant was transferred to a new tube. Internal standard and 600 ml of fresh buffer were added. The sample prep columns (Bond-Elut, 3 cc) were washed through the column, and the column was aspirated for 5 min. Next, 5 ml of methanol was washed through and the column was dried for 3 min. Finally, 1.5 ml of freshly prepared 2% ammonium hydroxide in ethyl acetate was slowly passed through the column and collected. A 500 ml portion of this material was placed in a glass vial along with 50 ml of 1% HCl in methanol, and the solution was dried. This was continued until all of the sample was dried. Next, 60 ml of heptafluorobutyl anhydride was added and the vial was sealed and heated at 70°C for 20 min. The vial was then cooled and dried, and 400 ml of ethyl acetate was added. The sample was sealed and vortexed, and 1 ml was injected in the GC/MS for analysis following the methods of Foltz et al. (1990).
Experiment 2
Assessment of stereotyped climbing induced by METH.In rodents, low doses of amphetamine (0.25–1.0 mg/kg) produce hypothermia and a general activation comprised of sniffing, locomotion, and rearing. Higher doses result in hyperthermia and a decreasing incidence of competing behaviors, as stereotyped behaviors begin to emerge (Randrup and Munkvad, 1967; Jellinek, 1971). Stereotypy may take many forms that depend in part on the animal species, dose of amphetamine, and time after injection, as well as environmental factors such as stress or previous pharmacological history with psychostimulants (Randrup et al., 1975). Rapid and repetitious paw-to-mouth movements (chewing) have been the most frequently studied behavior in rodents.
Protais et al. (1976) initially characterized and standardized climbing behavior as a stereotypy elicited by dopamine agonists in mice. After apomorphine administration, mice show increased locomotor activity and, if placed in a cage with mesh walls, enhanced climbing behavior (see, for example, Marcais et al., 1978; Balsara et al., 1981;Joshi et al., 1981; Cabib et al., 1995). A pilot study assessing the effects of several doses of METH determined that 16 mg/kg produced the most robust climbing response in both B6 and D2 mice. Therefore, 715 adult female mice from 24 of the BXD/Ty RI set—in addition to B6 and D2 mice—were tested on the climbing apparatus after intraperitoneal injection of either equivolume SAL or 16 mg/kg METH. Immediately after the injection, mice were placed in a single cage on corncob bedding surrounded by a 12-cm-diameter wire mesh chimney 24 cm high, constructed from 0.5 inch hardware cloth. The top of each chimney was covered to prevent escape. Animals were scored during a single 2 sec inspection each 5 min throughout a 60 min test period as follows (Protais, 1976): 4 paws on the floor (0), forefeet on the wall or grooming [i.e., 2 feet on floor (1)], 4 paws on the wall [climbing (2)].
Data analysis. Mean SAL strain values were subtracted from individual post-METH scores to obtain a measure of the METH-induced change in each behavior for Experiment 1. The effects of strain on METH sensitivity were then analyzed by one-way ANOVA at each dose for Experiment I. A single-factor (strain) ANOVA was also used to examine climbing in the RIs and their progenitors after 16 mg/kg METH in Experiment 2. Climbing scores were averaged over 15 min time periods to create four epochs. Because strains also differed in climbing behavior after SAL, METH sensitivity was again determined by subtracting each strain’s mean saline score for the corresponding time period from the individual post-METH scores within a strain. To determine the reliability of the measures assessed in these studies, split-half correlations were obtained for each trait and corrected using the Spearman–Brown method. Genetic correlations (i.e., among strain means) were estimated among the dependent measures examined in this study (Hegmann and Possidente, 1981).
QTL analysis. Over 1500 genetic markers (Mammalian Genome, 1996) have been genotyped in the BXD RI strains, mostly comprising microsatellite simple sequence repeat polymorphisms determined by PCR, so that the location and the B6- or D2-like repeat length of each of these alleles is known for each strain. By convention, D2 alleles are coded with values of 1, and B6 alleles with values of 0 in our database. Therefore, to identify the putative location of relevant QTLs, the strain mean values for chewing, body temperature, home cage activity, and climbing were correlated with the value (0 or 1) at each of the polymorphic loci in our data set. For body temperature changes and activity, each of the three doses was considered as a separate trait for QTL analysis. Thus, a total of eight behavioral traits were analyzed. A statistically significant correlation at p≤ 0.01 between the phenotype and the allelic state at a particular marker indicates that a gene in this region of the genome may influence the response to METH. A significant clustering of markers showing similar correlations on a particular region of the genome is more likely to point to a true locus of influence. Even so, primarily because of the large number of correlations calculated, there is a high probability of false positives. However, making the criteria more stringent would increase the likelihood of type II errors or false negatives (missing actual influential loci). Our strategy, therefore, is to use this analysis as an initial screen, suggesting provisional regions of influence that await confirmation or rejection in subsequent analyses. Our general approach for subsequent verification of the QTLs provisionally reported here has been discussed previously (Belknap et al., 1993, 1996a,b; Crabbe et al., 1994) (K. Buck, P. Metten, J. Belknap, and J. Crabbe, unpublished observations). We plan to examine other populations to follow up the results reported here with the goal of attaining aggregate levels of statistical significance required to accept linkage that have been proposed by Lander and Kruglyak (1995).
RESULTS
Dose-effect analysis in several METH-induced behaviors
Significant dose-related effects of METH on all responses assessed were evident (data not shown). Repetitive chewing showed a roughly linear dose–response function; the 16 mg/kg dose resulted in near 100% frequency for chewing for almost all strains. In contrast, body temperature and home cage activity responses evidenced biphasic curves. The lowest dose of METH, 4 mg/kg, caused statistically significant hypothermia, whereas the highest dose, 16 mg/kg, caused significant hyperthermia. Similarly, 4 mg/kg METH resulted in peak locomotor activation, whereas higher doses showed progressive decreases in home cage activity.
Strains differed widely with respect to METH-induced effects. For instance, B6 mice showed no effect of METH on body temperature at any dose, whereas others (such as BXD-1) showed consistent hyperthermia at all doses (consistent with Seale et al., 1986). Figure 1shows the METH-induced (METH − mean SAL) body temperature changes in all 27 strains at each of the three doses. At the 4 mg/kg dose (top panel), there was wide variation in body temperature changes (F(26,185) = 5.31,p < 0.0001). Eight strains showed significant hypothermia, and two were hyperthermic relative to saline control values. At the highest dose (16 mg/kg; bottom panel), there were also differences among strains (F(26,188) = 3.55, p < 0.0001): none of the strains were hypothermic, and ∼75% showed significant hyperthermia. The 8 mg/kg dose (middle panel) was intermediate but also variable (F(26,190) = 6.25, p < 0.0001). Overall, METH-induced body temperature changes were quite heritable, with genotype (strain) accounting for 43, 46, and 33% of the total observed variance (R2) at the 4, 8, and 16 mg/kg doses, respectively.
Changes in body temperature in 25 BXD RI strains and the two progenitor strains after 4 mg/kg (top), 8 mg/kg (middle), or 16 mg/kg (bottom) intraperitoneal METH injection. The strains are plotted in order from the most negative response (hypothermia) on the left to the most positive response (hyperthermia) on the right. Each bar shows the mean ± SEM for each dose. The strains along the x-axis refer to the BXD strain number; B6 and D2 mice are also labeled.
The typical dose–response curve with respect to home cage activity is an inverted U-shaped function. The individual strain means for home cage horizontal activity after the three doses of METH are shown in Figure 2. The 4 mg/kg dose (top panel) induced more activation than the other doses and varied significantly among strains (F(26,186) = 4.37, p < 0.0001). All strains were significantly activated by this dose of METH except for BXD-8 mice. Increases in home cage locomotor activity were also commonly observed after 8 mg/kg METH and, again, they differed significantly among strains (F(26,190) = 6.73,p < 0.0001; middle panel). In contrast, the 16 mg/kg dose (bottom panel) produced significant reductions in activity in nine strains and enhanced activity in only two strains (F(26,188)= 34.98, p < 0.0001). Again, METH-induced home cage activity was quite heritable, with genotype accounting for 38, 48, and 41% of the total variance for the 4, 8, and 16 mg/kg doses, respectively.
Mean ± SEM number of quadrant crossings/min (home cage activity) for B6, D2, and 25 BXD RI strains after 4 mg/kg (top), 8 mg/kg (middle), or 16 mg/kg (bottom) intraperitoneal METH. See legend to Figure1.
Repetitive chewing was almost never seen in animals treated with the 4 mg/kg dose of METH, whereas chewing was near maximal in most strains after the 16 mg/kg dose (data not shown). Thus, only the 8 mg/kg dose was particularly useful in assessing strain differences in sensitivity to this behavior (Fig. 3). As in the other behaviors, repetitive chewing showed large strain differences, from no response in three strains to continuous chewing throughout the observation period in three other strains (F(26,190) = 3.51,p < 0.0001). Genotype accounted for 33% of the total observed variance at this dose. Stereotypy was also measured by examining the climbing response in Experiment 2.
BXD RI strain chewing responses (mean ± SEM) after 8 mg/kg METH. See legend to Figure 1.
Spearman–Brown-corrected split-half reliability estimates indicated that all of these measures were reliable. For the temperature changesrSB = 0.81, 0.76, and 0.78 after 4, 8, and 16 mg/kg, respectively. Activity changes were also reliable:rSB = 0.75, 0.87, and 0.86 at the same doses. Finally, the corrected reliability estimate for chewing wasrSB = 0.85.
Brain METH concentrations
Although brain METH levels did vary significantly among strains at the 4 mg/kg dose (F26,183 = 1.724,p < 0.05), no significant differences in brain METH levels were evident after 8 mg/kg METH (p > 0.10) (Fig. 4). The correlation between the strain means for brain METH concentrations after these two doses was 0.61 (p < 0.001). However, brain concentrations of METH were not significantly correlated with any of the METH sensitivity measures taken in these same animals. Thus, because pharmacokinetic differences did not account for genotypic differences in drug response, we did not correct for brain METH concentrations when evaluating the effects of METH.
Brain METH levels (μg/gm) after 4 mg/kg (hatched bars) or 8 mg/kg (solid bars) METH in B6, D2, and 25 BXD RI strains. See legend to Figure 1.
Effect of 16 mg/kg METH on stereotyped climbing in BXD RIs
As with D2 mice, in many of the RI strains there was little evidence of stimulation during the initial 15 min period after injection. That is, for most RI strains, METH-treated groups did not differ initially from the SAL group and, therefore, data from the first epoch were excluded from further analysis. An initial analysis considered strain differences in SAL treatment across the latter three epochs using a two-way ANOVA. Baseline (SAL) scores were significantly variable across strains (F25,328 = 4.66,p < 0.0001), although SAL climbing did not vary across epochs (F2,656 = 0.449, p = 0.64) and the effect of strain did not depend on the time epoch examined (F50,656 = 0.92, p = 0.63). Therefore, data at each time point were indexed as the mean difference score of all the individual METH scores within a strain from the mean SAL score.
Because the strain-specific response to METH was highly correlated across epochs, a total climbing score was calculated. Figure5 shows that there was also wide variation within the BXD panel for the combined METH-induced difference scores (F25,334 = 13.52, p < 0.0001). Twelve strains responded with stimulated climbing, whereas six others showed decreased climbing and the remaining eight strains did not significantly alter their climbing activity after METH administration. The Spearman–Brown-corrected reliability estimate for the climbing change after 16 mg/kg METH was 0.79.
Total METH-induced climbing scores (±SEM) for all BXD RI strains for the period from 15 to 60 min after 16 mg/kg METH. See legend to Figure 1.
Genetic correlations among responses to METH
Table 1 shows correlations among strain responses to METH. Because they were based on strain means, they are predominantly genetic in origin (Hegmann and Possidente, 1981). Within the temperature and home cage activity measures, there was evidence for a high degree of genetic control in common across doses. For activity, the 4 and 8 mg/kg doses correlated r = 0.77, and the 8 and 16 mg/kg doses correlated r = 0.49. Similarly, temperature correlations across doses ranged from 0.49 to 0.84. The strains most activated by 4 or 8 mg/kg METH had greater increases (or lesser decreases) in METH-induced body temperature induced by all three doses. This suggests that some genes affect both thermal and locomotor responses to METH. The other behavioral responses were not generally correlated with each other, indicating that they are largely subserved by different mechanisms and/or neural substrates. An exception to this generalization was the significant genetic correlation between hyperthermia and climbing (r = −0.44) at the 16 mg/kg dose. Interestingly, the two measures thought to index forms of stereotyped behavior, chewing and climbing, tended to correlate negatively (r = −0.32, not significant), suggesting that chewing and climbing may be competitive behaviors.
Genetic correlations (among strain means) for behavioral measures taken in Experiments 1 and 2
In contrast to all other measures, climbing was evaluated in female mice (both genders were not evaluated for all phenotypes because of the limited availability of these strains). It is possible, therefore, that some of the putative common loci may in fact be gender-specific and that these correlations would change if one sex were examined for all measures. Combining data from mice of both sexes is perhaps most likely to result in an underestimation of common genetic etiology; furthermore, in pilot studies examining METH-induced climbing in both female and male B6 and D2 mice, strain differences did not depend on gender. Therefore, although it is certainly important to note this difference between the measures in Table 1, gender effects are unlikely to have a large impact on our interpretation of the data.
QTL analysis of METH sensitivity
To identify the putative location of relevant QTLs, the strain mean values for chewing, home cage activity, body temperature, and climbing were correlated with the value (0 or 1) at each of the polymorphic loci in our database. Multiple QTLs were provisionally identified (p ≤ 0.01; Table 2) for each measure of drug sensitivity. In several instances, possible candidate genes map to nearby locations. For each chromosome region, only the marker showing the highest correlation with each drug sensitivity measure is shown, but a list of all significant associations is available from the authors. The locations of provisionally mapped QTLs are shown schematically in Figure6.
Genetic markers, chromosome number, distance from centromere in cM (e.g., 1:8.7), and associations with METH responses for home cage activity, body temperature, chewing, and climbing
Schematic representation of genomic locations for putative QTLs underlying response to METH in the BXD RI strains and their progenitors. The 19 mouse chromosomes and the X chromosome are depicted proportionately; scale gives centimorgans (cM) from the centromere. For home cage activity (squares) and temperature changes (triangles), dose (mg/kg) is indicated by number inside the symbol. Markers from Silver et al. (1994).
DISCUSSION
Genetic factors mediate some of the variability between individuals in response to psychostimulants (c.f. Broadhurst, 1978;Seale, 1991). We have analyzed several responses to METH in the BXD RIs and their progenitors in order to identify putative genetic loci contributing to METH sensitivity. A comparison of the pattern of METH-induced strain mean effects with polymorphic genetic markers resulted in the identification of >70 provisional loci over all eight traits, or about eight to nine per trait. These can be used to generate hypotheses about candidate genes as well as to focus future testing in other genetic models. The next phase of the analysis will be to determine which of these provisional loci represent actual associations with high probability (Lander and Kruglyak, 1995). Computer simulations (Belknap et al., 1996a,b) and our experience (K. J. Buck, P. Metten, J. K. Belknap, and J. C. Crabbe, unpublished observations) suggest that approximately half of these loci will prove to represent true associations.
Candidate neurotransmitter modulation
Extensive pharmacological and neuroanatomical evidence that has been accumulating since the 1970s makes a strong case for the role of dopamine in psychostimulant-produced behaviors (see, for example,Groves and Rebec, 1976; Kuczenski, 1986; Hernandez et al., 1987). Congruent with these hypotheses, both temperature and activity changes induced by METH were associated with a region on chromosome 7 near the dopamine D4 receptor gene (Drd4). Locomotor activation in response to 4 mg/kg METH was also associated with a region on chromosome 16 that contains the gene for both the D3 receptor and the enzyme catechol-O-methlyltransferase (Comt), which is important for the catabolism of catecholamines. Even though genes for several dopamine receptor subtypes and the dopamine transporter have been mapped in mice, it is perhaps noteworthy that no QTLs suggested by analyses of either stereotyped behavior were found in the vicinity of these or other known dopamine-related genes. This finding does not necessarily imply that dopamine neurotransmission is not involved in METH-induced stereotypy but, rather, that genetic variability in dopamine activity in the BXD RIs may not underlie genetic differences in METH-induced stereotypy (as measured in the present studies). There is a dissociation between dopamine activation in the striatum induced by psychostimulants, and behavioral stereotypy, because there does not appear to be a strict temporal correlation between the amount of dopamine release in either the caudate nucleus or the nucleus accumbens with the stereotypic stimulation observed (Kuczenski et al., 1989, 1990, 1991). Our data are congruent with this interpretation and go further to suggest that the genetic variability in the chewing and climbing stereotypical responses depend more on other neural substrates which may, in turn, interact with dopamine mechanisms.
For instance, there is compelling evidence for modulation of METH-stimulated behaviors by acetylcholine and glutamate, neurotransmitters for which loci associated with relevant genes were identified in our analysis. Some of these areas are particularly attractive loci for attempts at verification because they overlap between studies of related phenotypes. For example, QTLs on mid-chromosome 1, associated with climbing, home cage activity, and temperature changes, are in a region where QTLs were provisionally mapped for open-field locomotor stimulation (horizontal movement) after 1 and 2 mg/kg METH and its sensitization (T. J. Phillips, M. Huson, C. Gwiazdon, and C. McKinnon, unpublished observations); this region contains genes coding for both the δ and the γ subunits of the nicotinic acetylcholine receptor. Moreover, a region on chromosome 9 that appears to be associated with the activity and temperature changes after METH in the present studies—also identified by the Phillips group as associated with cocaine-induced locomotor stimulation and METH sensitization—contains the gene for the a3 subunit of the acetylcholine receptor. The gene coding for the a2 subunit is in a region on proximal chromosome 14 with markers associated with a cluster of METH-related phenotypes (discussed below). There is also substantial biochemical and histochemical evidence suggesting that B6 and D2 mice differ with respect to cholinergic activity (Ebel et al., 1973; Durkin et al., 1977) and cell density (Albanese et al., 1985;Iacopino et al., 1986; Schwab et al., 1990). Furthermore, the genetic variability in cholinergic cell number in the striatum was shown byDains (1995) to be correlated with neuroleptic-induced catalepsy in the mouse, suggesting an interaction between acetylcholine and dopamine in the production of striatally mediated motor behavior. Although the relationship between striatal dopamine and acetylcholine has been well documented (Graybiel, 1990; Florin et al., 1992; Groves et al., 1995), the joint dopamine/acetylcholine responses to amphetamine do not appear to be sufficient to explain the behavioral changes observed, implying likely interactions with still other neurotransmitter systems (Kuczenski et al., 1977).
It is clear also that both glutamate and GABA are also involved in the locomotor activation resulting from dopamine agonist administration (Balsara et al., 1981; Liljequist and Karcz-Kubicha, 1993; Karler et al., 1995). For example, Karler et al. (1994, 1995) found that pharmacological blockade of dopamine, glutamate, or GABA receptors in the ventral lateral striatum could prevent the stereotypy induced by amphetamine or cocaine. The dopaminergic projections from the substantia nigra pars compacta are known to involve both of these other transmitters in the striatal–thalamic–cortical loop thought to underlie the locomotor activating effects of psychostimulants. Our data suggest that glutamate modulation by the enzyme glutamate dehydrogenase may play a role in determining individual differences in genetic sensitivity to activation by amphetamine, because variation in a marker near the Glud locus (and the acetylcholine a2subunit) on chromosome 14 was correlated with genetic variability in four observed traits (climbing and body temperature changes after 16 mg/kg METH and home cage activity after either 8 or 16 mg/kg METH). Given the consensus among genetic studies and the behavioral data suggesting a role for glutamate in responsivity to psychostimulants, further testing of this locus in an F2 intercross might be a fruitful endeavor.
QTLs affecting multiple phenotypes
There are three other published genetic analyses of psychostimulant responses in the BXD RIs. Data collected by Seale et al. (1985) were subject to QTL analyses by Gora-Maslak et al. in 1991. The hyperthermia induced by 20 mg/kg amphetamine was found to be associated with markers in four distinct regions of the mouse genome. These were from 10 RI strains and a very limited set of less than 200 markers, constraints that would imply that any associations observed would be likely to reflect fairly large genetic influences but would also have a higher probability of being spurious. In fact, all four of these regions overlap with those suggested for some traits by the current analyses, although we found QTLs for METH’s temperature effects in common with Gora-Maslak et al. (1991) only in the region from 27 to 31 cM on chromosome 6. Recently, Tolliver et al. (1994) andMiner and Marley (1995) mapped QTLs contributing to the locomotor activation in the BXD RIs after cocaine administration. Tolliver and colleagues evaluated the response to 32 mg/kg cocaine in 16 strains, while Miner and Marley examined locomotor stimulation to 10 mg/kg cocaine in 11 of the BXDs. Both Tolliver et al. (1994) and Miner and Marley (1995) identified the region of chromosome 9 associated with the locomotor response that overlaps with the large cluster of QTLs that we found associated with several responses to METH on the distal portion of this chromosome (discussed above). All three studies also found a common region of influence near the D17 Mit3 andHprt-ps1 loci on chromosome 17. As gene mapping studies of drug effects begin to accumulate, both discrepancies and commonalties between studies are potentially informative. The fact that common regions are found in studies of phenotypes that differ substantially (large dose ranges or different pharmacological agents) can lead to hypotheses about common mechanisms (see, for example, Crabbe et al., 1994). The Hprt-ps1 marker locus on chromosome 17, for instance, was associated with psychostimulant response in all of these studies. This is a pseudogene for the enzyme hypoxanthine-guanine phosphoribosyltransferase (HPRT) and, therefore, does not appear to translate any functional proteins (Isamat et al., 1988). It is interesting, however, that an HPRT enzyme deficiency is associated with the motor disorder Lesch-Nyhan syndrome in humans and with alterations in dopaminergic activity in motor nuclei and increased sensitivity to amphetamine-induced locomotion and stereotypy in mice (Jinnah et al., 1991, 1994), suggesting the possibility that polymorphisms in the retroposed HPRT pseudogene may differentially regulate (as endogenous antisense nucleotides) the HPRT gene itself (see, for example, Zhou et al., 1992; Akagi et al., 1994; Carlton et al., 1994; Manjanatha et al., 1994; Kas et al., 1995; Leedman et al., 1995).
A lack of concurrence between QTLs for different markers may be indicative of distinct genetic underpinnings for the different yet related phenotypes, or they may be attributable to other factors. It is possible, for instance, that more pleiotropic QTLs would have been identified in common if the same sex were tested for all phenotypes; as it happens, female mice were analyzed in the climbing study, whereas male mice were used to evaluate all other phenotypes. In a recent study by Melo et al. (1996) aimed at mapping genes associated with alcohol preference, two sex-specific QTLs were identified. Because these are the first gender-specific QTLs proposed for a pharmacological response, the prevalence of such sex-constrained loci is not yet known. However, many behavioral responses are known to vary significantly between females and males and, therefore, attempts at elucidating genetic substrates of such phenotypes should give this possibility serious consideration.
One of the ubiquitous effects of addictive drugs is behavioral stimulation (see Crabbe and Harris, 1991). This increased activity seen in both animals and humans after drug administration may reflect aspects of the reinforcing or rewarding effects of these drugs (Wise, 1978, 1987). The current use of QTL mapping in RI strains offers promise toward the identification of genes associated with locomotor activation in response to psychotropic drugs. The present findings suggest putative QTLs that may contribute to differences in METH sensitivity between individuals. The next step is to attempt to confirm these provisional QTL associations in an F2 (or other) population. In this way, offspring from B6 and D2 parents are individually phenotyped and then genotyped at areas of the genome suggested to reflect influences of QTLs from this RI study. When the correlations between behavior and allelic variation are confirmed in such a genetically segregating population, the case for a significant influence of these loci is strengthened further. Finally, congenic mice, isogenic at all loci except for a small region bracketing a QTL in question, can be created, contributing further to our understanding of the influence of a particular locus or gene product on a complex phenotype. In the interim, pharmacological or anatomical interventions can be used to test the hypothesis generated by preliminary QTL analysis in RIs (such as the role of glutamate or acetylcholine for genetic sensitivity to psychostimulants). These convergent methods of testing are likely to further the goal of understanding the precise molecular products that influence responding to psychostimulants and, therefore, susceptibility to drug abuse.
Footnotes
These studies were supported by Public Health Service Grants POl AA08621 and P50AA10760, National Institutes of Drug Abuse Contract 271-90-7405, and two grants from the Department of Veterans Affairs. J.E.G. was supported by Training Grant T32 AA07468. A preliminary report of this work has been given at the Special Conference on The Biology and Genetics of Complex Mammalian Traits, 1995, Bar Harbor, ME. We thank Dr. Steve R. Mitchell for his help with the QTL analyses, Dr. Jeffrey S. Mogil for his suggestions on the preparation of this manuscript, and Dr. Charlie Meshul for his helpful discussions.
Correspondence should be addressed to Judith E. Grisel, Research Service, 151W, Department of Veterans Affairs Medical Center, 3710 SW U.S. Veterans Hospital Road, Portland, OR 97201.
