Skip to main content

TSC patient-derived isogenic neural progenitor cells reveal altered early neurodevelopmental phenotypes and rapamycin-induced MNK-eIF4E signaling



Tuberous sclerosis complex (TSC) is a neurodevelopmental disorder with frequent occurrence of epilepsy, autism spectrum disorder (ASD), intellectual disability (ID), and tumors in multiple organs. The aberrant activation of mTORC1 in TSC has led to treatment with mTORC1 inhibitor rapamycin as a lifelong therapy for tumors, but TSC-associated neurocognitive manifestations remain unaffected by rapamycin.


Here, we generated patient-specific, induced pluripotent stem cells (iPSCs) from a TSC patient with a heterozygous, germline, nonsense mutation in exon 15 of TSC1 and established an isogenic set of heterozygous (Het), null and corrected wildtype (Corr-WT) iPSCs using CRISPR/Cas9-mediated gene editing. We differentiated these iPSCs into neural progenitor cells (NPCs) and examined neurodevelopmental phenotypes, signaling and changes in gene expression by RNA-seq.


Differentiated NPCs revealed enlarged cell size in TSC1-Het and Null NPCs, consistent with mTORC1 activation. TSC1-Het and Null NPCs also revealed enhanced proliferation and altered neurite outgrowth in a genotype-dependent manner, which was not reversed by rapamycin. Transcriptome analyses of TSC1-NPCs revealed differentially expressed genes that display a genotype-dependent linear response, i.e., genes upregulated/downregulated in Het were further increased/decreased in Null. In particular, genes linked to ASD, epilepsy, and ID were significantly upregulated or downregulated warranting further investigation. In TSC1-Het and Null NPCs, we also observed basal activation of ERK1/2, which was further activated upon rapamycin treatment. Rapamycin also increased MNK1/2-eIF4E signaling in TSC1-deficient NPCs.


MEK-ERK and MNK-eIF4E pathways regulate protein translation, and our results suggest that aberrant translation distinct in TSC1/2-deficient NPCs could play a role in neurodevelopmental defects. Our data showing upregulation of these signaling pathways by rapamycin support a strategy to combine a MEK or a MNK inhibitor with rapamycin that may be superior for TSC-associated CNS defects. Importantly, our generation of isogenic sets of NPCs from TSC patients provides a valuable platform for translatome and large-scale drug screening studies. Overall, our studies further support the notion that early developmental events such as NPC proliferation and initial process formation, such as neurite number and length that occur prior to neuronal differentiation, represent primary events in neurogenesis critical to disease pathogenesis of neurodevelopmental disorders such as ASD.


Tuberous sclerosis complex (TSC) is an autosomal dominant monogenic disorder with severe neurological manifestations including epilepsy, autism spectrum disorder (ASD), intellectual disability (ID), and hamartomas (benign tumor-like formations) in many organs. TSC is caused by mutations in the TSC1 or TSC2 gene, encoding tumor suppressor proteins hamartin (TSC1) and tuberin (TSC2) [1,2,3]. The TSC proteins form a functional complex that acts as a central hub relaying signals from diverse cellular pathways to inhibit mammalian/mechanistic target of rapamycin complex 1 (mTORC1) activity, which regulates cell growth and proliferation [4, 5]. In neuronal translation, mTORC1 signaling is a regulator of long-lasting synaptic plasticity and memory as it integrates signals from neuronal surface receptors/channels via MEK/ERK- and PI3K/AKT-mediated phosphorylation and inactivation of the TSC1-TSC2 complex [4,5,6]. The aberrant activation of mTORC1 in TSC has led to treatment with rapamycin analogs (rapalogs) as a lifelong therapy [7,8,9,10], with discontinuation leading to a rebound in growth of the TSC-associated lesions. Moreover, rapalog treatment has no significant effect on neurocognitive functioning or behavior in children with TSC [11]. Therefore, there is a clear need to identify novel therapeutics for treating TSC that are superior or complementary to rapalogs in terms of long-term effectiveness and efficacy toward various non-tumor CNS manifestations of TSC.

Several mouse models of TSC have provided valuable clues for neurological symptoms, but have limitations in faithfully recapitulating the human phenotypes [12]. Further, the inability to establish expandable human cell lines derived from various TSC-associated lesions, along with genetically matched control cell lines has made it difficult to define the precise pathogenic mechanisms involved in TSC. Patient-specific induced pluripotent stem cells (iPSCs) derived from somatic cells, followed by differentiation into specific cell types, are rapidly evolving to be powerful for disease modeling to study pathophysiology and to identify treatments [13,14,15,16,17]. More importantly, the emergence of powerful genome editing techniques has made it possible to generate isogenic pairs of disease and control human iPSCs that differ only with respect to disease-causing gene mutations [18,19,20]. Recent studies that employed either human embryonic stem cell lines with heterozygous or homozygous loss of TSC2 or TSC patient iPSCs examined post-mitotic neurons and focused on later developmental processes such as dendrite outgrowth and synapse formation [21,22,23,24,25,26]. However, modeling of neurodevelopmental disorders with patient-derived iPSCs from ASD patients suggests that ASD risk genes can affect early phenotypes of neurogenesis such as the proliferation of neural progenitor cells (NPCs). Recent studies reveal that earlier developmental events, such as NPC proliferation, process outgrowth, and migration that occur prior to neuronal differentiation are also critical in disease pathogenesis of ASD and schizophrenia [27,28,29,30,31].

In this study, we have generated isogenic sets of iPSCs from a TSC patient harboring a germline TSC1 mutation and have focused our efforts to examine the changes caused by either heterozygous or homozygous loss of TSC1 in NPCs. Both TSC1-Het and Null NPCs show enlarged cell size and mTORC1 activation when compared with the CRISPR-corrected WT, which are reversed by rapamycin. Further, we show basal activation of MEK-ERK signaling in TSC1-Het and Null NPCs, with further activation of ERK upon rapamycin treatment. Rapamycin also activates MNK-eIF4E signaling pathway, a regulator of 5’ cap-dependent translation. More importantly, we have observed genotype-dependent differences in early neurodevelopmental events such as aberrant NPC proliferation and neurite outgrowth, which are not affected by rapamycin treatment. Transcriptome analyses using RNA sequencing (RNAseq) revealed differential expression of genes related to ASD, ID, and epilepsy that were also altered in a genotype-dependent manner. Taken together, our results authenticate that both heterozygous and homozygous loss of TSC1 influence phenotypes, signaling, and gene expression in NPCs compared to the genetically matched control WT cells, supporting that heterozygous loss of TSC1/2 may indeed play a role in some of the neurological manifestations of TSC.


Human iPSC line derivation

TSC1 skin fibroblast samples were collected through the TSC Clinic at Massachusetts General Hospital (Boston, USA). Genomic DNA from fibroblasts was extracted, and a germline mutation in TSC1 exon 15 was identified using standard PCR and Sanger sequencing methods as previously described [32]. The TSC1 heterozygous fibroblasts were cultured in high glucose DMEM (Gibco) supplemented with 15% of fetal bovine serum (Sigma) and streptomycin/penicillin (Cellgro). Cells were reprogrammed using a synthetic modified mRNA-based method to obtain iPSC lines as described [33, 34]. Briefly, cells were transfected by nucleofection (Amaxa Nucleofector I) with in vitro transcribed mRNAs encoding OCT4, SOX2, KLF4, cMYC, and LIN28 (Stemgent). After picking clones, iPSC colonies were cultured in feeder-free culture conditions on Geltrex (ThermoFisher)-coated plates in Essential 8 medium (Gibco). Colonies were passaged every 4–6 days and the media was changed daily. We also utilized an unrelated wildtype control iPSC line 8330 that was originally generated from fibroblasts (GM08330) obtained from the Coriell Institute for Medical Research) as previously reported [35]. Karyotyping was performed by WiCell Cytogenetics Lab (Madison, WI). Sanger sequencing was performed by Eton Bioscience (San Diego, CA).

CRISPR/Cas9 method for generating isogenic iPSC lines

