What this is
- (AS) is a rare genetic disorder linked to the loss of UBE3A expression, leading to severe neurological deficits.
- This research investigates the electrophysiological and sleep characteristics in a mouse model of AS.
- Key findings include increased seizure susceptibility, elevated delta power, and reduced , which are critical for memory consolidation.
Essence
- Ube3a-del mice exhibit heightened seizure susceptibility, increased delta power, and significant sleep deficits compared to wild-type mice. Notably, they show reduced , which are vital for memory processes.
Key takeaways
- Ube3a-del mice show reduced latency to seizures after a high-dose pentylenetetrazole injection, indicating increased seizure susceptibility compared to wild-type (WT) mice.
- Electroencephalogram (EEG) recordings reveal that Ube3a-del mice have elevated delta power across various frequency bands, suggesting heightened excitability and abnormal brain activity.
- Ube3a-del mice demonstrate significantly less total sleep time and longer latencies to reach paradoxical sleep stages, along with fewer , which are critical for memory consolidation.
Caveats
- The study is limited to the exon 2 deletion mouse model and does not include other potential models, such as the rat model of AS.
- Further research is needed to correlate sleep spindle reduction with cognitive deficits observed in AS, as this study did not assess behavioral outcomes.
Definitions
- Angelman Syndrome: A genetic disorder characterized by severe intellectual disabilities, seizures, and sleep disturbances due to loss of UBE3A expression.
- Sleep Spindles: Short bursts of brain activity during sleep, believed to play a role in memory consolidation and learning.
AI simplified
Introduction
Angelman Syndrome (AS) is a rare (~ 1:15,000) neurodevelopmental disorder (NDD) characterized by impaired expressive communication skills, ataxia, motor and balance deficits, severe intellectual disabilities, recurring seizures and poor sleep [1, 2]. Among the most frequent (> 80%) clinical features of AS are seizures and abnormal encephalogram (EEG) patterns, where seizures often start early in life and are largely resistant to classic anti-epileptic drugs and distinct EEG signatures commonly precede most clinical features and are persistent throughout an AS individualâs lifetime [3, 4]. Seizures occur frequently and typically present across multiple seizure types including, but not limited to, absence, myoclonic and generalized clonic-tonic seizures and, while they are generally diagnosed early in life, they are consistent throughout an individualâs lifetime, contributing to a significantly higher burden of care [5â7]. Abnormal encephalogram (EEG) patterns are widespread in AS and include epileptiform discharges, increased rhythmic delta wave activity and intermittent elevated theta activity [8, 9]. In addition, sleep deficits are also common in AS (20â80%) and are one of the most difficult symptoms to manage, as reported by parents and caretakers [3, 10].
AS results from the loss of expression and function of the ubiquitin-protein ligase E6-AP (UBE3A) gene in neuronal cells [11]. UBE3A is maternally imprinted and, as such, is only expressed from the maternal allele in neuronal but not glial cells in the brain albeit biallelically expressed in all other cell types throughout the body. UBE3A is located on chromosome 15q11.2â13, and its protein product, Ube3a, mainly functions as a ligase responsible for polyubiquitinating chains to substrates, targeting them for degradation by the proteasome [12]. A number of studies indicate a role for Ube3a in synaptic and neural plasticity, which may underlie the imbalance of excitatory/inhibitory homeostasis and contribute to seizure phenotypes and irregularities in sleep [13â15]. Despite a broad understanding of the genetic etiology and some of the basic mechanistic function(s) of Ube3a, there is an unmet need for therapies for individuals with AS.
While we and others have previously described behavioral deficits in rodent models of AS [16, 17] and other neurodevelopmental disorders (NDDs) [18â26], an effort has been made to incorporate in vivo electrophysiology, expanding clinically analogous phenotypes that can be provided as proof of in vivo efficacy. Pursuant to this goal we sought to identify relevant functional phenotypes, including seizures, EEG signature, sleep patterns and sleep spindles using the most rigorously characterized mouse model of AS, with a deletion of Ube3a (Ube3a-del) inherited from the maternal allele resulting in Ube3amâ/p+ on the C57BL/6 J background generated from breeding Ube3am+/pâ females with Ube3am+/p+ males [11].
Many laboratories have used this AS mouse model of Jiang and Beaduet [11] and reported numerous characteristics that resemble core features of AS including seizure susceptibility, increased epileptiform activity, elevated delta and sleep deficits, though they vary widely on seizure induction methods, background strain and EEG collection techniques [15, 27â33]. In corroboration and extension, we quantified seizure susceptibility in Ube3a-del mice on the seizure-resistant C57BL/6J background strain [34]. Additionally, we investigated baseline epileptiform activity and its indication in hyperexcitability metrics. We compared spectral power signatures using our wireless telemetry devices that employed the use of surface electrodes to quantify cortical EEG while allowing free movement of the animal in its home cage. Finally, we extended earlier reports of sleep deficits and, for the first time, quantified sleep spindles in the Ube3a-del mouse model on the C57BL/6J background.
Materials and methods
Subjects
All animals were housed in a temperature-controlled vivarium maintained on a 12:12-h lightâdark cycle. All procedures were approved by the Institutional Animal Care and Use Committee (IACUC) at the University of California Davis and were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals. All experiments were performed on B6.129S7-Ube3atm1Alb/J (Ube3a) mice obtained from The Jackson Laboratory (Stock number 016590; Bar Harbor, ME, USA) and housed in a 24-h lightâdark cycle (7 a.m.â7 p.m.), temperature controlled room and fed a standard diet of Teklad global 18% protein rodent diets (Envigo, Hayward, CA, USA). To maintain the colony, Ube3am+/pâ male mice were paired with C57BL/6J wildtype females resulting in paternal transmission of the mutant allele that is silenced due to imprinting and litters with normal Ube3a expression. To create mice with maternal transmission of the mutant allele, Ube3am+/p+ (WT) male mice are paired with Ube3amâ/p+ females resulting in Ube3a-del mice and their WT littermate controls. To identify mice, pups were labeled by paw tattoo on postnatal day 2â3 using non-toxic animal tattoo ink (Ketchum Manufacturing Inc., Brockville, ON, Canada). At postnatal day 2â7, tails of pups were clipped (1â2 mm) for genotyping, following the UC Davis IACUC policy regarding tissue collection. Genotyping was performed with REDExtract-N-Amp (Sigma-Aldrich, St. Louis, MO, USA) using primers Wildtype Forward: TCA ATG ATA GGG AGA TAA AAC A, Common: GAA AAC ACT AAC ATG GAG CTC and Mutant Forward CTT GTG TAG CGC CAA GTG C.
Design
Three cohorts of mice were used in this study. All subjects were bred in our facility with a deletion of Ube3a (Ube3a-del) inherited from the maternal allele resulting in Ube3amâ/p+ mice generated from breeding Ube3am+/pâ females with Ube3am+/p+ males, termed Ube3a-del hereafter. Cohort 1 consisted of 20 WT and 8 Ube3a-del mice (5 Ube3a-del/9 WT males and 3 Ube3a-del/11 WT females) that were observed for 30 min before and after an administration of pentylenetetrazole (PTZ; 80 mg/kg; i.p.) for behavioral seizure characterization (SKU: P6500, Sigma-Aldrich, St. Louis, MO, USA). EEG data were acquired in Cohorts 2 and 3 where animals were anesthetized and implanted with a wireless telemetry device designed to measure electroencephalogram (EEG) and electromyogram (EMG) in freely moving animals (Data Science International, New Brighton, MN). Cohort 2 consisted of 7 Ube3a-del animals (3 males, 4 females) and 3 WT (2 males, 3 females) that were recorded for 24 h before administration of a lethal dose of pentylenetetrazole (80 mg/kg; i.p.) to observe EEG before and after seizure induction. Cohort 3 consisted of 7 Ube3a-del mice (4 males, 3 females) and 5 WT (4 males, 1 female) that were recorded for 72 h to collect sleep data. All animals were between 8 and 12 weeks old, and experimenters were blind to genotype. Subjects were implanted with the EEG device 7 days prior to data acquisition then, on day 8, began testing. All EEG recordings were collected in the subject animalâs home cage in a temperature-controlled testing room maintained on a 12:12-h lightâdark cycle. All animals were littermates and singly housed after EEG implantation to avoid possible device displacement due to cage-mate interactions.
PTZ administration
Seizure induction studies were conducted using 80Â mg/kg pentylenetetrazole delivered intraperitoneally. Before administration, subjects were observed for 30Â min and weighed to determine the appropriate solution volume. For those that were previously implanted, the weight of the implant (4.0Â g) was subtracted from their total weight. Dosing was conducted in the morning (10:00â11:00) in a dim (~â30Â lx) holding room. Directly after administration of the convulsant, subjects were placed in a clean, empty cage where subsequent seizure stages were live-scored for 30Â min. Seizure stages were scored using latencies to (1) first jerk/Straubâs tail, (2) loss of righting, (3) generalized clonic-tonic seizure, and (4) death. First jerk/Straubâs tail, previously described by Straub et al. (1911), was identified as a tonic dorsal extension of the tail usually accompanied by a jerk or jump of the animalâs entire body. Loss of righting was defined by the absence of both fore- and hindlimb paws from the surface of the cage bottom forâ>â1Â s. Generalized clonic-tonic seizures were identified as loss of righting followed by phases of rigidity and forelimb/hindlimb spasms. Time to each stage was taken in seconds and analyzed across genotype.
EEG implantation
Wireless EEG transmitters were implanted in anesthetized test animals using continuous isoflurane (2â4%). A 2â3Â cm midline incision was made over the skull and trapezius muscles and then expanded to expose the subcutaneous space. Implants were placed in the subcutaneous pocket lateral to the spine to avoid discomfort of the animal and displacement due to movement. Attached to the implant were 4 biopotential leads made of a nickel-cobalt-based alloy insulated in medical-grade silicone, making up two channels that included a signal and reference lead. These leads were threaded toward the cranial part of the incisions for EEG and EMG placement. The periosteum was cleaned from the skull using a sterile cotton-tip applicator and scalpel; then, two 1-mm-diameter burr holes were drilled (1.0Â mm anterior and 1.0Â mm lateral; ââ3.0Â mm posterior and 1.0Â mm lateral) relative to bregma. This lead placement allowed for measurement of EEG activity across the frontal cortical area. Steel surgical screws were placed in the burr holes, and the biopotential leads were attached by removing the end of the silicone covering and tying the lead to its respective screw. Once in place, the skull screws and lead connections were secured using dental cement. For EMG lead placement, the trapezius muscles of the animal were exposed, and each lead was looped through and sutured to prevent displacement. Finally, the incision was sutured using non-resorbable suture material and the animals were placed in a heated recovery cage where they received carprofen (5Â mg/kg; i.p.) directly after surgery and 24Â h post-surgery as an analgesic. Subjects were individually caged with ad libitum access to food and water for 1 week before EEG acquisition and monitored daily to ensure proper incision healing and recovery. Each implantation surgery tookâ<â45Â min, and no fatalities were observed.
EEG data acquisition, processing and analysis
After a 1-week recovery from surgical implantation, individually housed mice were assigned to PhysioTel RPC receiver plates that transmitted data from the EEG implants to a computer via the data exchange matrix using Ponemah software (Data Sciences International, St. Paul, MN). EEG and EMG data were collected at a sampling rate of 500Â Hz with a 0.1 Hz high-pass and 100 Hz low-pass bandpass filter. Activity, temperature and signal strength were collected at a sampling rate of 200Â Hz. Data acquired in Ponemah was read into Python and further processed with a bandpass filter from 0 to 50Â Hz to focus on our frequencies of interest.
Power spectral density analysis
Spectral analysis was performed in Python using MEG and EEG Analysis and Visualization (MNE) open-source software. Frequency bands were defined as delta 0.5â4Â Hz, theta 5â9Â Hz, alpha 9â12Â Hz, beta 13â30Â Hz, and gamma 30â50Â Hz. Spectral power was analyzed using the Welchâs method which windows over the signal and averages across spectral samples. For power spectral densities (PSD) investigated in Cohort 2, analysis started 3Â h into recording and finished 3Â h prior to the end of recording and PTZ administration, resulting in an 18-h sampling window. PSD analysis in Cohort 3 also began 3Â h into recording but continued over the three-day recording resulting in a 69-h sampling window. No statistical difference was detected in PSD within genotype between samples; therefore, both cohorts were combined. Total delta power was determined by adding the density data detected in the 0.5â4Â Hz frequency range, while total power summed all the power spectral density data in the 0.5â50Â Hz frequency range. Relative delta frequencies were calculated by dividing total delta power by total power per animal and averaging across genotype.
Spiking analysis
For spiking analysis, baseline EEG data were segmented into 30-s windows where the mean amplitude was calculated per window. Spiking analysis was conducted in data collected in the 24Â h prior to PTZ administration in Cohort 2 and all of the data collected in Cohort 3. In a first-pass assessment, potential spike data were demarcated as any point 2.5 standard deviations above or below the mean amplitude of a given window. To determine true spike events, the data were then filtered for peaks which were defined as points where the three data points prior to and following the peak were increasing and decreasing in amplitude, respectively, to the potential peak of interest. If activity was detected during a 30-s window, that data were not included in the spike count to avoid possible movement artifact. Similar to PSD analyses, the first and final 3Â h was removed from the spiking data for both cohorts. Spiking activity could not be combined between Cohorts 2 and 3 as the difference in recording time (24 versus 72Â h) greatly contributed to the number of spikes detected.
Sleep analysis
Sleep in mice was assessed using EEG/EMG signals and automatically binned with Neuroscore software (Data Sciences International, St. Paul, MN) into active wake, wake, slow-wave sleep or paradoxical sleep states. A wake state was characterized by a low-amplitude, high-frequency signal with low-EMG tone, while an active wake state was distinguished by high-EMG tone. Sleep was divided into either a slow-wave sleep state or a paradoxical sleep state. Slow-wave sleep was defined by having a high-amplitude, low-frequency signal with elevated delta power and low-EMG tone, while paradoxical sleep had a low-amplitude, low-frequency signal with elevated theta power and low-EMG tone. EEG data were segmented into 1-s windows, and the sleep stage was determined by Neuroscore. We defined a scoring epoch of 10-s and, if at least 50% of the epoch was predominantly one type of sleep stage, that epoch was marked with the majority stage. If a 50% criterion was not reached, then that epoch was not included in analysis. Sleep/wake stages were evaluated in Cohort 3 for the entirety of the acquisition period as this cohort did not conclude with seizure induction. Mean time in a sleep state was calculated by averaging the time spent in each bout of that state. Sleep latency was defined as the average latency to a sleep state from either active wake or wake. Total sleep time was summed across the entirety of the recording from sleep state bouts.
For sleep parameter analysis across lightâdark cycles, the first 24-h time period was sectioned into 2-h time bins starting at 12:00Â a.m. (0â2 time of day). Paradoxical sleep and slow-wave sleep were evaluated separately, while active wake and wake stages were combined into awake readouts. Raw counts of awake, paradoxical sleep, slow-wave sleep or artifact were taken every 10Â s per 2-h time bin, resulting in 720 possible scores. To determine percent time in each sleep stage, the total count of each sleep stage was divided by the total possible count (720) and multiplied by 100. To further quantify both sleep stage counts and percent sleep time between lightâdark cycles, sleep data were combined for the 7â19 time bins and 0â7 and 19â24 time bins. Counts and percent sleep time during the 7â19 time bins were considered the âlight cycle,â while those from the 0â7 and 19â24 time bins were the âdark cycle.â
Sleep spindle analysis
To identify and analyze sleep spindles, we developed a custom Python script modified from a study designed to validate automated sleep spindle detection [35]. Briefly, a bandpass filter with cutoff frequencies of 10 and 15 Hz was applied to include the mouse spindle peak frequency of 11 Hz [36]. Additionally, a Butterworth filter (3 Hz first stopband, 10 Hz first passband, 15 Hz second passband, 22 Hz second stopband, 24 dB attenuation level) was used to further filter for the frequency bands of interest. Next, the root-mean square (RMS) of the filtered signal was calculated with a 750-ms window to smooth the EEG trace before cubing the entire signal to amplify the signalânoise ratio. To detect spindles, a lower threshold (1.2 Ă mean-cubed RMS) was used to determine the start and end of a spindle, while an upper threshold (3.5 Ă mean-cubed RMS) was used to identify the peak of a spindle. Finally, a spindle had to be longer than 0.5 s and less than 10 s for detection. Spindle detection was analyzed for the entirety of the acquisition period of 72 h in Cohort 3.
Statistical analysis

