Introduction

Real-world listening occurs in environments cluttered by noise. To navigate these conditions, the auditory system relies on the integration of multiple sources of input. Broad integration of neural inputs, however, comes at a potential cost. Dense convergence can lead to inappropriate cross-talk between input streams, generating a processing “soup” whereby individual neural pathways cannot be independently modulated. Therefore, highly connected neural systems need to preserve the specificity of inputs in the face of massive integration (Fair et al., 2007; Sporns, Tononi, & Edelman, 2000; Tononi, Sporns, & Edelman, 1994)

Several brain systems have solved the integration/segregation problem by creating repeated structural motifs that allow independent modulation of different pathways. For example, in the basal ganglia, neurochemical periodicity of the striosome and matrix regions permits the retention of functional independence of input streams in the face of massive convergence. As a result, basal ganglia direct and indirect pathways may be independently modulated, while allowing controlled crosstalk between channels (Alexander, Crutcher, & DeLong, 1990; Brimblecombe & Cragg, 2017; Draganski et al., 2008; Gerfen, 1985)

In the auditory system, the integration/segregation challenge is particularly formidable in the inferior colliculus (IC). The IC serves as a specialized integration hub because of its heavy innervation by higher auditory centers as well as by nonauditory regions. Thus, the IC stands at the crossroads between multiple brain regions and is thus important for linking acoustic stimuli with non-acoustic inputs. The IC comprises a lemniscal division, the central nucleus (ICC), that receives primarily ascending auditory projections (Oliver, 2005; Willard & Ryugo, 1983; Willott, 2001) and two non-lemniscal divisions known as the dorsal (DC) and lateral (LC) cortices that receive massive descending auditory as well as non-auditory input (Bajo & Moore, 2005; Coleman & Clerici, 1987; Loftus, Malmierca, Bishop, & Oliver, 2008; Saldaña, Feliciano, & Mugnaini, 1996; Schreiner & Winer, 2005; Winer, Larue, Diehl, & Hefti, 1998) as reviewed in (A. M. H. Lesicko & Geffen, 2022).

The LC and DC appear to differ in their strategy for integrating these inputs. The DC does not contain known structural heterogeneities but has a tonotopic arrangement of neuronal responses on its surface (Barnstedt, Keating, Weissenberger, King, & Dahmen, 2015; A. B. Wong & Borst, 2019). In contrast, the LC contains repeated neurochemical motifs, here referred to as modules, that express the inhibitory neurotransmitter GABA as well as other metabolic markers (Chernock, Larue, & Winer, 2004). Areas outside the modules, referred to as matrix, stain strongly for calretinin (Dillingham, Gay, Behrooz, & Gabriele, 2017). We and others have determined that nearly all of the long- and short-range connectivity of the LC is governed by the module/matrix organization (Dillingham et al., 2017; Lamb-Echegaray, Noftz, Stinson, & Gabriele, 2019; A. M. H. Lesicko, Sons, & Llano, 2020; A. M. Lesicko, Hristova, Maigler, & Llano, 2016).

Unfortunately, how the LC integrates these inputs while maintaining distinct paths of information flow is not known. For example, LC neurons respond to multisensory stimuli (Jain & Shore, 2006; Zhou & Shore, 2006), but information from different modalities enters the LC in segregated channels (A. M. Lesicko et al., 2016). Thus, integration across the module/matrix border must occur, but the mechanisms of this integration are unknown. Traditional methods of studying the LC, such as in vivo electrophysiology, pose significant challenges in unambiguously assigning individually recorded neurons to residence in matrix or modules. Imaging approaches, such as two-photon (2P) microscopy, have proven to be enormously valuable in understanding the fine spatial structure of neural response properties in the auditory system (Bandyopadhyay, Shamma, & Kanold, 2010; Romero et al., 2020; Rothschild, Nelken, & Mizrahi, 2010). However, because of its superficial location on the dorsal surface of the mouse midbrain, the DC is the only part of the IC that has been characterized by 2P imaging (Barnstedt et al., 2015; Ito, Hirose, Murase, & Ikeda, 2014; A. B. Wong & Borst, 2019). Optically, it has been impossible to characterize the functions of the LC, which is laterally embedded deep in a brain sulcus.

Here, we describe a novel approach using a 45° microprism inserted into the sulcus that separates the LC from the cerebellum to investigate the fine-scale functional organization of the LC from its lateral surface through 2P imaging. The 2P images obtained from the microprism revealed the first-ever images of the LC in vivo showing its characteristic GABAergic modules. Compared to the DC imaged from the surface, the acoustically and somatosensory-evoked calcium signals obtained by the 2P imaging of the LC through the microprism revealed functional distinctions between the different cellular motifs of the LC. For instance, while modules were more responsive to somatosensory stimulation than matrix, the latter was more acoustically-responsive and had a lower detection threshold for complex sounds. Additionally, modules and matrix were found to differently process the spectral and temporal features of the acoustic information. These findings suggest that the non-lemniscal divisions of the IC contain rich but distinct platforms for the processing of different features of auditory and non-auditory information.

Results

The insertion of 45° microprism for imaging the LC from the side

GAD67-GFP knock-in mice on a Swiss Webster background were crossed with Thy1-jRGECO1a mice to produce a GAD67-GFPx jRGECO1a hybrid mouse, which was used to visualize GABAergic modules, to identify the dense network of GABAergic neurons and terminals in modules of the LC in the green channel (Tamamaki et al., 2003), and the calcium indicator, jRGECO1a, in the red channel (Dana et al., 2018). A standard craniotomy was performed (Goldey et al., 2014) with modifications to suit the specific characteristics of the IC in terms of its location, depth, and high density of the surrounding vasculature (Xiong et al., 2017). The craniotomy (∼ 3 x 4 mm) did not drill directly into the sagittal and lambdoid sutures which allowed minimal disruption of the superficial sinuses below the skull (Figure 1A) as indicated by the minimal bleeding near the bone flap (Figure 1B). During drilling, pressurized air and a cold stream of saline were used to remove the bone debris and dissipate heat. Because the IC sits deeper in the brain compared to the surrounding structures, the direct insertion of the 45° microprism was found to push the IC ventrally. Therefore, a fine diamond knife was used to separate the lateral surface of the IC from the surrounding tissues without damaging the transverse and sigmoid sinuses (Figures 1C and D). This step was followed by a gentle retraction of the knife medially to generate a pocket to secure the microprism (Figure 1E). The successful insertion of the microprism was assessed by examining the reflection of the diamond knife from the microprism’s top glass surface (Figures 1F and G). Holding the microprism down and retracting the knife upward simultaneously resulted in securing the microprism in its pocket (Figure 1H and I). Locking the microprism down in its pocket for 10 minutes with the diamond knife (Figure 1J) is critical to ensure that the microprism is stabilized and to prevent bleeding. Filling the gap between the skull surface and the microprism with saline was found to generate motion artifact due to the heart pulsation, so 1% agarose was used to fill the gap to reduce the motion artifact and add more stability to the microprism (Figure 1K). To compare response properties between the LC and DC, direct imaging (i.e., without a microprism) of the dorsal surface of the IC, where the DC is located, was also conducted by scanning the DC surface from the most medial to lateral regions following the convex curvature of the surface (Figures 1L and M).

Figure 1:

The 2P characterization of GABAergic and non-GABAergic neuronal distribution on the LC and DC surfaces with histological validation

