Skip to main content

Characterization of early markers of disease in the mouse model of mucopolysaccharidosis IIIB

Abstract

Background

Mucopolysaccharidosis (MPS) IIIB, also known as Sanfilippo Syndrome B, is a devastating childhood disease. Unfortunately, there are currently no available treatments for MPS IIIB patients. Yet, animal models of lysosomal storage diseases have been valuable tools in identifying promising avenues of treatment. Enzyme replacement therapy, gene therapy, and bone marrow transplant have all shown efficacy in the MPS IIIB model systems. A ubiquitous finding across rodent models of lysosomal storage diseases is that the best treatment outcomes resulted from intervention prior to symptom onset. Therefore, the aim of the current study was to identify early markers of disease in the MPS IIIB mouse model as well as examine clinically-relevant behavioral domains not yet explored in this model.

Methods

Using the MPS IIIB mouse model, we explored early developmental trajectories of communication and gait, and later social behavior, fear-related startle and conditioning, and visual capabilities. In addition, we examined brain structure and function via magnetic resonance imaging and diffusion tensor imaging.

Results

We observed reduced maternal isolation-induced ultrasonic vocalizations in MPS IIIB mice relative to controls, as well as disruption in a number of the spectrotemporal features. MPS IIIB also exhibited disrupted thermoregulation during the first two postnatal weeks without any differences in body weight. The developmental trajectories of gait were largely normal. In early adulthood, we observed intact visual acuity and sociability yet a more submissive phenotype, increased aggressive behavior, and decreased social sniffing relative to controls. MPS IIIB mice showed greater inhibition of startle in response to a pretone with a decrease in overall startle response and reduced cued fear memory. MPS IIIB also weighed significantly more than controls throughout adulthood and showed larger whole brain volumes and normalized regional volumes with intact tissue integrity as measured with magnetic resonance and diffusion tensor imaging, respectively.

Conclusions

Together, these results indicate disease markers are present as early as the first two weeks postnatal in this model. Further, this model recapitulates social, sensory and fear-related clinical features. Our study using a mouse model of MPS IIIB provides essential baseline information that will be useful in future evaluations of potential treatments.

Background

Mucopolysaccharidosis (MPS) IIIB, also known as Sanfilippo Syndrome B, is a devastating childhood disease caused by a known lysosomal enzyme deficiency [1, 2]. Specifically, a complete loss of the lysosomal enzyme alpha-N-acetylglucosaminidase (NAGLU) results in the accumulation of heparin sulfate and serious neurological issues [3]. Disease onset occurs around age five with developmental and speech delays, and the rise of somatic signs like hearing loss [1,2,3]. Throughout the disease progression, patients experience cognitive deterioration, sleep disturbances, and behavioral difficulties including hyperactivity and aggression, which are among the most distressing symptoms, as reported by parents [4, 5]. In addition, these children show autistic-like social and emotional abnormalities. Imaging studies revealed cerebellar atrophy and white matter thinning accompany the behavioral symptomatology yet do not track with disease severity [6, 7]. As the disease progresses, these children develop a lack of fear, motor decline, and eventual onset of severe dementia and early death.

Regrettably, there are currently no available treatments for MPS IIIB patients. Yet, animal models of lysosomal storage diseases, including for MPS IIIB, have been valuable tools in identifying promising avenues of treatment, including gene therapy, bone marrow transplant, and enzyme replacement therapy [8,9,10,11,12]. While not FDA approved yet, both gene therapy and enzyme replacement therapy have been ongoing clinical trials [13, 14]. The MPS IIIB murine model recapitulates many aspects of the human disease, including altered circadian rhythms, juvenile-onset hearing loss, later-age onset of balance issues and phenotypes suggestive of reduced fear [15,16,17,18,19]. However, many important markers of MPS IIIB disease remain unexplored in the mouse model, such as aggression, visual decline, sensorimotor gating, and fear conditioning. Thus, the emerging therapeutics have not yet been evaluated against these phenotypes.

In addition, a ubiquitous finding across rodent models of lysosomal storage diseases is that the best treatment outcomes resulted from intervention prior to symptom onset [20,21,22,23,24]. Therefore, identifying the earliest age at onset of behavioral abnormalities is necessary to determine the ideal age for treatment. Specifically, the MPS IIIB murine model has not been evaluated for neurobehavioral phenotypes prior to four weeks of age, missing possible early markers of developmental delay [15,16,17,18,19]. In addition, this model has not been well-characterized for behaviors related to motor function such as gait, and social and emotional abnormalities, which would be key markers of disease for evaluation of therapy effectiveness.

In the current study we assessed phenotypes in early postnatal development and expanded our assays to include clinically-relevant behavioral domains not yet explored in this model. Histopathological and biochemical analyses of disease require tissue collection and euthanasia of the animal, and therefore only a single time point can be examined. Neurobehavioral analyses are non-invasive and allow for the coupling of multiple areas of assessment in a single animal. Thus, in addition to early communicative behaviors and developmental motor trajectories, we assessed social and agonistic behaviors, visual acuity, sensorimotor gating, and fear conditioning in MPS IIIB mice. We coupled our behavioral assessment with neuroimaging, including magnetic resonance imaging to examine whole brain volumes, cerebellar structure and diffusion tensor imaging to interrogate structural and white matter integrity, and correlated these data with behavioral phenotypes. The findings from our study provide baseline phenotypic information for the MPS IIIB model that will be useful in future evaluations of potential treatments.

Methods

Animals

All experimental protocols were approved by and performed in accordance with the relevant guidelines and regulations of the Institutional Animal Care and Use Committee of Washington University in St. Louis and were in compliance with US National Research Council's Guide for the Care and Use of Laboratory Animals, the US Public Health Service's Policy on Humane Care and Use of Laboratory Animals, and Guide for the Care and Use of Laboratory Animals.

All mice used in this study were maintained and bred in the vivarium at Washington University in St. Louis from an existing colony at Washington University. For all experiments, adequate measures were taken to minimize any pain or discomfort. The colony room lighting was on a 12:12 h light/dark cycle (lights on 6am to 6 pm); room temperature (~ 20–22 °C) and relative humidity (50%) controlled automatically. Standard lab diet and water were freely available. Upon weaning at postnatal day (P)21, mice for behavioral testing were group-housed according to sex and genotype.

The MPS IIIB mouse model, maintained on a C57BL/6J background, containing a neomycin cassette disruption of exon 6 in the Naglu gene was used for this study. Generation of this model was described previously [19]. Naglu+/- (Control) mice were crossed with Naglu−/− (MPS IIIB) mice to generate male and female Control and MPS IIIB littermates for the study. MPS IIIB is an autosomal recessive disorder, therefore heterozygous mutants are reportedly normal [17] and serve as the appropriate littermate control. We examined the phenotypes of three independent cohorts. The first cohort included 10 litters (4 from heterozygous control sire x MPS IIIB dam pairings and 6 from heterozygous control sire x dam pairs) comprising 27 Control (13 female and 14 male) and 17 MPS IIIB (4 female and 13 male) mice. The second cohort included 8 litters (3 from MPS IIIB sire x heterozygous control dam pairings and 5 from heterozygous control sire x dam pairs) comprising 30 Control (10 female and 20 male) and 16 MPS IIIB (8 female and 8 male) mice. The third cohort included 7 litters (3 from MPS IIIB sire x heterozygous control dam pairings, 2 from heterozygous control sire x MPS IIIB dam pairings, and 2 from heterozygous control sire x dam pairs) comprising 20 Control (10 female and 10 male) and 19 MPS IIIB (9 female and 10 male) mice.

Behavioral tasks

All behavioral data were collected in the Animal Behavior Subunit of the Intellectual and Developmental Disabilities Research Center at Washington University in St. Louis. Cohorts 1 and 2 were assessed for developmental trajectories from early postnatal ages through young adulthood (Table 1; Figs. 1A, 2A). Cohort 3 was assessed for additional, clinically-relevant behavioral phenotypes beginning in adolescence through early adulthood (Table 1; Fig. 3A). For all tasks, the mice were acclimated to the testing room at least 30 min prior to the start of testing. All assays were conducted by female experimenters blinded to experimental group designations during testing, and all testing occurred during the light phase. Order of tests (Table 1) was chosen to minimize effects of stress.

Table 1 Order of and age at behavioral testing per cohort with related patient phenotype and specific assessment in mice
Fig. 1
figure 1

MPS IIIB postnatal pups exhibit reduced body temperature and USVs with altered spectrotemporal features. A Timeline of developmental assessments. B Representative spectrogram segment from a Control mouse and a MPS IIIB mouse. C-D Number of ultrasonic vocalizations (USVs) during maternal isolation-induced USV assay across ages (C) and collapsed for age (D). E Body weight across postnatal assessment ages. F Latency to exhibit righting reflex at P14. G Body temperatures recorded immediately prior to USV recordings at P6-10. H Covariate-adjusted means of number of USVs accounting for differences in body temperature. I Number of phrases of USV calls. J Average USV duration. K Average pause duration between USVs within a phrase. L Sound pressure level or volume of USV calls. M Mean USV frequency. N Frequency range of USV calls. O Percent of calls that included a frequency jump greater than 10 kHz. C, E-J, M, and N data are presented as means ± SEM. D, K, L, and O data are presented as boxplots with thick horizontal lines respective group medians, boxes 25th – 75th percentiles, and whiskers 1.5 × IQR. Individual data points are presented as open circles

