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

doi:10.21203/rs.3.rs-2074843/v1

Download PDF

Research Article

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

https://doi.org/10.21203/rs.3.rs-2074843/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this older preprint version

Read the latest preprint version →

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 WT mice were similar between TscD^+/− and Tsc2^+/− mice, and the overlapping genes whose expression was altered in mutant mice compared with wild type (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.

double mutations

Pdlim2

rapamycin

Tsc1

Tsc2

tuberous sclerosis complex

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–5]. One such pathway is the mechanistic target of rapamycin complex 1 (mTORC1) signaling cascade [6–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–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–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–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–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–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–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.

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℃ under a 12 h/12 h light/dark cycle (lights on 8:00 AM to 8:00 PM) 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

All behavioral experiments were performed with male and female mice at 14–46 weeks of age between 10:00 AM and 4:00 PM.

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 minutes. Total distance travelled 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-genome gene expression 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 8x60K 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.

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 Fig. 1 test compared with WT mice (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 (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 of 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).

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 Fig. 2 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.

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 log₂ (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 log₂ (fold-change) in TscD^+/− mice than in Tsc1^+/− and Tsc2^+/− mice (Fig. 3A).

In Tsc1^+/-, Tsc2^+/-, and TscD^+/- mice compared with WT mice, 1,523 transcripts, 1,225 transcripts, and 4,954 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 Fig. 3 analysis. Of the 1,119 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” (Table S1). Furthermore, the 148 transcripts that overlapped between all three DETs were enriched in the seven significant networks in the canonical Fig. 4 pathway modeling analysis in MetaCore (p < 0.05, zScore > 100; 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 (1,523 transcripts in Tsc1^+/- vs. WT, 1,225 transcripts in Tsc2^+/- vs. WT, and 4,954 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 Fig. 5
higher in mutant mice than in WT mice and reduced by rapamycin treatment (Fig. 4B).

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–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–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]. Tsc1^f/−;Tsc2^f/−;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 (Tsc2^GFAP1) mice exhibited a high frequency of seizures and shorter life span than Tsc1^GFAP1 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.

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, 46]. 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 [47]. 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 (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 (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 [48–50]. 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 [51–54]. 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 [55, 56]. Pdlim2 is reported to be associated with tumor suppression [57–59] and tumorigenesis [60, 61], but its function in the nervous system remains unknown.

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.

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.

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.

Ethics approval

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.

Availability of data and materials

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

Competing interests

The authors declare that they have no competing interests.

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.

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.

Acknowledgments

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.

Doernberg E, Hollander E. Neurodevelopmental Disorders (ASD and ADHD): DSM-5, ICD-10, and ICD-11. CNS Spectr. 2016;21:295-299.
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-384.
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-361.
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.
Manoli DS, State MW. Autism spectrum disorder genetics and the search for pathological mechanisms. Am J Psychiatry. 2021;178:30-38.
Santini E, Klann E. Reciprocal signaling between translational control pathways and synaptic proteins in autism spectrum disorders. Sci Signal. 2014;7:re10.
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-13842.
Sahin M, Sur M. Genes, circuits, and precision therapies for autism and related neurodevelopmental disorders. Science. 2015;350:aab3897.
Winden KD, Ebrahimi-Fakhari D, Sahin M. Abnormal mTOR activation in autism. Annu Rev Neurosci. 2018;41:1-23.
Datta AN, Hahn CD, Sahin M. Clinical presentation and diagnosis of tuberous sclerosis complex in infancy. J Child Neurol. 2008;23:268-273.
Napolioni V, Moavero R, Curatolo P. Recent advances in neurobiology of tuberous sclerosis complex. Brain Dev. 2009;31:104-113.
Buxbaum JD. Multiple rare variants in the etiology of autism spectrum disorders. Dialogues Clin Neurosci. 2009;11:35-43.
European Chromosome 16 Tuberous Sclerosis Consortium. Identification and characterization of the tuberous sclerosis gene on chromosome 16. Cell. 1993;75:1305-1315.
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-808.
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.
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-546.
Curatolo P, Moavero R, de Vries PJ. Neurological and neuropsychiatric aspects of tuberous sclerosis complex. Lancet Neurol. 2015;14:733-745.
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.
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.
Jozwiak J, Jozwiak S, Wlodarski P. Possible mechanisms of disease development in tuberous sclerosis. Lancet Oncol. 2008;9:73-79.
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.
Mizuguchi M, Ohsawa M, Kashii H, Sato A. Brain symptoms of tuberous sclerosis complex: pathogenesis and treatment. Int J Mol Sci. 2021;22:6677.
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-515.
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-1024.
Curatolo P, Moavero R, Roberto D, Graziola F. Genotype/phenotype correlations in tuberous sclerosis complex. Semin Pediatr Neurol. 2015;22:259-273.
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-1695.
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-2029.
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.
Gao X, Pan D. TSC1 and TSC2 tumor suppressors antagonize insulin signaling in cell growth. Genes Dev. 2001;15:1383-1392.
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-4059.
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-3831.
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-1211.
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-8767.
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-655.
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-848.
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.
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-772.
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.
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-327.
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.
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-410.
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-1143.
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-372.
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.
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-454.
Switon K, Kotulska K, Janusz-Kaminska A, Zmorzynska J, Jaworski J. Molecular neurobiology of mTOR. Neuroscience. 2017;341:112-153.
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-415.
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.
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-16163.
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-577.
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.
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-849.
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-719.
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.
Tanaka T, Grusby MJ, Kaisho T. PDLIM2-mediated termination of transcription factor NF-kB activation by intranuclear sequestration and degradation of the p65 subunit. Nat Immunol. 2007;8:584-591.
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.
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-7368.
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-1420.
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-613.
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-285.
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.

No competing interests reported.

Download PDF

Version 1

posted

You are reading this older preprint version

Read the latest preprint version →

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

Archived Versions:

Version 1

Abstract

Background

Figures

Background

Materials And Methods

Mice

Drug

Behavioral tests

Whole-genome gene expression analysis

Gene ontology and pathway enrichment analyses

Statistical analysis

Results

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

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

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

Discussion

Limitations

Conclusions

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Archived Versions:

Version 1

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

Archived Versions:

Version 1

Abstract

Background

Figures

Background

Materials And Methods

Mice

Drug

Behavioral tests

Whole-genome gene expression analysis

Gene ontology and pathway enrichment analyses

Statistical analysis

Results

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

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

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

Discussion

Limitations

Conclusions

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Archived Versions:

Version 1

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

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

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