Myelination Deficits in the Auditory Brainstem of a Mouse Model of Fragile X Syndrome

Auditory symptoms are one of the most frequent sensory issues described in people with Fragile X Syndrome (FXS), the most common genetic form of intellectual disability. However, the mechanisms that lead to these symptoms are under explored. In this study, we examined whether there are defects in myelination in the auditory brainstem circuitry. Specifically, we studied myelinated fibers that terminate in the Calyx of Held, which encode temporally precise sound arrival time, and are some of the most heavily myelinated axons in the brain. We measured anatomical myelination characteristics using coherent anti-stokes Raman spectroscopy (CARS) and electron microscopy (EM) in a FXS mouse model in the medial nucleus of the trapezoid body (MNTB) where the Calyx of Held synapses. We measured number of mature oligodendrocytes (OL) and oligodendrocyte precursor cells (OPCs) to determine if changes in myelination were due to changes in the number of myelinating or immature glial cells. The two microscopy techniques (EM and CARS) showed a decrease in fiber diameter in FXS mice. Additionally, EM results indicated reductions in myelin thickness and axon diameter, and an increase in g-ratio, a measure of structural and functional myelination. Lastly, we showed an increase in both OL and OPCs in MNTB sections of FXS mice suggesting that the myelination phenotype is not due to an overall decrease in number of myelinating OLs. This is the first study to show that a myelination defects in the auditory brainstem that may underly auditory phenotypes in FXS.


INTRODUCTION
Fragile X Syndrome (FXS) is the most common monogenic form of autism spectrum disorder (ASD) and occurs in 1:4,000-1:8,000 people in the United States (Hagerman and Hagerman, 2008). FXS has been commonly used as a model for studying ASD, especially because there are several commercially available animal models such as mouse and rat (The Dutch-Belgian Fragile X Consorthium et al., 1994;Till et al., 2015;Tian et al., 2017). There have been many proposed mechanisms for how loss of Fragile X Mental Retardation Protein (FMRP), the protein encoded by the gene Fmr1, leads to FXS and ASD phenotypes (Bear et al., 2004;Hagerman et al., 2009;Osterweil et al., 2013;Gantois et al., 2017;Rajaratnam et al., 2017, among others). However, drug therapies in animal models do not always "rescue" human phenotypes when extended to the clinic (Dahlhaus, 2018). Recent work has shown that myelination may underly some the phenotypes common in ASD and FXS (Pacey et al., 2013;Phan et al., 2020). Auditory phenotypes, such as auditory hypersensitivity, are often conserved between mouse and human FXS, and it is possible that myelination deficits are also conserved (McCullagh et al., 2020b;Phan et al., 2020). There is reduced or delayed myelination in FXS in many brain areas (Pacey et al., 2013;Phan et al., 2020) and one of the targets of FMRP is myelin basic protein (Jeon et al., 2017), suggesting that deficits seen in FXS may be caused at least in part by alterations to myelination. Additionally, research has shown FMRP expression in mature oligodendrocytes (OLs), potentially explaining how loss of FMRP impacts myelination (Wang et al., 2004;Giampetruzzi et al., 2013). The mechanisms by which binaural hearing is impaired in people with FXS and ASD are unknown, however, there is evidence the auditory brainstem is involved (Kulesza and Mangunay, 2008;Kulesza et al., 2011;Wang et al., 2014;Rotschafer et al., 2015;Garcia-Pino et al., 2017;McCullagh et al., 2017McCullagh et al., , 2020aZorio et al., 2017;Curry et al., 2018;El-Hassar et al., 2019;Lu, 2019).
Previous work has shown that FXS impacts the auditory brainstem and sound localization processing in a myriad of ways, from protein dysfunction to behavioral changes [reviewed in McCullagh et al. (2020b)]. FMRP expression is high in the brainstem, suggesting an important role in brainstem function (Wang et al., 2014;Zorio et al., 2017). Alterations to proteins which directly interact with FMRP can critically disrupt the maintenance of neuronal activity patterns in the brainstem Strumbos et al., 2010;El-Hassar et al., 2019). In addition, previous anatomical and physiological work has shown that excitation/inhibition balance within the auditory brainstem circuit is altered in FXS (Rotschafer et al., 2015;Garcia-Pino et al., 2017;McCullagh et al., 2017;Rotschafer and Cramer, 2017;Curry et al., 2018;Lu, 2019). This circuit relies on precise timing of excitatory and inhibitory inputs for computation of sound location information (reviewed in Grothe et al., 2010). Behavioral studies have shown subtle, but persistent, auditory specific behavioral phenotypes in FXS mice that likely originate in altered auditory brainstem processing (McCullagh et al., 2020a) however, the exact mechanism that underlies these alterations is unknown.
The ability to localize sound is dependent on a very precise comparison of sound input from the two ears, either through interaural timing differences (ITD) or interaural sound intensity/level differences (IID). Processing and interpretation of ITD and IID relies on extremely fast and temporally precise synaptic transmission, leading to some of the most heavily myelinated axons in the brain (Ford et al., 2015). In the auditory brainstem, the medial nucleus of the trapezoid body (MNTB) afferent axons stand out due to their large diameters, which are about ∼3 µm in healthy hearing animals (Ford et al., 2015). The large diameter is due to a high demand on speed and temporal precision by this synapse (Kim J.H. et al., 2013). Any alteration in this pathway that changes the precision and timing of this circuit will lead to substantial difficulties in the ability to separate competing auditory streams and localize sound, as occurs in ASD. Interestingly, increased changes to latency for both behavioral and physiological measures seem to be one of the more repeatable phenotypes in FXS mice McCullagh et al., 2020a, but see Rotschafer et al., 2015;El-Hassar et al., 2019). One way in which changes to myelination in the auditory brainstem would likely manifest is through changes in the speed of electrical propagation, and therefore latency of auditory brainstem timing.
Several studies have shown that myelination deficiencies contribute to ASD phenotypes in several brain areas. This suggests that an area of the brain most critically dependent on myelination for properly encoding sound location information, such as the auditory brainstem, should be especially impacted (Pacey et al., 2013;Phan et al., 2020). We use a combination of microscopy techniques [coherent anti-stokes Raman spectroscopy (CARS) and electron microscopy (EM)] to examine fine myelin microstructure in Fmr1 KO mice compared to wildtype controls (B6). Additionally, we measured the number of mature and precursor oligodendrocytes (OLs and OPCs) as a potential mechanism through which myelin microstructure is altered in FXS mice. We hypothesize that there will be reduced myelin (thickness and diameter) and concomitant changes in number of OLs and OPCs in FXS mice. The goal of this study is to provide a possible mechanism through which auditory phenotypes might arise in FXS.