Fig. 2
figure 2

MPS IIIB mice exhibit subtle yet persistent alterations in gait. A Schematic of gait data collection timeline. B Number of data collection attempts needed across developmental timepoints. C-D Weight (C) and body length (D) across developmental gait collection ages. E Duration of braking for forelimbs during developmental gait assessment. F Hindpaw angle across developmental gait assessment timepoints. G-H Adult body length and weight for both females and males. I Number of data collection attempts needed for adult mice at 30 cm/s and 40 cm/s belt speeds. J Duration of braking for forelimbs at 30 cm/s belt speed for adult mice. K Percent stride in hindlimb shared stance at 40 cm/s belt speed for adult mice. L Percent of stride in swing (left panel) and stance (right panel) for hindlimbs at 40 cm/s belt speed for adult mice. Observed p values reported in figures. B-D, G, and H Grouped data presented as means ± SEM with individual data points presented as open circles. E, F, and I-L Covariate-adjusted means ± SEM accounting for changes in body length presented

Fig. 3
figure 3

MPS IIIB mice weigh significantly more than controls with unaltered vision in adolescence. A Schematic of behavioral testing timeline for cohort 3. B-D Body weights as measured from adolescence through four months of age for all mice (left panel), females only (middle panel), and male mice only (right panel). E–F VOS measurements of mean visual acuity levels (E) and mean contrast levels (F). B-F Grouped data presented as means ± SEM with individual data points presented as open circles

Maternal isolation-induced ultrasonic vocalizations

During the first two weeks postnatal, we assessed the control and MPS IIIB littermates in cohort 1 for signs of gross developmental delay, communicative delay or motor delay (Fig. 1A; Table 1), following our previously published methods [25]. To evaluate gross development, the mice were weighed daily from P6—P14, and evaluated for physical milestones of development including pinna detachment by P6 and eye opening by P14. While human language cannot be explored in mice, vocal communication behavior is conserved across taxa [26]. Ultrasonic vocalizations (USVs) produced by the mouse pup is one of the earliest forms of social communication we can examine in the mouse, and these isolation calls elicit maternal care from the dam [27]. This behavior also has a developmental trajectory, beginning just after birth, peaking during the first week postnatal and disappearing around P14, making it useful for examining delay in early social circuits. Tissue was collected from all pups on P5 for genotyping. USVs were recorded on P6, P8, and P10 following our previously published methods [25, 28]. Briefly, the dam was removed from the nest and the litter placed in a warming cabinet to maintain body temperature while away from the dam. The surface temperature of each pup was recorded (HDE Infrared Thermometer) prior to placement in an empty polycarbonate cage (28.5 cm × 17.5 cm x 12 cm) in a sound-attenuating chamber. USVs were recorded for three minutes using an Avisoft UltraSoundGate CM16 microphone, Avisoft UltraSoundGate 116H amplifier, and Avisoft Recorder software (gain = 8 dB, 16 bits, sampling rate = 250 kHz). The pup was then weighed and returned to the nest. The dam was returned to the nest following recording of the last pup in the litter. Frequency sonograms were prepared from recordings using MATLAB (frequency range = 25 kHz to 120 kHz, FFT size = 512, overlap = 50%, time resolution = 1.024 ms, frequency resolution = 488.2 Hz). Individual syllables and other spectral features were identified and counted from the sonograms as previously described [29, 30]. To assess achievement of motor milestones, the surface righting reflex was assessed for each pup on P14 [25]. Each pup was placed on its back in an empty cage lined with a plastic bench pad and the time to return to a prone position was recorded up to 60 s. Three trials were averaged for analysis.

DigiGait

Developmental gait trajectories were assessed using the DigiGait Gait Analysis System (Mouse Specifics, Inc) following our previously published methods [31, 32]. Briefly, for all mice in cohort 2, gait metrics were assessed on P21, 24, 27, 30 at 20 cm/s and again at one time point after P60 (adulthood) at both 30 cm/s and 40 cm/s (Fig. 2A; Table 1). Prior to P21 testing, mice were habituated to the DigiGait on P20 by placing the mouse on the dark, motionless apparatus belt for two minutes, followed by two minutes in the light. Then, the belt was turned on starting at 5 cm/s and increased slowly to 20 cm/s until the mouse had about 30 s of uninterrupted running. The litters were weaned from the dam and sire prior to testing on P21. On test days, a ‘failure to run’ criterion of five attempts was implemented. If the animal failed to run on the belt after five attempts, it was returned to its homecage. Video processing and body length assessments were conducted as previously described [31, 32].

Social approach

Three chamber social approach task was used to assess sociability and preference for social novelty at P40 (Fig. 5A) using procedures adapted from our previous work [25]. Briefly, stimulus partner mice were sex-, age- and strain-matched. The full apparatus measured 60 cm x 39 cm x 22 cm with the center chamber (9.75 cm x 39 cm x 22 cm) and outer chambers (24.2 cm x 39 cm x 22 cm) separated by opaque walls. Stimulus cups were placed 2 cm from the external wall. Testing was completed in white lighting at 115 lx. Behavior of the animals was tracked using Any-maze tracking software (Stoelting, Co). Between animals, the stimulus cups were cleaned with 70% ethanol and the acrylic apparatus was cleaned with 0.02% unscented Chlorhexidine.

Virtual Optomotor System

Visual acuity and contrast testing was completed using the Virtual Optomotor System (VOS; Cerebral Mechanics Inc.) adapted from previous methods [33, 34]. Testing assessed the sharpness of vision and the ability to distinguish between the foreground and background at P47. The VOS system apparatus consisted of a chamber (52.5 cm × 52.2 cm x 95 cm) with a centered platform (5.5 cm diameter) raised 19 cm surrounded by 4 monitors. To measure acuity, vertical gratings rotated across the screens and varied sizes across testing until a threshold was met at which mice could no longer distinguish the moving lines. A response was measured by the reflexive head movement (optomotor reflex) in concert with the rotation, scored visually by experimenter. Visual contrast is the ability to distinguish the foreground from the background. Contrast assessment used a constant grating size (0.128 cycles/degree), with gratings varying in shade until the mouse could no longer distinguish movement and elicit the optomotor reflex. The direction of movement in the VOS system allowed for each eye to be tested independently. Threshold values of clockwise (CW) and counterclockwise (CCW) were analyzed and averaged to create variables for mean acuity and mean contrast. Testing was completed between 6am and 12 pm when mice are less active and thus the response more readily distinguished. Mice were habituated to the apparatus 1 day prior to testing for 10 min, with the screens on but no stimulus present. Acuity testing occurred on day 2 and contrast testing on day 3. Females and males were tested on alternating days. Fifty trials were completed per eye using a staircase method to determine threshold. Animals with extreme outlier values were retested at the end of the day for accurate assessment.

Acoustic startle response/prepulse inhibition

Sensorimotor gating and startle reflex was evaluated at P53 using the Acoustic Startle/Prepulse Inhibition (PPI) assay following our previously published methods [25]. Startle response to a 120 dB auditory stimulus pulse (40 ms broadband burst) and PPI (response to a prepulse plus the startle pulse) were measured concurrently in the mice. Beginning at stimulus onset, 1 ms force readings were averaged to obtain an animal's startle amplitude. A total of 20 startle trials were presented over a 20 min test period during which the first 5 min served as an acclimation period when no stimuli above the 65 dB white noise background were presented. The session began and ended by presenting 5 consecutive startle (120 db pulse alone) trials unaccompanied by other trial types. The middle 10 startle trials were interspersed with PPI trials (consisting of an additional 30 presentations of 120 dB startle stimuli preceded by pre-pulse stimuli of either 4, 8, or 16 dB above background (10 trials for each PPI trial type). A percent PPI score for each trial was calculated using the following equation: %PPI = 100*(ASRstartle pulse alone—ASRprepulse + startle pulse)/ASRstartle pulse alone. Startle responses to increasing decibel levels was used as an auditory function screen at the end of the PPI trial session, and a paired t-test for response to 10 120 dB startle trials vs non-startle trials was used to determine deafness.

Social dominance tube test

Laboratory mice acquire dominance behaviors to create social hierarchy ranks within their cage environments, or social groups, between six to eight weeks of age, which can be leveraged to examine normal social dominance behavior. Thus, we used the Tube Test to assess social dominance behavior at P62 between sex-matched control and MPS IIIB animal pairs following our previously published methods (Fig. 5J) [25].

Conditioned fear

