- Open Access
Identification and characterization of seven new exon 11-associated splice variants of the rat mu opioid receptor gene, OPRM1
Molecular Pain volume 7, Article number: 9 (2011)
The mouse mu opioid receptor (OPRM1) gene undergoes extensive alternative splicing at both the 3'- and 5'-ends of the gene. Previously, several C-terminal variants generated through 3' splicing have been identified in the rat OPRM1 gene. In both mice and humans 5' splicing generates a number of exon 11-containing variants. Studies in an exon 11 knockout mouse suggest the functional importance of these exon 11-associated variants in mediating the analgesic actions of a subset of mu opioids, including morphine-6β-glucuronide (M6G) and heroin, but not others such as morphine and methadone. We now have examined 5' splicing in the rat.
The current studies identified in the rat a homologous exon 11 and seven exon 11-associated variants, suggesting conservation of exon 11 and its associated variants among mouse, rat and human. RT-PCR revealed marked differences in the expression of these variants across several brain regions, implying region-specific mRNA processing of the exon 11-associated variants. Of the seven rat exon 11-associated variants, four encoded the identical protein as found in rMOR-1, two predicted 6 TM variants, and one, rMOR-1H2, generated a novel N-terminal variant in which a stretch of an additional 50 amino acids was present at the N-terminus of the previously established rMOR-1 sequence. When expressed in CHO cells, the presence of the additional 50 amino acids in rMOR-1H2 significantly altered agonist-induced G protein activation with little effect on opioid binding.
The identification of the rat exon 11 and its associated variants further demonstrated conservation of 5' splicing in OPRM1 genes among rodents and humans. The functional relevance of these exon 11 associated variants was suggested by the region-specific expression of their mRNAs and the influence of the N-terminal sequence on agonist-induced G protein coupling in the novel N-terminal variant, rMOR-1H2. The importance of the exon 11-associated variants in mice in M6G and heroin analgesia revealed in the exon 11 knockout mouse implies that these analogous rat variants may also play similar roles in rat. The complexity created by alternative splicing of the rat OPRM1 gene may provide important insights of understanding the diverse responses to the various mu opioids seen in rats.
Three families of opioid receptors were proposed from pharmacological studies [1, 2]. Of the three opioid receptor families, the mu opioid receptors are particularly important since they mediate the actions of most of the clinically relevant opioids, as well as those most widely abused such as heroin. Clinicians have observed a wide range of responses among patients, a variability confirmed among different strains of mice. These findings, along with receptor binding studies and the actions of selective antagonists, led us to propose the existence of multiple mu opioid receptor subtypes  long before the molecular characteristics of mu receptors were known.
The molecular cloning of the mu opioid receptor (MOR-1) [4–6] opened new opportunities to investigate the molecular underpinnings for the concept of multiple mu opioid receptors. A single mu opioid receptor gene (OPRM1) has been identified in mammals, raising questions on how to reconcile a single gene with the multiple pharmacologically defined mu opioid receptors. One possibility is alternative pre-mRNA splicing, which can provide enormous RNA and protein diversity. A number of G protein-coupled receptors undergo alternative splicing, such as dopamine D2 [7, 8], somatostatin 2 , prostaglandin EP3 , serotonin receptor subtypes , tachykinin NK(2) , metabotropic glutamate receptor [13, 14], and metabotropic muscarinic acetylcholine receptors [15, 16]. Antisense mapping studies provided early evidence suggesting alternative splicing of the mouse and rat OPRM genes [17, 18], which was further supported by the studies in an exon 1 knockout (KO) mouse model generated by Pintar and colleagues . In this mouse, loss of exon 1 eliminated all the full length variants, which contain exon 1. However, a series of exon 11-associated variants lacking exon 1 were still expressed. Pharmacologically, disrupting exon 1 in this mouse completely abolished morphine analgesia, but not that of either M6G or heroin, consistent with the possibility that alternatively spliced transcripts lacking exon 1 might be responsible for the residual M6G and heroin actions.
In recent years, alternative splicing of the OPRM1 genes has been extensively explored by our group and others [20–34]. In the mouse OPRM1 gene, over 28 alternatively spliced variants have been isolated. Of these splice variants, most are C-terminal variants that were generated through alternative splicing between exon 3 and 10 different downstream exons . These C-terminal variants from mice, rats and humans bound mu opioids with similar high affinities, but displayed marked differences in agonist-induced G protein coupling in both their potency, defined by the EC50 values, and efficacy, indicated by the maximal stimulation[25, 27, 36, 37]. Although it can be speculated that different C-terminal tails may alter interactions of receptor with different G proteins or other related proteins like regulator of G protein signaling (RGS) proteins based upon their intracellular location, the underlying mechanisms for these differences remain unclear. Morphine-induced internalization also varied among the C-terminal variants. For example, morphine given intracerebroventricularly in vivo internalized mMOR-1C in the mouse lateral septum, while mMOR-1 is not internalized by morphine . Mu agonist-induced interaction of phosphorylated mu opioid receptor with β-arrestins have been indicated to involve the receptor internalization and desensitization [39–41]. Various additional phosphorylation sites for β-adrenergic receptor kinase, protein kinase C, caseine kinase, tyrosine kinase and cAMP- and cGMP-dependent protein kinases have been predicted among different C-termini . Although highly speculative, these phosphorylation sites at different C-termini may differentially modulate the recruitment of β-arrestins and therefore contribute to the disparities in mu agonist-induced receptor internalization.
In the mouse, OPRM1 generates a set of splice variants associated with exon 11, located approximately 30 kb upstream of exon 1, under the control of a distinct exon 11 promoter [24, 28]. Of nine exon 11-associated variants, three variants encoded the original mMOR-1 protein, five variants lacked exon 1 and predicted a 6 TM receptor protein, and one variant predicted a protein with single TM. The functional relevance of exon 11-associated variants was established by studies in an exon 11 KO mouse model . Unlike the exon 1 KO mouse developed by Pintar , the exon 11 KO mouse retained full sensitivity towards morphine and methadone analgesia while the effects of M6G, fentanyl and heroin were greatly attenuated. This suggested that exon 11 associated variants mediated the actions of a subset of mu opioids, including M6G and heroin.
Previous studies also have reported many splice variants from the rat OPRM1 gene [20, 25] (Figures 1 & 2), as well as the human OPRM1 gene [21, 26, 28, 32, 33]. The current study reports the identification and characterization of the rat exon 11 homolog and seven exon 11-associated variants.
Cloning the rat exon 11-associated splice variants
To determine whether the rat OPRM1 gene contained an exon homologous to the mouse exon 11, we blasted the rat genome database in Ensembl using the mouse exon 11 sequence and found a highly homologous sequence with 86% identity, which included both coding and adjacent intron regions that included the splice site (Figure 3). We designated this sequence as the rat exon 11. The rat exon 11 was located at about 21 kb upstream of exon 1 in the rat OPRM1 locus of chromosome 1, a distance similar to the 30 kb seen in the mouse OPRM1 gene  and the 28 kb in the human OPRM1 gene . However, the sequence of the rat exon 11 predicted only seven amino acids before encountering a stop codon. To isolate potential rat exon 11-associated splice variants homologous to those identified in the mouse OPRM1 gene, we performed RT-PCR using sense primers designed from the rat exon 11 sequence together with two antisense primers from the 3'UTR of exon 4. We identified seven exon 11-associated splice variants, rMOR-1G1, rMOR-1G2, rMOR-1H1, rMOR-1H2, rMOR-1i1, rMOR-1i2 and rMOR-1i3, from rat brain (Figures 1B, 2 & 4).
The rat exon 11 contained an alternative splice site that divided the exon into two parts, a pattern similar to that seen in the human exon 11 , which we assigned as exon 11a and exon 11b. Alternative usage of the two splice sites in the exons 11a and 11b led to the series of splice variants. rMOR-1G1 contained exons 11a/11b/2/3/4. If the AUG in exon 11a was used, rMOR-1G1 only predicted a peptide with seven amino acids since the stop codon predicted from exon 11b terminated its translation. However, rMOR-1G1 still could use the first AUG from exon 2 as the translational start codon to yield a 6 TM protein, a situation similar to that human hMOR-1G1 , the mu3 receptor  and hMOR-1K . rMOR-1G2 had the same exon composition as rMOR-1G1 except that the stop codon in exon 11b was skipped due a downstream splice site of exon 11a. Thus, like both mMOR-1G and hMOR-1G2, translation of rMOR-1G2 can proceed from the exon 11a AUG to encode a 6 TM protein since the exon 11a reading-frame was in frame with that of exons 2/3/4 (Figures 2 & 3). The two AUGs in exon 11a were in the same reading-frame and separated by four amino acids. We arbitrarily assigned the first AUG as the translational start codon for rMOR-1G2, although it is not clear which AUG is actually used.
The other five variants, rMOR-1H1, rMOR-1H2, rMOR-1i1, rMOR-1i2 and rMOR-1i3, contained exons 11a, 1a, 2, 3 and 4, but with alternative splicing among exons 11a, 11b, 1a, 1b and 1c to generate different transcripts. Despite their differences in exon composition, rMOR-1H1, rMOR-1i1, rMOR-1i2 and rMOR-1i3 all predicted the same protein sequence as the original rMOR-1 when using AUG in exon 1a as translational start codon (Figure 1B, 2 & 4). Translation from the AUG of exon 11a predicted a short protein sequence due to early translation termination within exon 11b, exons 1b or 1c. The ability of four different exon11-containing transcripts to encoded the same protein as rMOR-1 protein mimics three mouse exon 11-containing variants, mMOR-1H, mMOR-1I and mMOR-1J .
The predicted protein sequence of rMOR-1H2 was intriguing. Splicing from exon 11a to exon 1a gave rise to a sequence that predicted an in-frame fusion protein from exon 11a to exons 1a/2/3/4 when the AUG in exon 11a was used as the translational start codon (Figures 2 & 4). Thus, rMOR-1H2 encoded a novel receptor protein containing the same amino acid sequence as rMOR-1, but with an additional 50 amino acids at the N-terminus. In vitro transcription coupled translation revealed a molecular weight for rMOR-1H2 that was approximately 5 kD higher than that of rMOR-1, suggesting the preferential usage of the AUG in exon 11a to initiate translation (Figure 5). The 50 amino acid sequence did not contain a predicted transmembrane domain, implying that rMOR-1H2 still encoded a 7 TM protein. Interestingly, the additional sequence did possess a potential N-glycosylation site.
Expression of the rat exon 11-associated variant mRNAs
The relative size and abundance of the variant mRNAs was assessed using Northern blot analysis (Figure 6). The exon 2/3a probe, designated to detect most of the variant mRNAs, hybridized several heavy and diffuse bands ranging from 2 - 15 kb, a band pattern similar to Northern blots using mouse and human brains with their respective exon 2/3a probes [24, 37]. The exon 11 probe detected a major strong band around 12 kb. A similar band with relatively same size was also seen in the blot with the exon 2/3a probe. Two weaker bands were seen around 1 - 1.5 kb and 2 - 4.5 kb, respectively.
We next examined the expression of the variant mRNAs in several brain regions using RT-PCR (Figures 7A, 7B, Table 1 & Additional files 1, 2 & 3). The rMOR-1 band was observed in all the regions with relatively equal abundance, except for lower levels in the cerebellum. However, the expression of the other mRNAs varied markedly among the regions. rMOR-1G1, rMOR-1G2 and rMOR-1H2 were highly expressed in the brain stem, hippocampus and spinal cord, but had very lower levels in the cerebellum and hypothalamus. In contrast, rMOR-1H1 was abundant in the cerebellum, hippocampus and spinal cord, but limited in the brain stem and hypothalamus. On the other hand, rMOR-1i1, rMOR-1i2 and rMOR-1i3 expression was mainly observed in the brain stem. The brain stem expressed all the variants at relatively high levels except for rMOR-1H1, whereas the hypothalamus expressed the most variants at very low levels. These results suggested region-specific alternative splicing of these variant pre-mRNAs.
Characterization of the rat exon 11-associated variants by receptor binding
Of the seven exon 11-associated variants, rMOR-1H1, rMOR-1i1, rMOR-1i2 and rMOR-1i3, predicted the same protein as rMOR-1, while rMOR-1H2 encoded a novel receptor protein with additional 50 amino acids extended at the N-terminal tip of rMOR-1. To examine the pharmacological binding profiles of these variants, we established CHO cell lines stably expressing these variants and examined [3H] DAMGO binding. Saturation studies demonstrated similar high affinities of [3H] DAMGO for all five variants (Table 2). Although the small difference in K D values between rMOR-1i3 and rMOR-1, rMOR-1H2, rMOR-1i1 and rMOR-1i2 were statistically significant, we believe that these small differences reflected differences in the assays rather than the receptor itself, particularly since they all predict receptors with identical amino acid sequences. While it theoretically might be due to the concurrent generation of the 16 amino acid fragment predicted from the methionine in exon 11a, there is, to date, nothing to indicate that this peptide is actually generated. Competition studies confirmed their mu selectivity (Table 3), with mu ligands such as morphine and M6G potently lowering binding while the kappa1-selective opioid U50,488H and the delta-selective ligand DPDPE did not. As expected, the variants with the same predicted protein as rMOR-1 displayed the similar binding characteristics as rMOR-1 itself. rMOR-1H2 bound both agonists and antagonists with affinities indistinguishable from other variants including rMOR-1, indicating that the additional 50 amino acids at N-terminus did not influence opioid binding.
Functional comparison of rMOR-1 with rMOR-1H2 in agonist-induced [35S]GTPγS binding
All five full-length variants with 7 TM contained exon 1. Of these, only rMOR-1H2 encoded a novel protein, differing from the others by the extended N-terminal sequence. Previously, we found that the C-terminal variants of the OPRM1 gene displayed differences in agonist-induced G protein activation despite their small differences in receptor binding profile . To investigate a possible functional effect of the additional N-terminal sequence in rMOR-1H2 on agonist-induced G protein activation, we compared the agonist-induce stimulation profiles of several agonists on [35S]GTPγS binding in stably transfected CHO cells expressing rMOR-1H2 and rMOR-1 (Table 4). All the drugs effectively stimulated [35S]GTPγS binding. However, we observed differences in both their potencies (EC50 value) and efficacies (% maximal stimulation) among the variants. For example, DAMGO and dynorphin A were more potent in rMOR-1H2 than in rMOR-1, as was β-endorphin. Maximal stimulation revealed that morphine was significantly more efficacious in rMOR-1H2 than in rMOR-1. There was little correlation between the EC50 and the maximal stimulation, as shown by the fact that dynorphin A was the most efficacious of the ligands tested, despite its lower potency.
To obtain a general indication of the intrinsic activity of the various ligands, we compared their EC50 values after normalizing for their receptor binding affinity (EC50/Ki). This provides an indication of the receptor occupancy needed to elicit the response. The lower the number, the greater is the intrinsic activity of the ligand. While the EC50 values for M6G and DAMGO were similar for rMOR-1, their EC50/Ki ratios differed by approximately 10-fold. A similar situation existed for rMOR-1H2. The potency of dynorphin A in stimulating [35S]GTPγS was far less than that of the other ligands, while its EC50/Ki ratio was the lowest, implying the greatest intrinsic activity. The rank order of the EC50/Ki values varied from their corresponding rank order of both the EC50 and the maximal stimulation values. Comparing the two variants, we also saw different rank-order ratios (Table 4). These results suggested that the extra 50 amino acids influence agonist-induced G protein activation.
Multiple mu opioid receptors were proposed in many years ago, mainly based upon pharmacological studies [3, 43–45]. However, to date only a single mu opioid receptor gene has been identified, raising the questions of how a single OPRM1 gene could explain the complex pharmacology of mu opioids in animals and humans. Our early antisense mapping studies suggested different exon combinations for the analgesic actions of the two mu agonists morphine and M6G, raising the possibility of alternative splicing in the OPRM1 gene [17, 46]. Since then, much effort has been devoted to identifying these OPRM1variants. To date, over 28 splice variants of the mouse OPRM1 gene have been isolated [22–24, 27, 35], some of which had been previously identified in humans  and rats .
The majority of variants were C-terminal variants, differing only at the C-terminal tip. These variants revealed marked differences in their regional distribution at both mRNA and protein level [22–24, 27, 47–52] and agonist-induced G protein activation and internalization [25, 36, 37]. We then identified a second set of variants associated with exon 11, a previously unknown exon located 30 kb upstream of exon 1 , and established their functional significance in an exon 11 KO mouse model . Disrupting exon 11 diminished M6G and heroin analgesia without affecting morphine or methadone actions, suggesting that exon 11 and its associated variants played an important role in the actions of a subset of mu opioids that include M6G and heroin. A number of C-terminal variants have been isolated from the rat [20, 25] and human OPRM1 genes [21, 26, 37]. We recently isolated a homolog exon 11 and three its associated variants in the human OPRM1 gene . The current studies have now extended a similar splicing pattern to the rat with the identification of a homologous exon 11 and seven associated variants in the rat OPRM1 gene. Additionally, the exon 11 sequence has been predicted from the OPRM1 genomic locus of six other mammalian species through NBCI and Ensembl databases, including chimpanzees, monkeys, guinea pigs, bats, cows and armadillos , but not in lower vertebrate species such as fish and amphibians that contain OPRM1 gene.
The nucleotide sequence and genomic location of the rat exon 11 were similar to those in the mouse and human. However, some differences exist between the rat and mouse exons 11. Whereas the rat contains an alternative splice site within exon 11 which splits it into exon 11a and exon 11b, a situation similar to the human exon 11, the mouse does not. Alternative usage of these two splice sites within exon 11, together with a choice of downstream exons, created a number of different variants. The rat exon 11b has a predicted stop codon when translated from its first AUG in exon 11a, leading to only seven amino acids in rMOR-1G1, rMOR-1H1 and rMOR-1i1. However, initiating translation from the first AUG of exon 2 in rMOR-1G1 predicts a 6 TM protein. Using the AUG in exon 1a of rMOR-1H1 and rMOR-1i1, rMOR-1i2, rMOR-1i3 leads to the same protein as the original rMOR-1.
The variants that skip exon 11b can translate through from the AUG in exon 11a, but differed in amino acid sequence depending upon their downstream exons. In rMOR-1i2 and rMOR-1i3, translation using the first AUG in exon 11a still predicted small proteins due to early termination of translation in exons 1b and exon 1c, respectively. However, both rMOR-1i2 and rMOR-1i3 can initiate translation from the AUG of exon 1a to generate the same protein as the original rMOR-1. Thus, together with rMOR-1H1 and rMOR-1i1, a total of four exon 11-associated transcripts can produce the identical rMOR-1 protein, a similar situation seen in the mouse exon 11-aasociated variants. This raises questions regarding why four different splice variants are needed to generate the same protein. It is interesting to speculate, that these differences may differentially regulate their cellular location and their ability to express the protein, but there is no evidence to date to support this possibility.
On the other hand, translation from the AUG of exon 11a also generated the 6 TM protein, rMOR-1G2. Skipping exon 11b maintained the reading frame from AUG of exon 11a through exons 2/3/4. Similarly, skipping exon 11b also enabled rMOR-1H2 to read through, yielding a novel receptor with extra 50 amino acids extended at the N-terminus of rMOR-1, a prediction that was supported by in vitro transcription coupled with translation. Thus, rMOR-1H2 is the first full length (i.e. 7 TM) rat splice variant isolated with a different protein sequence at the N-terminus.
The exon 11 mRNAs are relatively abundant in the brain, as illustrated by Northern blot analysis that displayed a major ~ 12 kb band with intensity comparable to that observed with the exons 2/3 probe. The expression of the exon 11-associated variant mRNAs differed markedly among brain regions, contrasting with the relatively homogenous expression levels of rMOR-1. This suggested that, like the mouse, there is region- and/or cell-specific RNA processing of the variant pre-mRNAs and/or varying levels of upstream promoter activity. Differential expression of the variant mRNAs among brain regions also raised questions regarding their functions. Recently, we observed high correlations between mRNA expression levels, including exon 11-associated variants, in selected brain regions with the degree of morphine and heroin dependence and tolerance among four inbred strains of mice (J Xu, B Kest and YX Pan, unpublished observations). These results suggest a possible contribution of alternative splicing of the OPRM1 gene in mu opioid tolerance and addiction in mice, although the relevance of these correlations needs to be further validated. It will be interesting to see if these correlations also exist in rat.
The genomic location of the rat exon 11 approximately 21 kb upstream of exon 1 suggested the existence of an upstream promoter controlling the expression of the exon 11-associated variants. Preliminary studies indicate that the 5' flanking region of the rat exon 11 has promoter activity, particularly in the neuroblastoma cell lines NIE115 and Be(2)C cells, assessed using a secreted alkaline phosphotase (SEAP) reporter assay (J Xu and XY Pan, unpublished observation).
rMOR-1H2 encoded a full length 7 TM mu opioid receptor with a unique, extended N-terminus. Its similar binding profile is consistent with the other variants was expected since it is believed that the binding pocket is contained within the transmembrane regions, which are identical among all the full length variants. However, the additional N-terminal 50 amino acids in rMOR-1H2 did influence agonist-induced G protein activation, a similar scenario seen in the human N-terminal variant, hMOR-1i . While similar results were observed with the C-terminal variants, but this was more easily understood because of the presumed ability of the C-terminus to influence coupling to transduction proteins. How the additional N-terminal sequence influences receptor activation is as yet unknown.
In the mouse, the exon 11-associated variants mMOR-1G, mMOR-1M and mMOR-1N predict 6 TM variant due to skipping of exon 1, which encodes the first TM. rMOR-1G2 predicted a similar 6 TM protein with translation of exon 11 that resembles mMOR-1G with exon 4 as the last coding exon. While rMOR-1G1 also predicts a 6 TM variant with a terminal exon 4, it requires using the AUG within exon 2 to initiate translation. Despite our efforts, we were unable to isolate rat homologs of the mouse mMOR-1M and mMOR-1N. While it is possible that they do not exist in rats, it also is possible that these homologs are localized to very specific brain regions with a low overall abundance. The functional relevance of the 6 TM mouse variants has been suggested by a range of studies. First, although they do not bind radiolabeled mu agonists with high affinity, the mouse 6 TM variants displayed a moderate binding affinity towards [3H]-diprenorphine (KD approximately 10 nM; J Xu, GW Pasternak and YX Pan, unpublished observation). Second, the 6 TM variants can physically associate with the regular 7 TM MOR-1 and modulate the expression of the 7 TM receptors on cell surface membrane (J Xu, GW Pasternak and YX Pan, unpublished observation). More importantly, disrupting exon 11 diminished M6G and heroin analgesia without affecting morphine and methadone, suggesting selective roles of the 6 TM exon 11-associated variants in the actions of M6G and heroin. Finally, the conservation of the exon 11 and exon 11-associated variants across species further supports their role.
We isolated a rat exon 11 and seven exon 11-associated splice variants from the rat OPRM1 gene, resembling splicing in both mice and humans and suggesting conservation of exon 11 and its associated variants in mammals. The rat OPRM1 gene now contains eleven exons spanning over 250 kb and whose combination by alternative splicing generates over sixteen variants. The functional significance of these rat exon 11-associated variants was suggested by the region-specific expression of their mRNAs and the influence of the novel N-terminal sequence on agonist-induced G protein coupling in the N-terminal variant, rMOR-1H2. The existence of the rat exon 11-associated variants raises questions regarding their potential role in mediating the actions of heroin and M6G in rat. The diversity and complexity created by alternative splicing of the rat OPRM1 gene may provide important insights of understanding the diverse responses to the various mu opioids seen in rat.
Genomic database searching
Alignment of the mouse exon 11 sequence in with the rat OPRM1 gene in the Ensembl human genome database revealed a sequence homologous to exon11. The rat exon 11 was mapped approximately 21 kb upstream of exon 1 in the rat OPRM1 locus. There is 86% identity at the nucleotide level between the rat exon 11 and the mouse exon 11 sequences (Figure 3).
Reverse transcription-polymerase chain reaction (RT-PCR) cloning
Total RNA was isolated from rat brain or selected brain regions by the guanidinium thiocyanate phenol-chloroform extraction method  and reverse transcribed with random primers and Superscript II reverse transcriptase (Invitrogen) as previously described [22, 55] or with an antisense primer from exon 4 (E4-AN1, 5'-CAT GTG CAG AGT GAA GTA GCC AGA G-3') and Superscript III reverse transcriptase (Invitrogen) following the manufacture's protocol. In a 20 μl of RT reaction with random primers and Superscript II, 5 μg of RNA together with 260 ng of random primer was first incubated at 70°C for 5 min and then quickly cooled on ice for 2 min. Following adding the reaction buffer together with 10 mM DTT and 1 mM dNTP and warming the mixture at 37°C for 2 min, 260 units of Superscript II were added. The mixture with the enzyme was incubated at room temperature for 10 min, then at 37°C for 5 min and finally at 42°C for 90 min. The reaction was terminated by heating at 75°C for 15 min. In a 20 μl of RT reaction with E4-AN1 primer and Superscript III, 5 μg of RNA together with 10 pmol of E4-AN1 primer was first incubated at 70°C for 5 min and then quickly cooled on ice. The reaction buffer together with 5 mM DTT, 1 mM dNTP and 200 units of Superscript III was added. The reaction was incubated at 53°C for 90 min and terminated by heating at 75°C for 15 min. Two-step or nested PCRs was used to amplify exon 11-associated full-length clones using Platinum Taq DNA polymerase (Invitrogen). The first-step PCRs were carried out using 5 μl of RT reaction as template with the appropriate primers (see below) for 39 cycles after 2 min at 94°C, each cycle consisting of a 20 sec denaturing step at 94°C, a 20 sec annealing step at 65°C and a 2 min extension at 72°C. In the second-step PCRs, 2 μl of the first-step PCR products was used as template with appropriate primers (see below) using the same PCR cycling conditions as the first-step PCR. The primers used for rMOR-1G1, rMOR-1G2, and rMOR-1H2 were: two sense primers from exon 11 sequence obtained from the genomic alignment (E11-SE1: 5'-CTT CCC ATA AGT CAT TTG CTG TCC TTG-3' and E11-SE2: 5'-GAA GAG GAA CAC CGA AAC TGG GAA GC-3') and two antisense primers from exon 4 (E4-AN1 and E4-AN2: 5'-GAC AGC AAC CTG ATT CCA CGT AGA TG-3'); for rMOR-1H1 and rMOR-1i1, an exon 11 sense primer (E11-SE3: 5'-GAA GGA TGG GAT CTG GTC CAA TGC TGT AAG CTT TCT CCA AGT CCG CAT TCC AAA AAC TGG ACA GGG AGA TAG AAA TCA AGA GGG GAA GTT ACC TCA G-3') and the two exon 4 antisense primers (E4-AN1 and E4-AN2); for rMOR-1i2, a exon 11 sense primer (E11-SE4: 5'-GAA GGA TGG GAT CTG GTC CAA TGC TTG CAT G-3') and E4-AN1 primer; and for rMOR-1i3, an exon 11 sense primer (E11-SE5: 5'-GCT TGA AGG ATG GGA TCT GGT CCA ATG CTA TAC GAG-3') and E4-AN1 and E4-AN2 primers. All PCR fragments were subcloned into pcDNA3.1/V5His-TOPO vector (Invitrogen) and sequenced with appropriate primers in both orientations.
Northern blot analysis
Northern blot analysis was performed as described [22, 55]. Briefly, 20 μg of total brain RNA/lane was separated on a 0.8% formaldehyde agarose gel, and transferred to GenePlus membrane. The membranes were hybridized with either a 257 bp 32P-labeled exon 11 probe generated by PCR with a sense primer (E11-SE1) and an antisense primer (E11-AN1: 5'-GAG GTA ACT TCC CCT CTT GAT TTC TAT CTC CC-3') from exon 11 or a 685 bp 32P-labeled exons 2 & 3a probe by PCR with a sense primer from exon 2 (E2-SE1: 5'-GAC TGC CAC CAA CAT CTA CAT TTT CAA C-3') and an antisense primer from exon 3 (E3-AN1: 5'-GTT CGT GTA ACC CAA AGC AAT GC-3').
Regional expression of the variant mRNAs
The selected brain regions were dissected from four separate groups of rats. Each group had one or two rats depending upon the size of the region. Total RNAs extracted from the selected brain regions were reverse-transcribed with an E4-AN1 primer and Superscript III as described in RT-PCR cloning (see above). The first-strand cDNAs were used as templates for two-step or nested PCRs. For exon 11-associated splice variants, the first-step PCRs were performed using 5 μl of RT reaction as template and a sense primer from exon 11 (E11-SE1) and an antisense primer from exon 4 (E4-AN1) for 35 cycles with the same PCR cycling conditions as described in RT-PCR cloning. In the second-step PCRs, 3 μl of the first-step PCR products was used as template with appropriate primers (see below). The primers used in the second-step PCRs were designed to specifically amplify each variant and listed as following: for rMOR-1G1, a sense primer (G1-SE: 5'-GAA GTT ACC TCA GAT ACA CCA AAA TGA-3') and an exon 3 antisense primer (E3-AN2: CAG CAG ACG ATA AAT ACA GCC ACG-3'); for rMOR-1G2, a sense primer (G2-SE: 5'-GGT CCA ATG CTA TAC ACC AAA ATG-3') and E3-AN2 primer; for rMOR-1H1, a sense primer (H1-SE: 5'-AAG TTA CCT CAG GGC TGG TCC-3') and an exon 2 antisense primer (E2-AN1: 5'-ATG TTC CCA TCA GGT AGT TGA CAC TC-3'); for rMOR-1H2, a sense primer (H2-SE: 5'-GGT CCA ATG CTG GCT GGT CC-3') and E2-AN1 primer; for rMOR-1i1, a sense primer (I1-SE: 5'-AAG TTA CCT CAG TGC ATG GAG ACC-3') and E2-AN1 primer; for rMOR-1i2, a sense primer (I2-SE: GGT CCA ATG CTT GCA TGG AGA C-3') and E2-AN1 primer; for rMOR-1i3, a sense primer (I3-SE: 5'-GGT CCA ATG CTA TAC GCG GA-3') and E2-AN1 primer. For detecting rMOR-1, the first-step PCRs were performed using the same PCR cycling conditions with 5 μl of RT reaction, an exon 1c sense primer (E1-SE1: 5'-CCC ACT TTA CAC TCG TTT ACA CGG-3') and E4-AN1 primer. In the second-step PCRs, 3 μl of the first-step PCR products was used as template with an exon 1 sense primer (E1-SE2: 5'-GAC AGC CTG TGC CCT CAG ACC-3') and E2-AN1 primer. The second-step PCRs were carried out for 35 cycles after 2 min at 94°C, each cycle consisting of a 20 sec denaturing step at 94°C, a 20 sec annealing step at 60 - 65°C and 45 - 90 sec extension at 72°C, depending upon melting temperature of primers and length of amplicons. The lengths of the PCR products were consistent with their predicted sizes: 606 bp for rMOR-1G1, 604 bp for rMOR-1G2, 539 bp for rMOR-1H1, 537 bp for rMOR-1H2, 794 bp for rMOR-1i1, 793 bp for rMOR-1i2, 1051 bp for rMOR-1i3 and 217 bp for rMOR-1. The sequences of the PCR products were confirmed by sequencing with appropriate primers. A negative control using ddH2O as template was included for each variant throughout the two-step PCRs. RNA loading was estimated by parallel one-step PCRs with a pair of primers for glyceraldehydes 3-phosphate dehydrogenase (G3PDH) (Clontech). The PCR products were separated on 1% agarose gel, stained with ethidium bromide. The agarose gel was photographed and analyzed using a FluorChem 8000 Image System (Alpha Innotech).
In vitro transcription coupled translation
The full-length cDNAs of rMOR-1 and rMOR-1H2 in the pcDNA3.1/V5His-TOPO vector were transcribed and translated in vitro with a TnT T7 coupled reticulocyte lysate system (Promega) following the manufacturer's protocol. Briefly, the plasmids were incubated with T7 RNA polymerase and reticulocyte lysate in the presence of 0.04 mCi of [35S]methionine (> 1000 Ci/mmol; PerkimElmer) at 25°C for 90 min. The translated products were separated on a 12% SDS-polyacrylamide gel, and the gel was treated with Amplify (GE Life), dried and exposed to Kodak BioMax MR film.
Expression of rMOR-1H1, rMOR-1H2, rMOR-1i1, rMOR-1i2 and rMOR-1i3 in Chinese hamster ovary (CHO) cells
The rMOR-1H1/pcDNA3.1-TOTO, rMOR-1H2/pcDNA3.1-TOPO, rMOR-1i1/pcDNA3.1-TOPO, rMOR-1i2/pcDNA3.1-TOPO, rMOR-1i3/pcDNA3.1-TOPO and rMOR-1/pcDNA3.1(-) plasmids were used to transfect Chinese Hamster Ovary (CHO) cells by LipofectAMINE reagent (Invitrogen). Stable transformants were obtained 10 - 14 days after selection with G418 and screened with a [3H]DAMGO binding assay.
Receptor binding assays
Membranes were prepared from stable transfectants as described previously . Saturation and competition binding assays were performed with [3H]DAMGO at 25°C for 60 min in 50 mM potassium phosphate buffer, pH 7.4, containing 5 mM magnesium sulfate. Specific binding was defined as the difference between total binding and non-specific binding, determined in the presence of 10 μM levallorphan. K D and K i values were calculated by non-linear regression analysis (GraphPad Prism 4.0, Carlsbad, CA). Protein concentrations were determined using the Lowry method as previously described using bovine serum albumin (BSA) as the standard [22, 28].
[35S]GTPγS binding assay
Membranes prepared from stable transfectants were incubated in the presence and absence of indicated opioids for 60 min at 30°C in the assay buffer (50 mM Tris-HCl, pH 7.7, 3 mM MgCl2, 0.2 mM EGTA, 10 mM NaCl) containing 0.05 nM [35S]GTPγS (> 1000 Ci/mmol, PerkimElmer) and 60 μM GDP, as previously reported [25, 36, 37]. Basal binding was determined in the presence of GDP and absence of drug. The reaction was terminated by rapid filtration under vacuum through glass fiber filters, followed by three washes with 3 ml of ice-cold 50 mM Tris-HCl, pH 7.4. Bound radioactivity was measured by liquid scintillation spectrophotometry in Liquid Scintillation Analyzer (TRI-CARB 2900TR, PerkimElmer) after overnight extraction in 5 ml liquiscint scintillation fluid (National Diagnostic Inc.).
mu opioid receptor
mu opioid receptor gene
polymerase chain reaction
Lord JAH, Waterfield AA, Hughes J, Kosterlitz HW: Endogenous opioid peptides: multiple agonists and receptors. Nature 1977, 267: 495–499. 10.1038/267495a0
Martin WR, Eades CG, Thompson JA, Huppler RE, Gilbert PE: The effects of morphine and nalorphine-like drugs in the nondependent and morphine-dependent chronic spinal dog. J Pharmacol Exp Ther 1976, 197: 517–532.
Wolozin BL, Pasternak GW: Classification of multiple morphine and enkephalin binding sites in the central nervous system. Proc Natl Acad Sci USA 1981, 78: 6181–6185. 10.1073/pnas.78.10.6181
Chen Y, Mestek A, Liu J, Hurley JA, Yu L: Molecular cloning and functional expression of a μ-opioid receptor from rat brain. Mol Pharmacol 1993, 44: 8–12.
Thompson RC, Mansour A, Akil H, Watson SJ: Cloning and pharmacological characterization of a rat μ opioid receptor. Neuron 1993, 11: 903–913. 10.1016/0896-6273(93)90120-G
Wang JB, Imai Y, Eppler CM, Gregor P, Spivak CE, Uhl GR: μ opiate receptor: cDNA cloning and expression. Proc Natl Acad Sci USA 1993, 90: 10230–10234. 10.1073/pnas.90.21.10230
Guiramand J, Montmayeur JP, Ceraline J, Bhatia M, Borrelli E: Alternative splicing of the dopamine D2 receptor directs specificity of coupling to G-proteins. J Biol Chem 1995, 270: 7354–7358. 10.1074/jbc.270.13.7354
Usiello A, Baik JH, Rouge-Pont F, Picetti R, Dierich A, LeMeur M, et al.: Distinct functions of the two isoforms of dopamine D2 receptors. Nature 2000, 408: 199–203. 10.1038/35041572
Vanetti M, Kouba M, Wang X, Vogt G, Hollt V: Cloning and expression of a novel mouse somatostatin receptor (SSTR2B). FEBS Lett 1998, 311: 290–294. 10.1016/0014-5793(92)81122-3
Namba T, Sugimoto Y, Negishi M, Irie A, Ushikubi F, Kakizuka A, et al.: Alternative splicing of C-terminal tail of prostaglandin E receptor subtype EP3 determines G-protein specificity. Nature 1993, 365: 166–170. 10.1038/365166a0
Lucas JJ, Hen R: New players in the 5-HT receptor field: genes and knockouts. Trends Pharmacol Sci 1995, 16: 246–252. 10.1016/S0165-6147(00)89034-3
Candenas ML, Cintado CG, Pennefather JN, Pereda MT, Loizaga JM, Maggi CA, et al.: Identification of a tachykinin NK(2) receptor splice variant and its expression in human and rat tissues. Life Sci 2002, 72: 269–277. 10.1016/S0024-3205(02)02240-3
Schulz HL, Stohr H, Weber BH: Characterization of three novel isoforms of the metabotrobic glutamate receptor 7 (GRM7). Neurosci Lett 2002, 326: 37–40. 10.1016/S0304-3940(02)00306-3
Malherbe P, Kew JN, Richards JG, Knoflach F, Kratzeisen C, Zenner MT, et al.: Identification and characterization of a novel splice variant of the metabotropic glutamate receptor 5 gene in human hippocampus and cerebellum. Brain Res Mol Brain Res 2002, 109: 168–178. 10.1016/S0169-328X(02)00557-0
Park YS, Lee YS, Cho NJ, Kaang BK: Alternative splicing of gar-1, a Caenorhabditis elegans G-protein-linked acetylcholine receptor gene. Biochem Biophys Res Commun 2000, 268: 354–358. 10.1006/bbrc.2000.2108
Park YS, Kim S, Shin Y, Choi B, Cho NJ: Alternative splicing of the muscarinic acetylcholine receptor GAR-3 in Caenorhabditis elegans. Biochem Biophys Res Commun 2003, 308: 961–965. 10.1016/S0006-291X(03)01508-0
Rossi GC, Pan YX, Brown GP, Pasternak GW: Antisense mapping the MOR-1 opioid receptor: Evidence for alternative splicing and a novel morphine-6β-glucuronide receptor. FEBS Lett 1995, 369: 192–196. 10.1016/0014-5793(95)00757-Z
Pasternak GW: Multiple opiate receptors: deja vu all over again. Neuropharmacology 2004,47(Suppl 1):312–323. 10.1016/j.neuropharm.2004.07.004
Schuller AG, King MA, Zhang J, Bolan E, Pan YX, Morgan DJ, et al.: Retention of heroin and morphine-6 beta-glucuronide analgesia in a new line of mice lacking exon 1 of MOR-1. Nat Neurosci 1999, 2: 151–156. 10.1038/5706
Zimprich A, Simon T, Hollt V: Cloning and expression of an isoform of the rat μ opioid receptor (rMOR 1 B) which differs in agonist induced desensitization from rMOR1. FEBS Lett 1995, 359: 142–146. 10.1016/0014-5793(95)00028-8
Bare LA, Mansson E, Yang D: Expression of two variants of the human μ opioid receptor mRNA in SK-N-SH cells and human brain. FEBS Lett 1994, 354: 213–216. 10.1016/0014-5793(94)01129-X
Pan YX, Xu J, Bolan E, Abbadie C, Chang A, Zuckerman A, et al.: Identification and characterization of three new alternatively spliced mu-opioid receptor isoforms. Mol Pharmacol 1999, 56: 396–403.
Pan YX, Xu J, Bolan E, Chang A, Mahurter L, Rossi G, et al.: Isolation and expression of a novel alternatively spliced mu opioid receptor isoform, MOR-1F. FEBS Lett 2000, 466: 337–340. 10.1016/S0014-5793(00)01095-4
Pan YX, Xu J, Mahurter L, Bolan E, Xu M, Pasternak GW: Generation of the mu opioid receptor (MOR-1) protein by three new splice variants of the Oprm gene. Proc Natl Acad Sci USA 2001, 98: 14084–14089. 10.1073/pnas.241296098
Pasternak DA, Pan L, Xu J, Yu R, Xu MM, Pasternak GW, et al.: Identification of three new alternatively spliced variants of the rat mu opioid receptor gene: dissociation of affinity and efficacy. J Neurochem 2004, 91: 881–890. 10.1111/j.1471-4159.2004.02767.x
Pan YX, Xu J, Mahurter L, Xu MM, Gilbert AK, Pasternak GW: Identification and characterization of two new human mu opioid receptor splice variants, hMOR-1O and hMOR-1X. Biochem Biophys Res Commun 2003, 301: 1057–1061. 10.1016/S0006-291X(03)00089-5
Pan YX, Xu J, Bolan E, Moskowitz HS, Xu M, Pasternak GW: Identification of four novel exon 5 splice variants of the mouse mu-opioid receptor gene: functional consequences of C-terminal splicing. Mol Pharmacol 2005, 68: 866–875.
Xu J, Xu M, Hurd YL, Pasternak GW, Pan YX: Isolation and characterization of new exon 11-associated N-terminal splice variants of the human mu opioid receptor gene. J Neurochem 2009, 108: 962–972. 10.1111/j.1471-4159.2008.05833.x
Doyle GA, Rebecca SX, Lin SS, Press DM, Grice DE, Buono RJ, et al.: Identification of three mouse micro-opioid receptor (MOR) gene (Oprm1) splice variants containing a newly identified alternatively spliced exon. Gene 2007, 388: 135–147. 10.1016/j.gene.2006.10.017
Doyle GA, Sheng XR, Lin SS, Press DM, Grice DE, Buono RJ, et al.: Identification of five mouse mu-opioid receptor (MOR) gene (Oprm1) splice variants containing a newly identified alternatively spliced exon. Gene 2007, 395: 98–107. 10.1016/j.gene.2007.02.004
Cadet P, Mantione KJ, Stefano GB: Molecular identification and functional expression of mu 3, a novel alternatively spliced variant of the human mu opiate receptor gene. J Immunol 2003, 170: 5118–5123.
Choi HS, Kim CS, Hwang CK, Song KY, Wang W, Qiu Y, et al.: The opioid ligand binding of human mu-opioid receptor is modulated by novel splice variants of the receptor. Biochem Biophys Res Commun 2006, 343: 1132–1140. 10.1016/j.bbrc.2006.03.084
Shabalina SA, Zaykin DV, Gris P, Ogurtsov AY, Gauthier J, Shibata K, et al.: Expansion of the human mu-opioid receptor gene architecture: novel functional variants. Hum Mol Genet 2009, 18: 1037–1051. 10.1093/hmg/ddn439
Kvam TM, Baar C, Rakvag TT, Kaasa S, Krokan HE, Skorpen F: Genetic analysis of the murine mu opioid receptor: increased complexity of Oprm gene splicing. J Mol Med 2004, 82: 250–255. 10.1007/s00109-003-0514-z
Pan YX: Diversity and complexity of the mu opioid receptor gene: alternative pre-mRNA splicing and promoters. DNA Cell Biol 2005, 24: 736–750. 10.1089/dna.2005.24.736
Bolan EA, Pan YX, Pasternak GW: Functional analysis of MOR-1 splice variants of the mouse mu opioid receptor gene Oprm. Synapse 2004, 51: 11–18. 10.1002/syn.10277
Pan L, Xu J, Yu R, Xu MM, Pan YX, Pasternak GW: Identification and characterization of six new alternatively spliced variants of the human mu opioid receptor gene, Oprm. Neuroscience 2005, 133: 209–220. 10.1016/j.neuroscience.2004.12.033
Abbadie C, Pasternak GW: Differential in vivo internalization of MOR-1 and MOR-1C by morphine. Neuroreport 2001, 12: 3069–3072. 10.1097/00001756-200110080-00017
Zhang J, Ferguson SS, Barak LS, Bodduluri SR, Laporte SA, Law PY, et al.: Role for G protein-coupled receptor kinase in agonist-specific regulation of mu-opioid receptor responsiveness. Proc Natl Acad Sci USA 1998, 95: 7157–7162. 10.1073/pnas.95.12.7157
Qiu Y, Law PY, Loh HH: Mu-opioid receptor desensitization: role of receptor phosphorylation, internalization, and representation. J Biol Chem 2003, 278: 36733–36739. 10.1074/jbc.M305857200
Bohn LM, Gainetdinov RR, Lin FT, Lefkowitz RJ, Caron MG: Mu-opioid receptor desensitization by beta-arrestin-2 determines morphine tolerance but not dependence. Nature 2000, 408: 720–723. 10.1038/35047086
Pan YX, Xu J, Xu M, Rossi GC, Matulonis JE, Pasternak GW: Involvement of exon 11-associated variants of the mu opioid receptor MOR-1 in heroin, but not morphine, actions. PNAS 2009, 106: 4917–4922. 10.1073/pnas.0811586106
Pasternak GW: Pharmacological mechanisms of opioid analgesics. Clin Neuropharmacol 1993, 16: 1–18. 10.1097/00002826-199302000-00001
Pasternak GW, Childers SR, Snyder SH: Opiate analgesia: Evidence for mediation by a subpopulation of opiate receptors. Science 1980, 208: 514–516. 10.1126/science.6245448
Pasternak GW, Snyder SH: Identification of a novel high affinity opiate receptor binding in rat brain. Nature 1975, 253: 563–565. 10.1038/253563a0
Rossi GC, Pan YX, Cheng J, Pasternak GW: Blockade of morphine analgesia by an antisense oligodeoxynucleotide against the mu receptor. Life Sci 1994, 54: PL375–379. 10.1016/0024-3205(94)90038-8
Abbadie C, Pan YX, Pasternak GW: Differential distribution in rat brain of mu opioid receptor carboxy terminal splice variants MOR-1C and MOR-1-like immunoreactivity: Evidence for region-specific processing. J Comp Neurol 2000, 419: 244–256. 10.1002/(SICI)1096-9861(20000403)419:2<244::AID-CNE8>3.0.CO;2-R
Abbadie C, Pan YX, Drake CT, Pasternak GW: Comparative immunhistochemical distributions of carboxy terminus epitopes from the mu opioid receptor splice variants MOR-1D, MOR-1 and MOR-1C in the mouse and rat central nervous systems. Neuroscience 2000, 100: 141–153. 10.1016/S0306-4522(00)00248-7
Abbadie C, Gultekin SH, Pasternak GW: Immunohistochemical localization of the carboxy terminus of the novel mu opioid receptor splice variant MOR-1C within the human spinal cord. Neuroreport 2000, 11: 1953–1957. 10.1097/00001756-200006260-00029
Abbadie C, Pasternak GW, Aicher SA: Presynaptic localization of the carboxy-terminus epitopes of the mu opioid receptor splice variants MOR-1C and MOR-1D in the superficial laminae of the rat spinal cord. Neuroscience 2001, 106: 833–842. 10.1016/S0306-4522(01)00317-7
Zhang Y, Pan YX, Kolesnikov Y, Pasternak GW: Immunohistochemical labeling of the mu opioid receptor carboxy terminal splice variant mMOR-1B4 in the mouse central nervous system. Brain Res 2006, 1099: 33–43. 10.1016/j.brainres.2006.04.133
Abbadie C, Rossi GC, Orciuolo A, Zadina JE, Pasternak GW: Anatomical and functional correlation of the endomorphins with mu opioid receptor splice variants. Eur J Neurosci 2002, 16: 1075–1082. 10.1046/j.1460-9568.2002.02173.x
Pan YX, Pasternak GW: Molecular Biology of Mu Opioid Receptors. In Opiate. Edited by: Pasternak GW. Humana Press; 2010. Ch.6.
Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 1987, 162: 156–159. 10.1016/0003-2697(87)90021-2
Pan YX, Cheng J, Xu J, Rossi G, Jacobson E, Ryan-Moro J, et al.: Cloning and fuctional characterization of a kappa 3 -related opioid receptor. Mol Pharmacol 1995, 47: 1180–1188.
This work was supported, in part, by research grants to GWP (DA02615 and DA00220) and Y-XP (DA13997 and DA02944) from the National Institute on Drug Abuse, a grant from the National Genetics Foundation, a grant from Mr. William H. Goodwin and Mrs. Alice Goodwin and the Commonwealth Foundation for Cancer Research and The Experimental Therapeutics Center of Memorial Sloan-Kettering Cancer Center to GWP and by a Core Grant to MSKCC from the National Cancer Institute (CA8748).
The authors declare that they have no competing interests.
JX and MX performed RT-PCR, Northern blot analysis, in vitro transcription coupled translation, receptor binding and [35S]GTPγS binding assays. GCR dissected brain regions for RNA extraction. Y-XP designed and coordinated the experiments and analyzed data. Y-XP and GWP oversaw the project and wrote the manuscript. All authors read and approved the manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Xu, J., Xu, M., Rossi, G.C. et al. Identification and characterization of seven new exon 11-associated splice variants of the rat mu opioid receptor gene, OPRM1. Mol Pain 7, 9 (2011). https://doi.org/10.1186/1744-8069-7-9
- Translational Start Codon
- Variant mRNAs
- Select Brain Region