del mice exhibited seizure susceptibility and increased spiking compared to WT during baseline EEG recordings. Seizure susceptibility measures were observed for 30 min before and after an i.p. injection of 80Â mg/kg PTZ.Reduced latencies to first jerk were observed in the-del mice compared to WT.Faster onset to generalized tonicâclonic seizure was observed in-del versus WT. To assess hyperexcitability,mice were implanted with a telemetric device that collected both EEG and EMG data andwas small enough to allow for untethered, unrestrained data collection from the home cage of the test animals.Representative EEG traces of both WT anddel animals sampled before convulsant administration. Quantification of spiking activity during baseline data acquisition in mice recorded for72 h and24 h indicated more spiking events in-del animals compared to their WT littermate controls. *<â0.05, Studentâstest between genotype Ube3a- Ube3a Ube3a Ube3a- Ube3a p t a b c d e, f g h

del mice exhibited elevated delta power. Surface EEG was collected in mice via a wireless telemetric device and analyzed for spectral power dissimilarities in delta (0.5â4Â Hz), theta (5â9Â Hz), alpha (9â12Â Hz), beta (13â30Â Hz) and gamma (30â50Â Hz) frequencies during baseline EEG recordings.,Power spectral density collected fromdel mice differed from WT, specifically in the delta frequency range wheredel mice had higher delta power. No significant differences were detected at theta, alpha, beta or gamma frequencies by genotype.Representative 10-min power distributions ofdel and WT mice illustrated the elevated delta phenotype observed in deletion animals. Power spectral analysis was also evaluated by total delta power (summation of power spectral densities across the delta frequency per animal), relative delta power (total delta power of a given animal divided by that animalâs total power) and total power (summation of power spectral densities across all frequencies per animal) between genotype. Quantification ofdelta power andrelative delta power was significantly increased indel mice compared to WT littermate controls.Total power was also significantly higher indel animals due to their elevated delta power. *<â0.05, Two-way RM ANOVA comparing genotypes across frequency and Studentâstest between genotype Ube3a- Ube3a- Ube3a- Ube3a- Ube3a- Ube3a- p t a c b d e c