Associative anxiety-related memory was assessed at P67 in the conditioned fear task (Fig. 4F). The procedure details were as we previously published [25]. Briefly, the animals were exposed to a 1.0 mA current paired with a 80 dB white noise tone three times on day 1. Contextual fear conditioning was assessed on day 2 by placing the animal back into the chambers from day 1 without shock or tone presentation. On day 3, the animals were placed into a contextually altered chamber and assessed for cued fear conditioning in the presence of the 80 dB tone from day 1. On all three days, freezing was quantified by FreezeFrame software (Actimetrics). Sensitivity to foot shocks was evaluated one week following the start of conditioned fear testing to confirm differences in freezing behavior were not confounded by differences in reactivity to the shock current, based on our previously published methods [35]. Mice were placed in the fear conditioning chambers and exposed to 3 s footshock every 20–30 s starting at 0.05 mA and increasing in intensity by 0.05 mA until a flinch, vocalization and escape behavior were elicited. The current at which each behavior occurred was recorded.

Fig. 4
figure 4

MPS IIIB mice exhibit altered responses to negative stimuli suggestive of reduced fear by early adulthood. A Weight measurements recorded during Acoustic Startle/PPI task. B Startle response to 120 dB acoustic stimulus across 4 trials. C Startle response to consecutive trials of acoustic stimuli presentations of increasing dB level. D Inhibition of startle response to 120 dB stimuli following prepulse tone 4, 8, or 16 dB above background. E Percent of mice per group that exhibited a significant difference in response between startle and non-startle trials. F Schematic of fear conditioning paradigm. G Percent time freezing during baseline and tone/shock conditioning on day one. H Percent time freezing in response to fear context on day two. I Percent time freezing on day three during baseline in a new context and during tone presentation. A, and G-I Grouped data are presented as means ± SEM with individual data points presented as open circles. B-D Grouped data are presented as covariate adjusted means ± SEM accounting for differences in body weight. E Data are presented as counts

Resident intruder assay

To evaluate aggressive and agonistic behaviors, we tested P96 male mice in the resident intruder paradigm based on our previous study [36]. Male mice were single housed and cages remained unchanged for ten days prior to testing to establish a territory. On day 6, a sexually naive strain- and age-matched female was added to each male cage and removed 24 h prior to testing to potentiate territorial behaviors. Intruders used were strain- and age-matched novel male mice bred in the same colony space. On the test day, under white light, the intruder mouse was added directly into the resident home cage in a sound-attenuating box, and allowed to freely interact for 10 min. Interactions were filmed with a digital camera (Sony HandyCam). Upon completion, the intruder was removed and placed in a clean holding cage until all cage-mates complete testing. This was repeated over two additional days with a novel intruder mouse. We monitored interactions for any attacks drawing blood, which would necessitate prematurely ending the assay. However, we did not observe any attacks at this level.

To process the videos, we used pose estimation software (DeepLabCut, version 2.2.1) to track body parts in space and time [37, 38], building upon a previously built model [39]. Specifically, we labeled 50 frames taken from 60 videos plus an additional 1700 frames of closely interacting animals, and 95% were used for training. We used a dlcrnet_ms5-based neural network with default parameters for 200,000 training iterations. We validated with two shuffles, and found the test error was 7.95 pixels, with the train error at 5.58 pixels (image size was 1920 by 1080). We then used a p-cutoff of 0.6 to condition the X,Y coordinates for future analysis. This network was then used to analyze all videos in the experiment.

The pose estimation was then fed through the Simple Behavior Analysis (SimBA) random forest classifier [40] to extract attack behavior. Training videos were annotated for attack behavior using the SimBA event-logger using 180 annotated behavior files, downloaded from https://osf.io/tmu6y/ in addition to four in-house annotated videos. All training files were annotated according to definitions found in the simBA preprint [39, 40]. Random forest classifiers were trained using default hyperparameters, and classifier performances were evaluated. We set the discrimination threshold for scramble attacks to 0.65. The minimum bout length was 200 ms. Once classifiers were validated, predictive classifiers were applied to all videos in the dataset. Statistics regarding outcome variables of total time, the number of frames, total number of instances, mean instance interval, and mean interval between behavior instances were generated.

Magnetic resonance imaging

All in vivo magnetic resonance imaging (MRI) experiments were performed on a 9.4 T Bruker BioSpec 94/20 MR system (Bruker, Ettlingen, Germany) equipped with an 11.4-cm inner-diameter gradient insert. Data were collected with a 4-channel-array receive-only mouse-brain CryoProbe in combination with a quadrature transmit-only coil. Mice were anesthetized using 1.0–1.5% isoflurane/O2 and were secured in the prone position. Respiratory rate and rectal temperature were monitored using a small-animal physiological monitoring system (SA Instruments; Stony Brook, NY) and mouse body temperature was maintained at 37 °C using a circulating warm water pad.

Diffusion tensor imaging

Diffusion tensor imaging (DTI) data were collected using a spin-echo diffusion-weighted sequence with 31 b-values, max b-value = 1000 s/mm2. Other MR acquisition parameters include repetition time (TR) = 1.5 s; echo time (TE) = 30 ms; diffusion-encoding times, Δ = 18 ms, δ = 6 ms; field of view = 16 mm × 16 mm; matrix size = 128 × 128; 31 axial slices of thickness 0.5 mm. The total acquisition time was approximately one hour.

DTI data processing

DTI datasets were analyzed via non-linear least squares modeling of the diffusion-weighted signal using MATLAB (The MathWorks, Natick MA). Firstly, the data were zero-padded to matrix size 256 × 256, resulting in an in-plane resolution of 62.5 mm × 62.5 mm. The Gibbs-ringing artifact was reduced using a previously published method [41] followed by Gaussian filtering (Sigma 0.75) to denoise and smooth images. Tensor eigenvalues and parametric maps of apparent diffusion coefficient (ADC) and fractional anisotropy (FA) were calculated according to standard equations [42, 43].

All parametric maps were converted to NIfTI format for subsequent analysis. Whole brain (15.5 mm in anterior–posterior direction) and brain-substructure regions of interest (ROI) were segmented manually using ITK-SNAP [44]. Total brain volumes were calculated by segmenting T2-weighted (b = 0) images. The cerebellum and corpus callosum (CC) were segmented based on T2-weighted images and maps of ADC and FA. Volumes and DTI parameters were extracted for subsequent statistical analysis. Cortical thickness measurements were conducted using ITK-SNAP based on previously published methods [45]. A ventral to dorsal measurement was taken from the peak of the cingulum bundle to the dorsal surface of the cortex. This was repeated in both hemispheres across two sections for each brain, which were averaged for analysis.

Statistical analyses

Statistical analyses and data visualization were conducted using IBM SPSS Statistics (v.28). Prior to analyses, data were screened for missing values and fit of distributions with assumptions underlying univariate analysis. This included the Shapiro–Wilk test on z-score-transformed data and qqplot investigations for normality, Levene’s test for homogeneity of variance, and boxplot and z-score (± 3.29) investigation for identification of influential outliers. Means and standard errors were computed for each measure. Linear modeling including mixed designs and analysis of variance (ANOVA), including repeated measures and covariate designs, were used to analyze data where appropriate, and simple main effects were used to dissect significant interactions. Sex was included as a biological variable in all analyses across all experiments. Where appropriate, the Greenhouse–Geisser or Huynh–Feldt adjustment was used to protect against violations of sphericity for repeated measures designs. Multiple pairwise comparisons were subjected to Bonferroni correction or the Benjamini–Hochberg method for false discovery rate, where appropriate. For data that did not fit univariate assumptions, non-parametric tests were used or data transformations were applied. Sex*genotype effects are reported where significant, otherwise data are discussed and visualized collapsed for sex, with details of all statistical tests and results reported in Table 2. Non-significant predictors were removed from final models to achieve the most parsimonious models. The critical alpha value for all analyses was p < 0.05 unless otherwise stated. Pearson’s correlation was used to assess relationships between continuous outcomes and one-sample t-tests were used to assess differences from chance levels. Figure illustrations were generated using BioRender. The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.

Table 2 Statistical analysis results

Results

MPS IIIB mice exhibited disrupted early communicative behaviors with subtle yet persistent changes to the developmental gait trajectories

Speech delay is a first symptom of MPS IIIB making it an important phenotype to mark age at disease onset. While human language cannot be explored in mice, vocal communication behavior is conserved across taxa [26]. Mouse pups produce ultrasonic isolation calls as a way to attract the dam for maternal care [27], thus it is one of the earliest forms of social communication we can examine in mice. In addition, developing mice exhibit a trajectory of gait development that establishes ratios of stride components similar to that observed in humans [31], which can be leveraged to identify sensitive markers of motor impairment [32]. Here, we characterize the trajectories of call production (Fig. 1A) and gait in the MPS IIIB mice to identify early markers of disease.

