Dysregulation of neuron differentiation in an autistic savant with exceptional memory
Molecular Brain volume 12, Article number: 91 (2019)
Autism spectrum disorder (ASD) is a heterogeneous group of complex neurodevelopmental disorders without a unique or definite underlying pathogenesis. Although savant syndrome is common in ASD, few models are available for studying the molecular and cellular mechanisms of this syndrome. In this study, we generated urinary induced pluripotent stem cells (UiPSCs) from a 13-year-old male autistic savant with exceptional memory. The UiPSC-derived neurons of the autistic savant exhibited upregulated expression levels of ASD genes/learning difficulty-related genes, namely PAX6, TBR1 and FOXP2, accompanied by hypertrophic neural somas, enlarged spines, reduced spine density, and an increased frequency of spontaneous excitatory postsynaptic currents. Although this study involved only a single patient and a single control because of the rarity of such cases, it provides the first autistic savant UiPSC model that elucidates the potential cellular mechanisms underlying the condition.
Autism spectrum disorder (ASD) is a constellation of early-onset neurodevelopmental disorders. It is characterized by social communication deficits and repetitive sensory-motor behaviors, and it is often accompanied by abnormalities in language development . According to recent surveys, the prevalence of ASD among 8-year-old children in the United Sates was approximately 16.8 per 1000 children; the estimated prevalence in China was lower than those in other countries [2, 3]. ASD is a highly heritable mental disorder. Approximately 80% of cases of ASD have no clear etiopathology or related model, although many genetic variants have been identified in syndromic and non-syndromic ASD . Many ASD risk gene variants have been linked to brain connectivity, the excitation-inhibition balance, and synaptic functioning [5,6,7]. According to reports on risk-genes, TBR1 and five other genes, namely CHD8, KYRK1A, GRIN2B, PTEN and TBL1XR1, may contribute to 1% of sporadic ASD cases . Recently, a group of 65 genes were identified with strongest statistical evidence, and this identification has expanded knowledge of the sets of highest-confidence risk genes associated with ASD . Among these genes, TBR1—a putative transcription factor (TF)—is highly expressed in glutamatergic early-born cortical neurons; it dictates the expression of other risk genes, controls cortical development, and is implicated in intellectual disability [10,11,12]. TBR1 expression requires PAX6 [13,14,15], and PAX6 is related to autistic behaviors and speech abnormalities [16, 17]. Following TBR1 dictation, FOXP2 is also related to some severe speech-language disorders and plays a role in cortical neurogenesis [18,19,20,21].
Savant syndrome is a condition in which prodigious talent can co-occur with developmental conditions . In some cases of ASD, special abilities are accompanied by deficits; moreover, according to parental reports and psychometric tests, a third of adults with ASD exhibit savant skills in different domains . Special isolated memory skills were the most frequently reported special abilities . ASD children with special abilities exhibit more autistic traits, and multiple talent genes also influence the differences across individuals with ASD . However, few models are available for studying the molecular and cellular pathogenesis of autistic savants. In the study of neurodevelopmental disease, the induced pluripotent stem cell (iPSC) approach has been particularly useful [26,27,28,29]. Human iPSC models can enable the analysis of neuronal phenotypes and the investigation of cellular mechanisms after the derivation of autistic savants’ somatic cells into neurons.
In this study, we generated a urinary iPSC (UiPSC) model of a 13-year-old autistic boy with a photographic memory and speech-language deficit. This idiopathic savant exhibited repetitive behaviors and impaired social communication. We discovered that compared with control neurons, upregulated transcription of ASD risk genes co-occurred with dysregulated cellular cortical development and synaptogenesis in the UiPSC-derived neurons of the autistic savant on day 42 after neural progenitor cells (NPCs) differentiation. Our study is the first to provide a UiPSC model of an autistic savant with a photographic memory.
Generation of UiPSC-derived neurons of the autistic savant
Exfoliated renal epithelial cells were isolated from the urine of the autistic savant and an unrelated healthy control; the cells were then cultured for expansion (Fig. 1a and b and Additional file 3: Figure S1a). Approximately 3 weeks after isolation, urinary cells were replicated in sufficient quantity for the subsequent infection and exome sequencing studies (Additional file 1: Table S1). Attached urinary cells were infected using a Sendai-virus delivery system carrying the human OCT4, SOX2, KLF4, and c-MYC TFs; the urinary cells were then reprogramed into human UiPSCs (Fig. 1c and Additional file 3: Figure S1b). The advantages of this procedure have been described elsewhere [30,31,32]. We obtained three clone lines of UiPSCs from the autistic savant and two lines from the healthy control; all lines positively expressed pluripotent markers such as OCT4, NANOG, SOX2, SSEA4 and TRA-1-60 (Additional file 3: Figure S1c). All five UiPSC lines used in this study maintained a normal karyotype (Additional file 3: Figure S1d). The pluripotency of the UiPSC clones in vivo was confirmed by the generation of teratomas in severe combined immunodeficiency (SCID) mice. Furthermore, all UiPSC clones could generate teratomas that contained all three embryonic germ layers: endoderm, mesoderm and ectoderm (Additional file 3: Figure S1e). The results of teratomas generation proved that the five UiPSC clones could differentiate in vivo. Otherwise, positive alkaline phosphatase (AP) activities were also observed, confirming that all UiPSCs could maintain an undifferentiated proliferative state in vitro (Additional file 3: Figure S1f). These results indicated that the reprogramming of UiPSCs was successful.
The procedure for forebrain-specific neuronal differentiation is outlined in Fig. 1d. We initiated neuron differentiation by forming embryoid bodies (EBs). After approximately 3 weeks of differentiation, neural tube-like rosettes were manually collected, dissociated, and re-plated. All vital phases of NPC formation are illustrated in Additional file 3: Figure S1g. The UiPSC-derived NPCs were confirmed to positively express early neural precursor markers, such as nestin, SOX2 and musashi1 (Mui1) (Additional file 3: Figure S1g). After further neuron differentiation, the NPCs were differentiated into MAP2+ neurons with abundant expression of VGLUT1 (Additional file 4: Figure S2a), this observation is consistent with a previous description . Six-week-differentiation constitutes a sufficient maturation process for detecting spontaneous synaptic activity in human iPSC-derived forebrain neurons . In this study, we mainly focused on analyses of UiPSC-derived neurons on day 42 after NPCs differentiation (Fig. 1e).
Upregulated expression of PAX6, TBR1, and FOXP2 in UiPSC-derived neurons of the autistic savant
On day 42, we performed RNA sequencing (RNA-seq) analyses in UiPSC-derived neurons, without astrocyte co-culture, from all five clones. We detected 794 differentially expressed genes (DEGs), with 408 upregulations and 486 downregulations in the group of neurons from the autistic savant (Fig. 2a and b). Moreover, we observed an enrichment of upregulated pathways associated with mental disorders, such as the glutamatergic synapse pathway, Wnt signaling pathway, calcium signaling pathway and pathway of long-term potentiation, in the UiPSC-derived neurons of the autistic savant (Fig. 2c). An analysis of upregulated gene ontology (GO) enrichment revealed an abundance of genes associated with organ development, particularly neural system development, such as regulation of nervous system development, forebrain development, neurogenesis and negative regulation of nervous system development (Fig. 2d).
Among the aforementioned DEGs, we observed 55 upregulated TFs and 43 downregulated TFs in UiPSC-derived neurons of the autistic savant (Additional file 4: Figure S2b). High-confidence ASD risk genes, such as PAX6, TBR1, and FOXP2 [7, 10, 19, 34, 35], were significantly upregulated in the UiPSC-derived neurons of the autistic savant (Fig. 2a and Additional file 4: Figure S2b). Both PAX6 and TBR1 constituted the vital TFs of the upregulated protein-protein interaction (PPI) network in the group of neurons from the autistic savant (Fig. 2e). Cortical connections are shaped by sequential maps of regional identity, propagated by the PAX6-TBR1 TF cascade [11, 15]. Moreover, TBR1 was reported to be a putative TF that was highly expressed in glutamatergic early-born cortical neurons  and was demonstrated to be a highest-confidence risk gene associated with ASD . In the present study, Western blot analysis results revealed that TBR1 in the UiPSC-derived neurons of the autistic savant was upregulated by approximately six times compared with that in the neurons of the control (Fig. 2f). This finding was unexpected because most studies have reported many impaired conditions in TBR1 in cases of ASD [7, 12, 18, 19, 34]. We also determined that PAX6 and FOXP2 in the UiPSC-derived neurons of the autistic savant were upregulated by approximately two and three times compared with those in the neurons from the control (Fig. 2f).
Although previous studies have revealed that NR2B (GRIN2B) expression was directly regulated by TBR1 [19, 36], we observed no significant differences in mRNA and protein expression levels between the groups (Additional file 2: Table S2 and Fig. 2f). Previous studies on iPSC-derived neurons have examined neuronal subtype differentiation by conducting quantitative analyses of specific neural markers for neurodevelopmental diseases [27, 29]. According to our RNA-seq analyses, markers II-V of the neocortex (BCL11B, POU3F2, SATB2 and CUX1) and markers of glutamatergic neurons/GABAergic neurons (SLC17A7 or CAMK2A/GAD1 or SLC32A1) did not differ between the two groups (Additional file 2: Table S2 and Additional file 6: Table S4). Our Western blot analysis results also showed that BRN2 (POU3F2), CUX1, GAD67 (GAD1), and VGAT (SLC32A1) in the UiPSC-derived neurons of the autistic savant were not different from those in controls (Fig. 2g and h). However, GAD65 was downregulated in the UiPSC-derived neurons of the autistic savant (Fig. 2h). A previous study demonstrated that TBR1 regulated the number of synapses in heterozygous mutant mice .. In the present study, we detected synapsin 1 (SYN1; a presynaptic marker) and PSD95 (a postsynaptic marker) expression in the UiPSC-derived neurons. Both SYN1 and PSD95 observed in the UiPSC-derived neurons of the autistic savant were not different from those observed in the neurons of the healthy control (Fig. 2i). This finding is consistent with the RNA-seq analysis results (Additional file 2: Table S2).
Hypertrophic neural somas with enlarged spines and reduced spine density in UiPSC-derived neurons of the autistic savant
Previous studies have reported dysregulations of cellular phenotypes and the number of synapses in iPSC-derived neurons of idiopathic autistic individuals [28, 37, 38]. To investigate whether the UiPSC-derived neurons of the autistic savant had any morphological alterations when compared with the control neurons, we specifically labeled the membranes of the UiPSC-derived neurons with Dil by using a gene gun [39, 40]. The UiPSC-derived neurons, co-cultured on confluent astrocyte layers for 42 days after NPCs differentiation, were labeled with Dil and analyzed (Fig. 3a and b). Morphological analyses of the UiPSC-derived neurons of the autistic savant revealed that the dendritic spine density in these neurons was significantly reduced when compared with that in the control neurons (Fig. 3c). However, the average spine area and average spine length of the UiPSC-derived neurons of the autistic savant were greater than those of the control neurons (Fig. 3d and f). The average spine area of the UiPSC-derived neurons of the autistic savant was approximately three times larger than that of the control neurons (Fig. 3d). We measured the soma sizes of UiPSC-neurons without astrocyte co-culture on day 42. We conducted such measurements by analyzing bright-field images, and the analysis results showed that the average soma size of the UiPSC-derived neurons of the autistic savant was significantly larger than that of the control neurons (Fig. 3e). Overall, these data indicate the presence of hypertrophic neuronal somas with enlarged spines and reduced spine density in the UiPSC-derived neurons of the autistic savant.
Increased frequency but not amplitude of spontaneous excitatory postsynaptic currents in UiPSC-derived neurons of the autistic savant
Previous studies on ASD have observed that the electrophysiological function of iPSC-derived neurons was dysregulated following cellular and synaptic morphological alternations [28, 37, 38]. In the present study, we conducted whole-cell recordings on UiPSC-derived neurons with astrocyte co-culture on day 42 after NPCs differentiation. The frequency, but not amplitude, of the spontaneous excitatory postsynaptic currents (sEPSCs) of the UiPSC-derived neurons of the autistic savant was significantly larger than that of the control neurons on day 42 (Fig. 4a and b). These electrophysiological results suggest a greater synaptic release probability in the excitatory synapses of the UiPSC-derived neurons of the autistic savant.
We studied an autistic savant with a photographic memory and speech-language deficit. Our study demonstrated that PAX6, TBR1, and FOXP2 (as vital TFs) were significantly upregulated in the UiPSC-derived neurons of the autistic savant when compared with the control neurons. In addition, our genomic data reveal that both the autistic savant and healthy control had amplifications in different loci of PAX6. No clear interpretations for the upregulation of PAX6 could be gleaned from the genomic data. This requires further investigation with regard to our models. Notably, in most previous reports on ASD and intellectual disability, PAX6, TBR1 and FOXP2 were either deficiently expressed or had lost functions [6, 8, 9, 16]. TBR1 mutations have been implicated in impairments in social interaction, vocalization, memory and cognition [11, 12, 19]. PAX6 and FOXP2 have also been observed to be involved in the regulation of higher cognitive functions, including speech and language [16, 17, 41, 42]. In the present study, all of them were upregulated in the UiPSC-derived neurons of the autistic savant. According to studies on songbirds, FOXP2 is upregulated in a striatal region essential for song learning . On the basis of our results, we could hypothesize that PAX6, TBR,1and FOXP2 upregulation explains the superior memory along with the speech-language deficit of the idiopathic autistic savant.
A decrease or increase in dendritic spine density is commonly observed in animal models of ASD-related genes, such as Tbr1, Mecp2, Ube3a, Shank3, Tsc1, Fmr1, and Grin2b mice . In a mouse model, Tbr1, which is a potential master regulator in ASD, promotes synapse numbers through Wnt7b . Knockdown of FoxP2 reduces spine density in Area X of the zebra finch . In mice carrying two alleles of the human FOXP2 gene, dendrite length in the striatum increased . The downstream TBR1 target genes CDH10, GPC6, LMO7, and CNTN2, which encode membrane proteins that regulate cell adhesion and axonal outgrowth , were upregulated in the UiPSC-derived neurons of the autistic savant. The upregulation of these genes may explain the hypertrophic neural somas with enlarged spines and reduced spine density in the UiPSC-derived neurons of the autistic savant. To further explain the changes in soma size and dendritic spine density and size, studies on autistic savants regarding TBR1 or FOXP2 knockdown in UiPSC-derived neurons or overexpression of control neurons must be conducted in the future.
We did not observe a difference in the expression levels of GAD1 and SLC32A1, which encode GAD67 and VGAT, between the two neuron groups in the RNA-seq analyses and Western blot analysis results, respectively. Notably, GAD65 levels were obviously decreased in the UiPSC-derived neurons of the autistic savant, although no difference in GAD2 was observed in RNA-seq analyses. The analyses indicated the relationship between protein and mRNA levels under various scenarios, such as steady state, spatial and temporal variations of mRNAs, and local availability of resources for protein biosynthesis . The mechanisms underlying decreased GAD65 are still unclear and require further investigations. Previous studies have determined that GAD65 and GAD67 differed with regard to their cellular and subcellular locations, transcription times, and functional purposes [49,50,51]. GAD65 is localized in nerve terminals and synthesizes GABA for neurotransmission in functions such as synaptogenesis and protection from neural injury [50, 52]. Moreover, GAD65 transcription occurs slightly later on in the development period when synaptic inhibition is more prevalent . GAD65-synthesized GABA was reported to play a crucial role in the dynamic regulation of neural network excitability . During neurotransmission, GAD65 forms a complex with VGAT; this helps to package GABA into vesicles for release . Our results suggest a deficit of presynaptic inhibition in the synaptic release of the UiPSC-derived neurons of the autistic savant. This is a potential mechanism underlying the increased frequency of sEPSCs in the UiPSC-derived neurons of the autistic savant. The mechanism underlying the increased frequency of sEPSCs in the UiPSC-derived neurons of the autistic savant is not clear. The study showed that surviving TG (rTg(tauP301L)4510 tau mutant mice) cells had a significantly reduced total spine density but that the frequency of sEPSCs was increased, suggesting that during progressive tauopathy, cortical pyramidal cells compensate for the loss of afferent input with increased excitability . The compensation might be the reason for the increased frequency of sEPSCs. However, further investigation is required on this topic, including recording mEPSCs, mIPSCs and sIPSCs parameters in the future.
Our results illustrated a possible link between dysregulated ASD genes, an aberrant cellular and synaptic morphology, and impaired electrophysiological function in cases of ASD accompanied by savant syndrome. Our UiPSC model of the autistic savant is useful for further research on the mechanisms underlying such conditions. Although this study generated UiPSC models of an autistic savant, it has some limitations that warrant consideration. The primary limitation is the inclusion of only one autistic savant and one control; it has been very hard to find other autistic children with a photographic memory and speech-language deficit. In addition, considering the high level of clonal variation in iPSCs, our study findings apply to only this autistic savant with a photographic memory and speech-language deficit.
The autistic savant subject
A 13-year-old male autistic savant with exceptional memory was selected from a Chinese SuperBrain Talent Pool. He was diagnosed as having ASD and a speech-language deficit by four hospitals when he was 6 years old. According to reports from his parents, his exceptional memory was noticed in that same year. We tested his memory when he was 13 years old. He exhibited exceptional recall upon quickly scanning the content of books that he had never seen previously.
Collection and expansion of urinary cells
We established UiPSC models of a 12-year-old unrelated healthy control and a 13-year-old autistic savant with exceptional memory, both of whom were boys. The procedure and appearances are illustrated in Fig. 1b and Additional file 3: Figure S1a. We collected urine samples (approximately 100–400 mL) by using a sterile bottle containing a 5-mlL penicillin-streptomycin solution (Gibco, cat. no. 15070063). The samples were subsequently transferred to 50-mL centrifuge tubes and centrifuged at 1000 rpm for 10 min at room temperature. The pellets were collected into a single tube and washed with washing buffer, which was Dulbecco’s phosphate-buffered saline (DPBS, Gibco, cat. no. A1285601) supplemented with a penicillin-streptomycin solution. Urinary cells were resuspended in UC medium and transferred to a single well of a 12-well plate that had been coated with 0.1% gelatin (Millipore, cat. no. ES-006-B). The UC medium was a 1:1 mixture of complete Dulbecco’s modified Eagle’s medium (DMEM) containing DMEM/high glucose supplemented with 10% FBS and the renal epithelial cell growth medium (REGM; Lonza, cat. no. CC-3190). For each day during the first 4 days, we added 1 mL of UC medium. At approximately 96 h after plating, the urinary cells were removed all the medium, washed with DPBS, and added fresh UC medium. Half of the UC medium was replaced daily for approximately 10 days until the attached cells reached a density of 90%. Subsequently, the cells were passaged with TrypLE Select enzyme (Gibco, cat. no. 12563011) onto a 10-cm dish that was coated with gelatin for further expansion. This was considered passage 1 (P1). We continued to passage cells if required. The generation of UiPSCs from urinary cells at passages 1–3 was of a higher efficiency .
The protocol and appearances are illustrated in Fig. 1c and Additional file 3: Figure S1b. Urinary cells were plated in a gelatin-coated 6-well plate and incubated for 24 h. The cells were then infected using the CytoTune-iPS 2.0 Sendai Reprogramming Kit (Invitrogen, cat. no. A16517). The time point at which the cells were infected was considered day 0. After 24 h, Sendai virus-infected cells were replaced with fresh UC medium and changed every other day. On day 7, the infected cells were reseeded on the dishes that had been cultured primary mouse embryonic fibroblasts (Millipore, cat. no. PMEF-CFX) in UC medium. After 24 h of incubation, the infected cells were transferred from UC medium to UiPSC induction medium (UI medium). The UI medium comprised DMEM/F12, 20% KnockOut serum replacement (KSR; Gibco, cat. no. 10828028), 1% GlutaMAX (Gibco, cat. no. 35050061), 1% non-essential amino acid solution (NEAA; Gibco, cat. no. 11140050), 1% insulin-transferrin-selenium-ethanolamine (ITS-X; Gibco, cat. no. 51500056), 100 μM β-mercaptoethanol (β-Me; Sigma, cat. no. M3148), 50 ng/mL β-FGF (Peprotech, cat. no. 100-18B), 50 mg/mL ascorbic acid (Sigma, cat. no. A8960), 0.25 mM sodium butyrate (Sigma, cat. no. B5887), and 100 μg/mL Primocin (InvivoGen, cat. no. ant-pm-1). The reseeded cells were incubated for approximately 14 days, with the medium being replaced daily. After approximately 5 days of reseeding, the UiPSC colonies could emerge (Additional file 3: Figure S1b). When the colonies were ready for transfer, they could be picked and subsequently transferred to Matrigel-coated (Corning, cat. no. 354277) 24-well plates in UI medium. This was considered P1. The concentration of β-FGF was gradually decreased in the UiPSC culture medium at every passage until the UI medium contained 10 ng/mL β-FGF (UiPSC medium) at P5.
UiPSCs were plated on the Matrigel-coated cell image dishes (Eppendorf, cat. no. 0030740.017) and incubated for 5 days. NPCs were plated on Matrigel-coated coverslips (Fisher Scientific, cat. no. 12–548-82) in the 12-well plates and incubated for approximately 7 days. The cells were fixed with 4% paraformaldehyde (PFA) for 15 min at room temperature, blocked with 0.25% Triton X-100 (Sigma, X100) and 10% goat serum in DPBS for 40 min at room temperature, and incubated with primary antibodies (Additional file 5: Table S3) at 4 °C overnight. They were subsequently incubated with secondary antibodies (Additional file 5: Table S3) for 1 h at room temperature and then mounted with antifade mountant with DAPI (Invitrogen, cat. no. P36935). Antibodies were prepared in 0.25% Triton X-100 and 10% goat serum in DPBS. Image were captured by a TCS SP8 STED 3X multiphoton confocal microscope (Leica).
On day 5 after passaging, UiPSCs were treated in a colchicine solution (BBI Life Sciences, cat. no. A600322–0100) with a final concentration of 0.12 μg/mL for 1 h at 37 °C. Subsequently, the cells were dissociated with 0.05 mM EDTA and Trypsin mixed at a 9:1 ratio; the dissociated cells were then washed with DPBS. The pellet was gently resuspended with a 75 mM KCl solution and maintained in a 37 °C water bath for 20 min. We added 8–10 drops of a previously warmed fixative comprising methanol and acetic acid at a 3:1 ratio, after which we fixed the cells for 5 min at room temperature. After centrifugation, the pellet was replaced with a fresh fixative and fixed for 30 min at room temperature. The cells were then washed twice with a fixative and resuspended with 1 mL of fixative. Cell suspensions were sent for analysis to the Shanghai Da An Inspection for Medical Laboratory.
Teratoma formation assays
UiPSCs were subcutaneously injected into the dorsal flank of SCID mice. Eight to 12 weeks after injection, the mice were monitored, and teratomas were dissected, fixed with 4% PFA in DPBS at 4 °C overnight, and gradually embedded in paraffin. Tissues were sectioned at 5 μm, followed by haematoxylin and eosin staining. Images were captured by a microscope (Nikon Eclipse 80i).
AP activity in UiPSC assessment
On day 5 after passaging, UiPSCs were treated with an AP Detection Kit (Millipore, cat. no. SCR004) to determine AP activity. The UiPSCs were fixed with 4% PFA for 1–2 min at room temperature and washed with rinse buffer (20 Mm Tris-HCl solution with 0.05% Twwen-20, pH 7.4). Subsequently, the cells were replaced with stain solution and incubated in the dark at room temperature for 15 min. After incubation, they were washed again. Images were captured by a microscope (OPTIKA XDS-2).
Neuron differentiation of UiPSCs
The procedure and appearances of UiPSC-derived neurons are illustrated in Fig. 1d and e and Additional file 3: Figure S1g. We differentiated the UiPSCs into forebrain-specific NPCs, as previously described . The UiPSCs were detached with 1 mg/mL collagenase solution (Gibco, cat. no. 17104019) for 1 h at 37 °C; they were resuspended in EB medium—which was DMEM/F12 medium with 20% KSR, 1% GlutaMAX, 1% NEAA, 2 μM A 83–01 (Sigma, cat. no. SML0788), 2 μM Dorsomorphin (Sigma, cat. no. P5499), 100 μM β-Me and 100 μg/mL Primocin in not-treated polystyrene 6-well plates (Corning, cat. no. 3736)—for 6 days, with the medium being replaced daily. EBs were then replaced with NPC induction medium (NPC medium); this medium comprised DMEM/F12 and Neurobasal Medium (Gibco, cat. no. 21103049) mixed at 1:1 ratio, in addition to 1% NEAA, 1% N2 supplement (Gibco, cat. no. 17502048), 1% B27 supplement (Gibco, cat. no. 12587010), 2 μM cyclopamine (Sigma, cat. no. C4116), 2 μg/mL heparin (Sigma, cat.no. H3149), and 100 μg/mL Primocin. Subsequently, the EBs were transferred to Matrigel-coated 6-well plates to form neural tube-like rosettes. Rosettes were incubated for 16 days, with the NPC medium being replaced every other day. The rosettes were then picked and transferred to ultra-low-attachment-surfaced plates (Corning, cat. no. 3471) to form neural spheres in NPC medium. Alternatively, the rosettes were picked, dissociated with Accutase (Gibco, cat. no. A1110501) into single cells, and transferred to Matrigel-coated plates and culture slides in NPC medium. If they subsequently differentiated into neurons, NPCs were replaced with neuron medium (NM); the NM comprised Neurobasal Medium, 1% GlutaMAX, 1% NEAA, 1% N2, 1% B27, 10 ng/mL BDNF, 10 ng/mL GDNF, and 100 μg/mL Primocin after attachment. This was considered to be day 0 for the UiPSC-derived neurons. UiPSC-derived neurons were replaced weekly with half the culture medium. This occurred during continuous culturing and stopped when they were tested.
On day 42, UiPSC-derived neurons without astrocyte co-culture were lysed with TRIZOL (Invitrogen) for gene expression analyses. The sample lysates of five lines were delivered at low temperature. RNA and library preparation, clustering, sequencing, and data analyses were performed by the Novogene Experimental Department. Sequencing libraries were generated using the NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB) according to the manufacturer’s protocol; index codes were added to attribute sequences to each sample. After cluster generation, the library preparations were sequenced on an Illumina platform, and 125 bp/150 bp paired-end reads were generated. A differential expression analysis between the two groups was performed using the DESeq2 package in R (1.16.1). Genes with an adjusted p value of < 0.05 (obtained by DESeq2) were considered to be differentially expressed. A corrected p value of 0.05 and absolute fold change of 2 were set as the threshold for significantly differential expression. The KEGG is a database resource for understanding high-level functions and utilities of a biological system. Data were analyzed using the clusterProfiler package in R to test the statistical enrichment of DEGs in KEGG pathways. A GO enrichment analysis was implemented with the same package. A protein-protein interaction (PPI) analysis of DEGs was based on the STRING database. Otherwise, we used the Cytoscape (3.7.0) to construct the networks of the DEGs.
Western blot analyses
On day 42, UiPSC-derived neurons were lysed with 250–300-μL ice-cold RIPA buffer (Sigma, cat. no. R0278) in each well of a 6-well plate; the RIPA buffer comprised a protease inhibitor and phosphatase inhibitor (Roche, cat. no. 4693132001 and cat. no. 4906837001). The cell lysis suspensions were incubated on ice for 30 min and then transferred to tubes. The UiPSC-derived neuron lysate was heated at 100 °C in a metal bath for 5 min; it was then immediately cooled in an ice bath. The sample matrixes were subsequently centrifuged at 12,000 g for 25 min at 4 °C. The protein lysate was aspirated, mixed with 5× loading buffer, dispensed into aliquots, and stored. The protein samples were electrophoresed and transferred to polyvinylidene fluoride (PVDF) membranes. The membranes were blocked with 5% BSA in Tris-buffered saline with Tween 20 (TBST; 0.05% Tween-20) for 2 h at room temperature, incubated overnight with primary antibodies (Additional file 5: Table S3) at 4 °C, incubated with secondary antibodies (Additional file 5: Table S3) for 1 h at room temperature, and washed with TBST (0.05% Tween-20). Antibodies were prepared in primary/secondary antibody dilution buffer (Beyotime Biotechnology, cat. no. P0256/P0258). Horseradish peroxidase signals were detected by a Tanon (5200) system. The derived images were analyzed using ImageJ (NIH).
Preparation of primary astrocytes
Primary astrocytes were prepared from the cortices of Institute of Cancer Research (ICR) mice on postnatal day 0 (P0). P0 cortices were dissected, cut into small pieces, and dissociated with papain (Thermo Scientific, cat. no. 88285) in Hank’s balanced salt solution (HBSS) at 37 °C for 30 min. Cells were washed with HBSS and filtered, and epithelial cells were removed by differential time attachment. The cells were resuspended in astrocyte medium (AM) that comprised DMEM, 10% FBS, 1% NEAA, 1% GlutaMAX and Primocin, and they were incubated at 37 °C for 2 h. Subsequently, cell suspension was collected and transferred to T25 flasks that have been coated with 0.1 mg/mL poly-D-lysine (PDL; Sigma, cat. no. P0296). Confluent astrocytes were treated with 20 μM cytosine β-D-arabinofuranoside (AraC; Sigma, cat. no. C6646) for 72 h to further eliminate other cell types. A total of 2 × 105 cells were seeded on PDL-coated tissue culture slides in 24-well plates with AM for further use in co-culture.
On days 41–43, patch clamp recordings from UiPSC-derived neurons with astrocyte co-culture were obtained at room temperature by using a Multiclamp 700B amplifier, Digidata 1550, and the corresponding acquisition software Clampex 10.3 from Molecular Devices (USA). Coverslips seeded with derived neurons were transferred into recording chambers and perfused continuously with an extracellular solution containing the following (in mM, all from Sigma-Aldrich): 128 NaCl, 5 KCl, 2 CaCl2, 1 MgCl2, 25 HEPES, and 10 glucose (added before use) with pH 7.4. All recordings were obtained with a glass patch electrode (3–6 MΩ tip resistance), which was filled with pipette solutions, in whole-cell recording mode and filtered at 1 kHz. In addition, sEPSCs were recorded for 5 min at a holding potential of − 60 mV by using Cs-based pipette solution containing the following (in mM, all from Sigma-Aldrich): 135 CsCl, 5 EGTA, 10 HEPES, 4 Mg-ATP, and 0.5 Na-GTP (pH 7.3). The series resistance was not compensated for but was monitored. Recordings with a > 25% change in series resistance were discarded. Data were analyzed using Clampfit 10.3 software.
UiPSC-derived neurons labeled by gene gun
On day 42, UiPSC-derived neurons with astrocyte co-culture were washed with DPBS and rapidly fixed with 4% PFA for 1–2 min at room temperature. The UiPSC-derived neurons were removed from the fixative and specifically labeled with a Helios Gene Gun (Bio-Rad, 297BR). Bullets were prepared using Dil Stain (1,1′-Dioctadecyl-3,3,3′,3′-Tetramethylindocarbocyanine Perchlorate; Invitrogen, cat. no. D282). A slice with neurons was shot three times by one bullet. We added 100 μL DPBS to the neurons f to prevent them becoming dry, and we incubated the neurons for 3 h in the dark at room temperature. Slices were mounted with antifade mountant. The excitation and emission wavelengths of Dil (D282) were 549 and 565 nm, respectively. Images were captured by a TCS SP8 STED 3X multiphoton confocal microscope and were analyzed using Fiji software (ImageJ).
Morphological structure of UiPSC-derived neurons analyses
We analyzed the morphological structure of the UiPSC-derived neurons; our analysis focused on spine density, spine area, spine length, and soma size. Illumination is presented in Fig. 3b. Images of neurons were captured and analyzed. Specifically, dendritic lengths were measured using a plug-in (Simple Neurite Tracer). Dendritic diameters, spine lengths, number of spines, and total areas of threshold (dendrite areas and spine areas) were measured using Fiji software. The average spine density was defined as the number of spines per 10 μm of dendritic length. The average spine area was defined as the total spine area divided by the number of spines. The area of total spines was defined as the difference between the total area of fluorescence and dendritic area. At least 950-μm-long dendrites from five or more neurons (from at least three cultures) were analyzed for each UiPSC line. Soma size was analyzed through ImageJ. At least 247 somas from three cultures without astrocyte co-culture were statistically analyzed for each UiPSC line.
Quantification and statistical analysis
All experiments were replicated at least three times, and data from parallel cultures were acquired, except the wWestern blot of PAX6. GraphPad Prism (6.07) software was used for all statistical analyses. All individual data are presented as the mean ± standard error of the mean (SEM). Statistical significance was accepted at p < 0.05. We used Student’s unpaired t test to compare two groups. If five lines were compared, we used a one-way analysis of variance (ANOVA) with post hoc tests between groups corrected for multiple comparisons (Tukey). Statistical analyses and scale bars are described in each figure legend.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
One-way analysis of variance
Autism spectrum disorder
Differentially expressed genes
Dulbecco’s modified Eagle’s medium
Dulbecco’s phosphate-buffered saline
Hank’s balanced salt solution
Institute of Cancer Research
induced pluripotent stem cell
Knock Out serum replacement
Non-essential amino acid solution
Neural progenitor cells
Postnatal day 0
Renal epithelial cell growth medium
Severe combined immunodeficiency
Standard error of the mean
Spontaneous excitatory postsynaptic currents
Tris-buffered saline with Tween 20
rTg(tauP301l)4510 tau mutant mice
- UI medium:
UiPSC induction medium
Urinary induced pluripotent stem cells
Lord C, Elsabbagh M, Baird G, et al. Autism spectrum disorder. Lancet. 2018;392:508–20.
Wang F, Lu L, Wang SB, et al. The prevalence of autism spectrum disorders in China: a comprehensive meta-analysis. Int J Biol Sci. 2018;14:717–25.
Baio J, Wiggins L, Christensen DL, et al. Prevalence of autism Spectrum disorder among children aged 8 years - autism and developmental disabilities monitoring network, 11 sites, United States, 2014. MMWR Surveill Summ. 2018;67:1–23.
Mariani J, Coppola G, Zhang P, et al. FOXG1-dependent Dysregulation of GABA/glutamate neuron differentiation in autism Spectrum disorders. Cell. 2015;162:375–90.
Ecker C, Spooren W, Murphy DG. Translational approaches to the biology of autism: false dawn or a new era? Mol Psychiatry. 2013;18:435–42.
De Rubeis S, He X, Goldberg AP, et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature. 2014;515:209–15.
Fazel Darbandi S, Robinson Schwartz SE, Qi Q, et al. Neonatal Tbr1 dosage controls cortical layer 6 connectivity. Neuron. 2018;100(4):831.
O’Roak BJ, Vives L, Fu W, et al. Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science. 2012;338:1619–22.
Sanders SJ, He X, Willsey AJ, et al. Insights into autism Spectrum disorder genomic architecture and biology from 71 risk loci. Neuron. 2015;87:1215–33.
Hevner RF, Shi L, Justice N, et al. Tbr1 regulates differentiation of the preplate and layer 6. Neuron. 2001;29:353–66.
Vegas N, Cavallin M, Kleefstra T, et al. Mutations in TBR1 gene leads to cortical malformations and intellectual disability. Eur J Med Genet. 2018;61(12):759.
McDermott JH, Study DDD, Clayton-Smith J, et al. The TBR1-related autistic-spectrum-disorder phenotype and its clinical spectrum. Eur J Med Genet. 2018;61:253–6.
Holm PC, Mader MT, Haubst N, et al. Loss- and gain-of-function analyses reveal targets of Pax6 in the developing mouse telencephalon. Mol Cell Neurosci. 2007;34:99–119.
Quinn JC, Molinek M, Martynoga BS, et al. Pax6 controls cerebral cortical cell number by regulating exit from the cell cycle and specifies cortical cell identity by a cell autonomous mechanism. Dev Biol. 2007;302:50–65.
Elsen GE, Hodge RD, Bedogni F, et al. The protomap is propagated to cortical plate neurons through an Eomes-dependent intermediate map. Proc Natl Acad Sci U S A. 2013;110:4081–6.
Davis LK, Meyer KJ, Rudd DS, et al. Pax6 3′ deletion results in aniridia, autism and mental retardation. Hum Genet. 2008;123:371–8.
Kim KC, Lee DK, Go HS, et al. Pax6-dependent cortical glutamatergic neuronal differentiation regulates autism-like behavior in prenatally valproic acid-exposed rat offspring. Mol Neurobiol. 2014;49:512–28.
den Hoed J, Sollis E, Venselaar H, et al. Functional characterization of TBR1 variants in neurodevelopmental disorder. Sci Rep. 2018;8:14279.
Deriziotis P, O’Roak BJ, Graham SA, et al. De novo TBR1 mutations in sporadic autism disrupt protein functions. Nat Commun. 2014;5:4954.
Enard W, Przeworski M, Fisher SE, et al. Molecular evolution of FOXP2, a gene involved in speech and language. Nature. 2002;418:869–72.
Sollis E, Deriziotis P, Saitsu H, et al. Equivalent missense variant in the FOXP2 and FOXP1 transcription factors causes distinct neurodevelopmental disorders. Hum Mutat. 2017;38:1542–54.
Hughes JEA, Ward J, Gruffydd E, et al. Savant syndrome has a distinct psychological profile in autism. Mol Autism. 2018;9:53.
Howlin P, Goode S, Hutton J, et al. Savant skills in autism: psychometric approaches and parental reports. Philos Trans R Soc Lond Ser B Biol Sci. 2009;364:1359–67.
Meilleur AA, Jelenic P, Mottron L. Prevalence of clinically and empirically defined talents and strengths in autism. J Autism Dev Disord. 2015;45:1354–67.
Happe F. Why are savant skills and special talents associated with autism? World Psychiatry. 2018;17:280–1.
Takahashi K, Tanabe K, Ohnuki M, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131:861–72.
Wen Z, Nguyen HN, Guo Z, et al. Synaptic dysregulation in a human iPS cell model of mental disorders. Nature. 2014;515:414–8.
Deshpande A, Yadav S, Dao DQ, et al. Cellular phenotypes in human iPSC-derived neurons from a genetic model of autism Spectrum disorder. Cell Rep. 2017;21:2678–87.
Marchetto MC, Belinson H, Tian Y, et al. Altered proliferation and networks in neural cells derived from idiopathic autistic individuals. Mol Psychiatry. 2017;22:820–35.
Zhou T, Benda C, Dunzinger S, et al. Generation of human induced pluripotent stem cells from urine samples. Nat Protoc. 2012;7:2080–9.
Zhou T, Benda C, Duzinger S, et al. Generation of induced pluripotent stem cells from urine. J Am Soc Nephrol. 2011;22:1221–8.
Afzal MZ, Strande JL. Generation of induced pluripotent stem cells from muscular dystrophy patients: efficient integration-free reprogramming of urine derived cells. J Vis Exp. 2015;95:e52032.
Halliwell RF. Electrophysiological properties of neurons derived from human stem cells and iNeurons in vitro. Neurochem Int. 2017;106:37–47.
Notwell JH, Heavner WE, Darbandi SF, et al. TBR1 regulates autism risk genes in the developing neocortex. Genome Res. 2016;26:1013–22.
Hansen DV, Lui JH, Parker PR, et al. Neurogenic radial glia in the outer subventricular zone of human neocortex. Nature. 2010;464:554–61.
Chuang HC, Huang TN, Hsueh YP. Neuronal excitation upregulates Tbr1, a high-confidence risk gene of autism, mediating Grin2b expression in the adult brain. Front Cell Neurosci. 2014;8:280.
Bagni C, Zukin RS. A synaptic perspective of fragile X syndrome and autism Spectrum disorders. Neuron. 2019;101:1070–88.
Zaslavsky K, Zhang WB, McCready FP, et al. SHANK2 mutations associated with autism spectrum disorder cause hyperconnectivity of human neurons. Nat Neurosci. 2019;22:556–64.
Gan WB, Grutzendler J, Wong WT, et al. Multicolor “DiOlistic” labeling of the nervous system using lipophilic dye combinations. Neuron. 2000;27:219–25.
Grutzendler J, Tsai J, Gan WB. Rapid labeling of neuronal populations by ballistic delivery of fluorescent dyes. Methods. 2003;30:79–85.
Lai CS, Fisher SE, Hurst JA, et al. A forkhead-domain gene is mutated in a severe speech and language disorder. Nature. 2001;413:519–23.
Mukamel Z, Konopka G, Wexler E, et al. Regulation of MET by FOXP2, genes implicated in higher cognitive dysfunction and autism risk. J Neurosci. 2011;31:11437–42.
Scharff C, Haesler S. An evolutionary perspective on FoxP2: strictly for the birds? Curr Opin Neurobiol. 2005;15:694–703.
Varghese M, Keshav N, Jacot-Descombes S, et al. Autism spectrum disorder: neuropathology and animal models. Acta Neuropathol. 2017;134:537–66.
Schulz SB, Haesler S, Scharff C, et al. Knockdown of FoxP2 alters spine density in area X of the zebra finch. Genes Brain Behav. 2010;9:732–40.
Enard W, Gehre S, Hammerschmidt K, et al. A humanized version of Foxp2 affects corticao-basal ganglia circuits in mice. Cell. 2009;137:961–71.
Chuang HC, Huang TN, Hsueh YP. T-Brain-1—a potential master regulator in autism spectrum disorders. Autism Res. 2015;8:412–26.
Liu Y, Beyer A, Aebersold R. On the dependency of cellular protein levels on mRNA abundance. Cell. 2016;165:535–50.
Erlander MG, Tillakaratne NJ, Feldblum S, et al. Two genes encode distinct glutamate decarboxylases. Neuron. 1991;7:91–100.
Pinal CS, Tobin AJ. Uniqueness and redundancy in GABA production. Perspect Dev Neurobiol. 1998;5:109–18.
Soghomonian JJ, Martin DL. Two isoforms of glutamate decarboxylase: why? Trends Pharmacol Sci. 1998;19:500–5.
Kaufman DL, Houser CR, Tobin AJ. Two forms of the gamma-aminobutyric acid synthetic enzyme glutamate decarboxylase have distinct intraneuronal distributions and cofactor interactions. J Neurochem. 1991;56:720–3.
Kash SF, Johnson RS, Tecott LH, et al. Epilepsy in mice deficient in the 65-kDa isoform of glutamic acid decarboxylase. Proc Natl Acad Sci U S A. 1997;94:14060–5.
Jin H, Wu H, Osterhaus G, et al. Demonstration of functional coupling between gamma -aminobutyric acid (GABA) synthesis and vesicular GABA transport into synaptic vesicles. Proc Natl Acad Sci U S A. 2003;100:4293–8.
Crimins JL, Rocher AB, Peters A, et al. Homeostatic responses by surviving cortical pyramidal cells in neurodegenerative tauopathy. Acta Neuropathol. 2011;122:551–64.
The authors are grateful to all those who participated in this study.
This work was supported by the National Natural Science Foundation of China (NSFC) 91632103, Program of Shanghai Subject Chief Scientist 17XD1401700, Shanghai Education Commission Research and Innovation Program 2019-01-07-00-02-E00037, “111” Program of Higher Education Discipline Innovation and “Eastern Scholar” (Shanghai Municipal Education Commission).
Ethics approval and consent to participate
All study procedures followed institutional medical ethics approval and animal protocols approved by the Bio-X Institutes and the Institutional Animal Care and Use Committee of Shanghai Jiao Tong University. Informed consent was obtained from members of the autistic savant’s family and the parent of the boy who was involved as the control.
Consent for publication
Written informed consent to publish was obtained from the patient, the healthy donor, and their respective parents.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Exome sequencing of family members and control _ CNVs and SNPs. Related to Fig. 1.
Table S2. Secondary RNA sequencing _ differentially expressed genes _ significant. Related to Fig. 2.
Figure S1. Sample images of UiPSCs generation and NPCs differentiation. Related to Fig. 1.
About this article
Cite this article
Song, J., Yang, X., Zhou, Y. et al. Dysregulation of neuron differentiation in an autistic savant with exceptional memory. Mol Brain 12, 91 (2019). https://doi.org/10.1186/s13041-019-0507-7
- Autistic savant
- Human induced pluripotent stem cells