Elevation in delta power is persistent across light and dark cycles indel mice. Power spectral densities across all frequency bands were analyzed by lightâdark cycles as well as delta power, relative delta power and total power between genotype.Power spectral density collected fromdel mice during the light cycle differed from WT, specifically in the delta frequency range wheredel mice had higher delta power.Trends toward elevated power spectral densities were also observed indel mice during the dark cycle. Quantification of,delta power and,relative delta power was significantly increased indel mice compared to WT littermate controls across light and dark cycles.,Total power was higher indel animals during the light cycle, with trends toward significance during the dark cycle (=â0.125), due to their elevated delta power. *<â0.05, Two-way RM ANOVA comparing genotypes across frequency and Studentâstest between genotype Ube3a- Ube3a- Ube3a- Ube3a- Ube3a- Ube3a- p p t a e b f c g d h

Sleepâwake cycles are altered indel mice.A schematic detailing how each sleep state is binned using representative samples. Awake state was characterized by a low-amplitude, high-frequency signal with low-EMG tone, while anactive wake state was distinguished by high-EMG tone and activity. Sleep was divided into either a slow-wave sleep state or a paradoxical sleep state.Slow-wave sleep was defined by having a high-amplitude, low-frequency signal with elevated delta power and low-EMG tone, whileparadoxical sleep had a low-amplitude, low-frequency signal with elevated theta power and low-EMG tone.del mice had significantly reducedmean time in paradoxical sleep andincreased latencies to paradoxical sleep states, suggesting that deep sleep in deletion mice was not only shorter, but delayed compared to WT littermate controls. Additionally,del mice exhibited decreased mean time inslow-wave sleep that when combined with the paradoxical sleep deficits resulted in lowertotal sleep time compared to WT controls. *<â0.05, Studentâstest between genotype Ube3a- Ube3a- Ube3a- p t a (i) (ii) (iii) (iv) b d c h