MPS IIIB mice produced fewer calls overall compared to controls during the first two weeks of life (Fig. 1B-D; Table 2). Additional files show the full spectrogram for a representative MPS IIIB and control mouse (see Additional Files 1 and 2). This was not due to gross developmental delay as MPS IIIB mice body weights and latency to exhibit righting reflex at P14 were comparable to controls (Fig. 1E,F), and no differences were observed in number of mice that exhibited pinnae detachment at P6 (Controls, n = 27/27; MPS IIIB n = 17/17) or eye opening at P14 (Controls, 12.2% closed; MPS IIIB 12.5% closed). MPS IIIB mice also displayed significantly lower body temperatures relative to controls (Fig. 1G). Lower temperature is typically reported to increase the call rate to elicit maternal care and warmth [46]. However, we observe a reduction in call rate in the MPS IIIB mice, not an increase. To further confirm temperature was not modulating call rate (despite the misaligned direction of effect), we conducted a covariate model regressing out temperature and the reduced call rate was maintained (Fig. 1H), indicating the altered body temperature is not driving the call phenotype.

We next examined the spectrotemporal features of the calls to determine if physical anomalies, such as altered physiology of the lungs or larynx, may be heavily driving the reduction in call rate. We observed a reduction in the number of phrases of calls (chunks of call with less than 922.3 ms pauses between; Fig. 1I). In addition, MPS IIIB calls were shorter in average duration (Fig. 1J) with longer pauses within a phrase (Fig. 1K), yet were comparable in sound pressure level (Fig. 1L). Because the difference in number of calls and call duration are in the same direction with pause duration in the opposite direction, we can conclude that the reduced call rate is likely not driven by differences in respiration timing. The MPS IIIB calls are also lower in mean frequency pitch (Fig. 1M) with a narrower frequency range (Fig. 1N) and fewer percent of calls including a jump in frequency > 10 kHz (Fig. 1O). These pitch differences may result from a vulnerability in the laryngeal muscles. It is also possible that there is a higher-level disconnect with patterning of calls. Together, these data reveal an early marker of disease in the form of reduced call rate and altered call features that is independent of gross developmental delay.

In a separate cohort of mice, we examined the trajectories of gait development beginning on P21. Gait was measured across four juvenile time points and again in early adulthood (Fig. 2A). During development, MPS IIIB mice were able to run on forced gait apparatus just as well as the controls and did not require an increased number of attempts to run (Fig. 2B; Table 2). All mice increased body weight (Fig. 2C) and body length (Fig. 2D) across the window of assessment (P21—P30) with no significant differences between groups. Because of this drastic change in body size across time, a covariate model was used to regress out the influence of body size change on gait metrics[31]. All mice were tested at the same belt speed (20 cm/s) across development to allow for appropriate comparisons of gait metrics as speed is a robust influence on gait [47, 48]. After correcting for multiple comparisons across 44 gait metrics (Table 3), two remained significantly affected in the MPS IIIB mice. Specifically, the duration of braking for the forepaws was significantly increased (Fig. 2E) and the absolute paw angle of the hindpaws was significantly wider (Fig. 2F). We examined the mice again during adulthood just after P60 to determine if the developmental gait anomalies persisted or resolved with age. By adulthood, body length and body weight were significantly greater in the MPS IIIB mice compared to controls (Fig. 2G,H). At 30 cm/s, MPS IIIB mice were able to run on the apparatus as well as controls (Fig. 2I). Duration to brake for the forepaws was again increased in the MPS IIIB mice (Fig. 2J). Further gait anomalies were apparent at a higher belt speed of 40 cm/s in the adult mice. MPS IIIB mice took significantly more attempts to successfully run on the apparatus compared to controls (Fig. 2I). In addition, this higher speed induced an increased percent of time that both hindlimbs are in stance phase in the MPS IIIB mice (Fig. 2K) as well as decreased percent of stride in swing phase and increased percent of stride in stance phase for the hindlimbs (Fig. 2L). Together, these data indicate that gait was mildly disrupted in the MPS IIIB mice, which persisted into adulthood and worsened when challenged at a higher speed.

Table 3 Gait metrics examined during development and adulthood in MPS IIIB mice compared to control mice

MPS IIIB mice showed significant weight gain during early adulthood with intact visual acuity and contrast

As an estimate of overall health and to identify possible early initiation of decline, the MPS IIIB mice in Cohort 3 were weighed throughout the testing period (Fig. 3A). As compared to controls, the MPS IIIB mice weighed significantly more starting in early adulthood (Fig. 3B; Table 2). This weight phenotype was observed in both females and males (Fig. 3C,D). Clinical evidence of visual acuity deficits have been revealed in some MPS IIIB patients [49]. To assess visual acuity and contrast thresholds, MPS IIIB mice and their littermate controls were tested using the Virtual Optomotor System (VOS). VOS allows for each eye to be tested independently through the innate rodent optomotor reflex. MPS IIIB mice showed comparable mean visual acuity (Fig. 3E) and visual contrast (Fig. 3F) to controls during late adolescence.

MPS IIIB demonstrated reduced startle responses and cued fear responses

To understand if the MPS IIIB mice exhibit normal startle responses and ability to filter irrelevant stimuli through proper sensorimotor gating, we examined these mice in the Acoustic Startle/PPI task. The increased weight phenotype was confirmed by measurement on the force plate used in this apparatus (Fig. 4A), therefore the data were analyzed using an ANCOVA model with weight as a covariate and covariate adjusted means are presented in Fig. 4. MPS IIIB mice demonstrated a decrease in startle response to the 120 dB startle stimulus (Fig. 4B; Table 2), yet show comparable sensitivity to gradually increasing startle stimulus levels (Fig. 4C). A greater inhibition of startle following a much lower dB prepulse tone was also exhibited by the MPS IIIB mice (Fig. 4D). Comparison of responses during non-startle trials to those in startle trials revealed all mice were significantly startled (Fig. 4E). and thus did not exhibit issues hearing the 120 dB full spectrum white noise. These phenotypes were independent of weight and sex, and suggest by late adolescence, the MPS IIIB mice already exhibit features possibly related to reduced fear.

In early adulthood, we examined fear responses directly using the conditioned fear task (Fig. 4F). Fear responses, which have a component of anxiety, are active avoidance behaviors quantified in situations where a threat is imminent and well-defined [50]. We used the fear conditioning task to measure freezing behavior to evaluate associative fear memory in the MPS IIIB mice. MPS IIIB mice showed normal acquisition of fear memory (Fig. 4G), indicating no impaired learning association of a tone with a foot shock. Testing of contextual fear memory indicated MPS IIIB mice had comparable performance relative to control mice (Fig. 4H). Cued fear memory testing revealed a mild impairment in the MPS IIIB mice, which exhibited less freezing on the final day of test (Fig. 4I). Together these results indicate while fear acquisition and contextual memory are similar to controls, cued fear memory was mildly impaired in MPS IIIB mice. Together with the Acoustic Startle/PPI results, these findings may represent early features of reduced fear in this model.

MPS IIIB exhibited intact social approach behavior with sex-specific differences in hierarchy behavior and aggression behavior

To understand sociability and socially-mediated preference behavior in adolescent MPS IIIB mice, we tested these animals in the 3-chamber social approach task on P40 (Fig. 3A, 5A). No differences were observed between groups in the total distance traveled during all four trials (Fig. 5B; Table 2), nor were any side biases observed during the habituation trial (Fig. 5C). During the sociability trial examining time with a novel social partner or an empty withholding cup, MPS IIIB mice showed comparable social preference (Fig. 5D) as well as total time spent interacting with the social stimulus comparable to controls (Fig. 5E). We also examined the social preference index across each minute of the test to see if a previous finding of increased preference during minute 1 in the MPS IIIB mice was replicated [18]. We found a non-significant decrease in social preference index during minute 1 for the MPS IIIB mice (Fig. 5F). During the social novelty trial examining the innate preference for a novel versus familiar social partner, the MPS IIIB mice exhibited comparable novelty preference (Fig. 5G) as well as time investigating the novel partner mouse during the full 10 min trial (Fig. 5H). We also observed a non-significant increase in novelty preference during the first minute of this trial for the MPS IIIB mice (Fig. 5I).

Fig. 5
figure 5

MPS IIIB mice exhibit normal social approach yet sex-specific altered dominance and aggressive behavior. A Schematic of social approach paradigm. B Distance traveled during all four trials of the social approach task. C Stimulus investigation time during full habituation. D-F Total social preference index (D), stimulus investigation time (E) and social preference index across the 10-min sociability trial (F). G-I Total novelty preference index (G), stimulus investigation time (H) and novelty preference index across the 10-min social novelty trial (I). J Schematic of social dominance tube test paradigm. K Percent of dominant bouts for controls and MPS IIIB mice split across males and females. L Schematic timeline of resident intruder paradigm. M Counts of intruder bouts with scramble attacks per group. N Duration of anogenital sniffing by the resident. O Number of total nose sniffs between residents and intruders. B-I Grouped data are presented as means ± SEM with individual data points presented as open circles. K,M Data are presented as counts. N,O Grouped data presented as boxplots with thick horizontal lines respective group medians, boxes 25th – 75th percentiles, and whiskers 1.5 × IQR. Individual data points are presented as open circles