Following the surgery, the anesthetized animals were taken immediately for imaging. Examination of the distribution of GABAergic modules in the LC was critical to validate the target area for functional calcium imaging. Therefore, the 920 nm laser beam was initially used to scan for GFP, which is expressed in the GABAergic cells of GAD67-GFP knock-in mice. The insertion of the microprism was initially viewed by low magnification fluorescence imaging of the GFP signals from the whole surface of the IC along with the inserted microprism (Figure 2A). This step validated the location of the microprism relative to the IC under the microscope and ensured there was no bleeding that could block the field of view. The 2P images of the GFP signals from the LC imaged via the microprism, or for simplicity LC_(microprism), across different animals revealed the LC characteristic GABAergic modules (Figures 2C and E, the irregular white lines), which were consistent with previous 3D reconstructions of modules from the coronal sections (Chernock et al., 2004; A. M. Lesicko et al., 2016). To further validate the microprism location and the ability to image the LC from the lateral surface, a laser lesion (see Methods) was created via the microprism (Figure 2G), which was validated histologically (Figure 2H). For comparison, the whole dorsal surface of the IC was scanned from medial to lateral borders (Figure 2B, the dotted white lines). Unlike the LC, most of the dorsal surface of the IC showed a non-clustered distribution of the GABAergic cells (Figures 2D and F, blue boxes). However, the lateral horizon of the IC dorsal surface showed clusters of GABAergic cells and terminals that were confirmed to be modules, which comprise the prominent feature of the LC (Figures 2D and F, red boxes). To further validate the presence of GABAergic modules in the dorsolateral surface of the IC, a laser lesion was created in a GABAergic module located in that area (Figures 2I and J, dotted red circle), which was confirmed histologically (Figure 2K), demonstrating that it is possible to image a portion of the LC directly from its dorsal surface. Using the 1040 nm laser, the GAD67-GFPx jRGECO1a hybrid mouse showed a non-clustered cellular expression of jRGECO1a on the surface of the DC (Figure 3A) or on the surface of the LC_(microprism) (Figure 3B), as expected. The IC of anesthetized animals in each experimental setup was then imaged to measure calcium signals over the course of pseudorandom presentations of either 500 ms pure tones of different combinations of frequencies (5-40 kHz, half octave intervals) and sound levels (40-80 dB SPL, 10 dB intervals) or broadband noise without or with amplitude modulation (AM-noise) of different combinations of AM frequencies (2-256 Hz, one-octave intervals) and sound levels (40-80 dB SPL, 10 dB intervals). As shown by two representative cells from the DC imaged from the dorsal surface or LC imaged via microprism, these acoustic stimulation paradigms successfully elicited transient cellular calcium signals that varied depending on the stimulus type (Figures 3C and D). The expression distributions of jRGECO1a were histologically examined for the 2P imaged areas of the DC and LC (Figure 3-figure supplement 1A). The DC and LC, have a similar fraction of cells expressing jRGECO1a (paired sample t-Test, t₍₉₎ = 0.98, p = 0.34, n = 10 coronal sections from 2 animals) (Figure 3-figure supplement 1B and C). Also, the cells expressing jRGECO1a in all regions of the IC were only nonGABAergic neurons as they did not express GFP which is expressed exclusively in the GABAergic neurons (Figure 3-figure supplement 1D).

Figure 2:

Figure 3:

Figure 3-figure supplement 1:

Acoustic responsiveness differs in the LC and DC

As indicated by the GFP signals and consistent with Figure 2, imaging the LC of the hybrid mice via microprism showed clusters of GABAergic modules, which take different irregular shapes on the surface of the LC (Figure 4A - irregular white lines). Based on the best pure tone frequency (BTF) of each responsive cell (Barnstedt et al., 2015), an inconsistent and fragmentary tonotopic organization was seen over the LC_(microprism) (Figure 4B) as indicated by a low R² of the best linear and quadratic regressions fit made between the BTFs of cells and their locations (Figure 4 - figure supplement 1A). Additionally, the best fit of these regressions was found to be along different axes across the tested animals (Figure 4 - figure supplement 1A). Despite this inconsistent topographical organization, it was found that the cells in the modules were tuned to relatively higher frequencies than those in the matrix (Mann-Whitney Test: z = 4.03, p = 5.6×10^-5, a median of best-tuned frequency = 10 and 7.1 Hz for modules (215 cells from 7 animals) and matrix, (230 cells from 7 animals) respectively).

Figure 4:

Figure 4-figure supplement 1:

For comparison, the surface of the DC was imaged directly from the dorsal surface using the routine imaging setup. To view the majority of the DC surface, imaging was done by scanning the DC from its medial to the lateral horizon to generate a composite map (see Methods). As indicated by the GFP signals, it was noted that the most lateral region of the surface of the IC showed clusters of the GABAergic cells and terminals (i.e., modules, Figure 4C – irregular white lines) and are therefore considered part of the LC. Based on the BTF of each responsive cell (Barnstedt et al., 2015), the DC showed a tonotopic organization of cells along the rostromedial to caudolateral axis tuned to lower or higher frequencies, respectively, (Figure 4D) as indicated by a high R² of the regression fits made between BTFs of cells and their locations along the rostromedial to caudolateral axis. (Figure 4-figure supplement 1B). It was noted that the R² of the nonlinear quadratic regression was higher than that of the linear fit. The additional quadratic fit allowed the regression model to include the cells tuned to higher frequencies at the rostromedial part of the DC and result in a better fit, which is consistent with the tonotopic organization that was previously described (Barnstedt et al., 2015; A. B. Wong & Borst, 2019). It was also noticed that a large area of the DC was tuned to the lowest tested frequency (5kHz), suggesting a low-frequency bias on the surface of the DC. To quantify tuning differences between DC and LC, the cumulative distribution function based on the cells’ BTFs showed a low-frequency bias in the DC, similar to previous work (Ito et al., 2014; Romand & Ehret, 1990; Shen, Xu, & Shen, 2003; Stiebler & Ehret, 1985). Although around 50% of the cells of LC_(microprism) showed also a low-frequency bias to 5 and 7.1 kHz, collectively, most of the LC cells were tuned to relatively higher frequencies than those of the DC (Figure 5A, Mann-Whitney Test: z = -11.34, p < 0.001).

Figure 5:

Detection of AM is critical to identify natural sounds as well as to aggregate sound sources during auditory scene analysis, which is important for communication (Joris, Schreiner, & Rees, 2004; Rees & Langner, 2005; Singh & Theunissen, 2003). The IC also is the first site in the ascending auditory pathway where a substantial number of neurons are rate-tuned to AM (Frisina, 2001; Joris et al., 2004; Langner, 1992; Nelson & Carney, 2007; Rabang et al., 2012). Even though calcium indicators do not provide a spike-by-spike readout of neural activity, their integrated signal does correlate with the number of spikes produced by a neuron over short time windows (Vogelstein et al., 2010; Yaksi & Friedrich, 2006). Thus, this integrated signal can be used to measure rate responses to AM. Using this integrated signal, we tested the ability of the cells of the LC _(microprism) to respond to broadband noise at different AM rates within each cellular motif. Within the fraction of responsive cells, the LC_(microprism) showed cells of different response profiles. While some cells responded only to pure tones or tone-selective (Tone-sel), other cells responded only to AM-noise or noise-selective (Noise-sel). However, some cells responded to both pure tones and AM-noise (Non-sel), when the two stimuli were individually presented. While the LC_(microprism) had similar fractions of Noise-sel and Non-sel cells, both fractions were significantly higher than that of the Tone-sel cells (Figure 5B, one-way-ANOVA, p = 1.6×10^-4). We also observed that modules and matrix differentially process the spectral information. For instance, the matrix had a significantly greater fraction of Noise-sel cells compared to the modules (Figure 5C, Mann-Whitney Test, z = -2.32, p = 0.02). In contrast, both cellular motifs had similar fractions of Tone-sel and Non-sel cells (Figure 5C, Mann-Whitney Test: z = 0.48, p = 0.62 for Tone-sel and z = 1.36, p = 0.17 for Non-sel). These findings suggest that neurons in the matrix are more responsive to sounds of greater spectral and temporal complexity than those in the modules. Similar to the LC, the DC showed significantly higher fractions of Noise-sel and Non-sel cells than Tone-sel cells, but the cellular fraction of Non-sel cells was higher than that of Noise-sel cells (Figure 5D, one way-ANOVA, f_(2,9) = 86.9, p = 1.2×10^-6).

By pooling the cells based on their response to either pure tones or AM-noise, animals showed a lower fraction of responsive neurons to the acoustic stimulation in the LC _(microprism) compared to the DC regardless of the sound type (two-way ANOVA, f_(1,1,18) = 9.1, 7.8, and 0.1, p = 0.007, 0.01, and 0.75 for region (LC vs DC), stimulus (pure tones vs AM-noise), and interaction, respectively, n = 7 animals for LC and 4 animals for DC, the % of responsive cells ± SEM = 48±7% for LC and 75±2% for DC). To evaluate if any optical artifacts due to imaging through a prism could be responsible for this finding, the number of responsive cells imaged from the lateral-most dorsal surface of the IC, which we consider part of the LC because of the presence of GABAergic modules (LC imaged from the dorsal surface or LC_{(dorsal surface)}), was quantified and compared with that obtained from the LC_(microprism). We observed that there were no significant differences between the number of tone or noise-responsive neurons between the LC_(microprism) or LC_{(dorsal surface)} (two-way ANOVA, f_{(1,1, 18)} = 0.07, 14.1, and 2.11, p = 0.79, 0.001, and 0.16 for the region, stimulus (pure tones vs AM-noise), and interaction, respectively, n = 7 animals for LC_(microprism) and 4 animals for LC_{(dorsal surface)}, the % of responsive cells ± SEM = 48±7% for LC_(microprism) and 45±11% for LC_{(dorsal surface)}).

