Severity of pilocarpine-induced DSE is independent of pup genotype
Pilocarpine injection at P12 caused diarrhea and seizures accompanied by automatic scratching, drumming motions, tremors, limb extensions, postural problems and wild running, which is consistent with signs described in the literature [47]. The signs appeared in characteristic order during seizure onset, and the time of their first occurrence was independent of pup genotype (Fig. 1a). Seizure severity was assessed by visual observation complemented by video recordings at two 15-min observation intervals, the first starting immediately after injection, and the second starting 4 hours after the injection. The signs were scored using a modification of the Racine scale. Clonic behaviors, such as shivering or forelimb clonus (“drumming “) counted as stage 3, while tonic behaviors (general stiffness, tail clonus, falls or failure of righting reflex) counted as stage 5. Behaviors which could represent both seizure-related automatisms and normal behavior (i. e. scratching) only counted as stage 1. The highest stage observed in the 15-min interval for each animal was then counted. Seizure severity was comparable in Tsc2+/- and wild type pups at both time points (Fig. 1b). This suggests that genotype did not play a role in the severity of the seizure. After 5–6 hours, the symptoms receded in most individuals, and the pups were then returned to the mothers. No mortality was observed.
Social behavior is altered by Tsc2+/- and pilocarpine-induced DSE during early development
We monitored body weight of the animals from DSE induction to adulthood (P11–P90). We found an obvious effect of age (F(1.261, 83.225) = 4207.379; p < 0.001; Ƞp2 = 0.985), but also an interaction between age and DSE treatment (F(1.261, 83.225) = 4.092; p = 0.037; Ƞp2 = 0.058). Genotype had no influence on weight. Simple effects showed significant differences between DSE and naïve groups at P14 (p < 0.001); P19 (p = 0.001); P27 (p = 0.015) and P40 (p = 0.010), with DSE animals exhibiting lower body weight. This indicates that DSE affected physical development of the pups, with persistent retardation of growth only normalizing at early adulthood (Fig. 1c, 1d). Moreover, the effects of DSE on social behavior at P19 were examined. Total duration of social play was selected as the parameter most relevant to eventual autism-like deficits. The data were Log-transformed because of non-normal distribution, and analyzed by t-test. In pairs consisting of DSE pups, social play was much less common: T(15) = 3.15, p = 0.007. Also, ultrasonic vocalizations (total call count) were much less common in DSE pups: T(15) = 3.035, p = 0.008 (Fig. 1e-g). Together, this confirms that DSE has a lasting impact on both physical and social development of the pups.
In experimentally naïve P7 pups (prior to DSE), we observed an altered isolation-induced ultrasonic vocalization (USV) profile in Tsc2+/- pups, suggesting that Tsc2+/- also leads to early changes in vocal communication patterns (for detailed description, see Supplementary results, Fig. S2 and Suppl. Tab. S1).
Everolimus reduces social behavior deficits in adult rats with Tsc2+/- but not after DSE
Young adult rats from all four experimental groups (wild-type naïve, Tsc2+/- naïve, wild-type DSE, Tsc2+/- DSE) were investigated in several behavioral paradigms under three different conditions: baseline before Everolimus treatment (T1), during treatment with Everolimus (T2), and after wash-out (T3).
In the social interaction test, the parameters of anogenital and non-anogenital social exploration were chosen for analysis (Fig. 2b, 2c, Supplementary Tab. S1). In anogenital social exploration, three-way ANOVA of T1-T3 revealed an effect of DSE, F(1, 32) = 23.909; p < 0.001; Ƞp2 = 0.428, with rats with a history of DSE being less explorative. There was also an effect of time F(2, 64) = 16.078; p < 0.001; Ƞp2 = 0.334. Importantly, there was a significant interaction between Tsc2+/- and time, F(2, 64) = 3.438; p = 0.038; Ƞp2 = 0.097. More detailed examination by simple effects analysis showed that Tsc2+/- genotype significantly decreased anogenital exploration at T1 (p = 0.036), but this deficit was no longer apparent at T2 and T3 (Fig. 2b). This indicates an effect of Everolimus on anogenital social exploration impaired by Tsc2+/-, but not by DSE. In non-anogenital social exploration, three-way ANOVA of T1-T3 was only significant for the effect of DSE F(1, 32) = 12.285; p = 0.001; Ƞp2 = 0.277. Two-way ANOVA of T1, T2 and T3 showed an effect only of DSE (Fig. 2c). Other parameters either did not differ between groups, or were too rare for analysis, and are not shown.
Locomotor activity in the EPM, measured by total arm visits, was affected neither by Tsc2+/- nor by DSE. In three-way ANOVA, the effect of time was significant, F(3,96) = 3.839; p = 0.012; Ƞp2 = 0.107, corresponding to a gradual decrease of activity across sessions, probably driven by habituation of the animals. (Fig. 3a). Open arm visits exhibited strongly non-normal distribution, which was not improved by data transformation. Therefore, we used Negative binomial model with log link function, and did not found a statistically significant effect of either DSE or genotype at T1, T2 or T3. However, number of open arm visits was decreased in Tsc2+/- animals in the habituation session: 0.168 (95% CI, 0.029 to 0.972) p = 0.046 (Fig. 3b). In looking down duration indicating risk assessment behavior, three-way ANOVA showed a significant effect of DSE, F(1, 32) = 7.780; p = 0.009; Ƞp2 = 0.196, with DSE rats spending more time looking down from the maze than their naïve counterparts. On the other hand, this behavior was reduced in Tsc2+/- rats, F(1, 32) = 6.205; p = 0.018; Ƞp2 = 0.162. There was also an effect of time, F(2.168, 69.377) = 3.71; p = 0.026; Ƞp2 = 0.104, suggesting gradual decrease. There was no interaction between the factors. Two-way ANOVA for individual sessions confirmed the effects of DSE at all time points. Effects of Tsc2+/- were present at habituation and T3, but not at T1 and T2 (Fig. 3c).
Locomotion of adult rats in the OF was neither affected by Tsc2+/- nor DSE. In three-way ANOVA, there was an effect of time (F(3, 87) = 18.233; p < 0.001; Ƞp2 = 0.386), again suggesting a gradual decrease of locomotion either due to habituation or age-dependent changes in activity (Fig. 3d).
In LD, the time spent in the light compartment, exhibited non-normal distribution and the data had to be log-transformed. In three-way ANOVA, there was an effect of time, F(1.443, 41.857) = 5.065; p = 0.019; Ƞp2 = 0.149, suggesting gradual decrease of anxiety in time. There was a significant interaction between DSE and time, F(1.443, 41.857) = 4.329; p = 0.030; Ƞp2 = 0.130. However, subsequent simple effect analysis didn’t find a specific significant difference on any testing day, only a trend at T3 (p = 0.071). Two-way ANOVA for individual sessions did not find specific differences at T1 and T2, however, animals that underwent DSE did spend less time in the light compartment at T3 (Fig. 3e). Bonferroni’s post-hoc tests showed that it was due to the Tsc2+/- DSE group exploring the light compartment significantly less than any other group (Supplementary Tab. S1).
In the Y-maze, locomotor activity indicated by the total number of arm visits, was only affected by time, F(2, 64) = 23.895; p < 0.001; Ƞp2 = 0.427, with gradual decrease evidencing habituation to the environment. Working memory, measured by spontaneous alternation, was decreased in DSE animals, according to three-way ANOVA, F(1, 32) = 4.317; p = 0.046; Ƞp2 = 0.119 (Fig. 3f, 3 g).
Recognition memory in NORT, measured as discrimination index, was neither affected by Tsc2+/-, DSE nor time (Fig. 3h).
Importantly, we verified that everolimus concentration in brain tissue at T2 was not significantly affected by either DSE or Tsc2+/- (Fig. 2d). It was also stable during the T2 behavioral testing period, and non-detectable after washing-out at T3 (Fig. 2d; Suppl. Fig. S3). Everolimus induced significant weight decrement, which was, however, quickly compensated after therapy discontinuation (Fig. 2e).
Tsc2 +/- and DSE induce epileptiform activity in adult rats which is not eliminated by everolimus
EEG in combination with video monitoring was recorded during T1 (session 1), T2 (session2), four weeks later (session 3) and during T3 (session 4). EEG records showed normal and symmetric background activity and no signs of any generalized tonic-clonic seizures. However, epileptiform activity patterns were apparent in the data (Fig. 4, Supplementary Tab. S1).
During all sessions, patterns with frequency of 15–20 Hz (mean duration 576 ± 11 ms) were identified in the EEG (Fig. 4d), corresponding to wet dog shake (WDS) behavior on video recordings. Statistical analysis by three-way ANOVA of sessions 1–4 showed a significant effect of Tsc2+/-, F(1, 11) = 6.850; p = 0.024; Ƞp2 = 0.384, with Tsc2+/- rats exhibiting more WDS. Simple effects showed that the difference between genotypes was significant in session 2 (p = 0.001) and session 3 (p = 0.019). The effect of DSE didn’t reach significance. There was an effect of time (F(3, 33) = 9.380; p < 0.001; Ƞp2 = 0.460) with a robust gradual decrease of WDS events across sessions and interaction between time and Tsc2+/- (F(3, 33) = 3.586; p = 0.024; Ƞp2 = 0.246). No specific effect of everolimus medication was apparent. Two-way ANOVA of individual sessions 1–4 mostly showed effects of both Tsc2+/- and DSE. Missing differences during session 4 could be due to a floor effect, as WDS decreased with time. We analyzed pooled data from sessions 1–4 to assess overall occurrence of WDS. Two-way ANOVA and additional Bonferroni post-hoc tests demonstrated individual and combined effects of both Tsc2+/- and DSE on increased WDS occurrence (Fig. 4a). To assess the temporal distribution of WDS within individual sessions, we pooled all sessions together and divided them into 5-min intervals. Three-way ANOVA showed a significant effect of Tsc2+/-, F(1, 14) = 13.007; p = 0.003; Ƞp2 = 0.482, a strong effect of time, F(7, 98) = 19.223; p < 0.001; Ƞp2 = 0.579, and interaction between time and Tsc2+/-, F(7, 98) = 2.202; p = 0.040; Ƞp2 = 0.136. The number of WDS was initially high, rapidly decreasing during the course of the session. Two-way ANOVA showed a disseminated pattern of differences caused by Tsc2+/- and DSE (Fig. 4b).
During periods of behavioral immobility we also noted long intervals (mean duration 3.5 ± 0.2 s) of rhythmic EEG activity (frequency 8–10 Hz) (Fig. 4e), which occurred predominantly in Tsc2+/- rats (27 events compared to 1 in wt).