Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Articles, Behavioral/Systems/Cognitive

Identification of Candidate Genes and Gene Networks Specifically Associated with Analgesic Tolerance to Morphine

Jenica D. Tapocik, Noah Letwin, Cheryl L. Mayo, Bryan Frank, Troung Luu, Ovokeraye Achinike, Carrie House, Russell Williams, Greg I. Elmer and Norman H. Lee
Journal of Neuroscience 22 April 2009, 29 (16) 5295-5307; DOI: https://doi.org/10.1523/JNEUROSCI.4020-08.2009
Jenica D. Tapocik
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Noah Letwin
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cheryl L. Mayo
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bryan Frank
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Troung Luu
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ovokeraye Achinike
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Carrie House
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Russell Williams
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Greg I. Elmer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Norman H. Lee
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Chronic morphine administration may alter the expression of hundreds to thousands of genes. However, only a subset of these genes is likely involved in analgesic tolerance. In this report, we used a behavior genetics strategy to identify candidate genes specifically linked to the development of morphine tolerance. Two inbred genotypes [C57BL/6J (B6), DBA2/J (D2)] and two reciprocal congenic genotypes (B6D2, D2B6) with the proximal region of chromosome 10 (Chr10) introgressed into opposing backgrounds served as the behavior genetic filter. Tolerance after therapeutically relevant doses of morphine developed most rapidly in the B6 followed by the B6D2 genotype and did not develop in the D2 mice and only slightly in the D2B6 animals indicating a strong influence of the proximal region of Chr10 in the development of tolerance. Gene expression profiling and pattern matching identified 64, 53, 86, and 123 predisposition genes and 81, 96, 106, and 82 tolerance genes in the periaqueductal gray (PAG), prefrontal cortex, temporal lobe, and ventral striatum, respectively. A potential gene network was identified in the PAG in which 19 of the 34 genes were strongly associated with tolerance. Eleven of the network genes were found to reside in quantitative trait loci previously associated with morphine-related behaviors, whereas seven were predictive of tolerance (morphine-naive condition). Overall, the genes modified by chronic morphine administration show a strong presence in canonical pathways representative of neuroadaptation. A potentially significant role for the micro-RNA and epigenetic mechanisms in response to chronic administration of pharmacologically relevant doses of morphine was highlighted by candidate genes Dicer and H19.

Introduction

The therapeutic value of opioids for pain relief is well established. Repeated opioid administration, however, can lead to a progressive decline in analgesic efficacy. Extensive research into the mechanisms involved in analgesic tolerance has provided substantial insight into its neurobiology (Williams et al., 2001). However, despite recent advances, the cellular and molecular alterations downstream of the primary effects at the μ-opioid receptor (μOR) are complex and not well understood (Law et al., 2000).

One strategy to discover the mechanisms of the behavioral effects of morphine is to identify regions of the genome associated with these traits through linkage analysis. Behavior genetic investigations of morphine tolerance demonstrate a high degree of heritability (Mas et al., 2000; Kest et al., 2002; Liang et al., 2006). Behavioral quantitative trait locus (bQTL) mapping assigns a significant portion of the genotype-dependent variance in morphine consumption and acute analgesia to the Mop2 locus on the proximal region of Chr10 (Berrettini et al., 1994a; Belknap et al., 1995; Bergeson et al., 2001; Doyle et al., 2006). This is a promising candidate region since the gene encoding the μOR (Oprm) is located in Mop2. Reciprocal congenic genotypes in which this region has been introgressed from the opposing background [C57BL/6J (B6) or DBA/2J (D2)] have been developed to further aid in the isolation of genetic factors (Bergeson et al., 2001; Ferraro et al., 2005). These studies have established a basic strategy with which to pursue specific genomic regions and candidate genes.

Large-scale gene expression profiling studies have also been used as a strategy to find candidate genes involved in tolerance. For instance, gene expression profiling in four opioid-naive (Korostynski et al., 2006) and morphine-treated mouse genotypes have identified basal and morphine-induced gene expression differences in the striatum (Korostynski et al., 2007) and nucleus accumbens (Grice et al., 2007) that may contribute to opioid-related phenotypes. However, none of the aforementioned studies anchored their gene expression results to within-subjects derived behavioral data or investigated multiple brain regions especially to include those more directly involved in opioid-induced analgesia [i.e., periaqueductal gray (PAG)] and to assess region by genotype and drug interactions. Characterizing several brain regions and anchoring expression data to the same experimental protocol and dosing schedule may improve the likelihood of focusing on genes directly relevant to the phenotype in question.

Reciprocal congenic genotypes have been an excellent tool for confirming regions of the genome hypothesized to be involved in a particular trait (Bergeson et al., 2001; Ferraro et al., 2005). The purpose of our study was to use reciprocal congenics to investigate the promising role of proximal Chr10 in nociception and the acute and chronic effects of pharmacologically relevant doses of morphine. A behavior genetics approach (Letwin et al., 2006) was used to identify changes in gene expression directly related to morphine analgesic tolerance. In this manner, morphine-induced changes in gene expression that are irrelevant to this trait were filtered from the pool of genes. Only those genes specifically correlated to the trait of interest were identified as candidate genes significantly influencing the development of tolerance.

Materials and Methods

Animals

Adult male C57BL/6J (B6) and DBA/J (D2) (The Jackson Laboratory), D2B6 (Mop2 locus from B6 introgressed into D2 background), and B6D2 (Mop2 locus from D2 introgressed into B6 background) mice 60–120 d of age and weighing ∼21–28 g at the start of the experiment were used. The D2B6 and B6D2 mice were shipped to the Maryland Psychiatric Research Center from the University of Pennsylvania in which brother by sister mating continued to propagate the reciprocal strains. All animals were experimentally naive, housed in a temperature-controlled room (21°C; 7:00 A.M. to 7:00 P.M. lights on), and given ad libitum access to chow and tap water during the entire experimental procedure.

The reciprocal congenic genotypes were generated in the following manner. The portion of the chromosome thought to account for a significant portion of variance in response to morphine was reciprocally introgressed between the parental genotypes, B6 and D2. The introgress protocol began with a B6 × D2 intercross followed by backcross of the resulting F1 generation progeny to both parental genotypes. The backcross consisted of F1 males and females being mated with progenitor mice. The N2 backcross progeny were genotyped using DNA extracted from tail snips at the time of weaning. Mice were then selected for subsequent rounds of backcross breeding based on heterozygosity at specific set of microsatellite DNA markers [B6D2: D10Mit75 (2.0 cM) proximally to D10Mit124 (15.0 cM) distally; D2B6: D10Mit75 (2.0 cM) proximally to D10Mit61 (32.0 cM) distally] (Ferraro et al., 2005). Ten serial backcross generations were created in a reciprocal manner. At generation N10, heterozygous mice were intercrossed and DNA from the N10 F1 progeny was analyzed to identify homozygous congenic mice D2B6 and reciprocal congenic genotype B6D2 (supplemental Fig. 1, available at www.jneurosci.org as supplemental material).

Behavioral testing: nociception and analgesia

Baseline sensitivity to pain and morphine-induced analgesia was measured by the hot-plate test (Harvard Hot Plate Analgesia Meter; Harvard Apparatus). The hot plate was maintained at a constant temperature of 55°C for all test conditions. Latency to paw lick (front or hind) or jump was used as the dependent measure. Under all experimental conditions, a cutoff time of four times the saline control value of each genotype was used to avoid tissue damage.

Baseline sensitivity.

To determine the baseline sensitivity for each genotype, we examined the latencies of each genotype naive to morphine. Each genotype received one injection of saline and was tested 15 min later on the hot plate (n = 6–9 animals/genotype). Data were analyzed using a one-way ANOVA with genotype as a main effect.

Acute analgesic sensitivity to morphine.

The subject was placed on the hot plate 15 min after subcutaneous injection of 0, 0.3, 1.0, 3.0, 10.0, or 30.0 mg/kg morphine to determine the relative potency of morphine to produce analgesia. Each dose for each genotype consisted of six to nine animals and was tested for only 1 d. A two-way ANOVA was performed for dose–response curves across genotype. Data were analyzed as the percentage of maximal analgesic response (%MPE) as determined by the following formula: 100 * [(test latency) − (saline baseline latency)] ÷ [(4 * baseline latency) − (baseline latency)]. The ED90 for each genotype was derived from the regression analysis of the linear portion of each dose–response curve.

Chronic dosing and the development of tolerance.

Analgesic tolerance to morphine was determined in three distinct experiments. In each experiment, the dose of morphine required to produce ∼90% (ED90 as determined above) of the maximal possible analgesic effect (%MPE) was used to induce tolerance in each genotype. The ED90 dose of morphine for the D2, D2B6, B6, and B6D2 mice was as follows: 7.7, 7.7, 31.2, and 31.2 mg/kg, respectively. Mice were tested on the hot plate 15 min after injection. Saline or the ED90 dose of morphine was administered subcutaneously between the hours of 8:00 and 11:00 A.M.