MATERIALS AND METHODS
All experiments complied with all applicable laws, National Institutes of Health guidelines, and were approved by the University of Colorado Anschutz Institutional Animal Care and Use Committee.

Animals
All experiments were conducted in either C57BL/6J (stock #000664, B6) background or hemizygous male and homozygous Fmr1 knock-out strain (B6.129P2-Fmr1 TM1Cgr /J stock #003025, Fmr1 KO) obtained from The Jackson Laboratory (Bar Harbor, ME, United States) (The Dutch-Belgian Fragile X Consorthium et al., 1994). Animals from both sexes were used in the experiments for both B6 and Fmr1 KO mice. There were no significant differences in any of the measures based on sex, so data for both sexes are combined for analysis (p-values shown are the main effect of sex; CARS p = 0.7454, OL count p = 0.9529, OPC count p = 0.2103, EM data was one animal of each sex so differences do not apply here). Exact number of animals used are listed in the figure legends for each experiment but ranged from 4-7 for each genotype and were between P72-P167 (postnatal day).

Tissue Preparation
Mice were overdosed with pentobarbital (120 mg/kg body weight) and transcardially perfused with phosphate-buffered saline (PBS; 137 mM NaCl, 2.7 mM KCl, 1.76 mM KH 2 PO 4 , 10 mM Na 2 HPO 4 Sigma-Aldrich, St. Louis, MO, United States) followed by 4% paraformaldehyde (PFA). After perfusion, the animals were decapitated, and the brains removed from the skull. Brains were kept overnight in 4% PFA before transferring to PBS. Brainstems were embedded in 4% agarose (in PBS) and sliced coronally using a Vibratome (Leica VT1000s, Nussloch, Germany) at either 200 µm thickness for myelination analysis with CARS or 70 µm thickness for labeling oligodendrocytes.
For EM imaging, mice were perfused with PBS followed by a solution containing 2.5% glutaraldehyde, 2% paraformaldehyde, and 0.1 M cacodylate buffer at pH 7.4. The brains were stored in the same solution for 24 h and sectioned at 500 µm using the same protocol as above. The sections were immersed in glutaraldehyde solution for a minimum of 24 h at 4 • C. The tissue was then processed by rinsing in 100 mM cacodylate buffer followed by immersion in 1% osmium and 1.5% potassium ferrocyanide for 15 min. Next, the tissue was rinsed five times in the cacodylate buffer and immersed again in 1% osmium for 1.5 h. After rinsing five times for 2 min each in cacodylate buffer and two times briefly in water, en bloc staining with 2% uranyl acetate was done for at least 1 h at 4 • C, followed by three rinses in water. The tissue was transferred to graded ethanols (50, 70, 90, and 100%) for 15 min each and then finally to propylene oxide at room temperature, after which it was embedded in LX112 and cured for 48 h at 60 • C in an oven. Ultra-thin parasagittal sections (55 nm) were cut on a Reichert Ultracut S (Leica Microsystems, Wetzlar, Germany) from a small trapezoid positioned over the tissue and were picked up on Formvar-coated slot grids (EMS, Hatfield, PA, United States).

Immunohistochemistry
For staining oligodendrocytes, six to eight free-floating sections from each brain were submerged in L.A.B. solution (Liberate Antibody Binding Solution, Polysciences, Inc., Warrington, PA, United States, Cat No. 24310) for ten min to help expose epitopes related to labeling mature and precursor oligodendrocytes. Next, the slices were washed two times in PBS and blocked in a solution containing 0.3% Triton-X (blocking solution), 5% normal goat serum (NGS) and PBS for 1 h on a laboratory shaker. After blocking, slices were stained with primary antibodies rabbit anti-Aspartoacylase (GeneTex, Irvine, CA, United States; 1:1,000) and mouse Sox-10 (A-2) (Santa Cruz, Dallas, TX, United States; sc-365692, 1:500) diluted in blocking solution with 1% NGS and incubated overnight ( Table 1). Slices were then washed three times (10 min each wash) in PBS and incubated in secondary antibodies ( Table 2) diluted in blocking solution with 1% NGS for 2 h. After three additional washes in PBS (5 min each wash), Nissl (Neurotrace 425/435 Blue-Fluorescent Nissl Stain, Invitrogen, Carlsbad, CA, United States, 1:100) in antibody media [AB media: 0.1 M phosphate buffer (PB: 50 mM KH 2 PO 4 , 150 mM Na 2 HPO 4 ), 150 mM NaCL, 3 mM Triton-X, 1% bovine serum albumin (BSA)] was applied for thirty min. Stained slices were then briefly washed in PBS and slide-mounted with Fluoromount-G (Southern Biotech, Cat.-No.: 0100-01, Birmingham, AL, United States). Slides were stored at 4 • C. Slices used for myelination analysis with CARS (four to five slices per brain) were stained with Nissl as described above, then stored free-floating at 4 • C until imaged. All antibody labeling was performed at room temperature.

Antibody Characterization
The primary antibody for mature oligodendrocytes (Aspartoacylase (ASPA), 1:1,000, GeneTex, Irvine, CA, United States; GTX113389; RRID:AB_10727411, Table 1) is a rabbit polyclonal antibody specific to mature oligodendrocytes (OLs). The ASPA antibody is specific to a recombinant protein encompassing a sequence within the center region of the human aspartoacylase. Aspartoacylase is responsible for the conversion of N-acetyl_L-aspartic acid (NAA) to aspartate and acetate. Hydrolysis by aspartoacylase is thought to help maintain white matter (Bitto et al., 2007). The ASPA antibody has been shown to label mature oligodendrocytes in mouse cerebral cortex and has been validated by protein overexpression and western blot analysis (Larson et al., 2018;Orthmann-Murphy et al., 2020). Additional western blot analyses with various whole cell extracts also showed that ASPA antibody detects aspartoacylase protein. The primary antibody Sox-10 [Sox-10 (A-2), 1:500, Santa Cruz Biotechnology, Heidelberg, Germany; sc-365692; RRID:AB_10844002] was used to label the entire oligodendrocyte lineage [oligodendrocyte precursor cells (OPCs) and OLs]. Sox-10 is a mouse monoclonal antibody specific to an epitope mapping between amino acids 2-29 at the N-terminus of the human Sox-10 gene. Sox-10 has been shown to specifically  label the oligodendrocyte lineage in the CNS and in mouse brains (Zuo et al., 2018;Barak et al., 2019;Imamura et al., 2020). Both primary antibodies were visualized using two fluorescent-conjugated secondary antibodies which are listed in Table 2.

Imaging
Brainstem slides for immunofluorescence were imaged using an Olympus FV1000 confocal microscope (Olympus, Tokyo, Japan) with lasers for 488, 543, and 635 nm imaging. Once MNTB was identified by the distinctive trapezoidal morphology and cell size, z-stacks were taken using a 20x objective (UPLSAPO20X, NA 0.75) so that the entire nucleus could be visualized and quantified. The nucleus was separated into medial, central, and lateral regions following the protocol in Weatherstone et al., 2016. Briefly, the MNTB was digitally extracted using FIJI software, and the tonotopic axis was estimated by drawing the longest possible dorsomedial-to-ventrolateral line through the nucleus (Weatherstone et al., 2016). This line was measured using FIJI and divided into thirds wherein perpendicular lines were drawn. These lines delineated the lateral, central and medial regions of MNTB. For CARS microscopy, brainstem sections were imaged using an Olympus FV-1000 (Olympus Corp., Tokyo, Japan) fitted with non-descanned detectors in both the forward and epi CARS directions. Sections were placed in a culture dish with coverslip (for inverted microscopy) and custom weight to keep tissue near coverslip. Z-stacks were taken with a 60X, 1.2NA infrared corrected water objective which served for collection of the CARS signal in the epi direction for medial and lateral MNTB. Signal in the forward direction was collected through an Olympus 0.55NA condenser. The APE picoEmerald laser set up contained an NKT fiber laser which provided the 80 MHz clock, and an OPO (Optical Parametric Oscillator) laser with a tunable range of 770-990 nm. Water absorption peak was at 1,388 cm −1 which implies a 901.1 nm pump/probe beam with the fixed Stokes beam at 1,031 nm and a CARS signal at 801.6 nm. To excite the CH 2 stretch at 2,845 cm −1 the APE picoEmerald lasers required a pump/probe beam set to 797.2 nm which resulted in a CARS signal at 649.8 nm, thus laser settings were always set to 797.2 nm prior to imaging. Pulse duration for the two lasers was approximately 2 ps and the polarization was linear and horizontal. The synchronization was primarily based on the OPO cavity length with the pulse pumped by the frequency doubled output from the NKT laser at 515.5 nm. The tunable pump/probe beam was used for the TPEF (two-photon excitation fluorescence), which was separated by a dichroic from the CARS signal. The NDD (non-descanned detector) epi-direction unit contained two PMTs of the same type. The CARS and fluorescence do not share PMTs. These settings were used for optimal lipid signal to selectively image myelination in this area.
As an additional control, one brain from a wild type mouse and one brain from a knockout were processed and imaged with electron microscopy. The main purpose of this control was to verify that CARS microscopy has sufficient resolution to discriminate between the small axon diameter differences that are studied here. Sections for EM were first imaged on a compound microscope, and MNTB was identified by the characteristic cell size and shape, as well as the location of the cells in the parasagittal slices. Once MNTB location was determined, slices were imaged on a FEI Tecnai G2 transmission electron microscope (Hillsboro, OR, United States) with an AMT digital camera (Woburn, MA, United States).

Cell Counting
Quantification of OL's and entire OL lineage [oligodendrocyte precursor cells (OPCs) and OLs] was performed using the optical fractionator approach and FIJI software (Schindelin et al., 2012). In a pilot study to determine the appropriate stereological parameters, five brains were sliced at 50 µm using a freezing microtome (SM2010R Leica, Nussloch, Germany) and z-stacks from five sections per brain (every other section) were taken with a 20x objective the same as above. Using FIJI each stack was brightened, the background subtracted, the image scaled and the MNTB was digitally extracted. Separate grid plugins were used to create the dissector and count frames for each image stack. To account for tissue shrinkage, the section thickness was measured at multiple points within MNTB using the Olympus FV1000 microscope and FV10-ASW microscopy software. Measurements were taken by recording the first z-level where tissue features came into focus, and then recording the last z-level where tissue features were in focus. The difference between the recorded levels for each measured area of MNTB was averaged and represented the section thickness for stereological calculations. Stacks were set to contain 20 slices regardless of the section thickness and the step size was kept constant within, but not between brains. Cells were counted by moving through the stack and using a cell counter plugin to mark the top of the cell when it first came into focus within the probe unless it contacted the exclusion lines of the count frame (Schmitz and Hof, 2005).
The pilot data indicated that selecting four 70 µm sections (every other section starting from a randomly selected slice) to give a section subfraction (SSF) of 0.280 was appropriate for calculating the total number of cells in a single MNTB. In addition, counts would use a probe area of (1,250 µm × 1,250 µm) and dissector (5,000 µm × 5,000 µm) giving an area subfraction (ASF) of 0.25. It should be noted that two brains had fewer than four slices (3.5 and three slices, respectively) which were suitable for quantification. OL analyses also included three brains from the pilot study that had adequate immunohistochemical staining for quantification. In these cases, the SSF was adjusted, or a fifth slice was counted to achieve an SSF that was closest to 0.280 (0.20-0.25). The height subfraction (HSF) was calculated separately for each brain using the mean dissector height (OL mean = 19.5-34.4 µm; OPC mean = 21.06-34.4 µm) and the mean section thickness (OL mean = 26.75-45.75 µm; OPC mean = 27.75-45.75 µm). Using these parameters, OL's and OPC's were counted and recorded in excel. OPC's were identified as cells that were labeled with Sox 10, but not co-labeled with ASPA (Figure 1). Calculation of the total cell count for each complete MNTB nucleus was done by multiplying the counted cells ( Q − ) with the reciprocal sampling fractions such that: Error in the stereological approach was calculated using the Gundersen-Jensen coefficient of error (CE) estimator and smoothness class of m = 1 (Gundersen et al., 1999). CE measures error as a function of the biological variance (noise) and the sampling variance from systemic uniformly random sampling. The coefficient of variance (CV) within and between brains was high (CV > 1) due to natural variations in the population of glial cells throughout MNTB. CE for OL and OPC estimates was greater than 10% (Mean CE = 0.51), but the sampling parameters chosen showed CE was negligible in terms of overall variation (CE 2 /CV 2 < 0.003) (Gundersen et al., 1999). Thus, the stereological parameters were determined to be sufficient for quantification estimates.

Myelination Analysis
Analysis of myelination morphology was done using FIJI software. Myelin diameter was measured from CARS images by using the line tool to measure the fiber diameter (Figure 2A) for 20 axons per image. The resolution was not sufficient in CARS images to quantify thickness or g-ratio of axons. Tissue shrinkage was not corrected for in this analysis. In addition, diameters of FIGURE 2 | Quantification of myelination morphology in CARS and EM images. CARS images were quantified using the FIJI line segment tool. A line, as indicated by white line (A), was drawn and the measure function used to measure the fiber diameter for 20 axons per image (purple indicates cell bodies stained with Nissl, CARS signal in cyan). For EM images, increased resolution allowed for more detailed analysis including fiber diameter (dark blue), axon diameter (cyan), myelin thickness (green), inner (orange) and outer (purple) fiber diameter (B). Inner and outer fiber perimeter were then used to calculate g-ratio, which is the inner fiber diameter/outer fiber perimeter.
Frontiers in Neuroscience | www.frontiersin.org 240 axons were also measured in the same way in EM images, however, myelin thickness was measured directly, and g-ratios were calculated from fiber and axon diameters ( Figure 2B; Ford et al., 2015). G-ratios were calculated from EM images using inner and outer fiber perimeter ratio (g-ratio = inner fiber perimeter /outer fiber perimeter). EM data is representative of one animal from each genotype and was performed primarily to confirm CARS diameter changes.

Statistical Analyses
Figures were generated in R (R Core Team, 2013) using ggplot2 (Wickham, 2016). Boxplots dislay the median and 25th-75th percentiles (or 1st and 3rd quartiles, respectively) the whiskers represent ± 1.5 times the interquartile range. Data that falls outside the range are plotted as individual points. Data for all experiments were analyzed using linear mixed-effects models [lme4; (Bates et al., 2014)]. To account for repeated measures and variability with animals, Animal was considered a random effect with fixed effects of genotype and location (medial, center, lateral, and total when appropriate). Additional analyses used fixed effects of age or sex to determine if these were contributing factors in the differences shown. It was expected that there may be differences between different regions of the MNTB (medial, center, lateral, total) therefore a priori, it was determined that estimated marginal means [emmeans; (Lenth, 2019)] would be used for pairwise comparisons between region and genotype.
To control for multiple comparisons, emmeans implements a Tukey method for these contrasts. Asterisks are used to indicate statistical significance between the two genotypes, as follows: * p < 0.05, * * p < 0.01, and * * * p < 0.001. Figures were prepared for publication using Photoshop and Illustrator (Adobe, San Jose, CA, United States).

RESULTS
We used several (EM and CARS) microscopy techniques to examine the myelin microstructure of putative globular bushy cell axons in the auditory brainstem that innervate the MNTB and form the Calyx of Held in Fmr1 KO mice compared to wildtype controls (B6). We also counted the number of both mature and precursor oligodendrocytes in the same region to show if structural changes are due to changes in number or type of oligodendrocytes which myelinate these axons.

Myelin Diameter as Measured by Coherent Anti-stokes Raman Spectroscopy
We first examined the diameter of axons in the medial and lateral MNTB using CARS microscopy (Figure 3). CARS microscopy for myelination imaging was used over electron microscopy (EM) due to its speed and compatibility with immunofluorescence (Wang et al., 2005), in this case we used a Nissl stain for simultaneous cell body staining. We measured axon fiber diameters in both medial and lateral MNTB because previous work has shown the potential for differences in myelin based on tonotopic location (medial or lateral) (Ford et al., 2015;Stange-Marten et al., 2017). Representative images for both location (medial or lateral columns) for B6 and Fmr1 KO (rows) are shown in Figures 3A-D. We found that there was a significant decrease in fiber diameter in Fmr1 KO compared to B6 for both medial and lateral MNTB ( Figure 3E). Not surprisingly when data were summed between medial and lateral locations, there was still a significant decrease in fiber diameter in Fmr1 KO mice.

Coherent Anti-stokes Raman Spectroscopy Myelin Characteristics Confirmed by Electron Microscopy
To confirm that the CARS technology was sensitive enough to determine differences between fiber diameters in Fmr1 KO and B6 mice, we also performed EM controls on one animal of each genotype (across multiple sections with multiple measurements, Figure 4). The advantage of EM, while costly and time consuming, is that it allows for higher resolution and thus more detailed description of changes to myelination than CARS microscopy techniques. Since alterations in fiber diameter seen in CARS data was not dependent on medial or lateral localization, and it is difficult to determine location with parasagittal EM sections, these data are not separated by location. Representative EM images for B6 and Fmr1 KO are shown in Figures 4A,B, respectively. Consistent with CARS analysis, EM data showed smaller fiber diameters in Fmr1 KO animals compared to B6 ( Figure 4C). Additionally, axon diameter and myelin thickness are also decreased in Fmr1 KO compared to B6 (Figures 4D,E). Lastly, g-ratio, an often-used index of optimal functional and structural myelination, was increased in Fmr1 KO animals compared to B6 (Figure 4F; Chomiak and Hu, 2009). These results are consistent with the relatively large axon diameters in the Fmr1 KO mice compared to fiber diameters (a measure for g-ratio) despite thinner myelin and smaller axons than B6 mice.

Number of Mature Oligodendrocytes
Oligodendrocytes (OLs) are the glial cells responsible for myelinating axons in the central nervous system. FMRP has been shown to be expressed in OLs, potentially providing a mechanism through which loss of FMRP can lead to changes in myelination (Wang et al., 2004;Giampetruzzi et al., 2013). For OL counts, the MNTB was divided into three sections, medial, central, and lateral. Figures 5A,B show representative images of the entire MNTB stained with ASPA (a marker for mature OLs) and sox-10 (a marker for the entire OL lineage) for B6 and Fmr1 KO mice, respectively. As shown by arrowheads, OLs are represented by staining with ASPA (which also coincides with sox-10 staining). Cells that are only stained for sox-10 (magenta) and not ASPA (yellow) are quantified as OPCs (see below and Figures 1, 5). Like CARS data, there was no difference in OL number based on specific tonotopy of the MNTB (medial, center, or lateral) however, when data was combined for all regions, there was a significant increase in OL number in Fmr1 KO mice compared to B6 (Figure 5C top panel). Using the sampling Q count technique, to estimate number of cells within the entire MNTB, similar results were found to raw counts.  Specifically, there was no difference in medial, center, or lateral estimates, but total combined regions showed a highly significant (p < 0.0001) increase in OLs in Fmr1 KO animals compared to B6. In addition to changes in number of mature OLs, another possible mechanism underlying alterations to myelination could be incomplete or impaired maturation of OLs.

Number of Precursor Oligodendrocytes
Other brain areas (cerebellum, neocortex, and others) show changes in myelination in FXS related to impaired maturation or function of oligodendrocyte precursor cells (Wang et al., 2005;Pacey et al., 2013;Jeon et al., 2017). Like OL count measures, the MNTB was divided into medial, central, and lateral subsections ( Figure 5D), with representative images of the entire MNTB shown in Figure 5 for B6 ( Figure 5B) and Fmr1 KO ( Figure 5A). As described above, OPCs were quantified as cells that were positive for the sox-10 (magenta) marker without ASPA (yellow) staining ( Figures 5A,B arrowheads, OPC label). Like OL counts, there was no significant difference between Fmr1 KO and B6 mice based on tonotopic location (medial, center, or lateral, Figure 5D), however, when data was summed for a total MNTB count, there was a significant increase in OPC number in Fmr1 KO mice compared to B6 (Figure 5D, top panel). Consistent with data from OLs, whole nucleus Q counts were consistent with raw counts such that there was no difference between genotypes based on tonotopy, but total estimates were significantly higher in Fmr1 KO animals compared to B6 (p < 0.0014). Together, these OPC and OL data show that changes in myelination observed through EM and CARS are not due to inherent decreases in numbers of mature or immature oligodendrocytes, however, there are changes in overall number of these cells compared to wildtype animals.

DISCUSSION
Mechanisms underlying auditory symptoms among patients with FXS are still poorly understood. The auditory brainstem is the first location along the ascending auditory pathway where sound information is processed, and it is likely key to identifying auditory difficulties since any changes at this level alter midbrain and cortical processing as well. Here we show decreased myelination in adult FXS mice (fiber diameter, axon diameter, thickness, and increased g-ratio) compared to wildtype in the region of the MNTB, which contains one of the largest and most heavily myelinated axons in the central nervous system. One possible mechanism for impaired myelination is changes in oligodendrocyte number and type. We show an increase in mature and precursor oligodendrocytes, suggesting perhaps an overcompensation for earlier reductions in OL numbers that others have shown (Pacey et al., 2013). Results from these measures suggest impaired myelination in the MNTB of Fmr1 KO mice, which could be a potential mechanism for auditory brainstem physiological and behavioral phenotypes. While this is the first study to examine myelination changes in the auditory brainstem in FXS, previous work has shown reduced myelin sheath growth, decreases in total myelin (Doll et al., 2020), fewer myelinated and thinner axons, reductions in myelin basic protein, and development impacts of loss of FMRP (Pacey et al., 2013). FMRP is present in both precursor and mature oligodendrocytes (Wang et al., 2004;Giampetruzzi et al., 2013) and located subcellularly within myelin sheaths (Doll et al., 2020). However, recent work suggests that FMRP is not acting directly to regulate myelin basic protein (MBP) or other myelin protein transcripts, suggesting a diversity of roles for FMRP in myelin and oligodendrocyte dynamics (Giampetruzzi et al., 2013;Doll et al., 2020). Interestingly, inconsistent with the literature, and unexpectedly based on reductions in myelin characteristics measured here, we saw an increase in number of both OLs and OPCs. These differences can potentially be explained by an overcompensation for OL number reductions during development in this brain area, or an increased number of OLs that are not necessarily myelinating efficiently (fewer myelin sheaths per OL, non-myelinating OLs, etc.). In addition, it is possible that there are changes in proliferation ability of OPCs and OLs that were not measured here with only single time point measurements (Bu et al., 2004). Lastly, we did not measure a direct effect of FMRP on myelination, and amount of myelination in the MNTB also depends on neuronal activity (Sinclair, 2017), we cannot rule out the possibility of loss of FMRP having secondary effects besides an impact on OLs. What consequences these reductions in myelination have on conduction speed and physiological properties of the Calyx of Held in synaptic transmission are potential areas for future research.
G-ratio is a common measure for not just structural but potential physiological properties of myelin conduction along an axon. The g-ratios we measured for Fmr1 KO and B6 animals are similar to previous work (particularly for wildtype B6 animals) in the MNTB (Ford et al., 2015;Stange-Marten et al., 2017). Conduction velocity is an important factor in reliable processing of sound location information, resulting in the Calyx of Held which is highly specialized for speed and temporal fidelity. Increased g-ratio suggests slower conduction velocities and thinner myelin in Fmr1 KO mice (Rushton, 1951;Berman et al., 2019). However, note that there may be other factors which contribute to conduction velocity, for example Ranvier node and internode spacing, length, and diameter (Ford et al., 2015). In addition, while we did not see differences in myelin properties based on location or tonotopic area, it is possible that node dynamics vary more with tonotopic area while myelin thickness and diameter are less impacted based on tonotopy. However, note that such differences may also be explained by the relative importance of ILD versus ITD processing in different animal models. The previous work showing tonotopic differences was performed in gerbils, a rodent model suitable to study both low and high-frequency hearing. In that animal model, one important role for MNTB neurons is to provide inhibitory input to ITD processing (Ford et al., 2015). By contrast, mice are primarily high-frequency specialists. Consistent with their hearing spectrum and the fact that mice primarily use ILD cues, previous work has shown that there are no differences in myelination properties based on tonotopy (Stange-Marten et al., 2017), similar to what we see in this study.
In the current study we did not directly label Calyx of Held axons. Rather, we made measurements within the MNTB, measuring axons that were projecting in the coronal plane. Although axons projecting to MNTB are significantly larger than any other passing fibers and can be discerned by eye relatively easily and reliably, we cannot rule out the possibility that some measured axons may be "passing through." If our analysis did include some passing fibers, these would likely decrease the observed differences between wild type and mutant, such that the results presented in this study may be a lower limit of the true differences between wild type and mutant mouse model. This would also limit the interpretation of conclusions about tonotopic location of the measurements and whether indeed there are differences based on frequency-coding. Indeed, we have no way of discerning if the medial MNTB measurements included lateral axons passing through to the lateral MNTB. Note that CARS nor EM techniques inherently discriminate between different types of axons. However, based on the size of the axons that we measured in both EM and CARS, fibers are consistent with the expected size of Calyx axons in the same area (Ford et al., 2015;Sinclair, 2017;Stange-Marten et al., 2017), therefore we are confident that at least the large majority of the axons are Calyceal projections. This is the first study to show myelination changes in the auditory brainstem sound localization pathway in FXS mice. Thinner and smaller diameter Calyx axons with increased g-ratio may underly sound localization and auditory hypersensitivity issues seen behaviorally in mice and described by patient's with FXS. Interestingly, we saw an increase in both mature and immature OLs suggesting that myelin deficits are not due simply to fewer myelinating OLs. Future studies elucidating when during development myelin deficits begin and how OLs develop from precursors into mature OLs in this brain area would be crucial to understanding the complete picture of auditory brainstem phenotypes in FXS.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The animal study was reviewed and approved by University of Colorado Anschutz Institutional Animal Care and Use Committee.

AUTHOR CONTRIBUTIONS
EM, SP, AL, and AK developed the ideas and methods. EM, SP, and AL collected the data for the manuscript. EM and AL created the statistical analyses and figures for the manuscript. All authors helped wrote and revised the manuscript.

FUNDING
Supported by NIH R01 DC 17924, R01 DC 18401 (Klug), and NIH 1R15HD105231-01, T32DC012280, and FRAXA (McCullagh). The CARS imaging was performed in the Advanced Light Microscopy Core part of the NeuroTechnology Center at the University of Colorado Anschutz Medical Campus supported in part by NIH P30 NS048154 and NIH P30 DK116073.