Sleep deficits observed indel mice are persistent across lightâdark cycles. Sleep parameters were analyzed across lightâdark cycles by sectioning the first 24-h time period into 2-h time bins starting at 12:00Â a.m. (0â2 time of day). Time in each sleep stage, represented as percent of a 2-h time bin, and bouts of sleep stages, represented as counts per 2-h time bin, were examined. Both wake and active wake stages were combined for awake analysis. Grey boxes highlighting hours 7â19 indicate the dark cycle, while unhighlighted hours ranging from 0 to 7 and 19â24 indicate light cycle.A trend toward increased percent time awake was detected indel mice compared to WT littermate controls (=â0.105),that was not specific to either the light or dark cycle.Awake bouts were not significantly different betweendel mice compared to WT, andno significant differences were observed between light and dark cycles. Interestingly,-del mice exhibited less percent paradoxical sleep time compared to WT controls that wassignificant across both the light and dark cycles.Likewise, paradoxical sleep bouts were significantly reduced in-del mice compared to WT littermate controls, againacross both light and dark cycles. Both percent time () and bouts () of slow-wave sleep were not significantly different between genotypes or across lightâdark cycles (,). *<â0.05, Two-way RM ANOVA comparing genotypes across time bins and mixed model effects between genotype and lightâdark cycle Ube3a- Ube3a- p Ube3a- Ube3a Ube3a p a b c d e f g h i k j l

