Taar1 gene variants have a causal role in methamphetamine intake and response and interact with Oprm1.
Free full text in Europe PMC
Abstract
Free full text
Taar1 gene variants have a causal role in methamphetamine intake and response and interact with Oprm1
Abstract
We identified a locus on mouse chromosome 10 that accounts for 60% of the genetic variance in methamphetamine intake in mice selectively bred for high versus low methamphetamine consumption. We nominated the trace amine-associated receptor 1 gene, Taar1, as the strongest candidate and identified regulation of the mu-opioid receptor 1 gene, Oprm1, as another contributor. This study exploited CRISPR-Cas9 to test the causal role of Taar1 in methamphetamine intake and a genetically-associated thermal response to methamphetamine. The methamphetamine-related traits were rescued, converting them to levels found in methamphetamine-avoiding animals. We used a family of recombinant inbred mouse strains for interval mapping and to examine independent and epistatic effects of Taar1 and Oprm1. Both methamphetamine intake and the thermal response mapped to Taar1 and the independent effect of Taar1 was dependent on genotype at Oprm1. Our findings encourage investigation of the contribution of Taar1 and Oprm1 variants to human methamphetamine addiction.
eLife digest
People who misuse drugs often do so partly in response to the environment they find themselves in, and partly because of their genetics. The genetic component of someone’s risk is influenced by many different genes, and most research has found that each gene has a small individual effect. A method called quantitative trait locus (QTL) analysis can help find parts of the genome that influence someone’s risk of misusing drugs. In 2013, researchers found one region on chromosome 10 in mice has a particularly large influence on how much methamphetamine an individual mouse will ingest if the drug was available in one of its two water bottles. A gene called Taar1 was particularly important in this region and another gene, called Oprm1, may also play a significant role.
When the Taar1 gene is switched off, mice consume larger amounts of methamphetamine, have a heightened reward response from the drug, and are insensitive to the adverse effects – such as hypothermia. But whether Taar1 directly caused these effects, and whether Taar1 and Oprm1 interact, had not yet been determined. If these genes played a causal role, they could be useful targets for treatment of methamphetamine-use disorder.
Stafford, Reed et al. – who include several of the researchers involved in the 2013 work – now report that when a particular variant of Taar1 was present in mice they consumed large amounts of methamphetamine. The variant codes for a faulty version of a receptor protein. When this variant was replaced with a working version using gene editing, the mice consumed less methamphetamine and also became sensitive to hypothermia induced by the drug. This confirms that this gene does play a causal role in methamphetamine consumption and hypothermia. Next, Stafford, Reed et al. tested mice with different combinations of variants of Oprm1 and Taar1 to see how the genes interacted. The results showed that the effects of Taar1 on both consumption of the drug and hypothermia depended on the Oprm1 variant present.
The findings suggest that variants of these two genes in humans could influence an individual’s risk of addiction to methamphetamine. It is possible that in future the disorder could be treated by drugs that modify the brain activity impacted by these receptors. But first, it will be important to find out if these genes play a similar role in humans as they do in mice.
Introduction
Amphetamine-like stimulants, including methamphetamine (MA), remain the second most common class of illicit drugs used worldwide, behind cannabis-containing drugs (United Nations Office on Drugs and Crime, 2016). MA is the most widely abused psychostimulant (Chomchai and Chomchai, 2015), and the consequences of MA addiction pose major societal and health concerns (National Institute on Drug Abuse, 2019). MA-associated deaths in the United States were relatively stable between 2005 and 2013, but have since been on an increasing trajectory, rising from about 3600 individuals in 2013 to almost 11,000 in 2017 (National Institute on Drug Abuse, 2018). Discovering and understanding genetic factors that contribute to MA addiction risk, and uncovering the linked molecular and cellular processes, may help to curb this trend and improve interventions and long-term treatment success.
Whereas use of large human cohorts is effective at mapping loci and nominating candidate genes, precise experimental control over genetic composition, development, and drug exposure history is practical only in animal models. Evidence that mouse models continue to contribute critical information to our understanding of pathological conditions was recently reviewed (Nadeau and Auwerx, 2019). Selected lines of rodents have been derived to model genetic risk factors for drug use. These include, particularly, those bred for various drug intake phenotypes, such as our mice bred for different amounts of voluntary MA intake (Hitzemann et al., 2019; Shabani et al., 2011; Wheeler et al., 2009). In previous studies, we demonstrated a strong genetic contribution to level of voluntary MA consumption using multiple replicate sets of mice that we bidirectionally, selectively bred for MA high drinking (MAHDR) and MA low drinking (MALDR), which for simplicity are also referred to here as the ‘high’ and ‘low’ lines. Multiple lines of evidence support the involvement of the trace amine-associated receptor 1 gene, Taar1, in differences in MA intake, including: quantitative trait locus (QTL) mapping in the high and low lines (Belknap et al., 2013); the direct agonist activity of MA at the trace amine-associated receptor 1 (TAAR1) (Bunzow et al., 2001; Wolinsky et al., 2007); the presence of a single nucleotide polymorphism (SNP) in Taar1 (Taar1m1J) that renders the receptor it expresses (TAAR1) non-functional in response to MA and to other direct agonists (Harkness et al., 2015; Shi et al., 2016); fixation of the Taar1m1J allele in high line mice of all five replicates, and retention of the wildtype, Taar1+ allele in all low line replicates (Harkness et al., 2015; Reed et al., 2018); greater level of MA intake in Taar1 knockout mice compared to wildtype mice (Harkness et al., 2015); strong association of Taar1 genotype with MA intake in multiple genetic models (Reed et al., 2018); and strong association of Taar1 genotype with additional MA-related traits that correspond with level of MA intake, including level of sensitivity to MA-induced hypothermia and conditioned taste aversion (Harkness et al., 2015; Reed et al., 2018).
In the current studies, to demonstrate a causal role for Taar1, we used CRISPR-Cas9 technology to replace the Taar1m1J allele found in high line mice with the Taar1+ allele, creating a knock-in with functional TAAR1. These mice were then tested for MA consumption in comparison to an unaltered high line control, with the prediction that knock-in of Taar1+ would reduce MA intake. In addition, we tested knock-in and control mice for MA-induced change in body temperature to determine if this gene has a pleiotropic effect on this MA trait. We predicted that insertion of the Taar1+ allele would increase sensitivity to MA-induced hypothermia in the knock-in mice.
Our existing data indicate that Taar1 genotype has a major impact on MA consumption, and gene mapping results demonstrate that a chromosome 10 QTL, at the location of Taar1, accounts for as much as 60% of the genetic variance in risk for MA intake (Belknap et al., 2013). Other undefined variants must account for the remaining genetic variance. We previously examined the µ-opioid receptor 1 gene, Oprm1, as a quantitative trait gene (QTG) residing at the proximal end of mouse chromosome 10, but it was eliminated as a risk factor for MA intake (Eastwood et al., 2018). Although Oprm1 variants do not predict risk, Oprm1 is involved in MA consumption via regulation by a top-ranked transcription factor network of genes that is differentially expressed between the MA drinking lines. We found that Oprm1 and its product, the µ-opioid receptor (MOP-r), are more highly expressed in the prefrontal cortex of low line mice, relative to high line mice (Belknap et al., 2013; Eastwood et al., 2018). In addition, when MOP-r agonists were given to high line mice to test the hypothesis that increased MOP-r-regulated activity reduces MA intake, the acquisition of MA intake was attenuated, simulating the voluntary intake phenotype of the low line (Eastwood et al., 2018; Eastwood and Phillips, 2014a). High line mice possess the same mutant Taar1m1J/m1J genotype as their progenitor DBA/2J (D2) strain and are more likely to also possess the Oprm1 variant from the D2 strain, due to the relatively close proximity of Taar1 and Oprm1 on chromosome 10 (17 Mb apart). Due to this linkage, we have been unable to explore independent and interactive effects of the Taar1 and Oprm1 genotypes on MA intake or other MA-related traits in the MA drinking lines. Instead, we explored this in the present studies, using a family of related recombinant inbred mouse strains derived from a C57BL/6J (B6) x D2 inbred strain intercross (collectively the BXD strains) in which the linkage is disrupted. We examined the independent and interactive effects of these genes on MA intake and sensitivity to MA-induced change in body temperature, and performed QTL mapping for both traits using the BXD strain data.
Results
Knock-in of Taar1+ converts high MA intake to low MA intake
We predicted that MA intake would be reduced in high line mice in which the Taar1m1J allele was removed and replaced with the Taar1+ allele. In the initial repeated measures analyses of variance (ANOVAs), there were no interactions involving sex for either mg/kg MA consumption or total volume consumed, so we performed additional analyses with data for the sexes combined. For MA consumption (Figure 1a; Figure 1—source data 1), there was a significant MA concentration x genotype interaction (F[1,36]=66.8, p<0.0001). Replacement of the Taar1m1J allele with the Taar1+ allele converted high MA intake to low MA intake for both MA concentrations. Control mice, which are homozygous Taar1m1J, consumed significantly more MA at the 40 mg/l concentration, compared to the 20 mg/l concentration, whereas low and comparable levels of MA consumption were found for both MA concentrations in knock-in mice. For total volume consumed from the water and MA tubes, during the time that MA was available (Figure 1b; Figure 1—source data 1), there was only a significant main effect of MA concentration (F[1,36]=21.4, p<0.0001). Total volume consumed was greater when mice were offered 40 mg/l MA, compared to 20 mg/l MA.
(a) MA consumption and (b) total volume consumed, when MA was offered vs. water for 18h/day at a concentration of 20 or 40 mg/l for control and knock-in mice. Presented are means (represented by bars) and individual data points (represented by circles) for data collapsed on sex. n = 20/genotype tested in a single cohort. Repeated measures ANOVA followed by simple main effects analysis, ***p<0.001 vs. control (MAHDR-Taar1m1J/m1J) at the same MA concentration; +++p < 0.001 vs. control (MAHDR-Taar1m1J/m1J) at 20 mg/l MA. Repeated measures ANOVA, ###p<0.001 for the main effect of 20 mg/l vs. 40 mg/l MA. MA, methamphetamine; MAHDR, MA high drinking mice; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele. The raw data represented in the graphs are available in Figure 1—source data 1.
Figure 1—source data 1.
MA consumption (mg/kg/18h) and total volume consumed (ml/18h) data for male and female, control MAHDR-Taar1m1J/m1J and knock-in MAHDR-Taar1+/+ mice.Knock-in of Taar1+ restores a hypothermic response to MA in high line mice
Thermal response to MA is a genetically-correlated response to selection for level of MA intake; thus, low MA drinking line mice display MA-induced hypothermia, whereas high line mice do not (Harkness et al., 2015). We tested the hypothesis that Taar1 genotype has a causal role. Knock-in of Taar1+ in high line mice restored a hypothermic response to MA that did not occur in control mice (Figure 2; Figure 2—source data 1). In the initial repeated measures ANOVA there were significant independent interactions of sex with genotype (F[1,90]=7.1, p=0.009), MA dose (F[1,90]=4.9, p=0.03) and time (F[4,360]=2.9, p=0.02). However, sex did not determine the pattern of response to saline vs. MA in each genotype across time. There was a significant 3-way genotype x MA dose x time interaction (F[4,360]=10.2, p<0.0001) that did not interact with sex. We next examined the effects of genotype and time within each dose group. In the saline-treated mice (Figure 2a), there were significant main effects of time (F[4,188]=35.9, p<0.0001) and genotype (F[1,47]=4.7, p=0.04). Temperature declined by 0.3°C over time, regardless of genotype, and knock-in mice had 0.25°C higher body temperature than control mice, regardless of time. In the MA-treated mice (Figure 2b), the genotype x time interaction was significant (F[4,188]=13.0, p<0.0001). Knock-in mice had a significant hypothermic response to MA at T30 compared to T0, whereas control mice were resistant to the hypothermic effects of MA, instead exhibiting significant MA-induced increases in body temperature at all post-MA treatment time points, compared to T0. Body temperatures were significantly higher in knock-in than control mice at baseline, by about 0.6°C, but significantly lower by about 1.2°C, at T30.
Body temperatures obtained immediately before (T0) or 30–120 min (T30-120) after (a) saline or (b) 2 mg/kg MA treatment for control and knock-in mice. Presented are means ± SEM for data collapsed on sex. n = 24–25/genotype/MA dose tested in 3 cohorts of 25–46 mice. Repeated measures ANOVA, #p<0.05 for the main effect of genotype. Repeated measures ANOVA followed by simple main effects analysis, ***p<0.001 for the effect of genotype at the indicated time point; Dunnett’s post hoc test,+++p < 0.001 compared to T0. MA, methamphetamine; MAHDR, MA high drinking mice; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele. The raw data represented in the graphs are available in Figure 2—source data 1.
Figure 2—source data 1.
Core body temperature (°C) data across time for saline and 2 mg/kg MA-treated male and female, control MAHDR-Taar1m1J/m1J and knock-in MAHDR-Taar1+/+ mice.MA consumption is impacted in BXD mice by a Taar1 X Oprm1 interaction
Initial analyses characterized 18 BXD strains for MA consumption (Figure 3; Figure 3—source data 1). There was a significant MA concentration x strain interaction (F[17,215]=17.2, p<0.0001). MA intake differed among the strains for both the 20 (Figure 3a; p<0.0001) and 40 (Figure 3b; p<0.0001) mg/l MA concentrations, and all Taar1m1J/m1J strains, but not Taar1+/+ strains, exhibited higher levels of MA intake from the 40 mg/l MA concentration than from the 20 mg/l concentration (p=0.0006 to<0.0001). For total volume consumed when mice were offered water vs. MA, there was a main effect of strain (F[17,215]=27.7, p<0.0001) and a main effect of MA concentration (F[1,215]=58.6, p<0.0001). Total volume consumed differed among the strains by as much as 4 ml (Figure 3c), but the strain distribution pattern for total volume did not correspond with the distributions for MA intake at either MA concentration. Total volume consumed was significantly less when water and 20 mg/l MA were available, compared to when water and 40 mg/l MA were available (mean ± SEM = 5.3±0.1 and 5.7 ± 0.1 for 20 and 40 mg/l, respectively).
MA consumption when MA was offered vs. water for 18h/day at a concentration of (a) 20 mg/l or (b) 40 mg/l. (c) Total volume consumed from the water and MA tubes during the 18h period when MA was offered vs. water, collapsed on MA concentration, because there was no concentration-dependent effect. Presented are means (represented by bars) and individual data points (represented by circles) for data collapsed on sex. n = 5–28/strain tested in 4 cohorts of 27–97 mice. Repeated measures ANOVA, ###p<0.001 for the main effect of strain. BXD, C57BL/6J (B6) x DBA2/J (D2); MA, methamphetamine; Oprm1, mu opioid receptor gene; Oprm1B6/B6, homozygous for B6 Oprm1 allele; Oprm1D2/D2, homozygous for D2 Oprm1 allele; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele. The raw data represented in the graphs are available in Figure 3—source data 1.
Figure 3—source data 1.
MA consumption (mg/kg/18h) and total volume consumed (ml/18h) data for male and female, BXD mice, according to Taar1 and Oprm1 genotype.We next examined the data from the BXD mice with regard to Taar1 and Oprm1 genotype to detect potential independent and epistatic effects on MA consumption. First, to identify potential independent effects, we calculated Pearson’s correlations of each genotype with individual 20 and 40 mg/l MA intake amounts. MA intake was significantly associated with Taar1 genotype (r = 0.74 and 0.81, ps <0.0001, for 20 and 40 mg/l MA, respectively), but not Oprm1 genotype (r = 0.06 and 0.10, p=0.36 and 0.13, for 20 and 40 mg/l MA, respectively). Next, we considered main and interaction effects by repeated measures ANOVA with MA concentration, Oprm1 genotype, Taar1 genotype, and sex as factors. There were no significant interactions involving sex, so further analyses were performed with data for the sexes combined. Data are summarized in Figure 4 (Figure 4—source data 1), with mice grouped according to their four possible Taar1/Oprm1 genotypes (Table 1). For MA consumption (Figure 4a), there was a significant Oprm1 genotype x Taar1 genotype x MA concentration interaction (F[1,225]=17.6, p<0.0001). For consumption of MA from the 20 mg/l MA concentration, main effects of the Oprm1 (F[1,229]=8.6, p=0.004) and Taar1 (F[1,229]=301.3, p<0.0001) genotypes indicated that both genes influenced MA intake, with greater MA intake in Taar1m1J/m1J and Oprm1D2/D2 mice, compared to Taar1+/+ and Oprm1B6/B6 mice, respectively. However, there was no significant Oprm1 genotype x Taar1 genotype interaction. In contrast, for MA consumption from the 40 mg/l MA concentration, there was a significant Oprm1 genotype x Taar1 genotype interaction (F[1,229]=13.9, p=0.0002). MA intake was significantly higher in Taar1m1J/m1J mice, compared to Taar1+/+ mice, regardless of Oprm1 genotype (ps <0.0001). However, Taar1m1J/m1J/Oprm1D2/D2 mice consumed significantly more MA than Taar1m1J/m1J/Oprm1B6/B6 mice. Oprm1 genotype did not significantly impact MA intake in Taar1+/+ mice.
(a) MA consumption and (b) total volume consumed, when MA was offered vs. water for 18h/day at a concentration of 20 or 40 mg/l for BXD strain mice with different combined Taar1/Oprm1 genotypes. Presented are means (represented by bars) and individual data points (represented by circles) for data collapsed on sex. n = 52–71/ genotype tested in 4 cohorts of 27–97 mice. Repeated measures ANOVA, ###p<0.001 for the main effect of Taar1+/+ vs. Taar1m1J/m1J. Repeated measures ANOVA followed by simple main effects analysis, ***p<0.001 compared to Taar1+/+ with the same Oprm1 genotype; +++p < 0.001 compared to Taar1m1J/m1J/Oprm1B6/B6 at 40 mg/l MA. BXD, C57BL/6J (B6) x DBA2/J (D2); MA, methamphetamine; Oprm1, mu-opioid receptor gene; Oprm1B6/B6, homozygous for B6 Oprm1 allele; Oprm1D2/D2, homozygous for D2 Oprm1 allele; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele. The raw data represented in the graphs are available in Figure 4—source data 1.
Figure 4—source data 1.
MA consumption (mg/kg/18h) and total volume consumed (ml/18h) data for male and female BXD mice, according to Taar1 and Oprm1 genotype.Table 1.
| Taar1/Oprm1 genotype combination | MA intake | MA-induced change in body temperature |
|---|---|---|
| Taar1+/+/Oprm1B6/B6 | BXD 154, 184, 196, 212, 218 | BXD 154, 184, 196, 218 |
| Taar1m1J/m1J/Oprm1B6/B6 | BXD 113, 160, 171, 194 | BXD 113, 160, 171, 194 |
| Taar1+/+/Oprm1D2/D2 | BXD 161, 199, 205, 215 | BXD 161, 199, 205 |
| Taar1m1J/m1J/Oprm1D2/D2 | BXD 178, 186, 204, 210, 216 | BXD 178, 186, 210, 216 |
Abbreviations: BXD, C57BL/6J (B6) x DBA2/J (D2); MA, methamphetamine; Oprm1, mu-opioid receptor gene; Oprm1B6/B6, homozygous for B6 Oprm1 allele; Oprm1D2/D2, homozygous for D2 Oprm1 allele; Taar1, trace-amine associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele.
For total volume consumed (Figure 4b), there were significant main effects of Oprm1 genotype (F[1,225]=24.8, p<0.0001), Taar1 genotype (F[1,225]=4.7, p=0.03), and MA concentration (F[1,225]=75.7, p<0.0001), but no significant interactions. Oprm1D2/D2 mice consumed more total volume than Oprm1B6/B6 mice, and Taar1m1J/m1J mice consumed more than Taar1+/+ mice. In addition, total volume consumed was greater when mice were offered 40 mg/l MA, compared to 20 mg/l MA.
MA-induced hypothermia is impacted in BXD mice by a Taar1 X Oprm1 interaction
We used a similar approach to analyze body temperature data that were collected after saline or MA treatment in 15 BXD strains (Figure 5; Figure 5—source data 1). For ease of viewing, the data in Figure 5 are separated by saline/MA treatment and Taar1 genotype for strain responses across time. A repeated measures ANOVA detected a significant strain x MA dose x time interaction (F[56,1172]=10.0, p<0.0001). In saline-treated mice (Figure 5a,b), the BXD strains differed at each time point (ps = 0.04 to<0.0001), except T0, and body temperature differed across time in strains 161 (Taar1+/+/Oprm1D2/D2), 171, 194 (both Taar1m1J/m1J/Oprm1B6/B6), 186, 210 (both Taar1m1J/m1J/Oprm1D2/D2), and 218 (Taar1+/+/Oprm1B6/B6) (ps = 0.004 to<0.0001). In MA-treated mice (Figure 5c,d), the BXD strains differed at each time point (ps <0.0001), except T0, and body temperature differed across time in all strains (ps = 0.007 to<0.0001), except 160, 171, 194 (all Taar1m1J/m1J/Oprm1B6/B6), 178, 186, and 210 (all Taar1m1J/m1J/Oprm1D2/D2). The maximum mean drop in body temperature in saline-treated mice was 0.9°C (Figure 5a,b), whereas there was a maximum mean drop of 4°C in MA-treated Taar1+/+ mice (Figure 5c), and a maximum increase of 1°C in MA-treated Taar1m1J/m1J mice (Figure 5d).
Body temperatures obtained just before (T0) or 30–120 min (T30-120) after injection for (a) saline-treated Taar1+/+ mice of either Oprm1 genotype, (b) saline-treated Taar1m1J/m1J mice of either Oprm1 genotype, (c) 2 mg/kg MA-treated Taar1+/+ mice of either Oprm1 genotype, and (d) 2 mg/kg MA-treated Taar1m1J/m1J mice of either Oprm1 genotype. Presented are means ± SEM for data collapsed on sex. n = 4–15/strain/dose tested in 16 cohorts of 8–36 mice. For clarity, significant findings are discussed in the text. BXD, C57BL/6J (B6) x DBA2/J (D2); MA, methamphetamine; Oprm1, mu-opioid receptor gene; Oprm1B6/B6, homozygous for B6 Oprm1 allele; Oprm1D2/D2, homozygous for D2 Oprm1 allele; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele. The raw data represented in the graphs are available in Figure 5—source data 1.
Figure 5—source data 1.
Core body temperature (°C) data across time for saline and 2 mg/kg MA-treated male and female, BXD mice, according to Taar1 and Oprm1 genotype.We next examined the BXD data with regard to Taar1 and Oprm1 genotype to detect potential independent and epistatic effects on the hypothermic effect of MA. Again, to identify potential independent effects, we calculated Pearson’s correlations of each genotype with MA-induced change in body temperature from T0 to T30. There was a significant correlation of body temperature change with Taar1 genotype (r = 0.71, p<0.0001), but not with Oprm1 genotype (r = 0.11, p=0.16). Next, we considered main and interaction effects by repeated measures ANOVA with time, Oprm1 genotype, Taar1 genotype, sex and MA dose as factors. There were significant independent interactions of sex with MA dose (F[1,307]=8.0, p=0.006) and time (F[4,1228]=5.0, p=0.0007). However, sex did not play a significant role in the pattern of response of each genotype to each dose across time. There was a significant Oprm1 genotype x Taar1 genotype x MA dose x time interaction (F[4,1228]=8.0, p<0.0001) that did not interact with sex. We next examined the effects of genotype and time within each dose group (Figure 6; Figure 6—source data 1). In the saline-treated mice (Figure 6a), the Oprm1 genotype x time interaction was significant (F[4,596]=2.6, p=0.04). Both Oprm1B6/B6 and Oprm1D2/D2 mice exhibited reductions in body temperatures at T60-120 (ps <0.001), compared to T0. However, these differences in body temperature were of no more than 0.3°C on average. In the MA-treated mice (Figure 6b), there was a significant 3-way interaction of Oprm1 genotype x Taar1 genotype x time (F[4,664]=11.4, p<0.0001). Taar1+/+/Oprm1D2/D2 mice had significantly lower body temperatures at all post-MA treatment time points, compared to Taar1+/+/Oprm1B6/B6 mice. Taar1+/+/Oprm1B6/B6 mice exhibited a significant hypothermic response to MA at T30-90, compared to T0, and recovered to baseline (T0) temperatures by T120, whereas Taar1+/+/Oprm1D2/D2 mice exhibited a significant hypothermic response to MA at T30-120, and remained below baseline for the duration of testing. However, the rate of increase after the initial decrease at T30 was similar for the two genotypes. Regardless of Oprm1 genotype, Taar1m1J/m1J mice were resistant to the hypothermic effects of MA, exhibiting increases in body temperature at all post-MA treatment time points, compared to T0.
Body temperatures obtained just before (T0) or 30–120 min (T30-120) after (a) saline or (b) 2 mg/kg MA treatment for BXD strain mice with different combined Taar1/Oprm1 genotypes. Presented are means ± SEM for data collapsed on sex. n = 30–47/genotype/MA dose tested in 16 cohorts of 8–36 mice. Repeated measures ANOVA followed by Dunnett’s post hoc test,+p < 0.05 compared to T0 of same genotype; +++p < 0.001 compared to T0 of same genotype. Repeated measures ANOVA followed by simple main effects analysis, ***p<0.001 compared to Taar1+/+/Oprm1B6/B6. Repeated measures ANOVA collapsed on Oprm1 genotype followed by Dunnett’s post hoc test, ###p<0.001 compared to T0 of Taar1m1J/m1J mice. BXD, C57BL/6J (B6) x DBA2/J (D2); MA, methamphetamine; Oprm1, mu-opioid receptor gene; Oprm1B6/B6, homozygous for B6 Oprm1 allele; Oprm1D2/D2, homozygous for D2 Oprm1 allele; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele. The raw data represented in the graphs are available in Figure 6—source data 1.
Figure 6—source data 1.
Core body temperature (°C) data across time for saline and 2 mg/kg MA-treated male and female, BXD mice, according to Taar1 and Oprm1 genotype.MA consumption and body temperature response to MA map to a region of chromosome 10 at the location of Taar1
The QTL results for both traits are displayed in Figure 7 (Figure 7—source data 1). A significant QTL in the same region of chromosome 10 emerged for both MA consumption (Figure 7a) and body temperature response to MA (Figure 7b). This QTL is at the location of the Taar1 SNP (23.9 Mb). The correlations between the strain means and the Taar1 SNP were r = 0.94, p<0.0001 for MA consumption and r = 0.82, p=0.02 for temperature response to MA. In addition, suggestive QTLs were detected on chromosomes 17 and 18 for MA effect on body temperature. There was a strong statistical trend for an increase in the chromosome 10 logarithm of the odds (LOD) score for MA intake (from 8.2 to 8.8; p=0.07), and a significant increase for the chromosome 10 LOD score for temperature response (from 3.7 to 4.2; p=0.01), when results were compared for interval mapping and composite interval mapping. Composite mapping identified a suggestive QTL on chromosome four for MA intake that had not reached the suggestive significance threshold in initial mapping; however, the LOD change from 2.5 to 2.7 was not statistically significant (Figure 7c). The LOD scores for the suggestive QTLs on chromosomes 17 and 18 for temperature response were comparable for initial and composite mapping (Figure 7d).
Means for 18 BXD strains were used for mapping of (a) consumption of MA from a 40 mg/l solution, and 15 strains for mapping of (b) body temperature response to 2 mg/kg MA at 30 min post-treatment. Composite interval mapping for MA consumption (c) and body temperature response (d) included Oprm1 genotype as a co-factor. Chromosome number and megabase position are represented along the x-axis (chromosome indicated at top of plot). The y-axis represents the LOD ratio, a measure of the strength of association between variation in the phenotype and genetic differences (alleles) at a particular chromosomal locus. The horizontal pink and gray lines indicate the significant and suggestive threshold levels, respectively. Mapping results were generated using the QTL mapping module of GeneNetwork (www.genenetwork.org). BXD, C57BL/6J (B6) x DBA2/J (D2); LOD, logarithm of the odds; MA, methamphetamine; QTL, quantitative trait locus. The raw data represented in the graphs are available in Figure 7—source data 1.
Figure 7—source data 1.
BXD means for QTL mapping of MA consumption (mg/kg/18h) from the 40 mg/l concentration and for QTL mapping of 2 mg/kg MA-induced body temperature change.Discussion
Using CRISPR-Cas9 technology, we confirmed that Taar1 is a QTG for MA consumption, with an apparent pleiotropic influence on sensitivity to MA-induced hypothermia. Data from the BXD mice indicate that Taar1 and Oprm1 interact in their effects on these traits. Furthermore, the BXD data confirmed the QTL on mouse chromosome 10 for MA consumption that we originally identified in the MA drinking lines (Belknap et al., 2013), and we mapped a QTL for MA-induced hypothermia to the same region in the BXD strains. Taar1 resides at this location. Overall, the findings for Taar1 are consistent with published results, including those demonstrating that Taar1 knockout mice consume higher amounts of MA and are insensitive to MA-induced hypothermia, compared to their wildtype littermates (Harkness et al., 2015); that Taar1 genotype is strongly correlated with MA consumption and MA-induced hypothermia in multiple genetic models on both homogeneous and heterogeneous backgrounds (Reed et al., 2018); and that TAAR1 agonists reduce MA self-administration (Cotter et al., 2015; Jing and Li, 2015). We conclude that TAAR1 activation has a causal role in protection against high levels of MA intake, and thus, has potential involvement in resistance to MA addiction.
Taar1 and MA consumption
We leveraged CRISPR-Cas9 technology to replace the mutant Taar1m1J allele with the common Taar1+ allele in selectively bred high line mice. The patterns of MA consumption in mice with the different Taar1 genotypes mimicked the patterns previously observed in the MA drinking lines (Harkness et al., 2015; Shabani et al., 2011; Wheeler et al., 2009). Thus, control mice escalated their MA consumption with increasing MA concentration, similar to previous findings in high line mice, whereas MA consumption was low and not concentration-dependent in the knock-in mice, similar to previous findings in the selectively bred low line mice.
Although oral is not the most common route of administration, some MA users do take MA orally, and when taken via this route, MA has a half-life similar to that for other routes of administration (Cruickshank and Dyer, 2009). Furthermore, oral consumption, like other routes of administration, can lead to dependence (Galloway et al., 2010). On average, mice with absent TAAR1 function consume about 6 mg/kg/18h from a 40 mg/l MA solution (Harkness et al., 2015; Hitzemann et al., 2019; Reed et al., 2018; Shabani et al., 2011; Wheeler et al., 2009). High line mice will consume an average of as much as 14 mg/kg/18h, when MA is offered as a 140 mg/l concentration (Shabani et al., 2016). Daily MA use in humans ranges from 300 to 800 mg/day on average (Cho et al., 2001; Simon et al., 2002), which translates to 3.9 to 10.4 mg/kg/day in a 77 kg individual. Whereas the half-life of MA is shorter in rodents compared to humans (Cho et al., 2001), relevant levels of MA are attained in our mouse model, with locomotor stimulant effects observed in the high line mice at consumed doses as low as 0.4 mg/kg in 1h (Shabani et al., 2012a).
In rats with functional TAAR1, TAAR1 agonists delivered systemically or to specific subregions of the mesocorticolimbic system, decreased MA seeking and cocaine seeking (Liu et al., 2017; Pei et al., 2017), suggesting that treatments that increase TAAR1 activity could be beneficial therapeutics. However, when TAAR1 is non-functional, receptor agonists are not an option, leaving us to consider downstream mechanisms of TAAR1 activation. Details of the mechanisms enlisted by TAAR1 activation remain elusive. Expressing TAAR1 using a heterologous expression system to identify its messenger system(s) has proved difficult, potentially due to the predominant intracellular localization of the receptor. Studying the function of TAAR1 is also complicated by the lack of good quality antibodies and selective antagonists (Liu and Li, 2018; Rutigliano et al., 2017). There is one available TAAR1 antagonist, EPPTB (Bradaia et al., 2009), but it has a high rate of clearance (Rutigliano et al., 2017; Stalder et al., 2011) and poor solubility, which limit in vivo use. Our attempts to develop better antagonists with chemist collaborators have not yet met with success, and other investigators in the field have had a similar experience (Lam et al., 2018).
One strategy that could be considered for reducing MA consumption is to target mechanisms that are impacted by the absence of TAAR1 function. For example, TAAR1 modulates monoaminergic neurotransmission (Espinoza et al., 2015a; Leo et al., 2014; Lindemann et al., 2008; Revel et al., 2011; Xie and Miller, 2008), and recently has been implicated in glutamatergic neurotransmission (Espinoza et al., 2018; Espinoza et al., 2015b). High line mice are hyperglutamatergic at baseline, compared to low line mice, in both the nucleus accumbens (NAc) and medial prefrontal cortex (mPFC). They also exhibit a larger glutamate response to MA in the NAc, but not mPFC (Lominac et al., 2016; Szumlinski et al., 2017). In addition, compared to low line mice, high line mice have blunted dopamine levels in both the NAc and mPFC, and exhibit a heightened dopaminergic response to MA in the mPFC, but not NAc (Lominac et al., 2014). Whether these differences are related to Taar1 genotype is unknown, but they suggest certain manipulations to be examined for their role in MA intake.
Genetic basis for MA-induced thermal response
Our previous findings in the MA drinking lines demonstrated that the thermal response to MA is genetically-correlated with level of MA consumption. Thus, high line mice are resistant to the hypothermic effect of MA, whereas low line mice are highly sensitive to this effect of MA (Harkness et al., 2015). Here, we performed QTL analysis using data from the BXD strains and identified a significant QTL on chromosome 10 for the effect of MA on body temperature. This QTL is in the SNP-poor region of chromosome 10 (Shi et al., 2016), where we mapped the QTL for MA consumption in the MA drinking lines and BXD strains. QTL mapping was previously performed in BXD strains for the effects of 4, 8 and 16 mg/kg MA on body temperature recorded 48 min after administration. Although a QTL was mapped on chromosome 10 for the 4 mg/kg dose, it was considerably distal to the current QTL (Grisel et al., 1997), and not likely to be the same one for several reasons. First, the MA doses and time after administration were considerably different from those in our studies. But, more importantly, the Grisel et al. (1997) paper was published well before the Taar1m1J mutation arose in D2 mice (Reed et al., 2018); therefore, all of the BXD strains in that study shared the common Taar1+/+ genotype. We have confirmed the Taar1+/+ genotype of many of the BXD strains that were derived prior to when the mutation appeared (Reed et al., 2018; Shi et al., 2016).
Relationship between MA consumption and body temperature
Sensitivity to MA-induced hypothermia corresponds with low levels of MA consumption in MAHDR-Taar1+/+ knock-in mice, MALDR line mice, the wildtype littermates of Taar1 knockout mice, and BXD strains with the Taar1+/+ genotype. These data could indicate that Taar1 has independent pleiotropic effects on the two traits. Another possibility is that the subjective experience of MA-induced hypothermia in mice with functional TAAR1 limits MA consumption, thereby playing a protective role. Not known is whether MA-induced hypothermia is subjectively unpleasant. However, hypothermia prolonged the associative period during which aversion could be conditioned in a conditioned taste aversion (CTA) procedure (Misanin et al., 1998; Misanin et al., 2002), and mice with functional TAAR1 form MA-induced CTA, whereas those without functional TAAR1 do not (Harkness et al., 2015; Shabani et al., 2012b; Wheeler et al., 2009). Low line and high line mice voluntarily consume similar amounts of MA on the first day that MA is offered, but low line mice then decrease consumption on the subsequent day, and remain at low intake levels from then on (Eastwood et al., 2014; Shabani et al., 2012a). It is possible that initial MA consumption induces hypothermia in low line mice, so that initial MA consumption is associated with this potentially unpleasant physiological effect, and reduces the desire to consume more MA. Additional research is needed to determine if MA drinking produces changes in body temperature, as occurs in response to an IP injection of MA.
Regardless of the chicken-egg question, the Taar1+ allele replacement in high line mice produced MA intake levels like those found in low line mice, indicating that this gene has a major role in determining level of MA intake. Similarly, both the pattern and magnitude of hypothermic response to MA in the knock-in mice resembled those of low line mice. Low line mice had a maximal response of 1–2°C (depending upon replicate line) after treatment with 2 mg/kg MA (Harkness et al., 2015), and the response in the knock-in mice was also maximal at 30 min and of about 1.2°C. Therefore, Taar1 appears also to have a primary role in the hypothermic response to MA. However, Taar1 genotype does not impact basal body temperature, nor the hypothermic response to ethanol (Harkness et al., 2015). Further, the Taar1 mutation does not appear to affect locomotor activity at baseline or in response to the 2 mg/kg dose of MA used here (Shabani et al., 2011; Wheeler et al., 2009), making it unlikely that differences in thermal response to MA in mice with different Taar1 genotypes are due to differences in movement after treatment.
In addition to the significant QTL on chromosome 10 for MA intake and MA-induced thermal response, we identified other genomic regions at the suggestive level of significance. None of these overlapped for the two traits, and the chromosome 4 MA intake QTL that reached the suggestive significance threshold with composite mapping was not identified in our previous studies of the MA drinking lines (Belknap et al., 2013). There may be additional alleles that account for smaller amounts of genetic variance that our analysis did not have the power to detect. The need for large sample sizes to detect rare variants and alleles with smaller independent effects for complex traits is recognized (Belknap, 1998; Belknap et al., 1996; Buchner and Nadeau, 2015; Flint, 2011; Solberg Woods, 2014).
Independent and epistatic effects of Taar1 and Oprm1
Previous findings indicate that Oprm1 regulation and opioid system activity contribute to MA consumption in mice (Belknap et al., 2013; Eastwood et al., 2018; Eastwood and Phillips, 2014a; Eastwood and Phillips, 2014b), and MA dependence/psychosis in humans (Ide et al., 2004; Ide et al., 2006). Our results in the BXD strains with different combinations of Taar1/Oprm1 genotypes support a potential epistatic interaction in the effects of these genes on MA consumption and on the hypothermic response to MA. The combined effects of Oprm1 and Taar1 genotype on each trait was non-additive. Oprm1 genotype impacted MA consumption in Taar1m1J/m1J mice, but not Taar1+/+ mice. The opposite pattern of effect emerged for MA-induced hypothermia; Oprm1 genotype impacted this trait in Taar1+/+ mice, but not Taar1m1J/m1J mice. For both traits, the Oprm1D2/D2 mice exhibited the more robust effects, compared to Oprm1B6/B6 mice, and Oprm1 genotype had an impact only in mice with the stronger MA trait. Thus, Taar1m1J/m1J/Oprm1D2/D2 mice consumed more MA than Taar1m1J/m1J/Oprm1B6/B6 mice, and Taar1+/+/Oprm1D2/D2 mice showed greater MA-induced hypothermia compared to Taar1+/+/Oprm1B6/B6 mice. Low levels or absence of the MA-related phenotype could preclude the Oprm1 genotype from having an effect. This interaction between Taar1 and Oprm1 was supported by composite interval QTL mapping for MA-induced body temperature change, in which controlling for Oprm1 genotype resulted in a significantly increased LOD score, specifically for the chromosome 10 QTL. For MA intake, the LOD score increase did not meet the significance threshold, but the strong statistical trend deserves further consideration, which could be accomplished by testing more BXD strains of appropriate genotypes, as they become available, to increase power.
When we assessed the independent associations of each genotype with these MA-related traits, significant correlations were found only with Taar1 genotype. That Oprm1 genotype alone did not correspond with amount of MA consumed is consistent with our previous findings, indicating that Oprm1 genotype is not a risk factor for MA intake, and its interactive effect with Taar1 may be associated with the regulation of Oprm1 by a significant MA intake risk gene network (Belknap et al., 2013). The interactive effect of Oprm1 with Taar1, in the absence of an independent effect of Oprm1, is consistent with other studies demonstrating epistatic interactions involving polymorphisms that do not have significant independent associations with phenotypes (Lehner, 2011).
Epistasis has been proposed as a reason for variability in disease occurrence among individuals with disease risk mutations. However, like the many intermolecular epistatic interactions with unknown mechanisms (Lehner, 2011), the mechanism(s) underlying the impact of Oprm1 genotype on these Taar1-associated traits is not yet known. One possible mechanism to explore relates to differential Oprm1 expression. A polymorphism in the promoter region of Oprm1 was related to promoter activity, in which D2 mice had greater Oprm1 promoter activity compared to B6 mice (Doyle et al., 2006). It was hypothesized that this might alter MOP-r expression. Doyle et al. (2014) found greater Oprm1 expression in the frontal cortex of D2, compared to B6 mice. It is possible that greater expression of Oprm1 could result in greater MA intake in Taar1m1J/m1J/Oprm1D2/D2 mice, and greater MA-induced hypothermia in Taar1+/+/Oprm1D2/D2 mice. However, Eastwood et al. (2018) examined MOP-r protein density in D2 vs. B6 mice in the prefrontal cortex and did not find a difference. On the other hand, our low line mice, which consumes little MA, is sensitive to MA-induced hypothermia and is largely of the Taar1+/+ genotype, had both greater Oprm1 expression and MOP-r density exclusively in the prefrontal cortex, compared to our high line mice, which has the opposite MA-related phenotypes and Taar1m1J/m1J genotype (Belknap et al., 2013; Eastwood et al., 2018). Thus, the hypothesis that greater Oprm1 expression results in greater MA-induced hypothermia in Taar1+/+ genotype mice is consistent with these findings, but the hypothesis that greater Oprm1 expression results in greater MA intake in Taar1m1J/m1J genotype mice is not consistent with the existing data. In fact, greater MA intake is associated with reduced prefrontal cortex Oprm1 and MOP-r expression in the high line, and peripheral application of drugs that increase MOP-r activity reduced MA consumption in high line mice (Eastwood et al., 2018; Eastwood and Phillips, 2014a). We are currently creating Taar1+/+ knock-in mice on a D2 strain background and Taar1m1J/m1J knock-in mice on a B6 background. Epistatic behavioral and molecular effects may be more approachable for study on these homogeneous backgrounds.
Effects of Taar1 and Oprm1 on total fluid intake
We observed genotype-dependent effects on total volume consumed from the MA and water tubes, but conclude that genotypic differences in MA consumption were not an artifact of overall fluid intake differences. The total volume consumption difference between Taar1m1J/m1J and Taar1+/+ mice was small, compared to their robust difference for MA intake. Furthermore, the interactive effect of Taar1 and Oprm1 was specific to MA intake. We previously reported greater total fluid intake in high line mice, compared to low line mice, of about the same magnitude (~0.5 ml) in some, but not all, replicate sets of lines (Hitzemann et al., 2019; Shabani et al., 2011; Shabani et al., 2019; Wheeler et al., 2009). Locomotor behavior has been found to be positively associated with MA intake in the low and high lines (Shabani et al., 2012a), and larger intake volumes in mice consuming more MA could be related to stimulant effects that impact fluid needs. But, because a difference in total fluid intake has not always been found, another interpretation is that, due to greater reinforcing effects of MA in mice of the Taar1m1J/m1J genotype, such as the high line mice (Shabani et al., 2012a), avidity for MA is greater, leading to larger intake volumes.
Body temperature differences unrelated to MA exposure
We observed some differences in body temperature between the knock-in and control mice that were unrelated to MA treatment. However, these baseline differences did not appear to account for the genotype-dependent MA response. Although the knock-in mice had higher baseline temperatures and therefore, a greater opportunity to experience a reduction in body temperature when treated with MA, the difference at baseline was of about 0.5°C, whereas the difference at 30 min post MA treatment was over 1°C. This was, in part, because the MA responses were qualitatively opposite, with the knock-in mice exhibiting decreases in body temperature, similar to MALDR mice, and the controls exhibiting a time-dependent increase, similar to MAHDR mice (Harkness et al., 2015).
Sex differences
In our previous studies, sex differences for MA consumption and the thermal response to MA were inconsistently found, and some are summarized in Reed et al. (2018). Generally, when there have been differences, females have consumed about 15% or 0.5 mg/kg more MA than males, but in no case were there differences between the MAHDR and MALDR lines or between mice with the different Taar1 genotypes in only one sex (Eastwood and Phillips, 2014a; Harkness et al., 2015; Reed et al., 2018; Shabani et al., 2011; Shabani et al., 2016; Wheeler et al., 2009). Likewise, sex did not have a significant impact on genotype-dependent differences in MA consumption in the present studies. Similarly, there were some sex effects in the body temperature studies, but sex did not determine the pattern of response to saline vs. MA in each genotype across time. Thus, these sex effects did not impact overall interpretations of data for either trait.
Potential limitations
One consideration regarding Taar1 involvement in MA consumption and MA thermal response is the potential for differential potency of agonists interacting with the two TAAR1 isoforms. However, our analyses indicate that the receptor expressed by Taar1m1J is non-functional, rather than sub-functional. Thus, the D2-like isoform of TAAR1, expressed by Taar1m1J, exhibits no cAMP response to a wide range of MA concentrations, or to TAAR1-specific agonists (Harkness et al., 2015; Shi et al., 2016). If the Taar1 mutation reduces potency of agonists, we would expect a shift in the agonist dose-response curve for cAMP accumulation. Further, at doses of MA up to 16 mg/kg, high line mice fail to exhibit MA-induced hypothermia, whereas low line mice exhibit hypothermia at doses as low as 1 mg/kg that reverses toward hyperthermia at higher doses (Harkness et al., 2015). Thus, the hypothermic effect of MA via TAAR1 may compete with the expression of MA-induced hyperthermia occurring via a non-TAAR1 mechanism in low mice.
The gRNA used to produce the MAHDR-Taar1+/+ knock-in mice was a perfect match for the sequence on either side of the chromosome 10 Taar1 SNP, with imperfect mapping to some additional chromosomes. In our original QTL analysis for MA drinking and in the current analysis (Figure 7), there were no significant or suggestive QTLs mapped to those additional chromosomes (Belknap et al., 2013; and data herein). This suggests that if off-target deletions or insertions occurred at these locations, they would be unlikely to have the profound impact on the MA-related traits that we have observed. The chromosome 10 QTL accounts for 60% of the genetic variance in MA intake, and Taar1 genotype-phenotype correlations range from r = 0.81 to 0.96 for 40 mg/L MA intake, and r = 0.71 to 0.82 for 2 mg/kg MA-induced hypothermia (Reed et al., 2018; and data herein). Therefore, the large effects we observed would be expected for replacement of Taar1m1J with Taar1+. The high line mice possess a genetically heterogeneous background; thus, sequencing the CRISPR knock-in and control mice would not provide specific information about off-target effects, because the existing genetic variation would not be separable from such effects. Further, there appears to be some consensus that off-target effects of CRISPR-Cas9 are rare and indistinguishable from the background rate of de novo mutation (e.g., Anderson et al., 2018; Ayabe et al., 2019; Iyer et al., 2018; Mianné et al., 2016; Nakajima et al., 2016; Willi et al., 2018).
The knock-in and control mice used in our studies were produced by independent breeding pairs. As noted in Materials and methods, this interbreeding of like-Taar1-genotype individuals is a consistent feature for all other existing populations in which the Taar1 SNP exists. The only exception has been for the B6 x D2 F2 mice that were produced to create the high and low selected lines, which existed as multi-Taar1-genotype littermates. An equally strong association of Taar1 genotype with MA intake was found in these mice (Reed et al., 2018).
Conclusions
In summary, these data demonstrate that Taar1 is a major contributor to MA intake and MA-induced hypothermia, and provide evidence for an interaction between Taar1 and Oprm1 in their effects on these MA traits. Genetic variation in the human Oprm1 gene has been associated with MA dependence/psychosis (Ide et al., 2004; Ide et al., 2006), but additional research results in this area have not appeared in the literature. The potential impact of human TAAR1 genetic variation on risk for MA use or on the magnitude of MA-related phenotypes is not yet known. There are over 200 non-synonymous SNPs with the potential for missense variants in human TAAR1 (Rutigliano et al., 2017) Initial examination of some of these human TAAR1 variants indicates that TAAR1 proteins with variable levels of function are expressed (Shi et al., 2016). The success of TAAR1 agonists in treatment is dependent upon at least partial receptor function that can be enhanced. Whether human TAAR1 variants are relevant to MA use and addiction in humans, and whether TAAR1 agonists are a viable treatment option, are important research questions to pursue. Further, research into potential interactive effects between human TAAR1 and OPRM1 variants for their impact on MA use risk and response would be valuable.
Materials and methods
Key resources table
| Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
|---|---|---|---|---|
| Strain, strain background (Mus musculus, females) | CD1/NCrl | Charles River | Strain code: 022 | foster dams for generation of knock-in mice |
| Strain, strain background (Mus musculus, females and males) | BXD recombinant inbred strains | University of Tennessee Health Sciences Center | original breeders from Dr. Robert Williams; offspring for research produced by Dr. Tamara Phillips at VA Portland Health Care System (VAPORHCS) | |
| Genetic reagent (Mus musculus, females) | MAHDR | VAPORHCS | methamphetamine high drinking mice; created by Dr. Tamara Phillips; all are homozygous for SNP rs33645709 | |
| Genetic reagent (Mus musculus, females and males) | MAHDR-Taar1m1J/m1J control | this paper | created by Oregon Health & Science University Transgenic Mouse Models Shared Resource; offspring for research produced by Dr. Tamara Phillips at VAPORHCS | |
| Genetic reagent (Mus musculus, females and males) | MAHDR-Taar1+/+ knock-in | this paper | created by Oregon Health & Science University Transgenic Mouse Models Shared Resource; offspring for research produced by Dr. Tamara Phillips at VAPORHCS | |
| Biological sample (Mus musculus) | tail snip | other | obtained from breeders and research animals to determine genotype | |
| Sequence-based reagent | Taar1 gRNA | this paper | designed by ThermoFisher Scientific; reagent for generation of knock-in mice (Mus musculus) | |
| Sequence-based reagent | 100b oligo single strand DNA for Taar1+ (ossDNA) | this paper | MGI:MGI2148258 | synthesized by ThermoFisher Scientific; reagent for generation of knock-in mice (Mus musculus) |
| Sequence-based reagent | Cas9 mRNA | TriLink | Catalog number: L-7606 | reagent for generation of knock-in mice (Mus musculus) |
| Sequence-based reagent | primers for genotyping Oprm1 | this paper | original reference for these primers: Ferraro et al. (2005); based on Oprm1 sequence NM_001304955 (Mus musculus) | |
| Sequence-based reagent | DdeI restriction enzyme | ThermoFisher Scientific | Catalog number: ER1882 | reagent for Oprm1 genotyping (Helicobacter pylori, RFL3) |
| Commercial assay or kit | QuickExtract DNA Extraction Solution | Lucigen | Catalog number: QE09050 | |
| Commercial assay or kit | Hotstart DNA Polymerase Kit | Qiagen | Catalog number: 203205 | |
| Commercial assay or kit | Taqman kit for determining Taar1 alleles | ThermoFisher Scientific | Custom order | probes based on Taar1 sequence: NM_053205.1; SNP rs33645709 |
| Chemical compound, drug | ethidium bromide | Sigma Aldrich | Catalog number: E1510 | reagent for Oprm1 genotyping |
| Chemical compound, drug | (+)-methamphetamine hydrochloride | Sigma Aldrich | Catalog number: M8750 | |
| Chemical compound, drug | sterile 0.9% saline | Baxter Healthcare Corporation | Catalog number: 2F7124 | vehicle for methamphetamine |
| Software, algorithm | Statistica | TIBCO Software Inc | ||
| Software, algorithm | GeneNetwork | University of Tennessee | RRID:SCR_002388 | www.genenetwork.org |
| Software, algorithm | R | The R Foundation for Statistical Computing | RRID:SCR_001905 | www.r-project.org |
Animals
All experiments were performed in accordance with the National Institutes of Health Guidelines for the Care and Use of Laboratory Animals (National Research Council, 2011) and were approved by the Institutional Animal Care and Use Committee of the VA Portland Health Care System (VAPORHCS). Male and female mice were group housed (2–5 per cage) in filtered polycarbonate shoebox cages (28.5 × 17.5×12 cm) that were lined with Bed-o’Cobs bedding (The Andersons, Inc, Maumee, OH) and fitted with wire tops, except during MA drinking studies. For these studies, mice were individually housed in the same type of caging and provided with a cotton fiber nestlet for enrichment. All mice were maintained on a 12:12h light:dark schedule (lights on at 0600h), and had ad libitum access to laboratory rodent block food (PicoLab Laboratory Rodent Diet 5LOD, 4.5% fat content; Animal Specialties, Woodburn, OR) and tap water.
Design of the CRISPR-Cas9 knock-in of Taar1m1J and CRISPR reagents
The Taar1m1J mouse SNP, rs33645709, encodes a non-synonymous proline to threonine mutation at amino acid position 77 that originally occurred in D2 mice in 2001–2003 as a spontaneous mutation (Reed et al., 2018). Proline is often highly conserved, because it is essential for conservation of a particular protein fold. In fact, this mutation renders TAAR1 non-functional, and the SNP is fixed (homozygous) in high line mice (Harkness et al., 2015; Reed et al., 2018; Shi et al., 2016). CRISPR-Cas9 technology (Jinek et al., 2012; Zhang, 2012) was used to replace the mutant Taar1m1J allele with the reference Taar1+ allele in high line mice to generate a homozygous Taar1+/+ knock-in on the high MA drinking line background at the Oregon Health & Science University Transgenic Mouse Models Shared Resource. CRISPR in vitro transcribed guide RNAs (gRNAs), targeting the specified region of Taar1, and donor (100b) oligo single strand DNA (ossDNA) for incorporation of the Taar1+ sequence were designed and synthesized by ThermoFisher Scientific (Carlsbad, CA, USA) (Figure 8a). During the development of the gRNA, in silico analysis with the basic local alignment search tool or BLAST was used to assure specificity of the sequence. The gRNA chosen for use was a perfect match for Taar1, with no mismatches, except at the SNP location. The gRNA sequence had similarity to nine additional regions on chromosomes 1, 3, 5, 9, 11 and 19, but with 2–3 mismatches in each case. This gRNA had the sequence tctgataatgAcctgcagcatgg. The location of the SNP is indicated by the capital A, encoding the mutant genotype present in the high line. The location of the SNP within the Taar1 ossDNA sequence, which is the reference genotype, is indicated below by the capital C, with the gRNA sequence underlined. The gene editing replaces A with C: cFFctggctccttcactccatggccattgtcgactttctgctgggctgtctgataatgCcctgcagcatggtgagaactgttgagcgctgttggtatZZt. The ossDNA was designed without a silent mutation and to be asymmetric to provide increased efficiency of the gRNA. Cas9 mRNA (catalog number L-7606) was purchased from TriLink (San Diego, CA, USA).
(a) Strategy for insertion of a point mutation into mouse Taar1 exon 1. The guide sequence is indicated in green text and the protospacer adjacent motif (PAM) is indicated by bold black text. The targeting template is indicated in black text, except for the single base pair to be inserted during the CRISPR-Cas9 process, which is indicated in red text. (b) A sequencing chromatograph of a homozygous Taar1m1J/m1J mouse that is representative of all MAHDR mice. The red box indicates the single nucleotide targeted for replacement during the CRISPR-Cas9 process. (c) A sequencing chromatograph of a homozygous MAHDR-Taar1+/+ edited knock-in mouse. The nucleotide that was successfully inserted during the CRISPR-Cas9 procedure is indicated by a red square. (d) An agarose genotyping gel example of Oprm1 genotypes for Oprm1B6/B6, Oprm1D2/D2 and Oprm1B6/D2 (from a B6 x D2 cross). All BXD mice were either homozygous Oprm1B6/B6 or Oprm1D2/D2, as well as homozygous for one or the other Taar1 allele. BXD, C57BL/6J (B6) x DBA2/J (D2); Oprm1, mu-opioid receptor gene; Oprm1B6/B6, homozygous for B6 Oprm1 allele; Oprm1D2/D2, homozygous for D2 Oprm1 allele; Oprm1B6/D2, heterozygous for B6 and D2 Oprm1 alleles; PAM, protospacer adjacent motif; Taar1, trace amine-associated receptor 1 gene; Taar1+/+, homozygous for reference Taar1+ allele; Taar1m1J/m1J, homozygous for mutant Taar1m1J allele.
Generation of MAHDR-Taar1+/+ knock-in mice
High line mice were embryo donors. Three-four week old, female mice were super-ovulated, mated with high line males, and 0.5 day old eggs were isolated as previously described (Hogan et al., 1994). A mixture of Cas9 mRNA, Taar1 gRNA, and donor ossDNA for Taar1+ (at concentrations of 100 ng/µl, 50 ng/µl, and 100 ng/µl, respectively) was injected into the pronuclei of one-cell fertilized eggs. The injected eggs were transferred into the oviducts of pseudo pregnant recipient CD1/NCrl foster dams (Charles River, San Diego, CA, USA), and the surviving offspring were sequenced (see Genotyping and sequencing) to determine if they had retained the Taar1m1J allele or if an alteration resulting in insertion of Taar1+ had occurred. Control mice (MAHDR-Taar1m1J/m1J) were derived at the same time, from those mice in which the Taar1m1J allele was not successfully altered. These offspring were transported to the VAPORHCS vivarium and bred to produce the MAHDR-Taar1+/+ knock-in and MAHDR-Taar1m1J/m1J control mice that participated in our experiments. Knock-in and control mice were maintained as within-line breeding pairs in a shared colony room. This breeding of individuals with the same Taar1 genotype is consistent with the maintenance of all other populations in which the Taar1 SNP exists (e.g., MAHDR vs. MALDR; BXD strains; D2 mice from The Jackson Laboratory vs. D2 mice from other suppliers; Reed et al., 2018). Representative chromatographs for the Taar1m1J (control) and Taar1+ (knock-in) sequences are displayed in Figure 8b and c, respectively.
Separate cohorts of knock-in and control mice were tested for two-bottle choice MA intake or MA-induced change in body temperature, as detailed below and in our previous publications (Harkness et al., 2015; Hitzemann et al., 2019; Shabani et al., 2011; Wheeler et al., 2009). For MA intake, 20 knock-in and 20 control mice were tested in a single cohort (n = 10/genotype/sex), and the mice were 87–90 days old at the start of testing. For body temperature, 48 knock-in and 50 control mice were tested in 3 cohorts of 25–46 mice (final n = 12–13/genotype/dose/sex), and the mice were 79–94 days old on the day of testing.
BXD mice
Breeding pairs that produced the BXD mice tested in these studies were obtained from RWW (University of Tennessee Health Science Center, Memphis, TN), and established within the VAPORHCS vivarium. Specific strains were chosen according to Taar1 and Oprm1 genotypes to allow identification of potential independent and interactive effects on MA-related phenotypes. Separate cohorts of offspring were tested for two-bottle choice MA intake or MA-induced body temperature change. Mice were homozygous for either the reference Taar1+ or mutant Taar1m1J allele, as well as for the B6 or D2 Oprm1 allele, in all possible combinations (Table 1). A total of 233 BXD mice (134 female and 99 male), ranging from 53 to 114 days of age, were tested for MA intake in 4 cohorts of 27–97 mice. A total of 325 BXD mice (171 female and 154 male), ranging from 49 to 124 days of age, were tested for MA-induced body temperature change in 16 cohorts of 8–36 mice. Numbers of BXD mice with the four potential combined Taar1/Oprm1 genotypes (Table 1) were n = 22–45/genotype/sex for MA intake, and n = 14–25/genotype/sex/MA dose for the body temperature study.
Genotyping and sequencing
We previously established and described genotyping methods for the Taar1 SNP (Harkness et al., 2015; Reed et al., 2018). Genomic DNA was extracted from ear punch or tail snip samples, using QuickExtract DNA Extraction Solution (Lucigen, Middleton, WI, USA). For samples from the BXD mice and the original mice created by the CRISPR-Cas9 procedure, the Taar1 region was PCR amplified using a Hotstart DNA polymerase kit (Qiagen, Valencia, CA, USA), with sequence-specific primers surrounding the region of interest (see Harkness et al., 2015 for details). PCR products were sequenced at the Oregon Health & Science University Sequencing Core and aligned with the mouse Taar1 sequence (NM_053205.1) to identify which allele(s) was present based on the rs33645709 SNP. All knock-in and control offspring of the original mice were genotyped with an updated protocol that uses a standard Taqman (ThermoFisher Scientific) with fluorescent probes to detect and differentiate the Taar1 alleles (Reed et al., 2018).
We performed Oprm1 genotyping for BXD mice as previously described (Ferraro et al., 2005). Sequence-specific primers for exon 10 of Oprm1 were used as the forward primer (5'-ggttatgcctctctggattag-3') and reverse primer (5'-tccatcgcttacatcttacca-3'). A SNP in exon 10 of the B6 Oprm1 gene (Oprm1B6) creates a DdeI restriction site (Ferraro et al., 2005), so that two bands (235 and 168 bp) are created when amplified DNA from mice with Oprm1B6 is digested with the DdeI restriction enzyme; DdeI does not excise the D2 Oprm1 gene (Oprm1D2), resulting in one band. After amplification, PCR products were digested with DdeI, run on a 4% agarose gel, and visualized using ethidium bromide staining. A representative gel is displayed in Figure 8d to indicate differentiation of homozygous (B6/B6, D2/D2) and heterozygous (B6/D2) Oprm1 genotypes. However, all BXD mice used in these studies were homozygous.
Drugs and reagents
(+)MA hydrochloride was purchased from Sigma (St. Louis, MO, USA) and dissolved in tap water for drinking or in sterile 0.9% saline (Baxter Healthcare Corp., Deerfield, IL, USA) for injection. For the body temperature studies, saline and MA were administered by intraperitoneal injection at a volume of 10 ml/kg body weight.
Two-bottle choice MA drinking
Methods were consistent with those we used to measure two-bottle choice MA intake during the production of the MA drinking lines (Hitzemann et al., 2019; Shabani et al., 2011; Wheeler et al., 2009). During the initial two days of the study, mice had access to two tubes filled with tap water to familiarize them with the novel drinking apparatus. Starting on the third day, one water tube was replaced with a tube containing MA dissolved in tap water to which mice had access for 18h/day, beginning 3h before the lights turned off. The MA tube was removed for the remaining 6h of each day. Unpublished data (Phillips) indicate that the 6h forced abstinence periods between periods of 18h MA access are important for higher MA intake levels, compared to 24h MA access periods. The MA concentration was 20 mg/l for 4 days, and was increased to 40 mg/l for an additional 4 days. The MA and water tube positions were alternated every other day to account for potential side bias. MA consumption was indexed as the average intake on the second and fourth days of access to each MA concentration, as these days represent the day after the tube positions were switched, when mice should be familiar with the relative locations of the water and MA tubes. MA consumption was measured in ml (accuracy = 0.2 ml), and then converted to mg/kg, based on body weight measured every two days. Total volume consumed (ml) from the water and MA tubes during the 18h MA access periods was also measured.
MA-induced change in body temperature
Methods for determining MA-induced change in body temperature were consistent with those we established previously (Harkness et al., 2015). Mice were tested after injection of saline or 2 mg/kg MA at an ambient temperature of 21 ± 2°C, for 120 min beginning at 0800h, two hours after lights on. The MA dose was chosen from previous results demonstrating that mice with functional TAAR1 exhibit a robust hypothermic response to 2 mg/kg MA 30 min after administration that is absent in mice with non-functional TAAR1 (Harkness et al., 2015; Reed et al., 2018). Treatment groups (saline or MA) were assigned so that males and females were equally represented, and mice from each family were distributed equally between the groups. Mice were weighed and then placed into individual perforated acrylic plastic cubicles (8 × 19×8 cm in W x H x L) that served to prevent huddling-associated effects on body temperature (Crabbe et al., 1987; Crabbe et al., 1989). Mice were undisturbed for 1h before a baseline rectal temperature was obtained at time 0 (T0) by inserting a glycerin-coated, 5 mm probe attached to a Thermalert model TH-8 digital thermometer (Sensortek, Clifton, NJ, USA) 2.5 cm into the rectum for 5 s. Saline or MA was then administered, and mice were returned to their cubicles. Rectal temperature was again measured at 30, 60, 90, and 120 min post-injection (T30-T120).
QTL analysis
QTL analyses using BXD strain means were conducted using the WebQTL mapping module of GeneNetwork (www.genenetwork.org; RRID:SCR_002388). The traits mapped were: 1) MA intake when the 40 mg/l MA concentration was offered, which is the phenotype used for selective breeding (Shabani et al., 2011; Wheeler et al., 2009), and 2) MA-induced body temperature change at 30 min post-treatment. The 30 min time point is when the hypothermic effect of the 2 mg/kg dose of MA is largest (Harkness et al., 2015; Reed et al., 2018). The change score was calculated by subtracting temperatures at 30 min from baseline temperatures taken immediately before treatment. QTL mapping was initially performed using the interval mapping function and then composite interval mapping was applied, with Oprm1 genotype at rs29351111 as a co-factor, to assess the potential interaction of Oprm1 and Taar1. QTL significance was assessed using the LOD score obtained after 1000 or 2000 permutations, for interval or composite interval mapping, respectively, and was considered significant if the genome-wide p-value was <0.05, and considered suggestive if the genome-wide p-value was <0.63, which yields one false positive per genome scan on average. These are the standard settings used for GeneNetwork QTL mapping. R version 3.6.0 (The R Foundation for Statistical Computing, https://www.r-project.org/foundation/; RRID:SCR_001905) was used to analyze changes in LOD scores computed from interval vs. composite interval QTL mapping with a Chi-squared test comparing additive vs. interactive models for the effects of Taar1 and Oprm1 genotype. MA consumption and MA-induced body temperature change means for the BXD strains are available in GeneNetwork. For additional information about using GeneNetwork for QTL mapping, see Mulligan et al. (2017).
Data analysis
Statistica 13.3 (TIBCO Software, Inc, Palo Alto, CA, USA) was used for statistical analyses other than QTL mapping. For MA drinking studies, MA consumption (mg/kg) and total volume consumed (ml) data were analyzed by repeated measures ANOVA, with MA concentration as the repeated measure, and genotype (or strain) and sex as possible independent variables. Body temperature data were analyzed by repeated measures ANOVA, with time as the repeated measure, and genotype (or strain), sex, and MA dose as possible independent variables. For BXD MA drinking and body temperature data, initial analyses characterizing the strains did not include sex as a factor, due to group sizes that were too small. However, analyses examining associations between Taar1/Oprm1 genotypes and MA-related phenotypes did examine the potential impact of sex. For these analyses, to examine potential interaction effects on the measured phenotypes, data for the entire BXD population were analyzed with sex, Taar1 genotype, and Oprm1 genotype coded as separate independent variables (Reed et al., 2018). Effects were considered statistically significant at p<0.05. Complex interactions were simplified by ANOVAs at levels of a particular factor. Two-way interactions were further interpreted using simple main effects analysis, with Bonferroni correction. To compare body temperature data at post-injection time points (T30-120) to T0 within a particular genotype, Dunnett’s post hoc test was used. Pearson’s r was used to calculate correlations between phenotypes and Taar1 or Oprm1 genotype. Sample sizes for these studies were based on our previous MA drinking and body temperature studies (e.g., Harkness et al., 2015; Reed et al., 2018).
Acknowledgements
We thank Zhen Zhu and Nicholas Varra for assistance with data collection, and Jason Erk for assistance with data collection and mouse colony maintenance. This study was supported by NIH NIDA grants R01DA046081 (TP), P50DA018165 (TP, AJ), U01DA041579 (TP), and P30DA044223 (RW), Department of Veterans Affairs I01B × 002106 (TP), I01B × 002758 (AJ), and I01B × 003279 (KN), Oregon Health & Science University-University Shared Resources Pilot Funding to the Transgenic Mouse Models Shared Resource (KN), University of Tennessee Center for Integrative and Translational Science (RW), and the VA Research Career Scientist program (TP, AJ). The contents do not represent the views of the US Department of Veterans Affairs or the United States Government.
Funding Statement
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. The contents do not represent the views of the U.S. Department of Veterans Affairs or the United States Government.
Contributor Information
Catherine Dulac, Harvard University, United States.
Jonathan Flint, University of California, Los Angeles, United States.
Funding Information
This paper was supported by the following grants:
National Institute on Drug Abuse R01 DA046081 to Tamara J Phillips.
National Institute on Drug Abuse P50 DA018165 to Aaron J Janowsky, Tamara J Phillips.
National Institute on Drug Abuse U01 DA041579 to Tamara J Phillips.
National Institute on Drug Abuse P30 DA044223 to Robert W Williams.
Department of Veterans Affairs I01BX002106 to Tamara J Phillips.
Department of Veterans Affairs I01BX002758 to Aaron J Janowsky.
Department of Veterans Affairs I01BX003279 to Kim A Neve.
Oregon Health & Science University University Shared Resource award to Kim A Neve.
U.S. Department of Veterans Affairs Research Career Scientist Program - Career Scientist award to Aaron J Janowsky, Tamara J Phillips.
University of Tennessee Center for Integrative and Translational Science Center support to Robert W Williams.
Additional information
Competing interests
No competing interests declared.
Author contributions
Conceptualization, Data curation, Formal analysis, Investigation, Writing—original draft, Writing—review and editing.
Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing—original draft, Project administration, Writing—review and editing.
Data curation, Validation, Investigation, Visualization, Methodology.
Conceptualization, Validation, Investigation, Visualization, Methodology, Writing—review and editing.
Data curation, Formal analysis, Investigation, Visualization.
Resources, Data curation, Formal analysis, Funding acquisition, Validation, Visualization, Methodology, Writing—review and editing.
Conceptualization, Funding acquisition, Methodology, Project administration, Writing—review and editing.
Conceptualization, Resources, Supervision, Visualization, Methodology, Project administration, Writing—review and editing.
Conceptualization, Funding acquisition, Writing—review and editing.
Conceptualization, Resources, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing.
Ethics
Animal experimentation: All experiments were performed in accordance with the National Institutes of Health Guidelines for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee of the VA Portland Health Care System (VAPORHCS), protocol numbers 3169-14, 3169-16 and 3140-17.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 2 through 7.
References
- Anderson KR, Haeussler M, Watanabe C, Janakiraman V, Lund J, Modrusan Z, Stinson J, Bei Q, Buechler A, Yu C, Thamminana SR, Tam L, Sowick MA, Alcantar T, O'Neil N, Li J, Ta L, Lima L, Roose-Girma M, Rairdan X, Durinck S, Warming S. CRISPR off-target analysis in genetically engineered rats and mice. Nature Methods. 2018;15:512–514. 10.1038/s41592-018-0011-5. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Ayabe S, Nakashima K, Yoshiki A. Off- and on-target effects of genome editing in mouse embryos. Journal of Reproduction and Development. 2019;65:1–5. 10.1262/jrd.2018-128. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Belknap JK, Mitchell SR, O'Toole LA, Helms ML, Crabbe JC. Type I and type II error rates for quantitative trait loci (QTL) mapping studies using recombinant inbred mouse strains. Behavior Genetics. 1996;26:149–160. 10.1007/BF02359892. [Abstract] [CrossRef] [Google Scholar]
- Belknap JK. Effect of within-strain sample size on QTL detection and mapping using recombinant inbred mouse strains. Behavior Genetics. 1998;28:29–38. 10.1023/A:1021404714631. [Abstract] [CrossRef] [Google Scholar]
- Belknap JK, McWeeney S, Reed C, Burkhart-Kasch S, McKinnon CS, Li N, Baba H, Scibelli AC, Hitzemann R, Phillips TJ. Genetic factors involved in risk for methamphetamine intake and sensitization. Mammalian Genome. 2013;24:446–458. 10.1007/s00335-013-9484-9. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Bradaia A, Trube G, Stalder H, Norcross RD, Ozmen L, Wettstein JG, Pinard A, Buchy D, Gassmann M, Hoener MC, Bettler B. The selective antagonist EPPTB reveals TAAR1-mediated regulatory mechanisms in dopaminergic neurons of the mesolimbic system. PNAS. 2009;106:20081–20086. 10.1073/pnas.0906522106. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Buchner DA, Nadeau JH. Contrasting genetic architectures in different mouse reference populations used for studying complex traits. Genome Research. 2015;25:775–791. 10.1101/gr.187450.114. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Bunzow JR, Sonders MS, Arttamangkul S, Harrison LM, Zhang G, Quigley DI, Darland T, Suchland KL, Pasumamula S, Kennedy JL, Olson SB, Magenis RE, Amara SG, Grandy DK. Amphetamine, 3,4-methylenedioxymethamphetamine, lysergic acid diethylamide, and metabolites of the catecholamine neurotransmitters are agonists of a rat trace amine receptor. Molecular Pharmacology. 2001;60:1181–1188. 10.1124/mol.60.6.1181. [Abstract] [CrossRef] [Google Scholar]
- Cho AK, Melega WP, Kuczenski R, Segal DS. Relevance of pharmacokinetic parameters in animal models of methamphetamine abuse. Synapse. 2001;39:161–166. 10.1002/1098-2396(200102)39:2<161::AID-SYN7>3.0.CO;2-E. [Abstract] [CrossRef] [3.0.CO;2-E&" target="pmc_ext" ref="reftype=other&article-id=5632511&issue-id=254822&journal-id=1950&FROM=Article%7CCitationRef&TO=Content%20Provider%7CLink%7CGoogle%20Scholar">Google Scholar]
- Chomchai C, Chomchai S. Global patterns of methamphetamine use. Current Opinion in Psychiatry. 2015;28:269–274. 10.1097/YCO.0000000000000168. [Abstract] [CrossRef] [Google Scholar]
- Cotter R, Pei Y, Mus L, Harmeier A, Gainetdinov RR, Hoener MC, Canales JJ. The trace amine-associated receptor 1 modulates methamphetamine's neurochemical and behavioral effects. Frontiers in Neuroscience. 2015;9:39. 10.3389/fnins.2015.00039. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Crabbe JC, Kosobud A, Tam BR, Young ER, Deutsch CM. Genetic selection of mouse lines sensitive (cold) and resistant (hot) to acute ethanol hypothermia. Alcohol and Drug Research. 1987;7:163–174. [Abstract] [Google Scholar]
- Crabbe JC, Feller DJ, Dorow JS. Sensitivity and tolerance to ethanol-induced hypothermia in genetically selected mice. The Journal of Pharmacology and Experimental Therapeutics. 1989;249:456–461. [Abstract] [Google Scholar]
- Cruickshank CC, Dyer KR. A review of the clinical pharmacology of methamphetamine. Addiction. 2009;104:1085–1099. 10.1111/j.1360-0443.2009.02564.x. [Abstract] [CrossRef] [Google Scholar]
- Doyle GA, Sheng XR, Schwebel CL, Ferraro TN, Berrettini WH, Buono RJ. Identification and functional significance of polymorphisms in the mu-opioid receptor gene (oprm) promoter of C57BL/6 and DBA/2 mice. Neuroscience Research. 2006;55:244–254. 10.1016/j.neures.2006.03.005. [Abstract] [CrossRef] [Google Scholar]
- Doyle GA, Schwebel CL, Ruiz SE, Chou AD, Lai AT, Wang MJ, Smith GG, Buono RJ, Berrettini WH, Ferraro TN. Analysis of candidate genes for morphine preference quantitative trait locus Mop2. Neuroscience. 2014;277:403–416. 10.1016/j.neuroscience.2014.07.020. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Eastwood EC, Barkley-Levenson AM, Phillips TJ. Methamphetamine drinking microstructure in mice bred to drink high or low amounts of methamphetamine. Behavioural Brain Research. 2014;272:111–120. 10.1016/j.bbr.2014.06.035. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Eastwood EC, Eshleman AJ, Janowsky A, Phillips TJ. Verification of a genetic locus for methamphetamine intake and the impact of morphine. Mammalian Genome. 2018;29:260–272. 10.1007/s00335-017-9724-5. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Eastwood EC, Phillips TJ. Morphine intake and the effects of naltrexone and buprenorphine on the acquisition of methamphetamine intake. Genes, Brain and Behavior. 2014a;13:226–235. 10.1111/gbb.12100. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Eastwood EC, Phillips TJ. Opioid sensitivity in mice selectively bred to consume or not consume methamphetamine. Addiction Biology. 2014b;19:370–379. 10.1111/adb.12003. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Espinoza S, Ghisi V, Emanuele M, Leo D, Sukhanov I, Sotnikova TD, Chieregatti E, Gainetdinov RR. Postsynaptic D2 dopamine receptor supersensitivity in the striatum of mice lacking TAAR1. Neuropharmacology. 2015a;93:308–313. 10.1016/j.neuropharm.2015.02.010. [Abstract] [CrossRef] [Google Scholar]
- Espinoza S, Lignani G, Caffino L, Maggi S, Sukhanov I, Leo D, Mus L, Emanuele M, Ronzitti G, Harmeier A, Medrihan L, Sotnikova TD, Chieregatti E, Hoener MC, Benfenati F, Tucci V, Fumagalli F, Gainetdinov RR. TAAR1 modulates cortical glutamate NMDA receptor function. Neuropsychopharmacology. 2015b;40:2217–2227. 10.1038/npp.2015.65. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Espinoza S, Leo D, Sotnikova TD, Shahid M, Kääriäinen TM, Gainetdinov RR. Biochemical and functional characterization of the trace Amine-Associated receptor 1 (TAAR1) Agonist RO5263397. Frontiers in Pharmacology. 2018;9:645. 10.3389/fphar.2018.00645. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Ferraro TN, Golden GT, Smith GG, Martin JF, Schwebel CL, Doyle GA, Buono RJ, Berrettini WH. Confirmation of a major QTL influencing oral morphine intake in C57 and DBA mice using reciprocal congenic strains. Neuropsychopharmacology. 2005;30:742–746. 10.1038/sj.npp.1300592. [Abstract] [CrossRef] [Google Scholar]
- Flint J. Mapping quantitative traits and strategies to find quantitative trait genes. Methods. 2011;53:163–174. 10.1016/j.ymeth.2010.07.007. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Galloway GP, Singleton EG, Buscemi R, Baggott MJ, Dickerhoof RM, Mendelson JE, Methamphetamine Treatment Project Corporate Authors An examination of drug craving over time in abstinent methamphetamine users. The American Journal on Addictions. 2010;19:510–514. 10.1111/j.1521-0391.2010.00082.x. [Abstract] [CrossRef] [Google Scholar]
- Grisel JE, Belknap JK, O'Toole LA, Helms ML, Wenger CD, Crabbe JC. Quantitative trait loci affecting methamphetamine responses in BXD recombinant inbred mouse strains. The Journal of Neuroscience. 1997;17:745–754. 10.1523/JNEUROSCI.17-02-00745.1997. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Harkness JH, Shi X, Janowsky A, Phillips TJ. Trace Amine-Associated receptor 1 regulation of methamphetamine intake and related traits. Neuropsychopharmacology. 2015;40:2175–2184. 10.1038/npp.2015.61. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Hitzemann, Iancu, Reed, Baba, Lockwood, Phillips Regional analysis of the brain transcriptome in mice bred for high and low methamphetamine consumption. Brain Sciences. 2019;9:155. 10.3390/brainsci9070155. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Hogan B, Beddington R, Costantini F, Lacy E. Manipulating the Mouse Embryo: A Laboratory Manual. Cold Spring Harbor Press; 1994. Recovery, culture, and transfer of embryos and germ cells; pp. 127–188. [Google Scholar]
- Ide S, Kobayashi H, Tanaka K, Ujike H, Sekine Y, Ozaki N, Inada T, Harano M, Komiyama T, Yamada M, Iyo M, Ikeda K, Sora I. Gene polymorphisms of the mu opioid receptor in methamphetamine abusers. Annals of the New York Academy of Sciences. 2004;1025:316–324. 10.1196/annals.1316.039. [Abstract] [CrossRef] [Google Scholar]
- Ide S, Kobayashi H, Ujike H, Ozaki N, Sekine Y, Inada T, Harano M, Komiyama T, Yamada M, Iyo M, Iwata N, Tanaka K, Shen H, Iwahashi K, Itokawa M, Minami M, Satoh M, Ikeda K, Sora I. Linkage disequilibrium and association with methamphetamine dependence/psychosis of mu-opioid receptor gene polymorphisms. The Pharmacogenomics Journal. 2006;6:179–188. 10.1038/sj.tpj.6500355. [Abstract] [CrossRef] [Google Scholar]
- Iyer V, Boroviak K, Thomas M, Doe B, Riva L, Ryder E, Adams DJ. No unexpected CRISPR-Cas9 off-target activity revealed by trio sequencing of gene-edited mice. PLOS Genetics. 2018;14:e1007503. 10.1371/journal.pgen.1007503. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna JA, Charpentier E. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science. 2012;337:816–821. 10.1126/science.1225829. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Jing L, Li JX. Trace amine-associated receptor 1: a promising target for the treatment of psychostimulant addiction. European Journal of Pharmacology. 2015;761:345–352. 10.1016/j.ejphar.2015.06.019. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Lam VM, Mielnik CA, Baimel C, Beerepoot P, Espinoza S, Sukhanov I, Horsfall W, Gainetdinov RR, Borgland SL, Ramsey AJ, Salahpour A. Behavioral effects of a potential novel TAAR1 antagonist. Frontiers in Pharmacology. 2018;9:953. 10.3389/fphar.2018.00953. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Lehner B. Molecular mechanisms of epistasis within and between genes. Trends in Genetics. 2011;27:323–331. 10.1016/j.tig.2011.05.007. [Abstract] [CrossRef] [Google Scholar]
- Leo D, Mus L, Espinoza S, Hoener MC, Sotnikova TD, Gainetdinov RR. Taar1-mediated modulation of presynaptic dopaminergic neurotransmission: role of D2 dopamine autoreceptors. Neuropharmacology. 2014;81:283–291. 10.1016/j.neuropharm.2014.02.007. [Abstract] [CrossRef] [Google Scholar]
- Lindemann L, Meyer CA, Jeanneau K, Bradaia A, Ozmen L, Bluethmann H, Bettler B, Wettstein JG, Borroni E, Moreau JL, Hoener MC. Trace amine-associated receptor 1 modulates dopaminergic activity. Journal of Pharmacology and Experimental Therapeutics. 2008;324:948–956. 10.1124/jpet.107.132647. [Abstract] [CrossRef] [Google Scholar]
- Liu JF, Siemian JN, Seaman R, Zhang Y, Li JX. Role of TAAR1 within the subregions of the mesocorticolimbic dopaminergic system in Cocaine-Seeking behavior. The Journal of Neuroscience. 2017;37:882–892. 10.1523/JNEUROSCI.2006-16.2016. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Liu JF, Li JX. TAAR1 in addiction: looking beyond the tip of the iceberg. Frontiers in Pharmacology. 2018;9:279. 10.3389/fphar.2018.00279. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Lominac KD, McKenna CL, Schwartz LM, Ruiz PN, Wroten MG, Miller BW, Holloway JJ, Travis KO, Rajasekar G, Maliniak D, Thompson AB, Urman LE, Phillips TJ, Szumlinski KK. Mesocorticolimbic monoamine correlates of methamphetamine sensitization and motivation. Frontiers in Systems Neuroscience. 2014;8:70. 10.3389/fnsys.2014.00070. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Lominac KD, Quadir SG, Barrett HM, McKenna CL, Schwartz LM, Ruiz PN, Wroten MG, Campbell RR, Miller BW, Holloway JJ, Travis KO, Rajasekar G, Maliniak D, Thompson AB, Urman LE, Kippin TE, Phillips TJ, Szumlinski KK. Prefrontal glutamate correlates of methamphetamine sensitization and preference. European Journal of Neuroscience. 2016;43:689–702. 10.1111/ejn.13159. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Mianné J, Chessum L, Kumar S, Aguilar C, Codner G, Hutchison M, Parker A, Mallon AM, Wells S, Simon MM, Teboul L, Brown SD, Bowl MR. Correction of the auditory phenotype in C57BL/6N mice via CRISPR/Cas9-mediated homology directed repair. Genome Medicine. 2016;8:16. 10.1186/s13073-016-0273-4. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Misanin JR, Wilson HA, Schwarz PR, Tuschak JB, Hinderliter CF. Low body temperature affects associative processes in long-trace conditioned flavor aversion. Physiology & Behavior. 1998;65:581–590. 10.1016/S0031-9384(98)00212-1. [Abstract] [CrossRef] [Google Scholar]
- Misanin JR, Anderson MJ, Christianson JP, Collins MM, Goodhart MG, Rushanan SG, Hinderliter CF. Low body temperature, time dilation, and long-trace conditioned flavor aversion in rats. Neurobiology of Learning and Memory. 2002;78:167–177. 10.1006/nlme.2001.4048. [Abstract] [CrossRef] [Google Scholar]
- Mulligan MK, Mozhui K, Prins P, Williams RW. GeneNetwork: a toolbox for systems genetics. Methods in Molecular Biology. 2017;1488:75–120. 10.1007/978-1-4939-6427-7_4. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Nadeau JH, Auwerx J. The virtuous cycle of human genetics and mouse models in drug discovery. Nature Reviews Drug Discovery. 2019;18:255–272. 10.1038/s41573-018-0009-9. [Abstract] [CrossRef] [Google Scholar]
- Nakajima K, Kazuno AA, Kelsoe J, Nakanishi M, Takumi T, Kato T. Exome sequencing in the knockin mice generated using the CRISPR/Cas system. Scientific Reports. 2016;6:34703. 10.1038/srep34703. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- National Institute on Drug Abuse Overdose death rates. [June 3, 2019];2018 https://www.drugabuse.gov/related-topics/trends-statistics/overdose-death-rates
- National Institute on Drug Abuse Methamphetamine. [June 3, 2019];2019 https://www.drugabuse.gov/publications/research-reports/methamphetamine/
- National Research Council . Guide for the Care and Use of Laboratory Animals. 8th ed. Washington, DC: The National Academies Press; 2011. [CrossRef] [Google Scholar]
- Pei Y, Asif-Malik A, Hoener M, Canales JJ. A partial trace amine-associated receptor 1 agonist exhibits properties consistent with a methamphetamine substitution treatment. Addiction Biology. 2017;22:1246–1256. 10.1111/adb.12410. [Abstract] [CrossRef] [Google Scholar]
- Reed C, Baba H, Zhu Z, Erk J, Mootz JR, Varra NM, Williams RW, Phillips TJ. A spontaneous mutation in Taar1 Impacts Methamphetamine-Related Traits Exclusively in DBA/2 Mice from a Single Vendor. Frontiers in Pharmacology. 2018;8:993. 10.3389/fphar.2017.00993. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Revel FG, Moreau JL, Gainetdinov RR, Bradaia A, Sotnikova TD, Mory R, Durkin S, Zbinden KG, Norcross R, Meyer CA, Metzler V, Chaboz S, Ozmen L, Trube G, Pouzet B, Bettler B, Caron MG, Wettstein JG, Hoener MC. TAAR1 activation modulates monoaminergic neurotransmission, preventing hyperdopaminergic and hypoglutamatergic activity. PNAS. 2011;108:8485–8490. 10.1073/pnas.1103029108. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Rutigliano G, Accorroni A, Zucchi R. The case for TAAR1 as a modulator of central nervous system function. Frontiers in Pharmacology. 2017;8:987. 10.3389/fphar.2017.00987. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Shabani S, McKinnon CS, Reed C, Cunningham CL, Phillips TJ. Sensitivity to rewarding or aversive effects of methamphetamine determines methamphetamine intake. Genes, Brain and Behavior. 2011;10:625–636. 10.1111/j.1601-183X.2011.00700.x. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Shabani S, Dobbs LK, Ford MM, Mark GP, Finn DA, Phillips TJ. A genetic animal model of differential sensitivity to methamphetamine reinforcement. Neuropharmacology. 2012a;62:2169–2177. 10.1016/j.neuropharm.2012.01.002. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Shabani S, McKinnon CS, Cunningham CL, Phillips TJ. Profound reduction in sensitivity to the aversive effects of methamphetamine in mice bred for high methamphetamine intake. Neuropharmacology. 2012b;62:1134–1141. 10.1016/j.neuropharm.2011.11.005. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Shabani S, Houlton SK, Hellmuth L, Mojica E, Mootz JR, Zhu Z, Reed C, Phillips TJ. A mouse model for Binge-Level methamphetamine use. Frontiers in Neuroscience. 2016;10:493. 10.3389/fnins.2016.00493. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Shabani S, Schmidt B, Ghimire B, Houlton SK, Hellmuth L, Mojica E, Phillips TJ. Depression-like symptoms of withdrawal in a genetic mouse model of binge methamphetamine intake. Genes, Brain and Behavior. 2019;18:e12533. 10.1111/gbb.12533. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Shi X, Walter NA, Harkness JH, Neve KA, Williams RW, Lu L, Belknap JK, Eshleman AJ, Phillips TJ, Janowsky A. Genetic polymorphisms affect mouse and human trace Amine-Associated receptor 1 function. PLOS ONE. 2016;11:e0152581. 10.1371/journal.pone.0152581. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Simon SL, Richardson K, Dacey J, Glynn S, Domier CP, Rawson RA, Ling W. A comparison of patterns of methamphetamine and cocaine use. Journal of Addictive Diseases. 2002;21:35–44. 10.1300/J069v21n01_04. [Abstract] [CrossRef] [Google Scholar]
- Solberg Woods LC. QTL mapping in outbred populations: successes and challenges. Physiological Genomics. 2014;46:81–90. 10.1152/physiolgenomics.00127.2013. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Stalder H, Hoener MC, Norcross RD. Selective antagonists of mouse trace amine-associated receptor 1 (mTAAR1): discovery of EPPTB (RO5212773) Bioorganic & Medicinal Chemistry Letters. 2011;21:1227–1231. 10.1016/j.bmcl.2010.12.075. [Abstract] [CrossRef] [Google Scholar]
- Szumlinski KK, Lominac KD, Campbell RR, Cohen M, Fultz EK, Brown CN, Miller BW, Quadir SG, Martin D, Thompson AB, von Jonquieres G, Klugmann M, Phillips TJ, Kippin TE. Methamphetamine addiction vulnerability: the glutamate, the bad, and the ugly. Biological Psychiatry. 2017;81:959–970. 10.1016/j.biopsych.2016.10.005. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- United Nations Office on Drugs and Crime . World Drug Report. United Nations Publication; 2016. [CrossRef] [Google Scholar]
- Wheeler JM, Reed C, Burkhart-Kasch S, Li N, Cunningham CL, Janowsky A, Franken FH, Wiren KM, Hashimoto JG, Scibelli AC, Phillips TJ. Genetically correlated effects of selective breeding for high and low methamphetamine consumption. Genes, Brain and Behavior. 2009;8:758–771. 10.1111/j.1601-183X.2009.00522.x. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Willi M, Smith HE, Wang C, Liu C, Hennighausen L. Mutation frequency is not increased in CRISPR-Cas9-edited mice. Nature Methods. 2018;15:756–758. 10.1038/s41592-018-0148-2. [Abstract] [CrossRef] [Google Scholar]
- Wolinsky TD, Swanson CJ, Smith KE, Zhong H, Borowsky B, Seeman P, Branchek T, Gerald CP. The trace Amine 1 receptor knockout mouse: an animal model with relevance to schizophrenia. Genes, Brain and Behavior. 2007;6:628–639. 10.1111/j.1601-183X.2006.00292.x. [Abstract] [CrossRef] [Google Scholar]
- Xie Z, Miller GM. Beta-phenylethylamine alters monoamine transporter function via trace amine-associated receptor 1: implication for modulatory roles of trace amines in brain. Journal of Pharmacology and Experimental Therapeutics. 2008;325:617–628. 10.1124/jpet.107.134247. [Abstract] [CrossRef] [Google Scholar]
- Zhang F, inventor. Systems methods and compositions for sequence manipulation. U.S. Provisional Patent Application 61/736,527. December 12, 2012; Published at US008697359B1 (Awarded) 2012 https://patents.google.com/patent/WO2015089486A3/en
Decision letter
Thank you for submitting your article "Taar1 gene variants have a causal role in methamphetamine intake and response and interact with Oprm1" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Catherine Dulac as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Amelie Baud (Reviewer #1) and Bryan Yamamoto (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
The manuscript describes experiments that demonstrate the role of a DBA/2J specific Taar1 allele on methamphetamine consumption and methamphetamine-induced hypothermia. Drawing upon previous work, the authors test the involvement of a Taar1 allele that arose spontaneously in the DBA/2J colony at the Jackson Laboratory. Using CRISPR-Cas9 technology they generate a mouse line that replaced the DBA/2J Taar1 allele (Taar1m1J) with a wildtype allele (Taar1+) in mice selectively bred for high consumption of methamphetamine. The authors confirm the effects of Taar1 on methamphetamine drinking as well as methamphetamine-induced hypothermia, and present evidence for an interaction between Taar1 and a locus on chromosome 10 (Oprm1) using a set of BxD RI strains
Essential revisions:
1) Please use a significance threshold that takes into account the number of tests employed. This needs to be stated and the threshold justified.
2) The control line MAHDR-Taar1m1J was derived from mice in which the allele was not successfully altered during CRISPR experiments. However, it's not clear whether these lines are maintained separately or bred together – please state whether the controls are littermates of the MAHDR-Taar1+ mice.
3) CRISPR can have off-target effects. The authors provide one sentence in the Materials and methods section saying there was an "initial" analysis of potential off-target sites but provide no further information. It's not clear whether any follow up studies of off-target effects were conducted. These details should be provided in the Materials and methods and Results and any potential caveats presented in the Discussion.
4) It seems a bit misleading to characterize the hypothermic effect of MA as prolonged in the Taar1+/+/Oprm1D2/D2 mice (Figure 7B). The rate of increase in temperature AFTER the initial decrease seems fairly consistent between these mice and the Taar1+/+/Oprm1B6/B6 mice. In comparison to T0, it may seem that the hypothermic effect is prolonged due to the significant difference in the temperature at T30, but the rate of increase wasn't measured so it's hard to conclude that Taar1+/+/Oprm1B6/B6 mice recover more quickly.
5) Please can you explain whether the interaction of Oprm1 genotype with Taar1 has been examined in the QTL studies for MA drinking and hypothermia – perhaps by using Oprm1 genotype as a covariate? Although there may not have been enough power to study gene x gene interactions due to sample size, it would be good to see this effect confirmed (or if not explain if this is due to low power)
6) The discussion about the relationship between MA consumption and body temperature (section “Relationship Between MA consumption and Body Temperature”) is confusing. It makes sense that MA-induced hypothermia could result in reduced MA drinking due to aversion, but the alternative hypothesis that MA drinking somehow alters MA-induced hypothermia is not so obvious. More explanation would be beneficial as to how Taar1 mutation might alter meth consumption. For example if these mutants cannot appropriately thermoregulate, do they have a depressed phenotype, or have an altered mesocorticolimbic pathway?
7) The authors make several comments in the discussion about the potential translational importance of Taar1 in human MA use and as a potential therapeutic. They should include a few sentences in the discussion about MA drinking in rodents as a model for human MA use. Are the pharmacokinetics of MA following oral consumption similar to normal human routes of MA administration?
8) One consideration that was not addressed was whether the genetic modifications were simply affecting the potency of METH. If METH acts on Taar1, then it seems reasonable that a non-functional Taar1 would require higher METH intake to produce a similar effect. Please comment on this suggestion.
Contributor Information
Jonathan Flint, University of California, Los Angeles, United States.
Amelie Baud, EMBL - EBI, United Kingdom.
Bryan Yamamoto, Indiana University School of Medicine, United States.
Author response
Essential revisions:
1) Please use a significance threshold that takes into account the number of tests employed. This needs to be stated and the threshold justified.
We have now used a Bonferroni correction to account for multiple tests in the revised manuscript, and have stated this in the “Data analysis” section; p-values have been adjusted throughout. These corrections for multiple tests did not change our conclusions. The only changes from the initial submission involved the saline group in the BXD hypothermia study, and the Taar1m1J/m1Jmice, which do not exhibit MA-induced hypothermia. Thus, none of the results for MA effects were impacted.
2) The control line MAHDR-Taar1m1J was derived from mice in which the allele was not successfully altered during CRISPR experiments. However, it's not clear whether these lines are maintained separately or bred together – please state whether the controls are littermates of the MAHDR-Taar1+ mice.
The CRISPR knock-in and control lines were maintained as within-line breeding pairs in a shared colony room. This breeding of individuals with the same Taar1 genotype is consistent with the maintenance of all other populations in which the Taar1 SNP exists (e.g., MAHDR vs. MALDR; BXD strains; D2 mice from The Jackson Laboratory vs. D2 mice from other suppliers; Reed et al., 2018). We now state this in the section “Generation of MAHDR-Taar1+/+ knock-in mice”. Please note that when an F2 was created between the DBA/2J and C57BL/6J mice, resulting in littermates with different Taar1 genotypes, an equally strong association of Taar1 genotype with MA intake was found (see Reed et al., 2018). We have added a discussion point in the section “Potential Limitations”.
3) CRISPR can have off-target effects. The authors provide one sentence in the Materials and methods section saying there was an "initial" analysis of potential off-target sites but provide no further information. It's not clear whether any follow up studies of off-target effects were conducted. These details should be provided in the Materials and methods and Results and any potential caveats presented in the Discussion.
We have now included more detailed methodological information in the paper. We have also included discussion regarding potential off-target effects. In our situation, in which the sequence change was created on a heterogeneous background, sequencing the rest of the genome of the CRISPR knock-in and control mice would not provide specific information about off-target effects. This is because the existing genetic variance would not be separable from CRISPR-induced effects. A study published in Nature Methods (Schaefer et al., 2017) reported extensive off-target effects of CRISPR and has since been retracted, due to significant concerns raised in multiple commentaries *(Kim et al., 2018; Lareau et al., 2018; Lescarbeau et al., 2018; Nutter et al., 2018; Wilson et al., 2018). There appears to be some consensus that off-target effects are rare and indistinguishable from the background rate of de novo mutation (Anderson et al., 2018; Ayabe et al., 2018; Iyer et al., 2018; Mianne et al., 2016; Nakajima et al., 2016; Willi et al., 2018). We have mentioned this consensus in the paper. In addition, we describe our use of the basic local alignment search tool (BLAST) to examine potential sequence identity for the guide RNAs generated by ThermoFisher and give details about the one that we chose that was a perfect match, except for the relevant SNP. The gRNA sequence had similarity to 9 additional regions on chromosomes 1, 3, 5, 9, 11 and 19, but with 2-3 mismatches in each case (see t Design of the CRISPR-Cas9 knock-in of Taar1m1J and CRISPR reagents”). As can be seen in Figure 8 of the current paper, and in our published work (Belknap et al., 2013), there are no mapped QTLs for MA drinking or MA-induced hypothermia on any of these chromosomes. Therefore, even if there was a genetic change in any of these locations, it would be extremely unlikely that it would produce the large changes in both MA intake and MA-induced hypothermia that we observed. On the other hand, because the chromosome 10 QTL accounts for 60% of the genetic variance in MA intake and Taar1 genotype-phenotype correlations range from r=0.81 to 0.96 for 40 mg/L MA intake and r=0.71 to 0.82 for 2 mg/kg MA-induced hypothermia (Reed et al., 2018; and data herein), the large effects we observed would be expected for replacement of Taar1m1J with Taar1+. Please see this added discussion in the section “Potential Limitations”.
*Schaefer and Commentary citations:
Kim ST, Park J, Kim D, Kim K, Bae S, Schlesner M, Kim JS. 2018. Response to “Unexpected mutations after CRISPR-Cas9 editing in vivo.” Nat Methods 15: 239-240. 10.1038/nmeth.4554
Lareau CA, Clement K, Hsu JY, Pattanayak V, Joung JK, Aryee MJ, Pinello L. 2018. Response to “Unexpected mutations after CRISPR-Cas9 editing in vivo.” Nat Methods 15: 238-239. 10.1038/nmeth.4541
Lescarbeau RM, Murray B, Barnes TM, Bermingham N. 2018. Response to “Unexpected mutations after CRISPR-Cas9 editing in vivo.” Nat Methods 15: 237. 10.1038/nmeth.4553
Nutter LMJ, Heaney JD, Lloyd KCK, Murray SA, Seavitt JR, Skarnes WC, Teboul L, Brown SDM, Moore M. 2018. Response to “Unexpected mutations after CRISPR-Cas9 editing in vivo.” Nat Methods 15: 235-236. 10.1038/nmeth.4559
Wilson CJ, Fennell T, Bothmer A, Maeder ML, Reyon D, Cotta-Ramusino C, Fernandez CA, Marco E, Barrera LA, Jayaram H, Albright CF, Cox GF, Church GM, Myer VE. 2018. Response to “Unexpected mutations after CRISPR-Cas9 editing in vivo.” Nat Methods 15:236-237. 10.1038/nmeth.4552
Schaefer KA, Wu WH, Colgan DF, Tsang SH, Bassuk AG, Mahajan VB. 2017. Unexpected mutations after CRISPR-Cas9 editing in vivo. Nat Methods 14:547-548. 10.1038/nmeth.4293
4) It seems a bit misleading to characterize the hypothermic effect of MA as prolonged in the Taar1+/+/Oprm1D2/D2 mice (Figure 7B). The rate of increase in temperature after the initial decrease seems fairly consistent between these mice and the Taar1+/+/Oprm1B6/B6 mice. In comparison to T0, it may seem that the hypothermic effect is prolonged due to the significant difference in the temperature at T30, but the rate of increase wasn't measured so it's hard to conclude that Taar1+/+/Oprm1B6/B6 mice recover more quickly.
We see the reviewers’ point and agree. We have edited the document regarding this point in the section “MA-induced Hypothermia is Impacted in BXD Mice by a Taar1 x Oprm1 Interaction”, and now only comment on the difference in the magnitude of effect between Taar1+/+/Oprm1D2/D2and Taar1+/+/Oprm1B6/B6mice, as this was a significant result. We note that the rate of return to baseline is similar between the two genotypes.
5) Please can you explain whether the interaction of Oprm1 genotype with Taar1 has been examined in the QTL studies for MA drinking and hypothermia – perhaps by using Oprm1 genotype as a covariate? Although there may not have been enough power to study gene x gene interactions due to sample size, it would be good to see this effect confirmed (or if not explain if this is due to low power)
Thank you for making this point! We now include composite interval QTL maps (Figure 8C-D) for MA consumption and hypothermia that used Oprm1 genotype as a co-factor. This resulted in increased LOD scores for both traits, that was significant for the thermal trait and a strong trend for the MA consumption traits. This additional analysis strengthens our conclusions that formerly were based on the behavioral results. We now describe the composite mapping in the Materials and methods (section “QTL Analysis”), Results (section”MA Consumption and Body Temperature Response to MA Map to a Region of Chromosome 10 at the Location of Taar1”), and Discussion (section “Independent and Epistatic Effects of Taar1 and Oprm1”).
6) The discussion about the relationship between MA consumption and body temperature (section “Relationship Between MA consumption and Body Temperature”) is confusing. It makes sense that MA-induced hypothermia could result in reduced MA drinking due to aversion, but the alternative hypothesis that MA drinking somehow alters MA-induced hypothermia is not so obvious. More explanation would be beneficial as to how Taar1 mutation might alter meth consumption. For example if these mutants cannot appropriately thermoregulate, do they have a depressed phenotype, or have an altered mesocorticolimbic pathway?
It is unlikely that the difference in MA effect on body temperature is due to inability to thermoregulate for several reasons. Mice with the Taar1m1J/m1Jgenotype do not differ in baseline temperatures or changes in those temperatures over time, compared to mice with the Taar1+/+genotype (Harkness et al., 2015, and data in the current manuscript). In addition, mice of the 2 genotypes exhibit comparable time-dependent patterns of response to other drugs with hypothermic effects, including ethanol (Harkness et al., 2015) and cocaine (unpublished data). Thus, the mutants display a hypothermic response, and the lack of hypothermia in response to MA is a MA-specific effect. Further, the Taar1 mutation does not appear to affect locomotor activity at baseline or in response to this dose of MA (Shabani et al., 2011; Wheeler et al., 2009), so it is not likely that Taar1 effects are due to differences in movement after treatment. We now discuss this in the section “Relationship Between MA Consumption and Body Temperature”. Also, we discuss other possible mechanisms by which TAAR1 may affect MA consumption (including alterations in the mesocorticolimbic pathway) in the last paragraph of the section “Taar1 and MA Consumption”.
We have not yet examined body temperature changes occurring during the time that MA is being consumed. We have observed that low line mice consume similar amounts of MA as high line mice on the first day that MA is offered, but then decrease their consumption on the subsequent day, remaining at low intake levels from then on (Shabani et al., 2012a; Eastwood et al., 2014). It is possible that once MA reaches a high enough dose during the first drinking session, that the low line mice experience hypothermia and that this plays a role in their reduced intake on subsequent days. We are in the process of obtaining equipment to allow us to examine body temperature simultaneously with locomotor activity during drinking sessions, which will allow us to consider these questions. We now discuss this in the section “Relationship Between MA Consumption and Body Temperature”.
7) The authors make several comments in the Discussion about the potential translational importance of Taar1 in human MA use and as a potential therapeutic. They should include a few sentences in the Discussion about MA drinking in rodents as a model for human MA use. Are the pharmacokinetics of MA following oral consumption similar to normal human routes of MA administration?
We discuss this in Shabani et al. (2012a), but agree that it is important to reiterate. Please see our discussion of the relevance of this mouse model in the section “Taar1 and MA Consumption”.
8) One consideration that was not addressed was whether the genetic modifications were simply affecting the potency of METH. If METH acts on Taar1, then it seems reasonable that a non-functional Taar1 would require higher METH intake to produce a similar effect. Please comment on this suggestion.
If the mutation produced a sub-functional receptor, this would of course be possible. However, it appears that the receptor produced is non-functional. The D2-like isoform of TAAR1, expressed by Taar1m1J, exhibits no cAMP response to a wide range of MA concentrations or to other TAAR1 agonists (Harkness et al., 2015; Shi et al., 2016). We have shown that, even at high doses of MA (up to 16 mg/kg), high line mice fail to exhibit hypothermia, whereas MALDR mice express hypothermia at doses as low as 1 mg/kg (Harkness et al., 2015). Thus, we are fairly certain that the genetic modification did not simply affect potency. We now discuss this in the section “Potential Limitations”.
Articles from eLife are provided here courtesy of eLife Sciences Publications, Ltd
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
Article citations
Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice.
Elife, 12:RP89096, 21 Feb 2024
Cited by: 3 articles | PMID: 38381482 | PMCID: PMC10942616
Discovery and validation of genes driving drug-intake and related behavioral traits in mice.
Genes Brain Behav, e12875, 02 Jan 2024
Cited by: 1 article | PMID: 38164795 | PMCID: PMC10780947
Neuron Navigator 1 (Nav1) regulates the response to cocaine in mice.
Commun Biol, 6(1):1053, 18 Oct 2023
Cited by: 0 articles | PMID: 37853211 | PMCID: PMC10584906
Rats have low motivation to self-administer oral methamphetamine across increasing response requirements.
Behav Brain Res, 455:114673, 16 Sep 2023
Cited by: 1 article | PMID: 37717660 | PMCID: PMC10591766
Robust aversive effects of trace amine-associated receptor 1 activation in mice.
Neuropsychopharmacology, 48(10):1446-1454, 13 Apr 2023
Cited by: 0 articles | PMID: 37055488 | PMCID: PMC10425385
Go to all (18) article citations
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
RefSeq - NCBI Reference Sequence Database (2)
- (2 citations) RefSeq - NM_053205.1
- (1 citation) RefSeq - NM_001304955
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Differential genetic risk for methamphetamine intake confers differential sensitivity to the temperature-altering effects of other addictive drugs.
Genes Brain Behav, 19(5):e12640, 15 Jan 2020
Cited by: 4 articles | PMID: 31925906 | PMCID: PMC7286770
Trace Amine-Associated Receptor 1 Regulation of Methamphetamine Intake and Related Traits.
Neuropsychopharmacology, 40(9):2175-2184, 05 Mar 2015
Cited by: 47 articles | PMID: 25740289 | PMCID: PMC4613607
Confirmation of a Causal Taar1 Allelic Variant in Addiction-Relevant Methamphetamine Behaviors.
Front Psychiatry, 12:725839, 26 Aug 2021
Cited by: 3 articles | PMID: 34512422 | PMCID: PMC8428522
An animal model of differential genetic risk for methamphetamine intake.
Front Neurosci, 9:327, 23 Sep 2015
Cited by: 16 articles | PMID: 26441502 | PMCID: PMC4585292
Review
Free full text in Europe PMC
Funding
Funders who supported this work.
BLRD VA (3)
Grant ID: I01 BX002758
Grant ID: I01 BX003279
Grant ID: I01 BX002106
Department of Veterans Affairs (3)
Grant ID: I01BX002106
Grant ID: I01BX003279
Grant ID: I01BX002758
NIDA NIH HHS (7)
Grant ID: P50 DA018165
Grant ID: R01 DA034388
Grant ID: T32 DA007262
Grant ID: U01 DA041579
Grant ID: P30 DA044223
Grant ID: R01 DA046081
Grant ID: T32 DA007027
National Institute on Drug Abuse (4)
Grant ID: R01 DA046081
Grant ID: U01 DA041579
Grant ID: P30 DA044223
Grant ID: P50 DA018165
Oregon Health & Science University (1)
Grant ID: University Shared Resource award
U.S. Department of Veterans Affairs (1)
Grant ID: Research Career Scientist Program - Career Scientist award
University of Tennessee Center for Integrative and Translational Science (1)
Grant ID: Center support
VA (3)
Grant ID: I01BX002758
Grant ID: I01BX003279
Grant ID: I01BX002106