Skip to main content

Tsc2 mutation rather than Tsc1 mutation dominantly causes a social deficit in a mouse model of tuberous sclerosis complex

Abstract

Background

Tuberous sclerosis complex (TSC) is an autosomal dominant disorder that is associated with neurological symptoms, including autism spectrum disorder. Tuberous sclerosis complex is caused by pathogenic germline mutations of either the TSC1 or TSC2 gene, but somatic mutations were identified in both genes, and the combined effects of TSC1 and TSC2 mutations have been unknown.

Methods

The present study investigated social behaviors by the social interaction test and three-chambered sociability tests, effects of rapamycin treatment, and gene expression profiles with a gene expression microarray in Tsc1 and Tsc2 double heterozygous mutant (TscD+/−) mice.

Results

TscD+/− mice exhibited impairments in social behaviors, and the severity of impairments was similar to Tsc2+/− mice rather than Tsc1+/− mice. Impairments in social behaviors were rescued by rapamycin treatment in all mutant mice. Gene expression profiles in the brain were greatly altered in TscD+/− mice more than in Tsc1+/− and Tsc2+/− mice. The gene expression changes compared with wild type (WT) mice were similar between TscD+/− and Tsc2+/− mice, and the overlapping genes whose expression was altered in mutant mice compared with WT mice were enriched in the neoplasm- and inflammation-related canonical pathways. The “signal transducer and activator of transcription 3, interferon regulatory factor 1, interferon regulatory factor 4, interleukin-2R α chain, and interferon-γ” signaling pathway, which is initiated from signal transducer and activator of transcription 4 and PDZ and LIM domain protein 2, was associated with impairments in social behaviors in all mutant mice.

Limitations

It is unclear whether the signaling pathway also plays a critical role in autism spectrum disorders not caused by Tsc1 and Tsc2 mutations.

Conclusions

These findings suggest that TSC1 and TSC2 double mutations cause autistic behaviors similarly to TSC2 mutations, although significant changes in gene expression were attributable to the double mutations. These findings contribute to the knowledge of genotype–phenotype correlations in TSC and suggest that mutations in both the TSC1 and TSC2 genes act in concert to cause neurological symptoms, including autism spectrum disorder.

Background

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder that is characterized by two core behavioral symptoms: deficits in social communication and social interaction, and restricted, repetitive patterns of behaviors, interests, or activities [1]. Despite its high clinical and etiological heterogeneity, recent genetic studies have revealed that ASD appears to share neurobiological mechanisms and be attributable to the dysregulation of some convergent pathways [2,3,4,5]. One such pathway is the mechanistic target of rapamycin complex 1 (mTORC1) signaling cascade [6,7,8,9].

Tuberous sclerosis complex (TSC), a prototypical mTORopathy, is an autosomal dominant disorder. Up to 50% of patients with TSC have comorbid ASD. The prevalence of TSC in overall ASD was estimated to be 0.01 [10,11,12]. Thus, TSC is one of the most common diseases that is responsible for syndromic ASD. Tuberous sclerosis complex is considered to be caused by pathogenic mutations in the TSC complex subunit 1 (TSC1) or TSC2 gene and is unique in the similarity of clinical phenotypes between pathogenic mutations in the two genes [13,14,15]. TSC1 and TSC2 proteins form a complex with complemental TBC domain family member 7, which suppresses a small guanosine triphosphate-binding protein, Rheb, and negatively regulates the mTORC1 signaling cascade [16]. The excessive activation of mTORC1 signaling that is caused by pathogenic mutations in the two genes is assumed to be the primary pathophysiological mechanism of TSC [17]. According to Knudson’s two-hit theory, TSC-associated tumors develop because of biallelic inactivation of the TSC1 or TSC2 gene [18]. Indeed, the loss of heterozygosity (LOH) of either gene is often detected in most TSC-associated tumors, and more than half of brain cortical tubers lack LOH [19,20,21], suggesting that the heterozygosity of TSC1 or TSC2 is sufficient to cause neuronal abnormalities [22].

Genotype–phenotype correlation studies suggest that TSC2 pathogenic mutations rather than TSC1 pathogenic mutations are associated with severe phenotypic manifestations, including early-onset seizures, a high rate of complications that are associated with ASD, low cognitive function, and a high number of renal angiomyolipomas [15, 23,24,25,26]. Somatic TSC1 mutations were identified in lesions from patients with a pathogenic germline TSC2 mutation, suggesting that TSC1/TSC2 double heterozygous mutations may promote tumorigenesis [21, 27, 28]. Therefore, evaluations of manifestations of double heterozygous pathogenic mutations compared with mutations of either gene alone will provide insights into the pathophysiology of TSC. However, previous studies exclusively focused on pathogenic mutations of either gene alone, and few studies have described TSC1/TSC2 double mutations [29,30,31].

The homozygous loss of either Tsc1 or Tsc2 in mice causes embryonic lethality [32, 33]. Conventional heterozygous mice recapitulate some neuropsychiatric phenotypes of TSC in humans, including learning and memory deficits and social impairments without apparent brain lesions [34,35,36]. Additionally, transient treatment with an mTORC1 inhibitor successfully rescued these deficits.

The present study focused on the impact of Tsc1 and Tsc2 double heterozygous mutations in mouse models. We analyzed social behaviors and gene expression profiles in Tsc1 and Tsc2 double heterozygous mice and assessed the therapeutic effects of the mTOR inhibitor rapamycin on autistic-like behaviors.

Materials and methods

Mice

Tsc1 heterozygous knockout (Tsc1+/−) and Tsc2 heterozygous knockout (Tsc2+/−) mice were generated and maintained on a C57BL6/J background as previously described [32, 33]. A Tsc1+/− mouse and a Tsc2+/− mouse were crossed to produce Tsc1+/+/Tsc2+/+ (WT) mice, Tsc1+/−/Tsc2+/+ (Tsc1+/−) mice, Tsc1+/+/Tsc2+/− (Tsc2+/−) mice, and Tsc1+/−/Tsc2+/− (TscD+/−) mice. C57BL/6J mice were purchased from CLEA Japan (Tokyo, Japan). After weaning at 3 weeks of age, all mice were maintained separately by sex (two to five mice per cage) in ventilated racks at 22 °C ± 2 °C under a 12 h/12 h light/dark cycle (lights on 8:00 A.M. to 8:00 P.M.) with access to a standard laboratory diet and water ad libitum.

Drug

Rapamycin (LC Laboratories, Woburn, MA, USA) was dissolved to 0.5 mg ml−1 with 10% dimethylsulfoxide (vehicle) and administered intraperitoneally in a volume of 10 ml kg−1 once daily for 2 days. The behavioral tests and brain collection for gene expression analysis began 24 h after the second rapamycin administration.

Behavioral tests

The hole board test, self-grooming test, social interaction test, three-chambered sociability test, elevated plus maze test, and social transmission of food preference test were performed with male and female mice at 14.4–21.6, 44.4–46.9, 17.1–43.4, 15.6–43.4, 15.6–26.1, and 16.0–21.7 weeks of age between 10:00 A.M. and 4:00 P.M., respectively. The number of WT, Tsc1+/−, Tsc2+/−, and TscD+/− mice that were used in the behavioral tests were 13–19, 14–20, 15–20, and 14–22, respectively.