LC and DC respond differently to spectral and temporal modulation of acoustic inputs

Given that more cells from LC _(microprism) and DC had responses to AM-noise compared to pure tones, we further quantified the selectivity to spectral or temporal modulation by computing the spectral modulation index (SMI) or the temporal modulation index (TMI), respectively, as well as the number of modulated cells for all the responsive cells (see Methods) across the tested sound levels.

For the LC_(microprism), the matrix has a greater proportion of spectrally modulated cells (cells responding only to pure tones and unmodulated noise) than modules. The higher fractions of modulated cells in matrix were most prominently at sound levels below 70 dB (Figure 5E, two-way ANOVA: p = 7.1×10^-7, 1.0, and 0.01 for the region, sound level, and interaction, respectively) suggesting that the cells in the matrix have a lower threshold to detect the spectral features of sound. Similarly, the matrix had more cells responsive to temporal modulation (cells responding to AM-noise) as a main effect (Figure 5F, two-way ANOVA, p = 2.7×10^-4, 1.0, and 0.32 for the region, sound level, and interaction, respectively), which was consistent with the above results (Figure 5C).

The DC and LC imaged via microprism were found to have different profiles of modulation of responses to spectral and temporal complexity. Although there is a general trend for SMI and TMI to decrease with increasing sound-pressure level, corresponding to a decrease in the response to pure-tone and AM-noise, respectively, overall there is a statistically-significant higher level of average SMI and TMI in the LC compared to the DC (Figures 5G and H, SMI: two-way ANOVA, p < 0.001 for all region, sound level, and interaction, TMI: two-way ANOVA, p < 0.001 for all regions, sound level, and interactions). These results suggested that neurons in the LC are more sensitive to spectral and temporal modulation than those in the DC.

LC modules and matrix differentially respond to acoustic and somatosensory stimuli

Somatosensory inputs from many cortical and subcortical structures have been identified in the non-lemniscal division of the IC, particularly in the LC (A. M. Lesicko et al., 2016; Lohse, Zimmer-Harwood, Dahmen, & King, 2022; Zhou & Shore, 2006) Previous work has shown that while somatosensory projections mostly target the modules, the auditory projections tend to target the matrix (A. M. Lesicko et al., 2016), which indicates that the LC is a multisensory integration site. Therefore, we examined the somatosensory and auditory cellular responses on the surface of the LC by coupling the 2P imaging with the microprism with a whisker or acoustic stimulation (see Methods). Acoustic stimulation was done with 500 ms of unmodulated broadband noise at 80 dB SPL. After the acoustic or somatosensory stimulation, evoked cellular responses were observed from the neurons of the LC_(microprism) as well as the DC imaged directly from the dorsal surface. According to these cellular responses, different types of neurons were found in modules and matrix of the LC_(microprism) as well as the DC (Figure 6A-D). Also, the auditory response evoked by auditory stimulation only was found to be modulated (suppressed vs enhanced) when the auditory stimulation was simultaneously presented with somatosensory stimulation (bimodal stimulation) (Figure 6E and F). The cells of every structure were then mapped according to their responses. Auditory and somatosensory responses were evoked in the neurons of both matrix and modules imaged via microprism (Figures 7A, B, and D). Within the pool of responsive cells, the cells were categorized into cells that respond only to auditory stimulation [auditory selective cells (Aud-sel)], only to somatosensory stimulation [somatosensory selective cells (Som-sel)], or nonselectively to both stimuli (Aud/Som-nonsel)] (Figure 7F). Both modules and matrix showed a higher fraction of Aud-sel cells than that of Som-sel and Aud/Som-nonsel cells (Figure 7H, Chi-Square test, χ² = 33.1, p = 3.5×10^-6). Pooling the cells based on their auditory vs somatosensory responses revealed that more auditory-responsive cells (Aud-sel + Aud/Som-nonsel) were found in the matrix, while more somatosensory-responsive cells (Som-sel + Aud/Som-nonsel) were clustered in the modules (Figures 7B, D, F, and I, Chi-Square test, χ² = 26.0, p = 9.3×10^-5), which was consistent with the reported difference in inputs to modules compared to matrix (A. M. Lesicko et al., 2016).

Figure 6:

Figure 7:

Additionally, while most of the acoustically-responsive cells in the matrix showed onset responses, the cells showing acoustically-driven offset responses were clustered in the modules (Figures 7C and J, Chi-Square test, χ² = 120.5, p <0.001). Moreover, a higher fraction of Aud/Som-nonsel cells in the modules had offset, not onset, auditory response than those located in the matrix (Chi-Square test, χ² = 13.2, p = 0.004, % of Aud/Som-nonsel cells = 57%, 15%, 8%, and 20%, for auditory offset cells in modules, auditory onset cells in modules, auditory offset cells in matrix, and auditory onset cells in the matrix, respectively), which suggests a multisensory integrative function by those cells as discussed later. Given the specific spatial location of the somatosensory-responsive cells in the modules which contain clusters of GABAergic cells, we asked if somatosensory stimulation could suppress the auditory response when both stimuli are presented simultaneously (bimodal stimulation). To test this hypothesis, we calculated and visualized the response index (RI), which was calculated as the ratio of the cellular response to bimodal stimulation to the cellular response to acoustic response alone for auditory-responsive cells. Accordingly, some cells showed suppressed (RI < 1) or enhanced (RI > 1) auditory responses in both modules and matrix following bimodal stimulation. We observed that matrix and modules have more cells with suppressed auditory responses than those having enhanced auditory responses. However, matrix had a higher fraction of cells with suppressed responses than modules, while modules had a higher fraction of cells with enhanced responses (Figures 7G and K, Chi-Square test, χ² = 12.9, p = 0.04).

Consistent with the above data, the LC_{(dorsal surface)} (Figures 8A-D, dorsolateral white rectangle) shared common response profiles with the LC_(microprism). For instance, both regions had a higher percentage of somatosensory responsive cells than the DC (Figures 8A-D, dorsal yellow and dorsomedial blue rectangles) as quantitatively shown (Figure 8H, one-way ANOVA, p = 0.02), which suggests that the LC is the main region of the IC shell that processes somatosensory information, consistent with the anatomical reports (Aitkin, Dickhaus, Schult, & Zimmermann, 1978; Jain & Shore, 2006; A. M. Lesicko et al., 2016; Liu et al., 2023). Similar to the LC_(microprism), the matrix and modules imaged from the dorsal surface had more percentage of Aud-sel cells than Som-sel and Aud/Som-nonsel cells (Figures 8E and I, Chi-Square test, χ² = 86.1, p <0.001). Collectively, the matrix had more auditory responsive cells than modules, while modules had more somatosensory responsive cells than matrix (Figure 8J, Chi-Square test, χ² = 66.5, p <0.001). Additionally, more cells with onset auditory response were found in the matrix, while the cells with offset auditory responses were clustered in the modules (Figures 8F and K, Chi-Square test, χ² = 500.2, p <0.001).

Figure 8:

Consistent with the LC_(microprism), some cells had suppressed or enhanced auditory responses following bimodal stimulation (Figure 8-figure supplement 1). The matrix imaged from the dorsal surface showed also more cells with suppressed auditory responses than those having enhancing auditory responses induced by bimodal stimulation. Additionally, the matrix had a higher fraction of cells with suppressed responses than modules, while modules had a higher fraction of cells with enhanced responses than matrix. In contrast to the LC_{(microprism),} the modules imaged directly from the dorsal surface had a higher fraction of cells with enhanced responses than those with suppressed responses, which may indicate that the ratio between suppressed and enhanced responses in modules could be altered based on the dorsoventral location (Figures 8G and L, Chi-Square test, χ² = 18.5, p = 1.6×10^-5). In general, these data confirmed the findings obtained from the LC imaged via microprism and emphasized the functional difference between the LC and DC.

LC and DC differ in their acoustic and multimodal response properties in awake animals