The tube test of social dominance was used to assess the social hierarchy behavior of MPS IIIB mice when paired with sex-matched controls (Fig. 5J). While males show no significant effect of genotype, female MPS IIIB mice were significantly more submissive to control females, demonstrated by only 18.5% of dominant bouts, which was significantly less than chance (50%; Fig. 5K). Finally, the resident intruder paradigm (Fig. 5L) was used to determine if male MPS IIIB mice exhibit increased aggressive behavior. Control mice failed to engage in any scramble attacks with the intruder mouse. However, MPS IIIB mice engaged in scramble attacks with the intruders during 20% of their bouts (Fig. 5M). Examination of non-aggressive social behaviors during this task revealed that MPS IIIB males spent less time engaged in anogenital sniffing of the intruder (Fig. 5N) and fewer nose to nose sniffing bouts with the intruder (Fig. 5O). Together, these data indicate female MPS IIIB mice may be more submissive while males are more aggressive compared to controls.

MPS IIIB mice have larger whole-brain and regional volumes compared to controls with intact tissue integrity

To determine if the MPS IIIB model recapitulates cerebellar atrophy and white-matter thinning reported in patients [4, 5], we performed both structural MR imaging and DTI experiments in cohort 3 animals following completion of behavioral testing at approximately four months of age. The MPS IIIB mice had significantly larger whole-brain volumes at this age compared to their control littermates (Fig. 6A,C-F; Table 2). While the MPS IIIB mice also exhibited larger body weights compared to controls (Fig. 6B), we observed differences in the magnitudes of volumetric differences (females = 5.8%, males = 7.0%) and body-weight difference (females = 13.8%, males = 19.4%) between genotypes. This difference in magnitude of effects suggests the brain size and body size differences may be independent of each other. However, a longitudinal assessment starting from a young age would be needed to confirm this.

Fig. 6
figure 6

MPS IIIB mice have larger whole brain and regional volumes with intact tissue integrity. A Whole brain volumes measured by MRI. B Body weight measured prior to MRI. C-F T2-weighted image and whole brain volume of representative control (C,E) and MPS IIIB mouse brains (D,F). G Ratio of cerebellar volume to whole brain volume. H Mean apparent diffusion coefficient for the cerebellum. I,K Representative T2-weighted image (I) and ADC map (K) of control mouse cerebellum (ADC unit, µm2/ms). J,L Representative T2-weighted image (J) and ADC map (L) of MPS IIIB mouse cerebellum (ADC unit, µm2/ms). M Ratio of corpus callosal volume to whole brain volume. N Fractional anisotropy of the corpus callosum as measured by DTI. O,Q Representative T2-weighted image (O) and FA map (Q) of control mouse ­­corpus callosum. P,R Representative T2-weighted image (P) and FA map (R) of MPS IIIB mouse ­­corpus callosum. S Ratio of cortical thickness measurement to whole brain volume. T Representative schematic of cortical thickness measurements indicated by the red line. Data are presented as means ± SEM with individual data points presented as open circles

Structural regional assessments were conducted by comparing ratios of regional gray matter volumes to whole brain volumes. In the cerebellum, this volumetric ratio was higher in the MPS IIIB mice compared with controls (Fig. 6G,I,J). The structural volumetric ratio also was higher in MPS IIIB mice for corpus callosal volumes (Fig. 6 M,O,P) Next, examining tissue integrity by DTI, the apparent diffusion coefficient for the cerebellum was unchanged between genotypes (Fig. 6H,K,L), as was the fractional anisotropy for the corpus callosum white-matter tract (Fig. 6N,Q,R). Finally, we examined cortical thickness using the T2-weighted images normalized to brain volume. We found comparable cortical thickness ratios for MPS IIIB and control mice (Fig. 6S-U).

To determine if the regional volume differences could be related to behavioral phenotypes, we conducted correlations for acoustic startle/PPI, contextual fear conditioning, and social dominance performance for the subset of cohort 3 mice that received both behavioral phenotyping and imaging. We found a low, negative correlation between the volumetric ratio for the corpus callosum and percent of dominance bouts (Table 4). No other significant relationships between behavioral performance and regional volumetric ratios were detected.

Table 4 Correlations between regional volumetric ratios and behavioral metrics with differential phenotypes across genotypes

Discussion

Characterization of early-onset and clinically-relevant novel phenotypes in the MPS IIIB mouse model will help to identify more sensitive markers to use in evaluation of potential therapies for Sanfilippo Syndrome B. Previous early markers include enhanced degradation of synaptophysin and priming of microglia by postnatal day 10 [51, 52], which require histopathological or biochemical analysis of collected tissue. In this report, we provide evidence of neurobehavioral features in this model suggestive of disease during the first two weeks postnatal as well as evidence of phenotypes during late adolescence and early adulthood.

To our knowledge, this is the first study to examine the trajectory of early communication in the preclinical MPS IIIB model to better understand if this model recapitulates aspects of clinical speech delay [2]. While human language cannot be explored in mice, vocal communication behavior is conserved across taxa [26]. USVs produced by the mouse pup is one of the earliest forms of social communication we can examine in the mouse, and these isolation calls elicit maternal care from the dam [27]. Late in the first week of life and into the second, MPS IIIB mice produced significantly fewer USVs compared to controls as well as alterations to temporal and spectral USV features. It is unknown if the rate of infant crying is diminished in the Sanfillippo population, which would be an interesting future clinical investigation. However, our findings suggest that communication circuits are already perturbed early in development, identifying a very early marker of disease. These very early phenotypes might be leveraged in future treatment studies as markers of in utero interventions, such as enzyme replacement therapy [53, 54].