The three experimental conditions were as follows: (1) “Three injection”: For each genotype, two groups of mice (saline and morphine, N = 12–16/group) were tested once every other day during a 5 d period (i.e., Monday, Wednesday, Friday). Brains were harvested for gene expression analysis only in this experimental condition. (2) “Five injection, phenotypic confirmation”: This test was performed as a complementary means to confirm the results of the three-injection schedule and to push the dosing to determine whether D2 mice would develop tolerance using a shorter injection interval and more injections. For each genotype, two groups of mice (saline and morphine, N = 6–9/group) were tested once every day during a 5 d period. (3) “Learning assessment”: A subset of genotypes were run in a manner that would determine to what degree, if any, the response to morphine of a genotype was influenced by repeated exposure to the hot plate. In this experiment, the ED90 dose of morphine was given in a manner exactly the same as the three-dose group (once every other day; Monday, Wednesday, Friday); however, in one group, the hot plate was turned off until the final day of testing, whereas the hot plate was turned on (55°C) for the other group as described above. Only B6D2 and B6 mice were run in this experiment since they were the two genotypes that develop tolerance to a degree in which a learning component, if present, might be discerned (n = 6–9/group per genotype). The mice were tested on the hot plate 15 min after injection.

In the three-injection and five-injection conditions, a one-way ANOVA was used to determine whether equianalgesic effects of morphine were produced on the first day of drug administration. A two-way (genotype by day) repeated-measures ANOVA was used to determine genetic differences in the analgesic effects of morphine across repeated administration. Morphine-treated mice at each respective dosing interval were compared with saline-treated mice tested under the same procedure. Data are presented as %MPE. The three-injection treatment was compared with the five-injection phenotype confirmation group using a three-way repeated-measures ANOVA (genotype by treatment paradigm by injection number). Two comparisons were made, one using the first three tests of each paradigm and second using the first, third, and fifth tests of the five-injection treatment paradigm (compared with the three tests of the three-injection paradigm). In the learning assessment paradigm, a one-way ANOVA (heated vs room temperature hot plate) was conducted on the last test session to determine whether learning occurred as a result of repeated exposure to the hot plate.

Dissection

Brain tissue from the three-injection tolerance study was harvested from all four genotypes and treatment conditions (n = 12–16 mice/group) 3 h after the end of tolerance testing on the last test day. Each mouse was killed by CO2 asphyxiation followed by cervical dislocation and decapitation. A Plexiglas brain mold (David Kopf Instruments) was used to slice the fresh whole brain into coronal slices appropriate for each section desired (see below).

Using the atlas of Paxinos and Franklin (2001) as a reference, the following sections were taken: prefrontal cortex (PFC), interaurally 5.90–3.70 mm, from the dorsal-most point to 2.25 mm ventrally at a medial-lateral width of 1 mm centered on midline; ventral striatum (VS), interaurally 5.90–3.70 mm, from a dorsal-ventral position of 3.75 mm ventrally to just short of the ventral surface (not including olfactory tubercle) and from a medial-lateral position 2.0 mm to midline (this region included regions of interest such as the nucleus accumbens, ventral pallidum, bed nucleus of the stria terminalis, and substantia innominata); temporal lobe (TL), interaurally 2.74–1.74 mm from the dorsal-most point at 3.80 mm to the ventral surface from a medial point of 1.80 mm laterally to the lateral surface; periaqueductal gray (PAG) containing dorsal raphe, interaurally 0.00 to −1.16 mm from the dorsal point 1.80–3.00 mm ventrally at a medial-lateral width of 1.5 mm centered on midline.

The PAG has consistently been implicated in nociception, antinociception, and tolerance (Tortorici et al., 2001). In addition, the VS (Manning et al., 1994), TL (Pavlovic et al., 1996; Nandigama and Borszcz, 2003), and PFC (Hardy and Haigler, 1985) have also been shown to play a role in nociception and the analgesic effects of opioids (Chudler and Dong, 1995; Zubieta et al., 2001). In addition to nociception and analgesia, these areas were of interest in terms of their role in conditioned reinforcement and the reinforcing properties of drugs of abuse (VS, PFC) (Wise, 1996; Shippenberg and Elmer, 1998; Everitt et al., 1999).

Microarray gene expression profiling

Tissue from the four brain regions was harvested for gene expression analysis on a 27,648 element mouse cDNA microarray using a common reference design (Letwin et al., 2006; Lee and Saeed, 2007). The microarray was comprised of the National Institute on Aging Ko set of 15,247 cDNA clones (Tanaka et al., 2000) and National Institute of Mental Health Brain Molecular Anatomy Project set of 11,136 cDNA clones, together representing ∼21,000 distinct mouse genes. Each brain region from two to three subjects was pooled to obtain adequate tissue samples, thereby forming three to five independent replicates. Total RNA was prepared using Trizol reagent (Invitrogen) and the RNeasy mini kit (QIAGEN) as per manufacturers' instructions. Each pooled RNA sample was further divided into two equal aliquots to enable a technical replication known as a dye swap hybridization (Lee and Saeed, 2007). Thus, expression analysis of a brain region for a particular genotype and treatment was derived from three to five independent hybridizations and three to five corresponding dye swap hybridizations. For statistical analysis, a biological replicate experiment was defined as the average of an independent hybridization and its corresponding dye swap hybridization (i.e., five independent and five corresponding dye swap hybridizations yield five biological replicates). A hybridization experiment consisted of Cy5-labeled cDNA that was reverse transcribed from 15 μg of pooled total RNA and cohybridized with Cy3-labeled cDNA synthesized from an equal amount of the Stratagene Universal Mouse Reference RNA. Dye swap hybridizations were performed by reversing the dyes for each of the RNA samples. Hybridizations were performed for 18–24 h at 42°C followed by washing in decreasing concentrations of SSC at room temperature and spun dry. Microarray image scanning, fluorescence intensity measurements, LOWESS data normalization, normalization across replicate experiments, experimental noise determination, and cluster analysis were performed as described previously (Teramoto et al., 2005; Letwin et al., 2006).

Microarray gene expression analysis

Gene expression values from saline- and morphine-treated subjects in four genotypes in four brain regions can be defined as cgsr or mgsr, where c is saline control, m is morphine, g is gene, s is genotype, and r is region. These data yield two base data sets: (A) control gene expression without morphine intervention (cgsr), and (B) gene expression after morphine administration (mgsr). From these two data sets, we derive an additional data set (C), the change in gene expression relative to the mean value of the genotype and region control expression (mgsr − xcgsr, or Δ). Since the morphine expression data are important mostly in relation to the control levels, the Δ data set was emphasized over mgsr in future analyses and discussion. Statistical analysis for genotype-dependent differences in data sets cgsr and Δ was accomplished by an ANOVA with 10% false discovery rate (FDR) to control for type I error resulting from multiple comparisons (Reiner et al., 2003). The implementation of a 10% FDR corresponded to an effective α value of ≤0.00787 (data not shown). Subsequently, the significant genes were filtered for a minimum 1.5-fold change across genotypes based on “self versus self” hybridizations in which aliquots of the same RNA sample were labeled separately with Cy3 and Cy5 dyes and cohybridized onto a microarray. Results from these control hybridizations (data not shown) and previous published studies (Yang et al., 2002; Lee et al., 2007; Liang et al., 2008) demonstrate that a Z-score of ∼2 to 3 corresponds to a 1.4-fold “difference,” providing a metric for inherent assay noise. Significant genes displaying a ≤1.4-fold change in the strain comparisons exhibited a low validation success rate (∼50%) as assessed by quantitative real-time reverse transcription (RT)-PCR (data not shown), whereas genes above this cutoff displayed a validation rate (>90%) in line with the chosen FDR. We will refer to these statistically significant and filtered data sets as follows: (1) “genotype signature genes,” cgsr genes that are differentially expressed under control conditions as a function of genotype in each brain region, and (2) “morphine signature genes,” Δ genes that are significantly different as a function of genotype in each brain region.

Association of gene expression with tolerance across genotypes

The phenotypic endpoint used to define tolerance was set as [100 − (%MPE on day 3)] for each genotype. This number reflects the diminution in analgesic efficacy by the end of the chronic treatment regiment and assigns a higher value to the genotypes showing the most tolerance by the end of the treatment. A two-tier analysis strategy was used to define candidate genes associated with the tolerance phenotype. First, genes had to exhibit significant genotype-dependent expression differences (i.e., genotype signature and morphine signature genes). Second, the genotype signature or morphine signature genes must be significantly correlated (p < 0.05) with the tolerance phenotype as defined by pattern matching/feature selection (Pavlidis and Noble, 2001; Letwin et al., 2006). The subset of genotype signature genes that were significantly associated with tolerance are defined as candidate “predisposition” genes. The subset of morphine signature genes that were significantly associated with tolerance are defined as candidate “tolerance” genes.

Identification of regional conservation of enriched gene expression

Biological themes associated with morphine signature genes were identified by using the three Gene Ontology (GO) categories of molecular function, biological process, and cellular component in the Expression Analysis Systematic Explorer (EASE) application (Hosack et al., 2003), which is executable using TIGR Multi Experiment Viewer (TMEV) (available at www.tigr.org/softlab).

Identification of gene networks and canonical pathways

Identification of candidate genes that are highly interconnected or overrepresent a biological pathway may provide insight into molecular events that predispose or are causally related to the development of tolerance. To provide a functional framework to assist interpretation and narrow the candidate gene list to high value targets, gene networks and canonical pathways were explored using Ingenuity Pathway Analysis (IPA), Pathway Studio, the GO Consortium (www.geneontology.org), GenMAPP (www.genmapp.org), and Kyoto Encyclopedia of Genes and Genomes (www.genome.jp/kegg).