To determine if the presence of ketamine/xylazine anesthesia influenced the response properties detailed above, additional experiments were performed in awake animals (n = 6 total, 3 for the LC_(microprism) and 3 for the DC). Similar to what was observed under anesthesia, the cells of the LC_(microprism) did not show clear evidence of tonotopy (Figure 9B), which was indicated by a low R² of the regressions fit between the BTFs of the cells and their locations over different anatomical axes (Figure 9-figure supplement 1A). Similarly, the cells of the modules were tuned to relatively higher frequencies than those in matrix (Mann-Whitney Test, z = 3.8, p = 1.0 x 10^-4, median of best frequency = 10 vs 7.1 kHz for modules (148 cells from 3 animals) vs matrix (287 cells from 3 animals)). Additionally, the data obtained from imaging the DC under awake preparations showed a similar tonotopic map (Figure 9D) as that observed under anesthesia indicated by a high R² of the regression fit between the BTFs of the cells and their locations over the rostromedial and caudolateral axis (Figure 9-figure supplement 1B). Yet, a greater proportion of high-frequency cells were visually observed in the mediodorsal portion of the DC, which was indicated by a higher R² of the nonlinear quadratic fit than that of the linear fit, leading to the observation of a narrow low-frequency zone similar to the previous work under awake preparations (A. B. Wong & Borst, 2019). The LC_(microprism) had also more cells tuned to higher frequencies compared to the DC (Figure 10A, Mann-Whitney test, p=1.6×10^-6). We also examined the width of the receptive field (RF) of the DC vs LC neurons at each tested sound level under anesthetized and awake preparations using the binarized receptive field sum (RFS) method (Bowen, Winkowski, & Kanold, 2020). Both DC and LC of the anesthetized animals had higher RFS mean and median than those of awake animals given that ketamine was known to broaden the RF (Guo et al., 2012). In both preparations (anesthetized and awake), the DC had a higher RFS mean than that of the LC (Figure 10-figure supplement 1, Kruskal-Wallis ANOVA, p < 0.001), which could be consistent with the finding that the DC had a lower SMI than LC.

Figure 9:

Figure 9-figure supplement 1:

Figure 10:

Figure 10-figure supplement 1:

Similar to what was shown in the anesthetized animals, the fraction of Noise-sel cells was higher than that of Tone-sel cells in both LC_(micrprism) and DC imaged directly from the top in awake animals (Figure 10B and D). In contrast to the anesthetized animals, the fraction of Non-sel cells had different profiles. For instance, while the percentage of Non-sel cells was higher than Noise-sel and Tone-sel cells in the LC_(microprism) (Figure 10B, one way ANOVA, f_(2,6) = 41.6, p = 3.0×10^-4), they were significantly lower than that of Noise-sel cells and similar to that of Tone-sel cells in the DC (Figure 10D, one way ANOVA, f_(2,6) = 12.8, p = 0.006). Despite that modules and matrix had similar fractions of Noise-sel, Tone-sel, and Non-sel cells (Figure 10C, Mann-Whitney Test (Modules vs matrix), z = -0.46, -0.46, and -1.85, p = 0.64, 0.64, and 0.06 for Noise-sel, Tone-sel, and Non-sel, respectively), matrix had more spectrally (Figure 10E, two-way ANOVA: p = 5.3×10^-10, 1.0, and 0.17 for the region, sound level, and interaction, respectively) and temporally (Figure 10F, two-way ANOVA: p = 4.6×10^-7, 1.0, and 0.11 for the region, sound level, and interaction, respectively) modulated cells than modules, which was similar to anesthetized animals. Additionally, the LC_(microprism) and DC showed different profiles than those found by anesthetized animals. While LC and DC had similar SMI (Figure 10G, two-way ANOVA: p = 0.22, 2.1×10^-12, and 1.1×10^-6 for the region, sound level, and interaction, respectively), the LC_(microprism) had a lower average of TMIs compared to the DC (Figure 10H, two-way ANOVA: p = 1.9×10^-0, 6.6×10^-24, and 0.01 for the region, sound level, and interaction, respectively).

Similar to the anesthetized preparations, the auditory and somatosensory responsive cells were topographically organized across modules and matrix (Figures 11A, C, and E). However, it was noticed that the somatosensory and auditory responsive cells were clustered within a pool of nonresponsive cells (Figure 11E). Some of these nonresponsive cells to unmodulated broadband noise (Figure 11-supplementary figures 1A) were found to respond to pure tones (Figure 11-supplementary figures 1B, blue arrows) indicating that these nonresponsive cells may be sensitive to a specific feature or type of sensory stimulations. Additionally, some of the cells that did not respond to any of the tested stimulations showed spontaneous activity (Figure 11-supplementary figure 1C) indicating the viability of the nonresponsive cells. Within the pool of responsive cells, the matrix and modules imaged via microprism had higher fractions of Aud-sel cells than those of Som-sel and Aud/Som-nonsel cells (Figures 11E and G, Chi-Square test, χ² = 42.5, p = 4.4×10^-8), which was similar to the data obtained from the anaethetized animals. Additionally, matrix had a higher fraction of auditory responsive cells than modules, while modules had a higher fraction of somatosensory responsive cells than matrix (Figures 11A, D, and H, Chi-Square test, χ² = 30.2, p = 1.2×10^-6). In contrast to the anesthetized animals, the offset cells were very rare (only 6 cells/3 animals), and were evenly distributed between modules and matrix resulting in no significance difference between matrix and modules (Figures 11C and I, Chi-Square test, χ² = 2.8, p = 0.41). Similar to anesthetized data, we found that matrix and modules had higher fractions of cells with auditory suppressed responses than those with auditory enhanced responses induced by bimodal stimulation indicating an overall suppression in the cellular auditory response. In contrast to the anesthetized animals, modules had a higher fraction of cells with suppressed response than matrix and a lower fraction of cells with enhanced responses compared to matrix (Figures 11G and J, χ² test = 31.9, p = 5.2×10^-7) indicating the possibility that the bimodal interactions in modules could be affected by the arousal state of the animal.

Figure 11:

Figure 11-figure supplement 1:

Discussion

Summary

As a new imaging approach, a microprism was used to image the LC surface, to provide the first-ever sagittal image of the LC in vivo and to examine its functional organization, which revealed the distinction between strategies used by LC and DC for processing acoustic and non-acoustic information. Imaging the LC either via microprism or from the dorsal surface revealed that it contains GABAergic modules that were embedded in a non-GABAergic matrix.

Although modules and matrix were found to be activated by somatosensory and auditory stimuli, they were functionally distinct. The matrix had more auditory-responsive cells with a lower threshold for complex sound detection and had more cells responsive to complex spectral and temporal auditory features than modules. In contrast, modules had more somatosensory responsive cells, consistent with anatomical differences in input to the modules compared to matrix (A. M. Lesicko et al., 2016). Somatosensory stimulation suppressed auditory responses in the matrix, suggesting that multisensory interactions exist between the matrix and modules.

The novelty of the imaging setup and technical considerations

Microprism-based 2P-microscopy has been successfully used to image the cortical layers in the somatosensory and visual cortices of mice (Andermann et al., 2013) and to image other inaccessible cortical surfaces like those folded in the brain sulci such as the medial prefrontal cortex (Low, Gu, & Tank, 2014) or entorhinal cortex (Beckmann et al., 2019; Heys, Rangarajan, & Dombeck, 2014; Low et al., 2014). In those studies, the glass window glued to the microprism was used, which requires significant skill and potential removal of brain regions to allow prism access (Heys et al., 2014). Another approach is to use a micro metal or diamond-made knife to create a pocket with minimal tissue damage (Andermann et al., 2013; Marks & Goard, 2021). Accordingly, we combined and modified some of these procedures to improve the microprism implantation to image the sagittal surface of the LC, which is deeply embedded into the brain tissue in a highly vascular area, with minimum damage to the IC or the surrounding vasculature.

Previous work using microprisms has been done chronically for up to 60 days (Heys et al., 2014; Low et al., 2014; Wenzel, Hamm, Peterka, & Yuste, 2017), with minimal tissue inflammation (Andermann et al., 2013). The current study was initially done on anesthetized animals to examine and validate the new approach, which later revealed similar organizational differences between LC and DC in awake animals. Although imaging through a microprism and 1-mm cover glass appears to have minimal impact on image quality compared to direct 2P imaging (Chia & Levene, 2009), it was important to compare LC imaging via microprism with the direct 2P imaging of the dorsal surface to validate the new imaging method. We found no significant differences between LC response properties when imaging via microprism or imaging directly from the dorsal surface. However, we did observe significant functional differences between the LC and the DC.