Along with USV reductions, we also noted a significant decrease in body temperature in MPS IIIB pups compared to controls. Temperature can be a significant modulator of USV rate with colder body temperatures increasing USV rates to elicit maternal care and warmth [46]. However, the temperature effects noted here are in the wrong direction, thus, temperature differences are not driving our USV phenotype. By late adolescence and early adulthood through four months of age, MPS IIIB mice exhibited significantly higher body weights compared to controls, which fall within the expected variability as reported by the Jackson Laboratory (https://www.jax.org/news-and-insights/jax-blog/2020/march/variability-in-aged-mice). However, our results are in contrast to a recent study reporting no differences in weight in MPS IIIB mice compared to controls during the same time period [55]. Hyperactivity has been repeatedly demonstrated in this mouse model [16, 17, 56]. Further, MPS IIIB do not show substantial physical decline and movement impairments until seven to nine months of life, immediately preceding death [57]. Thus it is unlikely the weight gain is a result of inactivity, however, it may reflect other metabolic or energy regulation issues. Previous work on multiple preclinical lysosomal storage disease models, including MPS IIIB, found adipose deficiencies that correlated with reduced leptin levels indicative of altered energy balance [58]. In addition, a recent report found neuronal damage in MPS IIIB mouse brains at five months of age was most notable in thalamic and hypothalamic regions [55]. It is also possible this weight gain is a result of storage material build up. Together, these findings suggest a possible hypothalamic dysfunction or energy regulation disruption, possibly due to the build up of storage material, that is present at very early ages.

In addition, when MR imaged at four months of age, MPS IIIB whole brain volumes as well as cerebellar and callosal volumetric ratios were significantly larger than controls. These were surprising results due to the reports of cerebellar atrophy and white matter abnormalities that accompany behavioral symptomatology [6, 7]. We performed imaging at an age prior to substantial physical decline and movement impairments, which occur around seven to nine months of life, immediately preceding death [57]. MR imaging at a later age nearer this period of decline may reveal features more representative of those observed clinically. Macrocephaly has been reported in a subset of Sanfilippo B patients, although this was coupled with MRI evidence of cortical atrophy in a higher portion of patients [59], which we did not observe here. Our imaging data did not correlate well with behavioral findings. This is consistent with previous study that found the rate of MRI changes were unrelated to the severity of clinical features [6]. However, more work is needed to understand if this is a recapitulation of clinical presentation. Further, it is unclear how larger brain volumes might relate to body size, as we were unable to obtain physical measurements such as femur length in cohort 3. Build up of storage materials occurs in several organ systems including the heart, spinal cord, and brain in both the murine and canine models [10, 60,61,62]. Thus, the larger brain volumes may reflect storage buildup. Further studies coupling structural MRI with quantification of disease-specific heparan sulfate and other storage materials is needed to test this hypothesis.

Loss of acquired functions, including motor functions, is a hallmark of MPS IIIB [63]. While speech decline presents early and frequently, loss of motor skills, including gait, occurs less severely and less frequently. Our data are consistent with these clinical findings. Specifically, we found largely intact developmental trajectories of gait in MPS IIIB mice with abnormalities not present until early adulthood, with the exception of brake duration in the forepaws and absolute hindpaw angle. These two features indicate that the MPS IIIB mice do not increase their speed with which they brake during the stance phase of a stride, a feature of gait maturation we observed previously [31]. We also observed a similar difference from controls for brake duration in two other forms of developmental disabilities genetic risk, specifically a mouse model of Neurofibromatosis Type 1 mutation and the complete deletion mouse model of WIlliams Syndrome [32]. In adulthood, the more robust differences in gait performance were not apparent unless the mice were challenged with a higher belt speed. Previously, we found MPS IIIB mice at 16 weeks of age performed comparably to controls in a standard accelerating rotarod procedure [18], while longitudinal work using a more challenging rocking rotarod observed motor performance declines in the MPS IIIB mice after 200 days of age [16, 17]. Together with our data here, these findings suggest more challenging procedures may be necessary to demonstrate the MPS IIIB motor phenotypes. In addition, gait may be a more sensitive motor domain than rotarod coordination measures to discern earlier motor decline in this model. Joint stiffness and bone abnormalities have been reported clinically [64, 65]. Thus it will be important for future work to investigate possible contributions of joint pathology and bone abnormalities to gait abilities in this mouse model.

Reduced fear is a later onset feature of MPS IIIB [15,16,17,18,19]. We observed mild impairment to cued fear memory in the fear conditioning task in young adult mice (P67). Previous work by others using this model demonstrated intact cued memory around the age studied here (2.5 months) with a deficit presenting by 4.5 months [19]. Li et al., presented the tone once for 10 s accompanied by an 0.5 mA aversive foot shock for the final 2 s, while we presented the tone three times for 20 s with a 1.0 mA aversive foot shock during the final 2 s. Thus, it may be the differences in our procedure allowed us to observe the mild cued fear impairment at this earlier age.

Startle responses, such as those measured in our Acoustic Startle/PPI task, can be potentiated by fear stimuli, such that increasing the fear response increases the amplitude of the acoustic startle response [66]. In addition, startle responses have been positively correlated with anxiety levels [67]. Here, we observed reduced acoustic startle responses in the MPS IIIB mice at P53. This may be an early marker of the clinically-relevant reduced fear phenotype in this model. Further, this may reflect previously described cholinergic disruptions in this model; specifically, decreased acetylcholinesterase activity [18]. Cholinergic neurons of the pedunculopontine tegmentum modulate startle responses [68], thus the reduced startle observed here may be driven by underlying emerging cholinergic pathway abnormalities. We also observed an accompanying increase in inhibition of startle following a prepulse in MPS IIIB mice. Increases in this type of inhibition have been associated with reduced glutamatergic and dopaminergic transmission and increased serotonergic transmission [69,70,71]. Thus, it is plausible these changes in startle inhibition reflect underlying changes to neural activity in the MPS IIIB brain at this early age. Future work in this area will be important for identifying the aberrant circuitry mediating this phenotype.

To determine if this model recapitulates social disruptions reported in clinical populations, we examined social approach, social hierarchy ranks, and interactions in a resident intruder paradigm. Consistent with previous work conducted at 16 weeks of age, we did not see robust changes to sociability in the social approach task [18]. In contrast to that work, MPS IIIB mice here failed to show an initial dip in social preference index. We tested the mice at P40, an age within the peak social reward learning critical period [72] when mice find social interaction highly rewarding. Thus, the increased rewarding value of social interactions at this age may be overcoming any effects of MPS IIIB pathology on social neural circuitry. This hypothesis is supported by observations in another IDD-related genetic risk model, Myt1l haploinsufficiency. Social behavior disruptions in the Myt1l model were observed during adulthood (2–3 months of age) but not during an age in the peak social reward period (P23) [73]. During early adulthood, the MPS IIIB mice showed sex-specific dominance effects. Specifically, female MPS IIIB mice were submissive to controls while the males' social hierarchy behavior was unchanged. Finally, when tested later in adulthood for direct social interactive behaviors, the male MPS IIIB mice demonstrated more aggressive behavior towards the intruders compared to control males, with less non-aggressive social sniffing, thus recapitulating the aggression reported in Sanfilippo B patients.

Conclusions

In the MPS IIIB mouse model, we observed neurobehavioral abnormalities as early as end of first week postnatal, as well as novel phenotypes recapitulating clinically-relevant features including reduced fear, gait disruption, and social behavior changes by late adolescence and early adulthood, respectively. In addition, physical changes and larger brain volumes were observed, potentially reflecting hypothalamic or metabolic dysfunctions and buildup of storage material, respectively, that necessitate further examination. Thus, we identified novel and early physiological and behavioral markers of disease in the MPS IIIB model that will serve as treatment outcome targets for future preclinical evaluations of potential therapeutics, such as gene or enzyme replacement therapies through possible in utero or very early postnatal life interventions [53, 54].

Availability of data and materials

All data and detailed protocols are available upon reasonable request.

Abbreviations

MPS:

Mucopolysaccharidosis

USV:

Ultrasonic vocalization

NAGLU:

Alpha-N-acetylglucosaminidase

VOS:

Virtual optomotor system

CW:

Clockwise

CCW:

Counterclockwise

PPI:

Prepulse inhibition

MRI:

Magnetic resonance imaging

DTI:

Diffusion tensor imaging

References

  1. Valstar MJ, Ruijter GJG, van Diggelen OP, Poorthuis BJ, Wijburg FA. Sanfilippo syndrome: a mini-review. J Inherit Metab Dis. 2008;31:240–52.

    Article  CAS  PubMed  Google Scholar 

  2. Wijburg FA, Węgrzyn G, Burton BK, Tylki-Szymańska A. Mucopolysaccharidosis type III (Sanfilippo syndrome) and misdiagnosis of idiopathic developmental delay, attention deficit/hyperactivity disorder or autism spectrum disorder. Acta Paediatr. 2013;102:462–70.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Shapiro EG, Jones SA, Escolar ML. Developmental and behavioral aspects of mucopolysaccharidoses with brain manifestations — Neurological signs and symptoms. Mol Genet Metab. 2017;122:1–7.

    Article  CAS  Google Scholar 

  4. Malcolm C, Hain R, Gibson F, Adams S, Anderson G, Forbat L. Challenging symptoms in children with rare life-limiting conditions: findings from a prospective diary and interview study with families. Acta Paediatr. 2012;101:985–92.

    Article  CAS  PubMed  Google Scholar 

  5. Hoffmann F, Hoffmann S, Kunzmann K, Ries M. Challenging behavior in mucopolysaccharidoses types I-III and day-to-day coping strategies: a cross sectional explorative study. Orphanet J Rare Dis. 2020;15:275.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Barone R, Nigro F, Triulzi F, Musumeci S, Fiumara A, Pavone L. Clinical and neuroradiological follow-up in mucopolysaccharidosis type III (Sanfilippo syndrome). Neuropediatrics. 1999;30:270–4.

    Article  CAS  PubMed  Google Scholar 

  7. Zafeiriou DI, Savvopoulou-Augoustidou PA, Sewell A, Papadopoulou F, Badouraki M, Vargiami E, et al. Serial magnetic resonance imaging findings in mucopolysaccharidosis IIIB (Sanfilippo’s syndrome B). Brain Dev. 2001;23:385–9.

    Article  CAS  PubMed  Google Scholar 

  8. Bradbury AM, Rafi MA, Bagel JH, Brisson BK, Marshall MS, Pesayco Salvador J, et al. AAVrh10 Gene Therapy Ameliorates Central and Peripheral Nervous System Disease in Canine Globoid Cell Leukodystrophy (Krabbe Disease). Hum Gene Ther. 2018;29:785–801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Dickson PI, Chen AH. Intrathecal enzyme replacement therapy for mucopolysaccharidosis I: translating success in animal models to patients. Curr Pharm Biotechnol. 2011;12:946–55.

    Article  CAS  PubMed  Google Scholar 

  10. Kan S-H, Aoyagi-Scharber M, Le SQ, Vincelette J, Ohmi K, Bullens S, et al. Delivery of an enzyme-IGFII fusion protein to the mouse brain is therapeutic for mucopolysaccharidosis type IIIB. Proc Natl Acad Sci U S A. 2014;111:14870–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Macauley SL, Wong AMS, Shyng C, Augner DP, Dearborn JT, Pearse Y, et al. An Anti-Neuroinflammatory That Targets Dysregulated Glia Enhances the Efficacy of CNS-Directed Gene Therapy in Murine Infantile Neuronal Ceroid Lipofuscinosis. J Neurosci. 2014;34:13077–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Shyng C, Nelvagal HR, Dearborn JT, Tyynelä J, Schmidt RE, Sands MS, et al. Synergistic effects of treating the spinal cord and brain in CLN1 disease. Proc Natl Acad Sci U S A. 2017;114:E5920–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Muschol N, Koehn A, von Cossel K, Okur I, Ezgu F, Harmatz P, et al. A phase I/II study on intracerebroventricular tralesinidase alfa in patients with Sanfilippo syndrome type B. J Clin Invest. 2023;133:e165076.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Deiva K, Ausseil J, de Bournonville S, Zérah M, Husson B, Gougeon ML, et al. Intracerebral Gene Therapy in Four Children with Sanfilippo B Syndrome: 5.5-Year Follow-Up Results. Hum Gene Ther. 2021;32:1251–9.

    Article  CAS  PubMed  Google Scholar 

  15. Cressant A, Desmaris N, Verot L, Bréjot T, Froissart R, Vanier M-T, et al. Improved Behavior and Neuropathology in the Mouse Model of Sanfilippo Type IIIB Disease after Adeno-Associated Virus-Mediated Gene Transfer in the Striatum. J Neurosci. 2004;24:10229–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Heldermon CD, Ohlemiller KK, Herzog ED, Vogler C, Qin E, Wozniak DF, et al. Therapeutic Efficacy of Bone Marrow Transplant, Intracranial AAV-mediated Gene Therapy, or Both in the Mouse Model of MPS IIIB. Mol Ther. 2010;18:873–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Heldermon CD, Hennig AK, Ohlemiller KK, Ogilvie JM, Herzog ED, Breidenbach A, et al. Development of sensory, motor and behavioral deficits in the murine model of Sanfilippo syndrome type B. PLoS ONE. 2007;2:e772.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Kan S-H, Le SQ, Bui QD, Benedict B, Cushman J, Sands MS, et al. Behavioral deficits and cholinergic pathway abnormalities in male Sanfilippo B mice. Behav Brain Res. 2016;312:265–71.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Li HH, Yu W-H, Rozengurt N, Zhao H-Z, Lyons KM, Anagnostaras S, et al. Mouse model of Sanfilippo syndrome type B produced by targeted disruption of the gene encoding α-N-acetylglucosaminidase. Proc Natl Acad Sci U S A. 1999;96:14505–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Escolar ML, Poe MD, Provenzale JM, Richards KC, Allison J, Wood S, et al. Transplantation of umbilical-cord blood in babies with infantile Krabbe’s disease. N Engl J Med. 2005;352:2069–81.

    Article  CAS  PubMed  Google Scholar 

  21. Reddy AS, Kim JH, Hawkins-Salsbury JA, Macauley SL, Tracy ET, Vogler CA, et al. Bone Marrow Transplantation Augments the Effect of Brain- and Spinal Cord-Directed Adeno-Associated Virus 2/5 Gene Therapy by Altering Inflammation in the Murine Model of Globoid-Cell Leukodystrophy. J Neurosci. 2011;31:9945–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Sands MS, Barker JE, Vogler C, Levy B, Gwynn B, Galvin N, et al. Treatment of murine mucopolysaccharidosis type VII by syngeneic bone marrow transplantation in neonates. Lab Invest. 1993;68:676–86.

    CAS  PubMed  Google Scholar 

  23. Sun A. Lysosomal storage disease overview. Ann Transl Med. 2018;6. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6331358/. cited 2020 Nov 23.

  24. Yeager AM, Brennan S, Tiffany C, Moser HW, Santos GW. Prolonged survival and remyelination after hematopoietic cell transplantation in the twitcher mouse. Science. 1984;225:1052–4.

    Article  CAS  PubMed  Google Scholar 

  25. Chen J, Lambo ME, Ge X, Dearborn JT, Liu Y, McCullough KB, et al. A MYT1L syndrome mouse model recapitulates patient phenotypes and reveals altered brain development due to disrupted neuronal maturation. Neuron. 2021;109:3775–3792.e14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Ehret G. Development of Sound Communication in Mammals. In: Rosenblatt JS, Hinde RA, Beer C, Busnel M-C, editors. Advances in the Study of Behavior. Academic Press; 1980. p. 179–225. Available from: http://www.sciencedirect.com/science/article/pii/S0065345408601187. cited 2020 Nov 23.

  27. Haack B, Markl H, Ehret G. Sound communication between parents and offspring. In: Williot, James F., editor. The auditory psychobiology of the mouse. New York, NY: Psychology Press; 1983. p. 57–97. Available from: https://oparu.uni-ulm.de/xmlui/handle/123456789/1201. cited 2017 Apr 24.

  28. Maloney SE, Chandler KC, Anastasaki C, Rieger MA, Gutmann DH, Dougherty JD. Characterization of early communicative behavior in mouse models of neurofibromatosis type 1. Autism Res. 2018;11:44–58.

    Article  PubMed  Google Scholar 

  29. Holy TE, Guo Z. Ultrasonic songs of male mice. PLoS Biol. 2005;3. Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1275525/. cited 2016 May 13.

  30. Rieger MA, Dougherty JD. Analysis of within subjects variability in mouse ultrasonic vocalization: pups exhibit inconsistent, state-like patterns of call production. Front Behav Neurosci. 2016;10. Available from: http://journal.frontiersin.org/article/https://doi.org/10.3389/fnbeh.2016.00182/abstract. cited 2017 Jan 3.

  31. Akula SK, McCullough KB, Weichselbaum C, Dougherty JD, Maloney SE. The trajectory of gait development in mice. Brain and Behavior. 2020;n/a:e01636.

  32. Rahn RM, Weichselbaum CT, Gutmann DH, Dougherty JD, Maloney SE. Shared developmental gait disruptions across two mouse models of neurodevelopmental disorders. J Neurodev Disord. 2021;13:10.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Prusky GT, Alam NM, Beekman S, Douglas RM. Rapid Quantification of Adult and Developing Mouse Spatial Vision Using a Virtual Optomotor System. Invest Ophthalmol Vis Sci. 2004;45:4611–6.

    Article  PubMed  Google Scholar 

  34. Brown JA, Emnett RJ, White CR, Yuede CM, Conyers SB, O’Malley KL, et al. Reduced striatal dopamine underlies the attention system dysfunction in neurofibromatosis-1 mutant mice. Hum Mol Genet. 2010;19:4515–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kopp ND, Nygaard KR, Liu Y, McCullough KB, Maloney SE, Gabel HW, et al. Functions of Gtf2i and Gtf2ird1 in the developing brain: transcription, DNA binding and long-term behavioral consequences. Hum Mol Genet. 2020;29:1498–519.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Kopp N, McCullough K, Maloney SE, Dougherty JD. Gtf2i and Gtf2ird1 mutation do not account for the full phenotypic effect of the Williams syndrome critical region in mouse models. Hum Mol Genet. 2019;28:3443–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Mathis A, Mamidanna P, Cury KM, Abe T, Murthy VN, Mathis MW, et al. DeepLabCut: markerless pose estimation of user-defined body parts with deep learning. Nat Neurosci. 2018;21:1281–9.

    Article  CAS  PubMed  Google Scholar 

  38. Nath T, Mathis A, Chen AC, Patel A, Bethge M, Mathis MW. Using DeepLabCut for 3D markerless pose estimation across species and behaviors. Nat Protoc. 2019;14:2152–76.

    Article  CAS  PubMed  Google Scholar 

  39. Nygaard KR, Maloney SE, Swift RG, McCullough KB, Wagner RE, Fass SB, et al. Extensive characterization of a Williams Syndrome murine model shows Gtf2ird1-mediated rescue of select sensorimotor tasks, but no effect on enhanced social behavior. bioRxiv; 2023. p. 2023.01.18.523029. Available from: https://www.biorxiv.org/content/https://doi.org/10.1101/2023.01.18.523029v1. cited 2023 Mar 16.

  40. Nilsson SR, Goodwin NL, Choong JJ, Hwang S, Wright HR, Norville ZC, et al. Simple Behavioral Analysis (SimBA) – an open source toolkit for computer classification of complex social behaviors in experimental animals. bioRxiv; 2020. p. 2020.04.19.049452. Available from: https://www.biorxiv.org/content/https://doi.org/10.1101/2020.04.19.049452v2. cited 2022 Jul 18.

  41. Kellner E, Dhital B, Kiselev VG, Reisert M. Gibbs-ringing artifact removal based on local subvoxel-shifts. Magn Reson Med. 2016;76:1574–81.

    Article  PubMed  Google Scholar 

  42. Horsfield MA, Jones DK. Applications of diffusion-weighted and diffusion tensor MRI to white matter diseases - a review. NMR Biomed. 2002;15:570–7.

    Article  PubMed  Google Scholar 

  43. Mori S, Zhang J. Principles of diffusion tensor imaging and its applications to basic neuroscience research. Neuron. 2006;51:527–39.

    Article  CAS  PubMed  Google Scholar 

  44. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. 2006;31:1116–28.

    Article  PubMed  Google Scholar 

  45. Laporte MH, Chatellard C, Vauchez V, Hemming FJ, Deloulme J-C, Vossier F, et al. Alix is required during development for normal growth of the mouse brain. Sci Rep. 2017;7:44767.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Harshaw C, Alberts JR. Group and Individual Regulation of Physiology and Behavior: A Behavioral, Thermographic, and Acoustic Study of Mouse Development. Physiol Behav. 2012;106:670–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Jordan K, Challis JH, Newell KM. Walking speed influences on gait cycle variability. Gait Posture. 2007;26:128–34.

    Article  PubMed  Google Scholar 

  48. Möckel G, Perka C, Labs K, Duda G. The influence of walking speed on kinetic and kinematic parameters in patients with osteoarthritis of the hip using a force-instrumented treadmill and standardised gait speeds. Arch Orthop Trauma Surg. 2003;123:278–82.

    Article  PubMed  Google Scholar 

  49. Valstar MJ, Bruggenwirth HT, Olmer R, Wevers RA, Verheijen FW, Poorthuis BJ, et al. Mucopolysaccharidosis type IIIB may predominantly present with an attenuated clinical phenotype. J Inherit Metab Dis. 2010;33:759–67.

    Article  PubMed  PubMed Central  Google Scholar 

  50. La-Vu M, Tobias BC, Schuette PJ, Adhikari A. To Approach or Avoid: An Introductory Overview of the Study of Anxiety Using Rodent Assays. Front Behav Neurosci. 2020;14:145.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Ausseil J, Desmaris N, Bigou S, Attali R, Corbineau S, Vitry S, et al. Early neurodegeneration progresses independently of microglial activation by heparan sulfate in the brain of mucopolysaccharidosis IIIB mice. PLoS ONE. 2008;3:e2296.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Vitry S, Ausseil J, Hocquemiller M, Bigou S, Dos Santos CR, Heard JM. Enhanced degradation of synaptophysin by the proteasome in mucopolysaccharidosis type IIIB. Mol Cell Neurosci. 2009;41:8–18.

    Article  CAS  PubMed  Google Scholar 

  53. Cohen JL, Chakraborty P, Fung-Kee-Fung K, Schwab ME, Bali D, Young SP, et al. In Utero Enzyme-Replacement Therapy for Infantile-Onset Pompe’s Disease. N Engl J Med. 2022;387:2150–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Wang F, Moen DR, Sauni C, Kan S, Li S, Le SQ, et al. Enzyme Replacement Therapy for Mucopolysaccharidosis IIID using Recombinant Human alpha-N-acetylglucosamine-6-sulfatase in neonatal mice. Mol Pharm. 2021;18:214–27.

    Article  CAS  PubMed  Google Scholar 

  55. Petrova R, Patil AR, Trinh V, McElroy KE, Bhakta M, Tien J, et al. Disease pathology signatures in a mouse model of Mucopolysaccharidosis type IIIB. Sci Rep. 2023;13:16699.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Langford-Smith A, Malinowska M, Langford-Smith KJ, Wegrzyn G, Jones S, Wynn R, et al. Hyperactive behaviour in the mouse model of mucopolysaccharidosis IIIB in the open field and home cage environments. Genes Brain Behav. 2011;10:673–82.

    Article  CAS  PubMed  Google Scholar 

  57. Fu H, Kang L, Jennings JS, Moy SS, Perez A, DiRosario J, et al. Significantly increased lifespan and improved behavioral performances by rAAV gene delivery in adult mucopolysaccharidosis IIIB mice. Gene Ther. 2007;14:1065–77.

    Article  CAS  PubMed  Google Scholar 

  58. Woloszynek JC, Coleman T, Semenkovich CF, Sands MS. Lysosomal Dysfunction Results in Altered Energy Balance*. J Biol Chem. 2007;282:35765–71.

    Article  CAS  PubMed  Google Scholar 

  59. Valle DA do, Santos MLSF, Telles BA, Cordeiro ML. Neurological, neurobehavioral, and radiological alterations in patients with mucopolysaccharidosis III (Sanfilippo’s syndrome) in Brazil. Frontiers in Neurology. 2022;13. Cited 11 Aug 2023. Available from: https://www.frontiersin.org/articles/https://www.frontiersin.org/articles/10.3389/fneur.2022.968297.

  60. Pearse Y, Clarke D, Kan S, Le SQ, Sanghez V, Luzzi A, et al. Brain transplantation of genetically corrected Sanfilippo B Neural Stem Cells induces partial cross-correction of the disease. bioRxiv; 2022. p. 2022.06.30.498131. Available from: https://www.biorxiv.org/content/https://doi.org/10.1101/2022.06.30.498131v1. cited 2023 Aug 11.

  61. Le SQ, Kan S, Roberts MS, Dearborn JT, Wang F, Li S, et al. Recombinant NAGLU-IGF2 prevents physical and neurological disease and improves survival in Sanfilippo B syndrome. bioRxiv; 2023. p. 2021.08.06.455469. Available from: https://www.biorxiv.org/content/https://doi.org/10.1101/2021.08.06.455469v4. cited 2023 Aug 11.

  62. Egeland MT, Tarczyluk-Wells MM, Asmar MM, Adintori EG, Lawrence R, Snella EM, et al. Central nervous system pathology in preclinical MPS IIIB dogs reveals progressive changes in clinically relevant brain regions. Sci Rep. 2020;10:20365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Kamp JJPVD, Niermeijer MF, Figura KV, Giesberts MA. Genetic heterogeneity and clinical variability in the Sanfilippo syndrome (types A, B, and C). Clin Genet. 1981;20:152–60.

    Article  PubMed  Google Scholar 

  64. de Ruijter J, Maas M, Janssen A, Wijburg FA. High prevalence of femoral head necrosis in Mucopolysaccharidosis type III (Sanfilippo disease): A national, observational, cross-sectional study. Mol Genet Metab. 2013;109:49–53.

    Article  PubMed  Google Scholar 

  65. Mitrovic S, Gouze H, Gossec L, Schaeverbeke T, Fautrel B. Mucopolysaccharidoses seen in adults in rheumatology. Joint Bone Spine. 2017;84:663–70.

    Article  CAS  PubMed  Google Scholar 

  66. Davis M, Falls WA, Campeau S, Kim M. Fear-potentiated startle: A neural and pharmacological analysis. Behav Brain Res. 1993;58:175–98.

    Article  CAS  PubMed  Google Scholar 

  67. Poli E, Angrilli A. Greater general startle reflex is associated with greater anxiety levels: a correlational study on 111 young women. Front Behav Neurosci. 2015;9:10.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Azzopardi E, Louttit AG, DeOliveira C, Laviolette SR, Schmid S. The Role of Cholinergic Midbrain Neurons in Startle and Prepulse Inhibition. J Neurosci. 2018;38:8798–808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Cano JC, Huang W, Fénelon K. The amygdala modulates prepulse inhibition of the auditory startle reflex through excitatory inputs to the caudal pontine reticular nucleus. BMC Biol. 2021;19:116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Ralph-Williams RJ, Lehmann-Masten V, Otero-Corchon V, Low MJ, Geyer MA. Differential Effects of Direct and Indirect Dopamine Agonists on Prepulse Inhibition: A Study in D1 and D2 Receptor Knock-Out Mice. J Neurosci. 2002;22:9604–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Fletcher PJ, Selhi ZF, Azampanah A, Sills TL. Reduced brain serotonin activity disrupts prepulse inhibition of the acoustic startle reflex. Effects of 5,7-dihydroxytryptamine and p-chlorophenylalanine. Neuropsychopharmacology. 2001;24:399–409.

    Article  CAS  PubMed  Google Scholar 

  72. Nardou R, Lewis EM, Rothhaas R, Xu R, Yang A, Boyden E, et al. Oxytocin-dependent reopening of a social reward learning critical period with MDMA. Nature. 2019;569:116–20.

    Article  CAS  PubMed  Google Scholar 

  73. Kim S, Oh H, Choi SH, Yoo Y-E, Noh YW, Cho Y, et al. Postnatal age-differential ASD-like transcriptomic, synaptic, and behavioral deficits in Myt1l-mutant mice. Cell Rep. 2022;40: 111398.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Steven Le for technical assistance, and Drs. Joshua Dearborn, David Wozniak and Mark Sands for study consultation. Some elements of the figures were created using BioRender.com.

Funding

This work was supported by the NICHD (P50HD103525, IDDRC@WUSTL), the Mallinckrodt Institute of Radiology at Washington University in St. Louis, the Engelhardt Family Foundation, and NIH CTSA (UL1 TR002345, ICTS@WUSTL).

Author information

Authors and Affiliations

Authors

Contributions

KBM, AT, KR, SC, XG, CMY and SEM collected data. KBM and SEM analyzed data. KBM, JRG and SEM wrote the manuscript. JDD, JRG, PD, CMY and SEM conceptualized, designed and oversaw the project. PD, CMY and SEM acquired funding. All authors reviewed and edited the manuscript.

Corresponding author

Correspondence to Susan E. Maloney.

Ethics declarations

Ethics approval and consent to participate

All experimental protocols were approved by and performed in accordance with the relevant guidelines and regulations of the Institutional Animal Care and Use Committee of Washington University in St. Louis and were in compliance with US National Research Council's Guide for the Care and Use of Laboratory Animals, the US Public Health Service's Policy on Humane Care and Use of Laboratory Animals, and Guide for the Care and Use of Laboratory Animals.

Consent for publication

Not applicable.

Competing interests

There are no conflicts of interest to report.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. 

Full spectrogram of ultrasonic vocalization recording for a representative MPS IIIB mouse.

Additional file 2. 

Full spectrogram of ultrasonic vocalization recording for a representative control mouse.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

McCullough, K.B., Titus, A., Reardon, K. et al. Characterization of early markers of disease in the mouse model of mucopolysaccharidosis IIIB. J Neurodevelop Disord 16, 16 (2024). https://doi.org/10.1186/s11689-024-09534-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s11689-024-09534-z

Keywords