Abstract
The lateral cortex of the inferior colliculus (LCIC) is a multimodal subdivision of the midbrain inferior colliculus (IC) that plays a key role in sensory integration. The LCIC is compartmentally-organized, exhibiting a series of discontinuous patches or modules surrounded by an extramodular matrix. In adult mice, somatosensory afferents target LCIC modular zones, while auditory afferents terminate throughout the encompassing matrix. Recently, we defined an early LCIC critical period (birth: postnatal day 0 to P12) based upon the concurrent emergence of its neurochemical compartments (modules: glutamic acid decarboxylase, GAD+; matrix: calretinin, CR+), matching Eph-ephrin guidance patterns, and specificity of auditory inputs for its matrix. Currently lacking are analogous experiments that address somatosensory afferent shaping and the construction of discrete LCIC multisensory maps. Combining living slice tract-tracing and immunocytochemical approaches in a developmental series of GAD67-GFP knock-in mice, the present study characterizes: (1) the targeting of somatosensory terminals for emerging LCIC modular fields; and (2) the relative separation of somatosensory and auditory inputs over the course of its established critical period. Results indicate a similar time course and progression of LCIC projection shaping for both somatosensory (corticocollicular) and auditory (intracollicular) inputs. While somewhat sparse and intermingling at birth, modality-specific projection patterns soon emerge (P4–P8), coincident with peak guidance expression and the appearance of LCIC compartments. By P12, an adult-like arrangement is in place, with fully segregated multimodal afferent arrays. Quantitative measures confirm increasingly distinct input maps, exhibiting less projection overlap with age. Potential mechanisms whereby multisensory LCIC afferent systems recognize and interface with its emerging modular-matrix framework are discussed.
Introduction
An essential role of the nervous system is to encode and integrate various sensory stimuli. Such processing requires highly-ordered network configurations that coordinate converging inputs from distinct systems. The inferior colliculus (IC) is a strategically positioned midbrain hub that receives an array of top-down and bottom-up afferents of multisensory origin. Within its lateral cortex (LCIC) single units respond to a variety of stimuli (Aitkin et al., , ; Gruters and Groh, ), suggesting an underlying neuroanatomical substrate capable of multisensory processing. In adult mice, input channels target discrete regions of the LCIC in a modality-specific manner (Lesicko et al., ). Somatosensory projections terminate as a series of discontinuous patches or modules that span LCIC layer 2. In contrast, auditory inputs exhibit complementary patterns, preferentially targeting the surrounding LCIC matrix (layers 1, 3, and intermodular zones). A recent study examining intrinsic LCIC circuits shows that local connectivity is largely confined to its own compartment (module or matrix), although unidirectional flow from matrix regions that receive auditory inputs into somatosensory-rich modular zones also occurs (Lesicko et al., ). Such intercompartmental connectivity likely underlies its reported multisensory capabilities that in turn inform precise response behaviors.
To date little is known about the development of multimodal LCIC afferent systems. Previous work from our lab defined an early postnatal LCIC critical period (birth: postnatal day 0, P0 through P12) based on the emergence of its characteristic compartmental framework (Chernock et al., ). A host of neurochemical stains reveal that while not apparent at birth, a discrete microarchitecture quickly emerges (Dillingham et al., ). Most noteworthy from the identified markers are GAD (glutamic acid decarboxylase) which reliably highlights LCIC modules, and calretinin (CR) which labels the extramodular matrix. Temporally correlated with the emergence of LCIC neurochemical compartments is the transient expression of matching Eph-ephrin guidance patterns (Gay et al., ; Stinson et al., 2021). Expression of both EphA4 and ephrin-B2 align with its developing modularity, while ephrin-B3 expression is complementary and restricted to the surrounding matrix.
We hypothesized that an early period of projection shaping exists whereby somatosensory and auditory afferents target appropriate LCIC compartments. Recently, we characterized the development of auditory inputs to the LCIC arising from its neighboring central nucleus (CNIC, Lamb-Echegaray et al., ). While initially diffuse at birth, an early projection specificity for the LCIC matrix soon develops. This mapping preference of LCIC auditory afferents appears adult-like by P12, thus sharing a similar critical period defined for its emerging micro-organization and matching Eph-ephrin guidance patterns. Here, we combine anterograde tract-tracing approaches in living preparations of developmental GAD67-GFP tissue with immunocytochemical methods to: (1) determine whether somatosensory inputs follow a similar developmental progression in targeting LCIC modular zones; and (2) visualize and quantify the segregation of two multimodal LCIC afferent patterns (somatosensory: corticocollicular; auditory: from the CNIC) with respect to each other. Results indicate a model in which initially overlapping projection distributions refine and segregate into modality-specific compartments during an early postnatal window. Potential mechanisms that may instruct discretely-organized LCIC multisensory maps are discussed.
Materials and Methods
Animals
Experiments were performed on neonatal GAD67-GFP mice (P0, P4, P8, and P12; n = 36, consisting of at least three mice at each age for both single- and dual tracer studies). The GAD67-GFP knock-in line has previously been validated in our lab (Gay et al., ) and allows for easy visualization of GAD-positive LCIC modules (Lesicko et al., , ; Gay et al., ; Lamb-Echegaray et al., ; Stinson et al., 2021; Brett et al., ). Specifics concerning the generation of the GAD67-GFP (Δneo) line are described elsewhere (Tamamaki et al., 2003; permission granted by Dr. Yuchio Yanagawa, Gunma University Graduate School of Medicine, Gunma, Japan). Heterozygous GAD-GFP males were crossed with C57BL/6J females and GFP-expressing progeny were identified before P4 using a Dark Reader Spot Lamp Goggle visualization system (Clare Chemical Research1, Dolores, CO, Cat# SL10S). Equal numbers of males and females were used in the described experimentation and no sex-specific differences were noted. All procedures were performed in keeping with the US National Research Council’s Guide for the Care and Use of Laboratory Animals and received prior approval by the Institutional Animal Care and Use Committee (Protocol No. 20-1421).
Anterograde Tracing in Living Preparations
Following an overdose of ketamine (200 mg/kg) and xylazine (20 mg/kg), mice at designated stages were perfused with chilled, oxygenated (95% O2, 5% CO2) artificial cerebrospinal fluid (aCSF; in mM; 126 NaCl, 3 KCl, 1.25 NaH2PO4, 10 dextrose, 20 NaHCO3, 1.2 MgSO4, 2.5 CaCl2, pH 7.4). Brains were blocked in the coronal plane (just rostral to somatosensory cortex and just caudal to the midbrain) and varied in total thickness depending upon age. Under a dissecting microscope, gross placements of biocytin crystal (B1758, Sigma-Aldrich, St. Louis, MO) were made in somatosensory cortex at varying depths depending on age to maximize labeling of layer 5/6 corticocollicular projections. In double-labeling studies, the cerebellum was then removed to facilitate deposition of a 10,000 MW dextran (D22904, AlexaFluor 647 direct conjugate, ThermoFisherScientic, Waltham, MA) in the CNIC ipsilateral to the cortical biocytin placement. Tissue blocks were bubbled at room temperature in aCSF for 18–24 h to allow for complete filling of terminal endings in the LCIC. Prior to sectioning, tissue was postfixed and cryoprotected in a series of 4% paraformaldehyde solutions with increasing concentrations (10%–30%) of sucrose.
Tissue Processing and Immunocytochemistry
A rostrocaudal series of sections were cut at 50 μm on a sliding freezing microtome and collected in 0.1 M phosphate buffered saline (PBS, pH 7.4). Free-floating sections were rinsed three times for 5 min in PBS. In experiments combining tract-tracing with a matrix marker, a 2.5 h DyLight 549 streptavidin step (1:200, SA-5549, Vector Laboratories, Burlingame, CA, RRID:AB_2336408) for visualizing biocytin always preceded subsequent CR immunostaining. Following another round of PBS rinses, sections were blocked in 5% normal donkey serum (NDS) in PBS for 30 min. Tissue was then incubated in anti-CR primary made in rabbit (1:250, CR 7697, Swant, Burgdorf, Switzerland, RRID:AB_2619710) for 40 min at room temperature and then overnight at 4°C. An Alexa Fluor 350 donkey anti-rabbit IgG (1:25, A10039, Thermo Fisher Scientific, Waltham, MA, RRID:AB_2534015) was then applied prior to a final series of PBS rinses. Sections were mounted on charged slides, coverslipped with Pro Long Diamond (P36970, Thermo Fischer Scientific, Waltham, MA), and stored in the dark until imaging.
Image Acquisition and Quantitative Methods
Aside from documentation of somatosensory cortex tracer placements, data collection focused exclusively on the ipsilateral IC, as corticofugal projections to the contralateral LCIC are sparse. Verification of accurate cortical and CNIC tracer deposits were made using age-matched sections in The Atlas of the Developing Mouse Brain (Paxinos et al., 2007). Sections throughout the rostrocaudal extent of the IC were imaged using epiflourescent microscopy (Nikon Eclipse Ti-2 microscope equipped with a monochrome, Hamamatsu ORCA-Flash 4.0 V3 sCMOS camera, Japan; PlanApo objectives 10× : NA = 0.30, 20× : NA = 0.75, and 40× : NA = 1.30). Filter sets (Chroma Technology, Bellows Falls, VT) were designed with careful attention to the spectra of the various fluorophores to ensure no cross-channel bleed through (Alexa Fluor 350 filter set: excitation 381–403 nm emission 417–477; GFP filter set: excitation 446–486 emission 500–550; biocytin-streptavidin DyLight 549 filter set: excitation 542–566 emission 582–636; dextran AlexaFluor 647 filter set: excitation 593–643 emission 663–733. An extended depth of field (EDF) algorithm was used to generate two-dimensional images from acquired Z-stacks (Elements Software; Nikon) using only focused portions of each optical slice. Monochrome channel acquisitions were pseudocolored accordingly (blue: CR, green: GAD, biocytin: red, dextran: cyan) and saved as lossless JPEG2000 files.
Quantification of somatosensory afferent pattern alignment with emerging LCIC compartments, as well as assessments of multisensory pattern segregation, focused on mid-rostrocaudal regions of the LCIC where the modular-matrix framework was most readily apparent. Data from a minimum of three mice were analyzed at each age and multiple sections (a minimum of two) were sampled to account for any potential regional LCIC variability. LCIC layer 2 sampling and generation of brightness plot profiles has been previously described in detail (Wallace et al., 2016; Dillingham et al., ; Gay et al., ; Lamb-Echegaray et al., ; Stinson et al., 2021; Brett et al., ). In brief, separate image channels were converted to grayscale and exported as uncompressed TIFFS prior to importing into ImageJ2 (NIH, Bethesda, MD, RRID:SCR_003070). A freehand tool (line thickness 100 μm) was used to sample LCIC layer 2 from ventral-to-dorsal bisecting GAD-positive modules. Sampling contours were approximated through presumptive LCIC layer 2 in P0 mice, as GAD-defined modules are not easily discerned at this earliest postnatal time point. An ROI function was employed to apply the same sampling contour to each subsequent channel sampled of the same image. Generated signals (single tracer studies: somatosensory vs. GAD/CR, dual tracer studies: somatosensory vs. auditory) were compiled as multi-channel brightness profiles to visualize relative signal overlap/non-overlap. Raw data values for each waveform were exported for cross-correlation and kurtosis analyses described below. The latter analysis used FS and FA to refer to sampled somatosensory and auditory fluorescence values, respectively.
Cross-correlation analyses (Microsoft Excel, Redmond, WA) were performed to numerically quantify the similarity of two signal patterns relative to each other (e.g., somatosensory and GAD, or somatosensory and CR). Cross-correlation function y-intercepts (i.e., no spatial shift or relative displacement) serve as an objective measure for assessing signal matching/mismatch. Values range from +1.0 to −1.0, with higher values indicating strong signal overlap, whereas lower values indicate waveforms that are non-overlapping or out-of-phase. Single factor ANOVA assessed similarity among developmental stages. Based on ANOVA high confidence dissimilarity (p < 0.05), subsequent comparisons between age pairs were made using independent, two-tailed Student’s t-tests assuming unequal variance (Type 3) with statistical significance (p < 0.05). Similar t-tests (albeit one-tailed) were performed comparing y-intercepts of somatosensory vs. GAD and somatosensory vs. CR signal pairings.
To quantify the extent of segregation of developing multimodal LCIC inputs, we build upon established unbiased methods applied to emerging retinogeniculate patterns (Torborg and Feller, 2004; Jaubert-Miazza et al., ; Torborg et al., 2005). These previous methods processed afferent data in four main stages: (1) a “rolling ball” algorithm correcting for uneven background; (2) logarithmic ratios (R values) of ipsilateral and contralateral fluorescence intensities, FI and FC respectively; R = log10(FI/FC); (3) R value histograms revealing the extent of segregation in their shape; and finally (4) variances of R-distributions quantifying the extent of segregation in a single metric. Typically, overlapping patterns produce unimodal R distributions, while decreased projection overlap (or segregation) produce bimodal R distributions. Testing for overlap/segregation therefore requires quantifying the shape of R distributions.
Accordingly, we also perform the same first three stages, but in the last stage employ a different R distribution statistic (kurtosis rather than variance) to quantify R distribution shape to assess relative input separation. Variance is the second central moment (i.e., average squared deviations from the mean) and a measure of distribution spread. Kurtosis is the fourth central moment (i.e., average fourth power deviations from the mean; normalized by variance) and a direct measure of distribution shape. Unimodal distributions yield more positive kurtosis values, while bimodal distributions yield more negative kurtosis values. In sum, we also process our dual tract-tracing data in four stages: (1) applying a “rolling ball” filter in ImageJ (diameter = 20 pixels) correcting for uneven background fluorescence (i.e., regions of the LCIC exhibiting no axonal labeling); (2) logarithmic ratios (R values) combining fluorescence intensities of somatosensory and auditory inputs, FS and FA respectively; R = log10(FS/FA); (3) histograms of R values showing extent of segregation in their shape; and (4) kurtoses of R-distributions to quantify extent of segregation in single statistics.
Results
Development and Quantification of Descending Somatosensory Projections Targeting LCIC Modules
To determine whether somatosensory inputs to the LCIC share the same critical period defined for auditory inputs (Lamb-Echegaray et al., ) and its emerging compartments, unilateral biocytin crystal placements were made in somatosensory cortex (Figure 1) in a developmental series of GAD67-GFP living preparations. Given the sparse nature of the contralateral corticocollicular pathway, analyses focused exclusively on the LCIC ipsilateral to the tracer placement. Terminals arising from somatosensory cortex were present in the LCIC at birth and diffusely distributed (Figure 2A), lacking any regional specificity. GAD-positive modules at this age are not readily apparent (Figure 2B), in keeping with previous reports (Dillingham et al., ; Gay et al., ; Lamb-Echegaray et al., ; Stinson et al., 2021; Brett et al., ). As layer 2 modular fields emerge at P4 and continue to develop through P12, a preference of somatosensory endings for these defined zones becomes increasingly apparent (Figures 2C–H). By P12, projection specificity appears adult-like, with elaborate arborizations within modular confines and little evidence of significant labeling in the surrounding matrix.
Figure 1
Qualitative observations of increased matching of somatosensory and GAD patterns with age (i.e., in-phase signal waveform oscillations, arrows in Figures 2D,F,H) were confirmed with cross-correlation analyses. The positively sloped linear regression of y-intercept values plotted by developmental stage supports the notion that initially dispersed somatosensory inputs become increasingly specific for LCIC modules (Figure 3A). Single factor ANOVA showed significant mean dissimilarity among age groups (p < 0.001). Thus, t-tests comparing age pairs were performed that showed cross-correlation y-intercept values statistically different for P0 vs. P4 (p = 0.019), P0 vs. P8 (p = 0.009), and P0 vs. P12 (p = 0.005), but not other age pairs (p > 0.05). Off-origin cross-correlation function maxima and minima occur at certain relative spatial shifts between the two series, signifying aligned and mis-aligned signals, respectively. A representative P12 cross-correlation function (Figure 3B) shows a distinctly positive y-intercept (i.e., strong signal matching with no spatial shift) and pronounced off-origin peaks and troughs that reflect the periodic nature of the LCIC modular-matrix framework at this age. Taken together, these data suggest a similar period of projection shaping for top-down somatosensory inputs that gradually sharpen and overlap emerging LCIC modular domains.
Figure 2
Figure 3
To further confirm the establishment of highly refined somatosensory inputs for LCIC modules by its critical period closure (P12), analogous labeling studies were performed coupled with immunostaining for the established matrix marker, calretinin. As anticipated, concentrated somatosensory terminal fields overlapping layer 2 GAD cell clusters filled voids in the encompassing CR-positive matrix (Figure 4, arrowheads). This complementary patterning was reliably observed in additional P12 mice (Figure 5A). Higher magnification reveals dense somatosensory terminal fields sculpted such that they nearly exclusively reside within modular confines (Figures 5B,C, dashed contours). Three-channel LCIC layer 2 sampling and corresponding brightness plot profiles illustrate oscillatory waveforms, with matching somatosensory and GAD signals that are mismatched with CR fluctuations (Figures 5D,E, arrows).
Figure 4
Figure 5
Comparisons of cross-correlation y-intercept values for somatosensory/GAD vs. somatosensory/CR signal combinations further demonstrate discrete afferent patterns at P12 that target LCIC modules (Figure 6). Tightly clustered positive y-intercepts for somatosensory/GAD labeling indicate reliable signal matching or registry. In contrast, more negative cross-correlation y-intercept values for the same axonal distributions with respect to CR matrix labeling confirm non-overlapping patterns or signal mismatch (Figure 6A). Such quantitative findings support our qualitative observations (Figures 4, 5) that somatosensory terminals spatially align with modular domains and are offset from the encompassing matrix. The y-intercept median and standard deviation for biocytin and GAD signals were + 0.45 and 0.12, respectively, whereas that for biocytin and CR signals were—0.07 and 0.22. Not surprisingly given these dissimilar medians and overall distributions, these two data sets proved to be statistically different from each other (p < 0.001). A representative P12 cross-correlation function for biocytin and GAD shows a positive y-intercept (i.e., strong overlap at zero shift) with an off-origin trough and peak indicative of strong signal mismatching and matching respectively at increasing relative shifts (Figure 6B). A contrasting P12 cross-correlation function for biocytin and CR exhibits a negative y-intercept (i.e., non-overlap at zero shift) with an off-origin peak and trough indicative of subsequent strong signal matching and mismatching with increasing relative shifts (Figure 6C).
Figure 6
Early Segregation of Multisensory Input Maps Into Discrete LCIC Compartments
These findings, taken together with previous work from our lab characterizing developing auditory patterns (Lamb-Echegaray et al., ), suggest multimodal LCIC inputs share a similar critical period of shaping that yields complementary projection arrangements. Furthermore, they point towards emergence of discrete maps from initially intermingling distributions that undergo a process of segregation. To directly test such notions, we performed dual tracer studies simultaneously labeling somatosensory (corticocollicular) and auditory (arising from the CNIC) inputs in GAD67-GFP mice (Figure 7A) over the defined critical period (P0–P12). In each case, coronal sections containing placement centers for both dyes were documented and matched with corresponding atlas plates to verify tracer location (Figures 7B,B’,C,C’; plates modified for simplicity from Paxinos et al., 2007). As in the single-labeling experiments, placements centered in somatosensory cortex were made sufficiently deep to include corticocollicular cells of origin located in layers 5/6 (Figures 7B,B’, see Slater et al., 2013; Stebbings et al., 2014). CNIC placements were appropriately large to fill this subdivision, while not bleeding into the neighboring LCIC or DCIC (Figures 7C,C’).
Figure 7
At birth, both somatosensory and auditory terminals are present in the LCIC, albeit sparse and unorganized (Figures 8A–D). Results at P4 were somewhat variable, with some cases still exhibiting intermingling patterns and inconclusive modules (Figures 8E–H), while others show clear indications of an emerging modularity and initial projection separation. These findings suggest P4 as a pivotal timepoint for the developing LCIC compartmental framework and its interfacing multimodal afferents. By P8, GAD-positive modules are easily discerned and projection segregation is in progress, as somatosensory inputs preferentially terminate within modular zones and auditory inputs the surrounding matrix (Figures 8I–L). By P12, the process of segregation is largely complete, with adult-like projection patterns that terminate almost exclusively within their respective compartments (Figures 8M–P, Supplementary Figure 1). Sampling of tracer channels at P12 consistently generate periodic waveforms with out-of-phase signal fluctuations, signifying fully segregated multimodal patterns (Figure 9).
Figure 8
Figure 9
Quantification Reveals Progressive Input Segregation Over Critical Period
The kurtosis statistic was used to quantify the shape of R distributions which were calculated from somatosensory/auditory fluorescence ratios. Generally higher kurtosis values at earlier ages, contrasting with negative kurtosis values at later ages, indicate gradual changes in distribution shape from unimodal to bimodal with increasing age (Figure 10). Since overlapping projections produce unimodal R distributions and segregated projections produce bimodal R distributions, this kurtosis trend demonstrates increasing input separation with age. Linear regression depicts this developmental trend (Figure 10A). Representative R value histograms at each age reveal the transition from unimodal (overlapping) to bimodal (non-overlapping) from P0 to P12 (Figure 10B).
Figure 10
Discussion
The current study provides evidence that mapping of somatosensory inputs to discrete LCIC compartments occurs alongside that previously described for auditory inputs (Lamb-Echegaray et al., ). During this early critical period a process of segregation takes place such that inputs of different modalities recognize and target complementary zones of its emerging patch-matrix-like organization. At birth, prior to the earliest indications of any developing modularity, somatosensory and auditory inputs intermingle throughout the LCIC. From P4 to P8, coincident with the appearance of its discrete neurochemical compartments (Dillingham et al., ) and matching Eph-ephrin guidance patterns (Gabriele et al., ; Cramer and Gabriele, ; Wallace et al., 2016; Gay et al., ; Stinson et al., 2021), inputs of the two modalities begin to separate. By P12, inputs appear adult-like and fully segregated, with highly specific terminal fields that recognize LCIC compartmental boundaries and exhibit little evidence of any remaining projection overlap. Quantification of LCIC afferent patterns confirm this developmental trend that initially overlapping projection distributions soon yield non-overlapping, modality-specific input maps (i.e., somatosensory: modules, auditory: matrix).
Beyond the two inputs examined in the present study, other top-down and bottom-up inputs of somatosensory and auditory origin exist in the adult and exhibit similar preference for distinct LCIC compartments (Wiberg and Blomqvist, 1984; Saldaña and Merchán, 1992; Saldaña et al., 1996; Winer et al., 1998; Zhou and Shore, 2006; Bajo et al., ; Torii et al., 2013; Stebbings et al., 2014; Lesicko et al., ). With such multi-level convergence arising from multiple sensory systems, several mechanisms likely work in concert to instruct LCIC topographic mapping. If pioneer work in unimodal systems is any indication (McLaughlin and O’Leary, ; Miller et al., 2006; Luo and Flanagan, ; Cang et al., ; Triplett et al., 2009, 2011; Imai et al., ; Cang and Feldheim, ; James et al., ), it is probable that precise IC map formation relies upon both activity-dependent and activity-independent mechanisms. Complex mapping along multiple axes normally involves a combination of molecular guidance cues, patterned neural activity, and axon-axon competition. Unlike the neighboring CNIC that exhibits continuous map features (Fathke and Gabriele, ), molecular gradients (Gabriele et al., ; Wallace et al., 2013, 2016; Cramer and Gabriele, ), and cochlear-generated spontaneous events, the LCIC is discretely-organized (Lesicko et al., ; Lamb-Echegaray et al., ) with matching guidance patterns (Gay et al., ; Stinson et al., 2021) and intrinsically generated activity likely arising from multiple sources (Tritsch and Bergles, 2010; Wang and Bergles, 2015; Nakazawa et al., 2020; Babola et al., ).
In lieu of gradients, Eph-ephrin guidance molecules are transiently expressed in the LCIC in discontinuous patterns that correlate temporally and spatially with its emerging modularity and developing afferent arrays. Just prior to, and throughout the period that somatosensory and auditory projections separate within the LCIC, EphA4 and ephrin-B2 expression is patchy and overlaps GAD-positive modular domains (Gay et al., ). Ephrin-B3 expression, on the other hand, is complementary and localized to the surrounding CR-positive matrix (Stinson et al., 2021). After inputs have fully segregated (P12), expression of these signaling molecules is sharply downregulated. It is likely that discrete Eph-ephrin patterns provide necessary positional information for establishing rough somatosensory vs. auditory terminal compartments. Additional experimentation including loss-of-function studies in mice are needed to tease how specific Eph-ephrin interactions instruct modality-specific LCIC compartments. Once a course map is molecularly established, neural activity often provides subsequent fine-tuning or refinement, such that inputs with correlated activity patterns target similar neuronal subsets, more so than those exhibiting unrelated activity (Butts and Rokhsar, ; Debski and Cline, ; Cang et al., ; Triplett et al., 2012). Applying these insights from other systems/structures to the LCIC model, one might hypothesize that rhythmic spontaneous bursting in the auditory pathway is key for aligning multiple auditory maps within its matrix. Similarly, evidence for early somatotopic activity triggered by spontaneous muscle twitches (Khazipov et al., ; Luhmann et al., ) may prove equally important for the registry of multiple somatosensory inputs that map to its modular domains. Finally, the juxtaposition of multimodal projection domains (modular/matrix) carrying immature activity presumed to be uncorrelated (i.e., of either somatosensory or auditory origin) may serve to further demarcate LCIC compartmental boundaries. Further investigations that replicate the present tracing experiments and quantify the degree of projection overlap in mice lacking necessary Eph-ephrin cues, structured activity, or both, are needed to determine how these distinct afferent sets achieve proper alignment. Such data will be instrumental in developing future computational models that make LCIC mapping predictions based on various experimental conditions, similar to those already being successfully employed in other analogous systems (Reber et al., 2004; Yates et al., 2004; Tsigankov and Koulakov, 2006; Owens et al., 2015; Savier et al., 2017; Savier et al., 2020).
Besides the potential mechanisms mentioned above, molecules and signaling canonically associated with immune function have been implicated of late as important regulators of developmental pruning and synaptic plasticity (Faust et al., ). In particular, microglia as the resident macrophages of the CNS have received significant attention for their presence and significant roles in sculpting newly formed networks. Microglial selective engulfment of extraneous or underutilized contacts has been described in several developing systems and appears to be activity-dependent (Paolicelli et al., 2011; Hoshiko et al., ; Schafer et al., 2012; Stephan et al., 2012; Pagani et al., 2015; Hong et al., ; Thion and Garel, 2017). Among others, the classical complement cascade (C3–CR3) and fractalkine (CX3CL1–CX3CR1) signaling figure prominently in microglial cell (MGC) recruitment to areas of active circuit assembly and tagging of less active connections for subsequent phagocytic removal. Disruption of MGC function/signaling often results in sustained deficits or delayed maturation of synaptic function, as well as impairments in projection refinement. Compromised complement signaling affects the normal segregation of ipsi- and contralateral retinogeniculate inputs, with significant maintenance of initial projection overlap (Stevens et al., 2007; Schafer et al., 2012). In somatosensory cortex, fractalkine signaling controls entry of MGCs into barrel fields and CX3CR1-deficiency delays functional maturation of its connectivity (Hoshiko et al., ). Microglia appear to also perform similar functions in the developing auditory brainstem, evidenced by significant pruning deficits and abnormal auditory brainstem responses when MGCs are pharmacologically depleted (Milinkeviciute et al., , ). While pruning at the calyx of Held is unaffected in fractalkine receptor mutants, these mice express increased glycinergic synaptic markers in the auditory brainstem, suggesting an important role for CX3CL1–CX3CR1 crosstalk in sculpting immature inhibitory connections (Milinkeviciute et al., ). Recent work from our lab shows microglia indeed invade the nascent IC and that compromised fractalkine signaling delays MGC occupancy of its emerging modules (Brett et al., ). The present study quantifying the normal development and segregation of multimodal LCIC inputs serves as important baseline for ongoing experimentation aimed at determining whether MGCs and their established signaling pathways influence selective pruning at the level of the midbrain, and if compromised, whether multisensory distributions fail to segregate. If true, such lack of refinement and sustained overlap might overwhelm LCIC circuits. Instead of first segregating sensory information so that it can then be shared locally (Lesicko et al., ) before being passed on for further integration at the level of the deep superior colliculus (SC, among other multisensory centers), unrefined network configurations may preclude the LCIC from performing its function as a staging area for meaningful merging of the senses.
In conclusion, our experiments demonstrate that modality-specific LCIC projection patterns segregate over an early postnatal critical period. Beyond sharing a similar window for circuit assembly and mapping adjustments defined for other uni- and multimodal structures (Hoshiko et al., ; Schafer et al., 2012; Triplett et al., 2012), it appears the spatial registry of converging LCIC maps may too be governed by a combination of chemical labeling, neural activity, and competitive influences. Understanding the mechanisms that instruct its early plasticity and functional maturation is vital for gaining insights into certain neurodevelopmental conditions, such as autism spectrum disorders and schizophrenia. Both are characterized by gross errors in pruning (under- and over-pruning, respectively; Lehrman et al., ; Faust et al., ) and multimodal processing deficits that likely contribute to a host of behavioral consequences (Kwakye et al., ; Stevenson et al., 2014a, b; Brandwein et al., ; Tseng et al., 2015; Robertson and Baron-Cohen, 2017). Through further study of multisensory map alignments and the dynamic processes that guide their selective refinement, we may better understand the factors underlying these debilitating disorders and perhaps uncover new means for their successful treatment.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by James Madison University’s Animal Care and Use Committee (approval number, 20-1421).
Author contributions
JW, EK, JC, and MG all contributed to the presented experiments. JW, EK, and JC performed all tissue processing, associated imaging, and data sampling. MG performed data quantification, prepared all figures, and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by the National Institutes of Health (1R15DC018885-01), the National Science Foundation (DBI-0619207 and DBI-1725855), and the JMU Department of Biology Light Microscopy and Imaging Facility.
Acknowledgments
We would like to sincerely thank Dr. Thomas Gabriele for his signal analysis consultations, Dr. Kristopher Kubow for his microscopy expertise, and Stephanie Atkins and Sarah Keegan for their invaluable assistance with breeding and maintenance of our mouse colony.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncir.2022.882485/full#supplementary-material.
Supplementary Figure 1Two additional P12 cases (A–D,E–H) showing complementary, modality-specific LCIC afferent patterns. Somatosensory terminals (red) are concentrated within GAD-positive modules (green, dashed contours), while auditory inputs (cyan) distribute throughout the surrounding matrix. Projections have essentially segregated at this age (i.e. minimal overlap) with patterns in keeping with those described in the adult, suggesting the end of an early critical period for shaping discrete, multisensory LCIC input arrays. Scale bars = 100 μm.
References
1
AitkinL. M.DickhausH.SchultW.ZimmermannM. (1978). External nucleus of inferior colliculus: auditory and spinal somatosensory afferents and their interactions. J. Neurophysiol.41, 837–847. 10.1152/jn.1978.41.4.837
2
AitkinL. M.KenyonC. E.PhilpottP. (1981). The representation of the auditory and somatosensory systems in the external nucleus of the cat inferior colliculus. J. Comp. Neurol.196, 25–40. 10.1002/cne.901960104
3
BabolaT. A.LiS.WangZ.KersbergenC. J.ElgoyhenA. B.CoateT. M.et al. (2021). Purinergic signaling controls spontaneous activity in the auditory system throughout early development. J. Neurosci.41, 594–612. 10.1523/JNEUROSCI.2178-20.2020
4
BajoV. M.NodalF. R.BizleyJ. K.MooreD. R.KingA. J. (2007). The ferret auditory cortex: descending projections to the inferior colliculus. Cereb. Cortex17, 475–491. 10.1093/cercor/bhj164
5
BrandweinA. B.FoxeJ. J.ButlerJ. S.FreyH. P.BatesJ. C.ShulmanL. H.et al. (2015). Neurophysiological indices of atypical auditory processing and multisensory integration are associated with symptom severity in autism. J. Autism Dev. Disord.45, 230–244. 10.1007/s10803-014-2212-9
6
BrettC. A.CarrollJ. B.GabrieleM. L. (2022). Compromised fractalkine signaling delays microglial occupancy of emerging modules in the multisensory midbrain. Glia7, 697–711. 10.1002/glia.24134
7
ButtsD. A.RokhsarD. S. (2001). The information content of spontaneous retinal waves. J. Neurosci.21, 961–973. 10.1523/JNEUROSCI.21-03-00961.2001
8
CangJ.FeldheimD. A. (2013). Developmental mechanisms of topographic map formation and alignment. Annu. Rev. Neurosci.36, 51–77. 10.1146/annurev-neuro-062012-170341
9
CangJ.WangL.StrykerM. P.FeldheimD. A. (2008). Roles of ephrin-As and structured activity in the development of functional maps in the superior colliculus. J. Neurosci.28, 11015–11023. 10.1523/JNEUROSCI.2478-08.2008
10
ChernockM. L.LarueD. T.WinerJ. A. (2004). A periodic network of neurochemical modules in the inferior colliculus. Hear. Res.188, 12–20. 10.1016/S0378-5955(03)00340-X
11
CramerK. S.GabrieleM. L. (2014). Axon guidance in the auditory system: multiple functions of Eph receptors. Neuroscience277, 152–162. 10.1016/j.neuroscience.2014.06.068
12
DebskiE. A.ClineH. T. (2002). Activity-dependent mapping in the retinotectal projection. Curr. Opin. Neurobiol.12, 93–99. 10.1016/s0959-4388(02)00295-7
13
DillinghamC. H.GayS. M.BehroozR.GabrieleM. L. (2017). Modular-extramodular organization in developing multisensory shell regions of the mouse inferior colliculus. J. Comp. Neurol.525, 3742–3756. 10.1002/cne.24300
14
FathkeR. L.GabrieleM. L. (2009). Patterning of multiple layered projections to the auditory midbrain prior to experience. Hear. Res.249, 36–43. 10.1016/j.heares.2009.01.004
15
FaustT. E.GunnerG.SchaferD. P. (2021). Mechanisms governing activity-dependent synaptic pruning in the developing mammalian CNS. Nat. Rev. Neurosci.22, 657–673. 10.1038/s41583-021-00507-y
16
GabrieleM. L.BrubakerD. Q.ChamberlainK. A.KrossK. M.SimpsonN. S.KavianpourS. M. (2011). EphA4 and eprhin-B2 expression patterns during inferior colliculus projection shaping prior to experience. Dev. Neurobiol.71, 182–199. 10.1002/dneu.20842
17
GayS. M.BrettC. A.StinsonJ. P. C.GabrieleM. L. (2018). Alignment of EphA4 and ephrin-B2 expression patterns with developing modularity in the lateral cortex of the inferior colliculus. J. Comp. Neurol.526, 2706–2721. 10.1002/cne.24525
18
GrutersK. G.GrohJ. M. (2012). Sounds and beyond: multisensory and other non-auditory signals in the inferior colliculus. Front. Neural Circuits6:96. 10.3389/fncir.2012.00096
19
HongS.Dissing-OlesenL.StevensB. (2016). New insights on the role of microglia in synaptic pruning in health and disease. Curr. Opin. Neurobiol.36, 128–134. 10.1016/j.conb.2015.12.004
20
HoshikoM.ArnouxI.AvignoneE.YamamotoN.AudinatE. (2012). Deficiency of the microglial receptor CX3CR1 impairs postnatal functional development of thalamocortical synapses in the barrel cortex. J. Neurosci.32, 15106–15111. 10.1523/JNEUROSCI.1167-12.2012
21
ImaiT.SakanoH.VosshallL. B. (2010). Topographic mapping–the olfactory system. Cold Spring Harb. Perspect. Biol.2:a001776. 10.1101/cshperspect.a001776
22
JamesS. S.KrubitzerL. A.WilsonS. P. (2020). Modelling the emergence of whisker barrels. eLife9:e55588. 10.7554/eLife.55588
23
Jaubert-MiazzaL.GreenE.LoF.BuiK.MillsJ.GuidoW. (2005). Structural and functional composition of the developing retinogeniculate pathway in the mouse. Vis. Neurosci.22, 661–676. 10.1017/S0952523805225154
24
KhazipovR.SirotaA.LeinekugelX.HolmesG. L.Ben-AriY.BuzsákiG. (2004). Early motor activity drives spindle bursts in the developing somatosensory cortex. Nature432, 758–761. 10.1038/nature03132
25
KwakyeL. D.Foss-FeigJ. H.CascioC. J.StoneW. L.WallaceM. T. (2011). Altered auditory and multisensory temporal processing in autism spectrum disorders. Front. Integr. Neurosci.4:129. 10.3389/fnint.2010.00129
26
Lamb-EchegarayI. D.NoftzW. A.StinsonJ. P. C.GabrieleM. L. (2019). Shaping of discrete auditory inputs to extramodular zones of the lateral cortex of the inferior colliculus. Brain Struct. Funct.224, 3353–3371. 10.1007/s00429-019-01979-6
27
LehrmanE. K.WiltonD. K.LitvinaE. Y.WelshC. A.ChangS. T.FrouinA.et al. (2018). CD47 protects synapses from excess microglia-mediated pruning during development. Neuron100, 120–134.e6. 10.1016/j.neuron.2018.09.017
28
LesickoA. M. H.HristovaT. S.MaiglerK. C.LlanoD. A. (2016). Connectional modularity of top-down and bottom-up multimodal inputs to the lateral cortex of the mouse inferior colliculus. J. Neurosci.36, 11037–11050. 10.1523/JNEUROSCI.4134-15.2016
29
LesickoA. M. H.SonsS. K.LlanoD. A. (2020). Circuit mechanisms underlying the segregation and integration of parallel processing streams in the inferior colliculus. J. Neurosci.40, 6328–6344. 10.1523/JNEUROSCI.0646-20.2020
30
LuhmannH. J.SinningA.YangJ. W.Reyes-PuertaV.StüttgenM. C.KirischukS.et al. (2016). Spontaneous neuronal activity in developing neocortical networks: From single cells to large-scale interactions. Front. Neural Circuits10:40. 10.3389/fncir.2016.00040
31
LuoL.FlanaganJ. G. (2007). Development of continuous and discrete neural maps. Neuron56, 284–300. 10.1016/j.neuron.2007.10.014
32
McLaughlinT.O’LearyD. D. M. (2005). Molecular gradients and development of retinotopic maps. Annu. Rev. Neurosci.28, 327–355. 10.1146/annurev.neuro.28.061604.135714
33
MilinkeviciuteG.ChokrS. M.CramerK. S. (2021a). Auditory brainstem deficits from early treatment with a CSF1R inhibitor largely recover with microglial repopulation. eNeuro8:ENEURO.0318-20.2021. 10.1523/ENEURO.0318-20.2021
34
MilinkeviciuteG.ChokrS. M.CastroE. M.CramerK. S. (2021b). CX3CR1 mutation alters synaptic and astrocytic protein expression, topographic gradients and response latencies in the auditory brainstem. J. Comp. Neurol.529, 3076–3097. 10.1002/cne.25150
35
MilinkeviciuteG.HenningfieldC. M.MuniakM. A.ChokrS. M.GreenK. N.CramerK. S. (2019). Microglia regulate pruning of specialized synapses in the auditory brainstem. Front. Neural Circuits13:55. 10.3389/fncir.2019.00055
36
MillerK.KolkS. M.DonoghueM. J. (2006). EphA7-ephrin-A5 signaling in mouse somatosensory cortex: developmental restriction of molecular domains and postnatal maintenance of functional compartments. J. Comp. Neurol.496, 627–642. 10.1002/cne.20926
37
NakazawaS.YoshimuraY.TakagiM.MizunoH.IwasatoT. (2020). Developmental phase transitions in spatial organization of spontaneous activity in postnatal barrel cortex layer 4. J. Neurosci.40, 7637–7650. 10.1523/JNEUROSCI.1116-20.2020
38
OwensM. T.FeldheimD. A.StrykerM. P.TriplettJ. W. (2015). Stochastic interaction between neural activity and molecular cues in the formation of topographic maps. Neuron87, 1261–1273. 10.1016/j.neuron.2015.08.030
39
PaganiF.PaolicelliR. C.MuranaE.CorteseB.Di AngelantonioS.ZuroloE.et al. (2015). Defective microglial development in the hippocampus of Cx3cr1 deficient mice. Front. Cell. Neurosci.9:111. 10.3389/fncel.2015.00111
40
PaolicelliR. C.BolascoG.PaganiF.MaggiL.ScianniM.PanzanelliP.et al. (2011). Synaptic pruning by microglia is necessary for normal brain development. Science333, 1456–1458. 10.1126/science.1202529
41
PaxinosG.HallidayG. M.KoutcherovY.WangH.WatsonC. (2007). Atlas of the Developing Mouse Brain: E17.5, P0 and P6.London: Academic Press.
42
ReberM.BurrolaP.LemkeG. (2004). A relative signalling model for the formation of a topographic neural map. Nature431, 847–853. 10.1038/nature02957
43
RobertsonC. E.Baron-CohenS. (2017). Sensory perception in autism. Nat. Rev. Neurosci.18, 671–684. 10.1038/nrn.2017.112
44
SaldañaE.FelicianoM.MugnainiE. (1996). Distribution of descending projections from primary auditory neocortex to inferior colliculus mimics the topography of intracollicular projections. J. Comp. Neurol.371, 15–40. 10.1002/(SICI)1096-9861(19960715)371:1<15::AID-CNE2>3.0.CO;2-O
45
SaldañaE.MerchánM. A. (1992). Intrinsic and commissural connections of the rat inferior colliculus. J. Comp. Neurol.319, 417–437. 10.1002/cne.903190308
46
SavierE. L.DunbarJ.CheungK.ReberM. (2020). New insights on the modeling of the molecular mechanisms underlying neural maps alignment in the midbrain. eLife9:e59754. 10.7554/eLife.59754
47
SavierE.EglenS. J.BathélémyA.PerrautM.PfriegerF. W.LemkeG.et al. (2017). A molecular mechanism for the topographic alignment of convergent neural maps. eLife6:e20470. 10.7554/eLife.20470
48
SchaferD. P.LehrmanE. K.KautzmanA. G.KoyamaR.MardinlyA. R.YamasakiR.et al. (2012). Microglia sculpt postnatal neural circuits in an activity and complement-dependent manner. Neuron74, 691–705. 10.1016/j.neuron.2012.03.026
49
SlaterB. J.WillisA. M.LlanoD. A. (2013). Evidence for layer-specific differences in auditory corticocollicular neurons. Neuroscience229, 144–154. 10.1016/j.neuroscience.2012.10.053
50
StebbingsK. A.LesickoA. M.LlanoD. A. (2014). The auditory corticocollicular system: molecular and circuit-level considerations. Hear. Res.314, 51–59. 10.1016/j.heares.2014.05.004
51
StephanA. H.BarresB. A.StevensB. (2012). The complement system: an unexpected role in synaptic pruning during development and disease. Annu. Rev. Neurosci.35, 369–389. 10.1146/annurev-neuro-061010-113810
52
StevensB.AllenN. J.VazquezL. E.HowellG. R.ChristophersonK. S.NouriN.et al. (2007). The classical complement cascade mediates CNS synapse elimination. Cell131, 1164–1178. 10.1016/j.cell.2007.10.036
53
StevensonR. A.SiemannJ. K.SchneiderB. C.EberlyH. E.WoynaroskiT. G.CamarataS. M.et al. (2014a). Multisensory temporal integration in autism spectrum disorders. J. Neurosci.34, 691–697. 10.1523/JNEUROSCI.3615-13.2014
54
StevensonR. A.SiemannJ. K.WoynaroskiT. G.SchneiderB. C.EberlyH. E.CamarataS. M.et al. (2014b). Evidence for diminished multisensory integration in autism spectrum disorders. J. Autism Dev. Disord.44, 3161–3167. 10.1007/s10803-014-2179-6
55
StinsonJ. P. C.BrettC. A.CarrollJ. B.GabrieleM. L. (2021). Registry of compartmental ephrin-B3 guidance patterns with respect to emerging multimodal midbrain maps. Front. Neuroanat.15:649478. 10.3389/fnana.2021.649478
56
TamamakiN.YanagawaY.TomiokaR.MiyazakiJ.ObataK.KanekoT. (2003). Green fluorescent protein expression and colocalization with calretinin, parvalbumin and somatostatin in the GAD67-GFP knock-in mouse. J. Comp. Neurol.467, 60–79. 10.1002/cne.10905
57
ThionM. S.GarelS. (2017). On place and time: microglia in embryonic and perinatal brain development. Curr. Opin. Neurobiol.47, 121–130. 10.1016/j.conb.2017.10.004
58
TorborgC. L.FellerM. B. (2004). Unbiased analysis of bulk axonal segregation patterns. J. Neurosci. Methods135, 17–26. 10.1016/j.jneumeth.2003.11.019
59
TorborgC. L.HansenK. A.FellerM. B. (2005). High frequency, synchronized bursting drives eye-specific segregation of retinogeniculate projections. Nat. Neurosci.8, 72–78. 10.1038/nn1376
60
ToriiM.HackettT. A.RakicP.LevittP.PolleyD. B. (2013). EphA signaling impacts development of topographic connectivity in auditory corticofugal systems. Cereb. Cortex23, 775–785. 10.1093/cercor/bhs066
61
TriplettJ. W.OwensM. T.YamadaJ.LemkeG.CangJ.StrykerM. P.et al. (2009). Retinal input instructs alignment of visual topographic maps. Cell139, 175–185. 10.1016/j.cell.2009.08.028
62
TriplettJ. W.PfeiffenbergerC.YamadaJ.StaffordB. K.SweeneyN. T.LitkeA. M.et al. (2011). Competition is a driving force in topographic mapping. Proc. Natl. Acad. Sci. U S A108, 19060–19065. 10.1073/pnas.1102834108
63
TriplettJ. W.PhanA.YamadaJ.FeldheimD. A. (2012). Alignment of multimodal sensory input in the superior colliculus through a gradient matching mechanism. J. Neurosci.32, 5264–5271. 10.1523/JNEUROSCI.0240-12.2012
64
TritschN. X.BerglesD. E. (2010). Developmental regulation of spontaneous activity in the mammalian cochlea. J. Neurosci.30, 1539–1550. 10.1523/JNEUROSCI.3875-09.2010
65
TsengH. H.BossongM. G.ModinosG.ChenK. M.McGuireP.AllenP. (2015). A systematic review of multisensory cognitive-affective integration in schizophrenia. Neurosci. Biobehav. Rev.55, 444–452. 10.1016/j.neubiorev.2015.04.019
66
TsigankovD. N.KoulakovA. A. (2006). A unifying model for activity-dependent and activity-independent mechanisms predicts complete structure of topographic maps in ephrin-A deficient mice. J. Comput. Neurosci.21, 101–114. 10.1007/s10827-006-9575-7
67
WallaceM. M.HarrisJ. A.BrubakerD. Q.KlotzC. A.GabrieleM. L. (2016). Graded and discontinuous EphA-ephrinB expression patterns in the developing auditory brainstem. Hear. Res.335, 64–75. 10.1016/j.heares.2016.02.013
68
WallaceM. M.KavianpourS. M.GabrieleM. L. (2013). Ephrin-B2 reverse signaling is required for topography but not pattern formation of lateral superior olivary inputs to the inferior colliculus. J. Comp. Neurol.521, 1585–1597. 10.1002/cne.23243
69
WangH. C.BerglesD. E. (2015). Spontaneous activity in the developing auditory system. Cell Tissue Res.361, 65–75. 10.1007/s00441-014-2007-5
70
WibergM.BlomqvistA. (1984). The projection to the mesencephalon from the dorsal column nuclei. An anatomical study in the cat. Brain Res.311, 225–244. 10.1016/0006-8993(84)90086-6
71
WinerJ. A.LarueD. T.DiehlJ. J.HeftiB. J. (1998). Auditory cortical projections to the cat inferior colliculus. J. Comp. Neurol.400, 147–174.
72
YatesP. A.HolubA. D.McLaughlinT.SejnowskiT. J.O’LearyD. D. M. (2004). Computational modeling of retinotopic map development to define contributions of EphA-ephrinA gradients, axon-axon interactions and patterned activity. J. Neurobiol.59, 95–113. 10.1002/neu.10341
73
ZhouJ.ShoreS. (2006). Convergence of spinal trigeminal and cochlear nucleus projections in the inferior colliculus of the guinea pig. J. Comp. Neurol.495, 100–112. 10.1002/cne.20863
Summary
Keywords
development, inferior colliculus, multisensory, GAD, anterograde, modules, matrix
Citation
Weakley JM, Kavusak EK, Carroll JB and Gabriele ML (2022) Segregation of Multimodal Inputs Into Discrete Midbrain Compartments During an Early Critical Period. Front. Neural Circuit 16:882485. doi: 10.3389/fncir.2022.882485
Received
23 February 2022
Accepted
18 March 2022
Published
07 April 2022
Volume
16 - 2022
Edited by
Jason W. Triplett, Children’s National Hospital, United States
Reviewed by
Guillermina Lopez-Bendito, Miguel Hernández University of Elche, Spain; Daniel Llano, University of Illinois at Urbana-Champaign, United States
Updates
Copyright
© 2022 Weakley, Kavusak, Carroll and Gabriele.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mark L. Gabriele gabrieml@jmu.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.