- Open Access
TSC patient-derived isogenic neural progenitor cells reveal altered early neurodevelopmental phenotypes and rapamycin-induced MNK-eIF4E signaling
Molecular Autism volume 11, Article number: 2 (2020)
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 . 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 . 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 . 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 . 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 . 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.
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.
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 . 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) (http://www.bioinformatics.babraham.ac.uk/projects/fastqc). 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’ . 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 (https://broadinstitute.github.io/picard/), RNASeQC , RSeQC , and SamTools , 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) , 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 , 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 org.Hs.eg.db (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 . 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  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).
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 . 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).
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 . 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 . 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).
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)  not only plays a role in cancer biology, but also in neurobiology by regulating 5’ cap-dependent translation of specific mRNAs in neuronal cells . 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).
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  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.
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 . 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 . 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 . 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.
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
Clustered regularly interspaced short palindromic repeats / CRISPR associated protein 9
Cytoplasmic fragile X protein-interacting protein 1
Eukaryotic translation initiation factor 4E
Extracellular signal-regulated kinase
Induced pluripotent stem cells
MAP kinase-interacting kinase
Mammalian/mechanistic target of rapamycin complex 1
Neural crest cells
Neural progenitor cells
Phospho-ribosomal protein S6
Phospho-ribosomal protein S6 kinase
Subependymal giant cell astrocytoma
Single-stranded oligo donor
Tuberous sclerosis complex
Caban C, Khan N, Hasbani DM, Crino PB. Genetics of tuberous sclerosis complex: implications for clinical practice. Appl Clin Genet. 2017;10:1–8.
Henske EP, Jozwiak S, Kingswood JC, Sampson JR, Thiele EA. Tuberous sclerosis complex. Nat Rev Dis Primers. 2016;2:16035.
Salussolia CL, Klonowska K, Kwiatkowski DJ, Sahin M. Genetic Etiologies, Diagnosis, and Treatment of Tuberous Sclerosis Complex. Annu Rev Genomics Hum Genet. 2019.
Han JM, Sahin M. TSC1/TSC2 signaling in the CNS. FEBS letters. 2011;585(7):973–80.
Huang J, Manning BD. The TSC1-TSC2 complex: a molecular switchboard controlling cell growth. Biochem J. 2008;412(2):179–90.
Kelleher RJ 3rd, Bear MF. The autistic neuron: troubled translation? Cell. 2008;135(3):401–6.
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.
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.
Julich K, Sahin M. Mechanism-based treatment in tuberous sclerosis complex. Pediatr Neurol. 2014;50(4):290–6.
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.
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.
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.
Anderson RH, Francis KR. Modeling rare diseases with induced pluripotent stem cell technology. Mol Cell Probes. 2018;40:52–9.
Connor B. Concise Review: The Use of Stem Cells for Understanding and Treating Huntington's Disease. Stem Cells. 2018;36:146–60.
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.
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.
Silva MC, Haggarty SJ. Human pluripotent stem cell-derived models and drug screening in CNS precision medicine. Ann N Y Acad Sci. 2019.
Musunuru K. Genome editing of human pluripotent stem cells to generate human cellular disease models. Dis Model Mech. 2013;6(4):896–904.
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.
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.
Blair JD, Hockemeyer D, Bateup HS. Genetically engineered human cortical spheroid models of tuberous sclerosis. Nat Med. 2018;24(10):1568–78.
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.
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.
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.
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.
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.
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.
Kaushik G, Zarbalis KS. Prenatal neurogenesis in autism spectrum disorders. Front Chem. 2016;4:12.
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.
Packer A. Neocortical neurogenesis and the etiology of autism spectrum disorder. Neurosci Biobehav Rev. 2016;64:185–95.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28(16):2184–5.
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.
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.
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.
Mandal PK, Rossi DJ. Reprogramming human fibroblasts to pluripotency using modified mRNA. Nat Protoc. 2013;8(3):568–82.
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.
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.
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.
Stead RL, Proud CG. Rapamycin enhances eIF4E phosphorylation by activating MAP kinase-interacting kinase 2a (Mnk2a). FEBS Lett. 2013;587(16):2623–8.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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).
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
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.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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).
About this article
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). https://doi.org/10.1186/s13229-019-0311-3
- Tuberous sclerosis complex
- Induced pluripotent stem cells
- Neural progenitor cells
- Early neurodevelopment
- Disease modeling