To generate isogenic iPSC lines (Corrected-WT and Null), CRISPR/Cas9 genome editing was performed employing TSC1 heterozygous patient-derived iPSCs. For the TSC1-Null lines, a single guide RNA (sgRNA) sequence was designed to target TSC1 exon 7 (sgRNA seq: GAGATAGACTTCCGCCACG). For CRISPR-mediated corrections, a sgRNA was designed to specifically target the mutant TSC1 exon 15 allele (sgRNA seq: GGGAGACTGTCTCAGTAAA) to correct the germline microdeletion mutation. sgRNAs were cloned into the pSpCas9(BB)-2A-Puro (PX459) vector and plasmids were prepared using an endo-free midi prep (Qiagen). To achieve a gene knock-in to correct the mutation on exon 15, we used a 99 base pair wildtype donor sequence designed for homologous recombination of the allele bearing the mutation. iPSCs were pre-incubated with ROCK-inhibitor at 10 μM for 2 h prior to nucleofection with 1 μg of vector using the human stem cell nucleofector kit I (Lonza) with the Amaxa Nucleofector I (program B-16). Cells were then plated and selected on the following day with 0.5 μg/ml of puromycin for 2 days. Selected cells were grown for 10 days allowing colonies to form. Colonies were manually isolated in Geltrex-coated 48 well plates and expanded. After 3–4 days, genomic DNA was extracted from a portion of each colony followed by Sanger sequencing using targeted primers (see Additional file 3: Table S2). Clones showing successful gene editing were then further expanded and used for subsequent experiments.

Differentiation of iPSCs into NPCs

TSC1-iPSC lines were differentiated using the directed monolayer differentiation protocol [36]. Briefly, iPSCs expressing the pluripotency marker TRA-1-60 were sorted and enriched using the MACS Microbead cell sorting technology (Miltenyi Biotec) and plated in feeder-free conditions at a density of 2 − 2.5 × 104 cells per cm2. Cells were cultured in neural induction media (neurobasal media supplemented with 1× neural induction supplements (Gibco)) for 7–9 days, after which they expressed polysialylated-neural cell adhesion molecule (PSA-NCAM). The PSA-NCAM-positive (+) cells were first isolated by MACS sorting, and then double sorted to enrich for NPCs representing CD271-/CD133+ cells. NPCs were then cultured in neural expansion media (50% neurobasal media and 50% advanced DMEM/F12 (Gibco) supplemented with 1× neural induction supplements (Gibco)) up to 15 passages and periodically assessed for expression of NPC markers, NESTIN, and SOX2.

Neurite outgrowth assay

NPCs were seeded on Poly-D-Lysine (0.1 mg/ml, Sigma) and Fibronectin (5 μg/ml, Corning) coated wells at 6250 cells per cm2 in growth factor depleted Neural Expansion Medium (30% NEM) containing 49.7% neurobasal media, 49.7% advanced DMEM (Gibco), 1× penicillin/streptomycin and 0.3× neural induction supplements. Cells were grown for 48 h and fixed with 4% paraformaldehyde (PFA) for 20 min prior to immunostaining. Four independent field images with approximately 50 cells per field were analyzed. Processes that were at least two times the length of the cell body were considered as neurites. The average neurite number per cell and the average longest neurite length per cell were analyzed using HCA-Vision software V2.2.0 (CSIRO).