Hole board test Each mouse was placed in a hole board apparatus (500 mm × 500 mm × 400 mm; Muromachi Kikai, Tokyo, Japan) with four holes (38 mm diameter) and allowed to freely explore it for 30 min. Behavior was automatically analyzed for counts of head dipping using a video tracking system.

Self-grooming test Each mouse was individually placed in a new cage without bedding material. After an initial 10-min habituation period, the animal was video-recorded for 10 min. The number of episodes of self-grooming and rearing behavior was measured.

Social interaction test Each mouse was left alone in its home cage for 15 min. A novel mouse (C57BL/6J) of the same sex and size was then introduced to the cage. The behavior of the test mouse was video recorded for 10 min. The total duration of active interaction (anogenital sniffing, allogrooming, close following, and mounting) was measured.

Three-chambered sociability tests A three-chamber arena (500 mm × 500 mm × 400 mm) was used to evaluate social approach, preference for social novelty, and social preference. After a 10 min habituation period, an unfamiliar C57BL/6J mouse of the same sex (stranger 1) was introduced to the columnar wire cage (100 mm diameter, 100 mm height) in a side chamber, and a metal block (object) was placed in another cage on the other side chamber. The test mouse was allowed to explore the entire arena for 10 min (Session 1: social approach). The block was then replaced with another unfamiliar mouse (stranger 2), and the test mouse was examined for an additional 10 min (Session 2: social novelty). The test mouse was then presented with another unfamiliar mouse (stranger 3) and a cagemate of the test mouse for 10 min (Session 3: social preference). The amount of time spent around each cage was automatically measured using a video tracking system. The approach-avoidance score was calculated as the following: (time exploring stranger 3) − (time exploring cagemate).

Elevated plus maze test Each mouse was individually placed in the center of the elevated plus maze apparatus (Muromachi Kikai) and allowed to freely explore the apparatus for 10 min. Total distance traveled and time spent on the open arms were automatically measured using a DV-Track Video Tracking System (Muromachi Kikai).

Social transmission of food preference test Twenty-four hours before testing, the test mouse was placed in a novel cage with a specially designed food jar that was filled with powdered chow. Eighteen hours before testing, the feeding jars were removed, and the mice were deprived of food until testing. To begin testing, a “demonstrator” mouse was given a food jar that was filled with powdered chow that was flavored by mixing it with either 1% ground cinnamon or 2% powdered cocoa. After 1 h, the food jar was removed and weighed to ensure that at least 0.2 g of food had been consumed. The demonstrator mouse was then immediately placed in the cage of the “observer” test mouse, and interactions between the demonstrator mouse and observer mouse were allowed for 30 min. The demonstrator mouse was then removed, and two jars of powdered food (one of each flavor) were placed in the cage with the observer mouse for 1 h. At the end of the 1 h choice session, the food jars were removed and weighed. The demonstrator mice were C57BL/6J mice that were matched to the test observer mouse according to weight and sex.

Whole transcriptome analysis

Total RNA was isolated from the right brain hemisphere of 12-week-old male mice using Sepasol-RNA I Super G (nacalai tesque, Kyoto, Japan) and purified with RNeasy Mini (Qiagen, Tokyo, Japan). Only RNA with an RNA Integrity Number above 8.0, determined by an Agilent Bioanalyzer (Agilent Technologies, Tokyo, Japan), was used for gene expression analysis. Cy3-labeled cRNA was prepared using the Low Input Quick Amp Labeling Kit (Agilent Technologies) according to the manufacturer’s protocol, hybridized to the SurePrint G3 Mouse Gene Expression v2 8 × 60 K Microarray (Agilent Technologies), and scanned using a SureScan Microarray Scanner (Agilent Technologies).

Microarray data were measured using Feature Extraction software (Agilent Technologies) and analyzed using GeneSpring GX software (Agilent Technologies). Differentially expressed transcripts (DETs) were identified by two-way analysis of variance (ANOVA) with p < 0.05. Common transcripts between more than two DETs were selected instead of multiple testing correction to minimize false positives. Clustering analysis was performed using hierarchical clustering analysis with the default Euclidean distance metric and default Ward’s linkage rule.

Gene ontology and pathway enrichment analyses

