Abstract
The mouse sensory neocortex is reported to lack several hallmark features of topographic organization such as ocular dominance and orientation columns in primary visual cortex or fine-scale tonotopy in primary auditory cortex (AI). Here, we re-examined the question of auditory functional topography by aligning ultra-dense receptive field maps from the auditory cortex and thalamus of the mouse in vivo with the neural circuitry contained in the auditory thalamocortical slice in vitro. We observed precisely organized tonotopic maps of best frequency (BF) in the middle layers of AI and the anterior auditory field as well as in the ventral and medial divisions of the medial geniculate body (MGBv and MGBm, respectively). Tracer injections into distinct zones of the BF map in AI retrogradely labeled topographically organized MGBv projections and weaker, mixed projections from MGBm. Stimulating MGBv along the tonotopic axis in the slice produced an orderly shift of voltage-sensitive dye (VSD) signals along the AI tonotopic axis, demonstrating topography in the mouse thalamocortical circuit that is preserved in the slice. However, compared with BF maps of neuronal spiking activity, the topographic order of subthreshold VSD maps was reduced in layer IV and even further degraded in layer II/III. Therefore, the precision of AI topography varies according to the source and layer of the mapping signal. Our findings further bridge the gap between in vivo and in vitro approaches for the detailed cellular study of auditory thalamocortical circuit organization and plasticity in the genetically tractable mouse model.
Introduction
The auditory cortex is functionally organized to contain spatial representations of basic stimulus attributes and stimulus-evoked response properties (for review, see Schreiner and Winer, 2007). Some of these spatially ordered representations, such as tonotopy, can be distributed over hundreds of thousands of neurons, while other features, such as the representations for pitch, sound level, spectral integration, and binaural or temporal modulation preference are clustered into a patchwork of spatially discrete modules grouped by laminae or cortical subregion (Middlebrooks et al., 1980; Schreiner and Mendelson, 1990; Polley et al., 2007; Atencio and Schreiner, 2010a; Bendor and Wang, 2010). The synaptic, neurochemical and cell-autonomous properties that distinguish one cortical module from the next have remained mysterious because this cellular level of analysis is most effectively performed in vitro. Yet, studies of specialized sensory modules require, by definition, an intact connection between the neuron and the sensory receptor epithelium.
Although in vivo whole-cell recordings have begun to uncover the interplay between synaptic excitation and inhibition that shape cortical receptive fields (Wehr and Zador, 2003; Wu et al., 2006; Montgomery and Wehr, 2010; Zhou et al., 2010), a detailed understanding of thalamocortical microcircuits, on the order of what has been achieved in the cochlear nucleus, would require a regional and cell type-specific analysis performed in vitro as well as in vivo (Godfrey et al., 1975; Young and Brownell, 1976; Manis, 1990; Reyes et al., 1994; Oertel et al., 2000). The mouse auditory thalamocortical slice offers a promising approach to understand the cellular bases for regional variations in the representational capacities of auditory cortex neurons. In this acute preparation, the connection between the ventral medial geniculate body (MGBv) and auditory cortex (AI) is preserved, permitting an in vitro analysis of thalamocortical synaptic transmission that can be mechanistically dissected with the genetic toolbox uniquely available in the mouse (Cruikshank et al., 2002, 2010; de la Rocha et al., 2008; Llano and Sherman, 2009). One of the principal aims of this study was to provide a Rosetta Stone for translating the spatial geometry of a sensory feature map delineated in AI and MGBv of intact mice onto the same neural circuitry present in the thalamocortical slice.
Tonotopy, the spatially ordered gradient of preferred sound frequency, is a reasonable first choice to be explored both in vivo and in vitro. Tonotopy is a highly conserved organizational feature documented in AI and MGBv of at least 20 mammalian species over the past 4 decades (for review, see Kaas, 2011). In the mouse auditory cortex, however, the existence of clearly organized tonotopic maps is controversial, with two studies reporting orderly gradients spanning AI and anterior auditory field (AAF) (Stiebler et al., 1997; Zhang et al., 2005) and two recent studies describing a heterogeneous organization of best frequency (BF) in the same fields (Bandyopadhyay et al., 2010; Rothschild et al., 2010). Moreover, the presence or absence of a tonotopic organization in MGBv of the mouse has never been reported.
The goals of the present study were to determine whether a precisely organized tonotopic map exists in MGBv and core fields of the auditory cortex, how the basic pattern of thalamocortical afferents give rise to such tonotopic order (or lack thereof) in the cortex, and finally whether an organizational pattern established in vivo could also be demonstrated in the thalamocortical brain slice. We were surprised to find that spatially organized activity patterns in the slice motivated additional experiments in the intact preparation, which may shed some light on the debate surrounding the existence of precise tonotopy in the mouse auditory cortex.
Materials and Methods
In vivo neurophysiology studies
Neurophysiological data collection.
All procedures were approved by Vanderbilt University and Harvard University animal care and use committees and followed the guidelines established by the National Institutes of Health for the care and use of laboratory animals. Female C57BL/6 mice aged 4–7 weeks were brought to a surgical plane of anesthesia using a combination of pentobarbital sodium (50 mg/kg followed by 10–15 mg/kg supplements as needed) and chlorprothixene (0.2 mg). Multiunit responses were recorded from the middle cortical layers of AI (420–440 μm from pial surface) with epoxylite-coated tungsten microelectrodes (2.0 MΩ at 1 kHz; FHC) and from MGB with 16-channel silicon probes (177 μm2 contact area, 50 μm intercontact separation; Neuronexus Technologies). Frequency response areas (FRAs) were measured with pseudo-randomly presented tone pips of variable frequency (5.5–45.3 kHz in 0.1 octave increments, 20 ms duration, 5 ms raised cosine onset/offset ramps, 600 ms intertrial interval) and level [0–60 dB sound pressure level (SPL) in 5 dB increments] delivered from a free field electrostatic speaker placed 12 cm from the contralateral ear.
Auditory core fields (AI and AAF) were identified by an unmistakable caudal-to-rostral mirror reversal in tonotopy bounded by sites with poor or abruptly shifted frequency tuning. The rostral-to-caudal tonotopic zone is restricted to a narrow swath of cortical tissue with a precise layout that cannot be reliably predicted by coordinates from bregma or vascular landmarks. For MGB recordings, the silicon probe was inserted through the auditory cortex at 15° off the horizontal plane under stereotaxic guidance to match the plane of the section used in tracer reconstruction and thalamocortical slice experiments. To avoid recording from the dorsal division of the MGB, the probe was initially inserted lateral to the auditory core fields, ∼3.5 mm caudal to bregma. The ventral edge of the MGB was identified by documenting the most lateral cortical insertion point that yielded driven responses from the MGB, some 2.5–3.0 mm from the cortical surface. Reconstruction of lesions and electrode tracks confirmed that this corresponded to the ventral shell of the MGB. To target the MGBv and the medial MGB (MGBm), the silicon probe was inserted 0.5 mm medial to this point, a trajectory that reliably corresponded to the center of the MGBv, as evidenced by histologic reconstruction of lesions and electrode tracks.
FRAs were reconstructed across the full rostral-to-caudal extent of the MGB by making successive penetrations rostral and caudal to the starting position (50–100 μm between penetrations) until the recording probe had advanced beyond the caudal and rostral poles of the MGB. FRAs were also compared along the full lateral-to-medial extent of the MGBv and MGBm by documenting variations in response properties across the linear array of contact sites spanning 0.75 mm. In some cases, electrolytic lesions were made at various rostral-caudal positions in the MGB identified with the silicon probe. FRAs were measured at different insertion depths with a tungsten microelectrode, and small lesions were made by passing 0.8 μA of current for 12 s at one or two points of interest along the lateral-to-medial penetration (e.g., the lateral or medial extremes of tone-driven recording sites or reversals in frequency tuning).
Analysis of individual FRAs.
We performed two traditional measurements of preferred frequency tuning: BF, the tone frequency that elicited the greatest number of spikes across all sound levels, and characteristic frequency, the frequency to which the neuron was responsive at threshold. BF values were distributed across a broader range of frequencies, produced more orderly maps, and were completely objective in their definition. BF was therefore judged to be the better choice for preferred tuning in the present study. Onset latency was defined as the first 1 ms poststimulus onset time bin that populated a high-intensity region of the FRA near the BF (>40 dB SPL, ±0.5 octaves). The low- and high-frequency borders of the FRA were then determined at each tone level, and the bandwidth of the tuning was expressed relative to the minimum response threshold.
Reconstruction of tonotopic maps in auditory cortex and thalamus.
Mapping the auditory cortex and thalamus was achieved by carefully marking the position of the recording electrode on a high-resolution photograph of the surface vasculature. For cortical maps, the topographic position was defined by fitting a line that connected the low-frequency BF region in AI and AAF, and calculating the position of each recording site along the line (see Fig. 1C). Topographic variation in BF was well fit with a quadratic function, with the peak of the fit line used to identify the boundary between AI and AAF. For MGB maps, caudal-to-rostral position was expressed relative to the caudal- and rostral-most penetrations that yielded a measurable FRA from two or more contiguous contacts on the silicon probe. In nearly all MGB penetrations, BFs increased, reached a plateau, and then decreased along the lateral-to-medial extent. The most medial recording site in the plateau portion of the function was defined as the BF reversal point. BF values lateral and medial to the reversal point were independently fit with linear functions, and the slope and range of the fit lines were calculated (see Fig. 5A).
Analysis of spatially distributed tone representations.
All AI and AAF recording sites for a given mouse were divided into 18 equally spaced regions of interest (ROIs) along the rostral-to-caudal extent of each field. Similarly, positions of all MGBv recording sites in a given mouse were assigned to one of nine spatial zones spanning the caudolateral-to-rostromedial tonotopic gradient. Next, for each recording site, we determined whether a single tone drawn from a set of evenly spaced frequencies (5.7, 8.6, 13.0, 19.7, 27.9, or 36.8 kHz) and levels (30, 40, 50, or 60 dB SPL) was contained within the FRA. In so doing, we were able to determine the probability that a given spatial zone of AI, AAF, or MGBv contained a given tone/level combination in their receptive fields. Response probability functions were derived from all recording sites within a given brain structure (see Fig. 7B–G). For the purpose of direct comparison between spatial maps generated with voltage-sensitive dye (VSD), the peak and half-maximum boundaries of the normalized response probability functions were calculated for each map individually (see Fig. 7H,I).
Retrograde tracer and lesion studies
Dual retrograde tracer injection.
For these experiments, 4–7-week-old C57BL/6 mice were anesthetized with ketamine (90 mg/kg with 30 mg/kg supplements as needed) and medetomidine (0.5 mg/kg). Following the initial microelectrode mapping, retrograde tracers, cholera toxin β subunit (CTB) conjugated to either Alexa 488 (green, catalog #34775, Invitrogen) or 594 (red, catalog #34777, Invitrogen) were injected into low and high domains of the BF map, respectively (n = 3). Glass pipettes with a 20–30 μm tip diameter containing CTB tracer (1% in 0.1 m phosphate buffer) and back-filled with mineral oil were attached to a Hamilton syringe. Small quantities (0.1–0.2 μl) of each tracer were pressure injected into the middle layers of each BF domain, which subsequently spread to all layers. Upon completion of the injections, the surgical area was closed and the animal allowed to recover for ∼80 h before transcardial perfusion.
Histology.
Animals were perfused through the heart with ∼100 ml of each of the following solutions prepared in 0.1 m phosphate buffer, pH 7.4: 0.9% saline, cold (4°C) 4% paraformaldehyde, and cold (4°C) 4% paraformaldehyde plus 10% sucrose. Following perfusion, the brains were removed and placed in 30% sucrose for 1–3 d. Blocks containing auditory cortex were cut frozen at 40 μm on a sliding microtome in the coronal plane (lesion cases) or in a modified horizontal plane (slice cases) that included A1 and MGB, as described by Cruikshank et al. (2002). Alternating series of sections were processed for Nissl using thionin, cytochrome oxidase (lesion cases), immunohistochemistry for parvalbumin (slice cases), and fluorescent tracers (slice cases) (Hackett and de la Mothe, 2009).
The relative position of retrogradely labeled cell bodies was quantified independently in MGBv and MGBm from each section by fitting a rectangle with each boundary set by the locations of the most caudal, rostral, lateral, and medial cell bodies. The Cartesian coordinates for the red- and green-labeled cells were then expressed relative to these normalized boundaries. For lesion studies, we calculated the minimum distance separating the center of each lesion from the medial boundary separating the MGBv and MGBm.
In vitro voltage-sensitive dye studies
Acute auditory thalamocortical slice preparation.
C57BL/6J mice (postnatal day 16–20) were rapidly decapitated without anesthesia and placed in ice-cold artificial CSF (ACSF). Slices (600 μm thick) were sectioned on a vibratome using previously established methods (Cruikshank et al., 2002). Slices incubated at 35°C for at least 15 min were then maintained at room temperature (20−22°C). The ACSF contained the following (in mm): 125 NaCl, 25 glucose, 25 NaHCO3, 2.5 KCl, 2 CaCl2, 1.25 NaH2PO4, and 1 MgCl2.
VSD imaging and analysis.
Voltage-sensitive dye Di-4-ANEPPS (catalog #D-1199, Invitrogen) dissolved at 0.01 mg/ml in dimethylsulfoxide stock solution was diluted (0.6 μl/ml) in ACSF, producing a final solution of 5 μg/l Di-4-ANEPPS. Slices were incubated (in 30 ml) before transfer to an ACSF (room temperature) chamber on an Olympus MVX10 microscope with 1× 0.25 numerical aperture objective. Excitation light from a shuttered 150 W halogen lamp (MHF-G150LR, Moritex) was bandpass filtered (515/35 nm) and reflected toward the sample by a 570 nm dichroic mirror. Emitted fluorescence was long-pass filtered (590 nm) and imaged using a MiCam Ultima CMOS-based camera (BrainVision and SciMedia).
Fluorescent signals were integrated across different ROIs through spatial averaging within a 125 × 125 μm area. Fluorescence change was normalized to resting fluorescence (ΔF/F0). Discrete sites in the MGBv were activated with an ACSF-filled patch pipette (1 ms long, 0.5 mA square pulse) at 100 ms post-recording onset. Response amplitude was defined as maximal fluorescence change (ΔF/F) per trial at a given ROI and then averaged across 10 trials. Variations between daily preparations were minimized by dividing responses by the maximum change in fluorescence measured in each set of experiments.
Results
Tonotopic organization in core fields of the auditory cortex
The spatial organization of BF was assessed across the caudal-to-rostral extent of AI and AAF with ultra-high-density mapping from 341 multiunit recording sites in 11 mice (N = 188 sites in 6 mice and 153 sites in 5 mice for right and left hemisphere, respectively). Recordings targeted the middle cortical layers (420–440 μm from the pial surface). We observed a clear tonotopic organization for BF with low-frequency BFs of ∼6 kHz forming the caudal boundary of AI, a high-frequency mirror reversal at ∼32 kHz forming the boundary between AI and AAF, and a progressive shift to lower frequency BFs at the rostral border of AAF (Fig. 1A). Note that mice can hear and vocalize up to 100 kHz, yet the highest frequency cortical BF recorded in our entire sample was 34.3 kHz. Our recordings were performed at 4–7 weeks of age, several months before the onset of high-frequency hearing loss in C57BL/6 mice (Willott, 1986). In fact, the paucity of BFs >40 kHz has been routinely reported in the inferior colliculus and auditory cortex of several inbred and outbred mouse strains, and it is thought that these acoustic frequencies may be represented in specialized fields (Stiebler et al., 1997) or, intriguingly, according to their subharmonic cochlear distortion products (Portfors et al., 2009).
Representative FRAs from AI and AAF recording sites show well demarcated V-shaped tuning profiles (Fig. 1B). BF changes were plotted across the caudal-to-rostral expanse of AI and AAF with an average sampling density of 40 recordings per millimeter, or one recording site every 25 μm. We observed an orderly parabolic shift in BF that was well fit by a quadratic function, with the peak used to define the high-frequency boundary separating AI and AAF (Fig. 1C). We then combined normalized maps from individual mice to represent 0 as the AI-AAF reversal point, −1 as the caudal low-frequency boundary of AI, and +1 as the rostral low-frequency boundary of AAF (Fig. 1D). Using this approach, we found a significant increase in BF across the caudal-to-rostral extent of AI (one-way ANOVA, F = 34.9, p < 1.10−6) and a significant decrease in BF along the caudal-to-rostral extent of AAF (F = 6.88, p < 0.005).
Tonotopic organization of MGBv projections to AI
As a next step, we sought to determine whether tonotopy in AI arose from topographically organized MGB projections. To address this possibility, fluorescent retrograde tracers (CTB-red and CTB-green) were injected into a low- and a high-frequency region of AI that had been identified through microelectrode mapping, as described above (Fig. 2A). Following a waiting period of ∼80 h for retrograde transport, the brain was sectioned along the same plane used in auditory thalamocortical slice experiments (Cruikshank et al., 2002), which permitted visualization of MGBv, the thalamocortical axon bundle, and AI (Fig. 2B). In the example shown in Figure 2C, CTB-green was injected into a 7.0 kHz BF region and CTB-red was injected into a 22.6 kHz BF region. Qualitatively, one can appreciate that well segregated AI injection sites were innervated by compartmentalized axon bundles that originated from spatially separable populations of MGB neurons.
The projection patterns were examined in greater detail from this case (Fig. 3A,B) as well as a second case with injections into the low (5.7 kHz) and high (32.0 kHz) BF extremes of the AI map (Fig. 3C,D). In both cases, we observed dense, spatially separable projections from MGBv to each BF region in AI accompanied by a sparser, comingled projection from MGBm (Fig. 3A,C). Cell body locations were plotted against the anatomical boundaries of MGBv and MGBm in multiple sections across the dorsal-to-ventral extent of MGBv. A consistent pattern emerged from all sections in both cases: in MGBv, low-frequency projection neurons were located caudally and laterally, and high-frequency projection neurons were rostral and medial; MGBm neurons projecting to AI were not spatially separable (Fig. 3B,D).
Further quantification revealed that 75% of retrogradely labeled cell bodies were from MGBv, 24% from MGBm, and the remaining 1% from the posterior thalamic complex. Low-frequency projection neurons in MGBv were significantly more lateral (unpaired t test, p < 1 × 10−6) and caudal (p < 1 × 10−6) than high-frequency projection neurons (Fig. 3E). No significant difference in medial-lateral position was observed among low- and high-frequency projection neurons in MGBm (p = 0.12), although high-frequency projection neurons were observed to be significantly more rostral (p < 0.005) (Fig. 3F). Together, these data demonstrate that cortical projections from MGBv are topographically ordered and further suggest that MGBv may be tonotopically organized along a low-to-high, caudolateral-to-rostromedial axis (Fig. 3G).
Neurophysiological dissociation of MGBv and MGBm
Findings from the tracer studies suggest three explicit predictions regarding tonotopic map orientation in the MGBv, as follows: (1) recordings from rostral-to-caudal positions along the lateral wall of the MGBv should yield a progressive downward shift in BF; (2) recordings across the lateral-to-medial extent of the MGBv should yield a progressive upward shift in BF; and (3) there should be no clear tonotopic organization within the MGBm. To test these predictions, we performed high-density multiunit mapping across the lateral-to-medial extent of the MGBv and MGBm using a multichannel silicon probe (50 μm separation between contacts). Penetrations were made throughout the rostral-to-caudal extent of the MGB (50–100 μm between penetrations) at the same angle as the plane of section in the thalamocortical slice (Fig. 4A). Example FRAs shown in Figure 4A illustrate the drop in BF along the lateral wall of the MGBv between rostral, intermediate, and caudal recording locations (left column). Furthermore, regardless of the BF at the most lateral recording position, BF increased across the first few medial recording locations. However, it was equally clear that BFs subsequently decreased as the recording location moved to deeper, more medial areas of the MGB.
Did this high-frequency mirror reversal in BF indicate that recordings had crossed into the MGBm, or did it reveal a BF organization that was not suggested from the tracer studies? To address this question, we identified positions along the rostral-to-caudal extent of the MGB that yielded clear mirror reversals with the silicon probe (Fig. 4B, gray arrow) and then inserted a conventional tungsten microelectrode along the same trajectory (Fig. 4B, white arrow) to make small electrolytic lesions lateral and/or medial to the BF reversal. An example of one such experiment is shown in Figure 4, B and C, where a lateral lesion was made at a well tuned high-frequency site 0.13 mm lateral to the high-frequency reversal (Fig. 4D, top), and a medial lesion was made at the deepest site with tone-evoked activity, 0.21 mm medial to the high-frequency reversal point (Fig. 4D, bottom, E, white crosses). Localization of lesions in coronal sections reacted for cytochrome oxidase revealed that the recording site lateral to the BF reversal was in MGBv and the lesion medial to the BF reversal straddled the border between MGBm and the posterior thalamic complex (Fig. 4C). Indeed, in 13 of 14 lesions, the BF reversal correctly predicted the anatomical boundary dividing MGBv from MGBm (Fig. 4E).
Tonotopy in MGBv and MGBm
Using the high-frequency reversal point as a neurophysiological marker for the boundary between the ventral and medial divisions, we were able to rigorously test each of the three predictions enumerated above. The BF reversal point was identified in each penetration (Fig. 5A, black square), and the slope of the linear fit lines was calculated independently from the rising and falling components of the BF functions (Fig. 5A). Returning to the first prediction, BFs measured along the lateral wall of the MGBv were found to decrease significantly along the rostral-to-caudal extent (N = 43, one-way ANOVA, F = 4.67, R2 = 0.96, p < 0.05) (Fig. 5B). For the second prediction, we found that BF increased significantly from lateral to medial positions within MGBv (F = 26.1, p < 1 × 10−6) (Fig. 5C, gray circles), and the slope of the increasing BF function did not differ according to where recordings were made within the rostral-to-caudal extent of MGBv (F = 1.0, p = 0.55) (Fig. 5D, gray circles).
With regard to the third prediction, despite the extensive comingling of low- and high-frequency AI projection neurons within the MGBm, neurophysiological recordings indicated that BF decreased significantly across the lateral-to-medial extent of the MGBm (F = 2.75, p < 0.005) (Fig. 5C, black squares), and this negative slope did not vary systematically across the rostral-to-caudal extent of the MGBm (F = 0.41, p = 0.97) (Fig. 5D, black squares). Therefore, tonotopy in MGBv closely matched our predictions from the retrograde tracer studies, yet we also found evidence for a mirror reversal decreasing BF gradient within MGBm, which was not predicted by the anatomical data.
Variations in spectral bandwidth and onset latency in auditory cortex and MGB
Analysis of response properties at individual recording sites proved to be another useful tool for distinguishing AI vs AAF and MGBv vs MGBm. FRA bandwidth measured 10 dB above threshold was significantly broader in AAF than AI (0.82 ± 0.03 vs 0.92 ± 0.04 octaves; unpaired t test, p < 0.05) and in MGBm compared with MGBv (1.12 ± 0.04 vs 0.89 ± 0.03 octaves; unpaired t test, p < 1 × 10−6) (Fig. 5E). Consistent with previous observations in mouse (Linden et al., 2003) and rat (Polley et al., 2007), we found that tone-evoked onset latency was significantly shorter in AAF than AI (12.39 ± 0.17 vs 14.79 ± 0.3 ms; unpaired t test, p < 0.0001). Interestingly, onset latencies decreased by 5.3 ms across the BF range in AI (one-way ANOVA, F = 5.98, p < 5 × 10−5) (Fig. 5F), but not in AAF (F = 1.32, p = 0.26). Similarly, onset latency was significantly shorter on average in MGBm than MGBv (9.75 ± 0.17 vs 10.75 ± 0.18 ms; unpaired t test, p < 0.0001) and decreased by 5.3 ms across the BF range in MGBv (F = 24.98, p < 1 × 10−6) without varying significantly in MGBm (F = 0.59, p = 0.7) (Fig. 5F). Note that the latency shift observed between high and low BF recording sites in AI and MGBv cannot be accounted for by the basilar membrane group delay. Tone-evoked onset spikes recorded in single nerve fibers from the basal (high-frequency) versus apical (low-frequency) regions of the mouse cochlea occur at 1.86 ± 0.07 and 2.89 ± 0.11 ms, respectively, which can only account for ∼20% of the shift observed in AI and MGBv (M. C. Liberman, personal communication) (Fig. 5F, black squares).
Topographic mapping in the thalamocortical slice
Using the BF gradients characterized in AI and MGBv of the intact mouse, we turned to the acute thalamocortical slice preparation to determine whether tonotopic organization corresponded to topographic connectivity. Specifically, we predicted that electrical stimulation of the MGBv along the caudolateral-to-rostromedial tonotopic gradient that emerged from tracer and mapping studies would yield a smooth shift in AI activity foci along the caudal-to-rostral, low-to-high BF gradient established from our cortical mapping studies. This prediction could be effectively tested through VSD imaging in the thalamocortical slice by applying a discrete stimulus to the MGBv and measuring subthreshold VSD signal changes across the entire caudal-to-rostral extent of AI simultaneously. VSD signal amplitude was measured across 18 ROIs (125 × 125 μm each) positioned either in layer IV, where VSD response amplitudes were greatest, or in an immediately dorsal region in layer II/III.
In support of our prediction, we observed that rostral areas of AI were most responsive to rostromedial MGBv stimulation (high BFs); whereas, more caudal areas of the AI map were more effectively activated by caudolateral MGBv stimulation (low/mid BFs) (Fig. 6A). Stimulation of MGBv at six loci along the optimized caudolateral-to-rostromedial orientation-induced, broad, overlapping, subthreshold activation profiles that nevertheless shifted significantly across the caudal-to-rostral extent of layer IV in AI (mean ± SEM ROI for site 1 vs 6, 7.6 ± 0.8 vs 12.8 ± 0.8; t test, p < 0.0001) (Fig. 6B). Average AI layer IV responses were normalized for each stimulation site and represented as a color map, with the upward diagonal band representing the point-to-point topographic mapping between MGBv and AI (Fig. 6C). Topographic shifts along the caudal-to-rostral axis were specific to the optimized stimulation orientation, as stimulating six loci along an orthogonal orientation did not significantly shift the activity profiles, consistent with the expected effects of stimulating the MGBv along an isofrequency contour (9.2 ± 1.9 vs 10.0 ± 1.1; p = 0.72) (Fig. 6D,E). As a final test, we returned to the optimized orientation but shifted the ROIs to layer II/III, rather than layer IV. Topographic shifts in layer II/III were weak, but not significant (7.6 ± 0.9 vs 10.4 ± 1.4; p = 0.1) (Fig. 6F,G), in agreement with recent reports of heterogeneous tuning in the upper layers of mouse AI (Bandyopadhyay et al., 2010; Broicher et al., 2010; Rothschild et al., 2010).
Linking maps of preferred frequency to spatially distributed frequency representations
In one view, VSD imaging data from layer IV are in close agreement with AI and MGBv mapping studies: moving the stimulating electrode across the caudolateral-to-rostromedial tonotopic map in MGBv yields an orderly progression of activity foci across the caudal-to-rostral tonotopic map in AI. This finding thus projects a basic sensory map onto the acute auditory thalamocortical slice preparation to enable high-resolution studies of synaptic physiology. However, from another perspective, the BF maps in vivo exhibited relatively precise point-to-point mapping between a tone frequency and its preferred representation in the cortical map; whereas, stimulating a single point in MGBv-evoked VSD activity across the majority of the AI map. In other words, positional shifts in peak VSD signals corroborate the in vivo mapping of AI, but the breadth and overlap of the VSD activity profiles do not. Why might this be?
The fact that VSD is sensitive to subthreshold changes in membrane potential, while microelectrode mapping reflects spiking activity can certainly explain some of the differences between the two datasets (Berger et al., 2007). On the other hand, the discrepancy may also be attributed to the nature of BF maps, which only represent the single preferred frequency for a given recording site, compared with VSD activation maps, which measure the entire spatial spread of population activity evoked by a point stimulus.
To address this difference in conceptual approach and to reconcile the basic differences in the degree of tonotopic mapping precision with these two methods, we reanalyzed our in vivo recordings according to the spatial spread of tone-evoked excitation rather than preferred frequency. Returning to the same mapping experiment shown in Figure 1, we surveyed which recording sites contained 8.6 and 19.7 kHz tones presented at 50 dB SPL within their FRAs (Fig. 7A). Compared with Figure 1A, one can appreciate that tone representations in AI appeared less spatially focused than the orderly BF map, and tone representations in AAF were even broader and more overlapping.
To facilitate direct comparison with the VSD imaging data, we grouped recording sites from each AI and AAF map into 18 ROIs that spanned the tonotopic gradient in each field. We then selected six tone frequencies homologous to the six MGBv stimulation sites and determined the probability that neurons in each ROI would contain a given 50 dB SPL tone frequency within their FRA. Spiking probability profiles in AI were broader than the BF preference maps, but not as broad as the subthreshold VSD maps (Fig. 7B). Tones of different frequency activated tonotopically appropriate regions within AI, yet the response profiles were quite broad, with frequencies <20 kHz activating approximately half of the map with a ≥0.5 probability.
Normalized color maps were created in a similar fashion to those presented in Figure 6 and further demonstrated that spatial representations of tone frequencies in AI were broad, but spatially organized with a narrow upward diagonal band that was fairly stable across tone levels (Fig. 7E). Compared with AI, spatial representations of tones in AAF were degraded and level intolerant (Fig. 7C,F), which was to be expected given the greater BF scatter (Fig. 1C,D), compressed range of tone frequency representation (Fig. 1C,D), and broader frequency tuning (Fig. 5E). We also investigated spatial activation profiles in the MGBv across nine ROIs positioned along the caudolateral-to-rostromedial BF gradient. Qualitatively, the spatial code for tone frequency representation in MGBv was less organized and less level tolerant than AI, but superior to AAF (Fig. 7D,G).
Spatial activation profiles were further quantified in AI on an individual mouse basis by measuring the peak and width at half maximum of the normalized response functions for each frequency summed across 30–60 dB SPL. We observed significant spatial shifts in the caudal and rostral boundaries at half maximum (Friedman nonparametric ANOVA, p < 0.005 for both) as well as the peak position (p < 0.01) with tones of increasing frequency (Fig. 7H). Overall, each tone activated one-third to one-half of AI ROIs. As a final step, we directly compared spiking half-height boundaries with the normalized VSD half-height boundaries and the average BF values across all ROIs (Fig. 7I).
The upward diagonal slant present in all three measures demonstrates that layer IV of AI is tonotopically organized regardless of whether the input signal to the mapping function is subthreshold VSD, suprathreshold spiking, or BF. However, the precision of the layer IV mapping function varies between techniques. Subthreshold activity profiles suggested a coarse mapping specificity wherein any point in the MGBv map can activate nearly any point in the AI map. Suprathreshold spiking patterns produce a much higher degree of mapping precision, wherein a given tone can recruit activity across a swath of the AI map, yet only a small cluster of neurons within this active zone will claim this tone as the BF.
Discussion
The aims of these experiments were to resolve the discrepancy surrounding the existence of tonotopy in the mouse auditory cortex, to provide the first evidence for or against the existence and orientation of a tonotopic map within the mouse MGB, and to determine whether the tonotopic organization in both brain areas could be captured within the thalamocortical slice in vitro. We determined that the middle layers of AI exhibited a clear tonotopic organization with BFs ranging from 5.6 to ∼34 kHz between the caudal and rostral boundaries, respectively. Retrograde tracer injections into low and high BF domains in AI revealed topographically separable MGBv projections combined with a weaker, topographically mixed projection from MGBm. The relative position of retrogradely labeled cell bodies suggested a low-to-high, caudolateral-to-rostromedial BF organization in MBGv, which was subsequently corroborated with dense two-dimensional receptive field mapping. Stimulating the thalamocortical slice along the same tonotopic gradient in MGBv confirmed the expected caudal-to-rostral shift in VSD activity foci along the AI tonotopic gradient in layer IV (Fig. 8A).
Toward a cellular analysis of regional specialization in the auditory thalamocortical circuit
The combination of in vivo and in vitro approaches will facilitate detailed cellular analyses of the mechanisms supporting specialized auditory signal processing within subdomains of the tonotopic map. For example, do variations in first spike latency and tuning bandwidth across the tonotopic gradient reflect differences in local circuit properties or cell-autonomous regulation of dendritic morphology, ligand-gated or voltage-gated ion channels, as have been reported in the auditory brainstem nuclei (Smith and Rubel, 1979; Sanes et al., 1987; Li et al., 2001; Sanchez et al., 2010; Strumbos et al., 2010)? The spatial extent of a tone representation within the AI or MGB tonotopic map can also be modified through learning, trauma, or normal age-related processes (for review, see Dahmen and King, 2007). It would be interesting to explore whether tonotopic remapping induced by high-frequency age-related hearing loss or focal lesions of the basilar membrane arise through regionally restricted alterations in excitatory–inhibitory balance, consistent with the global dysregulation following pervasive sensorineural hearing loss (Robertson and Irvine, 1989; Willott et al., 1993; Kamke et al., 2003; Kotak et al., 2005, 2008).
Similarly, numerous studies have found that passively rearing young animals in ambient sound environments dominated by a single tone frequency, or conditioning adult animals to associate a neutral tone frequency with a behaviorally relevant outcome, induces large-scale reorganization of the AI tonotopic map (Polley et al., 2006; de Villers-Sidani et al., 2007; Dorrn et al., 2010). Developmentally and spatially specific changes of intracellular and synaptic signaling that permit this reorganization in an otherwise stable map remain unclear. By establishing a tonotopic map in the thalamocortical slice, we can begin to explore how differences in cell morphology, gene expression, synaptic dynamics, and intrinsic physiology may give rise to regional specializations in the representation and plasticity of auditory stimuli.
What is a tonotopic map?
The history of tonotopic mapping studies in AI is reminiscent of the Indian parable of the three blind men who each examine a different part of an elephant's body with their hands (e.g., the trunk, a leg, and the tail). Each man infers the appearance of the entire elephant based on their restricted experience, but upon comparison discovers that his view is highly divergent and mutually incompatible with the others. The presence or absence of a tonotopic organization in AI has been in dispute since the earliest published accounts of electrophysiological recordings in the auditory cortex. Some initial studies reported that AI is tonotopically organized due to the similarity of BFs within a single vertical penetration and the smooth systematic variation of BF across the expanse of AI (Woolsey and Walzl, 1942; Tunturi, 1962; Merzenich and Brugge, 1973), while others have found that tonotopic organization in AI is weak or nonexistent based on observations that BFs could vary considerably within a single penetration (Evans et al., 1965) with no discernable gradient across AI (Evans and Whitfield, 1964; Goldstein and Abeles, 1975).
Contradictory findings continue to the present day, with a pair of elegant in vivo calcium imaging studies reporting that BFs between neighboring neurons in superficial layers of mouse AI can vary by several octaves and form a very coarse tonotopic organization only when compared over the entire extent of AI (Bandyopadhyay et al., 2010; Rothschild et al., 2010). The authors of these studies suggest that earlier reports of fine-scale tonotopy could be attributed to the comparatively coarse spatial sampling used by microelectrode mapping studies (at 100–200 μm intervals), pooling of heterogeneous local tuning through multiunit recording and a selection bias toward responsive recording sites. However, the fine-scale tonotopy reported here was acquired with higher spatial resolution (25 μm sampling density along the tonotopic axis) and would be virtually impossible to obtain through a random selection of locally heterogeneous tuning. Moreover, the electrode tip was inserted into the thalamorecipient layers (420–440 μm beneath the pia) without any systematic fine tuning to find auditory-driven sites. Last, even though AI FRAs were derived from multiunit clusters, the tuning was considerably narrower and more V shaped than what would be expected from the combination of scattered local tuning described in the imaging studies. The precisely organized neurophysiological maps and topographically organized anatomical projections from MGBv leave no doubt that mouse AI is tonotopically organized.
Much like the story of the blind men and the elephant, we would argue that the interpretations of both datasets are individually correct but ultimately fail to adequately describe the whole. In our own experiments, we found that tonotopic precision was considerably degraded when maps were derived from a population-based suprathreshold spiking analysis (Fig. 7), even further blurred when assessed with subthreshold VSD signals in layer IV (Fig. 6B,C) and not statistically significant based upon VSD measurements from layer II/III (Fig. 6F,G). It is our view, therefore, that AI is both tonotopically and nontonotopically organized, depending on the manner in which it is studied.
We would argue that tonotopic mapping precision depends on the following five principle factors: (1) whether recordings were made within AI versus other fields within the core, belt, parabelt/multimodal heterarchy (Hackett, 2011); (2) sound level used for the preferred frequency characterization; (3) anesthetic state; (4) spatial resolution of the measurement technique; and (5) laminae in which measurements are made (Fig. 8B). Studies that characterize preferred frequency from multiunit recordings in the middle layers of AI in barbiturate-anesthetized animals with near-threshold stimuli invariably find a high degree of tonotopic precision. Studies that deviate from this combination by, for example, recording calcium changes from isolated single units in the superficial layers of awake animals or measuring blood oxygenation-level dependent signals from awake subjects with a single high-intensity sound level are less likely to observe a tonotopic map.
In this scheme, tonotopically organized BF maps in AI, which predominantly reflect topographically organized MGBv projections into layers IIIb and IV, are the “tip of the iceberg” for frequency selectivity (Figs. 2, 3) (Winer et al., 1999; Lee et al., 2004; Budinger and Scheich, 2009). MGBv afferents provide the principal source of subcortical input to the developmentally mature AI (Romanski and LeDoux, 1993; de Venecia et al., 1998; Cetas et al., 1999; Huang and Winer, 2000; Kimura et al., 2003) by forming efficacious excitatory synapses onto stubby dendritic spines within 100 μm of the AI cell bodies (Richardson et al., 2009). The present and prior studies have shown that MGBv neurons project to isofrequency regions of the AI map and that AI neurons in layer IIIb/IV inherit their BF from a focal convergence of MGBv neurons with frequency preferences that fall within one-third of an octave (i.e., one critical band) (Imig and Morel, 1984; Winer et al., 1999). However, AI frequency preference can also reflect less tonotopically organized inputs arising from local and long-range corticocortical connections (Kaur et al., 2004, 2005; Liu et al., 2007; Lee and Sherman, 2008; Wu et al., 2008; Happel et al., 2010; Moeller et al., 2010). The influence of these sources on cortical frequency tuning may be accentuated when recordings are made from the upper cortical layers (i.e., the source of Ca2+ imaging signals), where response latencies are longer, spectrotemporal receptive fields more complex and labile, and horizontal corticocortical connections more abundant (Barth and Di, 1990; Ojima et al., 1991; Wallace et al., 1991; Kaur et al., 2005; Barbour and Callaway, 2008; Dahmen et al., 2008; Atencio and Schreiner, 2010a,b).
As suggested from the slice data presented here and demonstrated through whole-cell recordings in the intact animal, the range of tone-evoked subthreshold inputs is substantially broader than the spiking receptive field (Wehr and Zador, 2003; Zhang et al., 2003; Kaur et al., 2004). These unexpressed synaptic inputs provide a latent substrate for remodeling the frequency preference of AI neurons, changes that can occur within minutes and last for many hours to days (Bakin and Weinberger, 1990; Fritz et al., 2003; Froemke et al., 2007; Dorrn et al., 2010). Similarly, recordings made under anesthesia, particularly in the upper layers, are likely to significantly underestimate modulatory inputs that arise from local circuits and higher cortical areas and overexaggerate the normally limited contribution of tonotopically organized MGBv inputs (Kuhn et al., 2008; Fritz et al., 2010; Sadagopan and Wang, 2010). The diverse and dynamic set of potential input sources interacting within and between cortical columns in the unanesthetized cortex provides a means for neighboring neurons to assume substantially different, yet sharply tuned frequency preferences. Thus, the same set of AI neurons could be both tonotopically and nontonotopically organized, depending on the source of synaptic inputs that dominate under a given experimental approach.
Although tonotopically organized inputs can be de-emphasized in some cortical operating modes, it would be incorrect to conclude that topographically organized feature representations are merely an artifact of anesthesia and are unworthy of further study. Regionally specialized cortical modules in AI can highlight areas with strong mutual connectivity (Matsubara and Phillips, 1988; Read et al., 2001), specialized afferent input from the MGB (Middlebrooks and Zook, 1983; de la Mothe et al., 2006; Read et al., 2008), or precise corticofugal connections with neurons in the thalamus (Zhang and Yan, 2008), midbrain (Zhang et al., 1997), brainstem (Luo et al., 2008), and even indirectly to outer hair cells (Xiao and Suga, 2002). Furthermore, reorganization of cortical spatial maps and modules can be a sensitive index for developmental critical periods (Insanally et al., 2009; Popescu and Polley, 2010), perceptual learning (Polley et al., 2006), or functional status of the auditory periphery (Irvine and Rajan, 1996). Therefore, a deep understanding of AI circuitry will likely require a synthesis of complementary approaches.
Footnotes
This work was supported by National Institutes of Health Grants DC009836 (D.B.P.) and DC04318 (T.A.H.); Core Research for Evolutional Science and Technology, Japan Science & Technology Agency (T.K.H.); the Human Frontiers Science Program (T.K.H.); and the Harvard Society of Fellows (T.R.B.). We thank Michelle Young and Dr. Vivek Khatri for technical assistance and fruitful discussion.
- Correspondence should be addressed to Dr. Daniel B. Polley, Eaton-Peabody Laboratory, Massachusetts Eye and Ear Infirmary, 243 Charles Street, Boston, MA 02114. Daniel_polley{at}meei.harvard.edu