Abstract
Signaling between neurons and glia is necessary for the formation of functional neural circuits. A role for microglia in the maturation of connections in the medial nucleus of the trapezoid body (MNTB) was previously demonstrated by postnatal microglial elimination using a colony stimulating factor 1 receptor (CSF1R). Defective pruning of calyces of Held and significant reduction of the mature astrocyte marker glial fibrillary acidic protein (GFAP) were observed after hearing onset. Here, we investigated the time course required for microglia to populate the mouse MNTB after cessation of CSF1R inhibitor treatment. We then examined whether defects seen after microglial depletion were rectified by microglial repopulation. We found that microglia returned to control levels at four weeks of age (18 d postcessation of treatment). Calyceal innervation of MNTB neurons was comparable to control levels at four weeks and GFAP expression recovered by seven weeks. We further investigated the effects of microglia elimination and repopulation on auditory function using auditory brainstem recordings (ABRs). Temporary microglial depletion significantly elevated auditory thresholds in response to 4. 8, and 12 kHz at four weeks. Treatment significantly affected latencies, interpeak latencies, and amplitudes of all the ABR peaks in response to many of the frequencies tested. These effects largely recovered by seven weeks. These findings highlight the functions of microglia in the formation of auditory neural circuits early in development. Further, the results suggest that microglia retain their developmental functions beyond the period of circuit refinement.
Significance Statement
Auditory brainstem pathways are optimized for their special functions that are shaped during development, which relies on the functions of non-neuronal cells, such as microglia and astrocytes. When microglia were pharmacologically eliminated during the early postnatal period with a colony stimulating factor 1 receptor (CSF1R) inhibitor, excess calyces were not pruned and astrocytes did not mature properly in the auditory brainstem. Here, we show that once this drug is withdrawn, microglia gradually return to the auditory nuclei. After microglia reemerge in the medial nucleus of the trapezoid body (MNTB), synaptic pruning of calyces of Held resumes and maturation of astrocytes and auditory function recover. The findings suggest that the auditory brainstem pathways can be shaped by microglia even after their normal period of circuit development.
Introduction
Specialized neural circuits in the auditory brainstem carry out computations needed for auditory processing, including sound source localization (Grothe et al., 2010). Auditory nerve fibers synapse onto bushy cells in the anteroventral cochlear nucleus (AVCN; Held, 1893; Brawer and Morest, 1975; Ryugo and Sento, 1991; Lauer et al., 2013). Globular bushy cell (GBC) axon terminals form the elaborate calyx of Held (Held, 1893), which contacts principal neurons in the contralateral medial nucleus of the trapezoid body (MNTB; Morest, 1968; Smith et al., 1991; Kandler and Friauf, 1993). MNTB neurons inhibit neurons in the lateral superior olive (LSO; Kuwabara and Zook, 1991; Tollin, 2003). LSO neurons also receive tonotopically matched excitatory input from the ipsilateral side and use the balance between inhibition and excitation to compute interaural intensity differences (Boudreau and Tsuchitani, 1968; Sanes and Rubel, 1988; Grothe, 2003). This pathway, which allows for fast and reliable transmission of sound signals (Borst and Sakmann, 1996; Trussell, 1999; Taschenberger et al., 2002), relies on precision in developmental mechanisms.
Microglia mediate major functions in neural development (Reemst et al., 2016; Nelson and Lenz, 2017) and formation of neuronal circuits (Miyamoto et al., 2016; Basilico et al., 2019) through their roles in synapse elimination (Schafer and Stevens, 2010; Paolicelli et al., 2011; Schafer et al., 2012; Neniskyte and Gross, 2017), construction of new synapses (Parkhurst et al., 2013; Miyamoto et al., 2016), and synaptic remodeling (Tremblay et al., 2010; Schafer et al., 2012; Weinhard et al., 2018). When microglial function is disturbed, excess synapses persist in adulthood, resulting in anatomic abnormalities in neural circuits (Kettenmann et al., 2013).
Microglia are sparse in MNTB during the first postnatal week, then peak in density just after hearing onset at two weeks (Dinh et al., 2014). Microglia are found in close proximity to the developing calyx of Held (Holcomb et al., 2013; Dinh et al., 2014). From multiple small “protocalyces” contacting a single MNTB neuron (Kandler and Friauf, 1993), a single dominant calyx of Held emerges and monoinnervation of MNTB neurons is established (Morest, 1968; Sätzler et al., 2002; Hoffpauir et al., 2006; Holcomb et al., 2013). Microglial elimination before hearing onset results in excess polyinnervated neurons, signifying a microglial role in the development of neural circuits in the brainstem (Milinkeviciute et al., 2019).
Like microglia, astrocytes support synaptic development and remodeling. Despite differences between microglia and astrocytes in terms of their origin and structure, these cell types exhibit coordinated responses to insult (Schiweck et al., 2018; Vainchtein and Molofsky, 2020). Astrocytes and microglia communicate through bidirectional signaling and shared responses to environmental signals (Jha et al., 2019; Vainchtein and Molofsky, 2020). Microglia can influence expression of astrocytic proteins (Elmore et al., 2014; Spangenberg et al., 2016; Jin et al., 2017) and astrocytes can signal to microglia through several pathways (Bialas and Stevens, 2013; Vainchtein et al., 2018). In the auditory brainstem, microglial depletion prevented the emergence of astrocytic expression of glial fibrillary acidic protein (GFAP), a marker for mature astrocytes (Wofchuk and Rodnight, 1995; Gomes et al., 1999; Middeldorp and Hol, 2011), further highlighting the relationship between the two glial populations (Milinkeviciute et al., 2019).
Here, we tested whether microglia can return after postnatal pharmacological elimination, and whether microglial repopulation can restore development of auditory brainstem pathways. We further tested whether temporary postnatal microglial elimination has an effect on auditory function. We depleted microglia during the early postnatal period with an inhibitor of colony stimulating factor 1 receptor (CSF1R; Milinkeviciute et al., 2019), which is essential for microglial survival and proliferation (Stanley et al., 1997; Erblich et al., 2011). We then allowed microglia to repopulate by cessation of treatment after postnatal day (P)10. We found that microglia slowly recolonize the brainstem, first in lateral regions and later in progressively more medial regions. We found that defects present at P13 (Milinkeviciute et al., 2019), synaptic pruning and decreased expression of GFAP, were both corrected when microglia repopulated the MNTB. Additionally, we evaluated the overall auditory brainstem function of these animals using ABRs. We found elevated auditory thresholds, increased latencies and interpeak latencies, and decreased amplitudes in almost all ABR peaks in response to most of the frequencies tested at four weeks of age. ABRs were partially restored to normal by seven weeks of age. These results suggest that microglia have early developmental roles in synaptic pruning of calyces of Held, promoting astrocyte maturation, and maturation of auditory brainstem function. Our findings further suggest that microglia can act at later developmental times, at least in part, to restore normal formation of auditory brainstem circuits.
Materials and Methods
Animals and treatment
We used eight P14, seven P18, 18 three-week-old, 37 four-week-old, and 26 seven-week-old wild-type C57BL/6 mice of both sexes. All animal procedures were performed in accordance with the University of California, Irvine animal care committee’s regulations. Mice were housed in groups, with a maximum of five adult mice per cage. Mice were reared in a standard dark/light cycle, fed a standard diet and received food and water ad libitum. Whenever possible, mice that were used for ABRs were also used for histology with exception of pruning analysis. Each litter was divided into control and experimental mice. Experimental mice were repeatedly injected with BLZ945 (subcutaneous injections at P2, P4, P6, P8, and P10; MW: 398.48, MedChem Express HY-12768/CS-3971), a small molecule inhibitor of CSF1R (Pyonteck et al., 2013), to eliminate microglia as previously described (Milinkeviciute et al., 2019). BLZ945 was dissolved in dimethylsulfoxide (DMSO; D136-1; Fisher Scientific) and 0.01 ml of solution (200 mg/kg) was administered subcutaneously every 2 d starting from P2 with the last injection delivered at P10. Control mice were injected with DMSO following the same protocol. After cessation of treatment, mice were killed when they were 14 or 18 d or three, four, or seven weeks old (Fig. 1A).
Effects of BLZ945 treatment withdrawal on IBA1 expression in the auditory brainstem. A, The timeline of BLZ945 treatment. B, Areal coverage ratio of IBA1 labeling in MNTB. There was no difference between IBA1 area coverage ratio between control and BLZ945-treated mice at four weeks. With age, IBA1 coverage significantly decreased in control mice only. Row C, Images of IBA1 immunolabel (green) in AVCN, LSO, and MNTB in control mice at P14. LSO and MNTB are indicated by dashed lines. Rows D–G, Images of IBA1 immunolabel in AVCN, LSO, and MNTB in control mice at P18, three, four, and seven weeks. At seven weeks, microglia looked sparser in all three brainstem nuclei. Row H, Photographs of IBA1 immunolabel in AVCN, LSO, and MNTB in BLZ945-treated mice at P14. No IBA1-positive cells were observed throughout the auditory brainstem. Row I, Images of IBA1 immunolabel in BLZ945-injected mice at P18 showed only few IBA1-positive cells in AVCN but not LSO and MNTB. Row J, Photographs of BLZ945-treated mice at three weeks show AVCN populated by IBA1-positive cells and few of them were observed in LSO. MNTB was still devoid of microglia. Row K, Photographs of IBA1 immunolabel in AVCN, LSO, and MNTB of BLZ945-treated mice at four weeks. Microglia were spread throughout all the three nuclei. Row L, Photographs of IBA1 labeling in BLZ945-injected animals at seven weeks. IBA1 labeling appeared sparser but was still present in AVCN, LSO, and MNTB. Scale bar in AVCN panel in row C: 100 μm, applies to all AVCN panels across rows C–L. Scale bar in LSO panel in row C: 40 μm, applies to all LSO and MNTB panels across rows C–L; *p < 0.05.
Tissue preparation and immunofluorescence
Mice were weighed and perfused transcardially with 0.9% saline followed by 4% paraformaldehyde (PFA) in 0.1 m phosphate buffer, pH 7.3 (PBS) at three, four, or seven weeks (Fig. 1A). Brainstems were dissected and postfixed overnight in PFA solution. Brains were then transferred to a 30% sucrose solution in 0.1 m PBS overnight, embedded in OCT mounting medium, and sectioned coronally at 18 μm using a cryostat (CM 1850-3-1; Leica Microsystems). Tissue was mounted on chrome-alum-coated glass slides in a one-in-five series. Mounted sections were surrounded with a PAP pen hydrophobic barrier and rinsed in 0.1 m PBS for 10 min. Tissue was incubated for 5 min in 0.1% sodium dodecyl sulfate in 0.1 m PBS solution for antigen retrieval followed by 3 10 min washes in 0.1 m PBS. Sections were then blocked with normal goat blocking solution containing 5% normal goat serum (NGS; Vector Laboratories S-1000) and 0.3% Triton X-100 (Acros 9002-93-1) in 0.1 m PBS in a humidity chamber at room temperature. After 1 h, primary antibodies were applied (Table 1). The next day, the tissue was rinsed in 0.1 m PBS and incubated for 1 h in goat anti-rabbit or anti-chicken secondary antibody tagged with an Alexa (Invitrogen) fluorophore (Table 1). Sections were washed in 0.1 m PBS and incubated in blue or red fluorescent Nissl stain (NeuroTrace 435/455 or 530/615, Life Technologies N21482 or N21479) diluted in 1:200 in 0.3% Triton X-100 in 0.1 m PBS for 1 h. Tissue was then rinsed and coverslipped with Glycergel mounting medium (Dako C0563).
List of antibodies used in the project
Areal coverage analysis
We acquired 20× magnification images of sections throughout the rostro-caudal extent of the MNTB on both sides of the midline using a Zeiss Axioskop-2 microscope, an Axiocam camera, and Axiovision software. For each subject, a series of 20× multichannel fluorescent photographs spanning the rostro-caudal extent of the MNTB was taken and collected as an image stack in FIJI (Schindelin et al., 2012). The Nissl channel was used to outline the MNTB [a region of interest (ROI)] in every section in FIJI. A stack of ROIs from each animal was used to analyze the areal coverage of IBA1 or GFAP immunolabeling in the corresponding MNTB region as described in Milinkeviciute et al. (2019).
Neuronal tracing and confocal analysis
We used four control and four BLZ-treated mice and five control and five BLZ945-treated mice at three and four weeks, respectively, for neuronal tracing. After transcardial perfusion with artificial CSF (aCSF; 130 mm NaCl, 3 mm KCl, 1.2 mm KH2PO4, 20 mm NaHCO3, 3 mm HEPES, 10 mm glucose, 2 mm CaCl2, and 1.3 mm MgSO4 perfused with 95% O2 and 5% CO2), brains were quickly dissected and placed in a chamber with oxygenated aCSF. We temporarily placed the brain in a Petri dish with aCSF and filled calyces of Held in the MNTB with a rhodamine dextran amine (RDA; MW 3000, Invitrogen) solution (6.35% RDA with 0.4% Triton X-100 in PBS). Pulses of RDA were delivered through a pulled glass micropipette in the ventral acoustic stria (VAS) close to the midline using an Electro Square Porator (ECM830; BTX) at a rate of five pulses per second (pps) at 55 V for 50 ms. As a result, axonal projections from the AVCN terminating in calyces of Held in the MNTB were sparsely labeled with RDA on both sides of the brainstem. The brain was then transferred back into the aCSF chamber for ∼2 h under continuous oxygenation to allow dye transport along the GBC axons. The brain was then placed in 4% PFA overnight then transferred to 30% sucrose in 0.1 m PBS and later cryosectioned in the coronal plane at 18 μm. Sections were mounted on chrome-alum-coated glass slides in a one-in-five series and were immunolabeled with anti-VGLUT1/2 antibody (Table 1) following the protocol described above.
Stained slides were analyzed using confocal microscopy (Leica SP8, 63× oil objective, zoom: 1.5, pinhole: 1). Nissl, RDA, and VGLUT1/2 z-stack images of calyces of Held were acquired at a resolution of 1024 × 1024 pixels, with a z-step size of 0.5 μm (Grande et al., 2014; Wang et al., 2018; Milinkeviciute et al., 2019). Gain and offset were set for each fluorescent channel and each slide separately and adjusted if labeling intensity was noticeably different between sections on the same slide. Calyces throughout the MNTB were sparsely and randomly labeled with RDA, thus, the entire mediolateral extent of the MNTB was randomly sampled. If the number of labeled calyces from one slide was too small (<3), an additional slide was used (Table 2).
Numbers of calyces of Held analyzed and polyinnervated neurons
Image stacks were analyzed using the surface module in Imaris software (v9.5.1; Bitplane). Surfaces of RDA-filled calyces that were complete or near complete and that had adequate VGLUT1/2 immunofluorescence as well as a visible preterminal axon segment were reconstructed using a 0.4 surface detail and 0.4 background subtraction settings. These settings were determined before the analysis and were chosen so as to most reliably reconstruct the surface of the calyx of Held. If necessary, the rendering was adjusted manually to achieve an accurate RDA fill. Calyceal surface area as well as volume were recorded and used to calculate the shape factor of calyces from monoinnervated neurons. A shape factor was calculated as (SA1.5)/(6
V) where SA is the surface area of the calyx and V is its volume (DeHoff, 1978; Muniak et al., 2018). A shape factor of 1.0 characterizes a perfect sphere. A large shape factor corresponds with more complex structure, which, in the case of calyx of Held, represents its intricate branching. The analysis was performed blind to the animal treatment group.
Determination of monoinnervation and polyinnervation
The monoinnervation or polyinnervation status of MNTB neurons was evaluated during the reconstruction of RDA-filled calyces of Held. Neurons were classified as monoinnervated if the neuron was contacted only by a single RDA-labeled calyx, with no VGLUT1/2 immunolabel seen outside this calyx. Neurons were classified as polyinnervated if VGLUT1/2 immunolabel was present on the MNTB cell outside of the RDA-filled calyx. To be considered additional calyceal input, the VGLUT1/2 immunolabel outside of the RDA-filled calyx covered at least 25% of the MNTB neuron surface. The threshold of 25% was chosen because the smallest RDA-filled calyx detected across all the 10 analyzed animals contacted ∼25% of the neuronal surface area (Fig. 2). We then determined the percentage of analyzed MNTB neurons that were polyinnervated for the control and BLZ945-treated groups (Table 2).
Effects of BLZ945 treatment withdrawal on GFAP expression in the auditory brainstem. Row A, Images of GFAP (green) and IBA1 (magenta) immunolabel in AVCN, LSO, and MNTB in control mice at three weeks. Dashed lines indicate LSO and MNTB. GFAP expression was dense along outer edges of AVCN but sparse within the nucleus. GFAP expression was also sparse in LSO, but denser staining was observed in MNTB. Row B, Images of GFAP and IBA1 immunolabel in AVCN, LSO, and MNTB in BLZ945-treated mice at three weeks. GFAP immunolabel was seen along the outer edges of AVCN and was sparse in LSO. Some GFAP immunoloabel was seen in MNTB. Row C, Images of Nissl, GFAP immunolabel and an overlay of both in MNTB in control mice at four weeks. GFAP expression extended throughout the MNTB. Row D, Images of Nissl, GFAP immunolabel and an overlay of both in MNTB in BLZ945-treated mice at four weeks. GFAP expression was denser than at three weeks but was localized more to the boundaries of MNTB. Row E, F, Images of MNTB in control and BLZ945-injected mice, respectively, at seven weeks. GFAP immunolabeling was pronounced and found throughout MNTB. G, Areal coverage of GFAP expression in MNTB at four and seven weeks. GFAP expression was significantly decreased in BLZ945-treated animals at four weeks compared with age-matched controls. No difference was found between the two experimental groups at seven weeks, and there were no age-related changes in GFAP areal coverage ratio in control or BLZ945-injected mice. Scale bar in row A: 100 μm, applies to all panels in rows A, B. Scale bar in row C: 100 μm, applies to all panels in rows C–F; *p < 0.05.
Auditory brainstem recordings (ABRs)
ABRs were performed on 13 control and 14 BLZ945-treated mice at four weeks and on 11 control and 15 BLZ945 injected animals at seven weeks. Mice were anesthetized with an intramuscular injection of ketamine (75 mg/kg, KetaVed, VEDCO) and xylazine (15 mg/kg, AnaSed, NADA# 139-236). Body temperature was maintained at 35°C via a far infrared warming pad (Kent Scientific, RT-0501). Sterile ocular lubricant (Puralube Vet Ointment, 006PHM02-1-8) was applied on the eyes following ketamine/xylazine injection. Three pin electrodes were inserted subcutaneously. The positive, negative and ground electrodes were positioned at the vertex, the right cheek, and in the back near the right leg, respectively. Electrodes were connected to a Tucker-Davis Technologies (TDT) RA4PA 4-channel Medusa amplifier, which was connected to a TDT RA16 Medusa Base Station.
ABRs were performed in a sound-attenuating chamber (102 × 98 × 81 cm, Industrial Acoustics Company). Sound stimuli were generated using TDT SigGen software version 4.4 and presented 500 times at a rate of 21 stimuli per second using TDT MF1 Multi-Function Speaker through an ear tube inserted in the animal’s left ear. Acoustic stimuli were emitted via a TDT RP2.1 enhanced real time processor. We controlled the sound level via a TDT PA5 programmable attenuator. Recorded responses were amplified by a TDT SA1 stereo power amp, filtered with the control of BioSig software version 4.4. ABRs were recorded for 12 ms in response to 100-μs click stimuli or 3-ms pure tone stimuli (4, 8, 12, 16, 24, 32 kHz) at decreasing sound levels (5 dB SPL steps from 80 to 10 dB SPL). An averaged response was computed at each sound level and used for analysis.
ABR analysis
ABRs were analyzed for threshold, peak latency, interpeak latency, and peak amplitude. We computed the mean noise level from the first 20 data points, before the ABR trace. Threshold was defined as the lowest sound level at which the peak I level (μV) was ≥4 SDs above the noise level (Bogaerts et al., 2009). Peaks were manually determined by a blinded observer, and data were recorded to a spreadsheet in BioSig from which the following analyses were performed. Absolute peak latency was determined as the time period from the sound onset to the apex of the peak. Interpeak latency was measured by calculating the difference between absolute peak latencies between peaks I–II, II–III, III–IV, I–III, and I–IV. Peak amplitude was determined as the change in microvolts (μV) between the preceding trough and apex of the peak. Simple linear regressions of input-output functions of absolute peak latencies, interpeak latencies, and trough-to-peak amplitudes at recorded frequency and intensity levels were analyzed. Input-output functions were only analyzed if at least three animals per group were at or above threshold (Tables 3, 4, 5).
Descriptive statistics of immunostaining and confocal data
Two-way ANOVA and simple linear regressions
Sidak’s multiple comparisons
Experimental design and statistical analysis
We used multiple litters for each experimental group. Each litter contained pups injected with DMSO or BLZ945 to account for the variability between litters. Quantitative results are presented as mean scores ± SEM. Statistical analysis was done using Prism Software (v8.4.1; GraphPad Software). Comparisons between different treatment (control and BLZ945) and age groups (three and four weeks or four and seven weeks) were made using two-way ANOVA with Sidak’s multiple comparisons test unless otherwise indicated. Statistical significance was accepted at p < 0.05. Details of statistical analyses are presented in Tables 3-5.
Results
Microglia repopulate the MNTB after cessation of treatment
Treatment with BLZ945 early in development (repeated subcutaneous injections at P2, P4, P6, P8, and P10) was used to eliminate microglia in the brainstem of mice as previously described (Milinkeviciute et al., 2019). We first determined whether microglia return after treatment cessation and if so, to what extent and over what time course. We perfused mice at several posttreatment time points: P14, P18, three, four, and seven weeks (Fig. 1A). We used four DMSO and 4 BLZ945-injected mice at P14, three DMSO and 4 BLZ945-treated mice at P18, five DMSO and five BLZ945-treated mice at three weeks, six DMSO and six BLZ945-treated mice at four weeks and four DMSO and seven BLZ945-injected mice at seven weeks. Littermate controls injected with DMSO were used in each age group. Each cohort consisted of animals from a few different litters to account for the variability between litters. Microglia were identified using IBA1 immunofluorescence. The areal coverage of IBA1 immunolabel was used to evaluate the extent of microglia repopulation at four and seven weeks (Fig. 1B). There was no significant difference in the IBA1 areal coverage in the MNTB (four-week DMSO: 0.06 ± 0.004; four-week BLZ945: 0.05 ± 0.01; seven-week DMSO: 0.04 ± 0.01; seven-week BLZ945: 0.04 ± 0.003) between control and drug treated animals at both four (p = 0.0574, t = 2.360; df = 19) and seven (p = 0.9113, t = 0.3882; df = 19) weeks. IBA1 immunolabeling significantly decreased with age only in control animals (DMSO: p = 0.0166, t = 2.943; df = 19; BLZ945: p = 0.8425, t = 0.5286; df = 19), but there was no difference between DMSO controls and BLZ945-treated mice at seven weeks (p = 0.9113, t = 0.3882; df = 19; Fig. 1B).
Microglia densely populated VCN, LSO, and MNTB of control mice in each age group (Fig. 1C–G). In BLZ945-injected mice, microglia were not observed in the VCN, LSO or MNTB at P14 (Fig. 1H). A few microglial cells were seen in the VCN at P18 but no IBA1 staining was observed in the LSO or MNTB (Fig. 1I). By three weeks, microglia were seen throughout VCN and in LSO, but not in MNTB (Fig. 1J). At four and seven weeks, microglia were present in all three auditory brainstem nuclei (Fig. 1K,L). This gradual microglial repopulation, first in VCN, then LSO, then MNTB reflects an overall lateral-to-medial movement of microglial emergence that is seen throughout the brainstem in our sections, reflecting a similar trend to how microglia occupy the brainstem during the first two postnatal weeks (Dinh et al., 2014). Thus, microglial repopulation of MNTB takes place between the ages of three and four weeks, or 11–18 d after the treatment is discontinued.
GFAP expression in the MNTB returns to control levels after return of microglia
GFAP is regarded as a marker of mature (Wofchuk and Rodnight, 1995; Gomes et al., 1999; Middeldorp and Hol, 2011) or active astrocytes (Eng and Ghirnikar, 1994; Fuentes-Santamaria et al., 2017). In MNTB, GFAP expression changes during development, with no expression at P0, very sparse expression at P6, and increased density at P14. GFAP staining in VCN is first observed at P14 and some positive fibers are seen in MNTB at P23 (Dinh et al., 2014). We looked at 20× confocal images of sections through the VCN, LSO, and MNTB stained for IBA1 and GFAP at three weeks (11 d postcessation of treatment) in three DMSO and three BLZ945-treated mice. As noted above, microglia repopulated auditory brainstem sequentially, starting from the most peripheral nuclei. At three weeks, microglia were found throughout the VCN, LSO, and MNTB of control mice. In VCN, GFAP labeling was mostly confined to outer layer of the nucleus with sparse labeling in the deeper layers. In LSO, GFAP-positive cells were also sparse while MNTB appeared to be the most populated (Fig. 2A). Similar to control animals, in BLZ945-injected animals, microglia were spread throughout the VCN, however, in the LSO fewer cells were observed, and they appeared larger than in the VCN or LSO of control animals. MNTB was still devoid of IBA1-positive cells. In VCN GFAP was largely restricted to outer layers. GFAP staining in LSO was also very similar to that of control mice. In MNTB, GFAP-positive cells were observed throughout the nucleus (Fig. 2B). We showed that microglia fully repopulated MNTB by the age of four weeks, thus, next we focused on GFAP expression at four and seven weeks.
When microglia were depleted early postnatally, GFAP areal coverage was significantly reduced at P13 (Milinkeviciute et al., 2019). We examined GFAP areal immunofluorescence coverage ratio in 11 control and six BLZ945-treated mice at four weeks as well as six control and eight BLZ945-injected animals at seven weeks. Abundant GFAP was observed in four-week-old control mice (four-week DMSO: 0.05 ± 0.01; Fig. 2C). In drug treated four-week animals GFAP expression (four-week BLZ945: 0.04 ± 0.01; Fig. 2D) was localized more to the peripheral parts of the MNTB as reported in Milinkeviciute et al. (2019) at P13 while in control tissue GFAP labeling was found throughout the MNTB (Fig. 2A). At seven weeks, both, control and BLZ945-treated mice showed abundant GFAP labeling throughout the MNTB (seven-week DMSO: 0.04 ± 0.003; seven-week BLZ945: 0.04 ± 0.002 (Fig. 2E–G). The areal GFAP coverage was significantly reduced in BLZ945 injected mice at four weeks (p = 0.0258, t = 2.66; df = 27) but was comparable to control levels at seven weeks (p = 0.8450, t = 0.5215; df = 27; Fig. 2G). There was no significant age-related change in GFAP areal coverage ratio in any of the two experimental groups (DMSO: p = 0.2535, t = 1.537; df = 27; BLZ945: p = 0.8381, t = 0.5341; df = 27). These results suggest a relationship between microglia and the GFAP-positive astrocytic population in MNTB.
MNTB neurons are monoinnervated after microglia recovery
We labeled calyces of Held in the MNTB on both sides by injecting RDA into the VAS at the midline. Tissue was labeled with VGLUT1/2 antibody to identify excitatory inputs, which in MNTB mainly corresponds to calyces of Held (Billups, 2005). MNTB neurons contacted by an RDA labeled calyx were used for confocal imaging and subsequent reconstruction in 3D as previously described in Milinkeviciute et al. (2019). We reconstructed and analyzed calyces of Held from four DMSO and four BLZ945 at three weeks and five DMSO and five BLZ945-injected mice at four weeks. Monoinnvervated and polyinnervated MNTB neurons were identified in both control and BLZ945-treated mice at three weeks (Fig. 3A–D) and four weeks (Fig. 3E–H). We evaluated each MNTB neuron for VGLUT1/2 labeled inputs outside the RDA-labeled calyx (see Materials and Methods) to determine the percentage of neurons contacted by more than one calyx (three-week DMSO: 9.92 ± 6.75; three-week BLZ945: 38.78 ± 9.89; four-week DMSO: 5.08 ± 3.15; four-week BLZ945: 16.59 ± 5.82). At three weeks, the percentage of polyinnervated neurons was significantly increased in the BLZ945-treated mouse group compared with their age-matched controls (p = 0.0191, t = 2.998; df = 14). We found no significant difference in polyinnervation of MNTB neurons between control and drug treated animals at four weeks (p = 0.3645, t = 1.336; df = 14; Fig. 3I) Numbers of animals and calyces analyzed are presented in Table 2. These results suggest that neurons in MNTB remain polyinnervated when microglia are absent from MNTB.
Confocal analysis of MNTB neurons. A, Monoinnervated neuron from a three-week-old control mouse. Reconstructed RDA-labeled calyx of Held (magenta); MNTB neuron (blue); and VGLUT1/2 labeling (green). B, Monoinnervated neuron from a three-week-old BLZ945-treated mouse. C, D, Monoinnervated neuron from a four-week-old control and BLZ945-injected mouse, respectively. E, F, Polyinnervated neuron from a three-week-old control and BLZ945-injected mouse, respectively. Dashed line indicates the contour of the neuron while asterisks show the VGLUT1/2-positive terminal contacting the same neuron. G, H, Polyinnervated neuron from a four-week-old control and BLZ945-treated animal, respectively. An extra calyceal input was determined by a continuous VGLUT1/2 staining outside of the RDA-positive area (asterisks) surrounding the MNTB neuron. I, Comparison of calyx of Held surface area in control and BLZ945-treated mice at three and four weeks. There was no difference between the two animal groups at both ages tested. No age-related increase in calyceal area was observed. J, Comparison of polyinnervated MNTB neurons as a percentage of MNTB neurons analyzed revealed significantly more polyinnervated neurons in BLZ945-injected mice at three weeks when compared with controls. No difference between control and BLZ945 injected mice was found at four weeks. The age-related change in BLZ945-injected animals did not reach significance (p = 0.0574). Scale bars: 5 μm (A, B, E, F), 3 μm (C, G, H), and 4 μm (D).
We measured the surface area and volume of RDA-positive and reconstructed calyces from monoinnervated neurons. Calyces in the MNTB of control and drug injected mice did not differ in size at three weeks (p = 0.9895, t = 0.131; df = 14) or four weeks (p = 0.7005, t = 0.7724; df = 14). There was no significant age-related change in calyx surface area in control or BLZ945-treated mice (DMSO: p = 0.4941, t = 1.103; df = 14; BLZ945: p = 0.8527, t = 0.5126; df = 14; Fig. 3J). Complexity of calyces was measured by calculating a shape factor using the formula described in see Materials and Methods (Bogaerts et al., 2009; Muniak et al., 2018). Calyces did not become more complex with age (DMSO: p = 0.3715, t = 1.323; df = 14; BLZ945: p = 0.3717, t = 1.322; df = 14; Fig. 3K). There was no difference observed between the intricacy of calyceal branching between control and BLZ945-injected mice at three weeks (p = 0.3210, t = 1.425; df = 14) or at four weeks (p = 0.2488, t = 1.594; df = 14; Fig. 3K; Table 3). In summary, temporary microglial depletion and subsequent repopulation does not seem to affect the growth and morphologic maturation (in terms of branching complexity) of calyces of Held.
Temporary postnatal microglia elimination alters ABRs
Temporary microglia depletion results in elevated hearing thresholds
Auditory brainstem responses were measured in mice at four weeks (n = 13 DMSO, n = 14 BLZ945) and seven weeks (n = 11 DMSO, n = 15 BLZ945). Click and pure tone stimuli were presented to the left ear at decreasing intensities from 80 dB SPL to 10 dB SPL. Click and pure tone traces were analyzed for threshold, the lowest intensity at which peak I level (μV) was equal to or exceeded the 4 SDs above the mean noise level (Fig. 4A). BLZ945-treated mice had significantly higher thresholds in response to clicks at both four and seven weeks (four weeks: p < 0.0001, df = 47; seven weeks p = 0.0321, df = 47; one-way ANOVA with Tukey’s multiple comparisons; Fig. 4B). Age did not affect thresholds in both treatment groups (control: p = 0.1884, df = 47; BLZ945: p = 0.9999, df = 47; one-way ANOVA with Tukey’s multiple comparisons; Fig. 4B). At four weeks, BLZ945-treated mice displayed elevated thresholds in response to pure tones at 4, 8, and 12 kHz (4 kHz: p < 0.0001, df = 24.89; 8 kHz: p < 0.0001, df = 24.80; 12 kHz: p = 0.0007, df = 23.65; 16 kHz: p = 0.2228, df = 23.98; 24 kHz: p = 0.2179, df = 16.58; 32 kHz: p = 0.3584, df = 24.08; Fig. 4B). At seven weeks, BLZ945-treated mice displayed pure tone thresholds that were comparable to control, and there were no significant differences at any frequency tested (4 kHz: p = 0.1043, df = 24.0; 8 kHz: p = 0.4413, df = 23.63; 12 kHz: p = 0.8559, df = 24.0; 16 kHz: p = 0.4524, df = 20.73; 24 kHz: p = 0.8303, df = 20.54; 32 kHz: p = 0.9994, df = 22.56; Fig. 4C). Example traces indicating thresholds from each group at low and high frequencies are displayed in Figure 4D. These results indicate a temporary hearing impairment as a result of CSF1R inhibition during development and recovery after a prolonged period of microglial repopulation.
Effect of temporary BLZ945 treatment on auditory thresholds. A, Example trace for method of ABR analysis. Threshold was determined as the sound intensity for which peak I amplitude was greater than or equal to SDs above the averaged noise before the beginning of the ABR trace. Latency (blue) is the time of the apex of the peak from stimulus onset. Interpeak latency (red) is the difference in time between two respective peaks. Amplitude (black) is the difference in microvolts between the preceding trough to the apex of the peak. B, Averaged thresholds for responses to click stimuli at four and seven weeks. BLZ945 injected mice had significantly elevated thresholds compared with their age-matched controls at four and seven weeks. C, Averaged auditory thresholds for responses to pure tone stimuli at four weeks (black) and seven weeks (gray). BLZ945-treated animals had significantly elevated thresholds in response to 4, 8, and 12 kHz at four weeks. Treatment did not affect thresholds at seven weeks, compared with age-matched controls. D, Example traces in response to pure tones in control and BLZ945-treated mice at four and seven weeks of age. Thresholds are indicated with black arrows.
Peak latencies are increased in BLZ945-treated animals
Peaks in the ABR waveform reflect neuronal activity in different parts of the auditory pathway. Peak I indicates activity in the cochlea, spiral ganglion neurons (SGN), and VIIIth nerve, peak II represents, the cochlear nucleus (CN), peak III, the superior olivary complex (SOC), and peak IV, the lateral lemniscus (LL; Jewett et al., 1970; Jewett and Williston, 1971; Picton et al., 1974; Henry, 1979). Mice were tested at four and seven weeks to compare the auditory function at early and late time points after microglial repopulation. All ABR statistical analyses including input-output simple linear regressions are displayed in Tables 4, 5.
Peak latency was determined as the time (ms) from stimulus at which the apex of the peak was detected (Fig. 4A). At four weeks, BLZ945-treated mice showed significantly increased peak I latencies in lower and middle frequencies. Higher frequencies displayed significantly decreased peak I latencies (Fig. 5A). At seven weeks, treatment effects on peak I latencies were still observed at the lower and middle frequencies, but no differences were seen at the higher frequencies (Fig. 5A), indicating some recovery in peak I latencies. Peak II latency in BLZ945-treated mice was significantly increased in low and middle frequencies but significantly decreased in high frequencies at four weeks (Fig. 5B). By seven weeks, peak II latency remained significantly delayed in the lower and middle frequencies but recovered at high frequencies (Fig. 5B). BLZ945 treatment led to significantly increased peak III latency (Fig. 5C) in most frequencies tested at four weeks. These treatment effects remained at some frequencies at seven weeks, but multiple comparison tests did not show any intensity-specific differences (Tables 4, 5), suggesting an overall decrease in the effect of the treatment. In peak IV (Fig. 5D), BLZ945 treatment resulted in significantly increased latencies at four weeks in the lower and middle frequencies, but not the high frequencies. At seven weeks, latencies showed some recovery in the middle frequencies, but BLZ945-treated mice still showed significantly increased latencies at 32 kHz (Fig. 5D). Together, these results indicate that temporary CSF1R inhibition early in development affects peak latencies at four weeks. Some, but not all, of these latency changes recover after cessation of drug treatment by seven weeks.
Effect of temporary BLZ945 treatment on ABR peak latencies. A, Absolute peak I latency in response to 8-, 16-, and 32-kHz stimuli at four and seven weeks. Treatment with BLZ945 early in development resulted in significantly increased latency at the low and middle frequencies but significantly decreased latencies at the higher frequencies at four weeks. By seven weeks, the latencies at the lower frequencies remained elevated but latencies at the high frequencies were comparable to control. B, Peak II latency was significantly increased in the BLZ945-treated group. Latencies remained elevated by seven weeks at all frequencies except 32 kHz. C, BLZ945 treatment resulted in significantly increased peak III latency, with evidence of some recovery by seven weeks in all frequencies. D, Significantly increased latencies in peak IV were detected in BLZ945-treated mice in low and middle frequencies at four weeks. Some improvement was detected by seven weeks, but 32-kHz stimuli resulted in significantly longer latency in the BLZ945 group.
Increased interpeak latencies partially recover by seven weeks of age
We next measured interpeak latencies to test whether the observed latency differences in the BLZ945-treated mice is because of peripheral conduction defects, reflect any delays in central activity or is a result of both. Interpeak latency was determined as the difference in time between the apex of the peak and the apex of the preceding peak (ms). At four weeks, treatment effects on interpeak latencies were observed between peaks I and II at 8 and 24 kHz, where BLZ945-treated mice showed significantly decreased interpeak latencies. However, by seven weeks, there were no differences in peak I–II interpeak latencies, indicating that cochlea-SGN-VIIIth nerve signal transduction recovered (Fig. 6A). BLZ945 treatment led to significantly increased peak II–III interpeak latencies at four weeks at 8, 12, 16, 24, and 32 kHz (Fig. 6B). At seven weeks, BLZ945-treated mice showed normal peak II–III latencies at the lower and middle frequencies, but the higher frequencies remained delayed (Fig. 6B). Together, these data suggest that signal transmission from the CN to the SOC is delayed after temporary microglia depletion but recovers by at seven weeks. Peak III–IV interpeak latencies (Fig. 6C) were significantly elevated at 4, 8, 12, 16, and 24 kHz in BLZ945-treated mice at four weeks. Post hoc multiple comparison analyses showed that these interpeak latency differences were more apparent in the lower frequencies. Similar results were obtained at seven weeks, indicating that these increases in peak III–IV latencies do not recover with time (Fig. 6C).
Effect of temporary BLZ945 treatment on ABR interpeak latencies. A, ABR peak I–II interpeak latencies in response to 8 kHz were significantly decreased in BLZ945-treated mice at four weeks. By seven weeks, there were no differences at any frequency level. B, Peak II–III interpeak latency was significantly increased in BLZ945-treated mice at all frequency levels. By seven weeks, low and middle frequencies recovered, but high frequencies remained significantly elevated. C, Peak III–IV interpeak latencies for the BLZ945 group were significantly longer in low and middle frequencies in both ages, while the high frequencies were not affected. D, Peak I–III interpeak latency was significantly elevated in BLZ945-treated mice at all frequency levels at four weeks. By seven weeks, only the low frequency was comparable to the control group. E, BLZ945 treatment resulted in significantly elevated peak I–IV interpeak latency at all frequencies in both ages.
To investigate the overall effects of temporary CSF1R inhibition on central latencies, we measured interpeak latencies between peaks I–III and I–IV. BLZ945 treatment led to increased interpeak latencies between peaks I and III at four weeks at all frequencies tested when compared with control mice. At seven weeks, peak I–III latencies improved at most frequency levels except for 16 and 32 kHz (Tables 4, 5). Treatment effect comparisons for peak I–IV latencies showed a significant increase in the BLZ945-treated group at all frequencies tested at both four and seven weeks compared with controls (Fig. 6D). However, post hoc comparisons did not display any significant changes for any particular intensity at seven weeks, indicating some recovery in central latencies after prolonged microglia repopulation. Together, these findings demonstrate that BLZ945 treatment in the early postnatal period increases latencies in the auditory brainstem measured at four weeks, and that these effects are partially improved by seven weeks.
Peak amplitudes are diminished at four weeks but recover by seven weeks of age
We next examined the effect of BLZ945 treatment on ABR peak amplitudes, determined as the difference in μV from the preceding trough to the following peak. Peak I (Fig. 7A) amplitude showed that there was a significant decrease in treated animals across all frequencies except 8 and 32 kHz at four weeks. Post hoc multiple comparisons tests did not show significant differences at any specific intensity at all frequency levels. At seven weeks, peak I amplitude did not differ between BLZ945-treated mice and control mice. Peak II amplitude (Fig. 7B) was significantly diminished at 8, 12, 16, and 24 kHz at most intensities in BLZ945-treated mice at four weeks. By seven weeks, peak II amplitudes were comparable to control levels at all frequencies and intensities (Fig. 7B). BLZ945 treatment at four weeks resulted in diminished peak III amplitudes at the low and middle frequency levels (4, 8, 12, and 16 kHz; Fig. 7C). By seven weeks, peak III amplitudes in treated animals were comparable to controls at all frequencies and intensities, except for 32 kHz where a significant increase in amplitude was detected in the BLZ945-treated group (Fig. 7C). Peak IV (Fig. 7D) amplitudes were significantly decreased in BLZ945-trreated animals at four weeks at all frequency levels. By seven weeks, amplitudes at most frequencies recovered, except for 16 kHz where a significant decrease in amplitude was detected in BLZ945-treated mice (Fig. 7D). Multiple comparison analyses revealed no difference at any intensity at seven weeks, indicating that peak IV amplitudes were generally recovered. Overall, these data show that at four weeks, peak amplitudes are significantly decreased in BLZ945-treated animals compared with controls (see Tables 4, 5). After prolonged microglial repopulation, there is evidence of recovery of peak amplitudes throughout the auditory brainstem pathways.
Effect of temporary BLZ945 treatment on ABR amplitudes. A, Trough to peak amplitudes (μV) of ABR peak I in response to 8-, 16-, and 32-kHz stimuli at four and seven weeks. Amplitudes were comparable to control at most frequency levels except 16 kHz, where amplitudes were significantly diminished. B, Peak II amplitude was significantly decreased in the BLZ945-treated group at the low and middle frequencies, with some improvement by seven weeks. C, Peak III amplitudes were significantly lower in the BLZ945 mice at low and middle frequencies at four weeks, and recovered by seven weeks. The higher frequency showed that peak III amplitude was significantly elevated in the BLZ945-treated group at seven weeks. D, BLZ945-treated mice showed significantly decreased peak IV amplitudes at all frequency levels at four weeks. By seven weeks, peak IV amplitude was comparable to control at 8 and 32 kHz.
Discussion
In the present study, we used a CSF1R inhibitor to eliminate microglia early in development and then withdrew this drug to study microglia repopulation and its effects on auditory brainstem circuits. Previously, it was shown that such microglial depletion led to a decrease in GFAP expression as well as retention of polyinnervation in MNTB after the onset of hearing (Milinkeviciute et al., 2019). In this study, we showed that microglia recolonized the brainstem over an extended time course and polyinnervation was present when microglia were still absent from the MNTB at three weeks of age. The majority of MNTB neurons were no longer polyinnervated in BLZ945-treated mice at four weeks of age, a time when microglial numbers were comparable to control levels. Reduced GFAP expression took longer to reach that of controls; however, it was rectified by seven weeks. Functionally, early microglia depletion led to defects in auditory peak latencies and amplitudes, which largely recovered by the later age. This study emphasizes the plasticity of the brain after a challenge to the system and calls attention to the relationship between microglia and auditory function and between different glial cell types.
Microglial repopulation
Microglia can be eliminated using genetic methods (Bruttger et al., 2015; Wieghofer et al., 2015), pharmacological methods (Pyonteck et al., 2013; Gerber et al., 2018), or irradiation (Simard and Rivest, 2004; Han et al., 2019). We used subcutaneous BLZ945 injections early in development and achieved a complete elimination of IBA1-positive cells in the brainstem. Oral BLZ945 treatment in the five-week cuprizone mouse model resulted in reduced microglial numbers (up to 60%; Beckmann et al., 2018), highlighting the differences between the delivery method as well as age of the animal and region of the brain analyzed despite of the same inhibitory molecule used. Contrasting results between the two studies also extend to the timing of repopulation. After oral administration of BLZ945, active microglia appeared 3 d posttreatment and normal numbers were reached at 7 d (Beckmann et al., 2018). Similarly, microglial depletion with PLX3397 in adult mice resulted in re-emergence of immature microglia 3 d after the cessation of treatment. These cells gradually acquired adult-like phenotypes and reached control density after 14 d (Elmore et al., 2014). In our study, microglia did not populate the MNTB until after 18 d posttreatment, at four weeks of age. The difference in repopulation timing likely lies in our method in elimination of microglia early during postnatal development when microglial cells should be proliferating and populating the entire brain. BLZ945 treatment not only depleted microglia but, in part, prevented them from completing their initial brainstem colonization.
Origins of repopulating microglia
Microglia have a low turnover rate and can survive throughout the lifespan of an animal (Fuger et al., 2017). They can self-regulate and renew several times throughout life in rodents (Askew et al., 2017; Tay et al., 2017). We successfully eliminated microglia early in development and observed that, after cessation of treatment, microglia repopulate the brainstem. The origins of the reappearing brainstem microglia are not known. Studies implicate that repopulation can be achieved in one of three ways – proliferation of remaining microglia, stimulation of microglial precursors or infiltration of peripheral monocytes (for review, see Han et al., 2019). The process of repopulation likely depends on the method of depletion. It is possible that irradiation leads to an influx of macrophages from the periphery because of the disruption of the blood-brain barrier (BBB; Simard and Rivest, 2004). On the other hand, pharmacological depletion with drug-infused rodent diet could induce the proliferation from resident microglial progenitor cells such as nestin positive cells (Elmore et al., 2014). Additionally, proliferation of surviving microglial cells could support the microglial replenishment in the CNS. Genetic ablation of microglia (80%) using the Cx3cr1CreER-based system showed that repopulating microglial cells arose exclusively from an internal CNS-resident cell pool (Bruttger et al., 2015; Askew et al., 2017; Jakel and Dimou, 2017; Huang et al., 2018). Earlier it was shown that repeated subcutaneous BLZ945 injections resulted in complete microglial elimination in the brainstem (Milinkeviciute et al., 2019), thus, new cells do not likely emerge from surviving microglia. Peripheral macrophages could be another origin if the BBB was disrupted; however, despite the evidence that BLZ945 crosses the BBB (Pyonteck et al., 2013; Beckmann et al., 2018), there are no reports about disruption of the BBB. Consequently, CNS resident progenitors might be recruited to replenish the microglial cell pool after BLZ945 treatment is stopped. These studies support the view that the course of repopulation highly depends on the depletion approach used and the level of microglial elimination.
Microglia and calyceal pruning
Microglia are continuously surveying their environment by extending and retracting their processes (Nimmerjahn et al., 2005). During development, microglia have been shown to have a significant role in the formation of neural circuits (Miyamoto et al., 2016; Reemst et al., 2016; Nelson and Lenz, 2017; Basilico et al., 2019; Nonaka and Nakanishi, 2019) in health and pathologic conditions (Neniskyte and Gross, 2017; Wolf et al., 2017; Whitelaw, 2018). Microglia can eliminate synapses (Paolicelli et al., 2011; Schafer et al., 2012) or promote their formation (Schafer et al., 2012; Parkhurst et al., 2013; Miyamoto et al., 2016). In MNTB, microglia appear at the time coinciding with intense synaptic development. They peak in number at around P14 and a slight decline is observed at P23 (Dinh et al., 2014). We noticed a further significant reduction in microglial areal coverage between weeks 4 and 7. This change may be reflective of a shift in the function of microglia from a developmental role to a homeostatic role.
In the auditory brainstem, microglial elimination resulted in increased numbers of polyinnervated neurons remaining in the MNTB after hearing onset (Milinkeviciute et al., 2019). It was not clear, however, whether the calyceal pruning was delayed or directly related to the absence of microglia from MNTB. Here, we found that even when treatment is stopped at P10 and microglia are allowed to populate the brain, innervation of principal MNTB neurons remains abnormal until microglia inhabit the nucleus. It has been shown that microglia promote brain repair after injury (Willis et al., 2020). In our case there was no injury, and microglia appeared in the MNTB between weeks 3 and 4 for the first time. It appears that newly emerged microglia initialized the developmental process that normally occurs before the onset of hearing. Depletion and subsequent recolonization of microglia in the retina showed restoration of functions of endogenous microglia in homeostasis (Zhang et al., 2018). It is important to note that calyces of Held developed normally in terms of their surface area and complexity; however, the pruning of additional calyces was impaired, indicating that there may be different mechanisms that regulate these two processes. These findings suggest a correlation between microglial presence and pruning of calyces of Held in MNTB.
Microglia-astrocyte communication
Crosstalk between glial cells is essential for brain development and function in physiological and pathologic conditions. Bidirectional communication between microglia and astrocytes is established through cytokines, chemokines, ATP and growth factors early in development when they first populate the brain parenchyma (Jha et al., 2016; Liddelow et al., 2017; Vainchtein et al., 2018). Previously, it was shown that GFAP expression is reduced after microglia depletion (Milinkeviciute et al., 2019). Here, we found restoration of GFAP expression to control levels only after prolonged microglial repopulation. In contrast, oral BLZ945 treatment in the five-week cuprizone mouse model increased the number of GFAP-positive and ALDH1L1-positive astrocytes in the cortex, striatum, and corpus callosum (Beckmann et al., 2018). Similarly, PLX3397 treatment in adult mice resulted in a significant increase in GFAP mRNA (Elmore et al., 2014; Jin et al., 2017). This discrepancy may be attributed to the depletion method used; the delivery method; the developmental stage of the animal at the time of treatment; and/or the area of the brain analyzed. Here, microglia were depleted at the time when the brainstem should be undergoing the initial colonization by microglia (Dinh et al., 2014). This process was halted by injection of CSF1R inhibitor before microglia had a chance to occupy their positions and mature within the auditory brainstem. Given that GFAP is suggested to be a marker of mature (Wofchuk and Rodnight, 1995; Gomes et al., 1999; Middeldorp and Hol, 2011) astrocytes, it is possible that astrocytic maturation is impeded by microglial absence and when microglia are permitted to repopulate, GFAP expression is initiated and becomes comparable to that of control at seven weeks. It is yet to be determined which molecules may be involved in the process of astrocytic maturation.
Auditory brainstem function
Auditory function requires precision in the development of neuronal circuits as well as accurate organization of specialized synapses. Microglia are known to be necessary for the formation of some neural circuits (Miyamoto et al., 2016; Basilico et al., 2019; Nonaka and Nakanishi, 2019). We thus examined the overall effect of microglia depletion on the auditory function. We found that temporary depletion of microglia in the brainstem early in development resulted in deficits in the ABR. BLZ945 treatment during the first postnatal week caused a moderate hearing loss but thresholds recovered by seven weeks, while click thresholds did not show recovery at this later time period. We also found that BLZ945 treatment had a significant effect on absolute peak latencies as well as interpeak latencies at most frequencies tested, and at seven weeks, these latency defects were still mostly present. However, the diminished amplitudes observed at four weeks in BLZ945-treated mice showed almost full recovery by seven weeks.
Our observation of elevated ABR thresholds could indicate a peripheral defect in treated mice. Defective CSF1R signaling affects bone resorption and leads to abnormalities in the middle and inner ears (Aharinejad et al., 1999; Kanzaki et al., 2011; Okano and Kishimoto, 2019). In the cochlea, CSF1R-positive macrophages are found in the spiral ligament, spiral ganglion and stria vascularis (Okano et al., 2008; Okano and Kishimoto, 2019) while centrally, microglia rely on CSF1 signaling for survival (Ginhoux et al., 2010; Erblich et al., 2011). Thus, treatment with CSF1R inhibitor can affect many areas in the auditory pathway. Our ABR results are consistent with hearing deficits reported in Csf1op/op mice where auditory thresholds were significantly elevated likely because of bone deformities and subsequent conductive hearing loss (Aharinejad et al., 1999; Okano and Kishimoto, 2019). We observed delayed absolute latencies in subsequent peaks as well as interpeak latencies, consistent with brainstem auditory evoked potentials recorded in Csf1op/op mice (Michaelson et al., 1996). Central defects could result in part from impaired peripheral function.
In addition to potential peripheral effects, our ABR analyses suggest that BLZ945 may produce central effects independently through its elimination of microglia. We observed significant increases in interpeak latencies in treated animals. A minimal effect was seen for peak I–II, but greater effects were seen for all of the central intervals. These delays were partially corrected after microglia repopulation. Peak III corresponds to the superior olivary complex and reflects in part activity in MNTB. However, it is not known to what extent calyceal pruning directly influences the latency or amplitudes of this peak in our study. Nevertheless, the observation of increased interpeak latency, impaired pruning, and delayed astrocyte maturation in MNTB together support a central effect of microglia elimination on auditory function.
Microglia have been shown to communicate with other glial cells, including oligodendrocytes (Lloyd et al., 2017) and NG2 cells (Müller et al., 2009; Giera et al., 2018) and regulate oligodendrocyte differentiation (Miron and Franklin, 2014). Microglial depletion has revealed a microglial role in the development and maintenance of oligodendrocytes and their progenitors (Hagemeyer et al., 2017). Oligodendrocytes are responsible for myelination of axons and, thus regulate the precise signal conduction seen in auditory pathways (Seidl, 2014; Sinclair et al., 2017; Long et al., 2018). In the MNTB, myelination starts by P9 (Saliu et al., 2014) and increases around the time of ear canal opening (Sinclair et al., 2017). Taken together, microglial elimination may affect conduction velocity (Sinclair et al., 2017).
We showed that treatment with CSF1R inhibitor temporarily reduced peak amplitudes, especially the central peaks. Thus, there are both peripheral and central defects in the ABR. On one hand, bone deformities in the middle and inner ear (Aharinejad et al., 1999; Kanzaki et al., 2011; Okano and Kishimoto, 2019) could have contributed to these impairments which then propagated to the central auditory system, however, we did not test this hypothesis directly. On the other hand, centrally, oligodendrocytes are known to detect and respond to neuronal activity (Gibson et al., 2014) as well as secrete neurotrophic factors. In the auditory brainstem, lack of BDNF secreted by oligodendrocytes significantly reduced amplitudes of ABR peaks II–IV in response to click stimuli (Jang et al., 2019). NG2 cells near calyces of Held were suggested to participate in regulating the fast-signaling properties of calyceal inputs (Müller et al., 2009). Thus, microglial elimination could have affected oligodendrocytes or NG2 cells, leading to impaired neuronal activity and synaptic synchrony which resulted in the reduction of amplitudes we observed.
The observed decrease in ABR amplitudes may be also related to delayed maturation of GFAP positive astrocytes. Astrocytes are highly fenestrated and interspersed between the branches of calyces (Ford et al., 2009; Holcomb et al., 2013). GFAP astrocytes express glutamate transporters and regulate glutamate uptake activity (Bergles and Jahr, 1997; Matthias et al., 2003), and are crucial for its recycling (Pow and Robinson, 1994; Laake et al., 1995; Hertz et al., 1999). These glial cells respond to glutamate release from the calyx of Held, which leads to the uptake of excess glutamate and prevention of glutamate receptor saturation, especially in the immature calyx (Bergles and Jahr, 1997; Matthias et al., 2003). In this paper, we showed that temporary microglia depletion during postnatal development reduced GFAP+ astrocyte levels. It is thus possible that delayed astrocytic maturation led to impairment in glutamate uptake and allowed for glutamate receptor saturation, and subsequent decrease in ABR amplitudes of the central peaks.
At this point, we cannot predict what exactly causes deterioration of different aspects and parts of ABRs as microglia are found through the auditory system. It is possible that macrophages in the cochlea and microglia in different auditory nuclei have slightly different functions. For peak I, which reflects peripheral activity, frequency-specific effects could reflect tonotopic differences in the cochlea of molecules that interact with macrophages. Centrally, there may also be differences in the way tonotopic regions of individual nuclei interact with glial cells.
Because of the fact that ABR peak III encompasses the auditory response of the entire SOC, we cannot correlate anatomic abnormalities observed in MNTB with ABR results. Rather, we can conclude that temporary loss of CSF1R signaling affects conduction velocity and strength throughout the ascending auditory pathway.
In conclusion, microglia function in a myriad of processes ranging from immune responses to synaptogenesis in the CNS. We eliminated this major cell population and showed remarkable plasticity in the developing auditory brainstem. Not only do microglia completely recolonize the brainstem, but repopulation rectifies previously observed anatomic defects. We showed that microglial presence in MNTB is needed in order for calyceal pruning to take place and for full establishment of the GFAP-positive astrocytic population. Importantly, microglial absence postnatally has a significant impact on auditory function, which partially recovers after microglia return.
Acknowledgments
Acknowledgements: We thank Dr. Adeela Syed at the Optical Biology Core Facility of the Developmental Biology Center at the University of California, Irvine, for her advice and help with confocal imaging and analysis.
Footnotes
The authors declare no competing financial interests.
This work was supported by National Institutes of Health Grants R01 DC010796, R01NS083801 (National Institute of Neurological Disorders and Stroke), R01AG056768 (National Institute on Aging). This study was made possible in part through access to the Optical Biology Core Facility of the Developmental Biology Center, a shared resource supported by the Cancer Center Support Grant CA-62203 and Center for Complex Biological Systems Support Grant GM-076516 at the University of California, Irvine.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.