Cells were fixed with 4% paraformaldehyde for 20 minutes at room temperature. Non-specific labeling was blocked using 4% Normal Goat Serum in PBS with 0.1% of Triton-X-100 for 45 min at room temperature. Primary antibodies were diluted into 2% NGS-PBS-Triton-X 0.1% and incubated overnight at 4°C (see Additional file 2: Table S1). Secondary antibodies were diluted into 2% NGS-PBS-Triton-X 0.1% and incubated for 2 h in the dark at room temperature (see Additional file 2: Table S1). DAPI was used to stain nuclei (Invitrogen #D3571) at 5 μg/ml. Coverslips were mounted in ProLong Gold Antifade Mountant (Invitrogen #P36930) and images were captured using a Nikon Eclipse TE2000-U microscope and the NIS-Element BR 3.2 imaging software.

Immunoblot analyses

Cells were lysed in RIPA buffer as previously described [37, 38]. Protein lysates were resolved on 4–20% CriterionTM TGXTM gel (BioRad), transferred to nitrocellulose (Biorad) and then incubated with primary antibodies (see Additional file 2: Table S1). All immunoblotting data shown is a representative of 3 biological replicates.

Quantitative RT-PCR

Total RNA from iPSCs was isolated by lysis in TRIzol reagent (Ambion/Life Technologies; Grand Island, NY) according to the manufacturer’s instructions. Following lysis, RNA was rinsed in chloroform, and the aqueous layer was applied to Qiagen RNeasy kit (Qiagen) columns followed by purification according to manufacturer’s instructions. For cDNA synthesis, the Superscript VILO cDNA synthesis kit (Life Technologies) was used according to the manufacturer’s instructions, and quantitative RT-PCR (q-RT-PCR) was carried out using TaqMan according to manufacturer’s instructions.

RNA-seq library preparation and sequencing

Total RNA was isolated from NPC lines using TRIzol reagent. Briefly, pelleted cells were resuspended in TRIzol reagent and then extracted with chloroform, followed by isopropanol precipitation of RNA from the aqueous phase and three 70% ethanol washes. RNA pellets were solubilized in 30-50 μl of RNase-free water (Ambion, AM9937). RNA quality was assessed using the Agilent Bioanalyzer TapeStation 2200 (Agilent Technologies, Santa Clara CA). In total, 12 RNA-seq libraries were prepared in triplicate for each of the four NPC samples harboring WT, Het TSC1, and two clones of null TSC1 using the Illumina TruSeq Stranded mRNA Sample Prep Kit. Each library in this study included 1 μl of a 1:10 dilution of ERCC RNA Control Spike-Ins (Ambion) that were added from one of two mixes, each containing the same 92 synthetic RNA standards of known concentration and sequence. These synthetic RNAs cover a 106 range of concentration, as well as varying in length and GC content to allow for validation of dose-response and the fidelity of the procedure in downstream analyses [39]. Libraries were multiplexed, pooled, and sequenced on multiple lanes of an Illumina HiSeq2500, generating median 74.5 M paired-end reads per library of 100 bp.

RNA-seq data processing and analysis

Quality checking of sequence reads was assessed using fastQC (v.0.10.1) ( Sequence reads were aligned to the human reference genome (GRCh37, Ensembl build v. 75) using STAR (version 2.5.2a) with parameters ‘–outSAMunmapped Within –outFilterMultimapNmax 1–outFilterMismatchNoverLmax 0.1–alignIntronMin 21–alignIntronMax 0–alignEndsType Local–quantMode GeneCounts–twopassMode Basic’ [40]. STAR aligner also generated gene-level counts for all libraries relying on the human gene annotation provided for Ensembl GRCh37, build 75. Based on quality checking of alignments assessed by custom scripts utilizing PicardTools (, RNASeQC [41], RSeQC [42], and SamTools [43], no outlier sample was identified. Differentially expressed genes (DEGs) in three pair-wise comparisons including Het vs. WT, Null_A vs WT, and Null_B vs WT were identified by edgeR’s quasi-likelihood F test (v. 3.18.1) [44], testing differential expression between selected TSC1-mutant samples, Het or Null, and WT samples, which was run at the R platform (v. 3.4). In differential expression analysis, genes that passed the expression detection threshold, which was determined to be > 5 based on ERCC analysis as described in [45], in at least half of six analyzed samples in a given comparison were analyzed. Further comparison analysis revealed that 107 DEGs with Bonferroni adjusted p values < 0.05 overlapped among three pair-wise comparisons, of which 29 DEGs (9 upregulated, 20 downregulated) showed dosage effect in which significantly up- or downregulated DEGs with more than 2-fold changes in Het vs. WT comparison showed at least 1.5-fold up or downregulation in each Null vs. WT comparisons compared to their fold change in Het vs. WT. Gene ontology (GO) enrichment analysis was performed separately for 9 up-regulated and 20 downregulated DEGs represented by Ensembl gene IDs, using R/Bioconductor topGO package (v. 2.28) in R with “weight01” algorithm and “fisher” statistics and nodeSize = 10 relying on genome-wide annotations for Human in R/Bioconductor package (v. 3.4.1). In the GO enrichment analysis, union of all the analyzed genes in each of three comparisons, containing 22,235 genes were used as a background gene set.


Generation and characterization of an isogenic set of TSC1-iPSC lines

We established a TSC patient-derived iPSC line by reprogramming skin fibroblasts bearing a truncating nonsense mutation in exon 15 of TSC1 (1746C>T, Arg509X). We used a non-viral, synthetic modified mRNA-based method eliminating the risk of genomic integration and/or mutagenesis inherent to DNA-based approaches [46]. iPSCs showed typical colony formation in three independent clones (Fig. 1a), a normal karyotype (Fig. 1b), an expected decrease in the expression of pluripotency markers (OCT4, NANOG, and REX1) upon differentiation to embryoid bodies (EB). Further EB assays performed as recently described to assess the differentiation potential of iPSCs [47] revealed an increase in the expression of the three germ layer markers representing ectoderm (GFAP), endoderm (AFP), and mesoderm (Brachyury) (Fig. 1c, 3 biological replicates).

Fig. 1

Generation and characterization of isogenic TSC1-iPSCs. a Bright field images of a TSC1-iPSC colony (right) generated from skin fibroblasts of a TSC1 patient (left). b Representative normal karyotype of heterozygous TSC1-iPSCs. A total of 20 cells were counted to confirm normal diploidy of 46 and a total of 8 cells were analyzed in which the chromosomes were compared band by band to their homologues and a total of 4 were karyotyped. c Upon differentiation, real-time PCR in embryoid bodies (EB) from TSC1-iPS cells show decreased pluripotency markers (OCT4, NANOG, and REX1) at 8 days post-differentiation (8dEB) versus undifferentiated iPSCs and increased expression of the 3 germ layers markers (GFAP for ectoderm, AFP for endoderm, and Brachyury for mesoderm). Error bars represent standard deviation on 3 biological replicates. Data were normalized to the undifferentiated control. Mean values are shown, **p < 0.001 calculated with Student’s t test d CRISPR/Cas9-mediated correction of TSC mutation. Sequencing of original and corrected control iPSC lines derived from a TSC1 patient. e Quantitative RT-PCR shows increased TSC1 (left) expression in CRISPR-corrected iPSCs (Corr-WT) compared to heterozygous (Het) original iPSCs. Data were normalized to the Corr-WT. Mean values are shown, error bars represent standard deviation of 3 biological replicates, **p < 0.001 calculated with Student’s t test. No change in expression of TSC2 was noted. f Isogenic iPSCs heterozygous (Het) and corrected (Corr-WT), immunostained for pluripotency markers OCT4 (green) and NANOG (NAN, red). Bright field and nuclear DAPI (blue) are also shown. Scale bar = 200 μm. g Schematic representation of the two TSC1-Null clones of iPSCs generated using CRISPR/Cas9. h Immunoblot of TSC1 iPSC lines compared to an unrelated TSC iPSCs control (CT). Expression of TSC1 is reduced in TSC1-Het clones and completely lost after introducing a second somatic mutation, and mTORC1 is activated in both the Het and Null clones as shown by elevated expression of pS6. ac Independent iPSC clones of TSC1-Het (TSC1+/−) or TSC1-Null (TSC1−/−). Expression of TSC2, TSC1, phosphorylated S6 (pS6), and total amount of S6 (S6) were quantified from at least 3 independent experiments. Data were normalized to the Corr-WT. Mean values ± S.D. of three separate experiments are shown, *p < 0.01, **p < 0.001 calculated with Student’s t test

We then used the CRISPR/Cas9 technique to correct the mutation in TSC1-Het iPSCs by using a mutant specific sgRNA and a single-stranded oligo donor (ssODN) (Fig. 1d and Additional file 3: Table S2). The resulting wildtype iPSC line with the same genetic background as the original patient-derived TSC1-Het iPSCs is referred to as corrected wildtype (Corr-WT) and is used as a control. Real-time PCR analysis of Corr-WT TSC1-iPSCs showed increased expression of TSC1 compared to the original TSC1-Het line (Fig. 1e). We tested TSC2 mRNA expression since loss of TSC1 could affect the expression of TSC2 and observed no significant difference in expression of TSC2 mRNA between TSC1-Het and Corr-WT (Fig. 1e), all performed in three biological replicates. iPSC colonies were assessed for pluripotency by immunostaining and showed normal expression of OCT4 and NANOG (Fig. 1f). We also confirmed that the differentiation capacities of the Corr-WT iPSCs were intact by observing a decrease in pluripotent marker expression during embryoid body formation, and cytogenetic analysis revealed normal karyotyping (data not shown). To represent a second somatic mutation in TSC1, employing the CRISPR/Cas9 method again, we introduced mutations in the TSC1-Het iPSC line by targeting exon 7 using sgRNA cloned in pSpCas9(BB)-2A-Puro (PX459) vector (Additional file 3: Table S2). Two independent clones A and B with compound heterozygous mutations in TSC1 exon 7 (Fig. 1g) and showing normal OCT4 expression (Additional file 4: Figure S1) were chosen for further investigation. Immunoblotting for TSC1/hamartin protein in Corr-WT (TSC1+/+) showed expression level similar to another unrelated control iPSC line (CT), while the expression was reduced in three independent TSC1-Het iPSC clones (A, B, and C) and completely lost in TSC1-Null mutant iPSC lines. TSC2 expression was not changed significantly in either TSC1-Het or TSC1-Null iPSC lines when compared with the Corr-WT. As predicted, mTORC1 signaling was activated as shown by elevated expression of readout phosphorylated S6 (pS6) in both TSC1-Het and Null iPSC lines (Fig. 1h). These results were obtained on at least three independent experiments and confirmed the successful generation of an isogenic set of iPSCs representing TSC1 Corr-WT, Het, and Null, and we chose one clone for each genotype for all subsequent experiments.

Differentiation of isogenic TSC1-iPSC lines into NPCs

To generate a pure population of stable and expandable NPCs from iPSCs, we adapted a recently published protocol of directed monolayer differentiation using microbead sorting based on the expression of neural cell surface protein markers [36]. Here, we first pre-selected iPSCs based on the expression of an immature embryonic stem cell surface pluripotency marker (TRA-1-60) to maximize efficiency. We next selected PSA-NCAM-positive (PSA-NCAM+) cells to enrich for developing neuronal lineage cell populations and then sorted for CD133+/CD271 cells to isolate NPCs and eliminate the neural crest cell population. Following microbead purification, the resulting isogenic set of iPSC-derived TSC1-NPCs was confirmed by immunofluorescent staining for neural markers SOX2 and NESTIN at least 3 times for each genotype. (Fig. 2a).

Fig. 2

Characterization of TSC1 iPSC-derived NPCs. a All TSC1 NPC lines (Het, Null, Corr-WT) express expected neural progenitor markers SOX2 (upper panel) and NESTIN (panel below). DAPI in blue, SOX2 in green, and NESTIN in green. Scale bar = 100 μm. Immunostaining was performed at least 3 times. b TSC1-Het and Null NPCs display increased cell size compared to TSC1-Corr-WT as shown in bright field images (a) and by forward scatter FACS analysis; n = 3. c As expected, TSC1-Null and -Het NPCs show dose-dependent increased mTORC1 signaling (pS6 readout) compared to Corr-WT. Protein expression was quantified and normalized to the Corr-WT NPCs, n = 6, mean values ± s.e.m. are shown, *p < 0.01, **p < 0.001 calculated with Student’s t test. d Proliferation rate of NPC lines was quantified at day 0 (D0, equal cell seeding), and live cell numbers were assessed at D2, D3, and D5. Mutant TSC1 NPCs (Het and Null) revealed genotype-dependent increased proliferation compared to Corr-WT. Data was normalized to Corr-WT at D0, mean values ± S.D. of three separate experiments are shown,*p < 0.01, **p < 0.001 calculated with Student’s t test. e, f MAP2 immunostaining showed genotype-dependent increased neurite outgrowth (number and length) in TSC1 mutant NPCs, which were quantified using a custom image analysis pipeline and HCA Vision imaging software creating neurite segmentation (representative panel shown for the DMSO-treated NPCs). Analysis on n = 6 field images per treatment with approximately 50 cells per field. Data normalized to DMSO treated Corr-WT NPCs. Mean values + s.e.m. are shown. *p < 0.05, **p < 0.001, ****p < 0.0001, n.s. = not significant, calculated with Welch’s t test (GraphPad Prism 7.05). g TSC1 NPCs (Corr-WT, Het, and Null-clone B) were treated with 100 nM of rapamycin for 24 h or with DMSO and analyzed by flow cytometry using the forward scatter height (FSC-H) gating. For each cell line the blue curve represents the rapamycin-treated cells and the red curve represents the DMSO control cells. A shift in the curves shows a cell size difference. N = 3. h Proliferation rate of NPC lines after treatment with a vehicle control (DMSO) or rapamycin (100nM) was quantified at day 0 (D0, equal cell seeding), and live cell numbers were assessed at D2, D3, D4, and D5. No significant differences were observed between DMSO or rapamycin-treated NPCs in all cell types at all time point; n = 3. Mean values ± S.D. of three separate experiments are shown, data was normalized to the Corr-WT treated with DMSO at D0

TSC1 iPSC-derived NPCs show genotype-dependent phenotypes

We examined morphological differences between TSC1-CorrWT, TSC1-Het, and TSC1-Null NPCs, and found that TSC1-Het and Null NPCs were larger in size than the Corr-WT on 3 biological replicates and a representative image is shown in Fig. 2b. This is consistent with the reduction or loss of TSC1 in the Het and Null respectively, inducing a strong activation of mTORC1 as seen on at least six independent replicates by higher phospho-ribosomal protein S6 (pS6) expression levels in a dose-dependent manner (Fig. 2c). TSC2 expression is included as a control (Fig. 2c). We also observed that both the TSC1-Het and Null NPCs proliferate faster by day 5 when compared with the matched Corr-WT as determined by viable cell counts using trypan blue exclusion and automated cell counting on three independent experiments (Fig. 2d).

Previous studies of neurodevelopmental disorders including Rett Syndrome, Fragile-X Syndrome, TSC, ASD, and schizophrenia have used post-mitotic neurons derived from human iPSCs or mouse models to study morphological aspects such as dendrite outgrowth or synapse formation, or functional characteristics using neuronal electrophysiology, which represents late-stage neurodevelopmental events. However, NPC proliferation and initial process extension phenotypes, such as neurite number and length represent earlier events in neurogenesis and are strongly implicated in ASD [29, 30]. Therefore to explore the potential early neurodevelopmental deficits, we examined the isogenic NPC set using phalloidin and MAP2 immunostaining to quantify neurite length, number and branching points and to compare the genotype-specific differences. Interestingly, TSC1-Het and Null NPCs revealed a significant increase in average neurite number as well as average neurite length per cell (n = 6 field images for each genotype with approximately 50 cells per field). (Fig. 2e, f). Rapamycin treatment did not have an effect on neurite number or length (n = 3) (Fig. 2f). Rapamycin treatment (100 nM) for 24 h reduced cell size in TSC1-Het and Null (Fig. 2g) but had no effect on proliferation as shown by viable cell counts (n = 3) (Fig. 2h) and confirmed by flow cytometric analysis of cell cycle using quantification of propidium iodide DNA staining or a Cell-Titer-Glo viability assay that assess cellular ATP levels (Additional file 1: Supplementary Materials and Methods; Additional file 4: Figure S2).

These data suggest that early neurodevelopmental phenotypes, such as proliferation and neurite outgrowth that occur prior to neuronal differentiation are altered in TSC1-Het and Null NPCs when compared with the matched WT and that mTORC1 inhibition though rapamycin treatment does not affect these phenotypes.

TSC1-mutant NPCs reveal activation of MEK-ERK signaling

We have previously reported MEK-ERK signaling to be aberrantly elevated in TSC patient subependymal giant cell astrocytoma cells (SEGAs) compared to normal brain [48]. Therefore, we investigated whether TSC1-mutant NPCs had similar dysregulation of this pathway. Interestingly, we observed at least in five independent experiments an increase in pERK1/2 in TSC1-Het and TSC1-Null NPCs compared to Corr-WT (Fig. 3a), which is consistent with our previous report and differs from Tsc1/2-Null MEFs where ERK signaling is downregulated due to feedback inhibition by Grb10 [49]. Our results suggest that in addition to mTORC1 signaling, MEK-ERK signaling is upregulated in TSC1-mutant human neuronal cell types. Rapamycin treatment, as expected, blocked mTORC1 activation as detected by reduced phospho-ribosomal protein S6 kinase (pS6K) and pS6 levels, (Fig. 3a). Rapamycin however significantly increased pERK1/2 in TSC1-Het and Null NPCs (Fig. 3a).

Fig. 3

Activation of MEK/ERK and MNK-eIF4E pathways in TSC1 mutant NPCs. a In addition to mTORC1/2 activation, immunoblotting also showed Increased phosphorylation of ERK1/2 (pERK1/2) in TSC1-Het and Null NPCs. As predicted, rapamycin (100 nM, 24 h) inhibited mTORC1 signaling (pS6K and pS6 readouts). Rapamycin increased pERK1/2 in TSC1-Null NPCs; n = 5. b Treatment with rapamycin (100 nM, 24 h) leads to increased phosphorylation of translational subunit eIF4E (p-eIF4E) in TSC1-NPCs; n = 6. c Treatment with a MEK inhibitor tramatinib (250 nM, 24 h) alone or combined with rapamycin (100 nM, 24 h) inhibits phosphorylation of ERK1/2 (pERK1/2) but has no effect on the translational subunit eIF4E; n = 3. d Increased p-eIF4E in TSC1-Het and Null NPCs is blocked using MNK inhibitor eFT508 (50 nM, 24 h), alone or combined with rapamycin (100 nM, 24 h); n = 3. Mean values ± S.D. of three separate experiments are shown, *p < 0.01, **p < 0.001 calculated with Student’s t test

Recent studies demonstrated that a single, conserved Ser residue (S209) in the Eukaryotic Translation Initiation Factor 4E (eIF4E), phosphorylated by MAP kinase-interacting kinase (MNK) [50] not only plays a role in cancer biology, but also in neurobiology by regulating 5’ cap-dependent translation of specific mRNAs in neuronal cells [51]. Rapamycin is known to increase eIF4E phosphorylation at S209 in many cancer cell types [50, 52, 53]. To our knowledge, the phospho-status of eIF4E has not been assessed with TSC1/2-deficiency. In our TSC1 NPCs, we found no significant basal activation; however, upon treatment with rapamycin, we observed upregulation of p-eIF4E S209 (n = 6) (Fig. 3b). These results suggest that mTORC1 inhibition in TSC1-mutant NPCs leads to increases in both MEK-ERK and MNK-eIF4E signaling pathways, both of which are known to regulate protein translation [6, 51, 54].

To understand the mechanism of rapamycin-induced activation of ERK1/2 and eIF4E in TSC1 NPCs, we examined the effects of a MEK inhibitor trametinib, and a MNK inhibitor eFT508 (tomivosertib), which is currently in clinical development, either alone or in combination with rapamycin. Trametinib treatment inhibited phosphorylation of ERK1/2 as expected, but had no effect on p-eIF4E (n = 3) (Fig. 3c). Conversely, eFT508 treatment completely inhibited p-eIF4E but not pERK1/2 (n = 3) (Fig. 3d). These results suggest that phosphorylation of eIF4E by MNK is not regulated by MEK-dependent ERK signaling in TSC1 NPCs.

Rapamycin induced activation of eIF4E is partly dependent on PI3K activation

We observed that the increase in p-eIF4E after rapamycin treatment was more pronounced in TSC1-Null NPCs that have strong activation of mTORC1 (n = 3) (Fig. 3c, d), suggesting that inhibition of mTORC1 may relieve the negative feedback regulation on PI3K signaling, resulting in an increase in p-eIF4E. Further, studies performed in human cancer cells have shown that MNK-dependent eIF4E phosphorylation is regulated by PI3K signaling [52, 53]. Therefore, we examined the effects of a PI3K inhibitor wortmannin and observed that increased eIF4E phosphorylation after rapamycin treatment is partly dependent on PI3K signaling (n = 3) (Fig. 4a). Similarly, rapamycin-induced increase in pERK1/2, particularly in TSC1-Null NPCs was partially inhibited by wortmannin supporting a PI3K-ERK signaling axis in TSC1 mutant NPCs (n = 3) (Fig. 4b).

Fig. 4

Rapamycin induced activation of eIF4E and ERK1/2 is partly dependent on PI3K activation. a Treatment with wortmannin (500 nM or 1 μM for 3 h) combined with rapamycin (100 nM for 24 h) in TSC1-Null NPCs decreases phosphorylation of eIF4E; n = 3. b Treatment with wortmannin (1 μM for 3 h) alone or in combination with rapamycin (100 nM for 24 h) inhibits pAKT (s473) and decreases phosphorylation of ERK1/2 in TSC1 mutant NPCs; n = 3. pS6 serves as a control. For each panel, protein expression was quantified and normalized to the Corr-WT NPCs treated with DMSO. Mean values ± S.D. of three separate experiments are shown, *p < 0.01, **p < 0.001 calculated with Student’s t test

Transcriptome analyses of TSC1 isogenic NPCs

To further characterize the NPCs at the transcriptome level, we performed RNA-seq of TSC1 NPCs representing Corr-WT, Het, and two independent Null clones (Null-A and Null-B) in triplicates. Differential expression analyses of these samples identified distinct and shared upregulated or downregulated genes between TSC1-Het versus Corr-WT, and TSC1-Null (A and B) versus Corr-WT (Fig. 5a–c). We observed a total of 107 differentially expressed genes that overlapped between Het and Null, when compared with the Corr-WT at a high stringency of analyses (Bonferroni adjusted p < 0.05) (Fig. 5d). It is noteworthy that of the 107 genes, 29 displayed a genotype-dependent linear response, i.e., genes upregulated or downregulated in TSC1-Het NPCs that were further increased or decreased in TSC1-Null NPCs, respectively. In particular, genes such as ANXA1, CNTN6, HLA-B, PCDH19, and PCDH10, have been linked to ASD, epilepsy, ID, and other neuropsychiatric disorders [55,56,57,58,59,60], were significantly upregulated or downregulated and warrant further investigation (Fig. 5e). In addition, several genes encoding members of the Zinc-finger protein (ZNF) family of sequence-specific, DNA binding transcription factors [61] were downregulated in TSC1-Het and Null NPCs raising the possibility of their role in neurodevelopment. Gene ontology analyses for TSC1-Het and Null NPCs revealed enriched categories for upregulated genes including endosome membrane, endoplasmic reticulum (ER), and ER-to-Golgi transport vesicle membrane, whereas downregulated genes revealed enrichment related to DNA binding and regulation of transcription (Fig. 5f). Taken together, these data support the notion that loss of TSC1 causes aberrant changes in the transcriptome and in changes in pathways critical for shaping the neural proteome relevant for early neurodevelopment.

Fig. 5

Transcriptome analysis of TSC1 NPCs reveals genotype-dependent changes in gene expression. a–c Differential expression analysis identified distinct and shared upregulated and downregulated genes in the isogenic TSC1 NPC set for Het or 2 independent Null clones (a, b) versus Corr-WT. d Of the 107 genes shared between Het and Null NPCs, 29 showed a genotype-dependent linear response (showed on the heat map: e). This subset included ASD-related genes (red, up; blue, down) as well as multiple members of the Zinc-finger gene family of transcription factors (*). f Gene ontology analysis for the genes that showed a linear response


Neurodevelopmental syndromes including ASD, Fragile X, and TSC are commonly considered disorders of synaptic homeostasis and therefore, several prior studies have focused on differentiated neurons and synaptic defects. However, emerging studies reveal early neurodevelopmental events such as NPC proliferation, neurite outgrowth, and migration that precede synaptogenesis also play a vital role in disease pathogenesis of ASD and other neuropsychiatric disorders [27,28,29,30,31, 35, 62]. Further, exome sequencing studies in ASD as well as network analyses of large numbers of ASD-implicated genes also indicate that defects in neural progenitor cell division could be a shared phenomenon regulated by these genes [63,64,65,66]. Recent studies in TSC that employed either human embryonic stem cell lines with heterozygous or homozygous loss of TSC2 or TSC patient iPSC-derived neurons with heterozygous loss of TSC2 confirmed that mTORC1 inhibition corrects synaptic defects [22, 25]. However, while the treatment of TSC patients with an mTORC1 inhibitor manages tumor growth, its effectiveness in treating TSC-associated neuropsychiatric defects has remained equivocal [11]. Clinical studies and mouse models provide strong evidence to support the hypothesis that the NPC lineage is the cell of origin for the CNS manifestations of TSC, and the neural crest cell (NCC) lineage is responsible for other aspects of TSC [67]. To further elucidate these mechanisms, here we have generated a cellular model using NPCs to study early neurodevelopmental aspects of TSC. With an isogenic set of NPCs (Corr-WT, Het, and Null) derived from TSC patient iPSCs with heterozygous loss of TSC1, we focused our characterization on neurogenesis phenotypes such as NPC proliferation and neurite outgrowth. Our results convincingly show enhanced proliferation in TSC1-Het and Null NPCs when compared with the isogenic control (Corr-WT), which is consistent with previous reports [23, 25]. However, unlike the previous studies, rapamycin treatment did not have an effect on enhanced proliferation and neurite outgrowth in TSC1-Het and Null NPCs when compared to the Corr-WT, suggesting that early neurodevelopmental phenotypes seen upon loss of TSC1 are not solely dependent on mTORC1 activation. Thus, our isogenic TSC1 NPC cell models provide an opportunity to screen for drugs that could reverse the early neurodevelopmental phenotypes such as NPC proliferation and neurite outgrowth, which may ultimately lead to better treatment for TSC-associated epilepsy and neuropsychiatric defects.

MNK-dependent phosphorylation of eIF4E S209 plays a role in neurobiology by regulating the translation of specific mRNAs in neuronal cells [51]. It is well established that the mTORC1 signaling pathway is critically involved in protein translation through regulation of initiation. mTORC1 activation results in hyperphosphorylated 4E-BP, which in turn releases the bound eIF4E and facilitates the assembly of eukaryotic initiation factor (eIF) 4F, a heterotrimeric complex composed of eIF4E, a cap-binding protein; eIF4A, an RNA helicase; and eIF4G, a large scaffolding protein that recruits 40S ribosomes to mRNA templates (Fig. 6a). Our results suggest that in TSC1 mutant NPCs, while rapamycin inhibits translation regulated by mTORC1 through phospho-4E-BP, it possibly enhances translation mediated by MEK-ERK and MNK-eIF4E pathways, which couples synaptic activity to the translational machinery and plays an essential role in neuronal translation. Our results also reveal that rapamycin-induced activation of MNK-eIF4E could be mediated by PI3K (Fig. 6b), consistent with prior studies in other cancer cells [52, 53]. Interestingly, the MNK-eIF4E axis controls the translation of mRNAs and synaptic plasticity through regulation of the translational repressor, Cytoplasmic Fragile X Protein-Interacting Protein 1 (CYFIPI) [51, 68, 69]. Taken together, we believe that neuronal translation mediated by MEK-ERK and MNK-eIF4E could play a role in CNS aspects of TSC. Moreover, treatment with rapalogs, the current standard of care for TSC patients, may lead to adaptive changes in neuronal cells via signaling feedback mechanisms, and future translatome (referring to mRNAs being actively translated) studies performed before and after mTOR and MNK inhibitor treatment in TSC patient-derived NPCs could be valuable in providing in-depth information regarding distinct translational changes.

Fig. 6

Model depicting downstream signaling pathways in untreated and rapamycin-treated TSC1-mutant NPCs. a In untreated TSC1-mutant NPCs, loss of the TSC protein complex results in activated MEK-ERK and mTORC1 signaling. Activated mTORC1 leads to phosphorylation and release of the inhibitory 4EBP1 from the eIF4E subunit enabling the formation of an active eIF4F complex. b Upon rapamycin (Rapa) treatment, while mTORC1-dependent eIF4E activation is inhibited, an alternate mechanism of MNK-mediated eIF4F phosphorylation/activation is enhanced, which is blocked by co-treatment with the MNK inhibitor eFT508

It is noteworthy that our transcriptome studies on the isogenic sets of TSC1 NPCs revealed genotype-dependent alterations in the expression of 29 genes suggesting they may play an essential role in TSC pathogenesis. Genes implicated in ASD such as Annexin 1 (ANXA1) (50), HLA-B [55, 70] and CNTN6 with deletion or duplication in a spectrum of neurodevelopmental disorders and ID [71, 72] were upregulated in TSC1-Het and Null NPCS when compared with the WT. The protocadherin genes, PCDH10 and PCDH19 that are linked to ASD and epilepsy [58, 73, 74] were downregulated along with several members of Zinc-finger gene family of transcription factors. Follow-up studies are necessary to define the roles of these genes in neurodevelopmental phenotypes that we observe in NPCs and to understand whether expression changes correlate at the protein level and could be dependent on mTORC1 activation.


A limitation of this study is not understanding the role of genes linked to ASD, epilepsy, and ID, which were identified as significantly upregulated or downregulated in NPCs with heterozygous and homozygous loss of TSC1. The undertaking of such studies is beyond the scope of this work and follow up studies are necessary.


Our results clearly establish that both heterozygous and homozygous loss of TSC1 influence early neurodevelopmental phenotypes, signaling, and gene expression in NPCs compared to the genetically matched WT cells. Our approach of using TSC patient’s iPSC-derived NPCs will provide a useful platform for large-scale omics and drug screening studies that may identify drugs that could be superior to rapamycin or effective in combination with rapamycin to ultimately treat TSC-associated neuropsychiatric symptoms.

Availability of data and materials

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



Autism spectrum disorder


Central nervous system


Corrected wildtype


Clustered regularly interspaced short palindromic repeats / CRISPR associated protein 9


Cytoplasmic fragile X protein-interacting protein 1


Eukaryotic translation initiation factor 4E


Endoplasmic reticulum


Extracellular signal-regulated kinase




Intellectual disability


Induced pluripotent stem cells


MAP kinase-interacting kinase


Mammalian/mechanistic target of rapamycin complex 1


Neural crest cells


Neural progenitor cells


Phosphoinositide 3-kinase


Phospho-ribosomal protein S6


Phospho-ribosomal protein S6 kinase


RNA sequencing


Subependymal giant cell astrocytoma


Single-stranded oligo donor


Tuberous sclerosis complex


Zinc-finger protein


  1. 1.

    Caban C, Khan N, Hasbani DM, Crino PB. Genetics of tuberous sclerosis complex: implications for clinical practice. Appl Clin Genet. 2017;10:1–8.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  2. 2.

    Henske EP, Jozwiak S, Kingswood JC, Sampson JR, Thiele EA. Tuberous sclerosis complex. Nat Rev Dis Primers. 2016;2:16035.

    PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Salussolia CL, Klonowska K, Kwiatkowski DJ, Sahin M. Genetic Etiologies, Diagnosis, and Treatment of Tuberous Sclerosis Complex. Annu Rev Genomics Hum Genet. 2019.

  4. 4.

    Han JM, Sahin M. TSC1/TSC2 signaling in the CNS. FEBS letters. 2011;585(7):973–80.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  5. 5.

    Huang J, Manning BD. The TSC1-TSC2 complex: a molecular switchboard controlling cell growth. Biochem J. 2008;412(2):179–90.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  6. 6.

    Kelleher RJ 3rd, Bear MF. The autistic neuron: troubled translation? Cell. 2008;135(3):401–6.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  7. 7.

    Bissler JJ, McCormack FX, Young LR, Elwing JM, Chuck G, Leonard JM, et al. Sirolimus for angiomyolipoma in tuberous sclerosis complex or lymphangioleiomyomatosis. N Engl J Med. 2008;358(2):140–51.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  8. 8.

    Davies DM, de Vries PJ, Johnson SR, McCartney DL, Cox JA, Serra AL, et al. Sirolimus therapy for angiomyolipoma in tuberous sclerosis and sporadic lymphangioleiomyomatosis: a phase 2 trial. Clin Cancer Res. 2011;17(12):4071–81.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  9. 9.

    Julich K, Sahin M. Mechanism-based treatment in tuberous sclerosis complex. Pediatr Neurol. 2014;50(4):290–6.

    PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Krueger DA, Care MM, Holland K, Agricola K, Tudor C, Mangeshkar P, et al. Everolimus for subependymal giant-cell astrocytomas in tuberous sclerosis. N E J Med. 2010;363(19):1801–11.

    Article  CAS  Google Scholar 

  11. 11.

    Krueger DA, Sadhwani A, Byars AW, de Vries PJ, Franz DN, Whittemore VH, et al. Everolimus for treatment of tuberous sclerosis complex-associated neuropsychiatric disorders. Ann Clin Transl Neurol. 2017;4(12):877–87.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  12. 12.

    Sahin M, Henske EP, Manning BD, Ess KC, Bissler JJ, Klann E, et al. Advances and Future Directions for Tuberous Sclerosis Complex Research: Recommendations From the 2015 Strategic Planning Conference. Pediatr Neurol. 2016;60:1–12.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Anderson RH, Francis KR. Modeling rare diseases with induced pluripotent stem cell technology. Mol Cell Probes. 2018;40:52–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  14. 14.

    Connor B. Concise Review: The Use of Stem Cells for Understanding and Treating Huntington's Disease. Stem Cells. 2018;36:146–60.

    PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Li L, Chao J, Shi Y. Modeling neurological diseases using iPSC-derived neural cells : iPSC modeling of neurological diseases. Cell Tissue Res. 2018;371(1):143–51.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  16. 16.

    Liu J, Koscielska KA, Cao Z, Hulsizer S, Grace N, Mitchell G, et al. Signaling defects in iPSC-derived fragile X premutation neurons. Hum Mol Genet. 2012;21(17):3795–805.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  17. 17.

    Silva MC, Haggarty SJ. Human pluripotent stem cell-derived models and drug screening in CNS precision medicine. Ann N Y Acad Sci. 2019.

  18. 18.

    Musunuru K. Genome editing of human pluripotent stem cells to generate human cellular disease models. Dis Model Mech. 2013;6(4):896–904.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  19. 19.

    Reinhardt P, Schmid B, Burbulla LF, Schondorf DC, Wagner L, Glatza M, et al. Genetic correction of a LRRK2 mutation in human iPSCs links parkinsonian neurodegeneration to ERK-dependent changes in gene expression. Cell Stem Cell. 2013;12(3):354–67.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  20. 20.

    Tai DJ, Ragavendran A, Manavalan P, Stortchevoi A, Seabra CM, Erdin S, et al. Engineering microdeletions and microduplications by targeting segmental duplications with CRISPR. Nat Neurosci. 2016;19(3):517–22.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  21. 21.

    Blair JD, Hockemeyer D, Bateup HS. Genetically engineered human cortical spheroid models of tuberous sclerosis. Nat Med. 2018;24(10):1568–78.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  22. 22.

    Costa V, Aigner S, Vukcevic M, Sauter E, Behr K, Ebeling M, et al. mTORC1 inhibition corrects neurodevelopmental and synaptic alterations in a human stem cell model of tuberous sclerosis. Cell Rep. 2016;15(1):86–95.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  23. 23.

    Li Y, Cao J, Chen M, Li J, Sun Y, Zhang Y, et al. Abnormal neural progenitor cells differentiated from induced pluripotent stem cells partially mimicked development of TSC2 neurological abnormalities. Stem Cell Reports. 2017;8(4):883–93.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  24. 24.

    Nadadhur AG, Alsaqati M, Gasparotto L, Cornelissen-Steijger P, van Hugte E, Dooves S, et al. Neuron-glia interactions increase neuronal phenotypes in tuberous sclerosis complex patient iPSC-derived models. Stem Cell Reports. 2019;12(1):42–56.

    PubMed  Article  PubMed Central  Google Scholar 

  25. 25.

    Sundberg M, Tochitsky I, Buchholz DE, Winden K, Kujala V, Kapur K, et al. Purkinje cells derived from TSC patients display hypoexcitability and synaptic deficits associated with reduced FMRP levels and reversed by rapamycin. Mol Psychiatry. 2018;23:2167–83.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Winden KD, Sundberg M, Yang C, Wafa SMA, Dwyer S, Chen PF, et al. Biallelic Mutations in TSC2 Lead to Abnormalities Associated with Cortical Tubers in Human iPSC-Derived Neurons. J Neurosci. 2019;39(47):9294–305.

    PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Brennand K, Savas JN, Kim Y, Tran N, Simone A, Hashimoto-Torii K, et al. Phenotypic differences in hiPSC NPCs derived from patients with schizophrenia. Mol Psychiatry. 2015;20(3):361–8.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  28. 28.

    Kaushik G, Zarbalis KS. Prenatal neurogenesis in autism spectrum disorders. Front Chem. 2016;4:12.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  29. 29.

    Marchetto MC, Belinson H, Tian Y, Freitas BC, Fu C, Vadodaria K, et al. Altered proliferation and networks in neural cells derived from idiopathic autistic individuals. Mol Psychiatry. 2017;22(6):820–35.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  30. 30.

    Packer A. Neocortical neurogenesis and the etiology of autism spectrum disorder. Neurosci Biobehav Rev. 2016;64:185–95.

    PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Mellios N, Feldman DA, Sheridan SD, Ip JPK, Kwok S, Amoah SK, et al. MeCP2-regulated miRNAs control early human neurogenesis through differential effects on ERK and AKT signaling. Mol Psychiatry. 2018;23(4):1051–65.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  32. 32.

    Niida Y, Lawrence-Smith N, Banwell A, Hammer E, Lewis J, Beauchamp R, et al. Analysis of both TSC1 and TSC2 for germline mutations in 126 unrelated patients with tuberous sclerosis. Hum Mutat. 1999;14:412–22.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  33. 33.

    Schlaeger TM, Daheron L, Brickler TR, Entwisle S, Chan K, Cianci A, et al. A comparison of non-integrating reprogramming methods. Nat Biotechnol. 2015;33(1):58–63.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  34. 34.

    Warren L, Manos PD, Ahfeldt T, Loh YH, Li H, Lau F, et al. Highly efficient reprogramming to pluripotency and directed differentiation of human cells with synthetic modified mRNA. Cell Stem Cell. 2010;7(5):618–30.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  35. 35.

    Sheridan SD, Theriault KM, Reis SA, Zhou F, Madison JM, Daheron L, et al. Epigenetic characterization of the FMR1 gene and aberrant neurodevelopment in human induced pluripotent stem cell models of fragile X syndrome. PloS One. 2011;6(10):e26203.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Sellgren CM, Gracias J, Watmuff B, Biag JD, Thanos JM, Whittredge PB, et al. Increased synapse elimination by microglia in schizophrenia patient-derived models of synaptic pruning. Nat Neurosci. 2019;22(3):374–85.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  37. 37.

    James MF, Han S, Polizzano C, Plotkin SR, Manning BD, Stemmer-Rachamimov AO, et al. NF2/merlin is a novel negative regulator of mTOR complex 1, and activation of mTORC1 is associated with meningioma and schwannoma growth. Mol Cell Biol. 2009;29(15):4250–61.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  38. 38.

    Han S, Kim S, Bahl S, Li L, Burande CF, Smith N, et al. The E3 ubiquitin ligase protein associated with Myc (Pam) regulates mammalian/mechanistic target of rapamycin complex 1 (mTORC1) signaling in vivo through N- and C-terminal domains. J Biol Chem. 2012;287(36):30063–72.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  39. 39.

    Jiang L, Schlesinger F, Davis CA, Zhang Y, Li R, Salit M, et al. Synthetic spike-in standards for RNA-seq experiments. Genome Res. 2011;21(9):1543–51.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  Google Scholar 

  41. 41.

    DeLuca DS, Levin JZ, Sivachenko A, Fennell T, Nazaire MD, Williams C, et al. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics. 2012;28(11):1530–2.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  42. 42.

    Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28(16):2184–5.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  43. 43.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  44. 44.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  45. 45.

    Blumenthal I, Ragavendran A, Erdin S, Klei L, Sugathan A, Guide JR, et al. Transcriptional consequences of 16p11.2 deletion and duplication in mouse cortex and multiplex autism families. Am J Hum Genet. 2014;94(6):870–83.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  46. 46.

    Mandal PK, Rossi DJ. Reprogramming human fibroblasts to pluripotency using modified mRNA. Nat Protoc. 2013;8(3):568–82.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  47. 47.

    International Stem Cell I. Assessment of established techniques to determine developmental and malignant potential of human pluripotent stem cells. Nat Commun. 2018;9(1):1925.

    Article  CAS  Google Scholar 

  48. 48.

    Han S, Santos TM, Puga A, Roy J, Thiele EA, McCollin M, et al. Phosphorylation of tuberin as a novel mechanism for somatic inactivation of the tuberous sclerosis complex proteins in brain lesions. Cancer Res. 2004;64(3):812–6.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  49. 49.

    Yu Y, Yoon SO, Poulogiannis G, Yang Q, Ma XM, Villen J, et al. Phosphoproteomic analysis identifies Grb10 as an mTORC1 substrate that negatively regulates insulin signaling. Science. 2011;332(6035):1322–6.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  50. 50.

    Stead RL, Proud CG. Rapamycin enhances eIF4E phosphorylation by activating MAP kinase-interacting kinase 2a (Mnk2a). FEBS Lett. 2013;587(16):2623–8.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  51. 51.

    Bramham CR, Jensen KB, Proud CG. Tuning Specific Translation in Cancer Metastasis and Synaptic Memory: Control at the MNK-eIF4E Axis. Trends Biochem Sci. 2016;41(10):847–58.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  52. 52.

    Sun SY, Rosenberg LM, Wang X, Zhou Z, Yue P, Fu H, et al. Activation of Akt and eIF4E survival pathways by rapamycin-mediated mammalian target of rapamycin inhibition. Cancer Res. 2005;65(16):7052–8.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  53. 53.

    Wang X, Yue P, Chan CB, Ye K, Ueda T, Watanabe-Fukunaga R, et al. Inhibition of mammalian target of rapamycin induces phosphatidylinositol 3-kinase-dependent and Mnk-mediated eukaryotic translation initiation factor 4E phosphorylation. Mol Cell Biol. 2007;27(21):7405–13.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  54. 54.

    Kelleher RJ 3rd, Govindarajan A, Jung HY, Kang H, Tonegawa S. Translational control by MAPK signaling in long-term synaptic plasticity and memory. Cell. 2004;116(3):467–79.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  55. 55.

    Al-Hakbany M, Awadallah S, Al-Ayadhi L. The Relationship of HLA Class I and II Alleles and Haplotypes with Autism: A Case Control Study. Autism Res Treat. 2014;2014:242048.

    PubMed  PubMed Central  Google Scholar 

  56. 56.

    Correia CT, Conceicao IC, Oliveira B, Coelho J, Sousa I, Sequeira AF, et al. Recurrent duplications of the annexin A1 gene (ANXA1) in autism spectrum disorders. Mol Autism. 2014;5(1):28.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  57. 57.

    De Rubeis S, He X, Goldberg AP, Poultney CS, Samocha K, Cicek AE, et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature. 2014;515(7526):209–15.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  58. 58.

    Morrow EM, Yoo SY, Flavell SW, Kim TK, Lin Y, Hill RS, et al. Identifying autism loci and genes by tracing recent shared ancestry. Science. 2008;321(5886):218–23.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  59. 59.

    Smith L, Singhal N, El Achkar CM, Truglio G, Rosen Sheidley B, Sullivan J, et al. PCDH19-related epilepsy is associated with a broad neurodevelopmental spectrum. Epilepsia. 2018;59(3):679–89.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  60. 60.

    Yao BZ, Yu SQ, Yuan H, Zhang HJ, Niu P, Ye JP. The role and effects of ANXA1 in temporal lobe epilepsy: A Protection Mechanism? Med Sci Monit Basic Res. 2015;21:241–6.

    PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Cassandri M, Smirnov A, Novelli F, Pitolli C, Agostini M, Malewicz M, et al. Zinc-finger proteins in health and disease. Cell Death Discov. 2017;3:17071.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  62. 62.

    Williams M, Prem S, Zhou X, Matteson P, Yeung PL, Lu CW, et al. Rapid detection of neurodevelopmental phenotypes in human neural precursor cells (NPCs). J Vis Exp. 2018;133.

  63. 63.

    Cotney J, Muhle RA, Sanders SJ, Liu L, Willsey AJ, Niu W, et al. The autism-associated chromatin modifier CHD8 regulates other autism risk genes during human neurodevelopment. Nat Commun. 2015;6:6404.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  64. 64.

    Pinto D, Delaby E, Merico D, Barbosa M, Merikangas A, Klei L, et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am J Hum Genet. 2014;94(5):677–94.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  65. 65.

    Sugathan A, Biagioli M, Golzio C, Erdin S, Blumenthal I, Manavalan P, et al. CHD8 regulates neurodevelopmental pathways associated with autism spectrum disorder in neural progenitors. Proc Natl Acad Sci U S A. 2014;111(42):E4468–77.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  66. 66.

    Wilkinson B, Grepo N, Thompson BL, Kim J, Wang K, Evgrafov OV, et al. The autism-associated gene chromodomain helicase DNA-binding protein 8 (CHD8) regulates noncoding RNAs and autism-related genes. Transl Psychiatry. 2015;5:e568.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  67. 67.

    Delaney SP, Julian LM, Stanford WL. The neural crest lineage as a driver of disease heterogeneity in tuberous sclerosis complex and lymphangioleiomyomatosis. Front Cell Dev Biol. 2014;2:69.

    PubMed  PubMed Central  Article  Google Scholar 

  68. 68.

    Genheden M, Kenney JW, Johnston HE, Manousopoulou A, Garbis SD, Proud CG. BDNF stimulation of protein synthesis in cortical neurons requires the MAP kinase-interacting kinase MNK1. J Neurosci. 2015;35(3):972–84.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  69. 69.

    Napoli I, Mercaldo V, Boyl PP, Eleuteri B, Zalfa F, De Rubeis S, et al. The fragile X syndrome protein represses activity-dependent translation through CYFIP1, a new 4E-BP. Cell. 2008;134(6):1042–54.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  70. 70.

    Puangpetch A, Suwannarat P, Chamnanphol M, Koomdee N, Ngamsamut N, Limsila P, et al. Significant association of HLA-B alleles and genotypes in Thai children with autism spectrum disorders: a case-control study. Dis Markers. 2015;2015:724935.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  71. 71.

    Hu J, Liao J, Sathanoori M, Kochmar S, Sebastian J, Yatsenko SA, et al. CNTN6 copy number variations in 14 patients: a possible candidate gene for neurodevelopmental and neuropsychiatric disorders. J Neurodev Disord. 2015;7(1):26.

    PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Kashevarova AA, Nazarenko LP, Schultz-Pedersen S, Skryabin NA, Salyukova OA, Chechetkina NN, et al. Single gene microdeletions and microduplication of 3p26.3 in three unrelated families: CNTN6 as a new candidate gene for intellectual disability. Mol Cytogenet. 2014;7(1):97.

    PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

    Kolc KL, Sadleir LG, Scheffer IE, Ivancevic A, Roberts R, Pham DH, et al. A systematic review and meta-analysis of 271 PCDH19-variant individuals identifies psychiatric comorbidities, and association of seizure onset and disease severity. Mol Psychiatry. 2019;24(2):241–51.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  74. 74.

    Liu A, Yang X, Yang X, Wu Q, Zhang J, Sun D, et al. Mosaicism and incomplete penetrance of PCDH19 mutations. J Med Genet. 2019;56(2):81–8.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

Download references


Not applicable


This work was supported by US Army Medical Research Tuberous Sclerosis Complex Research program, W81XWH-15-1-0414 (VR, SJH, and PM), NIH R01 NS109540 (VR, SJH, and PM) MGH Executive Committee on Research Interim Support Funding (VR), and Stuart & Suzanne Steele MGH Research Scholars Program (SJH).

Author information




VR contributed to the conceptualization. PM, VW, SJH, and VR contributed to the methodology. SE and MT contributed to the formal analysis. PM, VW, RLB, and GS helped in the investigation. SAR, SDS, ET, SJH contributed to the resources. PM, SE, and VR helped in writing—original draft. PM, SE, RLB, SDS, SJH, and VR contributed to the writing review and editing. PM, VW, RLB, SAR, and SE contributed to the visualization. VR supervised the study and was responsible for the project administration. SJH and VR acquired the funding. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Vijaya Ramesh.

Ethics declarations

Ethics approval and consent to participate

This study of using patient-derived iPSCs was determined as not human subjects research by the Institutional Review Board.

Consent for publication

Not applicable

Competing interests

SJH is a member of the scientific advisory boards of Frequency Therapeutics, Psy Therapeutics, and Souvien Therapeutics, none of who were involved in this study. SDS is a scientific advisor for Outermost Therapeutics, Inc., which played no part in the present study. The other authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1. Additional Material and Methods.

Additional file 2: Table S1. List of antibodies used.

Additional file 3: Table S2. List of sgRNA or PCR primers used.

Additional file 4: Figure S1. Both TSC1-Null iPSC clones (A and B) are pluripotent. Immunostaining of TSC1-Null iPSC clones after introducing a second somatic mutation in the TSC1 gene. DAPI in blue, OCT4 in red. Scale bar=500μm. Figure S2. Rapamycin reversed TSC1 NPC cell size but not proliferation. A. Bright field images of TSC1-Het and Null NPCs showing a decreased cell size after 24h treatment with 500 nM of rapamycin. Scale bar=100μm; n=3. B. TSC1 NPCs (Corr-WT, Het and Null-clone B) were treated with increasing concentrations of rapamycin in a 6 point, 5-fold serial dilution series (0 - 10 μM) for 72h. Cell viabilities were assessed using CellTiter-Glo assays and plotted as % viability (relative to DMSO). Dose response curve data is presented as +SEM (3 replicates/dose). C. TSC1-Null NPCs were treated with DMSO (left pannel) or 100 nM of rapamycin for 24h (right pannel) and stained with propidium iodide for cell cycle analysis. Proliferating cells are represented in the S phase of the cell cycle (red).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Martin, P., Wagh, V., Reis, S.A. et al. TSC patient-derived isogenic neural progenitor cells reveal altered early neurodevelopmental phenotypes and rapamycin-induced MNK-eIF4E signaling. Molecular Autism 11, 2 (2020).

Download citation


  • Tuberous sclerosis complex
  • TSC1
  • mTORC1
  • Induced pluripotent stem cells
  • Neural progenitor cells
  • Early neurodevelopment
  • Disease modeling
  • CRISPR/Cas9
  • MEK-ERK1/2
  • MNK1/2-eIF4E