Abstract
A conserved feature of sound processing across species is the presence of multiple auditory cortical fields with topographically organized responses to sound frequency. Current organizational schemes propose that the ventral division of the medial geniculate body (MGBv) is a single functionally homogenous structure that provides the primary source of input to all neighboring frequency-organized cortical fields. These schemes fail to account for the contribution of MGBv to functional diversity between frequency-organized cortical fields. Here, we report response property differences for two auditory fields in the rat, and find they have nonoverlapping sources of thalamic input from the MGBv that are distinguished by the gene expression for type 1 vesicular glutamate transporter. These data challenge widely accepted organizational schemes and demonstrate a genetic plurality in the ascending glutamatergic pathways to frequency-organized auditory cortex.
Introduction
The ability to discriminate sound sources in natural scenes relies on distinct sound processing properties of different frequency-organized auditory cortices in mammals (Polley et al., 2006; Lomber and Malhotra, 2008; Nodal et al., 2010), but the mechanisms underlying these functional differences remain a mystery. By definition, “core” auditory cortices have a fine-grain topographic organization of responses to sound frequency and receive their primary thalamic input from the ventral division of the medial geniculate body (MGBv) of the thalamus (Kaas and Hackett, 2000). Several mammals including humans have multiple “core” fields outside the primary (A1) auditory cortical field (Formisano et al., 2003; Schreiner and Winer, 2007; Hackett, 2011). Gerbils, guinea pigs, and rats have overlapping hearing ranges and similar frequency organization of A1 relative to the anterior auditory field (AAF) and a third ventral field with distinct frequency organization (Wallace et al., 2000; Kalatsky et al., 2005; Donishi et al., 2006; Budinger and Scheich, 2009). In the rat, this ventral field is called the suprarhinal auditory field (SRAF) (Polley et al., 2007). Neurons in caudal SRAF (cSRAF) respond more selectively to sound frequency and position cues than their frequency-matched counterparts in A1 (Polley et al., 2007; Higgins et al., 2010). These sound-encoding differences could stem from physiologic and anatomic differences evident along the caudal-rostral axis of MGBv (Redies et al., 1989; Rodrigues-Dagaeff et al., 1989; Cant and Benson, 2007; Storace et al., 2011). However, at present there are no known anatomic or molecular indices that subdivide MGBv along its caudal-to-rostral axis (Morest, 1964; Clerici and Coleman, 1990; Winer et al., 1999a, 2005; Storace et al., 2010).
The majority of the principal neurons in MGBv likely release glutamate neurotransmitter as they express the gene and corresponding mRNA encoding type 2 vesicular glutamate transporter (VGLUT2) which in turn is critical for loading glutamate into synaptic vesicles (Fremeau et al., 2004a; Barroso-Chinea et al., 2007; Hackett et al., 2011; Ito et al., 2011). This expression pattern is consistent with a VGLUT2 protein-dense terminal field found within thalamic recipient layers of sensory cortex (Kaneko et al., 2002; Nahmani and Erisir, 2005; Graziano et al., 2008; Hackett and de la Mothe, 2009). However, recent studies find that type 1 glutamate transporter (VGLUT1) mRNA is coexpressed with VGLUT2 mRNA in a subset of MGBv neurons in rodents (Barroso-Chinea et al., 2007; Ito et al., 2011). This raises the question of whether MGBv pathways to different frequency-organized cortical fields have different gene expression patterns for glutamate transporter type.
Here, we identify sound response differences between adjacent frequency-organized cortices A1 and cSRAF in the rat. We confirm that the MGBv provides the majority of inputs to a targeted injection layer in both fields. Furthermore, we find that MGBv neurons that project to A1 express VGLUT1 mRNA, whereas those that project to cSRAF do not. These findings offer a novel explanation and molecular marker for functional diversity of ascending sensory pathways to frequency-organized auditory cortices.
Materials and Methods
Surgery and imaging in adult rats.
Animals were housed and handled according to a protocol approved by the Institutional Animal Care and Use Committee (University of Connecticut, Storrs). Intrinsic imaging was performed on 28 male rats with a mean age and weight of 71 ± 4 d and 214 ± 14 g, respectively. Each of these rats was used in one or more additional experiments that include retrograde tracer injections, in situ hybridization, unit recordings, and immunohistochemical analysis. Anesthesia was induced with a mixture of ketamine (50 mg/kg), xylazine (5 mg/kg), and acepromazine (2.5 mg/kg), and was maintained in a stage III, plane II of general anesthesia with sodium pentobarbital (25 mg/kg, i.p.) for surgical and experimental procedures. Anesthesia was supplemented as needed to maintain areflexia, and anesthetic depth was monitored periodically via pedal reflex and heart rate (electrocardiogram, GE Healthcare Life Sciences). Tracheotomy and cisterna magnum puncture were performed and atropine sulfate (0.1 mg/kg) and dexamethasone (0.25 mg/kg) were administered to minimize respiratory noise and cerebral edema. The skull and dura were removed to expose the temporal cortex. The exposure was covered with agar and sealed with a glass coverslip for optical imaging of intrinsic metabolic responses.
Sound stimuli, delivery, and design.
All sounds were presented to both ears through hollow ear-bars via custom modified (Beyer DT 770) speakers calibrated for frequencies in the range of 1–45 kHz (±3 dB) and linearity (input/output coherence > 0.95) in a closed system with a 400-tap FIR inverse filter implemented on a Tucker Davis Technologies (TDT) RX6 multifunction processor. Auditory stimuli used for multiunit recordings were flat spectrum noise tokens (<5 ms rise time, 50 ms duration) delivered via a TDT RX6 multifunction processor at 96 kHz (Higgins et al., 2010). Auditory stimuli used for Fourier intrinsic imaging were a continuous sequence of 16 pure tones (50 ms duration, 5 ms rise and decay time, 250 ms interstimulus interval) from 2 to 32 kHz (0.25 octave steps) with a 4 s period in ascending or descending order (Kalatsky et al., 2005; Polley et al., 2007; Higgins et al., 2010; Storace et al., 2010, 2011) delivered via a Lynx Studio Technology professional audio card at 98 kHz.
Fourier optical imaging.
Fourier optical imaging was used to map tone frequency response organization in auditory cortex. Images were acquired with a Dalsa 1M30 CCD camera with a 512 × 512 pixel array covering a 4.6 × 4.6 mm2 area (Kalatsky et al., 2005; Polley et al., 2007; Higgins et al., 2010; Storace et al., 2010, 2011). Surface vascular patterns were imaged with a green (546 nm) interference filter at a 0 μm plane of focus. Intrinsic activity was imaged at a plane of focus of 600 μm below the surface blood vessels using a red (610 nm) interference filter. Ascending and descending tone sequences yielded similar frequency topographies in all regions; thus, hemodynamic delay was corrected by subtracting the ascending and descending frequency phase maps to generate a difference phase map (Kalatsky et al., 2005). The optical best frequency (BF) was the value obtained from single pixels at designated positions in the optical difference phase map. The locations of A1, cSRAF, rostral SRAF (rSRAF), and ventral (VAF), posterior (PAF), and anterior auditory field (AAF), were determined physiologically by a change in the direction of the frequency gradient observed in the optical frequency map as described in detail previously (Higgins et al., 2010). Post hoc histochemistry confirmed that the dorsal border of A1 and the ventral border of cSRAF fall within the cytoarchitectonically defined temporal area 1, Te1 (LeDoux et al., 1985; Paxinos and Watson, 1997).
Temporal multiunit response analysis.
Multiunit responses were recorded in A1 and cSRAF (NA1 = 159; NcSRAF = 91) in response to white noise bursts presented at varying interaural level differences (ILDs) and average binaural levels (ABLs). The dataset range was limited in the ABL dimension from threshold to 40 dB above threshold, and best ILD ± 5 dB. Analyses for the calculation of best ILD and ABL thresholds have been previously detailed (Higgins et al., 2010). The poststimulus time histogram (PSTH) was generated from the limited response range and used to find the maximum amplitude, time to onset and time to peak. Rising and falling slopes were identified by fitting the PSTH with an α function of the form: where A is the magnitude, t is time, t0 is the time of significant onset, Trise is the time from onset to peak, and Tdecay is the time from peak to offset. Times to onset and offset were determined as the first and second crossings of a set significance threshold (α = 0.05) (Escabí et al., 2007), and were used to determine the significant portion of the PSTH to model. The best fit function for each unit was calculated using the lsqcurvefit MATLAB function (MathWorks), and cross validation was used to assess quality of fit. Each stimulus was presented 16 times, therefore the model was based on eight trials, then validated by comparison to the other eight trials. This procedure was repeated 500 times for each unit with a random selection of trials designated in each iteration. A subset of the data was analyzed over 1000 iterations of the cross validation procedure and no significant differences were observed.
Retrograde tracer and VGLUT2 histology.
Following optical imaging, the retrograde tracers Cholera toxin B (CTB, 10 μg/μl List Biological Laboratories) and CTB conjugated with gold (CTBG, 10 μg/μl List Biological Laboratories) were injected (Nanoliter 2000, World Precision Instruments) using a glass capillary (tip diameter 10–12 μm) into frequency-matched (NA1 = 8, BFA1 = 8.4 ± 0.6 kHz, NcSRAF = 8, BFSRAF = 7.3 ± 0.4 kHz, t(16) = 1.56, p = 0.14) contours in A1 and cSRAF. Injections targeted layer IV in both areas and had a cortical surface separation of 1.8 ± 0.05 mm. Tracer transport occurred during an ∼24 ± 0.5 h survival time under anesthesia. For experiments in which VGLUT2 protein expression was analyzed, a similar procedure was used except that green and red fluorescent beads (undiluted, Lumafluor) were injected into A1 and cSRAF (N = 4) for identification purposes. Rats were given a lethal dose of pentobarbital and perfused transcardially with 4% paraformaldehyde. Brains were removed, blocked (Model PA 001 Rat Brain Blocker, Kopf Instruments), and postfixed in 4% paraformaldehyde. Following postfixation, the brains were cut at 50 μm on a vibratome (N = 2; retrograde tracer experiments) (VT1000, Leica Microsystems), or cryoprotected in diethylpyrocarbonate-treated 30% sucrose in phosphate buffer for at least 3 d, and cut at 40 μm (N = 7; retrograde tracer and in situ experiments) or 20 μm (N = 4; VGLUT2 immunohistochemistry experiments) on a freezing microtome (HM440E, Microm International). CTB label (CTB-IGG, List Biological Laboratories) and VGLUT2 protein (Synaptic Systems, catalog #135 402) were detected using immunohistochemistry, and CTBG was detected via silver enhancement (KPL Labs). All rats that received dual CTB and CTBG tracer injections had at least one series of tissue double-reacted for both tracers. In cortex, locations of Te1 and temporal areas 2 and 3 (Te2 and Te3, respectively) were determined from cytoarchitectonics of Nissl-reacted and immunoreacted sections (LeDoux et al., 1985; Paxinos and Watson, 1997). Photomicrographs and drawings of histological sections were prepared using Adobe Photoshop and Illustrator (Adobe Systems) and were contrast enhanced, brightened, and color balanced.
Retrograde tracer analysis.
Retrograde tracer injection deposit volumes and locations did not overlap and were measured in Stereo Investigator (MBF Bioscience) using the Cavalieri estimator. Injection deposits had an average caudal-rostral spread of 600 ± 52 μm (N = 14). Injection deposit volumes for A1 and cSRAF were 0.13 ± 0.03 (N = 7) mm3 and 0.2 ± 0.02 (N = 7) mm3, respectively, and were not significantly different (t(12) = 0.27, p = 0.79). The distance from the pia to the center of each injection site was 512 ± 27 (N = 13) μm. The injection deposits fell within the boundaries of A1 and cSRAF, and were located 3.4 ± 0.09 (N = 6) and 1.7 ± 0.06 (N = 7) mm dorsal to the nadir of the rhinal fissure, respectively. Two experiments were excluded from injection volume and distance analysis as their cortex was removed and cut differently for a separate experiment. All histological distance and area measurements were taken directly from the sections.
Bright-field in situ hybridization histochemistry.
Digoxigenin (DIG)-labeled sense and antisense riboprobes were made from the cDNAs of mouse VGLUT1 (nucleotides of 152–1085, GenBank accession number NM_182993.2). The specificity of this riboprobe for VGLUT1 has been previously established in the mouse (Nakamura et al., 2007) and in both the rat and mouse auditory pathway (Ito et al., 2011). Briefly, this sequence for mouse VGLUT1 is 96% homologous to the corresponding rat cDNA, while very low for other members of the VGLUT family. The mouse VGLUT1 riboprobe shares only 75% and 71% homology with the mouse VGLUT2 and VGLUT3 mRNA, respectively. As a result, expression of this VGLUT1 probe is totally absent in the rat inferior colliculus and other structures that lack VGLUT1 expression in mouse (Ito et al., 2011).
In a subset (N = 6) of the brains that received dual retrograde tracer injections, neurons expressing VGLUT1 mRNA were visualized via in situ hybridization in a separate series of tissue as described previously (Ito et al., 2011). Sections in this series were enhanced with silver to visualize neurons retrogradely labeled with CTBG that was injected into either A1 (N = 4), or cSRAF (N = 2) using light microscopy (Basbaum, 1989). One experiment (A1 injection) was excluded from colocalization analysis due to poorly visualized retrograde label, and another (cSRAF injection) was excluded from caudal-to-rostral analysis due to having been cut at a different angle, however it was otherwise included in our colocalization analysis.
Quantifying labeled neurons and estimating MGB division borders.
All labeled neurons (CTB, CTBG, CTB/CTBG, VGLUT1 and VGLUT1/CTBG) positions and total counts were measured using Stereo Investigator (MBF Bioscience) for a series of sections through the MGB. Soma position corresponded to the nucleolus center as measured with a 100× objective lens. Labeled neurons were counted if they contained a single distinct nucleolus, and colocalization was identified as neurons containing the black granular reaction product of silver enhanced CTBG confined within the boundaries of immunohistochemical (DAB, brown) or in situ hybridization (NBT/BCIP, blue) reaction product. Data were collected from noncontiguous sections that were separated by at least 160 or 200 μm, although alternate columns were substituted in the case of damaged or missing sections. For each case, the section corresponding to the bregma position of −6.30 mm was identified and coordinates for all other sections are relative to that reference in 160 μm increments. Sections separated by 200 μm were placed into the nearest bin. Stereo Investigator data were imported into MATLAB to generate cell position plots.
Borders of MGBv, dorsal (MGBd), medial (MGBm), and suprageniculate (SG) nuclei of MGB were based on Nissl-, immuno-, and in situ hybridization-reacted sections, and previous studies in the MGB (LeDoux et al., 1985; Clerici and Coleman, 1990; Paxinos and Watson, 1997; Winer et al., 1999a,b). Quantification of cell counts in each MGB division was performed in MATLAB and measured using cell position plots aligned and overlaid with manually drawn cytoarchitectonic divisions. An automatic routine grouped individual cells into the appropriate nuclear division based on their location within cytoarchitectonic MGB division borders. CTB-, CTBG-, CTB/CTBG-, VGLUT1-, and VGLUT1/CTBG-labeled neurons were analyzed separately from each other.
To control for spurious single-labeled neuron “outliers” in the spatial distributions of labeled neurons, additional bootstrapping statistics were computed for each section using a modified version of our previous technique (Storace et al., 2010, 2011). One thousand random cell distributions were generated from each MATLAB cell position plot marked within a 1500 × 1500 μm2 area of the MGB. Each sample was divided into 60 × 60 μm2 bins to generate two-dimensional density histograms. The mean bootstrapped histogram was computed and bins were selected for further analysis if they qualified as significant by being three or more SDs above the overall mean count for the section. The areal boundaries of significant bins were marked with a contour, and drawn onto the original cell position plots. Additional cell count measurements were made based on whether labeled neurons were located within these contours. These contours selected for 69.5% of the original plotted neurons within the boundaries of the MGBv, MGBd, MGBm, and SG subdivisions of the MGB. Original cell counts for neurons projecting to A1 and cSRAF were 3485 and 3781, respectively, and bootstrapped counts were 2533 and 2523, respectively. The areal pattern of significant clusters was also analyzed using the same coefficient of variation analysis and significance criterion described previously but not included in detail here, and provided additional support for our cytoarchitectonic boundaries in this and prior studies (Polley et al., 2007; Storace et al., 2010, 2011).
Results
Glutamate transporter types covary with response magnitude and timing (Fremeau et al., 2004b); hence, we examined the possibility that A1 and cSRAF have distinctions in these two response properties (Fig. 1). Frequency organization and regional boundaries of A1 and cSRAF were determined by mapping the BF responses to tones using Fourier optical imaging of intrinsic metabolic signals (Materials and Methods; Fig. 1A). Subsequently, multiunit responses to noise were recorded in both cortical regions, and the time course of the PSTH was fit with an α function as shown for PSTH responses from a site in A1 (black) and a site in cSRAF (gray, Fig. 1B). Response magnitudes and slope were greater in A1 than cSRAF, and the response delays were shorter (Fig. 1C–E). Furthermore, the time or delay to peak response increased with dorsal to ventral shifts in cortical positions from A1 to cSRAF (Fig. 1F). These results indicate that A1 and cSRAF have distinct frequency organization and sound response properties.
The pattern and relative density of VGLUT2 protein expression has not been demonstrated for physiologically identified A1 or cSRAF in the rat. Here we used optical imaging to map the frequency response organization and to guide injections of fluorescent beads into frequency-matched regions of A1 and cSRAF (Fig. 2A). The brains were subsequently processed to examine the laminar organization of cortical cell bodies and VGLUT2 protein-expressing terminal fields in A1 and cSRAF. Both cortical fields showed a similar Nissl profile (Fig. 2B,E) and a high density of VGLUT2 protein expression in the middle cell layers (Fig. 2C,F). In contrast, VGLUT2 protein expression dropped markedly in the transition between rSRAF and the ventral “auditory association” area, Te3 (Fig. 2C, Te3 vs rSRAF). A plexus of VGLUT2 putative terminals with bouton-like swellings was evident at high magnification in A1 and cSRAF (Fig. 2D,G). This resembles that seen in other auditory and sensory cortices (Nahmani and Erisir, 2005; Graziano et al., 2008; Hackett and de la Mothe, 2009) and suggests that VGLUT2 protein expression alone does not readily distinguish the thalamocortical pathways to frequency-organized fields, A1 versus cSRAF.
We examined thalamic input sources to A1 and cSRAF as a potential source of response property differences in the two fields (Fig. 3). Retrograde tracers CTB and CTBG were deposited in layer IV of frequency-matched 8 kHz contours in A1 and cSRAF (Fig. 3A). Nissl-stained sections illustrate the distribution of all MGB neurons in caudal (Fig. 3B) and rostral (Fig. 3E) MGB from one dual-injected brain hemisphere. Adjacent sections processed to visualize both tracers illustrate independent and complementary caudal and rostral locations of the subsets of thalamic neurons projecting to cSRAF (Fig. 3C) and A1 (Fig. 3F), respectively. This caudal-to-rostral organization is demonstrated in a series of photomicrographs of sections through the MGB in a separate experiment (Fig. 4A–H, Nissl; I–P, immunoreacted). The complementary caudal-to-rostral distribution is also evident in the population average of retrogradely labeled cell counts in MGB for each section in a series of 9 sections per animal in 8 animals (Fig. 3I). While the entire division of MGBv contained Nissl-stained neurons (Figs. 3B,E, 4A–H), retrogradely labeled cells only occupied a small portion of the MGBv area in each section (Fig. 4I–P). Rostral MGBv neurons that lack retrograde tracer and yet stain for Nissl in adjacent sections (Figs. 3E,F, 4G,O) likely project to neighboring isofrequency contours in A1, as demonstrated in a previous study with multiple retrograde tracers into 4 neighboring cortical isofrequency contours that span over 3 octaves of representation (Storace et al., 2011). An average of 8.5 (±3.2) neurons per animal were double-labeled for both retrograde tracers, accounting for <1% of all labeled neurons across all animals. These data establish that frequency-matched regions of A1 and cSRAF receive input from separate populations of thalamic neurons.
Neurons located outside of the MGBv project more densely to auditory cortical fields outside of A1 in the rat (Arnault and Roger, 1990). Hence, we considered the possibility that the response property differences in cSRAF could be attributed to input from these non-MGBv thalamic divisions. To address this question, the positions of retrogradely labeled neurons were plotted and analyzed for a caudal-to-rostral series of sections spanning ∼1.5 mm, as illustrated for one dual-injected brain hemisphere (Fig. 5A–H). The mean number and percentage of retrogradely labeled neurons in the MGB (Fig. 6A,B, respectively) was similar for injections in A1 (black bars) or cSRAF (gray bars). Most (70%) of the labeled thalamic neurons with injections into A1 and cSRAF combined were located in MGBv, with the remaining located in non-MGBv divisions. Slightly more MGBd neurons were connected to cSRAF; whereas, slightly more SG neurons were connected to A1 (Fig. 6A,B). This confirms that MGBv is the primary source of connections to the target layer of both cortical fields as it typical of a core field, and that non-MGBv connections may serve to distinguish the two thalamocortical pathways.
Previously, VGLUT1 mRNA expression was seen in the cell bodies of the majority of neurons in the MGBv of the rat (Ito et al., 2011). Here, we quantified this pattern of gene expression and asked whether the thalamic neurons expressing VGLUT1 projected to A1 as well as cSRAF. VGLUT1 gene expression was minimal in the caudal portions of the MGBv, the area that contained the greatest number of neurons that projected to cSRAF (Fig. 3D). In contrast, neurons expressing the VGLUT1 gene were highly concentrated in the rostral portion of MGBv (Fig. 3G), the area that contained the greatest number of neurons that projected to A1. The number of neurons projecting to A1 increased ∼100-fold, on average, over the ∼1.5 mm distance from the caudal to rostral MGB (Fig. 3I). Likewise, the number of VGLUT1+ neurons increased by ∼1000-fold, on average, over the same distance (Fig. 3J). The caudal-rostral transition of VGLUT1 mRNA expression is demonstrated in three different animals as photomicrographs (Fig. 4Q–X) and two series of plotted sections (Fig. 5I–P,Q–X). The caudal-to-rostral relationship between MGB neurons that project to A1, and those that express VGLUT1 mRNA is further supported by a significant positive correlation (r = 0.93, p < 0.001), and the presence of <2% of VGLUT1+ neurons in the three sections containing most cells (75%) that projected to cSRAF (Fig. 3I,J, sections 2–4). These data confirm glutamate transporter gene expression changes along the caudal-to-rostral axis of the MGBv.
To determine whether the population of neurons that project to A1 versus cSRAF could be distinguished by the presence of VGLUT1 mRNA, tissue was dual processed to reveal neurons containing retrograde tracer and expressing VGLUT1 mRNA. Dual processed sections are shown from two animals in caudal (Fig. 7A,B) and rostral (Fig. 7C,D) MGB following injections into cSRAF (Fig. 7A,C) or frequency-matched A1 (Fig. 7B,D). Neurons in caudal MGB had low VGLUT1 mRNA expression, and were densely connected to cSRAF (Fig. 7A), but not to A1 (Fig. 7B). In contrast, neurons in rostral MGB had high VGLUT1 mRNA expression, and were densely connected to A1 (Fig. 7D), but not to cSRAF (Fig. 7C). High-magnification images reveal neurons that expressed VGLUT1 mRNA only (Fig. 7E), neurons labeled with retrograde tracer only (Fig. 7F, CTBG), and neurons double-labeled with VGLUT1 mRNA and retrograde tracer (Fig. 7G).
Quantification of labeling patterns in 3 animals revealed that 81% (±5.2) of the neurons that projected to A1 expressed VGLUT1 mRNA (Fig. 7H), and were positioned in rostral MGBv (Fig. 7I). In contrast, quantification of labeling patterns in 2 animals revealed that <3% of the neurons that projected to cSRAF expressed VGLUT1 mRNA (Fig. 7H). The caudal-rostral distributions for experiments double-reacted for VGLUT1 mRNA, and silver enhanced for tracer injection into A1 or cSRAF are plotted (Fig. 5I–P,Q–X, respectively). These results indicate that the majority of neurons projecting to the 8 kHz contour of A1 express VGLUT1, while those that project to frequency-matched cSRAF do not.
Discussion
It is common for mammals including human to have multiple cortical areas with distinct sound frequency response organization and processing and yet the etiology and purpose of these parallel representations is still not fully known (Schreiner et al., 2000; Schreiner and Winer, 2007; Hackett, 2011). Here we find the MGBv division of thalamus is the primary source of connections to layer IV for two frequency-organized fields, A1 and cSRAF, in the rat (Fig. 6). VGLUT2 gene expression is high throughout the caudal-to-rostral extent of MGBv (Barroso-Chinea et al., 2007; Ito et al., 2011) and a high VGLUT2 protein expression in cortical layers III and IV (Fig. 2) suggests both thalamocortical pathways use glutamate as a neurotransmitter (Fig. 8). However, we find that A1 and cSRAF have completely independent representations of sound frequency (Fig. 1A) and are connected to rostral and caudal halves of the MGBv, respectively, with <1% of MGBv neurons projecting to both regions. Furthermore, the majority (81%) of rostral MGBv neurons that project to A1 express the VGLUT1 gene, whereas few (3%) of the caudal MGBv neurons that project to the frequency-matched region of cSRAF express the VGLUT1 gene. Together this indicates functionally independent and parallel MGBv pathways to frequency-matched regions of neighboring frequency-organized auditory cortical fields (Fig. 8, caudal vs rostral pathways).
This study offers a promising new approach to identifying functionally distinct auditory pathways and potential molecular markers for specialization of auditory pathways across development and mammalian species. VGLUT2 gene expression is observed early in development throughout auditory thalamus and expressed minimally in cortical neurons (Fremeau et al., 2004b; Hackett et al., 2011; Ito et al., 2011) and corresponding VGLUT2 protein expression is likely situated on MGB terminals synapsing in cortex. In primate, the density of VGLUT2 protein expression is higher in A1 than surrounding belt and parabelt fields (Hackett and de la Mothe, 2009). Similarly, in the rat we find a higher density of VGLUT2 protein expression in A1 compared with the non-frequency-organized auditory “association” field Te3 located dorsal to the rhinal fissure and ventral to SRAF (Fig. 2C). A finer grain quantitative analysis may reveal additional laminar pattern differences for VGLUT2 protein expression in the rat. In primate, convergent inputs from non-MGBv (VGLUT2 expressing) thalamic pathways is another source of differences across cortical fields (Kaas and Hackett, 2000). Here, we find similar percentages and absolute numbers of MGBv neurons projecting to A1 and cSRAF (Fig. 6). Slightly more MGBd neurons projected to cSRAF; whereas, slightly more SG neurons projected to A1 (Fig. 6). The non-MGB connections with auditory cortex may be underestimated with the approach used here, as deposits were centered in layer IV, not layer III where MGBd synapses are highly concentrated (Huang and Winer, 2000; Smith et al., 2012). The highest density of VGLUT1 gene expression is in rostral MGBv in adult rat and mouse, indicating phenotypic conservation across rodent species (Ito et al., 2011). Preliminary studies suggest a possible phylogenetic difference and lack of VGLUT1 expression in MGBv of adult primate (Ito and Takada, 2011). The latter may be related to other phylogenetic differences including higher concentrations of inhibitory interneurons in MGBv (Winer et al., 1996; Peruzzi et al., 1997) and higher resolution of frequency responses to sound in cats and primates versus rats (Schreiner et al., 2000; Cheung et al., 2001; Read et al., 2002; Polley et al., 2007). These and many other factors likely contribute to variation in vesicular transporter gene and protein expression patterns during development as well as across species.
What might be the physiologic advantage to expressing two glutamate transporters in pathways to A1? Prior studies suggest coexpression of VGLUT1 and VGLUT2 allow for high response magnitudes and rates at corresponding synapses. VGLUT2 is expressed early in development and current theories suggest the highly reliable glutamate release associated with VGLUT2 is needed to establish thalamocortical sensory pathways (Nakamura et al., 2005; Liguz-Lecznar and Skangiel-Kramska, 2007). The sensitivity and magnitude of thalamocortical sensory responses increase during development (Zhang et al., 2001). Corresponding increased neural activity shapes the VGLUT1 and VGLUT2 expression patterns in several systems, indicating a refinement in synaptic processing (De Gois et al., 2005; Yoshida et al., 2009). In adult rats, the majority of MGBv cells that express VGLUT1 also coexpress VGLUT2 (Ito et al., 2011); hence, the population of cells identified here as projecting to A1 and expressing VGLUT1, likely coexpresses VGLUT2 (Fig. 8A). VGLUT1 and VGLUT2 transporters can colocalize to the same synaptic vesicles (Herzog et al., 2006). Furthermore, VGLUT1 and VGLUT2 bind to distinct vesicle trafficking proteins (De Gois et al., 2006), and the combined expression of both transporters results in more transmitter per vesicle, bigger and faster postsynaptic responses and higher spike rate output (Fremeau et al., 2004b; Wojcik et al., 2004; Wilson et al., 2005). In theory, the higher magnitude and faster time course of responses in A1 that we report here could require the coexpression of VGLUT1 and VGLUT2 at corresponding thalamocortical synapses (Fig. 8). The present study demonstrates a gene expression pattern that marks differences in parallel core ascending auditory pathways, and calls for future studies to more directly relate glutamate transporter type to dynamic properties of cortical sound processing.
Footnotes
This work was supported by National Institute of Child and Health Development Grant HD2080 (H.L.R.), National Institutes of Health Grant DC000189 (D.L.O.), and an intramural University of Connecticut Grant (H.L.R.). We thank A. Burkhalter and A. Reyes for comments on an early draft of this manuscript.
- Correspondence should be addressed to Heather L. Read, Department of Psychology, Behavioral Neuroscience Division, University of Connecticut, 406 Babbidge Road, Unit 1020, Storrs, CT 06269. heather.read{at}uconn.edu