Responsiveness to spectral and temporal features

While the LC showed inconsistent and fragmented tonotopy, the DC showed a consistent one, as previously reported (Barnstedt et al., 2015; A. B. Wong & Borst, 2019). These data were supported by a lower R² of the best regression fit between BTFs of the LC cells and their locations compared to a higher R² obtained for the DC. The direction of the best regression fit shown by the LC was inconsistent across animals, which could be related to the fact that the microprism did not permit a full view of the LC surface. Given that microprism permits imaging of only a single area per animal, it is expected that each animal may have a different profile of cellular organization based on the location of the imaged area on the LC.

The surface of the DC was biased to lower frequencies compared to LC, consistent with previous imaging and recording studies (Ito et al., 2014; Romand & Ehret, 1990; Shen et al., 2003; Stiebler & Ehret, 1985). Despite their response to pure tones, the DC and LC, particularly in the matrix, were more responsive to AM-noise, which fits their function as nonlemniscal divisions (Chen, Cheng, Ito, & Song, 2018; Faye-Lund & Osen, 1985; Loftus et al., 2008).

The spectral and temporal features of animal vocalization and speech are not independent (Singh & Theunissen, 2003). We observed that LC neurons, particularly in the matrix, had greater responsiveness to spectral and temporal modulation than the DC. Given the role of inhibition in spectral and temporal processing (Cai et al., 2018; Chang, Bao, Imaizumi, Schreiner, & Merzenich, 2005; Wehr & Zador, 2003), GABAergic cells of the modules could play an important role in increasing the spectral and temporal modulation indices, especially at lower sound levels, which could be critical for localizing objects in noisy environments given the importance of spectral and temporal cues for sound localization (Davis, Ramachandran, & May, 2003; Goodman, Benichoux, & Brette, 2013; Sinex, 2005). Therefore, future work is critical to image the functional responses of the GABAergic cells of the modules and compare them with the non-GABAergic cells inside and outside the modules.

Somatosensory vs. auditory responses

Consistent with the previous anatomical findings (A. M. Lesicko et al., 2016), modules showed more somatosensory responsive cells than matrix. Since modules receive somatosensory inputs from different brain regions (A. M. Lesicko et al., 2016; Lohse et al., 2022; Zhou & Shore, 2006), the somatosensory response shown by modules could be modulated by different types of somatosensory stimulations, which is an open area for future studies. In contrast to modules, we observed that most auditory-responsive neurons, particularly onset-responsive neurons, were found in the matrix, which receives mostly auditory inputs (A. M. Lesicko et al., 2016). However, a subpopulation of acoustically-driven offset-responsive neurons preferably resides in the modules compared to the matrix. Given that the offset response could be generated as a rebound to the termination of the ongoing inhibition during the sound presentation (Xu, 2014), GABAergic cells in the modules could play an important role in the coding of sound termination as found before in the cortex (Solyga & Barkat, 2021). Because some of the offset acoustically-driven cells inside the modules were also activated by the somatosensory stimulation, this fraction of cells could then encode the termination of somatosensory-associated sounds through the local inhibitory circuits of the modules. Under in-vivo recording of the mouse IC, the offset responses were previously reported from cells that were non-GABAergic and located in the caudal and shallower parts of the IC (Kasai, Ono, & Ohmori, 2012). These data suggest that the module/matrix system permits the preservation of distinct multimodal response properties, which is consistent with their anatomical distinction, despite the massive integration of inputs in the LC.

Consistent with previous reports (Aitkin et al., 1978; Aitkin, Kenyon, & Philpott, 1981; Jain & Shore, 2006), we observed that combined auditory-somatosensory stimulation generally suppressed neural responses to auditory stimuli and that this suppression was the most prominent in the LC matrix of anesthetized animals. Mechanistically, it is not known if such suppressive responses are created in the LC or are inherited from other brain regions since auditory-somatosensory interactions occur as early as in the cochlear nucleus, and are also found in the IC, thalamus, and cortex (Aitkin et al., 1981; Haenggeli, Pongstaporn, Doucet, & Ryugo, 2005; Kimura, Imbe, Donishi, & Tamai, 2007; Li & Mizuno, 1997; Lohse, Dahmen, Bajo, & King, 2021; Shore, 2005; Shore, Vass, Wys, & Altschuler, 2000; C. Wong, Chabot, Kok, & Lomber, 2015), reviewed in (Brunelle & Llano, 2023; Lohse et al., 2022). We note that a systematic exploration of the stimulus space was not conducted here such that other acoustic or somatosensory stimuli, different relative amplitudes, or timings could have produced enhancement rather than suppression. Future work will be needed to further explore this stimulus space to more fully understand multisensory processing in the LC.

Effect of ketamine anesthesia

The data obtained from anesthetized animals were similar to those obtained from awake animals in the findings related to the tonotopy, the preference of the cells to respond to AM-noise in the LC and DC, the increased sensitivity of matrix to detect the spectral and temporal features of sound, and the bimodal response properties of the matrix vs modules. However, anesthetized and awake animals showed some differences that could be related to the effect of ketamine. For instance, ketamine was found to broaden the RF of the auditory-responsive neurons (Guo et al., 2012). Consistently, the width of the RF of all DC and LC neurons was higher under ketamine anesthesia compared to awake preparation. However, the DC neurons had a wider RF than those in the LC. The relatively narrower RF of the LC could be consistent with their higher SMI (discussed below) compared to that of the DC. Both findings could be related to the degree of inhibition that may shape the inhibitory bands of the RF in the LC neurons given that the LC has prominent GABAergic modules. Also, while the average of the SMIs was negative with no difference between the LC_(microprism) and the DC of awake animals, it was positive with a higher value in the LC_(microprism) compared to the DC of the anesthetized animals indicating the preference of cells to respond more to pure tones under anesthesia. Given that ketamine was reported to increase the duration of the tone-evoked response in the auditory cortex (Guo et al., 2012), the same effect could increase the AUC for the tone-evoked response that might result in a positive average of SMI. Additionally, while LC had a higher average of TMIs than that of the DC of anesthetized animals, the LC had a lower mean of TMIs than that of the DC of awake animals indicating that ketamine may specifically modulate the response properties of the LC via its effects on GABAergic cells (Behrens et al., 2007; Deane et al., 2020; Weckmann et al., 2019), given that GABAergic inhibition reportedly contribute to the processing of the temporally modulated signals (Alluri et al., 2021; Burger & Pollak, 1998; Cai et al., 2018; Gourévitch, Mahrt, Bakay, Elde, & Portfors, 2020). In both anesthetized and awake animals, matrix neurons had always a consistent ratio between the cells with suppressed (∼ 70%) and enhanced (∼30%) auditory responses induced by bimodal stimulation. However, modules had variable ratios across anesthetized and awake animals. While modules had a lower fraction of cells with enhanced responses than matrix in awake animals, they had a higher fraction of cells with enhanced responses in anesthetized animals, which could be consistent with the direct inhibitory effect of ketamine on GABAergic cells (Behrens et al., 2007; Deane et al., 2020; Weckmann et al., 2019). Given that modules are clusters of GABAerigic cells and terminals, modules could be specifically modulated by ketamine anesthesia.

Conclusions

Here, using a novel microprism approach to image the LC in anesthetized and awake mice, we observed that 1) DC displayed tonotopic organization which was not consistently present in the LC, 2) DC neurons were generally more sound-responsive than LC neurons but that neurons in both regions had preferential responses to spectrally- and temporally-complex stimuli, 3) that matrix regions of the LC were more responsive to spectrally- and temporally-complex stimuli than modules and that modules were more responsive to somatosensory stimuli and 4) that bimodal stimulation led to suppressive responses that were more prominent in the matrix. In general, these data highlight significant organizational differences between the LC and DC and between modules and matrix in the LC. Future work using microprisms in the LC will help to reveal the potential differential role of matrix vs. module neurons in multisensory integration tasks in behaving animals.

Materials and methods

Animal subjects

