Abstract
In recent years, the distribution of dopamine receptor subtypes among the principal neurons of the neostriatum has been the subject of debate. Conventional anatomical and physiological approaches have yielded starkly different estimates of the extent to which D1 and D2 class dopamine receptors are colocalized. One plausible explanation for the discrepancy is that some dopamine receptors are present in physiologically significant numbers, but the mRNA for these receptors is not detectable with conventional techniques. To test this hypothesis, we examined the expression of DA receptors in individual neostriatal neurons by patch-clamp and RT-PCR techniques. Because of the strong correlation between peptide expression and projection site, medium spiny neurons were divided into three groups on the basis of expression of mRNA for enkephalin (ENK) and substance P (SP). Neurons expressing detectable levels of SP but not ENK had abundant mRNA for the D1a receptor. A subset of these cells (∼50%) coexpressed D3 or D4receptor mRNA. Neurons expressing detectable levels of ENK but not SP had abundant mRNA for D2 receptor isoforms (short and long). A subset (10–25%) of these neurons coexpressed D1aor D1b mRNAs. Neurons coexpressing ENK and SP mRNAs consistently coexpressed D1a and D2 mRNAs in relatively high abundance. Functional analysis of neurons expressing lower abundance mRNAs revealed clear physiological consequences that could be attributed to these receptors. These results suggest that, although colocalization of D1a and D2 receptors is limited, functional D1 and D2 class receptors are colocalized in nearly one-half of all medium spiny projection neurons.
The signaling pathways activated by dopamine in the neostriatum have been the subject of intense study since it was discovered that the loss of dopamine leads to the psychomotor symptoms of Parkinson’s disease (Hornykiewcz, 1973). Subsequently, several other common psychomotor disorders, including schizophrenia and Tourette’s syndrome, have been linked to alterations in neostriatal dopaminergic signaling (Nemeroff and Bissette, 1988; Erenberg, 1992).
In recent years, significant progress has been made in characterizing the membrane receptors transducing the signals of dopamine in the neostriatum and the brain in general. Molecular cloning work has shown that there are at least five dopamine receptor genes (D1a, D2, D3, D4, D1b;Sibley, 1995). These receptors can be grouped into D1(D1a, D1b) and D2 classes (D2, D3, D4) on the basis of their affinities for classical dopaminergic ligands and their coupling to intracellular signaling cascades.
How these receptors are distributed among the principal neuronal cell types in the neostriatum has been the subject of debate (Surmeier et al., 1993). This controversy stems primarily from discrepancies in the results obtained from functional studies on the one hand and anatomical studies on the other. The most compelling anatomical data are in situ hybridization studies suggesting that D1a and D2 mRNA are segregated primarily in the two major efferent neostriatal populations (Gerfen, 1992; LeMoine and Bloch, 1995). In particular, D1a receptor mRNA is found in substance P-expressing (SP) neurons projecting to the substantia nigra, whereas D2 receptor mRNA is found in enkephalin-expressing (ENK) neurons projecting exclusively to the globus pallidus. More recent immunocytochemical work supports this conclusion (Hersch et al., 1995), although others have reported significant degrees of receptor protein colocalization (Ariano et al., 1995).
Functional studies, on the other hand, repeatedly have observed responses to D1 and D2 class agonists that are difficult to explain if these receptor classes are not colocalized (Uchimura et al., 1986; Akaike et al., 1987; Cepeda et al., 1993; for review, see Surmeier et al., 1993). The most compelling evidence comes from patch-clamp studies of acutely isolated neostriatal neurons—where synaptic interactions have been removed—showing neuromodulatory effects of both D1 and D2 class agonists in the same cell (Surmeier et al., 1992). In this study, it was also shown that neurons projecting axons to the substantia nigra coexpressed D1a and D2 receptor mRNAs. However, the sample of neurons in which the mRNA analysis was performed was small, making it unclear whether the molecular identification of the receptors mediating the physiological responses could be generalized.
An implicit assumption in this debate has been that the presence of D1a and D2 receptors can be inferred from the responses to D1 and D2 classligands. The underpinnings for this line of argument are the apparent low level of D3, D4, and D1b mRNA expression in the dorsal neostriatum. In the experiments described below, the expression of dopamine receptor mRNAs in medium spiny neurons is reexamined with a newly developed technique—single-cell RT-PCR—that obviates some of the difficulties inherent in more conventional molecular techniques. Of particular interest to us was how dopamine receptor mRNA expression is coordinated in individual, identified neurons. Functional assays for particular receptor mRNAs were performed by applying whole-cell voltage-clamp techniques to neurons subjected to PCR-based analysis. Our results show that, although D1a and D2 mRNAs are coexpressed in only a subset of neurons, the presence of other dopamine receptor mRNAs (D3, D4, D1b) results in the coexpression of D1 and D2 class mRNAs in nearly one-half of all medium spiny neurons. Taken together, our results provide a framework within which previous anatomical and functional studies of dopaminergic signaling in the neostriatum can be reconciled.
MATERIALS AND METHODS
Acute dissociation procedure. Neostriatal neurons from young adult (>3 wks) rats were acutely dissociated by procedures similar to those we have described previously (Bargas et al., 1994;Surmeier et al., 1995a). All solutions coming into contact with the tissue during the dissociation procedure were prepared aseptically and nominally RNase-free. In brief, rats were anesthetized with methoxyflurane and decapitated; brains were removed quickly, iced, and then blocked for slicing. The blocked tissue was cut in 400 μm slices with a Vibroslice (Campden Instruments, London, UK) while bathed in a low Ca2+ (100 μm) HEPES-buffered salt solution (in mm): 140 Na isethionate, 2 KCl, 4 MgCl2, 0.1 CaCl2, 23 glucose, and 15 HEPES, pH 7.4, 300–305 mOsm/l. Then slices were incubated for 1–6 hr at room temperature (20–22°C) in an NaHCO3-buffered saline bubbled with 95% O2/5% CO2 (in mm): 126 NaCl, 2.5 KCl, 2 CaCl2, 2 MgCl2, 26 NaHCO3, 1.25 NaH2PO4, 1 pyruvic acid, and 10 glucose, pH 7.4, with NaOH, 300–305 mOsm/l. In later experiments ascorbic acid (200 μm),NG-nitro-l-arginine (100 μm), and kynurenic acid (1 mm) were added to the holding solutions, because their addition improved the apparent viability of dissociated neurons. All reagents were obtained from Sigma Chemical (St. Louis, MO). Then slices were removed into the low Ca2+ buffer, and, with the aid of a dissecting microscope, regions of the dorsal neostriatum were dissected and placed in an oxygenated Cell-Stir chamber (Wheaton, Millville, NJ) containing Pronase (1–3 mg/ml) in HEPES-buffered HBSS (Sigma Chemical) at 35°C. Dissections were limited to tissue rostral to the anterior commissure to reduce the possibility of contamination from pallidum. After 20–30 min of enzyme digestion, tissue was rinsed three times in the low Ca2+, HEPES-buffered saline and mechanically dissociated with a graded series of fire-polished Pasteur pipettes. Then the cell suspension was plated into a 35 mm Lux Petri dish mounted on the stage of an inverted microscope containing HEPES-buffered HBSS saline. After allowing the cells to settle, the solution bathing the cells was changed to our normal recording external solution.
Whole-cell recordings. Whole-cell recordings used standard techniques (Hamill et al., 1981; Bargas et al., 1994). Electrodes were pulled from Corning 7052 glass and fire-polished before use. The internal solution consisted of (in mm): 180 N-methyl-d-glucamine (NMG), 40 HEPES, 4 MgCl2, 0.1 BAPTA, 12 phosphocreatine, 2 Na2ATP, 0.2 Na3GTP, and 0.1 leupeptin, pH 7.2–7.3, with H2SO4, 265–270 mOsm/l. The pH of NMG solutions was measured with a Corning Model 476570 probe. The external solution consisted of (in mm): 135 NaCl, 20 CsCl, 1 MgCl2, 10 HEPES, 0.001 TTX, 5 BaCl2, 10 glucose, pH 7.3, with NaOH, 300–305 mOsm/l. All reagents were obtained from Sigma Chemical except ATP and GTP (Boehringer Mannheim, Indianapolis, IN), BAPTA and leupeptin (Calbiochem, La Jolla, CA), and sulfuric acid (Fluka, Ronkonkoma, NY).
Recordings were obtained with an Axon Instruments 200 patch-clamp amplifier and controlled and monitored with a PC 486 clone running pCLAMP (v. 6.0) with a 125 kHz interface (Axon Instruments, Foster City, CA). Electrode resistances were typically 3–6 mΩ in the bath. After seal rupture, series resistance (7–15 MΩ) was compensated (80–90%) and monitored periodically. Because whole-cell currents never exceeded 2 nA (usually <600 pA), errors in voltage attributable to inadequate compensation never should have exceeded a few mV. The adequacy of voltage control was assessed after compensation by examining the tail currents generated by strong depolarizations. Cells in which tail currents did not decay rapidly and smoothly at subthreshold potentials were discarded.
The D1 class agonist (+)-6-chloro-7,8-dihydroxy-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine hydrobromide (6-chloro-PB or SKF 81297) and the D2 class agonisttrans-(−)-4aR-4,4a,5,6,7,8,8a,9-octahydro-5-propyl-1H-pyrazolo[3,4-g]quinoline hydrochloride (quinpirole or LY-171,555) and 5-hydroxytryptamine hydrochloride (5-HT) were made up as concentrated stocks in deoxygenated water containing 0.1% ascorbic acid or sodium metabisulfite. Stocks were thawed and diluted immediately before use; solutions were protected from ambient light and contained 0.01% metabisulfite. Drugs were applied with a gravity-fed “sewer pipe” system. The array of application capillaries (∼150 μm inner diameter) was positioned a few hundred microns from the cell under study. Solution changes were effected by altering the position of the array with a DC drive system controlled by a microprocessor-based controller (Newport-Klinger, Irvine, CA). Solution changes were complete within <1 sec; the time constant of Cd2+ block of Ca2+ current was 400 msec. Potentials were not corrected for the liquid junction potential with the NMG-SO4 internal solution, which was estimated to be 7 mV.
Poisson statistics were used to estimate blank probability as a function of cDNA abundance. After Feller’s (1950) treatment of the distribution of r balls in n boxes, the blank probability (PB) is given by: in which λ = n · exp(−r/n), r = copy number (× PCR threshold), n = (fraction of total cellular cDNA)−1. Simulations were run with Igor Pro (Wavemetrics, Lake Oswego, OR).
Single-neuron RT-PCR analysis. After recording, cells were lifted into a stream of control solution and aspirated into the electrode by applying negative pressure. External solutions were made with nominally RNase-free water (Milli-Q PF, Millipore, Bedford, MA). Electrodes contained ∼5 μl of sterile recording solution (see above). Some cells were harvested without recording with electrodes filled with water. The capillary glass used for making electrodes had been heated to 200°C for 4 hr. Sterile gloves were worn during the procedure to minimize RNase contamination.
After aspiration, the electrode was broken and the contents were ejected into a 0.5 ml Eppendorf tube containing 5 μl of diethyl pyrocarbonate (DEPC)-treated water, 0.5 μl RNAsin (28,000 U/ml), and 0.5 μl dithiothreitol (DTT; 0.1 m). Either oligo dT (0.5 μg/ml or random hexanucleotides (50 ng/μl) primers were added, and the mixture was heated to 70°C for 10 min and then incubated on ice for >1 min. Single-strand cDNA was synthesized from the cellular mRNA by adding SuperScript II reverse transcriptase (1 μl, 200 U/μl) and buffer [4 μl, 5× first-strand buffer containing (in mm): 250 Tris-HCl, 375 KCl, and 15 MgCl2], RNAsin (0.5 μl, 28,000 U/ml), DTT (1.5 μl, 0.1 m), and mixed dNTPs (1 μl, 10 mm). The reaction mixture (20 μl) was incubated at 42°C for 50 min. The reaction was terminated by heating the mixture to 70°C for 15 min and then icing. Then the RNA strand in the RNA–DNA hybrid was removed by adding 1 μl of RNase H (2 U/μl) and incubating for 20 min at 37°C. All reagents except for RNAsin (Promega, Madison, WI) were obtained from Life Technologies (Grand Island, NY). No attempt was made to quantify the cDNA from individual cells. The cDNA from the reverse transcription (RT) of RNA in single neostriatal neurons was amplified with the polymerase chain reaction (PCR).
PCR amplification was performed with a thermal cycler (MJ Research, Watertown, MA) with thin-walled plastic tubes. Amplification was performed with one of two protocols. In our initial experiments, a fraction of the single-cell cDNA was used as a template for conventional PCR amplification. Reaction mixtures contained 2–2.5 mm MgCl2, 0.5 mm of each of the dNTPs, 0.8–1 μm primers, 2.5 U Taq DNA polymerase (Promega), 5 μl of 10× buffer (Promega), and one-tenth to one-fifth (2–4 μl) of the cDNA template made from the single-cell RT reaction (see above). The thermal cycling program for all primer sets was 94°C for 1 min, 58°C for 1 min, and 72°C for 1.5 min for 45 cycles.
The second protocol was designed to maximize our ability to detect low abundance dopamine receptor mRNAs (in this protocol, peptide mRNAs were amplified as above with 2 μl of template). To this end, a two-stage amplification was designed. In the first step, targeted cDNAs were amplified with three-fourths of the single-cell cDNA (15 μl) as a template in a multiplex PCR reaction. All dopamine receptor primers were added to a reaction mixture containing the same concentration of reagents as with conventional PCR, except for slightly elevated MgCl2 (4.0 mm) and dNTPs (1.0 mm) (Chamberlain and Chamberlain, 1994). Twenty cycles were performed with the parameters of 94°C for 1 min, 58°C for 1 min, and 72°C for 3 min. Then an aliquot of this PCR product (2 μl) was used as a template for a second round (30 cycles) of “touchdown” PCR amplification with each pair of specific primers. In this round, after the first ten cycles, the annealing temperature was decreased one degree every other cycle, resulting in a final annealing temperature of 48°C.
The PCR primers were developed from GenBank sequences for dopamine receptors and peptides (ENK and SP) with commercially available software (OLIGO, National Biosciences, Plymouth, MN). Primers were synthesized either by the St. Jude’s Molecular Resource Center (Memphis, TN) or by Life Technologies. The DA receptor primers targeted regions, including the putative third cytoplasmic domain, where the greatest divergence in base sequences was present. For the D1a cDNA, the primers were D1aR1–5′-CAG TCC ATG CCA AGA ATT GCC AGA-3′ and 5′-AAT CGA TGC AGA ATG GCT GGG TCT-3′ (225 bp product) or D1aR2–5′-GAC AAC TGT GAC ACA AGG TTG AGC-3′ and 5′-ATT ACA GTC CTT GGA GAT GGA GCC-3′ (609 bp product; GenBank accession M35077). For the D2 cDNA, the primers were D2R–5′-GCA GTC GAG CTT TCA GAG CC-3′ and 5′-TCT GCG GCT CAT CGT CTT AAG-3′ (404 and 317 bp; GenBank accession M36831). For the D3 cDNA, the primers were D3R1–5′-CAT CCC ATT CGG CAG TTT TCA A-3′ and 5′-TGG GTG TCT CAA GGC AGT GTC T-3′ (201 bp) or D3R2–5′-AGC ATC TGC TCC ATC TCC AAC CC-3′ and 5′-A GGA GTT CCG AGT CCT TTC CAC G-3′ (461 bp product) (GenBank accession M69189). For the D4 cDNA, the primers were D4R1–5′-CTA CTC AGG GTC CCC TCT TC-3′ and 5′-TGA TCT TGG CGC CTC TCT TTC-3′ (189 bp product) or D4R2–5′-TC ATG CTA CTG CTT TAC TGG GCC A-3′ and 5′-T CTG AGA GAG GTC TGA CTC TGG TC-3′ (223 bp product; GenBank accession M84009). For the D1b cDNA, the primers were D1bR1–5′-AGT CGT GGA GCC TAT GAA CCT GAC-3′ and 5′-GCG TCG TTG GAG AGA TTT GAG ACA-3′ (517 bp product; GenBank accession M69118).
The primers for ENK cDNA (GenBank accession M28263; Howells, 1986) were 5′-AAC AGG ATG AGA GCC ACT TGC-3′ (nucleotides 314–334) and 5′-CTT CAT CCG AGG GTA GAG ACT-3′ (nucleotides 770–790). This set yielded a product of 476 bp. The primers for SP (Carter and Krause, 1990) were 5′-TGA GCA TCT TCT TCA GAG AA TCGC-3′, which targets exon 3, and 5′-ATC GCT GGC AAA CTT GTA CAA CTC-3′, which targets exon 7. This set should yield products of 513 or 468 bp, depending on splicing of exon 4. The product for genomic DNA is expected to be 5741 bp, which was never observed.
PCR products (amplicons) were visualized by staining with ethidium bromide and separated by electrophoresis in 1.5–2% agarose gels. Amplicons were sequenced with a dye termination procedure by the University of Tennessee Molecular Resource Center (Memphis, TN) and found to closely match published sequences.
Care was taken to insure that the PCR amplicons arose from cellular mRNA. RT and PCR reactions were performed following procedures designed to minimize the chances of cross contamination (Cimino et al., 1990). Negative controls for contamination from extraneous and genomic DNA were run for every batch of neurons. To insure that genomic DNA did not contribute to the PCR products, neurons were aspirated and processed in the normal manner, except that the reverse transcriptase was omitted. Contamination from extraneous sources was checked by replacing the cellular template with water. Both controls were consistently negative in these experiments.
Tissue RT-PCT. Dorsal neostriatum (see Fig.1A) was dissected from rat brain slices (400 μm) and homogenized in 1 ml of Trizol reagent (Life Technologies) per 10–30 mg of tissue. After 5 min of incubation at room temperature, 0.2 ml of chloroform was added and mixed with the homogenized sample. The tube was incubated at room temperature for 2–3 min and then centrifuged for 15 min at 4°C. The upper aqueous phase that contained RNA was transferred to a fresh tube. Then RNA was precipitated from the aqueous phase by mixing with 0.5 ml of isopropyl alcohol, incubating at room temperature for 10 min, and then centrifuging for 10 min at 4°C. The supernatant was removed, and the RNA pellet was washed with 75% ethanol. The RNA pellet was air-dried and then dissolved in DEPC-treated water. The isolated RNA had anA260/A280ratio of 1.6–1.8, and the RNA yield was ∼1.5 μg/mg of brain tissue.
Before RT-PCR, the isolated RNA was treated with amplification grade DNase I (Life Technologies) to eliminate genomic DNA. The reaction mixture (10 μl) contained 1 μg of RNA, 1 μl of 10× DNase I reaction buffer, 1 μl of DNase I (1 U/μl), and 8 μl of DEPC-treated water. The tube was incubated at room temperature for 15 min. The reaction was terminated by adding 1 μl of 25 mmEDTA and heating for 10 min at 65°C. Conventional PCR analysis of dopamine receptor cDNAs was performed as described above. After 40 cycles, products were separated by electrophoresis in a 1.5% agarose gel stained with ethidium bromide. As a control for genomic contamination, samples were prepared as described above, except that the reverse transcriptase was omitted. These experiments consistently yielded negative results.
RESULTS
All five DA receptor mRNAs are expressed in dorsal neostriatum
As a first step in the study of cellular localization of DA receptor mRNAs, tissue expression was examined. Conventional RT-PCR analysis of dorsal neostriatal mRNA isolated from coronal slices revealed that all five receptor mRNAs were present (Fig.1A,B). Although a quantitative analysis was not attempted, D1a and D2 mRNAs seemed to be the most abundant of the five in these assays, as suggested by other techniques (Sibley, 1995).
The isolation of single cells for patch-clamp analysis requires that slices be held for up to 1 hr before dissociation. During this period of time, it is possible that changes in the transcription of DA receptor genes take place. To test for this possibility, tissue was placed immediately after slicing into an inhibitor of RNA polymerase II, α-amanitin (10 μg/ml; Lewin, 1990), and held for 1 hr. As shown in Figure 1C, the RT-PCR profile of tissue in which RNA polymerase II had been inhibited was indistinguishable from untreated tissue (Fig. 1B; n = 3). As a consequence, it seems unlikely that significant alterations in the transcription of DA receptor mRNAs occur during tissue preparation or in the first hours after slicing.
Over one-half of all acutely isolated neostriatal neurons respond to both D1 and D2 class agonists
Previous electrophysiological studies of DA receptor function have found that substantial fractions of medium spiny neurons respond to both D1 and D2 class agonists (Uchimura et al., 1986; Akaike et al., 1987; Hu et al., 1990; Surmeier et al., 1992;Cepeda et al., 1993), suggesting that D1 and D2class receptors are colocalized, but much of this work is open to the criticism that direct and indirect effects of dopaminergic agonists were not convincingly discriminated, making strong inferences about receptor colocalization difficult.
To attempt to understand the extent to which these earlier estimates of receptor colocalization were a consequence of indirectly mediated effects, we studied acutely dissociated neostriatal neurons. We have shown earlier that sodium currents in acutely isolated, retrogradely labeled striatonigral neurons commonly were modulated by both D1 and D2 class agonists (Surmeier et al., 1992). To follow up these observations, we patch-clamped randomly sampled, medium-sized neurons (<10 pF) and determined the ability of D1 and D2 class agonists to modulate Ba2+ current through Ca2+ channels (Surmeier et al., 1995a,b). As shown previously, D1 and D2class agonists frequently were able to modulate currents in the same cell. A plot of peak current as a function of time and ligand application in one such cell is shown in Figure2A. In this cell, both quinpirole (Fig. 2B) and 6-chloro-PB (Fig. 2C) reversibly reduced peak currents. Of our sample (n = 26), ∼50% responded to both agonists, whereas ∼20% responded only to D1 class agonists and ∼20% responded only to D2 class agonists (Fig. 2D).
Single-cell RT-PCR analysis of DA receptor mRNAs
To determine the molecular identity of the receptors mediating the responses to the dopaminergic ligands, single-cell RT-PCR techniques were used. Neurons initially were divided into three groups on the basis of ENK and SP mRNA expression. The expression of these releasable peptides is correlated strongly with the axonal projection pattern of medium spiny neurons (Graybiel, 1990; Gerfen, 1992). Neurons expressing ENK mRNA alone project primarily to the globus pallidus, whereas neurons expressing either SP alone or both ENK and SP project to both the globus pallidus and the substantia nigra (Kawaguchi et al., 1990; our unpublished observations). As described below, the DA receptor mRNAs expressed by these three groups differed markedly.
An obvious problem in attempting to determine how the expression of several mRNAs is coordinated in the same cell is that low abundance mRNAs may be missed inadvertently. For example, to assay for five different mRNAs, we must divide the total compliment of reverse-transcribed cDNA into at least five separate aliquots using conventional PCR techniques. If there were only one or two times the number of templates needed to reach the detection threshold for the PCR reaction, then there would be a significant chance of missing that particular cDNA. To generate a model of the probability of such an occurrence, we used Poisson statistics (Feller, 1950). Shown in Figure3 are plots of blank probability as a function of cDNA abundance (expressed in PCR threshold units). As the fraction of the total cellular cDNA used in the detection reaction increased (from one-twentieth to one-half), the probability of detection increased (as expected). This simulation suggests that the relationship between detection probability and the fraction of total cellular cDNA used in the detection reaction can be used to draw inferences about mRNA abundance at an ordinal scale of measurement.
Enk+/SP− neurons express primarily D2 class DA receptor mRNAs
Medium spiny neurons having detectable levels of enkephalin, but not substance P mRNA (ENK+/SP−), expressed high levels of D2 receptor mRNA. With a one-tenth aliquot, D2 cDNA was detected in nearly 90% of this group. A photomicrograph of an agarose gel in which the PCR products from an individual cell have been separated by electrophoresis is shown in Figure 4A. The other DA receptor mRNAs were rarely detected with this small aliquot. The pattern of DA receptor mRNA expression in ENK+/SP− neurons studied in this way is summarized in Figure 4B. In this panel (and in subsequent summary diagrams), the extent to which particular mRNAs were found together is coded by the extent to which their lanes are shaded at similar points along the abscissa. For example, two mRNAs that were always found together would have lanes that were shaded at the same points along the horizontal axis; mRNAs that did not colocalize would have lane shading that did not overlap.
As the fraction of the total cellular cDNA used in the PCR reactions increased, other DA receptor mRNAs emerged. Doubling the template concentration for D1a, D1b, D3, and D4 reactions led to a modest increase in the detection frequency of D1b mRNA (Fig. 4C,D). The short isoform of the D2 mRNA was seen consistently in this condition (Fig. 4C,D). Using a two-stage multiplex procedure to maximize the detectability of low abundance mRNAs increased the detection rate of D3 and D4 mRNAs, primarily. An example of an ENK+/SP− cell in which D2, D3, and D4 mRNAs were seen is shown in Figure 4E. Using this approach, we found D3 and/or D4 mRNAs in >50% of ENK+/SP− neurons. The increasing detection frequency of these mRNAs with aliquot size suggests that they were present in low abundance, as compared with D2 mRNAs. The variation in the detection frequency of D1a and D1b mRNAs suggests that they are present in moderate abundance in a relatively small percentage (10–25%) of these neurons.
ENK−/SP+ neurons express D1amRNA in addition to D3 or D4 receptor mRNAs
Medium spiny neurons having detectable levels of substance P, but not enkephalin mRNA (ENK−/SP+), expressed high levels of D1a receptor mRNA. With one-tenth aliquots, D1a mRNA was detected in 80% of the sample (Fig.5A,B). D3 mRNA was detected in one-half of this sample, but the other receptor mRNAs either were not seen or were seen rarely. Doubling the aliquot size increased the D1a detection rate to 100% and did not change the detection rate for the D3 mRNA (Fig. 5C,D).
With the multiplex protocol, the detection frequencies for D1a and D3 mRNAs were similar to those seen with smaller aliquots (Fig. 5F). However, the detection frequency of D4 mRNA increased to 25% (as seen with one-fifth aliquots; Fig. 5E,F). In addition, D2 mRNA, which had not been seen with one-tenth or one-fifth aliquots, was seen in ∼20% of the sample, suggesting that D2 mRNAs were of relatively low abundance in this group.
Enk+/SP+ neurons coexpress D1aand D2 mRNAs
Neurons coexpressing detectable levels of substance P and enkephalin (ENK+/SP+) were the only group in which a large percentage of neurons coexpressed D1a and D2 mRNAs. With one-tenth aliquots, D2 mRNA was detected most commonly in this sample (78%; Fig.6A,B). D1a receptor amplicons were seen less frequently, and both amplicons were detected in approximately one-fifth of the sample. The other receptor mRNAs either were not detected or were seen rarely with the small aliquot. Doubling the aliquot did not change the D1a and D2 detection frequencies significantly (Fig.6D).
However, with the multiplex protocol, the detection rate for D1a mRNA increased to 90% (Fig. 6F). As a consequence, D1a and D2 mRNAs were detected in 70% of this group (Fig. 6E,F). The detection frequency of D3, D4, and D1b mRNAs also increased with the size of the aliquot used, suggesting that they were expressed at relatively low levels in 20–30% of ENK+/SP+ neurons.
Low abundance DA receptor mRNAs give rise to functional protein
Although the RT-PCR technique was capable of detecting the presence of mRNAs for D3, D4, and D1b receptors, there is no guarantee that functionally significant levels of receptor protein were generated from these transcripts. In an attempt to determine whether these transcripts could give rise to functional protein, whole-cell recording was combined with conventional single-cell RT-PCR. Several neurons were found to express detectable levels of D2 and D1b mRNA, but not D1a mRNA (Fig. 7A, see above). In two of three neurons with this expression profile, the application of the D1 class agonist 6-chloro-PB brought about a robust modulation of Ca2+ currents, suggesting that D1b receptors were responsible (Fig. 7B). As shown above, in other cells D3 or D4 mRNA was detected, but not D2 mRNA (Fig. 7C). In two of three neurons with this expression profile, the application of the D2 class agonist quinpirole reversibly decreased Ca2+ currents, suggesting that D4 receptors were responsible for the modulation (Fig. 7D). It is possible in both of these experimental series that undetected levels of D1a or D2 receptor mRNA were responsible for the protein mediating the response, but this seems unlikely. Providing evidence that D3 receptors possessed a similar coupling proved difficult, in that all of the cells studied that expressed D3 receptors also had detectable levels of another D2 class receptor mRNA (n = 5). Although these results do not prove that D3, D4, and D1b receptor mRNAs invariably give rise to functional protein, they do suggest that these mRNAs are capable of being translated and properly processed.
D2 receptor splice variants are colocalized in the majority of ENK+ neurons
In the whole-striatum RT-PCR analysis and in the single-cell analyses, the long-splice variant of the D2 receptor mRNA (D2l) was the most abundant. This agrees with previous estimates of D2 isoform expression in the neostriatum (Sibley, 1995). Although less abundant, the short-splice variant (D2s) was seen in a significant number of neurons, usually in conjunction with D2l mRNA. For example, in the ENK+/SP− group, over one-half of the neurons coexpressed detectable levels of short and long isoforms (as seen with the multiplex protocol; n = 36; Fig.8A). Neurons expressing detectable levels of only the short form were a small fraction of the total (8%). In the ENK+/SP+ group, there was a similar pattern of expression (approximately one-half of the neurons coexpressed D2s and D2l) with a slightly higher percentage of neurons (20%, n = 15) expressing only the short isoform (Fig. 8B). Although the number of ENK−/SP+ neurons expressing D2isoforms was small (n = 3/54), the short and long isoforms were never found together in this group (Fig. 8C).
DISCUSSION
The principal goal of this study was to understand why anatomical and functional studies of DA receptor localization in the neostriatum have been so discordant (Gerfen, 1992; Surmeier et al., 1993). The issue is obviously an important one, having a direct bearing on treatment strategies for a variety of psychomotor disorders of the basal ganglia, most prominently Parkinson’s disease.
It was our working hypothesis that at the heart of the dispute was the assumption that only D1a and D2receptors were functionally significant in the workings of the dorsal neostriatum. This inference allowed the equation of responsiveness to D1 class agonists with the presence of D1areceptors and responsiveness to D2 class agonists with the presence of D2 receptors. In agreement with a number of other studies (Ariano and Sibley, 1994; Bergson et al., 1995; Booze and Wallace, 1995; Choi et al., 1995; Larson and Ariano, 1995; Schoots et al., 1995; Vrana et al., 1995), our results argue that, in fact, this assumption is incorrect. Significant levels of D3, D4, and D1b receptor mRNAs are present in medium spiny neurons, and their expression is correlated with that of the releasable peptides SP and ENK.
Strengths and weaknesses of single-cell RT-PCR
The view that only D1a and D2 receptors are important in the dorsal neostriatum derives primarily from in situ hybridization studies (Bouthenet et al., 1991; Gerfen, 1992;Meador-Woodruff et al., 1992; Diaz et al., 1995). Although this approach has proven to be an important tool in localizing mRNAs, it cannot be used definitively to exclude the presence of particular mRNAs. Several factors can compromise the ability of this technique to accurately localize mRNAs, such as secondary structure and protein binding (Chesselet et al., 1995). In addition, mRNA must be localized in relatively high abundance within the soma to be detected. A number of studies have shown that mRNAs coding for proteins found in dendritic synapses (like DA receptors) are moved to these dendritic sites for translation (Steward and Banker, 1992; Miyashiro et al., 1994).
RT-PCR techniques are not encumbered by the same limitations. The detection threshold for these techniques is in principle only a few mRNA copies (Mullis et al., 1994). However, conventional RT-PCR techniques with mRNA pooled from all the cell types in a brain region obviously lack the ability to resolve heterogeneity within a population of cells. This shortcoming is of obvious importance not only in distinguishing between mRNAs originating from neuronal and non-neuronal sources but also in distinguishing between mRNAs originating from different neuronal phenotypes within a region. Single-cell RT-PCR (Monyer and Lambolez, 1995; Yan and Surmeier, 1996) and RT-aRNA (Van Gelder et al., 1990; Eberwine et al., 1992; Surmeier et al., 1992) techniques overcome this problem. However, these approaches are not without shortcomings. RNAs that are exclusively or primarily dendritic may be missed. This is potentially a problem in our experiments and may have led to an underestimate of colocalization.
However, the converse argument can also be made, at least for the somatic compartment. Although in situ hybridization techniques may miss particular mRNAs, the exquisite sensitivity of RT-PCR techniques makes them capable of detecting mRNAs with very low copy numbers. At the extreme, genomic DNA may be amplified. This was a major concern in our initial studies. Although there were obvious differences in cellular PCR profiles (arguing against genomic amplification), we tried explicitly to design primer sets that spanned a splice site when possible (D2 class receptor genes all contain introns). If genomic DNA were amplified, then a high-molecular-weight, unspliced product would be present in the amplified material. We have never seen any evidence for this in our experiments, including those using 80 rounds of amplification. Johansen et al. (1995) have also found that genomic material is not amplified in single-cell preparations.
Given that the RT-PCR procedure amplifies transcripts generated from mRNA, can inferences be drawn about relative mRNA abundance? A number of studies have attempted to generate abundance estimates at the interval or ratio scale of measurement with PCR (Ramakrishnan et al., 1994). Although the PCR amplification stage is capable of accurate quantitation, the reverse transcription (RT) step is not, at present. Although there is no way to escape the limitations this uncertainty poses completely, restricting the level of quantitation to an ordinal scale (rank order) mitigates the consequences of deficiencies in the RT step. The simplest, least-assumption-bound way of generating ordinal estimates of mRNA abundance was to look at correlations between the probability of detection and the fraction of the total cellular cDNA used as template. This approach showed that, when the cellular analysis was restricted to only high abundance mRNAs, a pattern very similar to that found with in situ hybridization was seen. Namely, D1a mRNA was detected in SP-expressing neurons and D2 mRNA was detected in ENK-expressing neurons, whereas D3 mRNA was seen in a small population of cells (∼20%) and D4 and D1b mRNAs were virtually undetected. Increasing the sensitivity of our mRNA assay revealed more colocalization of D1a and D2 mRNAs but also more widespread expression of D3, D4, and D1b mRNAs. The frequency with which medium spiny neurons responded to both D1 and D2 class agonists can be explained only if these lower abundance receptors give rise to functionally significant receptors.
The expression patterns of D3, D4, and D1b mRNAs
Previous work had shown that D3 receptor mRNA and protein are present in the dorsal neostriatum, albeit at lower levels than in ventral neostriatum (Ariano and Sibley, 1994; Booze and Wallace, 1995; Larson and Ariano, 1995; Vrana et al., 1995). Our work extends these studies in showing that D3 mRNA is clearly present and most abundant in neurons expressing SP but not ENK. Approximately one-half of these neurons had detectable levels of D3 mRNA. This proportion did not change significantly with increasing amounts of cellular cDNA, suggesting that it was relatively abundant in this subset of cells and not expressed or expressed at very low levels in other SP+/ENK− neurons. In contrast, the D3 mRNA detection probability steadily increased with the amount of cellular cDNA used as starting material in the ENK+/SP− population, suggesting that D3 mRNA was expressed at lower levels in these neurons. Even lower levels of D3 mRNA were present in the group of neurons coexpressing SP and ENK.
Although the presence of D4 mRNA and receptors in dorsal neostriatal neurons has been controversial (O’Malley et al., 1992;Murray et al., 1995; Schoots et al., 1995; Vrana et al., 1995), our results suggest that they are present at low-to-moderate levels in a subset of medium spiny neurons. To ensure that the D4transcript was properly identified, we used two different primer sets [one that spanned a splice site (O’Malley et al., 1992)], and the PCR product was sequenced (as for the other amplicons). D4mRNA was found to some extent in all three groups of neurons but seemed to be the most abundant in neurons expressing SP, but not ENK. Our results also suggest that these transcripts give rise to functional receptor protein. In many of these neurons, D3 mRNA was also present. Because all of the neurons in this group expressed D1a mRNA, the presence of D3 and D4mRNAs argues that D1 and D2 class receptors should be colocalized in nearly 60% of ENK−/SP+ neurons (Fig.9A).
In contrast, presumptive striatopallidal neurons (ENK+/SP−) did not commonly have detectable levels of both D1 and D2 class mRNAs (Fig.9B). All of the neurons in this group had D2mRNA, frequently coexpressing short and long isoforms, but D1a and D1b mRNAs seemed to be present in only a small subset of neurons (∼10-25%). Previous in situhybridization work had suggested that D1b mRNA was not present in the rat dorsal striatum (Meador-Woodruff et al., 1992; cf.Laurier et al., 1994). Given its restriction to a small subset of neurons, this is not surprising. However, within this group, activation of the receptor does have functionally significant consequences, as our results show. Furthermore, in another numerically small but functionally significant group of neurons—the cholinergic interneuron—D1b mRNA is found in nearly every cell tested (n = 34; Yan, Song and Surmeier, unpublished observations). The pattern of D1b expression seen here and its presence in cholinergic neurons agree well with the distribution of mRNA (Choi et al., 1995) and protein (Bergson et al., 1995) found in the primate brain.
The only group of neurons that consistently coexpressed D1aand D2 mRNAs also coexpressed SP and ENK mRNAs. Nearly all of these neurons (80%) had both D1 and D2class mRNAs (Fig. 9C). The restriction of high levels of D1a and D2 mRNA to this group of cells is in agreement with previous in situ hybridization work showing a strong correlation between the expression of ENK and D2mRNA on the one hand and SP and D1a mRNA on the other (Gerfen, 1992; LeMoine and Bloch, 1995).
An unresolved issue is how large a population coexpresses these peptides (and receptors). In a sample of acutely isolated medium-sized neurons profiled just for peptide mRNA (n = 64), 17% had detectable levels of SP and ENK, whereas 38% were ENK+/SP−, 42% ENK−/SP+, and 3% were ENK−/SP−. Both ENK−/SP+ and ENK+/SP+neurons (but not ENK+/SP− neurons) are retrogradely labeled by rhodamine bead injections into the substantia nigra (data not shown). Previous peptide colocalization estimates from immunocytochemical studies have varied widely—from >30 to 1–2% (Penny et al., 1986; Anderson and Reiner, 1990; Besson et al., 1990). The variance in the estimates of D1a and D2mRNA colocalization is surprisingly similar (cf. Gerfen, 1992; Lester et al., 1993; LeMoine and Bloch, 1995). The lack of agreement in these studies undoubtedly turns on the substantial technical obstacles facing attempts to judge protein or mRNA coexpression. These obstacles are considerably less significant in our situation. If our estimate of the relative sizes of the three efferent populations is taken at face value (Fig. 9), then it suggests that nearly 70% of the medium spiny neurons in the “direct” efferent pathway (projecting to substantia nigra/entopeduncular nucleus) colocalizes D1 and D2 class receptors—in stark contrast to neurons presumably projecting exclusively to the globus pallidus (ENK+/SP−). A schematic summary of our results is shown in Figure 10, in which we have assumed the well established correlation between peptide expression and projection site. An important question is whether these groups can be subdivided further on the basis of membership in other functionally significant classes or compartments (most notably the striosomal and matrix compartments).
Functional implications
The demonstration that D1 and D2class receptors are coexpressed by a significant percentage of medium spiny neurons, particularly those within the striatonigral pathway, establishes a cellular and molecular foundation for the well known interactions between these signaling pathways that does not require the postulation of poorly defined circuit level interactions. An understanding of the functional consequences of the differences in receptor expression hinges on a more complete characterization of the subcellular distribution of the lower abundance receptors (particularly D3 and D4 receptors) and how their intracellular linkages differ from D1a and D2receptors. It is evident, however, that the striatopallidal and striatonigral systems cannot be viewed as functionally homogeneous. This is particularly true of the striatonigral system in that there are at least three neuronal subsets (all of approximately the same size) that can be distinguished on the basis of their expression of D2 class receptors. It is tempting to speculate that these basal ganglia subsystems may represent components of parallel processing modules identified by other criteria (Alexander et al., 1986; Graybiel et al., 1994).
Footnotes
This work was supported by U.S. Public Health Service Grants NS28889 and NS34696 to D.J.S., NS26487 to D.J.S. and S. T. Kitai, and a Parkinson’s Disease Foundation fellowship to W.J.S. We thank Drs. Jim Eberwine, Marie-Francoise Chesselet, and Kent and Sheila Vrana for many helpful discussions, Drs. Kent and Sheila Vrana for help in the initial PCR experiments, and Dr. Lorina Dudkin for her excellent technical assistance.
Correspondence should be addressed to Dr. D. James Surmeier, Department of Anatomy and Neurobiology, University of Tennessee, Memphis, TN 38163.