Ingenuity Systems is an online-based bioinformatics resource that algorithmically identifies gene networks and canonical pathways based on the gene lists. The gene lists uploaded into the database consisted of unique gene identifier (e.g., RefSeq, GenBank accession, Affymetrix ID), expression and correlation values. To identify gene networks, IPA rank orders the gene list in order of interconnectedness (triangle connectivity), and then adds genes to the top focus gene based on specific connectivity (how much the connections of a new gene overlap with the seed gene), and then combines smaller networks through the addition of linker genes found either in the uploaded list or IPA's Global Molecular Network. A p value, the number of focus molecules (f), and top functional potential is given for each network. The p value is defined as the probability of finding f or more Focus Genes in a set of n genes randomly selected from the Global Molecular Network database using a right-tailed Fisher's exact test. All significant biological themes associated with these networks were subsequently cross-validated using Pathway Studio, GenMAPP, GO, and KEGG pathways. To identify canonical pathways, IPA quantifies the likelihood that the number of uploaded genes that overlap with any one of 163 established canonical pathways occurs at a level greater than chance using a right-tailed Fisher's exact test. The p value and ratio value (number of molecules in a given pathway divided by total number of molecules that make up the pathway) are provided for each pathway.

Association of predisposition and tolerance genes with expression QTLs

Predisposition and tolerance genes in the PAG, PFC, TL, and VS were associated with expression QTLs (eQTLs) using the WebQTL module available at www.genenetwork.org/home.html. INIA Brain mRNA M430 (Jun06) RMA, VCU BXD PFC Sal M430 2.0 (Dec06) RMA, INIA Brain mRNA M430 (Jun06) RMA, HBP Rosen Striatum M430V2 (Apr05) RMA served as the databases for eQTL searches in the PAG, PFC, TL, and VS, respectively. Significant linkage of eQTL markers to correlated genes was defined by marker regression plots with 1000 permutations. The results from the permutation tests provide LRS (likelihood ratio statistics) scores that are suggestive, significant, or highly significant.

Real-time RT-PCR validation of genes

Expression levels of a portion of the genes were confirmed by real-time RT-PCR on an ABI Prism 7700 Sequence Detection System as previously described (Malek et al., 2002). Total RNA from PFC, TL, VS, and PAG in all genotypes was reverse transcribed using random primers following the manufacturer's instructions. The resulting cDNA was diluted and used as a template for quantitative PCR (qPCR). PCR primers were selected for specificity by the National Center for Biotechnology Information BLAST (basic local alignment search tool) of the mouse genome, and amplicon specificity was verified by first-derivative melting curve analysis using software provided by PerkinElmer and Applied Biosystems. Quantitation and normalization of relative gene expression were accomplished by using the comparative threshold cycle method as described previously (Joe et al., 2005). The following genes were identified by microarray analysis as differentially expressed across genotypes and subjected to RT-PCR validation: protein kinase Cγ (NM_011102), protein kinase A (NM_008854), cannabinoid receptor 1 (NM_007726), adenosine receptor 1 (NM_001008533), and dicer1 (NM_148948). The expression of the “housekeeping” genes carnitine acetyltransferase (NM_007760), peptidylprolyl isomerase-like 4 (NM_026141), and succinate-CoA ligase α subunit (NM_019879) were used for normalization as these genes did not exhibit differential expression in our microarray assays. Sequence of primers used for the real-time RT-PCR is given in supplemental Table 1 (available at www.jneurosci.org as supplemental material).

Results

Behavioral testing: nociception and analgesia

Introgression of the Oprm region did not alter baseline sensitivity or analgesic effects of morphine

Baseline sensitivity.

There was a significant main effect of genotype [F(genotype) df(3,47) = 8.67; p < 0.0001] on baseline nociception (Fig. 1A). The B6 and B6D2 genotypes were significantly more sensitive than the D2 and D2B6 mice. The B6 and B6D2 genotypes did not differ from each other nor did the D2 and D2B6 genotypes. Therefore, baseline pain sensitivity was not affected by the introgression of the proximal region of Chr10.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Baseline and acute analgesic effects of morphine. A, Baseline nociception. B, Dose–response curve at 0, 1, 5.6, 10, 30 mg/kg of morphine. Results are the mean ± SE of seven to nine animals per genotype. *p < 0.05.

Acute analgesic sensitivity.

There was a significant main effect of genotype on acute sensitivity to the analgesic effects of morphine [F(genotype) df(3,23) = 35.40; p < 0.0001]; the B6 and B6D2 genotypes were significantly less sensitive than the D2 and D2B6 mice (Fig. 1B). ED50 values for morphine were 4.2 and 5.1 mg/kg for the D2 and D2B6 mice and 11.8 and 12.2 mg/kg for the D2 and D2B6 mice, respectively. As was true for baseline nociception, introgression of the proximal Chr10 region did not alter acute sensitivity to the analgesic effects of morphine; the B6 and B6D2 genotypes did not differ from each other nor did the D2 and D2B6 genotypes. This was a surprising finding considering bQTL analyses have mapped the proximal chromosome region of 10 to morphine analgesia (Belknap et al., 1995) and other morphine-related behaviors (Berrettini et al., 1994a; Ferraro et al., 2005).

Reciprocal introgression of the Oprm region altered the development of tolerance

Figure 2A (three-injection experiment) shows the analgesic effects of morphine given every other day for 5 d as a percentage of analgesia produced by the first injection. There was no significance difference in the analgesic effects of morphine after the first injection; therefore, tolerance assessment was performed at approximately equianalgesic doses. Genotype significantly influenced the development of tolerance [F(genotype) df(3,50) = 17.25; p < 0.0005]. Tolerance developed most rapidly in the B6 followed by the B6D2 strain. Tolerance did not develop in the D2 mice and only slightly in the D2B6 animals. Pairwise comparisons representing statistical significance after the last injection are delineated by the arrows.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Tolerance to the analgesic effects of morphine. A, When the ED90 dose for each genotype is given once every other day for 5 d (mean ± SE of 12–16 animals per genotype). B, When given once a day for 5 d (mean ± SE of 6–9 animals per genotype). The average of the last 3 d is given to the right of the graph for comparison. Arrows indicate significant pairwise comparisons.

Figure 2B (five injection, phenotype confirmation) shows the analgesic effects of morphine given everyday for 5 d as a percentage of that produced by the first injection. Again, there was no significant difference in the analgesic effects of morphine after the first injection; therefore, tolerance assessment was performed at approximately equianalgesic doses. Genotype significantly influenced the development of tolerance [F(genotype) df(3,23) = 19.74; p < 0.0001]. Tolerance developed most rapidly in the B6 followed by the B6D2 and D2B6 mice. Tolerance did not develop in the D2 mice. Pairwise comparisons (genotype by day) for statistical significance after the last injection are delineated by the arrows. Overall, the shorter dose interval increased the rate of tolerance that developed after three injections [F(treatment schedule) df(1,71) = 7.19; p < 0.009] as did the additional injections of morphine [F(treatment schedule) df(1,71) = 14.27; p < 0.0003] (Fig. 2, compare A, B). However, regardless of the test paradigm, the D2 genotype did not develop tolerance to morphine and there was no genotype by treatment schedule interaction (p < 0.31; NS); decreasing the interinjection interval and increasing the number of injections produced a parallel change in the development of tolerance across genotypes. Thus, the five-injection phenotype confirmation experiment confirms the results of the three-injection schedule and extends the generality to schedules using a shorter injection interval and more injections.

The degree of tolerance was not influenced by learning factors associated with repeated hot-plate testing. There was no difference between the experimental groups exposed to a heated plate versus those exposed to room temperature plate in the B6D2 or B6 genotypes (p < 0.86, p < 0.34, respectively; both NS) (supplemental Fig. 2, available at www.jneurosci.org as supplemental material).

Microarray gene expression analysis: genotype-dependent differences

The expression profiles of ∼21,000 cgsr genes were measured in each of four brain regions (PAG, PFC, TL, VS) across four saline-treated strains (B6, D2, B6D2, D2B6), and a one-way ANOVA (main factor genotype) at 10% FDR was performed per brain region. Together, a total of 1070, 2550, 1758, and 2684 genotype signature genes were identified (cgsr genes exhibiting significant strain-dependent differences in expression of at least 1.5-fold or greater) in the PAG, PFC, TL, and VS, respectively (Fig. 3A; supplemental Table 2, available at www.jneurosci.org as supplemental material), representing 4562 unique genes across the four brain regions. To define genotype-dependent differences in morphine-induced changes in gene expression, the Δ value for each gene was computed from the morphine and saline treatment data and subjected to a one-way ANOVA at 10% FDR per brain region. A total of 1182, 1642, 2041, and 2074 morphine signature genes were detected in the PAG, PFC, TL, VS, respectively (Fig. 3B; supplemental Table 2, available at www.jneurosci.org as supplemental material), representing 3942 unique genes across the four brain regions. The number of unique morphine signature genes is not unexpected given the well documented effects of morphine on the expression and/or activity of transcription factors such as ΔFosB and CREB (Han et al., 2006; Zachariou et al., 2006). The maximum Δ/minimum Δ ratio for each gene across the four genotypes ranged from 1.5- to 6-, 1.5- to 8-, 1.5- to 11-, and 1.5- to 10-fold in the PAG, PFC, TL, and VS, respectively (data not shown).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Number of significant gene expression differences across four genotypes in four brain regions. A, Genotype signature genes (significant cgsr genes in saline-treated animals). B, Morphine signature genes (significant mgsr − xcgsr genes) (see Materials and Methods section for explanation). Statistical significance was defined by ANOVA per brain region with 10% FDR.

Biological themes associated with morphine signature genes