Male mice of 8-12 weeks of age were used. Males were used as our preliminary data suggest male/female differences in DC organization (Ibrahim, Shinagawa, Xiao, Asilador, & Llano, 2022) and will be published separately. GAD67-GFP knock-in mice (developed and shared with permission from Dr. Yuchio Yanagawa at Gunma University and obtained from Dr. Douglas Oliver at the University of Connecticut), where GFP is exclusively expressed in GABAergic cells (Tamamaki et al., 2003), were used to visualize the GABAergic cells and LC modules. To simultaneously monitor calcium signals and visualize the GFP+ cells, Tg(Thy1-jRGECO1a)GP8.20Dkim/J mice (Jackson Laboratory, Stock# 030525) were crossed with GAD67-GFP knock-in mice to generate GAD67-GFPxThy1-jRGECO1a hybrid mice. The mice were housed and bred at the animal research facility of The Beckman Institute for Advanced Science and Technology. The animal care facilities are approved by the Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC). To phenotype the GFP-positive animals, a single-photon fluorescence microscope was used for the transcranial examination of the pups at postnatal day 4 using excitation: 472/30 nm and emission 520/35 nm filters along with dichroic 505 nm long pass at low power magnification (2.5x/0.08, Olympus objective, MPlanFL N, Japan). The positive pups exhibited a green fluorescence in the cerebellum, the cerebral cortex, and the olfactory bulb and were only kept with their parents until weaning. The genotyping for jRGECO1a was done using a probe composed of a forward primer: GCCGCCGAGGTCAAGA & a reverse primer: TCCAACTTGATGTCGACGATGTAG) by Transnetyx™ (Transnetyx, USA). Samples of the mice’s tail snips around one month of age were used for genotyping. All applicable guidelines for the care and use of animals were followed. All procedures were approved by the Institutional Animal Care and Use Committee (IACUC).

Surgery

The general procedures of craniotomy were previously described (Goldey et al., 2014), with modifications to place craniotomy over the IC. Before surgery, mice were anesthetized with a mixture of ketamine, xylazine, and acepromazine (100 mg/kg, 3 mg/kg, and 3mg/kg, respectively) delivered intraperitoneally. The anesthesia was maintained during the surgery and imaging using only ketamine (100 mg/kg). To prevent neural edema during or after the craniotomy, an intramuscular injection of dexamethasone sodium (4.8 mg/kg) was given just before the surgery using an insulin syringe. After placing the animal in the stereotaxic apparatus (David Kopf Instruments, USA), both eyes were protected by applying Optixcare lubricant eye gel (Aventix Animal Health, Canada). The hair on the scalp was then removed by massaging the scalp with a depilatory cream (Nair) using a cotton-tipped applicator and leaving the cream on the scalp for 4-5 minutes. The cream was then removed by a thin plastic sheet (flexible ruler) to leave a hair-free area on the scalp. The remaining tiny hairs were then removed by alcohol swab and the area was then sterilized by applying 10% Povidone Iodine (Dynarex, USA) using a sterile cotton-tipped applicator. The medial incision was made with a scalpel blade #10, and 0.2 ml of 0.5% lidocaine was injected intradermally into the scalp. The skin extending from the medial line to the temporalis muscle was completely removed using a pair of microscissors to produce a wide skinless area above the skull. Number 5/45 forceps were used to remove any remaining periosteum. The remaining dried or firmly attached pieces of periosteum were removed with a scalpel blade #10. The skull was cleaned with sterile saline and dried with gently pressurized air. Using the stereotaxic apparatus, a wide area of ∼3 x 4 mm above the left IC was made. A micro drill bit (size #80, Grainger, USA) was used to drill through the skull starting from the rostrolateral (Figure 1A, yellow circle) region to lambda (Figure 1A, black circle) following the border in a clockwise direction (Figure 1A, dotted black rectangle). To prevent overheating of the superficial brain tissues and to mitigate the occasional spurts of skull bleeding during the drilling, ice-cold sterile saline was used to intermittently irrigate the surface. A stream of pressurized air was also applied during the drilling procedures to prevent overheating and remove the debris produced by the drilling. Caution was taken not to pierce the dura when performing the craniotomy while crossing the sagittal or the lambdoid sutures to avoid damaging the underlying sinuses. After drilling, the skull was irrigated in sterile saline and the bone flap (the undrilled bone over the craniotomy area) was gently examined for complete separation from the rest of the skull (Figure 1B). Using a pair of no. 5/45 forceps, the bone flap was gently removed. To control the bleeding if it occurred, a piece of sterile hemostatic gel (Adsorbable Gelatin Sponge USP, Haemosponge, GBI, India) which was pre-soaked in ice-cold saline, was applied to the bleeding spot. Once the bleeding ceased, the brain was kept covered in sterile saline. In some surgeries, the dura was peeled off from the surface of the IC while removing the bone flap. In the surgeries where the dura remained intact, a Bonn microprobe (F.S.T, Germany, Item # 10032-13) was used to pierce the dura in an area that is not above the IC and devoid of cortical vasculature (e.g., a part of exposed cerebellum). After piercing, the dura was carefully and gently lifted, and a pair of no. 5/45 forceps were used to grab the dura to gently tear it to the extent of the transverse sinus to avoid bleeding of the major venous structures.

Using the stereotaxic apparatus, a 30° diamond knife (F.S.T, Germany, Item # 10100-30) was placed in its vertical position. From the medial line, the knife was placed at the lateral horizon of the IC. The knife was then inserted 1.5 mm ventrally. Caution was taken to avoid hitting any major blood vessels to avoid bleeding. From its position, the knife was retracted 0.3-0.5 mm medially to create an initial pocket for the insertion of the microprism (Figures 1C-E). The microprism (1.5 mm-silver coated, Tower Optical Corporation) was then inserted, so one of its sides was parallel to the lateral surface of the diamond blade, the other side was facing up, and its hypotenuse was facing ventrally (Figures 1F and G). The wooden piece of a cotton swab was carefully driven down to push the microprism until its top surface was at the same level as the IC dorsal surface. Keeping the microprism in its place with the wooden end of a swab, the diamond knife was carefully retracted upward leaving the microprism behind (Figure 1H). If the retraction of the diamond knife pulled the prism up, the wooden piece would be used to push it down simultaneously with the knife movement. To control the bleeding, a piece of sterile hemostatic gel was used as before. The tip of the diamond knife was used to secure the microprism in its position for at least 10 minutes (Figure J). During that time, the solution of the 1% agarose was prepared for use. A pair of forceps were then used to gently place a sterile 5 mm cover glass #1 (Thomas Scientific, USA) over the skull covering the opened area above the IC. The cover glass was secured by a wooden trimmed piece of sterile cotton swab by gently pressing the cover glass from the top. Since the surface of the IC was located below the surfaces of both the cerebellum and cerebral cortex, the cover glass usually traps about a 0.7-1.0 mm thick sheet of sterile saline. The movement of the liquid medium caused by heart pulsation was found to be a source of motion artifact during imaging, so the gap between the surface of the IC and the cover glass was filled with 1% agarose gel that was made in saline at a temperature range of 33-35°C (Figure 1K). The extra agarose gel pieces were cut off and removed by the scalpel blade and the wet skull was then dried out using pressurized air. A titanium headpost as described before (Goldey et al., 2014) was glued carefully on the top of the skull to be at the same level as the cover glass. Following the manufacturer’s instructions, the C&B Metabond (Parkell, Japan) was used to secure the headpost in its place. Excluding the steps of the microprism insertion, the same procedures were followed for DC imaging from the surface (Figure 1L). For awake preparations, the same procedures were followed under isoflurane anesthesia starting with 4% of isoflurane as an induction dose with 1-2% isoflurane during the surgery. The animal was transferred under the microscope objective and left to recover from isoflurane anesthesia for 1 hour before imaging. To ensure the animal’s recovery, the animal was visually inspected for its breathing rate (Ewald, Werb, & Egeblad, 2011) and ongoing limb movements.

Acoustic stimulation

Using a custom-made MATLAB (The MathWorks, Natick, MA, USA) code, either 500 ms pure tone or AM-noise was generated. Thirty-five (5×7) pure tones with a series of sound pressure levels (40, 50, 60, 70, 80 dB SPL) and a series of carrier frequencies (5000-40000 Hz with a half-octave gap) were presented with a cosine window. Forty-five (5×9) 100% AM white noise bursts with a series of sound pressure levels (40, 50, 60, 70, 80 dB SPL) and modulation frequencies (0, 2, 4, 8, 16, 32, 64, 128, 256 Hz) were generated. The stimuli of either pure tone or AM-noise combinations were played in random sequence to the mice with a 600 ms interstimulus interval (ISI) by a TDT RP2.1 processor (Tucker-Davis Technologies, US) and delivered by a TDT ES1 speaker (Tucker-Davis Technologies, US). In other sets of experiments, 500 ms of unmodulated broadband noise was presented 5 times with ISI of 600 ms at 80 dB SPL to examine the acoustic activity between modules and matrix compared to somatosensory stimulation.

The output of the TDT ES1 speaker was calibrated using a PCB 377A06 microphone, which feeds a SigCal tool to generate a calibration file for all tested frequencies (5–40 kHz). To enable the custom-made MATLAB code to read this calibration file, the values were first processed by MATLAB signal processing toolbox (sptool) to generate a 256-tap FIR filter to apply the calibration using the following parameters [arbitrary magnitudes, least square, order: 256, sampling rate: 97,656.25, frequency vector (5–40 kHz), amplitude vector (40–80 dB SPL), and weight vector ones (1128)].

Somatosensory stimulation

The somatosensory stimulation was done by deflecting the right whiskers of the animals at an average excursion of 1.2 mm with a 50 Hz rate over 500 ms to match the acoustic stimulus. A total of five stimuli were presented with 600 ms ISI. The whisker deflection was performed using a Brüel & Kjær vibrator (V203 10/32 UNF-CE, United Kingdom) that was driven by a wave generator (4063, BK precision, USA) after amplifying the signals with an SLA1 100 W power amplifier (ART, USA).

Two-Photon imaging

Immediately after surgery, the anesthetized animal was taken and secured under the microscope objective by clamping the arms of the head post to two perpendicular metal posts mounted on the microscope stage. A custom-built 2P microscope was used. The optical and the controlling components were supplied from Bruker, Olympus, and Thorlabs. The imaging of the DC was made using a 20x water-immersion objective (LUMPlanFI/IR, 20X, NA: 0.95, WD: 2 mm; Olympus Corporation, Tokyo, Japan), while the lateral surface of the IC was imaged by the long working distance 16x water-immersion objective (N16XLWD-PF - CFI LWD Plan Fluorite Objective, 16X, NA: 0.8, WD: 3.0 mm; Nikon, Tokyo, Japan) to be able to reach the focal point of the microprism at its hypotenuse. For imaging both the GFP or jRGECO1a signals, the excitation light was generated by InSight X3 laser (Spectra-Physics Lasers, Mountain View, CA, USA) tuned to a wavelength of 920 or 1040 nm, respectively. A layer of a 1:1 mixture of wavelength multipurpose ultrasound gel (National Therapy, Canada) with double-deionized water was used to immerse the objective. This gel was able to trap the water and reduce its evaporation during imaging. The emitted signals were detected by a photomultiplier tube (Hamamatsu H7422PA-4, Japan) following a t565lp dichroic and a Chroma barrier et525/70m filter for GFP and et595/50m filter for jRGECO1a signals. Images (512×512 pixels) were collected at a frame rate of 29.9 Hz in the resonant galvo mode. Data were collected from the dorsal surface of the IC by scanning the surface of the IC based on the GFP and jRGECO1a signals through the medial and lateral horizons of the IC. Generally, the scanning was started by moving the 20x objective to the most ventromedial position, where there is no ability to see any cells (Figure 2B, medial horizon). Then, the objective was driven to slowly move dorsolaterally to get cells in focus (Figure 1M), where the average depth of the first field of view was 169 μm. The lateral movement of the objective was intended to be aligned and parallel to the curvature of the DC surface indicated by getting the superficial cells in focus to avoid any possibility of imaging the upper layers of the central nucleus of the IC, so the average depth of the most superficial point of the imaging was 92 μm. The scanning was ended by moving the objective to the most ventrolateral position of the IC, where there is no ability to see GFP laterally (Figure 2B, lateral horizon), where the last field of view was collected at an average depth of 225 μm. Generally, each field of view was selected based on the expression of GFP signals as an indicator for the cells in focus and being acoustically active by using a search stimulus that was 500 ms broadband noise with zero modulation at 80 dB SPL. The z-axis of the objective was determined relative to the pia surface. The relative position of each region was tracked through the micromanipulator (MP-285, Sutter Instruments) that controlled the microscope objective. The positions of imaged areas were further aligned across animals to a common coordinate using the midline, lateral extremes, and the vasculature landmarks taken via imaging of the GFP signals of the cranial window with low magnification (4x Olympus objective, 4×4 binning, 50 ms frame time) using a Lumen 200 bulb (Prior Scientific Inc, USA) and CoolSNAP MYO camera (Photometrics, USA, excitation: 488 nm and emission: 515–550 nm). Imaging of the lateral surface of the IC via microprism was done similarly. In both cases, the frame timing of the scanner and the sound stimuli were both digitized and time-locked using a Digidata 1440A (Molecular Devices, Sunnyvale, CA, USA) with Clampex v. 10.3 (Molecular Devices, Sunnyvale, CA, USA).

Laser lesion

The laser lesion was produced using a 720 or 820 nm laser beam at 74-80 mW power at the level of the objective for 100 to 250 seconds depending on the size of the lesioned area, the speed of laser scanning (frame rate), and the depth of the lesion (z-axis length). In brief, the area of interest and the dimension of the laser lesion were initially determined using 920 nm as indicated by the GFP signals through scanning one of the GABAergic modules. The depth of the GABAergic modules was then determined by counting the optical slices starting from their superficial surface. The tunable laser was then switched from 920 nm to the shorter wavelength (720 or 820 nm) to scan the same area and depth at a rate (1-2 frames/ second).

Data processing

Data collection

The data were collected as separate movies (512×512 pixels) for each pure tone or AM-noise runs in a resonant galvo mode. Depending on the amplitude and frequency combinations for each type of acoustic stimulus, 40- or 52-second periods were assigned as a movie’s length for pure tone (35 stimulus combinations), or AM-noise (45 stimulus combinations), respectively. Using ImageJ software, the z-projection was used to compute one single image representing either the sum, the standard deviation, or the median of all the image sequences in the movie. Based on these single images, the region of interest (ROI) was manually drawn around each detectable cell body. Figures 3A and B represent the sum of the image sequences of the movie collected from the dorsal surface, or the lateral surface of the IC via microprism, respectively. For the subsequent processing steps, Python open-source libraries were used (listed below).

Motion Correction and Filtering

The imread function from the OpenCV library was used in grayscale mode to import images into a numpy array from a single folder containing TIFF images. The array was then exported as a single TIFF stack into a temporary folder using the mimwrite function from the imageio library, and the process was repeated for each folder. The NoRMCorre algorithm (Pnevmatikakis & Giovannucci, 2017) embedded in the CaImAn library (Giovannucci et al., 2019) was used to apply motion correction to each of the TIFF files. The data were then imported into a numpy array, rounded, and converted to 16-bit integers. The images were filtered using a 2D Gaussian filter with a sigma value of 1 (surface View) / 2 (prism View) in each direction, then a 1D Gaussian temporal filter with a sigma value of 2 was applied using the ndimage.gaussian_filter and ndimage.gaussian_filter1d function from the scipy library, respectively.

Data Extraction

The ROI sets, which were manually created using ImageJ, were imported using the read_roi_zip function from the read_roi library. The sets were then used to create two masks; one mask was used as a replica of the ROIs and the second mask was made around the original ROI (roughly four times larger in area). The smaller mask was applied to find the average pixel value within each ROI, while the larger mask was applied to find the average pixel value of the neuropil. The neuropil correction was applied using the following equation (Akerboom et al., 2012; Kerlin, Andermann, Berezovskii, & Reid, 2010);

To identify the calcium signals, Δf/f was calculated by using the following equation.

where the background value is the slope estimating the background levels with fluorescence bleaching factored in. The data were then reorganized so that all segments with the same stimulus frequency and stimulus amplitude were grouped. The area under the curve (AUC) of the calcium signals of the excitatory responses only was then used as a metric to determine the magnitude of the response and for subsequent analysis (Berens et al., 2018; Overk, Rockenstein, Florio, Cheng, & Masliah, 2015).

Cell Flagging

The correlation coefficient between each of the trials was calculated using the stats. Pearson function from the scipy library. The average correlation coefficient was calculated for each stimulus frequency and stimulus amplitude. Similar to previous work (A. B. Wong & Borst, 2019), if the average correlation coefficient was above the threshold of 0.6, the cell was flagged as being responsive to that stimulus combination of frequency and amplitude, and the best tone and modulation frequencies were calculated for every cell (Barnstedt et al., 2015). Knowing BTF and the best modulation frequency (BMF) for each neuron enabled us to calculate the SMI and TMI at each sound level following these equations;

The SMI has values ranging from (-1 to 1). While neurons with values closer to 1 are more responsive to pure tone, neurons with values closer to -1 are more responsive to noise with no modulation.

Similarly, the TMI has values ranging from (-1 to 1). While neurons with values closer to 1 are more responsive to AM-noise, neurons with values closer to -1 are more responsive to noise with no modulation.

For the comparison between somatosensory and auditory responses, the threshold of 0.4 was used as a cutoff threshold to determine the responsive cells to each stimulation. The average response of each somatosensory, acoustic (500 ms of unmodulated noise), or simultaneous somatosensory and acoustic (bimodal) stimulation was then calculated for each responsive cell. The correlation method was validated as a tool to determine the responsive cells to somatosensory stimulation by generating an ROC curve comparing the automated method to a blinded human interpretation. The AUC of the ROC curve was 0.88. This high AUC value indicates that the correlation method can rank the random responsive cells than the random nonresponsive cells. At the correlation coefficient (0.4), which was the cutoff value to determine the responsive cells for somatosensory stimulation, the specificity was 87% and the sensitivity 72%, the positive predictive value was 73%, and the negative predictive value was 86%. Additionally, all the false-responsive cells were excluded from the analysis. To examine the effect of somatosensory stimulation on auditory responses, the response index (RI) was calculated for the auditory-responsive cells using the following equation;

Reurons with RI values > 1 have enhanced or higher responses induced by simultaneous auditory and somatosensory stimulation, while neurons with RI values < 1 have downregulated or suppressed auditory responses induced by auditory and simultaneous somatosensory stimulation.

Map Generation

The average radius for all ROIs was calculated to ensure that all cells on the tonotopic map had uniform radii. A color key was also generated for every corresponding map. A blank canvas was generated using the Image module from the pillow library, and a circle was created for each cell and filled with a corresponding color to its value.

Receptive field sum (RFS) analysis

The binary RFS was calculated for every cell by determining whether the cell was responsive to a particular combination of frequency and amplitude (Bowen et al., 2020). A value of 1 was given to the responsive cell, if it responds to a single frequency/amplitude combination. In contrast, a value of 1 was added to the tally for each responsive combination if the cell was broadly responsive. Accordingly, the higher and lower values of RFS indicate a broad and narrow RF, respectively. To examine the differences between the RFS of the cells at different areas under different preparations, the cumulative distribution function was used.

Assessment of tonotopic organization

The tonotopy was assessed by examining the best fit of linear or nonlinear quadratic polynomial regressions between the BTFs of cells and their locations along different anatomical axes. For the LC_(microprism), the regression fit was examined along the main four axes such as dorsal to ventral, rostral to caudal, dorsocaudal to ventrorostral, and dorsorostral to ventrocoudal. Based on the previous reports showing the tonotopic organization over the DC (Barnstedt et al., 2015; A. B. Wong & Borst, 2019), the tonotopy was examined along the rostromedial to caudolateral axis. The x₂ and y₂ coordinates of each cell were determined. At the starting point of each axis, the x₁ and y₁ coordinates were assigned values as (0, 0). Then the distance (d) between each cell and the starting point at each axis was calculated based on the following equation;

For each 50 um (bin size), the geometric average of the BTFs of all cells as well as their distances from the starting point at each axis were computed. The linear or quadratic polynomial regressions were then examined between the geometric averages of BTFs and the geometric averages of distances of cells for each animal.

Detecting the spontaneous activity from the nonresponsive cells

Movies have been recorded from the nonresponsive cells of the LC_(microprism) when no stimulation was presented. The time trace of jRGECO1a signals was extracted from each cell. A polynomial fitting of the time trace was used to correct for photobleaching and to adjust for the baseline for the whole trace. The sliding average of the baseline was then calculated for 1 second (∼29 frames), and the z score was calculated for the subsequent frames. The spontaneous activity was detected when the z-score was > 3.0.

Statistics

All statistical analyses and graphs were made using Origin Pro software. The normality of the distributions of the data was examined using a Kolmogorov-Smirnov test. The Mann-Whitney test as well as two-sample or paired t-test were used to examine the difference in data distribution or the difference between samples of non-normal and normal distribution, respectively. Either one-way or two-way ANOVA as well as Fisher as a post hoc test were used for the mean comparison between different groups of normal distribution. The Chi-square (χ²) test or Kruskal-Wallis ANOVA tests were used to compare the samples based on cell numbers of their cellular population. The cumulative frequency function was computed to calculate the fraction of cells tuned at each frequency based on their best frequency. Differences were deemed significant when p-value <0.05. The upper and lower borders of each box plot are at the 25^th and 75^th percentile, with whiskers at 95% spread.

Histology

After the 2P imaging, the mouse was transcardially perfused with 4% paraformaldehyde (PFA) solution in PBS solution. The isolated brain was then placed overnight in PFA for a complete fixation. The brain was then moved from the PFA to a graded sucrose solution (10, 20, and 30%) until it sunk. 50 μm coronal sections were taken at the level of the IC using a Leica cryostat. The sections were then mounted on a slide, coverslipped, and imaged using a confocal microscope for imaging the GFP signals as an indicator for the GABAergic cells in the GAD67-GFP knock-in mouse (excitation: 488 nm and emissions: 515–550 nm).

For the examination of jrgeco1a expression in the IC, the mouse was transcardially perfused with saline with a high concentration of potassium and a high ratio of calcium to magnesium [(in mM) 3.5 KCl, 25 glucose, 125 NaCl, 1.2 KH2PO4, 2 CaCl2, 0.5 MgSO4] under ketamine anesthesia. The brain was then fixed by a subsequent perfusion of 4% PFA. Selected 40 μm coronal sections at the level of the IC were immunostained for NeuN using (NeuN (D4G40) XP® Rabbit mAB, Cell Signaling, Cat# 62994, 1:100) primary antibody and (goat anti-rabbit IgG conjugated to Alexa 405, 1:200) secondary antibody following the immunostaining procedures that were previously described (Issa, Sekaran, & Llano, 2023; Vaithiyalingam Chandra Sekaran et al., 2021). Confocal images were obtained using a Confocal Zeiss LSM 710 Microscope using 10 and 40x objectives and tile scan functions. For NeuN, Alexa 405 signals were obtained by a 405 nm laser and emission at 415-466 nm. The GFP signals were obtained by a 488 nm laser and emission at 501-551 nm. The jRGECO1a signals were obtained by 561 nm laser and emission at 600-700 nm. All cells with somas expressing a red fluorescence of jrgeco1a were quantified in 2P-imaged areas of the DC and LC. The density of the cells expressing jrgeco1a was then calculated based on the cells expressing NeuN.

Artwork

The cartoon of the mouse brain in Figures 1 C, F, and I was obtained from Brain Explorer 3D viewer Allen Mouse Brain (https://connectivity.brain-map.org/3d-viewer?v=1) and modified to be suitable for illustration. All figures were designed and made using Adobe Illustrator (Adobe, San Jose, CA). To keep working within the Adobe environment to avoid losing the resolution of the figures, Adobe Photoshop (Adobe, San Jose, CA) was used to crop the borders of some images to save space. The color pallets were chosen to avoid the combination of the green and red colors for all maps, so they will be suitable for color-blind readers. To induce the same effect in the microscope images combining the red (jRGECO1a) and green (GFP) fluorescence, we kept the green signals of the GFP unchanged, but we substituted the red fluorescence of jRGECO1a with magenta as a pseudocolor.

Acknowledgements

We thank Dr. Kush Paul (Bruker) for initial training on the craniotomy surgery, Dr. Jason Maclean (The Department of Neurobiology, University of Chicago) for assistance with 2P imaging, Dr. Michael Jacobs Goard (Neuroscience Research Institute, The University of California-Santa Barbara) for his guidance in the use of microprisms, Dr. Taher Saif (Department of Mechanical Science and Engineering, The University of Illinois Urbana-Champaign) for measuring the deflecting distance of the somatosensory stimulating probe using a high-speed camera, Dr. Pengfei Song (Department of Electric and Computer Engineering, The University of Illinois Urbana-Champaign) for his assistance with the somatosensory stimulating probe, and Dr. Nathiya Vaithiyalingam Chandra Sekaran (Department of Molecular and Intergartive physiology, The University of Illinois Urbana-Champaign) for her assistance in immunohistochemistry and confocal imaging. We also thank the NIDCD for funding this work under grants # R01DC013073 and R01DC016599.

Competing interests

On behalf of all authors, there are no financial or non-financial competing interests to be declared.