Gene ontology (GO) and pathway enrichment analyses were performed using MetaCore version 6.31 (https://portal.genego.com/cgi/data_manager.cgi; accessed July 22, 2022). Map Folders and GO Process Networks are mainly based on biochemical and signaling cascades and cellular processes, respectively. The statistical significance value and false discovery rate (FDR) were used to rank GO. Canonical pathway modeling was used to identify pathways with statistical significance, zScore, and gScore by considering the number of objects of datasets and canonical pathways in the networks.

Statistical analysis

The statistical analyses of behavioral data were performed using SPSS Statistics 24 software (IBM Japan Ltd., Tokyo, Japan). Motor, sensory, and social behaviors between mutant and WT mice were statistically evaluated using one-way ANOVA followed by the Tukey Honestly Significant Difference (HSD) test. Unpaired t-tests were used to evaluate differences between vehicle and rapamycin treatment in the social interaction test. Paired t-tests were used to evaluate differences between familiar and stranger mice in the three-chambered sociability tests and between cued food and non-cued food in the social transmission of food preference test. Values of p < 0.05 were considered statistically significant.

Results

Autistic-like behaviors in Tsc1+/−, Tsc2+/−, and TscD+/− mice

The homozygous knockout of either the Tsc1 or Tsc2 gene causes embryonic lethality in mice, whereas both Tsc1+/− and Tsc2+/− mice are apparently healthy, do not die prematurely, and lack apparent cerebral lesions and spontaneous seizures [32, 33]. TscD+/− mice were born at the expected Mendelian frequency and exhibited no apparent spontaneous seizures, no early death, and no delay of physical development (data not shown). Adult Tsc1+/− and Tsc2+/− mice have been reported to exhibit normal motor and sensory behaviors [36]. Similar to Tsc1+/− and Tsc2+/− mice, TscD+/− mice did not display any significant differences in motor or sensory behaviors in the hole board test (Fig. 1A), elevated plus maze test, or self-grooming test compared with WT mice (Additional file 1: Fig. S1A–C). TscD+/− mice normally discriminated the scent of a food that was presented by stranger mice in the social transmission of food preference test, similar to Tsc1+/− and Tsc2+/− mice, indicating that these mutant mice can differentiate novel nonsocial stimuli that are transmitted through social interaction (Additional file 1: Fig. S1D). TscD+/− mice exhibited more frequent self-grooming behavior than WT mice, similar to Tsc1+/− and Tsc2+/− mice (p = 0.024, 0.031, and 0.042 for Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, Tukey HSD test; Fig. 1B). TscD+/− mice also exhibited a reduction in interaction time compared with WT mice, similar to Tsc1+/− and Tsc2+/− mice, in the social interaction test (p = 0.001, 0.004, and 0.004 for Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, Tukey HSD test; Fig. 1C). In the social approach session of the three-chambered sociability test, TscD+/− mice exhibited normal social preference, with a significant difference between the novel mouse and non-social object, similar to WT, Tsc1+/−, and Tsc2+/− mice (Fig. 1D). In the social novelty session, WT mice spent significantly more time exploring the novel mouse than the familiar mouse, whereas TscD+/− mice exhibited no significant preference for social novelty, similar to Tsc1+/− and Tsc2+/− mice (p = 0.005, 0.581, 0.099, and 0.845 for WT, Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, paired t-test; Fig. 1E). In the social preference session, TscD+/− mice spent a similar amount of time exploring both the novel mouse and the cagemate, similar to Tsc2+/− mice, whereas WT and Tsc1+/− mice preferred the novel mouse significantly more than the cagemate (p < 0.001, p < 0.001, p = 0.725, and p = 0.422 for WT, Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, paired t-test; Fig. 1F). The calculated approach-avoidance score revealed a severe reduction of social preference for a novel mouse in Tsc2+/− and TscD+/− mice rather than Tsc1+/− mice (p = 0.045, p < 0.001, p < 0.001 for Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, compared with WT mice; p = 0.015, Tsc1+/− mice compared with Tsc2+/− mice, Tukey HSD test; Fig. 1G).

Fig. 1
figure 1

Social behaviors in Tsc1+/−, Tsc2+/−, and TscD+/− mice. A–C Frequencies of head-dipping (A) and self-grooming (B) and time spent in active interaction behaviors (C) in Tsc1+/−, Tsc2+/−, and TscD+/− mice in the hole board test, self-grooming test, and social interaction test, respectively. D Time spent exploring (social approach) a stranger mouse (S1) and non-social object (O) in the three-chamber test. E Time spent exploring a familiar mouse (S1) and stranger mouse (S2) in the social-novelty preference test. F Time spent exploring a cagemate (C) and stranger mouse (S3) in the social preference test. G Approach avoidance score calculated as the following: (time exploring S3) – (time exploring C). *p < 0.05; **p < 0.01; ***p < 0.001, one-way ANOVA followed by Tukey HSD test (AC, G) and paired t-test (D–F). Each bar indicates the mean ± SEM

Rapamycin ameliorated impairments in social behaviors in TscD+/− mice similarly to Tsc1+/− and Tsc2+/− mice

We previously reported that 5–10 mg kg−1 rapamycin ameliorated impairments in social interaction in Tsc1+/− and Tsc2+/− mice in the social interaction test [36]. The time spent engaged in active interaction was increased by 5 mg kg−1 rapamycin in TscD+/− mice in the social interaction test similarly to Tsc1+/− and Tsc2+/− mice, with no effect in WT mice (p = 0.559, 0.012, 0.016, and 0.001 for WT, Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, unpaired t-test; Fig. 2A). Rapamycin also increased social novelty preference in Tsc1+/− mice (p = 0.600 and 0.007 for vehicle and rapamycin treatment, respectively, paired t-test; Fig. 2B), Tsc2+/− mice (p = 0.513 and 0.008 for vehicle and rapamycin treatment, respectively, paired t-test; Fig. 2B), and TscD+/− mice (p = 0.749 and 0.003 for vehicle and rapamycin treatment, respectively, paired t-test; Fig. 2B) and social preference in the three-chambered sociability tests in Tsc2+/− mice (p = 0.106 and 0.007 for vehicle and rapamycin treatment, respectively, paired t-test; Fig. 2C) and TscD+/− mice (p = 0.314 and 0.042 for vehicle and rapamycin treatment, respectively, paired t-test; Fig. 2C), which were similar to Tsc2+/− mice but not Tsc1+/− mice.

Fig. 2
figure 2

Effects of rapamycin on social behaviors in Tsc1+/−, Tsc2+/−, and TscD+/− mice. A–C Effects of rapamycin and vehicle on the time spent in active interaction in the social interaction test (A), time spent exploring a familiar mouse (S1) and stranger mouse (S2) in the social novelty session of the three-chambered sociability test (B), and time spent exploring a stranger mouse (S3) and cagemate (C) in the social preference session of the three-chambered sociability test (C). *p < 0.05; **p < 0.01; ***p < 0.001, unpaired t-test (A) and paired t-test (B, C). Each bar indicates the mean ± SEM. Veh vehicle, Rap rapamycin

Changes in gene expression in Tsc1+/−, Tsc2+/−, and TscD+/− mice

To evaluate similarities and differences in autistic-like behaviors in mutant mice, gene expression levels were profiled in the brain in mutant mice. After excluding control and low-intensity transcripts from a total of 56,745 transcripts that were included in the microarray, 45,313 transcripts were analyzed. The standard deviations of log2 (fold-change) were 0.156, 0.165, and 0.226 for Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, compared with WT mice, and the approximate curve showed a larger variation of log2 (fold-change) in TscD+/− mice than in Tsc1+/− and Tsc2+/− mice (Fig. 3A).

Fig. 3
figure 3

Gene expression profiles in the brain in Tsc1+/−, Tsc2+/−, and TscD+/− mice. A Histograms of the log2 (fold-change) (n = 45,313, bin width = 0.02) for each mutant mouse compared with WT (cyan: Tsc1+/− vs. WT, magenta: Tsc2+/− vs. WT, blue: TscD+/− vs. WT). Solid lines represent approximate normal distribution curves. SD standard deviation. B Venn diagram of three DETs between mutant and WT mice. Among overlapping transcripts (yellow area), transcripts that changed in different directions (▲: upregulated, : downregulated) were excluded from each entity for further analysis. C Heat map of valid 1119 transcript profiles in vehicle-treated samples for four genotypes (n = 4 for each genotype) that were in common with at least two DETs and met the criteria. Hierarchical clustering was based on the relative expression value after baseline correction to the median of vehicle-treated WT samples (red: upregulated, blue: downregulated). D Top 20 rank map folders enriched for differentially expressed 1119 transcripts in mutant mice. The blue and red folders indicate cancer/neoplasm and immunity/inflammation, respectively. FDR false discovery rate

In Tsc1+/−, Tsc2+/−, and TscD+/− mice compared with WT mice, 1523 transcripts, 1225 transcripts, and 4954 transcripts, respectively, were identified as DETs (two-way ANOVA, p < 0.05; Fig. 3B). To minimize false-positive changes in gene expression, only transcripts that were in common with at least two DETs (yellow area) were the focus of the subsequent analysis. Of the 1119 transcripts that were obtained after excluding transcripts that were regulated adversely between DETs, 289 and 830 were expressed upward and downward, respectively. Hierarchical clustering analysis showed that the expression patterns of these significant transcripts in TscD+/− mice were similar to Tsc2+/− mice but not Tsc1+/− mice (Fig. 3C).

Among the top 20 canonical processes that were identified in Map Folders, two major categories were included: cancer or neoplasms (e.g., lung cancer, prostatic neoplasm, hepatocellular carcinoma, and stomach neoplasm; nine processes) and diseases or pathways related to immunity or inflammation (e.g., systemic lupus erythematosus, allergic dermatitis, and asthma; eight processes; Fig. 3D). This result was consistent with the Tsc1 and Tsc2 gene functions of suppressing tumors and regulating immune and inflammatory responses via the mTOR pathway. Among GO process networks, the top three enrichment processes included “cellular component organization or biogenesis,” “regulation of cation channel activity,” and “nervous system development” (Additional file 3: Table S1). Furthermore, the 148 transcripts that overlapped between all three DETs were enriched in the seven significant networks in the canonical pathway modeling analysis in MetaCore (p < 0.05, zScore > 100; Additional file 4: Table S2).

Differentially expressed transcripts that were reversed by rapamycin treatment in the “STAT3, IRF1, IRF4, IL-2R alpha chain, IFN-γ” pathway initiated from PDLIM2

The effect of rapamycin on transcriptional changes was investigated, and we assessed whether the effects corresponded to the results of the behavioral tests in mutant mice. Among the DETs that were obtained to compare mutant and WT mice (1523 transcripts in Tsc1+/− vs. WT, 1225 transcripts in Tsc2+/− vs. WT, and 4954 transcripts in TscD+/− vs. WT), the expression of 180, 79, and 174 transcripts was significantly reversed by rapamycin treatment in Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively (two-way ANOVA, p < 0.05; Fig. 4A). The pathway enrichment analysis was performed with canonical pathway modeling in MetaCore to elucidate signal transduction pathways that are related to DETs between mutant and WT mice and are reversed by rapamycin. These DETs were enriched in seven networks, five networks, and six networks in Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively, and “signal transducer and activator of transcription 3 (STAT3), interferon regulatory factor 1 (IRF1), IRF4, interleukin-2R (IL-2R) α chain, and interferon-γ (IFN-γ)” was the common network in all mutant mice (Fig. 4A). Most DETs in the network were regulated by STAT3, which was downstream of IRF1, IRF4, IL-2R α chain, and IFN-γ, as well as STAT4 (Fig. 5). The network was initiated from PDZ and LIM domain protein 2 (PDLIM2), which was the only gene that overlapped among DETs in all mutant mice. The gene expression of Pdlim2 was higher in mutant mice than in WT mice and reduced by rapamycin treatment (Fig. 4B).

Fig. 4
figure 4

Rapamycin-induced changes in gene expression in the brain in Tsc1+/−, Tsc2+/−, and TscD+/− mice. A The cyan-colored circle indicates transcripts whose expression significantly changed in Tsc1+/− mice and was reversed by treatment with rapamycin ([Tsc1+/− vs. WT] ∩ [Tsc1+/− Veh vs. Rap]). Magenta- and blue-colored circles indicate gene expression changes in Tsc2+/− and TscD+/− mice, respectively. Canonical pathway modeling analysis for these three groups of transcripts identified seven, five, and six significant networks in Tsc1+/−, Tsc2+/−, and TscD+/− mice, respectively (p < 0.05, zScore > 100). B Box-and-whisker plot of normalized log2 signal intensity for the gene expression of Pdlim2 in vehicle- or rapamycin-treated mutant and WT mice. **p < 0.01; ***p < 0.001 (two-way ANOVA). Veh vehicle, Rap rapamycin

Fig. 5
figure 5

Signaling pathway associated with DETs in Tsc1+/−, Tsc2+/−, and TscD+/− mice with rapamycin treatment. Arrows indicate directional edges, with main actions shown in green. All pathways are commonly initiated from PDLIM2 to STAT3 via STAT4, IRF1, IFNγ, and IRF in Tsc1+/−, Tsc2+/−, and TscD+/− mice (symbol definitions can be found at https://portal.genego.com/legends/MetaCoreQuickReferenceGuide.pdf; accessed July 14, 2022)

Discussion

Autism spectrum disorder is a polygenic and genetically heterogeneous disorder for which many genes and genetic loci have been identified as causal factors in human genetic studies. The contribution of each gene or genetic locus to ASD is low (up to a few %), indicating that no genetic model covers all or most all ASD cases. The mTOR signaling pathway is a convergent pathway that is disturbed by genetic variations in the TSC1 and TSC2 (TSC), PTEN (Bannayan-Riley-Ruvalcaba syndrome, Lhermite-Duclos syndrome, and Cowden syndrome), FMRP (fragile X syndrome), and NF1 (neurofibromatosis type 1) genes. All of these disorders are associated with ASD, which together account for approximately 3–5% of overall ASD cases [12]. Thus, Tsc1/2 heterogeneous knockout mice are considered models for at least 3–5% of ASD cases.

Tuberous sclerosis complex manifests with such neurological symptoms as epilepsy, intellectual disability, and syndromic ASD, the symptoms of which do not differ from non-syndromic ASD [37]. Therefore, TSC is a disease model of the pathogenesis of ASD. Tuberous sclerosis complex is considered to be caused by pathogenic germline mutations (e.g., single nucleotide variations, deletions, insertions, and copy number variations) in either the TSC1 or TSC2 gene [38,39,40]. These mutations were reported to be mostly structural mutations that induce missense (7–33%) and frameshift (24–29%) mutations and affect splice sites (13–40%), whereas 18–28% are nonsense mutations. A total of 9–15% of TSC patients were found to have no identifiable mutation, and genotype–phenotype causal relationships have not yet been elucidated [39,40,41]. Furthermore, somatic TSC1 mutations were identified in lesions from patients with a pathogenic germline mutation of TSC2 [21, 27, 28], suggesting that TSC1 and TSC2 mutations may work in concert to elicit phenotypic differences among TSC patients. The present study was performed with TscD+/− mice as a model of TSC patients with both TSC1 and TSC2 mutations.

Tsc1+/− mice were previously reported to exhibit lower active interaction in the social interaction test [34, 36]. Tsc2+/− mice and Ekar (Tsc2+/−) rats were also reported to exhibit lower active interaction in the social interaction test [36, 42, 43] and impairments in social behavior in the social novelty session of the three-chambered test [42], but they did not exhibit abnormalities in the social approach session of the three-chambered test [42, 44]. The present results indicated that heterozygous Tsc2 mutation-induced social impairment was slightly more severe than the heterozygous Tsc1 mutation. TscD+/− mice exhibited impairments in social behaviors that were similar to Tsc2+/− mice and were associated with the mTOR pathway. Tsc1 and Tsc2 proteins form a heterodimeric complex that constitutively inactivates mTORC1 through inhibition of the small guanosine triphosphatase (GTPase) protein Rheb. Tsc2 possesses the GTPase-activating protein (GAP) domain for Rheb in the C terminus, and complex formation with Tsc1 stabilizes Tsc2 and enhances its GAP activity [16]. TscD+/− mice were expected to exhibit additive or synergistic phenotypes of both Tsc1+/− and Tsc2+/− mice. Indeed, the number of astrocytes in the hippocampus was reported to be higher in TscD+/− mice than in Tsc1+/− and Tsc2+/− mice [30]. Tsc1f/−; Tsc2f/−; hGFAP-Cre mice have a shorter lifespan than mice with mutations of any of these genes alone [31]. Additionally, glial-fibrillary acidic protein (GFAP)-positive cell-specific Tsc2 conditional knockout (Tsc2GFAP1) mice exhibited a high frequency of seizures and shorter life span than Tsc1GFAP1 conditional knockout mice [45]. Taken together with the present results, Tsc2 rather than Tsc1 appears to be dominant for causing neurological symptoms, including impairments in social behaviors that are caused by heterozygous Tsc1 and/or Tsc2 mutations, and impairments in social behaviors in TscD+/− mice would be partially mediated through different mechanisms from astrocytic abnormalities in an mTOR-dependent manner.

Unlike humans, which have marked individual differences, inbred mouse strains on homogeneous genetic backgrounds can be analyzed with much smaller sample sizes than humans. Indeed, a microarray analysis with 6–8 mice or 57 humans reported that standard deviations of the log (sample/reference) were approximately 0.3 and 0.8, respectively [46]. The Agilent microarray has superior detection compared with other microarray platforms [47]. Our transcriptome data that were based on Agilent microarrays showed that standard deviations of log (mutant/WT) in Tsc1+/−, Tsc2+/−, and TscD+/− mice were as low as 0.158, 0.165, and 0.226, respectively (Fig. 3A). Therefore, because the mice that were used in the present study were on a homogeneous genetic background and the standard deviation for each mutant group was low, the transcriptome analysis in which four mice were analyzed for each genotype and treatment group was considered to have the necessary power to detect differences between groups. Brain structures are left–right asymmetrical and have different functions that are controlled by each hemisphere [48], and these asymmetries of brain structures and functions have been reported to be broken in ASD [49,50,51]. Left hemisphere impairments and microstructural abnormalities were found in high-functioning ASD [52, 53], whereas the right hemisphere failed to respond to temporal novelty in ASD [54]. Social deficits in ASD have been reported to be associated with gamma-band electroencephalographic asymmetry in the right hemisphere [55], and social processing is right hemisphere-dominant in vertebrates [56, 57]. Unknown are the precise brain regions and cellular subtypes that are responsible for social deficits in ASD that are associated with TSC. Therefore, whole transcriptome analysis was performed using the right brain hemisphere rather than the left hemisphere and specific brain regions in the present study.

Over 1,000 genes with alterations of gene expressions in each mutant mice compared with WT mice were identified in the study. The magnitude of changes in gene expression in TscD+/− mice was higher than in Tsc1+/− and Tsc2+/− mice. The number of transcripts whose expression was significantly different from WT mice was also higher in TscD+/− mice than in Tsc1+/− and Tsc2+/− mice. Our results indicated that unlike social behaviors, gene expression in the brain in TscD+/− mice was affected by synergistic effects of heterozygous Tsc1 and Tsc2 mutations. In the present study, few genes in the mTOR1 signaling pathway had alterations of expression in addition to the Tsc1 and Tsc2 genes. Similarly, genes in the mTOR signaling pathway were not enriched in the gene expression analysis of cortical tubers from TSC patients, suggesting that transcription does not necessarily feedback from changes in translation, whereas mTOR1 signaling regulates translation via ribosomal protein S6 kinase, polypeptide 1 (S6K1), and eukaryotic translation initiation factor 4E binding protein 1 (4E-BP1) [21, 58]. Transgenic mice in which eukaryotic translation initiation factor 4E (Eif4e) gene expression was driven by the β-actin promoter exhibited autistic-like behaviors and a ~ 30% increase in cap-dependent translation, but the enhanced translation was limited to only a part of proteins [59]. Our results indicated that all mutant mice exhibited impairments in social behaviors despite large differences in the number of genes with alterations of gene expression between mutant and WT mice, suggesting that impairments in social behaviors that were caused by Tsc1 and Tsc2 mutations would be attributable to changes in specific molecular pathways and not a general increase in transcription or translation.

The hierarchical clustering analysis indicated that gene expression patterns in overlapping genes with alterations of expression between mutant and WT mice were similar in Tsc2+/− and TscD+/− mice, which corresponded to the severity of impairments in social behavior in mutant mice. Furthermore, genes whose expression was different between Tsc1+/− and Tsc2+/− mice were mostly enriched in the group of genes that were altered in Tsc2+/− and TscD+/− mice but not in Tsc1+/− mice compared with WT mice, also indicating the similarity of changes in gene expression between Tsc2+/− and TscD+/− mice (Additional file 2: Fig. S2A). The GO analysis of overlapping genes showed that they were enriched in cancer/neoplasm- and immunity/inflammation-related processes. These cancer/neoplasm- and immunity/inflammation-related processes showed the up- and downregulation of genes that were altered in Tsc2+/− and TscD+/− mice (Additional file 2: Fig. S2B), respectively, indicating that the severity of symptoms that are related to cancer/neoplasm and immunity/inflammation would be different between TSC patients with underlying TSC1 and TSC2 mutations.

STAT3 is phosphorylated by mTOR1 and activated in TSC2-deficient cells, which are inhibited by rapamycin [60,61,62]. The expression of inflammatory gene markers, including IL-1β, IL-6, and Toll-like receptors, and phosphorylated STAT3 increased in brain lesions in TSC patients [63,64,65,66]. These previous studies support our findings that the “STAT3, IRF1, IRF4, IL-2R α chain, IFN-γ” pathway network was associated with impairments in social behaviors in mutant mice and suggest that suppressing inflammatory responses in TSC patients may improve social deficits. This pathway network is initiated by Pdlim2, which is ubiquitously expressed throughout the body of mice but especially in the genitourinary, nervous, and visceral organ systems. Pdlim2 encodes an E3 ubiquitin ligase and ubiquitinates STAT3 and STAT4 to lead to their degradation [67, 68]. Pdlim2 is reported to be associated with tumor suppression [69,70,71] and tumorigenesis [72, 73], but its function in the nervous system remains unknown.

Limitations

The “STAT3, IRF1, IRF4, IL-2R α chain, IFN-γ” pathway network initiated from PDLIM2 and STAT4 was associated with impairments in social behaviors in Tsc1+/−, Tsc2+/−, and TscD+/− mice, but the regulatory mechanisms of the mTOR pathway and social behaviors by the pathway network were not elucidated in the present study. It is also unclear whether the pathway network also plays a critical role in impaired social behaviors in autism spectrum disorders not caused by Tsc1 and Tsc2 mutations.

Conclusions

In the present study, detailed analyses of social behaviors and transcriptomes in Tsc1+/−, Tsc2+/−, and TscD+/− mice revealed the genetic pathology of impairments in social behaviors. Mutations in the TSC2 gene are likely to have dominant effects on social behaviors rather than TSC1 gene if these mutations have similar potency. However, rapamycin has therapeutic effects on impairments in social behaviors caused by the mutations in both genes. It is also suggested that social behavioral abnormalities in patients with mutations in both the TSC1 and TSC2 genes are similar to those in patients with mutations in the TSC2 gene. The Pdlim2-mediated signaling pathway is suggested to be associated with altered brain function and impaired social behaviors due to mutations in the TSC1 and TSC2 genes. The present findings further our understanding of genotype–phenotype relationships and may contribute to the development of medications that can ameliorate social deficits in ASD.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

4E-BP1:

4E binding protein 1

ANOVA:

Analysis of variance

ASD:

Autism spectrum disorder

DETs:

Differentially expressed transcripts

FDR:

False discovery rate

GAP:

GTPase-activating protein

GFAP:

Glial-fibrillary acidic protein

GO:

Gene ontology

GTPase:

Guanosine triphosphatase

HSD:

Honestly significant difference

IFN-γ:

Interferon-γ

IL-2R:

Interleukin-2R

IRF1:

Interferon regulatory factor 1

LOH:

Loss of heterozygosity

mTORC1:

Mechanistic target of rapamycin complex 1

PDLIM2:

PDZ and LIM domain protein 2

S6K1:

S6 kinase, polypeptide 1

STAT3:

Signal transducer and activator of transcription 3

TSC:

Tuberous sclerosis complex

TSC1 :

TSC complex subunit 1

Tsc1 +/− :

Tsc1 Heterozygous knockout

Tsc2 +/− :

Tsc2 Heterozygous knockout

TscD +/− :

Tsc1 And Tsc2 double heterozygous mutant

WT:

Wild type

References

  1. Doernberg E, Hollander E. Neurodevelopmental disorders (ASD and ADHD): DSM-5, ICD-10, and ICD-11. CNS Spectr. 2016;21:295–9.

    Article  Google Scholar 

  2. Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, et al. Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011;474:380–4.

    Article  CAS  Google Scholar 

  3. de la Torre-Ubieta L, Won H, Stein JL, Geschwind DH. Advancing the understanding of autism disease mechanisms through genetics. Nat Med. 2016;22:345–61.

    Article  Google Scholar 

  4. Masini E, Loi E, Vega-Benedetti AF, Carta M, Doneddu G, Fadda R, et al. An overview of the main genetic, epigenetic and environmental factors involved in autism spectrum disorder focusing on synaptic activity. Int J Mol Sci. 2020;21:8290.

    Article  CAS  Google Scholar 

  5. Manoli DS, State MW. Autism spectrum disorder genetics and the search for pathological mechanisms. Am J Psychiatry. 2021;178:30–8.

    Article  Google Scholar 

  6. Santini E, Klann E. Reciprocal signaling between translational control pathways and synaptic proteins in autism spectrum disorders. Sci Signal. 2014;7:re10.

    Article  Google Scholar 

  7. Huber KM, Klann E, Costa-Mattioli M, Zukin RS. Dysregulation of mammalian target of rapamycin signaling in mouse models of autism. J Neurosci. 2015;35:13836–42.

    Article  CAS  Google Scholar 

  8. Sahin M, Sur M. Genes, circuits, and precision therapies for autism and related neurodevelopmental disorders. Science. 2015;350:aab3897.

    Article  Google Scholar 

  9. Winden KD, Ebrahimi-Fakhari D, Sahin M. Abnormal mTOR activation in autism. Annu Rev Neurosci. 2018;41:1–23.

    Article  CAS  Google Scholar 

  10. Datta AN, Hahn CD, Sahin M. Clinical presentation and diagnosis of tuberous sclerosis complex in infancy. J Child Neurol. 2008;23:268–73.

    Article  Google Scholar 

  11. Napolioni V, Moavero R, Curatolo P. Recent advances in neurobiology of tuberous sclerosis complex. Brain Dev. 2009;31:104–13.

    Article  Google Scholar 

  12. Buxbaum JD. Multiple rare variants in the etiology of autism spectrum disorders. Dialogues Clin Neurosci. 2009;11:35–43.

    Article  Google Scholar 

  13. European Chromosome 16 Tuberous Sclerosis Consortium. Identification and characterization of the tuberous sclerosis gene on chromosome 16. Cell. 1993;75:1305–15.

    Article  Google Scholar 

  14. van Slegtenhorst M, de Hoogt R, Hermans C, Nellist M, Janssen B, Verhoef S, et al. Identification of the tuberous sclerosis gene TSC1 on chromosome 9q34. Science. 1997;277:805–8.

    Article  Google Scholar 

  15. Au KS, Williams AT, Roach ES, Batchelor L, Sparagana SP, Delgado MR, et al. Genotype/phenotype correlation in 325 individuals referred for a diagnosis of tuberous sclerosis complex in the United States. Genet Med. 2007;9:88–100.

    Article  CAS  Google Scholar 

  16. Dibble CC, Elis W, Menon S, Qin W, Klekota J, Asara JM, et al. TBC1D7 is a third subunit of the TSC1-TSC2 complex upstream of mTORC1. Mol Cell. 2012;47:535–46.

    Article  CAS  Google Scholar 

  17. Curatolo P, Moavero R, de Vries PJ. Neurological and neuropsychiatric aspects of tuberous sclerosis complex. Lancet Neurol. 2015;14:733–45.

    Article  Google Scholar 

  18. Hino O, Kobayashi T. Mourning Dr. Alfred G. Knudson: the two-hit hypothesis, tumor suppressor genes, and the tuberous sclerosis complex. Cancer Sci. 2017;108:5–11.

    Article  CAS  Google Scholar 

  19. Niida Y, Stemmer-Rachamimov AO, Logrip M, Tapon D, Perez R, Kwiatkowski DJ, et al. Survey of somatic mutations in tuberous sclerosis complex (TSC) hamartomas suggests different genetic mechanisms for pathogenesis of TSC lesions. Am J Hum Genet. 2001;69:493–503.

    Article  CAS  Google Scholar 

  20. Jozwiak J, Jozwiak S, Wlodarski P. Possible mechanisms of disease development in tuberous sclerosis. Lancet Oncol. 2008;9:73–9.

    Article  CAS  Google Scholar 

  21. Martin KR, Zhou W, Bowman MJ, Shih J, Au KS, Dittenhafer-Reed KE, et al. The genomic landscape of tuberous sclerosis complex. Nat Commun. 2017;8:15816.

    Article  CAS  Google Scholar 

  22. Mizuguchi M, Ohsawa M, Kashii H, Sato A. Brain symptoms of tuberous sclerosis complex: pathogenesis and treatment. Int J Mol Sci. 2021;22:6677.

    Article  CAS  Google Scholar 

  23. van Eeghen AM, Black ME, Pulsifer MB, Kwiatkowski DJ, Thiele EA. Genotype and cognitive phenotype of patients with tuberous sclerosis complex. Eur J Hum Genet. 2012;20:510–5.

    Article  Google Scholar 

  24. Kothare SV, Singh K, Hochman T, Chalifoux JR, Staley BA, Weiner HL, et al. Genotype/phenotype in tuberous sclerosis complex: associations with clinical and radiologic manifestations. Epilepsia. 2014;55:1020–4.

    Article  Google Scholar 

  25. Curatolo P, Moavero R, Roberto D, Graziola F. Genotype/phenotype correlations in tuberous sclerosis complex. Semin Pediatr Neurol. 2015;22:259–73.

    Article  Google Scholar 

  26. Overwater IE, Swenker R, van der Ende EL, Hanemaayer KB, Hoogeveen-Westerveld M, van Eeghen AM, et al. Genotype and brain pathology phenotype in children with tuberous sclerosis complex. Eur J Hum Genet. 2016;24:1688–95.

    Article  CAS  Google Scholar 

  27. Tyburczy ME, Wang JA, Li S, Thangapazham R, Chekaluk Y, Moss J, et al. Sun exposure causes somatic second-hit mutations and angiofibroma development in tuberous sclerosis complex. Hum Mol Genet. 2014;23:2023–9.

    Article  CAS  Google Scholar 

  28. Avgeris S, Fostira F, Vagena A, Ninios Y, Delimitsou A, Vodicka R, et al. Mutational analysis of TSC1 and TSC2 genes in tuberous sclerosis complex patients from Greece. Sci Rep. 2017;7:16697.

    Article  Google Scholar 

  29. Gao X, Pan D. TSC1 and TSC2 tumor suppressors antagonize insulin signaling in cell growth. Genes Dev. 2001;15:1383–92.

    Article  CAS  Google Scholar 

  30. Uhlmann EJ, Apicelli AJ, Baldwin RL, Burke SP, Bajenaru ML, Onda H, et al. Heterozygosity for the tuberous sclerosis complex (TSC) gene products results in increased astrocyte numbers and decreased p27-Kip1 expression in TSC2+/− cells. Oncogene. 2002;21:4050–9.

    Article  CAS  Google Scholar 

  31. Mietzsch U, McKenna J 3rd, Reith RM, Way SW, Gambello MJ. Comparative analysis of Tsc1 and Tsc2 single and double radial glial cell mutants. J Comp Neurol. 2013;521:3817–31.

    Article  CAS  Google Scholar 

  32. Kobayashi T, Minowa O, Kuno J, Mitani H, Hino O, Noda T. Renal carcinogenesis, hepatic hemangiomatosis, and embryonic lethality caused by a germ-line Tsc2 mutation in mice. Cancer Res. 1999;59:1206–11.

    CAS  Google Scholar 

  33. Kobayashi T, Minowa O, Sugitani Y, Takai S, Mitani H, Kobayashi E, et al. A germ-line Tsc1 mutation causes tumor development and embryonic lethality that are similar, but not identical to, those caused by Tsc2 mutation in mice. Proc Natl Acad Sci U S A. 2001;98:8762–7.

    Article  CAS  Google Scholar 

  34. Goorden SM, van Woerden GM, van der Weerd L, Cheadle JP, Elgersma Y. Cognitive deficits in Tsc1+/− mice in the absence of cerebral lesions and seizures. Ann Neurol. 2007;62:648–55.

    Article  Google Scholar 

  35. Ehninger D, Han S, Shilyansky C, Zhou Y, Li W, Kwiatkowski DJ, et al. Reversal of learning deficits in a Tsc2+/− mouse model of tuberous sclerosis. Nat Med. 2008;14:843–8.

    Article  CAS  Google Scholar 

  36. Sato A, Kasai S, Kobayashi T, Takamatsu Y, Hino O, Ikeda K, et al. Rapamycin reverses impaired social interaction in mouse models of tuberous sclerosis complex. Nat Commun. 2012;3:1292.

    Article  Google Scholar 

  37. Jeste SS, Varcin KJ, Hellemann GS, Gulsrud AC, Bhatt R, Kasari C, et al. Symptom profiles of autism spectrum disorder in tuberous sclerosis complex. Neurology. 2016;87:766–72.

    Article  Google Scholar 

  38. Liu J, Zhao W, Ou X, Zhao Z, Hu C, Sun M, et al. Mutation spectrums of TSC1 and TSC2 in Chinese women with lymphangioleiomyomatosis (LAM). PLoS ONE. 2019;14: e0226400.

    Article  CAS  Google Scholar 

  39. Lin S, Zeng JB, Zhao GX, Yang ZZ, Huang HP, Lin MT, et al. Tuberous sclerosis complex in Chinese patients: phenotypic analysis and mutational screening of TSC1/TSC2 genes. Seizure. 2019;71:322–7.

    Article  Google Scholar 

  40. Tyburczy ME, Dies KA, Glass J, Camposano S, Chekaluk Y, Thorner AR, et al. Mosaic and intronic mutations in TSC1/TSC2 explain the majority of TSC patients with no mutation identified by conventional testing. PLoS Genet. 2015;11: e1005637.

    Article  Google Scholar 

  41. Peron A, Vignoli A, Briola F, Morenghi E, Tansini L, Alfano RM, et al. Deep phenotyping of patients with tuberous sclerosis complex and no mutation identified in TSC1 and TSC2. Eur J Med Genet. 2018;61:403–10.

    Article  Google Scholar 

  42. Tang G, Gudsnuk K, Kuo SH, Cotrina ML, Rosoklija G, Sosunov A, et al. Loss of mTOR-dependent macroautophagy causes autistic-like synaptic pruning deficits. Neuron. 2014;83:1131–43.

    Article  CAS  Google Scholar 

  43. Waltereit R, Japs B, Schneider M, de Vries PJ, Bartsch D. Epilepsy and Tsc2 haploinsufficiency lead to autistic-like social deficit behaviors in rats. Behav Genet. 2011;41:364–72.

    Article  Google Scholar 

  44. Ehninger D, Sano Y, de Vries PJ, Dies K, Franz D, Geschwind DH, et al. Gestational immune activation and Tsc2 haploinsufficiency cooperate to disrupt fetal survival and may perturb social behavior in adult mice. Mol Psychiatry. 2012;17:62–70.

    Article  CAS  Google Scholar 

  45. Zeng LH, Rensing NR, Zhang B, Gutmann DH, Gambello MJ, Wong M. Tsc2 gene inactivation causes a more severe epilepsy phenotype than Tsc1 inactivation in a mouse model of tuberous sclerosis complex. Hum Mol Genet. 2011;20:445–54.

    Article  CAS  Google Scholar 

  46. Wei C, Li J, Bumgarner RE. Sample size for detecting differentially expressed genes in microarray experiments. BMC Genom. 2004;5:87.

    Article  Google Scholar 

  47. van Iterson M, t Hoen PA, Pedotti P, Hooiveld GJ, den Dunnen JT, van Ommen GJ, Boer JM, Menezes RX. Relative power and sample size analysis on gene expression profiling data. BMC Genom. 2009;10:439.

    Article  Google Scholar 

  48. Güntürkün O, Ströckens F, Ocklenburg S. Brain lateralization: a comparative perspective. Physiol Rev. 2020;100:1019–63.

    Article  Google Scholar 

  49. Mundorf A, Peterburs J, Ocklenburg S. Asymmetry in the central nervous system: a clinical neuroscience perspective. Front Syst Neurosci. 2021;15: 733898.

    Article  Google Scholar 

  50. Kong XZ, Postema MC, Guadalupe T, de Kovel C, Boedhoe PSW, Hoogman M, et al. Mapping brain asymmetry in health and disease through the ENIGMA consortium. Hum Brain Mapp. 2022;43:167–81.

    Article  Google Scholar 

  51. Sha Z, van Rooij D, Anagnostou E, Arango C, Auzias G, Behrmann M, et al. Subtly altered topological asymmetry of brain structural covariance networks in autism spectrum disorder across 43 datasets from the ENIGMA consortium. Mol Psychiatry. 2022;27:2114–25.

    Article  CAS  Google Scholar 

  52. Perkins TJ, Stokes MA, McGillivray JA, Mussap AJ, Cox IA, Maller JJ, et al. Increased left hemisphere impairment in high-functioning autism: a tract based spatial statistics study. Psychiatry Res. 2014;224:119–23.

    Article  Google Scholar 

  53. Peterson D, Mahajan R, Crocetti D, Mejia A, Mostofsky S. Left-hemispheric microstructural abnormalities in children with high-functioning autism spectrum disorder. Autism Res. 2015;8:61–72.

    Article  Google Scholar 

  54. Orekhova EV, Stroganova TA, Prokofiev AO, Nygren G, Gillberg C, Elam M. The right hemisphere fails to respond to temporal novelty in autism: evidence from an ERP study. Clin Neurophysiol. 2009;120:520–9.

    Article  CAS  Google Scholar 

  55. Van Hecke AV, Stevens S, Carson AM, Karst JS, Dolan B, Schohl K, et al. Measuring the plasticity of social approach: a randomized controlled trial of the effects of the PEERS intervention on EEG asymmetry in adolescents with autism spectrum disorders. J Autism Dev Disord. 2015;45:316–35.

    Article  Google Scholar 

  56. Forrester GS, Todd BK. A comparative perspective on lateral biases and social behavior. Prog Brain Res. 2018;238:377–403.

    Article  Google Scholar 

  57. Karenina K, Giljov A. Mother and offspring lateralized social behavior across mammalian species. Prog Brain Res. 2018;238:115–41.

    Article  Google Scholar 

  58. Switon K, Kotulska K, Janusz-Kaminska A, Zmorzynska J, Jaworski J. Molecular neurobiology of mTOR. Neuroscience. 2017;341:112–53.

    Article  CAS  Google Scholar 

  59. Santini E, Huynh TN, MacAskill AF, Carter AG, Pierre P, Ruggero D, et al. Exaggerated translation causes synaptic and behavioural aberrations associated with autism. Nature. 2013;493:411–5.

    Article  CAS  Google Scholar 

  60. Yokogami K, Wakisaka S, Avruch J, Reeves SA. Serine phosphorylation and maximal activation of STAT3 during CNTF signaling is mediated by the rapamycin target mTOR. Curr Biol. 2000;10:47–50.

    Article  CAS  Google Scholar 

  61. Zhou J, Wulfkuhle J, Zhang H, Gu P, Yang Y, Deng J, et al. Activation of the PTEN/mTOR/STAT3 pathway in breast cancer stem-like cells is required for viability and maintenance. Proc Natl Acad Sci U S A. 2007;104:16158–63.

    Article  CAS  Google Scholar 

  62. Weichhart T, Costantino G, Poglitsch M, Rosner M, Zeyda M, Stuhlmeier KM, et al. The TSC-mTOR signaling pathway regulates the innate inflammatory response. Immunity. 2008;29:565–77.

    Article  CAS  Google Scholar 

  63. Boer K, Jansen F, Nellist M, Redeker S, van den Ouweland AM, Spliet WG, et al. Inflammatory processes in cortical tubers and subependymal giant cell tumors of tuberous sclerosis complex. Epilepsy Res. 2008;78:7–21.

    Article  CAS  Google Scholar 

  64. Shu HF, Zhang CQ, Yin Q, An N, Liu SY, Yang H. Expression of the interleukin 6 system in cortical lesions from patients with tuberous sclerosis complex and focal cortical dysplasia type IIb. J Neuropathol Exp Neurol. 2010;69:838–49.

    Article  CAS  Google Scholar 

  65. Boer K, Crino PB, Gorter JA, Nellist M, Jansen FE, Spliet WG, et al. Gene expression analysis of tuberous sclerosis complex cortical tubers reveals increased expression of adhesion and inflammatory factors. Brain Pathol. 2010;20:704–19.

    Article  CAS  Google Scholar 

  66. Prabowo AS, Anink JJ, Lammens M, Nellist M, van den Ouweland AM, Adle-Biassette H, et al. Fetal brain lesions in tuberous sclerosis complex: TORC1 activation and inflammation. Brain Pathol. 2013;23:45–59.

    Article  CAS  Google Scholar 

  67. Tanaka T, Grusby MJ, Kaisho T. PDLIM2-mediated termination of transcription factor NF-κB activation by intranuclear sequestration and degradation of the p65 subunit. Nat Immunol. 2007;8:584–91.

    Article  CAS  Google Scholar 

  68. Tanaka T, Yamamoto Y, Muromoto R, Ikeda O, Sekine Y, Grusby MJ, et al. PDLIM2 inhibits T helper 17 cell development and granulomatous inflammation through degradation of STAT3. Sci Signal. 2011;4:ra85.

    Article  Google Scholar 

  69. Sun F, Xiao Y, Qu Z. Oncovirus Kaposi sarcoma herpesvirus (KSHV) represses tumor suppressor PDLIM2 to persistently activate nuclear factor κB (NF-κB) and STAT3 transcription factors for tumorigenesis and tumor maintenance. J Biol Chem. 2015;290:7362–8.

    Article  CAS  Google Scholar 

  70. Zhao L, Yu C, Zhou S, Lau WB, Lau B, Luo Z, et al. Epigenetic repression of PDZ-LIM domain-containing protein 2 promotes ovarian cancer via NOS2-derived nitric oxide signaling. Oncotarget. 2016;7:1408–20.

    Article  Google Scholar 

  71. Wurster KD, Hummel F, Richter J, Giefing M, Hartmann S, Hansmann ML, et al. Inactivation of the putative ubiquitin-E3 ligase PDLIM2 in classical Hodgkin and anaplastic large cell lymphoma. Leukemia. 2017;31:602–13.

    Article  CAS  Google Scholar 

  72. Kang M, Lee KH, Lee HS, Park YH, Jeong CW, Ku JH, et al. PDLIM2 suppression efficiently reduces tumor growth and invasiveness of human castration-resistant prostate cancer-like cells. Prostate. 2016;76:273–85.

    Article  CAS  Google Scholar 

  73. Bassiri K, Ferluga S, Sharma V, Syed N, Adams CL, Lasonder E, et al. Global proteome and phospho-proteome analysis of Merlin-deficient meningioma and schwannoma identifies PDLIM2 as a novel therapeutic target. EBioMedicine. 2017;16:76–86.

    Article  Google Scholar 

Download references

Acknowledgements

The authors thank the members of the animal facility at the Tokyo Metropolitan Institute for Medical Science. The authors thank Michael Arends for proofreading the manuscript.

Funding

This research was supported by Grants-in-Aid for Scientific Research from MEXT/JSPS 24659490, 26670491, 26860836, 15H01303, 16K15565, 16K08924, JP22H04922 (AdAMS), and AMED JP16ek0109015.

Author information

Authors and Affiliations

Authors

Contributions

HK and SK were involved in the gene expression analysis, gene ontology and pathway analyses, data interpretation, and manuscript drafting. AS was involved in the behavioral analysis. YH was involved in the preparation of mice. YN was involved in the microarray procedure. TK and OH provided heterozygous knockout mice. MM and KI conceived the study. All authors contributed to the manuscript and approved the submitted version.

Corresponding author

Correspondence to Kazutaka Ikeda.

Ethics declarations

Ethics approval and consent to participate

All animal procedures were approved by the Institutional Animal Care and Use Committee of the Tokyo Metropolitan Institute of Medical Science (Approval No. for 5 years: 18028, 19026, 20-019, 21-071, and 22-012).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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

. Fig. S1. Anxiety and sensory behaviors in Tsc1+/−, Tsc2+/−, and TscD+/− mice.

Additional file 2

. Fig. S2. Differences in changes in gene expression between Tsc1+/− and Tsc2+/− mice.

Additional file 3

. Table S1. GO enrichment analysis of significant DETs in the brain in Tsc1+/−, Tsc2+/−, and TscD+/− mice.

Additional file 4

. Table S2. Canonical model pathways commonly enriched in the brain in Tsc1+/−, Tsc2+/−, and TscD+/− mice.

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

Kashii, H., Kasai, S., Sato, A. et al. Tsc2 mutation rather than Tsc1 mutation dominantly causes a social deficit in a mouse model of tuberous sclerosis complex. Hum Genomics 17, 4 (2023). https://doi.org/10.1186/s40246-023-00450-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40246-023-00450-2

Keywords