Morphine signature genes (Δ genes) were subjected to EASE analysis to identify potential biological themes associated with the four brain regions. Genes were classified by using GO categories belonging to the ontologies of molecular function, biological process, and cellular component. A total of 73, 166, 139, and 144 GO categories were significantly overrepresented (Fisher's exact test, p < 0.05) by morphine signature genes in the PAG, PFC, TL, and VS, respectively. Particularly noteworthy were the identification of epigenetic-related categories such as chromatin remodeling complex (GO:0016585) in the PFC, RNA interference (GO:0016246) and regulation of cell differentiation (GO:0045595) in the TL, and plasticity-related categories such as neuron development (GO:0048666) and cell projection organization and biogenesis (GO:0030030) in the VS.

Association of gene expression with tolerance across genotypes

Gene expression associations with tolerance in different brain regions

An associative analysis was performed between signature genes and the behavioral endpoint tolerance for candidate gene discovery. Genotype signature genes that correlate with tolerance are postulated to predispose animals to morphine-induced tolerance (such gene associations are termed “predisposition genes”). A similar approach has been used to identify “naive” genes predisposing inbred mice to ethanol-related behaviors (Letwin et al., 2006). The tolerance value [100 − (%MPE on day 3)] for each genotype was used as the template to search for genotype signature genes (Fig. 3A); the expression value of a gene across the four genotypes was required to match the tolerance template across the four genotypes to be identified. This was done in each of the four brain regions. Using pattern matching as a correlative search tool (also known as feature selection) (Pavlidis and Noble, 2001), 64, 53, 86, and 123 predisposition genes (positive and negative) were identified in the PAG, PFC, TL, and VS, respectively (pattern match p value, <0.05) (Fig. 4A; supplemental Table 3, available at www.jneurosci.org as supplemental material). Tspyl4 and 6330407J23Rik were the only predisposition genes physically located within the congenic region demarcated by markers D10Mit75 and D10Mit61. Using the GeneNetwork eQTL resource (http://www.genenetwork.org), Tspyl4 was found to be regulated in trans by an eQTL locus associated with marker rs3682996 on chromosome 1, whereas 6330407J23Rik was determined to be regulated in cis (supplemental Table 3, available at www.jneurosci.org as supplemental material).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Gene expression correlations with tolerance in different brain regions. A, Saline correlations in naive animals. B, Δ correlations in morphine-treated animals. Cluster images of template matched genes with positive (top panel) and negative correlations (bottom panel). Highly expressed genes are indicated in red, intermediate in black, and weakly expressed in green. Mouse strains are indicated above each column. The black graph above each cluster represents tolerance template values, whereas the red graph represents the mean expression value of all matched genes within each column. Data presented are the mean of 4–12 hybridizations.

Next, the tolerance value for each genotype was used as the template to search for “tolerance genes” (Fig. 3B) within the morphine signature data set; again, the expression value of a gene across the four genotypes was required to match the tolerance template across the four genotypes to be identified as a candidate gene. Using feature selection, 81, 96, 106, and 82 tolerance genes (positive and negative) were identified in the PAG, PFC, TL, and VS, respectively (pattern match p value, <0.05) (Fig. 4B; supplemental Table 3, available at www.jneurosci.org as supplemental material). These genes are postulated to represent potential mediators of tolerance whose expression is modulated by morphine.

We also performed a two-way ANOVA (genotype and treatment serving as the main factors) with 10% FDR correction in each brain region and subjected the resulting significant genotype by treatment interactions to pattern matching to identify the tolerance genes. Despite the more conservative nature of the two-way ANOVA approach compared with the one-way ANOVA strategy used above, a substantial overlap of pattern-matched genes was found between the two approaches. For example, 79% of the PFC tolerance genes originally identified by the one-way ANOVA approach were also captured in the two-way ANOVA approach. Overall, the overlap between the two approaches across the four brain regions ranged from 77% in the PAG to 80% in the VS. Given the paucity of information pertaining to genes associated with tolerance behavior, all subsequent analyses will be based on the results from the one-way ANOVA.

A phenotype–gene association (PGA) clustergram was generated by hierarchical clustering of the 369 tolerance genes across the four brain regions (Fig. 5). This plot is analogous to hierarchical clustering of gene expression values typically seen in microarray experiments (Eisen et al., 1998); however, rather than clustering one-dimensional expression values, PGA analysis clusters biological information pertaining to the direction and significance level of gene–phenotype correlations (Reiner-Benaim et al., 2007). Tolerance genes exhibiting significant positive correlations are indicated in red, significant negative correlations in green, and nonsignificant correlations in black. Noteworthy in Figure 5 is the paucity of genes that are correlated in more than one brain region (<2%). A similar finding was noted for the predisposition set (<2%) (supplemental Table 3, available at www.jneurosci.org as supplemental material). These observations suggest that networks comprised of different combinations of predisposition and/or tolerance genes in each brain region may be working in concert to promote tolerance susceptibility and acquisition. Tolerance genes could be grouped into a number of interesting categories such as gene expression, second messenger-mediated signaling, neuronal synaptic plasticity, response to oxidative stress, signal transduction/G-protein-coupled receptor and transport (Fig. 5).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

PGA plot. Cluster image represents morphine signature genes (by row) and the strength of its correlation (color coded) with the behavioral phenotype tolerance in each of the four brain regions (column). Positively correlated morphine signature genes are indicated in red, noncorrelated in black, and negatively correlated in green.

In silico mapping of tolerance and predisposition genes to behavioral and expression QTLs

The mapping of correlated genes to behavioral intervals provides additional supportive evidence when associating a gene to a behavior. Physical map locations of our correlated genes were obtained from ENSEMBL (http://www.ensembl.org) and morphine-related bQTL information was acquired from the literature. Provisional bQTLs mapping to chromosomes 1, 6, and 10 for morphine preference (Berrettini et al., 1994a), chromosomes 3, 7, 8, 9, 10, 12, and 18 for morphine consumption (Belknap and Crabbe, 1992), chromosomes 2, 3, 4, 7, 9, 10, 11, and 19 for morphine analgesia (Belknap and Crabbe, 1992; Belknap et al., 1995; Bergeson et al., 2001), and chromosomes 1, 5, and 10 for morphine withdrawal (Kest et al., 2004) have been identified through crosses of the B6 and D2 strains, representing ∼15% of the mouse genome. For our analysis, a bQTL interval was defined as two logarithm of odds units around the linkage peak (Rapp, 2000), provided that information on multiple genetic markers in the vicinity of the peak was given (Berrettini et al., 1994a). Otherwise, the interval was defined by 10–15 Mb on either side of a lone genetic marker defining the bQTL. Of the 369 tolerance genes, 316 had unambiguous mapping information, and of these mapped genes 73 (23%) were found to reside in bQTLs that influence morphine-related behaviors (supplemental Table 3, available at www.jneurosci.org as supplemental material). Similar results were obtained in predisposition gene set in which 70 of 258 mapped genes (27%) were found to reside in bQTLs (supplemental Table 3, available at www.jneurosci.org as supplemental material). To ascertain whether or not the mapping of correlated genes in bQTLs may have occurred by chance, we compared the number of correlated and noncorrelated genes residing in bQTLs for each brain region. There was a significant occurrence not caused by chance alone of tolerance genes mapping to bQTLs for the TL (χ2, p = 0.0069) and PAG (χ2, p = 0.0216); and predisposition genes mapping to bQTLs for the VS (χ2, p = 0.0059) and PFC (χ2, p = 0.0221). A number of interesting tolerance genes residing in bQTLs associated with morphine analgesia included the glycine receptor β subunit (positively correlated in PAG), signal transducer and activator of transcription 5 (positively correlated in PFC), ADP-ribosylation factor interacting protein 1 (negatively correlated in TL), and tumor suppressor candidate 2 (negatively correlated in VS).

Of the tolerance and predisposition genes residing inside bQTLs but outside of the congenic region, 13 and 6 genes, respectively, were regulated in trans by eQTLs located in the congenic region (supplemental Table 3, available at www.jneurosci.org as supplemental material). Moreover, 39 tolerance and 25 predisposition genes residing outside bQTLs were found to be regulated in trans by eQTLs in the congenic region (supplemental Table 3, available at www.jneurosci.org as supplemental material). There was a significant difference in the occurrence of tolerance versus noncorrelated genes mapping inside eQTLs for the TL (χ2, p = 0.0469), VS (χ2, p = 0.0381), and PAG (χ2, p = 0.0053); and predisposition versus noncorrelated genes for the TL (χ2, p = 0.0392) and PAG (χ2, p = 0.0032). Notable trans-regulated examples included synaptic vesicle membrane proteins (synaptotagmin XI, syntaxin 3), transcription factors (Stat5a) and extracellular matrix/cell adhesion molecules (laminin B1, integrin β3 binding protein).

Identification of gene networks and canonical pathways

To address whether the correlated genes are potentially interlinked into a gene regulatory network, all predisposition and tolerance genes were computationally analyzed using microarray bioinformatics resources available through Ingenuity Pathway Analysis, Pathway Studio, the GO Consortium (www.geneontology.org), GenMAPP (www.genmapp.org), and Kyoto Encyclopedia of Genes and Genomes (www.genome.jp/kegg). Three network “types” can be generated in a region-specific manner: a network of genes based solely from the predisposition set, a network derived solely from the tolerance set, and a network derived from the combined consideration of predisposition and tolerance genes. The top three networks in each network type for each brain region are provided in supplemental Table 4 (available at www.jneurosci.org as supplemental material). Each network type affords a slightly different viewpoint. The predisposition network may interact in such a way as to provide a genotype-dependent starting “set point” that predisposes the subject to tolerance. These genes may not necessarily change expression in response to morphine but set the stage for the effects of morphine on other genes. The tolerance networks may interact in such a way as to cause (or prevent via compensatory mechanisms) tolerance. These genes do not necessarily differ at baseline; rather they change expression in response to morphine administration in a genotype-dependent manner. The combined consideration networks define potential preexisting interactions of predisposition and tolerance genes during the acquisition of tolerance. Promoter analysis of the three network types using the TRANSFAC motif database and MAST (Malek et al., 2006; Lee et al., 2007) did not identify any statistically significant overrepresented transcription factor binding motifs (data not shown).

Figure 6 (right) shows an example network derived from the tolerance gene list in the PAG. Nineteen of the 34 genes in the PAG tolerance network were either positively (12) or negatively (7) associated with the development of tolerance. Eleven of the tolerance genes reside in bQTLs previously associated with morphine-related behaviors (supplemental Table 3, available at www.jneurosci.org as supplemental material). The probability of having 11 of 19 network tolerance genes residing in morphine-related bQTLs was significantly greater than would be expected by chance alone (χ2, p < 0.0001). Despite having <2% overlap overall between the predisposition and tolerance genes (across all brain regions), seven genes in the PGA tolerance network had corresponding predisposition gene counterparts. As a means to consider the potential alterations in network associations that may occur after chronic drug treatment, Figure 6 (left) shows a simplified network comprised of the seven predisposition gene counterparts (H19, Ttk, Hba2, B4galt1, Hmox1, Casp3, cFos). A comparison of Figure 6 (left and right) was made to demonstrate how prevailing gene networks may undergo alterations in network associations after chronic drug treatment. Some predisposition gene counterparts undergo a reversal in direction of their association with tolerance after morphine treatment (H19, Ttk, Hba2, B4galt1). For instance, Ttk was positively correlated in saline-treated animals but negatively correlated in animals injected with morphine. Hba2 in the same network was negatively correlated in saline-treated animals but oppositely correlated on treatment with morphine. Overall, the genes modified by chronic morphine administration show a strong presence in canonical pathways representative of neuroadaptation (see next section) and illustrate potential alterations in network associations after chronic drug administration.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Tolerance gene network in the PAG of saline-treated animals (left) and morphine-treated animals (right). A gene regulatory network with lipid metabolism as the associated biological theme was constructed from the tolerance gene set (correlated morphine signature genes). The correlation status of the same genes in saline-treated animals provides a reference point for the transitioning of drug-naive animals to morphine tolerant animals. Edges (lines) connecting nodes (genes) represent regulatory interactions such as regulation of gene expression, protein–protein interactions, protein–nucleic acid interactions, protein–hormone interactions. Positively correlated genes are indicated in red, noncorrelated in black, and negatively correlated in green. See supplemental Tables 4 and 5 (available at www.jneurosci.org as supplemental material) for definitions of gene symbol abbreviations.

A complementary analysis to the de novo network construction using candidate genes is to determine their representation in molecular pathways with known function. Supplemental Table 5 (available at www.jneurosci.org as supplemental material) presents the top 15 canonical pathways associated with genotype signature, morphine signature, predisposition, and tolerance genes in each brain region along with the number of genes found in each pathway and its corresponding p value. As noted previously, the pathways identified using the genotype and morphine signature lists show a strong presence in canonical pathways representative of neuroadaptation. In this regard, synaptic long-term potentiation (LTP) was identified in the top five genotype and morphine signature lists in the PAG and TL, and PAG and PFC (supplemental Fig. 3, available at www.jneurosci.org as supplemental material). Axonal guidance signaling was identified in the top five canonical pathways in all brain regions within genotype signature list and in four of the five brain regions in the morphine signature list. Additional pathways within the morphine signature list suggesting a major role of synaptic plasticity include the ephrin receptor, Huntington's disease, and neuregulin signaling. The smaller number of genes in the predisposition and tolerance lists limit the statistical probability that they will populate a canonical pathway. Nevertheless, there were several pathways in the tolerance list suggestive of neural plasticity [neuregulin, notch, and ERK (extracellular signal-regulated kinase)/MAPK (mitogen-activated protein kinase) signaling] and cell death [NRF-2 mediated oxidative stress response, GM-CSF (granulocyte-macrophage colony-stimulating factor), and p53 signaling].

Real-time RT-PCR

Five tolerance genes (protein kinase Cγ, protein kinase A α subunit, cannabinoid receptor 1, adenosine A1 receptor, dicer-1 protein) involved in six correlations were selected for microarray validation. These genes exhibited the smallest fold change (1.5- to 2-fold) of maximum Δ/minimum Δ across the four strains, presumably representing the lower limits of our analysis to define significant correlations. Protein kinase Cγ, protein kinase A catalytic α subunit, cannabinoid receptor 1, and adenosine A1 receptor have been implicated in morphine-related behaviors including tolerance, albeit these studies are typically performed in a single mouse strain (Kaplan and Leite-Morris, 1997; Narita et al., 2001; Smith et al., 2006; Trang et al., 2007). It should be noted that none of the aforementioned genes had, before this study, been correlated with tolerance across multiple genetic strains. The fifth gene chosen for validation, dicer-1, has not previously been linked to drug tolerance. Finally, protein kinase A catalytic α subunit and adenosine A1 receptor genes were of interest because they mapped to bQTLs associated with morphine consumption and preference, respectively (supplemental Table 3, available at www.jneurosci.org as supplemental material). Quantitative real-time RT-PCR validation of the five genes in six correlations is shown in Figure 7. A significant correlation was found between the microarray and real-time RT-PCR results (p < 0.05) and between the RT-PCR results and tolerance behavior for all genes, supporting the reliability of our gene expression measurements.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Real time RT-PCR validation of microarray data. Validation of individual template matched genes plotted against tolerance (left panels) and microarray data (right panels). Prkcc VS, Protein kinase Cγ in ventral striatum (correlation with tolerance, r = −0.99, p = 0.0042; correlation with microarray data, r = 0.99, p = 0.0019); Prkaca PFC, protein kinase A catalytic α subunit in prefrontal cortex (correlation with tolerance, r = −0.99, p = 0.0033; correlation with microarray data, r = 0.97, p = 0.028); Prkaca TL, protein kinase A catalytic α subunit in temporal lobe (correlation with tolerance, r = −0.96, p = 0.038; correlation with microarray data, r = 0.99, p = 0.013); Cnr1 PAG, cannabinoid receptor 1 in periaqueductal gray (correlation with tolerance, r = 0.99, p = 0.0023; correlation with microarray data, r = 0.98, p = 0.025); Adora1 TL, adenosine A1 receptor in temporal lobe (correlation with tolerance, r = −0.96, p = 0.041; correlation with microarray data, r = 0.93, p = 0.049); Dicer1 PFC, dicer 1 in prefrontal cortex (correlation with tolerance, r = 0.93, p = 0.041; correlation with microarray data, r = 0.98, p = 0.024). Data are represented as the mean ± SE. •, B6; ♦, B6D2 congenic; ▴, D2B6 congenic; ■, D2.

Discussion

One of our goals was to investigate the role of the proximal region of Chr10 in the acute and chronic effects of morphine. Previous bQTL studies suggest that this region between markers D10Mit28-D10Mit3 (4–21 cM) accounts for a significant amount of variance in morphine two-bottle choice preference and analgesia (Berrettini et al., 1994a,b; Belknap et al., 1995; Bergeson et al., 2001) (supplemental Fig. 1, available at www.jneurosci.org as supplemental material). A comparison of the B6 and D2 strains in the Mouse Genome Database (http://www.informatics.jax.org) reveals the presence of 21 nonsynonymous coding and 207 noncoding single-nucleotide polymorphisms (SNPs) between D10Mit28 and D10Mit3, which were thought to account for the observed variance in behavioral responses. Oprm, the gene responsible for the primary receptor target of morphine (μOR), is located within this region at 8.0 cM. Reciprocal introgression of this promising region, however, did not alter baseline nociception or the acute analgesic effects of morphine in our study. These results are in contrast to those found by Bergeson et al. (2001) using a similar congenic approach in which the Chr10 region between markers D10Mit28 and D10Mit3 (4–21 cM) from the donor B6 strain was introgressed into the D2 background (supplemental Fig. 1, available at www.jneurosci.org as supplemental material). In that case, acute sensitivity to 10 mg/kg morphine (intraperitoneal) was altered in the reciprocal congenic strains; baseline nociception was not reported nor was tolerance assessed. Differences in stimulus intensity, route of administration, endpoint, or alterations in baseline sensitivity (Elmer et al., 1998) could contribute to the discrepancy. Equally important is the fact that the introgressed regions differ substantially in size between our congenics and the ones studied by Bergeson et al. (2001). Specifically, our D2B6 strain carries an additional 11 cM portion or 37.4 Mb of Chr10 derived from the B6 strain (supplemental Fig. 1B, available at www.jneurosci.org as supplemental material). Consequently, genetic interactions (i.e., networks) likely diverge between the two D2B6 lines because of the presence of 162 SNPs residing within 73 genes (many encoding transcription factors) found in this 11 cM region. In fact, 16 of the 162 SNPs have been identified as eQTLs predicted to trans-regulate a subset of our tolerance and predisposition genes. The consequences of different introgressed genetic footprints are analogous to the “congenic footprint” phenomenon, which can lead to profound differences in gene expression and behavioral responses (Lee et al., 2007; Schalkwyk et al., 2007). Also noteworthy is our finding that, although reciprocal introgression was without effect on acute morphine analgesia, presumably because of the genetic footprint, profound effects were observed for chronic morphine analgesia. Although differences in environmental context and protocol may ultimately explain a large part of the differences between studies, the difference in introgressed footprint may fortuitously have explanatory significance worth pursuing.

Tolerance to the analgesic effects of morphine develops at a rate directly related to the magnitude of the initial pharmacological effect (Cox and Tiffany, 1997). If the same dose of morphine is used to induce tolerance in subjects that differ significantly in the acute analgesic potency of morphine, tolerance could appear to develop more rapidly in sensitive genotypes given the greater magnitude of the initial pharmacological effect. This confound would affect the tolerance outcome and interpretation. The current study accounts for this by using genotype-specific doses within a “therapeutic” range (ED90). Under these conditions, reciprocal introgression of proximal Chr10 altered the development of tolerance. Thus, the acute analgesic effects of morphine are determined in part by genetically separable mechanisms from those that are involved after chronic administration. This conclusion is supported by previous behavior genetic (Oliverio and Castellano, 1974; Kest et al., 2002; Liang et al., 2006), molecular, and pharmacological investigations (Kaneto et al., 1985; Ocaña et al., 2007). A potential confound to this approach is that the genotype-specific ED90 dose is sufficient to produce analgesia acutely in all genotypes but does not reach a threshold for inducing the molecular mechanisms involved in tolerance in all genotypes. If the morphine dose was increased sufficiently to produce tolerance in the D2 mice (i.e., ∼10-fold dose increase produces tolerance in D2) (Kest et al., 2002), a different set of genes may be revealed to produce tolerance. However, this approach also runs the risk of confound by nonspecific neuropharmacological mechanisms engaged at supramaximal dosing.

Given the strong bQTL evidence, a straightforward explanation for differential response to chronic morphine in the B6 and D2 mice would be found in the variations in Oprm expression or protein regulation. In this regard, B6 and D2 mice differ qualitatively in VS μOR regulation after chronic high-dose morphine administration; μOR concentration decreases (−50%) in the B6 mice and increases (+20%) in the D2 mice (Petruzzi et al., 1997). There are two promoter regions in the Oprm gene, one upstream (5′ of exon 11) and one downstream (5′ of exon 1) (Liang and Carr, 1997; Pan, 2002). Five SNPs have been identified within a 1.7 kb region of the Oprm downstream promoter region between the B6 and D2 strains; three were near or within potential transcription binding sites (Doyle et al., 2006). A difference in the basal Oprm promoter activity was seen in vitro between the B6 and D2 constructs, suggesting a role of the polymorphisms in gene transcription (Doyle et al., 2006). However, our current study and those using B6 and D2 constructs in μOR-positive BE(2) cell lines (Doyle et al., 2006) suggest other factors are involved in tolerance in addition to Oprm expression since μOR was not differentially expressed in the cell lines, nor did we find differential expression of the μOR or differential expression of genes known to promote or repress Oprm (Law et al., 2004; Contet et al., 2008). Given the complex and oftentimes controversial evidence surrounding μOR regulation in tolerance, systematic large-scale gene expression analysis may yield alternative molecular mechanisms not found by investigating prototypical candidate mechanisms.

Large-scale gene expression studies were performed in conjunction with the behavioral studies to identify candidate genes associated with tolerance. A stumbling block in expression studies has been the lack of a systematic method to iteratively focus on the most promising candidates. We used canonical pathway and network analyses to help identify high-value candidate predisposition and tolerance genes. Canonical pathway analysis of the morphine-responsive genes suggests a significant role for neuroadaptive processes in response to chronic morphine administration (i.e., LTP, axonal guidance, ephrin, neuregulin pathways) (for LTP example, see supplemental Fig. 3, available at www.jneurosci.org as supplemental material). Network analysis was used to identify highly interconnected genes to provide insight into molecular events that predispose or are causally related to the development of tolerance. It is possible to analyze and view gene networks from several different vantage points. One viewpoint would be to determine network associations comprised of tolerance genes (Fig. 6, right panel) and determine the correlation status of the same genes before morphine administration (Fig. 6, left panel). An alternative permutation would consist of predisposition genes (supplemental Fig. 3, left panel, available at www.jneurosci.org as supplemental material) and their status after chronic morphine administration (supplemental Fig. 3, right panel, available at www.jneurosci.org as supplemental material). Changes in the correlation status of genes transitioning from baseline to morphine treatment in the network permutations may be of significant interest. This is nicely exemplified by the top network involving tolerance genes in the PAG. Nineteen of the 34 genes in this network were either positively (12) or negatively (7) associated with the development of tolerance. Protein kinase C, cannabinoid receptor 1, and calmodulin are known to play a prominent role in the current theories of morphine tolerance and neuroadaptation (Law et al., 2000; Williams et al., 2001; Trang et al., 2007; McClung and Nestler, 2008). Of the 19 tolerance genes contributing to this network, more than one-half (11) remarkably resided in bQTLs previously associated with morphine-related behaviors. The occurrence of these correlated genes in bQTLs was significantly greater than would be expected by chance alone, presumably reflecting their importance in tolerance acquisition. Before morphine administration, 7 of the 34 genes making up this network were found to be associated with tolerance. Three of these genes (c-Fos, Casp, Pld) were not correlated after tolerance developed which could be traced to a parallel morphine response (downregulation) across the genotypes. Strikingly, four genes (B4galt1, Hba2, Ttk, H19) exhibited a reversal in the direction of correlation after chronic morphine administration. For example, high H19 expression is predictive of tolerance development (direction of predisposition association), and a negative change in expression (i.e., downregulation) is associated with tolerance (direction of tolerance association). H19 expression appears to decrease in the genotypes that become tolerant to morphine and increase in those that do not.

Several recent studies provide intriguing insight into the potential role of H19 in response to chronic morphine administration. First, it has been suggested that H19 is a primary micro-RNA precursor (to miR-675) (Cai and Cullen, 2007). Coincidently, our data reveal that Dicer1, encoding the pre-micro-RNA processing enzyme (Bernstein et al., 2001), is associated with the development of tolerance in the PFC, and EASE analysis of morphine-responsive genes identified the RNA interference pathway as playing a role in the TL. Second, analysis of 12,000 genes and ESTs (expressed sequence tags) between control undifferentiated and dedifferentiation derived cell clones establish only one gene as differentially regulated, H19 (Scott et al., 2005). In a separate study, Ayesh et al. (2002) identified 48 differentially expressed genes modulated by H19 in a manner that may promote metastasis; interestingly, we found that morphine administration altered the expression of 50% of these genes, but not all were directly associated with the development of tolerance. Third, the H19 imprinting region has been implicated in dopamine neuron differentiation (Freed et al., 2008) and gene regulation via interchromosomal interactions involving the CCCTC-binding factor (Ling et al., 2006). Interestingly, our findings reveal the CCCTC-binding factor to be a genotype signature gene as well. Moreover, five dopamine-related genes (dopamine D3 and D4 receptor, dopamine transporter, DARPP-32, and calcyon) were all found to be differentially regulated by morphine in this study. Although none of these genes were found to be correlated with analgesic tolerance, they may be involved in traits such as the rewarding effects of morphine. Overall, our findings suggest that a form of transdifferentiation or metaplasia-like activity may be an important component in the neuroadaptive response to chronic morphine. Thus, although this imprinting region has been recognized as important in neurodevelopment and cellular differentiation, our study suggests that it may play a significant role in response to chronic administration of pharmacologically relevant doses of morphine and warrants future investigation.

In summary, epigenetic and micro-RNA mechanisms and neuroadaptation appear to play a significant role in the nervous systems response to chronic morphine administration as corroborated by converging lines of evidence from EASE/GO, network, and canonical pathway analyses. Specific components of the neuroadaptive response are likely related to specific behavioral endpoints.

Footnotes

  • This work was funded by National Institute on Drug Abuse–National Institutes of Health Grant DA015087 (G.I.E., N.H.L.). We also acknowledge Dr. Wade H. Berrettini and Thomas N. Ferraro for their valuable contribution of mice and consultation.

  • Correspondence should be addressed to either of the following: Greg I. Elmer, Maryland Psychiatric Research Center, Department of Psychiatry, University of Maryland, Baltimore, Maple and Locust Street, Baltimore, MD 21228, gelmer{at}mprc.umaryland.edu; or Norman H. Lee, Department of Pharmacology and Physiology, George Washington University Medical Center, 2300 I Street, Washington, DC 20037, E-mail: phmnhl{at}gwumc.edu

References

  1. ↵
    1. Ayesh S,
    2. Matouk I,
    3. Schneider T,
    4. Ohana P,
    5. Laster M,
    6. Al-Sharef W,
    7. De-Groot N,
    8. Hochberg A
    (2002) Possible physiological role of H19 RNA. Mol Carcinog 35:63–74.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Belknap JK,
    2. Crabbe JC
    (1992) Chromosome mapping of gene loci affecting morphine and amphetamine responses in BXD recombinant inbred mice. Ann N Y Acad Sci 654:311–323.
    OpenUrlPubMed
  3. ↵
    1. Belknap JK,
    2. Mogil JS,
    3. Helms ML,
    4. Richards SP,
    5. O'Toole LA,
    6. Bergeson SE,
    7. Buck KJ
    (1995) Localization to chromosome 10 of a locus influencing morphine analgesia in crosses derived from C57BL/6 and DBA/2 strains. Life Sci 57:PL117–PL124.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bergeson SE,
    2. Helms ML,
    3. O'Toole LA,
    4. Jarvis MW,
    5. Hain HS,
    6. Mogil JS,
    7. Belknap JK
    (2001) Quantitative trait loci influencing morphine antinociception in four mapping populations. Mamm Genome 12:546–553.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bernstein E,
    2. Caudy AA,
    3. Hammond SM,
    4. Hannon GJ
    (2001) Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature 409:363–366.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Berrettini WH,
    2. Ferraro TN,
    3. Alexander RC,
    4. Buchberg AM,
    5. Vogel WH
    (1994a) Quantitative trait loci mapping of three loci controlling morphine preference using inbred mouse strains. Nat Genet 7:54–58.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Berrettini WH,
    2. Alexander R,
    3. Ferraro TN,
    4. Vogel WH
    (1994b) A study of oral morphine preference in inbred mouse strains. Psychiatr Genet 4:81–86.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Cai X,
    2. Cullen BR
    (2007) The imprinted H19 noncoding RNA is a primary microRNA precursor. RNA 13:313–316.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Chudler EH,
    2. Dong WK
    (1995) The role of the basal ganglia in nociception and pain. Pain 60:3–38.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Contet C,
    2. Filliol D,
    3. Matifas A,
    4. Kieffer BL
    (2008) Morphine-induced analgesic tolerance, locomotor sensitization and physical dependence do not require modification of mu opioid receptor, cdk5 and adenylate cyclase activity. Neuropharmacology 54:475–486.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Cox LS,
    2. Tiffany ST
    (1997) Associative and nonassociative tolerance: the effects of dose and interdose interval. Pharmacol Biochem Behav 57:31–36.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Doyle GA,
    2. Sheng XR,
    3. Schwebel CL,
    4. Ferraro TN,
    5. Berrettini WH,
    6. Buono RJ
    (2006) Identification and functional significance of polymorphisms in the mu-opioid receptor gene (Oprm) promoter of C57BL/6 and DBA/2 mice. Neurosci Res 55:244–254.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Eisen MB,
    2. Spellman PT,
    3. Brown PO,
    4. Botstein D
    (1998) Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 95:14863–14868.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Elmer GI,
    2. Pieper JO,
    3. Negus SS,
    4. Woods JH
    (1998) Genetic variance in nociception and its relationship to the potency of morphine-induced analgesia in thermal and chemical tests. Pain 75:129–140.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Everitt BJ,
    2. Parkinson JA,
    3. Olmstead MC,
    4. Arroyo M,
    5. Robledo P,
    6. Robbins TW
    (1999) Associative processes in addiction and reward. The role of amygdala-ventral striatal subsystems. Ann N Y Acad Sci 877:412–438.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Ferraro TN,
    2. Golden GT,
    3. Smith GG,
    4. Martin JF,
    5. Schwebel CL,
    6. Doyle GA,
    7. Buono RJ,
    8. Berrettini WH
    (2005) Confirmation of a major QTL influencing oral morphine intake in C57 and DBA mice using reciprocal congenic strains. Neuropsychopharmacology 30:742–746.
    OpenUrlPubMed
  17. ↵
    1. Freed WJ,
    2. Chen J,
    3. Bäckman CM,
    4. Schwartz CM,
    5. Vazin T,
    6. Cai J,
    7. Spivak CE,
    8. Lupica CR,
    9. Rao MS,
    10. Zeng X
    (2008) Gene expression profile of neuronal progenitor cells derived from hESCs: activation of chromosome 11p15.5 and comparison to human dopaminergic neurons. PLoS ONE 3:e1422.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Grice DE,
    2. Reenilä I,
    3. Männistö PT,
    4. Brooks AI,
    5. Smith GG,
    6. Golden GT,
    7. Buxbaum JD,
    8. Berrettini WH
    (2007) Transcriptional profiling of C57 and DBA strains of mice in the absence and presence of morphine. BMC Genomics 8:76.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Han MH,
    2. Bolaños CA,
    3. Green TA,
    4. Olson VG,
    5. Neve RL,
    6. Liu RJ,
    7. Aghajanian GK,
    8. Nestler EJ
    (2006) Role of cAMP response element-binding protein in the rat locus ceruleus: regulation of neuronal activity and opiate withdrawal behaviors. J Neurosci 26:4624–4629.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Hardy SG,
    2. Haigler HJ
    (1985) Prefrontal influences upon the midbrain: a possible route for pain modulation. Brain Res 339:285–293.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hosack DA,
    2. Dennis G Jr.,
    3. Sherman BT,
    4. Lane HC,
    5. Lempicki RA
    (2003) Identifying biological themes within lists of genes with EASE. Genome Biol 4:R70.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Joe B,
    2. Letwin NE,
    3. Garrett MR,
    4. Dhindaw S,
    5. Frank B,
    6. Sultana R,
    7. Verratti K,
    8. Rapp JP,
    9. Lee NH
    (2005) Transcriptional profiling with a blood pressure QTL interval-specific oligonucleotide array. Physiol Genomics 23:318–326.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Kaneto H,
    2. Yamazaki A,
    3. Kihara T
    (1985) Evidence for the dissociation of morphine analgesia, tolerance and dependence. J Pharm Pharmacol 37:507–508.
    OpenUrlPubMed
  24. ↵
    1. Kaplan GB,
    2. Leite-Morris KA
    (1997) Up-regulation of adenosine transporter-binding sites in striatum and hypothalamus of opiate tolerant mice. Brain Res 763:215–220.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Kest B,
    2. Hopkins E,
    3. Palmese CA,
    4. Adler M,
    5. Mogil JS
    (2002) Genetic variation in morphine analgesic tolerance: a survey of 11 inbred mouse strains. Pharmacol Biochem Behav 73:821–828.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Kest B,
    2. Palmese CA,
    3. Juni A,
    4. Chesler EJ,
    5. Mogil JS
    (2004) Mapping of a quantitative trait locus for morphine withdrawal severity. Mamm Genome 15:610–617.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Korostynski M,
    2. Kaminska-Chowaniec D,
    3. Piechota M,
    4. Przewlocki R
    (2006) Gene expression profiling in the striatum of inbred mouse strains with distinct opioid-related phenotypes. BMC Genomics 7:146.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Korostynski M,
    2. Piechota M,
    3. Kaminska D,
    4. Solecki W,
    5. Przewlocki R
    (2007) Morphine effects on striatal transcriptome in mice. Genome Biol 8:R128.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Law PY,
    2. Wong YH,
    3. Loh HH
    (2000) Molecular mechanisms and regulation of opioid receptor signaling. Annu Rev Pharmacol Toxicol 40:389–430.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Law PY,
    2. Loh HH,
    3. Wei LN
    (2004) Insights into the receptor transcription and signaling: implications in opioid tolerance and dependence. Neuropharmacology 47:300–311.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Lee NH,
    2. Saeed AI
    (2007) Microarrays: an overview. Methods Mol Biol, 265–300.
  32. ↵
    1. Lee NH,
    2. Haas BJ,
    3. Letwin NE,
    4. Frank BC,
    5. Luu TV,
    6. Sun Q,
    7. House CD,
    8. Yerga-Woolwine S,
    9. Farms P,
    10. Manickavasagam E,
    11. Joe B
    (2007) Cross-talk of expression quantitative trait loci within 2 interacting blood pressure quantitative trait loci. Hypertension 50:1126–1133.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Letwin NE,
    2. Kafkafi N,
    3. Benjamini Y,
    4. Mayo C,
    5. Frank BC,
    6. Luu T,
    7. Lee NH,
    8. Elmer GI
    (2006) Combined application of behavior genetics and microarray analysis to identify regional expression themes and gene-behavior associations. J Neurosci 26:5277–5287.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Liang DY,
    2. Guo T,
    3. Liao G,
    4. Kingery WS,
    5. Peltz G,
    6. Clark JD
    (2006) Chronic pain and genetic background interact and influence opioid analgesia, tolerance, and physical dependence. Pain 121:232–240.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Liang M,
    2. Lee NH,
    3. Wang H,
    4. Greene AS,
    5. Kwitek AE,
    6. Kaldunski ML,
    7. Luu TV,
    8. Frank BC,
    9. Bugenhagen S,
    10. Jacob HJ,
    11. Cowley AW Jr.
    (2008) Molecular networks in Dahl salt-sensitive hypertension based on transcriptome analysis of a panel of consomic rats. Physiol Genomics 34:54–64.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Liang Y,
    2. Carr LG
    (1997) Transcription of the mouse mu-opioid receptor gene is regulated by two promoters. Brain Res 769:372–374.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Ling JQ,
    2. Li T,
    3. Hu JF,
    4. Vu TH,
    5. Chen HL,
    6. Qiu XW,
    7. Cherry AM,
    8. Hoffman AR
    (2006) CTCF mediates interchromosomal colocalization between Igf2/H19 and Wsb1/Nf1. Science 312:269–272.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Malek RL,
    2. Irby RB,
    3. Guo QM,
    4. Lee K,
    5. Wong S,
    6. He M,
    7. Tsai J,
    8. Frank B,
    9. Liu ET,
    10. Quackenbush J,
    11. Jove R,
    12. Yeatman TJ,
    13. Lee NH,
    14. et al.
    (2002) Identification of Src transformation fingerprint in human colon cancer. Oncogene 21:7256–7265.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Malek RL,
    2. Wang HY,
    3. Kwitek AE,
    4. Greene AS,
    5. Bhagabati N,
    6. Borchardt G,
    7. Cahill L,
    8. Currier T,
    9. Frank B,
    10. Fu X,
    11. Hasinoff M,
    12. Howe E,
    13. Letwin N,
    14. Luu TV,
    15. Saeed A,
    16. Sajadi H,
    17. Salzberg SL,
    18. Sultana R,
    19. Thiagarajan M,
    20. Tsai J
    (2006) Physiogenomic resources for rat models of heart, lung and blood disorders. Nat Genet 38:234–239.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Manning BH,
    2. Morgan MJ,
    3. Franklin KB
    (1994) Morphine analgesia in the formalin test: evidence for forebrain and midbrain sites of action. Neuroscience 63:289–294.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Mas M,
    2. Sabater E,
    3. Olaso MJ,
    4. Horga JF,
    5. Faura CC
    (2000) Genetic variability in morphine sensitivity and tolerance between different strains of rats. Brain Res 866:109–115.
    OpenUrlCrossRefPubMed
  42. ↵
    1. McClung CA,
    2. Nestler EJ
    (2008) Neuroplasticity mediated by altered gene expression. Neuropsychopharmacology 33:3–17.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Nandigama P,
    2. Borszcz GS
    (2003) Affective analgesia following the administration of morphine into the amygdala of rats. Brain Res 959:343–354.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Narita M,
    2. Aoki T,
    3. Ozaki S,
    4. Yajima Y,
    5. Suzuki T
    (2001) Involvement of protein kinase Cgamma isoform in morphine-induced reinforcing effects. Neuroscience 103:309–314.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Ocaña M,
    2. Entrena JM,
    3. Baeyens JM,
    4. Del Pozo E
    (2007) The antinociceptive effect of morphine is reversed by okadaic acid in morphine-naive but not in morphine-tolerant mice. Pharmacol Biochem Behav 86:21–26.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Oliverio A,
    2. Castellano C
    (1974) Genotype-dependent sensitivity and tolerance to morphine and heroin: dissociation between opiate-induced running and analgesia in the mouse. Psychopharmacologia 39:13–22.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Pan YX
    (2002) Identification and characterization of a novel promoter of the mouse mu opioid receptor gene (Oprm) that generates eight splice variants. Gene 295:97–108.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Pavlidis P,
    2. Noble WS
    (2001) Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol 2:RESEARCH0042.
    OpenUrlPubMed
  49. ↵
    1. Pavlovic ZW,
    2. Cooper ML,
    3. Bodnar RJ
    (1996) Opioid antagonists in the periaqueductal gray inhibit morphine and beta-endorphin analgesia elicited from the amygdala of rats. Brain Res 741:13–26.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Paxinos G,
    2. Franklin KBG
    (2001) The mouse brain in stereotaxic coordinates (Academic, San Diego), Ed 2.
  51. ↵
    1. Petruzzi R,
    2. Ferraro TN,
    3. Kürschner VC,
    4. Golden GT,
    5. Berrettini WH
    (1997) The effects of repeated morphine exposure on mu opioid receptor number and affinity in C57BL/6J and DBA/2J mice. Life Sci 61:2057–2064.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Rapp JP
    (2000) Genetic analysis of inherited hypertension in the rat. Physiol Rev 80:135–172.
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Reiner A,
    2. Yekutieli D,
    3. Benjamini Y
    (2003) Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics 19:368–375.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Reiner-Benaim A,
    2. Yekutieli D,
    3. Letwin NE,
    4. Elmer GI,
    5. Lee NH,
    6. Kafkafi N,
    7. Benjamini Y
    (2007) Associating quantitative behavioral traits with gene expression in the brain: searching for diamonds in the hay. Bioinformatics 23:2239–2246.
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Schalkwyk LC,
    2. Fernandes C,
    3. Nash MW,
    4. Kurrikoff K,
    5. Vasar E,
    6. Kõks S
    (2007) Interpretation of knockout experiments: the congenic footprint. Genes Brain Behav 6:299–303.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Scott RE,
    2. Gao S,
    3. Kim CK,
    4. Cox R
    (2005) De-differentiation-derived mesenchymal stem cells demonstrate selective repression in H19 bioregulatory RNA gene expression. Differentiation 73:294–302.
    OpenUrlPubMed
  57. ↵
    1. Shippenberg TS,
    2. Elmer GI
    (1998) The neurobiology of opiate reinforcement. Crit Rev Neurobiol 12:267–303.
    OpenUrlCrossRefPubMed
  58. ↵
    1. Smith FL,
    2. Javed RR,
    3. Smith PA,
    4. Dewey WL,
    5. Gabra BH
    (2006) PKC and PKA inhibitors reinstate morphine-induced behaviors in morphine tolerant mice. Pharmacol Res 54:474–480.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Tanaka TS,
    2. Jaradat SA,
    3. Lim MK,
    4. Kargul GJ,
    5. Wang X,
    6. Grahovac MJ,
    7. Pantano S,
    8. Sano Y,
    9. Piao Y,
    10. Nagaraja R,
    11. Doi H,
    12. Wood WH 3rd.,
    13. Becker KG,
    14. Ko MS
    (2000) Genome-wide expression profiling of mid-gestation placenta and embryo using a 15,000 mouse developmental cDNA microarray. Proc Natl Acad Sci U S A 97:9127–9132.
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Teramoto H,
    2. Castellone MD,
    3. Malek RL,
    4. Letwin N,
    5. Frank B,
    6. Gutkind JS,
    7. Lee NH
    (2005) Autocrine activation of an osteopontin-CD44-Rac pathway enhances invasion and transformation by H-RasV12. Oncogene 24:489–501.
    OpenUrlCrossRefPubMed
  61. ↵
    1. Tortorici V,
    2. Morgan MM,
    3. Vanegas H
    (2001) Tolerance to repeated microinjection of morphine into the periaqueductal gray is associated with changes in the behavior of off- and on-cells in the rostral ventromedial medulla of rats. Pain 89:237–244.
    OpenUrlCrossRefPubMed
  62. ↵
    1. Trang T,
    2. Sutak M,
    3. Jhamandas K
    (2007) Involvement of cannabinoid (CB1)-receptors in the development and maintenance of opioid tolerance. Neuroscience 146:1275–1288.
    OpenUrlCrossRefPubMed
  63. ↵
    1. Williams JT,
    2. Christie MJ,
    3. Manzoni O
    (2001) Cellular and synaptic adaptations mediating opioid dependence. Physiol Rev 81:299–343.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Wise RA
    (1996) Neurobiology of addiction. Curr Opin Neurobiol 6:243–251.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Yang IV,
    2. Chen E,
    3. Hasseman JP,
    4. Liang W,
    5. Frank BC,
    6. Wang S,
    7. Sharov V,
    8. Saeed AI,
    9. White J,
    10. Li J,
    11. Lee NH,
    12. Yeatman TJ,
    13. Quackenbush J
    (2002) Within the fold: assessing differential expression measures and reproducibility in microarray assays. Genome Biol 3:research0062.
  66. ↵
    1. Zachariou V,
    2. Bolanos CA,
    3. Selley DE,
    4. Theobald D,
    5. Cassidy MP,
    6. Kelz MB,
    7. Shaw-Lutchman T,
    8. Berton O,
    9. Sim-Selley LJ,
    10. Dileone RJ,
    11. Kumar A,
    12. Nestler EJ
    (2006) An essential role for DeltaFosB in the nucleus accumbens in morphine action. Nat Neurosci 9:205–211.
    OpenUrlCrossRefPubMed
  67. ↵
    1. Zubieta JK,
    2. Smith YR,
    3. Bueller JA,
    4. Xu Y,
    5. Kilbourn MR,
    6. Jewett DM,
    7. Meyer CR,
    8. Koeppe RA,
    9. Stohler CS
    (2001) Regional mu opioid receptor regulation of sensory and affective dimensions of pain. Science 293:311–315.
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

The Journal of Neuroscience: 29 (16)
Journal of Neuroscience
Vol. 29, Issue 16
22 Apr 2009
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Identification of Candidate Genes and Gene Networks Specifically Associated with Analgesic Tolerance to Morphine
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Identification of Candidate Genes and Gene Networks Specifically Associated with Analgesic Tolerance to Morphine
Jenica D. Tapocik, Noah Letwin, Cheryl L. Mayo, Bryan Frank, Troung Luu, Ovokeraye Achinike, Carrie House, Russell Williams, Greg I. Elmer, Norman H. Lee
Journal of Neuroscience 22 April 2009, 29 (16) 5295-5307; DOI: 10.1523/JNEUROSCI.4020-08.2009

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Identification of Candidate Genes and Gene Networks Specifically Associated with Analgesic Tolerance to Morphine
Jenica D. Tapocik, Noah Letwin, Cheryl L. Mayo, Bryan Frank, Troung Luu, Ovokeraye Achinike, Carrie House, Russell Williams, Greg I. Elmer, Norman H. Lee
Journal of Neuroscience 22 April 2009, 29 (16) 5295-5307; DOI: 10.1523/JNEUROSCI.4020-08.2009
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Articles

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Articles

Behavioral/Systems/Cognitive

  • Choice Behavior Guided by Learned, But Not Innate, Taste Aversion Recruits the Orbitofrontal Cortex
  • Maturation of Spontaneous Firing Properties after Hearing Onset in Rat Auditory Nerve Fibers: Spontaneous Rates, Refractoriness, and Interfiber Correlations
  • Insulin Treatment Prevents Neuroinflammation and Neuronal Injury with Restored Neurobehavioral Function in Models of HIV/AIDS Neurodegeneration
Show more Behavioral/Systems/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2022 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.