del mice exhibited reduced sleep spindle production that was significantly correlated with delta rhythmicity. Sleep spindles were automatically detected using a custom written Python script.del mice had significantly less sleep spindles when compared to their wildtype littermate controls. When taking spindles as a fraction of total sleep time,-del animals had more spindle counts given their reduced sleep time, suggesting increased frequency of spindle counts during their infrequent sleep periods. To further analyze the-del spindle count by reported sleep deficitscorrelational analyses were run between spindles and total sleep time, mean time in slow-wave sleep, mean time in paradoxical sleep and latency to paradoxical sleep. No significant correlations were detected between spindle count and sleep. Additionally,spindle production was correlated with spiking events to see if elevated spiking indel mice could predict spindle count, but no relationship was found. Finally,spindle count was negatively correlated with delta power indel mice, where animals with enhanced delta had less spindles than animals with lower delta power. *<â0.05, Studentâstest between genotype and linear regression for correlation across behaviors Ube3a- Ube3a- Ube3a Ube3a Ube3a- Ube3a- p t a b câf f g
Results
Increased seizure susceptibility and spiking activity indel mice. Ube3a-
Seizure susceptibility was evaluated in Cohort 1 by latency to (1) first myoclonic jerk and/or Straubâs tail, (2) loss of righting reflex, (3) generalized clonic-tonic seizure and (4) death after administration of PTZ (80 mg/kg; i.p.). A reduction in latency indicated susceptibility, while an increase in latency indicated resistance. Ube3a-del animals exhibited seizure susceptibility via reduced latency to first jerk and generalized clonic-tonic seizure (Fig. 1a, b: t(26) = 2.287, p = 0.031; t(26) = 2.627, p = 0.014), compared to WT. No differences were detected in loss of righting reflex or death (Additional file 1: Fig. 1a, b: t(26) = 1.461, p = 0.156; t(26) = 1.333, p = 0.194). Hyperexcitability and seizure susceptibility were further analyzed by implanting Cohort 2 subjects with a wireless telemetric device that captured continuous EEG (Fig. 1c, d). Epileptiform activity, such as spiking, is sufficient in the detection and diagnosis of epilepsy. As expected, Ube3a-del mice in both Cohort 2 and Cohort 3 exhibited increased spiking activity during baseline recording compared to WT controls (Fig. 1g, h: t(10) = 3.435, p = 0.006; t(8) = 3.132, p = 0.014).
Elevated delta power indel mice Ube3a-
Rhythmic delta activity is a consistent spectral signature of AS [37, 38] and has been previously reported in mouse models [31, 32], but not via a single-channel electrode placed in the skull via a wireless telemetric system. To test whether our methodology was able to capture EEG signal sensitivity that would also detect increased delta activity, we evaluated spectral power dynamics of WT and Ube3a-del mice. When comparing power spectral densities across 0â50 Hz frequencies, Ube3a-del mice displayed robust elevated delta power (Fig. 2a: F(1, 20) = 6.432, p = 0.020; t(20) = 2.763, p = 0.012 for delta). When densities were binned into delta, theta, alpha, beta and gamma frequency groups, delta power in Ube3a-del mice was significantly higher than WT littermate controls (Fig. 2c: F(1, 20) = 5.862, p = 0.025; t(20) = 2.763, p = 0.012 for delta).
Delta-specific dynamics were further examined by looking at total delta, relative delta and total power across the entire recording. Ube3a-del mice exhibited higher total and relative delta power compared to WT across the entire recording (Fig. 2d, e: t(20) = 2.763, p = 0.012; t(20) = 4.089, p = 0.0001). Interestingly, total power was also higher in Ube3a-del animals, likely a result of the robust increase in delta (Fig. 2f: t(20) = 2.564, p = 0.019).
As delta activity is elevated in sleep stages, particularly slow-wave sleep, delta activity was also analyzed across lightâdark cycle to control for the possible influence of sleepâwake cycles. When looking at the light cycle, elevated delta activity was detected across all 0â50 Hz power spectral densities (Fig. 3a: F(1, 20) = 8.263, p = 0.009; t(20) = 3.020, p = 0.007 for delta). Furthermore, both relative and total delta power were significantly higher in Ube3a-del animals as well as total power compared to WT (Fig. 3bâd: t(20) = 3.042, p = 0.006; t(20) = 2.278, p = 0.034; t(20) = 2.989, p = 0.007). Similarly, in the dark cycle, Ube3a-del mice displayed elevated delta activity (Fig. 3e: t(20) = 2.178, p = 0.042 for delta), with trends toward significantly increased PSDs across the 0â50 Hz frequency bins (Fig. 3e: F(1, 20) = 4.027, p = 0.059). Unsurprisingly, higher total and relative delta power was observed in Ube3a-del animals; however, no change in total power was detected (Fig. 3fâh: t(20) = 2.191, p = 0.041; t(20) = 2.823, p = 0.011; t(20) = 1.60, p = 0.125).
Sleep deficits and abnormal sleepâwake cycles detected indel mice Ube3a-
Sleep data were acquired in Cohort 3 animals and subsequently parsed into 4 distinct sleepâwake stages: active wake, wake, slow-wave sleep, and paradoxical sleep. A wake state was characterized by a low-amplitude, high-frequency signal with low-EMG tone, while an active wake state was distinguished by high-EMG tone. Sleep was divided into either a slow-wave sleep state or a paradoxical sleep state. Slow-wave sleep was defined by having a high-amplitude, low-frequency signal with elevated delta power and low-EMG tone, while paradoxical sleep had a low-amplitude, low-frequency signal with elevated theta power and low-EMG tone. Ube3a-del mice displayed reduced mean time in paradoxical sleep (Fig. 4b: t(10) = 2.91, p = 0.016) and took longer to reach paradoxical sleep stages compared to WT (Fig. 4d: t(10) = 3.691, p = 0.004). Furthermore, trends toward reduced mean time in slow-wave sleep were detected in Ube3a-del mice (Fig. 4d: t(10) = 1.931, p = 0.0823), contributing to a significantly reduced total sleep time compared to WT littermate controls (Fig. 4e: t(10) = 2.741, p = 0.021).
Sleep parameters were further analyzed by lightâdark cycles. The first 24-h time period was sectioned into 2-h time bins starting at 12:00 a.m. (0â2 time of day) where counts of each sleep stage were scored as well as the percent time of a given sleep stage. Both wake and active wake stages were combined for awake analysis. A trend toward increased percent awake time was detected in Ube3a-del mice compared to WT littermate controls (Fig. 5a: F(1,10) = 3.173, p = 0.105) that was consistent across both the light and dark cycles (Fig. 5b: F(1, 6) = 4.105, p = 0.089). No change in awake counts was observed across time bins or lightâdark cycles between genotypes (Fig. 5c: F(1,10) = 0.465, p = 0.511; Fig. 5d: F(1,6) = 465, p = 0.521). Interestingly, Ube3a-del animals had reduced percent paradoxical sleep time across time bins (Fig. 5e: F(1,10) = 12.72, p = 0.005) that was observed for both the light and dark cycle compared to WT littermate controls (Fig. 5f: F(1,6) = 28.02, p = 0.002; WT and Ube3a-del by light cycle p = 0.009; WT and Ube3a-del by dark cycle p = 0.014). Expectedly, Ube3a-del mice showed decreased paradoxical sleep counts across time bins (Fig. 5g: F(1,10) = 12.71, p = 0.005) that was reduced in both light and dark cycles compared to WT (Fig. 5h: Genotype effect F(1,6) = 28.01, p = 0.002; WT and Ube3a-del by light cycle p = 0.009; WT and Ube3a-del by dark cycle p = 0.014). Finally, the percent time in slow-wave sleep was not significantly different between genotypes (Fig. 5i: F(1,10) = 0.399, p = 0.542), or across lightâdark cycle (Fig. 5j: F(1,6) = 1.310, p = 0.296). No significant differences in slow-wave sleep counts were detected across both time bins (Fig. 5k: F(1,10) = 0.397, p = 0.543) and between lightâdark cycles (Fig. 5l: Genotype effect F(1,6) = 1.302, p = 0.297).
Sleep spindles are reduced indel mice and negatively correlated with the elevated delta phenotype Ube3a-
We wanted to identify, validate and quantify sleep spindles in the Ube3a-del mouse, adding another clinically relevant functional phenotype for therapeutic testing. First, a 10â15 Hz bandpass filter was first applied, followed by a Butterworth filter (3 Hz first stopband, 10 Hz first passband, 15 Hz second passband, 22 Hz second stopband, 24 dB attenuation level) for the frequency bands of interest. Next, the root-mean square (RMS) of the filtered signal was calculated with a 750-ms window to smooth the EEG trace before cubing the entire signal to amplify the signalânoise ratio. To detect spindles, a lower threshold (1.2 Ă mean-cubed RMS) was used to determine the start and end of a spindle, while an upper threshold (3.5 Ă mean-cubed RMS) was used to identify the peak of a spindle. Ube3a-del mice exhibited less sleep spindles total compared to WT littermate controls (Fig. 6a: t(10) = 2.357, p = 0.04). To control for sleep time we first divided spindle count by total sleep time. Curiously, Ube3a-del mice had significantly higher spindle production during their shorter sleep periods compared to WT controls (Fig. 6b: t(10) = 2.361, p = 0.04). As additional controls for sleep, spindle count was correlated with total sleep time (Fig. 6c: F(1,3) = 1.322, p = 0.334), mean time in slow-wave sleep (Fig. 6d: F(1,3) = 4.723, p = 0.118), mean time in paradoxical sleep (Fig. 6e: F(1,3) = 5.422, p = 0.102) and latency to paradoxical sleep (Fig. 6f: F(1,3) = 0.913, p = 0.410) in Ube3a-del animals. No significant correlations were detected between spindle count and any other sleep metric. Spindle count and spiking count correlations were also analyzed, though no significant correlation was detected (Fig. 6h: F(1,3) = 0.077, p = 0.800). Interestingly, spindle count was negatively correlated with elevated delta power, where animals with higher delta power had lower spindle counts (Fig. 6g: F(1,3) = 59.44, p = 0.005). Additional correlation studies showed no significant relationship between spiking, delta and sleep phenotypes (Additional file 2: Fig. 2aâj).
Discussion
Novel therapies in development for genetic precision medicine for AS that could be âcurativeâ have resulted in (1) unsilencing of paternal Ube3a; (2) molecular reversal of Ube3a expression levels; and (3) some degree of functional phenotypic rescue by a wide variety of molecular therapies. These include gene therapy by antisense oligonucleotides (ASO) [39], viral vector delivery [40], and artificial transcription factors (ATFs) [41]. In fact, two ASO compounds are in Phase I clinical trials (GeneTx NCT04259281â; Roche NCT04428281â). Outcome measures are required to demonstrate the utility of these innovative therapeutic designs as well as to validate other traditional medicinal therapies that may be in the drug discovery pipeline by biotechnology and pharmaceutical companies for AS. Previous research has comprehensively characterized the AS mouse model line of Jiang and Beaduet both behaviorally and biochemically and discovered motor abnormalities [42], learning and memory deficits in fear conditioning [39, 40], strain-dependent seizure susceptibility [15, 32], elevated delta spectral power [31] and abnormal sleep signatures [43]. This study rigorously reproduced findings on hyperexcitability and seizure susceptibility and extended studies using a unique chemoconvulsant mechanism, GABAB antagonist, pentylenetetrazole. Innovation reported herein is the detection of similar effect sizes and phenotypes using a translational approach of acquisition of EEG in the home cage over several days via a skull screw and wireless telemetry over hippocampal depth recordings of local field potentials requiring a head mount and tethered system. Previous work from Sahin and Rotenberg has established this technology as useful for models of NDDs, by their extensive studies in Tuberous Sclerosis Complex and PhelanâMcDermid Syndrome models [24, 44]. Both of these genetic NDDs share common phenotypes with AS, including motor difficulties, intellectual disabilities, recurring seizures and sleep difficulties, in addition to diagnoses of autism spectrum disorders (ASD).
It is interesting, although not surprising, that no behavioral seizures were observed in either the mouse or rat model of AS, given it has been recognized by our laboratory and others that mouse models on C57BL6/J backgrounds are resistant to seizures [22, 34]. Since multiple genetic NDDs that are syndromic forms of ASD have high seizure comorbidity, our laboratory adapted to uncovering subthreshold behavioral seizures with wireless, untethered, telemetry implants that acquire the EEG signal. This is the first report to identify with translationally relevant methods, global (non-hippocampal; non-depth local field potentials; freely moving non-tethered) neurophysiology and sleep in AS mice.
In addition to observations of epileptiform and spike trains via EEG in AS subjects, we observed elevated spectral power in the delta frequency band in mice with Ube3amâ/p+ deletions. Our laboratory works closely with patient advocacy groups and gathers observations of clinical AS from leading epileptologists (Thibert, Anderson) providing us numerous deidentified clinical examples for which EEG signatures in individuals with mutations or excess expression of Ube3a are unique and distinguishable. For example, Drs. Thibert and Anderson have carefully explained their reports of elevated delta power in AS clinics [31, 32]. Further, we reproduced Sidorov et al. earlier work [45] that highlighted delta power as translational biomarker [31, 33]. This phenotypic reliability paves the pathway to power spectral signatures as therapeutic windows for precision treatment within AS and beyond including other NDDs, such as Dup15q and genetic forms of ASD [46â49].
Sleep is greatly affected in individuals with AS [50â52]. Sleep disturbances reported include reduced overall time sleeping, higher number of nighttime awakenings, and longer onset latencies to falling asleep are the most common [53, 54]. Sleep analysis is highly translational, since a majority of genetic NDDs, including AS, have clinical sleep disruption. This report defined rodent sleep cycles over a 36-h acquisition period. We measured alterations by quantifying time in: (a) wake, defined as high-frequency, low-amplitude signals with without EMG or video activity; (b) active wake, similar to wake but with detected activity; (c) time in slow-wave sleep, defined as low-frequency, high-amplitude signal with elevated delta; (d) time to sleep onset and time spent in paradoxical (~ REM) sleep, characterized by a low-frequency, low-amplitude signal with elevated theta.
We reproduced Ehlen and Philpotâs earlier work [45] that maternal Ube3a loss had a striking effect on the architecture of sleepâwake cycles. As reported, Ube3a-del had reduced REM sleep bouts when employing 2-h time bins to establish percent time scores in Ube3a-del mice compared to WT controls during the dark phase. We reproduced these findings by our discoveries of reduced percent time in paradoxical sleep using 2-h time bins and counts/10 s of paradoxical sleep. In contrast, we observed these deficits across both light and dark cycles. Methodologies and coding algorithms may account for the differences. We are continuing to examine this discrepancy. This earlier work also described shorter NREM bout durations throughout the dark cycle in Ube3a-del mice, and specifically NREM frequency of bouts was shorter at the beginning of the dark cycle and then increased over time. The number of NREM bouts also remained high during the early dark in Ube3a-del mice. These results concluded that Ube3a-del mice had a âfragmentedâ NREM sleep in the dark cycle compared to WT mice. While we did not observe changes in percent slow-wave sleep time or count across our 24-h sampling window, we did observe trends toward reduced slow-wave sleep time across the total recording. Furthermore, Ube3a-del mice showed impairments in total sleep time, namely through reductions in mean paradoxical sleep time and trends toward reduced slow-wave sleep compared to WT animals. Ube3a-del mice took longer to reach paradoxical sleep compared to WT, suggesting difficulty reaching deep sleep stages and, once there, difficulty remaining in those stages. We believe our data corroborate the impaired sleep phenotype described by Ehlen and Philpot although our signal collection methods and output metrics differed.
For the first time, in a preclinical model of AS, we have defined sleep spindles with custom-built automation, following manually filtering data and processed via our custom algorithm. Sleep spindles are thalamocortical oscillations ranging from 11 to 16 Hz and are theoretically thought to mediate memory consolidation. There have been reports of sleep spindle deficits in the AS population, but not in a mouse model of AS [43]. Ube3a-del mice exhibited less spindles, as we hypothesized given the reduction in sleep spindles observed clinically in AS [43]. Interestingly, when accounting for total sleep time, spindle count was significantly higher in Ube3a-del mice compared to WT, but not directly correlated with total sleep time, mean time in paradoxical sleep and mean time in slow-wave sleep. From these data we posit that spindle production is abnormal in Ube3a-del as observed by reduced total amount of spindles and increased spindle frequency during shorter sleep periods. Future studies will be required to investigate the nuances of these phenotypes.
Fewer, abnormal sleep spindles have been historically implicated in intellectual disabilities [55, 56] and for numerous genetic NDDs, such as PhelanâMcDermid and PraderâWilli syndromes [57, 58] as well as childhood epilepsies [59] 60 and ASDs [61, 62]. Sleep spindles have been reported altered when mGluR5 is deficient [63], and mGluR5 dysregulation has been postulated as underlying Fragile X Syndrome phenotypes, another NDD with ASD and poor sleep regulation [64, 65]. Interestingly, there was a significant negative correlation between spindles and delta power, where animals with higher delta power tended to have lower spindle counts. It is unclear whether elevated delta power directly correlates with any core clinical features of AS, but these data offer a promising link between delta rhythmicity and possible cognitive outcome measures.
Limitations to this work were our lack of inclusion of anatomical substrate analysis either by gross MRI scan or histopathology in regions to identify neuronal integrity or lack thereof. While less regionally specific than typical depth recording techniques, we sought to use a more clinically analogous EEG recording method that collected signals from the surface of the skull and was wireless to allow for free movement of test animals. We hypothesized that quantifying EEG measures in subjectsâ home cage environment using wireless telemetry devices would be sufficient to assess seizures, sleep deficits and spectral band abnormalities associated with models of AS and NDDs as well as offer translational observation of global neuronal activity, as previously described [44]. Further limitations were our lack of correlation of the reduction of spindles with lower cognitive behavior observed by others in AS mice, a popular theory of spindle function [55, 56].
In summary, our data enhance rigor and translatability of EEG readouts as having biomarker potential for preclinical testing of therapeutics. Our study provides important corroboration of previous reports of rodent epileptiform with EEG that illustrated that Ube3a-del mice had increased epileptiform spiking activity and elevated delta power, which corroborates reported work and recapitulates clinical reports in AS. This is the first report to use a cortical surface-based recording by a wireless telemetry device over tethered/fixed head-mount depth recordings. Less time in both paradoxical and slow-wave sleep, longer latencies to paradoxical sleep stages and total less sleep time in Ube3a-del mice were observed compared to WT. For the first time, we detected fewer sleep spindles in the AS mouse model, a critical marker of memory consolidation during sleep, in a preclinical model of AS. This study was limited to the exon 2 deletion mouse model, and future work will investigate the rat model of AS, containing a complete Ube3a deletion and pair EEG with behavior. Our data enhance rigor and translatability of EEG, spectral power and sleep as beneficial outcomes for therapeutic testing.
Supplementary Information
Additional file 1: Supplemental Figure 1. Ube3a-del mice exhibited similar latencies to loss of righting and death after convulsant administration. Seizure susceptibility measures were observed for 30 min after an i.p. injection of 80 mg/kg PTZ. While reduced latencies to both first jerk and generalized clonictonic seizures were observed in Ube3a-del mice Figure 1 (a, b), no genotype differences were detected in either A latency to loss of righting or B death. *p < 0.05, Studentâs t-test between genotype.Additional file 1: Supplemental Figure 2. No significant correlations were detected across spiking events and delta rhythmicity, sleep metrics and delta rhythmicity, or spiking and sleep in Ube3a-del mice. Spiking events and delta power were correlated in A Cohort 3 Ube3a-del mice and B Cohort 2 Ube3a-del mice, with no significant relationship detected. Similarly, when analyzing CâF mean time in paradoxical sleep, mean time in slow-wave sleep, latency to paradoxical sleep and total sleep time with delta power and GâJ mean time in paradoxical sleep, mean time in slow-wave sleep, latency to paradoxical sleep and total sleep time with spiking events in Cohort 3 animals, no significant correlations were found. *p < 0.05, linear regression for correlation across behaviors.