|Home | About | Journals | Submit | Contact Us | Français|
Abuse of heroin and prescription opiate medications has grown to disturbing levels. Opioids mediate their effects through mu opioid receptors (MOR), but minimal information exists regarding MOR-related striatal signaling relevant to the human condition. The striatum is a structure central to reward and habitual behavior and neurobiological changes in this region are thought to underlie the pathophysiology of addiction disorders.
We examined molecular mechanisms related to MOR in postmortem human brain striatal specimens from a homogenous European Caucasian population of heroin abusers and control subjects and in an animal model of heroin self-administration. Expression of ets-like kinase 1 (ELK1) was examined in relation to polymorphism of the MOR gene OPRM1 and drug history.
A characteristic feature of heroin abusers was decreased expression of MOR and extracellular regulated kinase (ERK) signaling networks, concomitant with dysregulation of the downstream transcription factor ELK1. Striatal ELK1 in heroin abusers associated with the polymorphism rs2075572 in OPRM1 in a genotype dose-dependent manner and correlated with documented history of heroin use, an effect reproduced in an animal model that emphasizes a direct relationship between repeated heroin exposure and ELK1 dysregulation. A central role of ELK1 was evidenced by an unbiased whole transcriptome microarray that revealed ~20% of downregulated genes in human heroin abusers are ELK1 targets. Using chromatin immuneprecipitation, we confirmed decreased ELK1 promoter occupancy of the target gene Use1.
ELK1 is a potential key transcriptional regulatory factor in striatal disturbances associated with heroin abuse and relevant to genetic mutation of OPRM1.
Opioids have become the second most abused class of drugs in the United States and are prevalent among teens and young adults.(1) Given that opioid abuse is associated with high rates of drug-related mortality, morbidity and criminality, the ramifications of this trend are cause for public health concern.(2) The rewarding effects of heroin are mediated by the activation of the G-protein coupled μ-opioid receptor (MOR).(3, 4) MOR stimulation leads to modulation of multiple effectors, including the inhibition of adenylyl cyclase,(5) regulation of voltage-gated calcium (Ca2+) channels,(6, 7) and control of the canonical mitogen-activated protein kinases (MAPK) extracellular regulated kinase (ERK)-1 and -2.(8) Active ERK1/2 phosphorylates the transcription factors ELK1 and cyclic AMP responsive element binding protein (CREB), which translocate to the nucleus to regulate expression of target genes,(9) the physiological consequences of which result in synaptic plasticity and aberrant behaviors that maintain the cycle of drug abuse.(10) Experimental models have led to the postulation that long-term repeated exposure to opioid drugs results in receptor desensitization and/or down-regulation events linked to opioid tolerance and dependence.(11, 12) However, the signaling and transcriptional events that underlie adaptation within the human brain as a result of repeated MOR activation remain poorly understood.
An important consideration in human studies is genetics which confer an increased vulnerability to heroin addiction (13, 14) and gene variants encoding opioid receptors have been shown to be associated with this disorder.(15–18) Numerous single nucleotide polymorphisms (SNPs) have been identified in the human MOR gene (OPRM1) and while the functional relevance of some variants to MOR signaling(19–21) and expression(22) has been described, the downstream transcriptional events have not been well characterized.
In the present study, we investigated signaling molecules directly associated with MOR function and downstream transcription factors in the striatum, a central structure involved in mediating goal-directed behavior, habit formation and reward.(23–25) We report dysregulation of striatal MOR and ERK signaling pathways in human heroin abusers and identify ELK1 as a key transcriptional regulator for differentially expressed genes in the striatum of heroin abusers and in association with OPRM1 variants. These findings highlight ELK1 as a potentially new molecular target in the cellular dysregulation of opioid addiction.
Postmortem brains samples from heroin abusers and control subjects were collected at the Department of Forensic Medicine at Semmelweis University, Hungary, and the National Institute of Forensic Medicine at Karolinska Institutet, Stockholm, Sweden under the ethics guidelines approved by each institution. Heroin subjects died from heroin intoxication, were predominantly heroin users not receiving methadone or buprenorphine treatment, and were negative for human immunodeficiency virus infection. No information was known about their physical dependence. The control group had negative blood levels of opiates and other drugs of abuse. Demographic information can be found in Supplemental Table 1. Striatal tissue was dissected as previously described.(18) Due to the low abundance of nucleus accumbens (NAc) tissue, putamen tissue was used for the majority of protein experiments. Putamen punches (~200 mg) were pulverized on dry ice and stored at −80°C. For laser captured samples, striatal regions were dissected from 20μm sections from human subjects using an Arcturus XT Laser Capture Microdissection Instrument (Molecular Devices, Sunnyvale, CA).
A heroin self-administration paradigm was performed in adult male Long Evans rats (Charles River, Wilmington, MA) for two weeks as previously described.(26–28) This dosing paradigm does not induce physical dependence.(29) Animals were sacrificed either 1 hour or 24 hours after the last self-administration session by CO2 followed by decapitation, brains frozen and striata dissected using a 2mm micro-punch.
Genomic DNA was purified from human cerebellar tissue using DNeasy columns (Qiagen, Valencia, CA). Genotyping for the OPRM1 rs1799971 (Taqman Assay ID: C_8950074_1_) and rs2075572 (Taqman Assay ID: C_1691815_1_, Applied Biosystems, Foster City, CA) SNPs was performed as previously described.(18, 30) The PLINK 1.07 genetic association analysis program(31) was used to verify SNP data quality, test for departure from Hardy-Weinberg equilibrium, and test individual SNPs for statistical association with significance set at p<0.05.
Western blots were performed on a subset of heroin (N=32) and control (N=15) human subjects, and on dorsal striata from the rodent self-administration study, as previously described(32) and probed with primary antibodies directed against the following proteins: MOR (GeneTex, Inc., Irvine, CA); ERK1/2, ELK1 and phospho-ELK1(Ser383)(Cell Signaling Technology, Danvers, MA for human samples and Santa Cruz Biotechnology, Santa Cruz, CA for rat); βarrestin2 (Santa Cruz Biotechnology); and MEK1 (Invitrogen, Carlsbad, CA). Proteins were analyzed using ImageJ software(33, 34) and normalized to total protein levels(35) obtained with the Memcode Reversible Protein Stain Kit (Thermo Scientific). Logarithmic transformations were performed to render normal data distribution. A general linear stepwise regression was used to calculate statistical significance and identify covariates (brain pH, gender, age) using JMP software (SAS Institute, Cary, NC). Two-tailed t-tests were used when no covariates were found. Spearman correlations were calculated to assess the relationship between protein levels and heroin toxicology.
mRNA microarray experiments were performed on NAc tissue from a subset of heroin (N=22; 19 males and 3 females) and control (N=27; 22 males and 5 females) subjects at the Purdue Pharma, L.P. (Cranbury, NJ). Total RNA isolation, assessment of RNA integrity, and microarray hybridization were performed as previously described.(36) Expression data were collected using HG-U133A Affymetrix Chips (Affymetrix, Santa Clara, CA). Raw data were normalized using the Robust Multichip Average (RMA)(37) from the Affymetrix Expression Console. To control for RNA quality, actin 3′/5′ ratios of less than 3 and % present calls of 50% or greater were used. To identify Differentially Expressed Genes (DEGs), two-tailed t-tests were performed with a multiple-test correction (p<0.05) (38). Hierarchical clustering was performed on log2 expression levels of DEGs with Euclidian distance and average-linkage for constructing dendograms. Common transcription factors for DEGs were identified with Lists2Networks, a web-based software system for performing gene list enrichment analysis.(39) The transcription factor gene set library in List2Networks was created from putative targets of transcription factors.(40) All other computational analyses were performed in MATLAB. (Mathworks, Natick, MA)
Microarray findings were verified with Nanostring Technologies (Nanostring, Seattle, WA) using 100ng NAc RNA from human heroin (N=20) and control (N=16) subjects. Design and synthesis of probes was performed according to company protocol and hybridization reactions were carried out as previously described.(41) Expression values were normalized to ACTB and GAPDH, which were stable in the human population. Probe sequences are listed in supplementary table 2. Statistical analyses were performed as above using JMP software.
Chromatin immuneprecipitation (ChIP) was performed as previously described (42) using the NAc from rats that underwent the self-administration paradigm. 30ug of purified chromatin per animal was immuneprecipitated with 1.25ug of ELK1 antibody (Epitomics, Burlingame, CA) or nonspecific rabbit IgG (Santa Cruz Biotechnology, Santa Cruz, CA). Chromatin was analyzed by QPCR at genomic regions of BAG1 and USE1. Primer sequences are listed in supplementary table 3. Statistical analyses as described above were performed using JMP software.
Protein expression levels of MOR (89.55±3.42%, p<0.005; F2,43=6.73) and β-arrestin2 (87.11±4.26%, p<0.05; F3,41=8.1), which regulates MOR desensitization and internalization,(43) were significantly reduced in heroin abusers compared to controls (Fig 1A–B), indicating that the tone of MOR-signaling cascades may be altered. One of the key intracellular pathways activated upon MOR stimulation is the MAPK pathway, depicted in Figure 1C. We focused on the core components of the pathway including MAPK ERK kinases (MEK)1/2 and ERK1/2 because they regulate phosphorylation status and activity of transcription factors. Levels of both MEK1 (91.28±2.55%, p<0.05; F1,41=4.33) and its substrate ERK1 (84.78 ± 3.83%, p<0.05; F1,46= 4.40) were significantly reduced. We observed an opposing increase in the protein levels of ERK2 (111.78±2.32%. p<0.005; F1,44=9.33). The levels of phosphorylated ERK1 (105.67 ± 13.01%) and phosphorylated ERK2 (120.1 ± 11.4%) were not different from controls due to the strong confound of brain pH which significantly impacted pERK1 (p<0.0005) and pERK2 (p< 0.00001) levels (data not shown), suggesting that phosphorylation of ERK1/2 is sensitive to the perimortem state (data not shown).
To assess further downstream events of MAPK signaling, we examined expression of transcription factors ELK1 and CREB, which are nuclear targets for ERK1/2.(44, 45) CREB is activated in response to ERK1/2 through signaling by kinases such as p90RSK.(46–48) Both p90RSK (82.01±10.2% control) and CREB (79.07±4.69% control) were reduced in the putamen, but this was confounded by sensitivity to pH (p<0.0001). However, heroin abusers had significantly higher expression of ELK1 (133.26±8.21%, p<0.01; F1,46=7.38). In order to gain entry into the nucleus, ELK1 must be phosphorylated by ERK1/2(49) and the levels of activated pELK1 were significantly reduced (86.92±3.35%, p<0.05; F1,41=5.21) in heroin subjects. Furthermore, the ratio of pELK1 to total ELK1 was significantly decreased in heroin abusers (66.37±5.46, p<.001; F1,46=12.01; Fig 1D), indicating that repeated heroin use leads to an imbalance in total versus phosphorylated ELK1 protein.
The NAc plays a key role in heroin abuse vulnerability,(50, 51) thus we next examined protein levels of ELK1 in the NAc to determine whether alterations in the putamen were also evident in the mesolimbic component of the striatum. Laser capture microdissected samples were taken from the NAc, putamen and caudate nucleus from glass-mounted striatal sections from the same subjects studied above (Fig. 2A). Similar to the observations in putamen tissue punches, ELK1 expression was significantly elevated in all striatal regions of the heroin abusers, but was most pronounced in the NAc (p<0.05; F2,41=4.923, Fig. 2B).
Given the contribution of genetics to opiate abuse, we determined if individual polymorphisms of OPRM1 are related to ELK1 expression. The rs1799971 SNP (A118G variant) in exon 1 is the most common mutation in the coding region as well as the most studied OPRM1 polymorphism.(52) Unfortunately, only two control subjects carried the minor G allele which rendered it challenging to interpret the results in regard to non-drug protein expression. Nevertheless, there was an overall group-genotype effect (F2,42=4.789, p<0.05) with ELK1 levels significantly higher in heroin subjects than controls for the A/A genotype (p<.001) and heroin subjects carrying the G allele had the highest ELK1 levels (p<0.01; Fig. 3A).
Considering the low frequency of rs1799971, we also evaluated the rs2075572 SNP, which is located in intron 2(53) and has been associated with initiation of smoking.(54) There was not an overall significant association of rs2075572 with heroin abuse, but a greater number of G/G subjects tended to be heroin abusers, while the control group tended to have higher numbers of C/C individuals (Supplementary Table 4). We observed a pronounced genotype-dose effect association of rs2075572 with striatal ELK1 protein levels (F4,46=4.338, p<0.01;Fig. 3B). Heroin carriers of the minor C allele had significantly higher ELK1 levels than heroin subjects homozygous for the G allele (Fig 3B). A similar trend was observed for changes in ELK1 protein expression between control subjects, indicating that individual variation in the OPRM1 gene may result in heterogeneous levels of ELK1. Importantly, the mean value of ELK1 expression for each individual heroin genotype (C/C- 170.8%; C/G- 129.2%; G/G- 106.6%) was higher than the mean value for all control subjects combined (100%).
Assessment of toxicology in heroin abusers revealed a strong positive correlation between ELK1 protein expression in the putamen and urine concentrations of 6-MAM, the rapid metabolite of heroin (r2=0.70, p<0.05), but this was strongly influenced by one subject (Figure 4A). Interestingly, documentation of drug use history available for some subjects showed that levels of ELK1 mRNA tended to be negatively correlated with age of onset of heroin use for subjects carrying the G allele of rs2075572 (r2=0.37, p<0.05), suggesting elevated ELK1 with long-term heroin use, but this was predominantly driven by those with the G/G genotype (r2=0.95, p<0.01; Fig 4B). Importantly, this subset of G allele heroin abusers was representative of the larger population of subjects as their ELK1 mRNA expression levels (124.00±4.04% of controls, p<0.01; F1,28=10.328, data not shown) were comparable to that in the larger cohort.
Given the challenge of human postmortem studies that lack a complete drug history about the subjects, we examined the regulation of ELK1 in a rat model of addiction using a heroin self-administration paradigm where animals controlled their individual drug intake.(26–28) Since the post-mortem interval for many of the human heroin subjects is relatively short, we examined ELK1 protein 1 hour after the last self-administration session as well as 24 hours later to obtain insights regarding the stability of the observations. ELK1 expression was tightly correlated with each animal’s history of heroin intake. We observed a positive correlation 1 hour after the final heroin self-administration session (r2=0.79, p<0.05; Fig 4C) that was reversed 24 hours later with a negative correlation to ELK1 levels (r2=0.72, p<0.001; Fig 4D).
To assess the functional state of ELK1 in relation to heroin exposure, we also measured levels of pELK1 in animals that self-administered the drug. While no changes in pELK1 were observed at the 1 hour time-point, we observed a significant increase in pELK1 relative to both total protein and total ELK1 protein at 24 hours after the last heroin session (pELK1/total protein: 201.31 ± 32.13%, p<0.05; F1,10= 6.03; pELK1/total ELK1: 158.14 ± 16.10%, p<0.05; F1,10= 6.50; Figure 5A.). Correlations to heroin intake revealed that at both time points, pELK1 negatively correlated to heroin intake (1 hour: r2=0.86, p<0.05; Figure 5B; 24 hours: r2=0.79, p<0.05; Figure 5C), indicating that high levels of heroin were associated with decreased active ELK1. These data support the notion that repeated heroin use dysregulates ELK1 expression relative to the overall history of drug exposure which may also account for individual differences observed in humans.
We next addressed whether ELK1 could play a role in transcriptional regulation relevant to heroin abuse. We chose an unbiased strategy in which downstream gene targets could be identified to guide future studies. As such, we utilized a NAc microarray dataset from a cohort of our heroin postmortem population and examined transcription factor binding motifs common to the observed transcriptome alterations. Interestingly, evaluation of common transcription factors for the significantly regulated genes using the Lists2Networks program (55) revealed that ELK1 (p<0.05) and ETS1 (p<0.01) had the most putative target genes overlapping with down-regulated genes (36/195) (Fig. 6A and Supplementary Table 5), confirming the relevance of ELK1 as a transcriptional regulator of gene expression in heroin abuse. Hierarchical clustering on a subset (N=32) of subjects in the two groups that were most different from one another revealed differential expression of 1518 genes, of which 115 were identified as downregulated ELK1 target genes (p<.000001) (Supplementary Fig 1 and Supplementary Table 6). These findings were validated by measuring mRNA expression levels of the ELK1 target genes BAG1 (80.46 ± 3.61%; p<0.01) and USE1 (92.09 ± 2.76%; p<0.05; Fig 6B) in the NAc using Nanostring technology. In addition, Nanostring analyses revealed a significant increase in ELK1 mRNA in heroin subjects compared to controls (120.67 ± 3.72%; p<0.01; Fig 6B), which is in agreement with the increased protein expression of ELK1 previously observed.
ELK1 binding sites are located slightly upstream of the transcriptional start site of target genes.(56) To determine if the gene expression changes could be attributed to decreased transcriptional activity of ELK1 directly at the target loci, we used ChIP to assess the occupancy of ELK1 at the promoter regions of downstream targets in the NAc of rats that underwent the self-administration paradigm. We examined published ChIP sequencing data obtained using ELK1, CREB, and SRF antibodies to guide our selection of ELK1 target genes.(56, 57) Such studies suggested Bag1 and Use1 are regulated by ELK1 in a CREB- and SRF-independent manner.(56, 57) We observed no enrichment of ELK1 binding over nonspecific IgG controls at any region of Bag1 studied (data not shown). In animals that were sacrificed 1 hour after the last heroin self-administration session, ELK1 occupancy was significantly decreased at the promoter of Use1 versus saline animals (0.034± 0.007% versus 0.058±0.006; p<.05; Fig 6C), indicating a direct functional regulation of ELK1 on downstream gene targets as a consequence of heroin use. At the 24-hour time-point, ELK1 binding at the Use1 promoter was not significantly different from controls, indicating that heroin did not induce a prolonged decrease in ELK1 occupancy.
The current study provides convergent lines of evidence identifying ELK1 disturbances in association with heroin use history and variants of the OPRM1 gene. The findings suggest that repeated activation of MOR, and subsequently the ERK pathway, efficiently controls the availability of ELK1 to regulate striatal transcriptional machinery in a dose-dependent manner. The ELK1 transcription factor has been generally overlooked in the context of addiction, but the current results implicate it as a potential substrate in the pathophysiology of heroin abuse.
Limited insights are available regarding the role of ELK1 in neuropsychiatric disorders, but its molecular regulation has been described. Phosphorylation of ELK1 by ERK1/2 is critical for its translocation into the nucleus where it can exert effects on chromatin and gene expression.(58) The fact that total ELK1 was increased in the putamen of human heroin abusers concomitant with a decrease in pELK1 levels suggested that less ELK1 is trafficked to the nucleus and transcription of ELK1 target genes may be decreased. Indeed, total ELK1 was positively correlated to heroin at 1 hour after heroin intake in animals that directly self-administered heroin, similar to human abusers, and pELK1 levels were negatively correlated with heroin use, irrespective of the time of death following the last drug intake. These findings emphasize an important contribution of the history of repeated heroin use on ELK1 phosphorylation. Interestingly, the absolute levels of pELK1 were elevated 24 hours after heroin use, indicative of an acute withdrawal effect or a potential compensation to counter reduction of pELK1 as a consequence of repeated heroin use. Acute opiate withdrawal 24 hours following heroin self-administration is generally associated with enhanced stress response that involves the activation of the hypothalamic-pituitary axis and glucocorticoid receptor system.(59, 60) Interestingly, glucocorticoid receptor stimulation increases pELK1 (61), suggesting that increased pELK1 24 hours after heroin may be a result of withdrawal-induced stress. Such time course changes could also have distinct functional relevance depending on the cellular localization of pELK1 at these periods since ELK1 localization in dendritic, cytoplasmic and nuclear compartments has significant implications regarding its effects on dendritic sprouting and transcriptional regulation. (62–64)
Although ELK1 has not been widely investigated in addiction, other animal studies have shown that acute administration of cocaine (65, 66) or delta-9-tetrahydrocannabinol (psychoactive component of cannabis)(67) increases pELK1 immediately after drug exposure. Whether these noted differences in pELK1 relate to acute versus chronic drug exposure, type of drugs, cellular localization or other factors in experimental designs of these studies needs to be further explored. However, the fact that ELK1 is altered by several drugs of abuse is intriguing considering the role of ELK1 phosphorylation in neurite outgrowth (62–64) and that dysregulation of structural plasticity is a feature of drugs of abuse.
ELK1 is a member of the ternary complex factor subfamily of ETS transcription factors and intriguingly, our unbiased computational analysis of the NAc microarray data also identified ETS1, the prototype of the ETS1 family, as enriched in the promoter region of down-regulated genes of heroin abusers. In addition to cellular remodeling linked to synaptic plasticity, ETS1 family members are known to play important roles in various cellular functions such as proliferation,(68) differentiation,(69) migration,(70) and apoptosis.(71) Other transcriptional regulators implicated in addiction such as CREB also target some genes common to ELK1. By focusing on gene targets regulated by ELK1 in a CREB- and SRF-independent manner, we confirmed reduction of the genes BAG1 and USE1 identified in the microarray. BAG1 has antiapoptotic effects that confer prosurvival signals to enhance neurite outgrowth activity and is activated downstream of RAF1 and MEK activity.(72–74) Unfortunately, Bag1 shares an ELK1 promoter binding site with another gene in close proximity, Chmp5, which makes the assessment of transcription factor binding difficult to distinguish with ChIP.(56) We were, nevertheless, able to detect ELK1 at the promoter of Use1, a gene known to be antiapoptotic (75), and to act as an E2 ubiquitin-conjugating enzyme to aide in the proteasomal degradation of proteins.(76) ELK1 occupancy at Use1 was decreased 1 hour after heroin self-administration, validating direct gene transcription regulation in the NAc after repeated heroin exposure. ELK1 occupancy at the Use1 promoter was normalized 24 hours later, in line with the dynamic nature of transcription factors. Clearly, the time course of ELK1 regulation is an important question. Moreover, while the current study focused on SRF-independent downstream targets, SRF-dependent ELK1 transcriptional activity is also relevant to addiction and warrants future studies in heroin abuse.
The importance of individual factors to striatal ELK1 expression was also evident in our study, both in relation to drug use history as supported by our animal model and by mutation of the OPRM1 gene. The rs2075572 OPRM1 variant has not been studied in opioid abuse and the functional consequences of this variant on MOR activity are not known. However, a similar genotype-dose effect pattern was evident for the well-known A118G SNP in heroin abusers, suggesting that these mutations relate to general functioning of MOR and its downstream signaling. Although mutations of OPRM1 contributed to variability in the individual ELK1 expression levels, it is important to emphasize that heroin use still enhanced ELK1 at each respective genotype and thus, history of drug use plays a significant role in ELK1 regulation.
While we concentrated on the novel ELK1 findings, there is also clear dysregulation of other components of striatal MOR signaling, including MOR and β-arrestin2. Such impairments may relate to altered drug sensitivity or tolerance to heroin, as β-arrestin2 modulates MOR internalization and desensitization in rodent models.(77) These events contribute significantly to drug tolerance as evidenced by the fact that β-arrestin2 knockout mice do not develop tolerance.(43) A lack of MOR desensitization or proper recycling due to deficits in β-arrestin2 may lead to heightened reward response with repeated heroin use.(78) Human opiate abusers also have reduced expression in the frontal cortex of MOR, β-arrestin2, and GRK2,(79) a kinase that also contributes to MOR desensitization and internalization, which together with the current findings would suggest similar MOR-related pathology in brain areas critically involved in addictive behaviors.
Another key finding indicating compromised integrity of MOR signaling in heroin abuse was the dysregulation of the MAPK/ERK pathway, including differential alterations of ERK1 and ERK2. Most studies have normally assumed that measures of ERK1 are indicative of ERK2 alterations since MEK phosphorylates both these kinases. However, ERK1 inhibits ERK2.(80) Thus increased ERK2 levels observed in the heroin abusers would be in line with compensation for the decreased ERK1 expression. Altogether such aberrant signaling activity of MEK1 and ERK1/2 would be expected to contribute to subsequent downstream changes in ELK1 and gene expression that mediate cellular adaptations to repeated heroin exposure. In summary, the current results demonstrate that the striatum of human heroin abusers is characterized by dysregulation of MOR signaling and the MAPK pathway linked to disturbances of the downstream transcription factor, ELK1. The importance of ELK1 impairment is highlighted by the fact that it targeted the promoter of a large proportion of the repressed genes in the NAc of heroin abusers. The significant contribution of heroin intake history and OPRM1 polymorphisms to individual differences in ELK1 expression also emphasizes the strong link between MOR and ELK1. Overall, these findings suggest that ELK1 may be a key transcriptional substrate in pathophysiological disturbances relevant to drug addiction.
Figure S1: Heatmap of gene expression changes in the NAc of human heroin abusers. The hierarchical clustering of differentially expressed genes detected between the control (N=19) and heroin (N=13) groups in the NAc using mRNA microarrays. Dark colors represent normalized low expression and light colors high expression. The x-axis lists subjects represented by a three digit ID number; C-control, H-heroin user; m-male, f-female; age.
Table SI: Demographic data
Table S2: Probes used with Nanostring technologies.
Table S3- Primer pair sequences used in qPCR.
Table S4: Frequency of the rs2075572 OPRM1 polymorphism in human subjects.
Table S5: ELK1 target genes identified in 49 subjects.
Table S6: ELK1 target genes identified in 32 subjects.
We thank James Sperry and Nayana Patel for technical assistance as well as Dr. Gang Jin previously at Purdue Pharma L.P. for microarray analysis. This work was funded by grants from the National Institute on Drug Abuse DA015446 (YLH), DA013997 (YXP), and T32DA007135 (MMJ), as well as grant ETT/098/2009 (EK).
Disclosure of conflicts of interest
The authors declare that we have no competing interests.