- Open Access
Differentiation of multipotent neural stem cells derived from Rett syndrome patients is biased toward the astrocytic lineage
Molecular Brainvolume 8, Article number: 31 (2015)
Rett syndrome (RTT) is one of the most prevalent neurodevelopmental disorders in females, caused by de novo mutations in the X-linked methyl CpG-binding protein 2 gene, MECP2. Although abnormal regulation of neuronal genes due to mutant MeCP2 is thought to induce autistic behavior and impaired development in RTT patients, precise cellular mechanisms underlying the aberrant neural progression remain unclear.
Two sets of isogenic pairs of either wild-type or mutant MECP2-expressing human induced pluripotent stem cell (hiPSC) lines were generated from a single pair of 10-year-old RTT-monozygotic (MZ) female twins. Mutant MeCP2-expressing hiPSC lines did not express detectable MeCP2 protein during any stage of differentiation. The lack of MeCP2 reflected altered gene expression patterns in differentiated neural cells rather than in undifferentiated hiPSCs, as assessed by microarray analysis. Furthermore, MeCP2 deficiency in the neural cell lineage increased astrocyte-specific differentiation from multipotent neural stem cells. Additionally, chromatin immunoprecipitation (ChIP) and bisulfite sequencing assays indicated that anomalous glial fibrillary acidic protein gene (GFAP) expression in the MeCP2-negative, differentiated neural cells resulted from the absence of MeCP2 binding to the GFAP gene.
An isogenic RTT-hiPSC model demonstrated that MeCP2 participates in the differentiation of neural cells. Moreover, MeCP2 deficiency triggers perturbation of astrocytic gene expression, yielding accelerated astrocyte formation from RTT-hiPSC-derived neural stem cells. These findings are likely to shed new light on astrocytic abnormalities in RTT, and suggest that astrocytes, which are required for neuronal homeostasis and function, might be a new target of RTT therapy.
Rett syndrome (RTT; MIM 312750) is a representative X-linked neurodevelopmental disease distinguished by repetitive and stereotypic hand movements replacing purposeful hand use, accompanied by gait ataxia, seizures, and autistic features. Patients present with apparently normal psychomotor development during the first 6–18 months of life, followed by a short period of developmental stagnation, and then rapid disease progression with the characteristic features described above . The incidence of RTT is 1 in 10,000–15,000 female births [2,3]. A linkage analysis demonstrated that the region harboring a causative gene is located between Xq27 and Xqter . Meanwhile, a mutation analysis identified an X-linked gene encoding methyl-CpG-binding protein 2 (MeCP2) in the Xq28 region as the cause of RTT, where > ~95% of patients with a classic RTT phenotype carry the MECP2 mutation [5,6].
MeCP2 is a transcription repressor that inhibits transcription by binding to methylated CpG dinucleotides, and also by recruiting co-repressors and chromatin remodeling proteins . Thus, mutant MeCP2 affects large-scale chromatin organization , resulting in the mal-regulation of a number of genes, including neural and synaptic genes [9-13]. Furthermore, because RTT is an X-linked dominant disorder, phenotypic differences between female RTT patients have generally been attributed to variances in X chromosome inactivation (XCI) patterns, with skewing in favor of the mutant allele for the more severe clinical phenotypes [14-16].
Nevertheless, it remains unknown how developmental defects occur in the RTT brain at the cellular level. Recently, human induced pluripotent stem cell (hiPSC) technology has facilitated the modeling of neurological diseases by permitting the reprogramming of somatic cells into pluripotent cells . So far, several studies have been performed with hiPSCs derived from patients with RTT and other neurological and neurodevelopmental diseases [18-26]. Previous reports of differentiated cells derived from RTT patient-specific hiPSCs demonstrated several abnormal in vitro phenotypes, such as diminished cell soma and nuclear sizes, reduced expression of neuronal markers, and attenuated dendritic spine density [19,21-24].
We recently reported a rare monozygotic (MZ) case of RTT in twins in which the genomic sequences were identical, including a MECP2 frame-shift mutation (G269AfsX288) . Interestingly, the patients (designated the RTT-MZ twins) showed divergent symptom severity regarding impaired neurological development, despite an identical genomic structure. Taking advantage of the nonrandom pattern of XCI in female hiPSCs  and the shared genetic background of the RTT-MZ twins , we aimed to generate two sets of isogenic pairs of wild-type MECP2- and mutant MECP2-expressing hiPSC lines from the RTT-MZ individuals. We then went on to show detailed defects in neural cells related to MeCP2 protein deficiency-induced developmental defects in RTT. Furthermore, we set out to clarify whether epigenetic and environmental cues participate in the abnormal neural development of RTT patients by comparison of the RTT-hiPSC lines derived from the RTT-MZ twin patients, who presented with varying symptom severity.
Mosaic expression patterns in fibroblasts procured from RTT-MZ twins
We recently established fibroblast cell lines from the RTT-MZ twins (patients RS1 and RS2). These RTT twins have a de novo frame-shift mutation in exon 4 (c.806delG) that truncates the MeCP2 protein within the transcriptional repression domain (Fig. 1A). We also reported that fibroblasts generated from both patients exhibited random XCI patterns , which were detected by the methylation-specific polymerase chain reaction (PCR)-based HUMARA (human androgen receptor) XCI assay . To examine the expression patterns of MeCP2 in RS1 and RS2 fibroblasts, immunostaining was performed with a specific primary antibody against MeCP2. Consequently, the fibroblast lines derived from the RTT-MZ twins included both MeCP2-positive and MeCP2-negative cells (Fig. 1B). Such mosaic expression patterns for the MeCP2 suggests that the fibroblasts comprise MeCP2-positive cells with the X chromosome harboring wild-type MECP2 as the active MeCP2 species, and MeCP2-negative cells with the X chromosome harboring mutant MECP2 as the active MeCP2 species. The fractions of MeCP2-positive cells among the RS1 and RS2 fibroblasts were 0.64 and 0.60, respectively (Fig. 1C).
Generation and characterization of RTT-MZ hiPSC lines
We utilized standard methods and transduction of OCT4-, SOX2-, KLF4-, and c-MYC-containing retroviruses to reprogram fibroblasts derived from the RTT-MZ twins into RTT-hiPSCs. Each RTT-hiPSC clone was clonally isolated and selected by morphological criteria and transgene silencing. We also verified pluripotency of the stem cells by immunostaining of the hiPSCs with primary antibodies against pluripotency markers (NANOG, OCT4, and TRA-1-81) (Fig. 2A and Additional file 1A).
To evaluate in vivo pluripotency, we injected the RTT-hiPSCs into the testes of immunodeficient mice, and confirmed the formation of teratomas containing derivatives of all three embryonic germ layers (Additional file 1B). No abnormalities were found in the karyotypes of any of the hiPSC lines (Additional file 1C). Notably, most of the selected hiPSC clones were either all MeCP2-positive or MeCP2-negative, and putatively originated from a single MeCP2-positive or MeCP2-negative fibroblast (Fig. 2A). Therefore, we isolated both wild-type MECP2-expressing (RS1-52 M and RS2-65 M) and mutant MECP2-expressing (RS1-61P and RS2-62P) clones from each patient (RS1 and RS2) to generate four hiPSC lines.
We next analyzed the XCI patterns in the hiPSC lines, and found that the maternal X chromosome was active in the MeCP2-positive (RS1-52 M and RS2-65 M) clones, while the paternal X chromosome was active in the MeCP2-negative (RS1-61P and RS2-62P) clones (Fig. 2B). Although we isolated several clones that partially include iPS cells with two active X chromosomes, we only used clones with an inactive X chromosome (i.e., the “standard” XCI status of undifferentiated hiPSCs ) in the present study.
Our previous study  revealed that the maternally-derived X chromosome carries the wild-type MECP2 allele, whereas the paternally-derived X chromosome carries the mutant MECP2 allele. These results were shown by sequencing the MECP2 gene in somatic hybrid cell clones carrying either the maternal or the paternal X chromosome of the RTT-MZ twins. Accordingly, the RS1-52 M and RS2-65 M hiPSC lines, in which maternal wild-type MECP2 was preferentially active, exhibited MeCP2 expression in the nuclei, whereas the RS1-61P and RS2-62P hiPSC lines, in which paternal mutant MECP2 was preferentially active, did not (Fig. 2A).
MeCP2 expression in neural cells differentiated from RTT-hiPSCs
We next differentiated the four RTT-hiPSC lines into neural cells and examined MeCP2 expression in the progeny. Immunostaining revealed that all of the cells derived from the RS1-52 M and RS2-65 M hiPSC lines (containing preferentially active maternal wild-type MECP2) expressed MeCP2 in the nucleus, like the parental hiPSCs. In particular, MAP2-positive neurons expressed MeCP2 more strongly than glial fibrillary acidic protein (GFAP)-positive astrocytes (Fig. 3A) as previously shown in mice [30,31]. However, no MeCP2 expression was found in the progeny of the RS1-61P and RS2-62P hiPSC lines (containing preferentially active paternal mutant MECP2), or in the nucleus of RS1-61P/RS2-62P hiPSC-derived MAP2-positive neurons or GFAP-positive astrocytes (Fig. 3A). Therefore, neural cells differentiated from the RTT-hiPSC lines apparently retain the XCI status of the undifferentiated cells.
Global gene expression in RTT-hiPSCs and differentiated neural cells
We next investigated whether the loss of MeCP2 protein affects the transcriptional network in undifferentiated hiPSCs and neural progeny differentiated from RTT-hiPSCs. Global gene expression levels in wild-type MECP2- and mutant MECP2-expressing cells were examined by comparative microarray analyses of undifferentiated hiPSCs and differentiated neural cells. The gene expression data were normalized and subjected to principal component analysis (PCA) and hierarchical clustering, with the exception of those genes with low expression in all samples.
The samples were clustered into hiPSC or neural cell (N) groups depending on the MECP2 expression pattern. In the PCA analysis, the original hiPSCs were clustered tightly into one region, signifying few differences from one hiPSC clone to the next. On the other hand, the neural cells were roughly clustered into two groups that reflected the status of the activated X chromosome (Fig. 3B and Additional file 2). Next, we constructed scatter plots to compare wild-type MECP2- and mutant MECP2-expressing cells. In the undifferentiated hiPSCs, the correlation coefficient (R2) was > 0.99 for any pair of hiPSC clones. However, in the differentiated neural cells, higher correlations (R2 > 0.98) were only found in comparisons of hiPSCs with the same XCI pattern (i.e., paternal X vs. paternal X or maternal X vs. maternal X), while lower correlations (R2 < 0.98) were only found in comparisons of hiPSCs with different XCI patterns (paternal X vs. maternal X) (Fig. 3C). The clustering dendrogram (Additional file 3A) also showed that global gene expressions for all four RTT-hiPSC clones during the undifferentiated stage. After this time, mutant MECP2-expressing hiPSCs became distinguishable from wild-type MECP2-expressing hiPSCs due to neural differentiation. Nevertheless, Tanaka et al.  recently characterized five patient-speific cell lines with different mutations to show that hiPSCs with MECP2 gene mutations are distinguishable from normal hiPSCs and ESCs even during the undifferentiated stage. By contrast, our isogenic cell lines might largely exclude the effect of donor divergence, permitting the elucidation of MECP2 function during neural differentiation. While Tanaka and colleagues  demonstrated that the gene expression of a key mitochondrial transcription factor, NR3C1, was upregulated in mutant MECP2-expressing cells after neuronal differentiation, we failed to detect such a difference in the current study (Additional file 3B). These results suggest that MeCP2 plays a more important role as a transcriptional regulator in differentiated neural cells than in undifferentiated hiPSCs.
Enhanced astrocytic differentiation of mutant MECP2-expressing neural stem cells derived from RTT-hiPSCs
Previous studies revealed that MeCP2 is involved in the regulation of astroglial gene expression [32-34]. Gfap and S100β are expressed at significantly higher levels in astrocytes derived from MeCP2-null mouse embryonic stem cells than in those derived from wild type mouse embryonic stem cells [35,36]. Moreover, a truncated form of RTT-associated MeCP2 (R168X) is reported to be unable to promote neuronal differentiation in mice, and instead promotes an abnormally high degree of astrocytic differentiation . We characterized neurospheres derived from the four RTT-hiPSC lines (RS1-61P, RS2-62P, RS1-52 M, and RS2-65 M), and found no significant differences in neurosphere number, size, or expression of neural stem cell markers between the two groups derived from mutant vs. wild-type MECP2-expressing hiPSCs (Additional file 4). Next, the neurospheres were differentiated into neural cells by using an adhesion culture method without fibroblast growth factor-2 (FGF-2) at ~30 days in vitro (Fig. 4A). The differentiated cells were classified by their expression of βIII-tubulin and GFAP. Immunocytochemical analysis revealed that the neural cells originating from MeCP2-negative hiPSC lines (RS1-61P and RS2-62P) contained significantly higher proportions of GFAP-positive cells than those originating from MeCP2-positive hiPSC lines (RS1-52 M and RS2-65 M) (Fig. 4B). Most of the GFAP-positive cells failed to express neuronal markers; therefore, we concluded that these cells were in fact astrocytes.
During the neural development of the mammalian central nervous system (CNS), neural stem cells initially differentiate into neurons, followed by astrocytes and oligodendrocytes at a later stage. We and others reported that MeCP2 binds to the promoter region of astrocyte-specific genes, including S100β and Gfap, to inhibit the conversion of neurons into astrocytes in mammals [33,34,37-39]. To determine how MeCP2-regulated astrocyte-specific genes are influenced by the absence of MeCP2, we next examined the expression of S100β and GFAP in MeCP2-negative and McCP2-positive neural cells. As a result, MeCP2-negative RS1-61P and RS2-62P neural cells showed significantly enhanced expression of the astrocyte-specific genes, GFAP and S100β (Fig. 4C). However, decreased expression levels of the neuronal genes, TUBB3 and MAP2, were observed in several samples (i.e., TUBB3 in RS2-62P cells and MAP2 in RS1-61 cells) (Fig. 4C). These observations suggest that MeCP2 deficiency in the neural cell lineage increases astrocytic differentiation from multipotent neural stem cells.
Lack of MeCP2 binding to the GFAP gene in RTT-hiPSC-derived neural stem cells
MeCP2 was previously reported to bind to the highly methylated exon 1 region of the Gfap gene in mouse neurons [33,38]. To confirm these results in wild-type or mutant MECP2-expressing neurospheres, we performed chromatin immunoprecipitation (ChIP) assays with a specific primary antibody against MeCP2. Binding data were quantitated via real-time PCR by using primers spanning from the STAT3 binding site within the GFAP promoter to exon 1, relative to the transcription start site of GFAP (Fig. 5A, indicated by blue bidirectional arrows). For both the STAT3 binding site and GFAP exon 1, precipitated genomic fragments were scarcely detectable by quantitative PCR (qPCR) in neurospheres derived from mutant MECP2-expressing hiPSC lines (RS1-61P, RS2-62P), while significant precipitation of genomic fragments was detected in neurospheres derived from wild-type MECP2-expressing hiPSC lines (RS1-52 M and RS2-65 M) (Fig. 5B).
Based on these results, we hypothesized that the increased expression of GFAP was induced by lack of direct MeCP2 binding to the GFAP gene in the mutant MECP2-expressing cells. However, it is possible that MeCP2 dysfunction indirectly enhanced the maturation of neural stem cells with truncated MeCP2, accompanied by decreased methylation of gfap , which might in turn modify MeCP2 protein recruitment. Thus, we performed bisulfite sequencing at the STAT3 binding site within the GFAP locus (Fig. 5A, indicated by black bidirectional arrows) by using wild-type or mutant MECP2-expressing neurospheres. Nevertheless, the MeCP2 binding sites in the GFAP gene were similarly hypermethylated in both wild-type and mutant MECP2-expressing neural stem cells differentiated from the four hiPSC lines (Fig. 5C). These findings signify that aberrant GFAP gene expression in neural stem cells containing truncated MeCP2 is not caused by changes in DNA methylation status (i.e., hypomethylation) within the genomic region encompassing the STAT3 binding site. Instead, the aberrant gene expression may be due to alterations in the amount of MeCP2 that binds to the STAT3 binding site within the GFAP promoter and the hGFAP exon1.
Female RTT patients display somatic mosaicism, because the MECP2 gene is located on the X chromosome and undergoes random XCI. Here, we successfully generated two kinds of isogenic RTT-hiPSC lines that expressed only the wild-type or the mutant MECP2 allele from the same RTT patient. So far, several groups have demonstrated disease modeling of RTT by using patient-specific hiPSCs [19,21-24]. However, some of these studies were conducted by comparisons between mutant hiPSCs derived from RTT patients and control hiPSCs derived from healthy subjects .
Regarding hiPSC technologies, clonal variations caused by donor differences and cell types of origin reportedly determine the differentiation potential of hiPSC clones. Hence, isogenic controls genetically engineered by using ZFN, TALEN or CRISPR-Cas9 systems are required for precise disease modeling with patient-specific hiPSCs . Our isogenic RTT-hiPSC system has two significant advantages over these previous systems. First, our system allows clonal generation of hiPSCs expressing either wild-type or mutant MeCP2 protein. Second, the generated hiPSCs are isogenic without genetic engineering. By utilizing our homogenous RTT-hiPSC populations, we were able to show that global gene expression patterns were distinctly different between wild-type MeCP2-cells and mutant MeCP2-expressing cells in terms of both neurospheres and differentiated neural stem cells. Although MeCP2 was weakly expressed in undifferentiated wild-type hiPSCs, significant differences in gene expression profiles were not observed among the four hiPSC lines. These results suggest that the MeCP2 primarily functions in cells already committed to a neural lineage.
We also found increased proportions of astrocytes among mutant MeCP2-expressing neural cells compared with wild-type MeCP2-expressing neural cells. Other investigators similarly documented MeCP2 expression in glial cells, including astrocytes , implying that glial MeCP2 might play a distinct role in the maintenance of neuronal functions, such as neuronal maturation and/or dendritic arborization . Although previous studies showed that a conditional MeCP2 knockout mouse (which lacked MeCP2 expression in mature neurons) displayed a RTT-like neurological phenotype, mutant astrocytes derived from this model mouse failed to support normal dendritic arborization of either wild-type or mutant neurons . Furthermore, mutant MECP2-expressing astrocytes differentiated from RTT-hiPSCs apparently had significant effects on wild-type neurons . Intriguingly, another study in mice showed that the MeCP2-deficient state can spread between astrocytes via gap junctions .
The current investigation did not examine dendritic phenotypes in neurons derived from RTT-hiPSCs. However, our model may be useful for mixed culture experiments (i.e., wild-type MECP2-expressing neurons and mutant MECP2-expressing astrocytes) to distinguish cell-autonomous from non-cell-autonomous effects in RTT patient-derived neural cells. Additionally, we demonstrated that astrocytic marker genes, particularly GFAP, were aberrantly expressed under conditions of MeCP2 deficiency due to the deregulation of transcription, and we employed ChIP and bisulfite sequencing analyses to show induction of GFAP expression in RTT-hiPSC-derived neural cells caused by the absence of direct binding of mutant MeCP2 to the GFAP gene. Hence, mutant MeCP2 may contribute to the excessive differentiation of GFAP-positive astrocytes.
Several groups reported that human female RTT patients present with microcephalic brains, as do male and female Rett syndrome model mice [44-46]. Moreover, analysis of postmortem brain tissue from RTT patients revealed increased GFAP protein expression in the RTT brain . Consistent with these results, we anticipate that the brains of RTT patients will exhibit a larger than normal percentage of astrocytes. Thus, astrocytes differentiated from RTT patient-derived, mutant MECP2-expressing hiPSC lines may be useful for drug development for RTT therapy, because some drugs for mental disorders (e.g., antidepressants) affect not only neurons, but also activate astrocytes. The activated astrocytes then go on to carry out specific functions that result in the reactivation of cortical plasticity, potentially leading to the readjustment of abnormal neural networks [47,48].
This study aimed to clarify the involvement of epigenetic and environmental cues in the neural development of RTT-MZ twins (patients RS1 and RS2) with differing symptom severity. Although we successfully generated hiPSCs from both patients, no significant differences in mutant MeCP2 phenotypes were observed between RS1 and RS2 cells. These results suggest that epigenetic and environmental cues might work together to affect the neural phenotypes of RTT patients. Another possibility is that skewed XCI occurs in the CNS in RTT patients. Although we previously reported that biased XCI does not take place in the other somatic cells (i.e., skin, blood, and hair) of these patients , biased XCI in the CNS may underlie the discrepancy in symptom severity between isogenic RTT-MZ twin patients.
In summary, we successfully generated two sets of isogenic RTT-hiPSC lines that expressed only the wild-type or the mutant MECP2 allele from the same RTT patient, and revealed that abnormal astrocytic differentiation is involved in RTT pathogenesis. Specifically, we used hiPSC technology to show that (1) MeCP2 protein deficiency results in dysregulation of GFAP expression in mutant MECP2-expressing neural cells, and (2) an increased number of astrocytes are differentiated from RTT-hiPSC-derived neural stem cells carrying mutant vs. wild-type MECP2.
Materials and methods
Isolation of skin fibroblasts and generation of hiPSCs
A skin-punch biopsy from a single pair of 10-year-old Japanese MZ twins with RTT (patients RS1 and RS2)  was performed with written informed consent from their parents. The obtained fibroblasts were used to establish the four hiPSC lines discussed above. The study protocol was reviewed and approved by the Research Ethics Committees of Keio University (Approval No. 20080016) and the University of Yamanashi (Approval Nos. 523 and 699). This study was conducted in accordance with the principles expressed in the Declaration of Helsinki.
A standard retrovirus method was used to establish the hiPSCs. The maintenance of fibroblasts, lentiviral production, retroviral production, infection, stem cell culture, and characterization of hiPSCs was conducted as described previously .
Sample collection and DNA isolation
All peripheral blood, skin fibroblasts, and hiPSC samples used in this study were obtained from RTT-MZ patients whose parents gave informed consent for medical research at the University of Yamanashi and Keio University Hospital. DNA from lymphocytes, skin fibroblasts, and hiPSCs was extracted with the DNeasy Blood and Tissue Kit (QIAGEN) according to the manufacturer’s instructions.
XCI patterns were obtained in the extracted DNA by using a methylation-specific PCR-based method. Bisulfite treatment of genomic DNA followed by PCR was performed as previously described . Briefly, DNA was treated with sodium bisulfite, and the treated DNA was amplified via PCR by using two primer sets (one for the methylated inactive X chromosome, and the other for the unmethylated active X chromosome). The primer sets were designed within the CpG island (exon 1) of the HUMARA gene located at Xq11-a12. The PCR products were mixed with a size standard and separated on an ABI 310 DNA Sequencer equipped with Genescan Software (Applied Biosystems).
Teratoma formation was assessed by injecting hiPSCs into the testes of 8-week-old NOD.CB17-Prkdc scid/J mice (OYC International, Inc.), as previously described . Eight weeks after hiPSC transplantation, tumors were dissected and fixed with phosphate buffered saline (PBS) containing 4% formaldehyde. Paraffin-embedded tissue was sectioned and stained with hematoxylin & eosin. Images were obtained by using a BZ-9000 microscope (Keyence).
Karyotyping of hiPSCs
Standard G-banding analysis was performed in the four hiPSC lines (RS1-52 M, RS1-61P, RS2-65 M, and RS2-62P) to rule out the possibility of abnormal karyotypes that can occur during the generation of hiPSCs.
In vitro differentiation of hiPSCs
Neural differentiation of hiPSCs was performed as described previously , except that epidermal growth factor and Noggin were omitted from the culture system (Matsumoto et al., manuscript in submission). The hiPSCs were plated onto the bottom of T75 flasks (Nunclon) and maintained for 14 days. Neurospheres were repeatedly passaged by dissociation into single cells, followed by culture in the same manner. Neurospheres were differentiated into neural cells by using an adhesion culture method without FGF-2, and then subjected to qPCR and immunostaining analysis. Typically, neurospheres between passages 3 and 10 were used for the analysis. RNA samples were collected on day 14 of neurosphere culture from an arbitrary number of cells. For terminal differentiation, neurospheres were cultured in 1 × B27 medium, dissociated neural precursor cells allowed to adhere to poly-L-ornithine-/fibronectin-coated coverslips and then cultured for 30–32 days. RNA was then isolated from terminally differentiated neurospheres. Alternatively, the differentiated cells were immunostained with primary antibodies against βIII-tubulin and GFAP. Approximately 1 × 105 neural stem cells were plated onto 14-mm glass coverslips for immunostaining of differentiated neural cells, and ~ 1.2 × 106 neural stem cells were plated onto a well using 6-well plate for RNA isolation from differentiated neural cells. For long-term neural differentiation, mouse astrocytes obtained from E16.5 ICR mice were utilized with a cell culture insert system (Falcon® product #353102). RNA samples were obtained from iPSCs, neurospheres, and neuronal cells at passages 12–20, passages 3–9 (day 14), and day 30–32, respectively.
Immunocytochemical analysis of hiPSCs and neural cells differentiated from hiPSCs
Cells were fixed with PBS containing 4% paraformaldehyde for 30 min at room temperature and incubated with primary antibodies against the following proteins: MeCP2 (1:200, Cell Signaling Technology), phalloidin (1:2000, Dyomics), NANOG (1:500, CosmoBio), OCT4 (1:200, Santa Cruz Biotechnology, Inc.), TRA-1-60 (1:1000, Millipore), TRA-1-81 (1:1000, Millipore), βIII-tubulin (1:2000, Sigma Chemical Co.), MAP2 (1:250, Sigma Chemical Co.), and GFAP (1:1000, Invitrogen). They were then washed with PBS and incubated with an Alexa Fluor 488-, Alexa Fluor 555-, or Alexa Fluor 647-conjugated secondary antibody (1:500, Invitrogen), as appropriate. Images were obtained using an Axioplan 2 microscope (Carl Zeiss). The number of GFAP-positive cells was counted among 100 Hoechst-positive cells for each experiment (n = 5).
For microarray analysis, RNA quality was assessed by using a 2100 Bioanalyzer (Agilent Technologies). Total RNA (100 ng) was reverse transcribed, biotin- labeled, and hybridized to a GeneChip® Human Genome U133 plus 2.0 Array (Affymetrix). The array was subsequently washed and stained in a Fluidics Station 450 (Affymetrix), as instructed by the manufacturer. The microarrays were scanned by using a GeneChip® Scanner 3000 (Affymetrix), and the raw image files were converted into normalized signal intensity values by using the MAS 5.0 algorithm.
Next, targets were selected that (a) were “present” in at least one of the eight arrays analyzed, and (b) had a probe intensity of ≥ 50. In total, 23,463 targets were identified from this initial screen. The normalized logs were hierarchically clustered based on uncentered correlation with complete linkage by using Cluster 3.0 Software . The PCA plot was generated by using TIGR MeV (Multiple Experimental Viewer) Software . The three-dimensional image of the PCA plot and the scatter plot were generated by using GeneSpring GX Software (ver. 12.6.1).
Quantitative reverse transcription PCR (qPCR) assay
Total cellular RNA was extracted by using the TRIzol Reagent® (Life Technologies), the RNase-Free DNase Set (QIAGEN), and the RNeasy® Mini Kit (QIAGEN). Next, cDNA synthesis was performed by using the SuperScript III First-Strand Synthesis System for RT-PCR (Life Technologies) with oligo-dT primers according to the manufacturer’s guidelines. Real-time RT-PCR was then performed on a ABI PRISM® 7900HT Sequence Detection System (Applied BioSystems) by using SYBR Premix ExTaq Tli RNaseH Plus (Takara). The qPCR amplification was performed by using the following primers: GFAP-forward (5’-TGTGAGGCAGAAGCTCCAGGATGA-3’) and GFAP-reverse (5’- AGGGTGGCTTCATCTGCTTCCTGT-3’); S100β-forward (5’-GTGGCCCTCATCGACGTTTT-3’) and S100β-reverse (5’-ACCTCCTGCTCTTTGATTTCCTCT-3’); TUBB3-forward (5’-ATTTCATCTTTGGTCAGAGTGGGGC-3’) and TUBB3-reverse (5’-TGCAGGCAGTCGCAGTTTTCAC-3’); MAP2-forward (5’-GGCCCAAGCTAAAGTTGGTTCTC-3’) and MAP2-reverse (5’-GCAGTGACATCCTCAGCCAAAG-3’); and ACTB-forward (5’-TGAAGTGTGACGTGGACATC-3’) and ACTB-reverse (5’-GGAGGAGCAATGATCTTGAT-3’). Relative gene expression levels of GFAP, S100β, TUBB3, and MAP2 were normalized to ACTB expression and standardized by the RS1-52 M value, which was set to 1.
Nuclei were obtained from neurospheres derived from each hiPSC line. ChIP was performed with antibodies against MeCP2 (Cell Signaling Technology), as described previously , with the following modifications. Neurospheres were crosslinked with 1% formaldehyde for 10 min to yield crosslinked chromatin. They were then incubated with glycine at a final concentration of 200 mM for another 5 min, and stored at −80°C until use. The lysed nuclear pellets were sonicated six times with a 30-sec ON, 60-sec OFF cycle by using a Bioruptor® Sonicator (Diagenode, Inc.). The crosslinked chromatin was subsequently eluted from magnetic beads (Dynabeads® M-280 Sheep anti-Rabbit IgG; Life Technologies). Co-immunoprecipitated DNA was detected via qPCR by using the following primers: STAT3 binding site-forward (5’-TCATGCCCAGTGAATGACTC-3’) and STAT3 binding site-reverse (5’-AGATGCCAGGCTGTCAGG-3’); and hGFAP exon1-forward (5’-AGAGCAGGATGGAGAGGAGA-3’) and hGFAP exon1-reverse (5’- CCTTGAAGCCAGCATTGAGT-3’).
Bisulfite sequencing analysis
Sodium bisulfite treatment of genomic DNA was performed by using the Methylamp DNA Modification Kit (Epigentek) according to the manufacturer’s instructions. A region in the GFAP promoter region containing the STAT3 binding site was amplified from bisulfite-treated genomic DNA via PCR with the following forward and reverse primers: hGFAP promoter-forward (5’-TTGGGGAGGAGGTAGATAGTTAGGTTTT-3’) and hGFAP promoter-reverse (5’-CATCCCCTAATCCCCTTTCCTAAA-3’). PCR products were cloned into the pT7Blue vector (Novagen), and at least 12 randomly selected clones were sequenced.
All quantifiable data are expressed as the means ± the standard deviation of the mean. The statistical significance of differences between conditions was analyzed by using Student’s t-test and Welch’s t-test. In all cases, p values of < 0.05 were considered statistically significant.
Fabio RA, Giannatiempo S, Antonietti A, Budden S. The role of stereotypies in overselectivity process in Rett syndrome. Res Dev Disabil. 2009;30:136–45.
Hagberg B, Aicardi J, Dias K, Ramos O. A progressive syndrome of autism, dementia, ataxia, and loss of purposeful hand use in girls: Rett's syndrome: report of 35 cases. Ann Neurol. 1983;14:471–9.
Bienvenu T, Philippe C, De Roux N, Raynaud M, Bonnefond JP, Pasquier L, et al. The incidence of Rett syndrome in France. Pediatr Neurol. 2006;34:372–5.
Webb T, Clarke A, Hanefeld F, Pereira JL, Rosenbloom L, Woods CG. Linkage analysis in Rett syndrome families suggests that there may be a critical region at Xq28. J Med Genet. 1998;35:997–1003.
Amir RE, Van den Veyver IB, Wan M, Tran CQ, Francke U, Zoghbi HY. Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. Nat Genet. 1999;23:185–8.
Monteggia LM, Kavalali ET. Rett syndrome and the impact of MeCP2 associated transcriptional mechanisms on neurotransmission. Biol Psychiatry. 2009;65:204–10.
Ng HH, Zhang Y, Hendrich B, Johnson CA, Turner BM, Erdjument-Bromage H, et al. MBD2 is a transcriptional repressor belonging to the MeCP1 histone deacetylase complex. Nat Genet. 1999;23:58–61.
Agarwal N, Becker A, Jost KL, Haase S, Thakur BK, Brero A, et al. MeCP2 Rett mutations affect large scale chromatin organization. Hum Mol Genet. 2011;20:4187–95.
Chen WG, Chang Q, Lin Y, Meissner A, West AE, Griffith EC, et al. Derepression of BDNF transcription involves calcium-dependent phosphorylation of MeCP2. Science. 2003;302:885–9.
Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, et al. DNA methylation-related chromatin remodeling in activity-dependent BDNF gene regulation. Science. 2003;302:890–3.
Horike S, Cai S, Miyano M, Cheng JF, Kohwi-Shigematsu T. Loss of silent-chromatin looping and impaired imprinting of DLX5 in Rett syndrome. Nat Genet. 2005;37:31–40.
Chahrour M, Jung SY, Shaw C, Zhou X, Wong ST, Qin J, et al. MeCP2, a key contributor to neurological disease, activates and represses transcription. Science. 2008;320:1224–9.
Miyake K, Hirasawa T, Soutome M, Itoh M, Goto Y, Endoh K, et al. The protocadherins, PCDHB1 and PCDH7, are regulated by MeCP2 in neuronal cells and brain tissues: implication for pathogenesis of Rett syndrome. BMC Neurosci. 2011;12:81.
Amir RE, Van den Veyver IB, Schultz R, Malicki DM, Tran CQ, Dahle EJ, et al. Influence of mutation type and X chromosome inactivation on Rett syndrome phenotypes. Ann Neurol. 2000;47:670–9.
Young JI, Zoghbi HY. X-chromosome inactivation patterns are unbalanced and affect the phenotypic outcome in a mouse model of rett syndrome. Am J Hum Genet. 2004;74:511–20.
Gibson JH, Williamson SL, Arbuckle S, Christodoulou J. X chromosome inactivation patterns in brain in Rett syndrome: implications for the disease phenotype. Brain Dev. 2005;27:266–70.
Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131:861–72.
Hotta A, Cheung AY, Farra N, Vijayaragavan K, Seguin CA, Draper JS, et al. Isolation of human iPS cells using EOS lentiviral vectors to select for pluripotency. Nat Methods. 2009;6:370–6.
Marchetto MC, Carromeu C, Acab A, Yu D, Yeo GW, Mu Y, et al. A model for neural development and treatment of Rett syndrome using human induced pluripotent stem cells. Cell. 2010;143:527–39.
Muotri AR, Marchetto MC, Coufal NG, Oefner R, Yeo G, Nakashima K, et al. L1 retrotransposition in neurons is modulated by MeCP2. Nature. 2010;468:443–6.
Cheung AY, Horvath LM, Grafodatskaya D, Pasceri P, Weksberg R, Hotta A, et al. Isolation of MECP2-null Rett Syndrome patient hiPS cells and isogenic controls through X-chromosome inactivation. Hum Mol Genet. 2011;20:2103–15.
Kim KY, Hysolli E, Park IH. Neuronal maturation defect in induced pluripotent stem cells from patients with Rett syndrome. Proc Natl Acad Sci U S A. 2011;108:14169–74.
Ananiev G, Williams EC, Li H, Chang Q. Isogenic pairs of wild type and mutant induced pluripotent stem cell (iPSC) lines from Rett syndrome patients as in vitro disease model. PLoS One. 2011;6:e25255.
Cheung AY, Horvath LM, Carrel L, Ellis J. X-chromosome inactivation in rett syndrome human induced pluripotent stem cells. Front Psychiatry. 2012;3:24.
Tanaka Y, Kim KY, Zhong M, Pan X, Weissman SM, Park IH. Transcriptional regulation in pluripotent stem cells by methyl CpG-binding protein 2 (MeCP2). Hum Mol Genet. 2014;23:1045–55.
Williams EC, Zhong X, Mohamed A, Li R, Liu Y, Dong Q, Ananiev GE, Mok JC, Lin BR, Lu J, Chiao C, Cherney R, Li H, Zhang SC, Chang Q. Mutant astrocytes differentiated from Rett syndrome patients-specific iPSCs have adverse effects on wild-type neurons. Hum Mol Genet. 2014;23:2968-80.
Miyake K, Yang C, Minakuchi Y, Ohori K, Soutome M, Hirasawa T, et al. Comparison of Genomic and Epigenomic Expression in Monozygotic Twins Discordant for Rett Syndrome. PLoS One. 2013;8:e66729.
Kubota T, Nonoyama S, Tonoki H, Masuno M, Imaizumi K, Kojima M, et al. A new assay for the analysis of X-chromosome inactivation based on methylation-specific PCR. Hum Genet. 1999;104:49–55.
Tchieu J, Kuoy E, Chin MH, Trinh H, Patterson M, Sherman SP, et al. Female human iPSCs retain an inactive X chromosome. Cell Stem Cell. 2010;7:329–42.
Ballas N, Lioy DT, Grunseich C, Mandel G. Non-cell autonomous influence of MeCP2-deficient glia on neuronal dendritic morphology. Nat Neurosci. 2009;12:311–7.
Skene PJ, Illingworth RS, Webb S, Kerr AR, James KD, Turner DJ, et al. Neuronal MeCP2 is expressed at near histone-octamer levels and globally alters the chromatin state. Mol Cell. 2010;37:457–68.
Colantuoni C, Jeon OH, Hyder K, Chenchik A, Khimani AH, Narayanan V, et al. Gene expression profiling in postmortem Rett Syndrome brain: differential gene expression and patient classification. Neurobiol Dis. 2001;8:847–65.
Setoguchi H, Namihira M, Kohyama J, Asano H, Sanosaka T, Nakashima K. Methyl-CpG binding proteins are involved in restricting differentiation plasticity in neurons. J Neurosci Res. 2006;84:969–79.
Kohyama J, Kojima T, Takatsuka E, Yamashita T, Namiki J, Hsieh J, et al. Epigenetic regulation of neural cell differentiation plasticity in the adult mammalian brain. Proc Natl Acad Sci U S A. 2008;105:18012–7.
Okabe Y, Kusaga A, Takahashi T, Mitsumasu C, Murai Y, Tanaka E, et al. Neural development of methyl-CpG-binding protein 2 null embryonic stem cells: a system for studying Rett syndrome. Brain Res. 2010;1360:17–27.
Okabe Y, Takahashi T, Mitsumasu C, Kosai K, Tanaka E, Matsuishi T. Alterations of gene expression and glutamate clearance in astrocytes derived from an MeCP2-null mouse model of Rett syndrome. PLoS One. 2012;7:e35354.
Tsujimura K, Abematsu M, Kohyama J, Namihira M, Nakashima K. Neuronal differentiation of neural precursor cells is promoted by the methyl-CpG-binding protein MeCP2. Exp Neurol. 2009;219:104–11.
Namihira M, Nakashima K, Taga T. Developmental stage dependent regulation of DNA methylation and chromatin modification in a immature astrocyte specific gene promoter. FEBS Lett. 2004;572:184–8.
Forbes-Lorman RM, Kurian JR, Auger AP. MeCP2 regulates GFAP expression within the developing brain. Brain Res. 2014;1543:151–8.
Ran FA, Hsu PD, Wright J, Agarwala V, Scott DA, Zhang F. Genome engineering using the CRISPR-Cas9 system. Nat Protoc. 2013;8:2281–308.
Nagai K, Miyake K, Kubota T. A transcriptional repressor MeCP2 causing Rett syndrome is expressed in embryonic non-neuronal cells and controls their growth. Brain Res Dev Brain Res. 2005;157:103–6.
Kishi N, Macklis JD. MeCP2 functions largely cell-autonomously, but also non-cell-autonomously, in neuronal maturation and dendritic arborization of cortical pyramidal neurons. Exp Neurol. 2010;222:51–8.
Maezawa I, Swanberg S, Harvey D, LaSalle JM, Jin LW. Rett syndrome astrocytes are abnormal and spread MeCP2 deficiency through gap junctions. J Neurosci. 2009;29:5051–61.
Chen RZ, Akbarian S, Tudor M, Jaenisch R. Deficiency of methyl-CpG binding protein-2 in CNS neurons results in a Rett-like phenotype in mice. Nat Genet. 2001;27:327–31.
Duncan Armstrong D. Neuropathology of Rett Syndrome. J Child Neurol. 2005;20:747–53.
Chahrour M, Zoghbi HY. The story of Rett syndrome: from clinic to neurobiology. Neuron. 2007;56:422–37.
Sanacora G, Banasr M. From pathophysiology to novel antidepressant drugs: glial contributions to the pathology and treatment of mood disorders. Biol Psychiatry. 2013;73:1172–9.
Czeh B, Di Benedetto B. Antidepressants act directly on astrocytes: evidences and functional consequences. Eur Neuropsychopharmacol. 2013;23:171–85.
Ohta S, Imaizumi Y, Okada Y, Akamatsu W, Kuwahara R, Ohyama M, et al. Generation of human melanocytes from induced pluripotent stem cells. PLoS One. 2011;6:e16182.
Chaddah R, Arntfield M, Runciman S, Clarke L, van der Kooy D. Clonal neural stem cells from human embryonic stem cell colonies. J Neurosci. 2012;32:7771–81.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95:14863–8.
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34:374–8.
Kimura H, Hayashi-Takanaka Y, Goto Y, Takizawa N, Nozaki N. The organization of histone H3 modifications as revealed by a panel of specific monoclonal antibodies. Cell Struct Funct. 2008;33:61–73.
We thank the RTT patients and their parents for their cooperation in this study. We also thank J. Kohyama for technical advice on the ChIP assay; M. Isoda, H. Ebise, and H. Nakagawa (Sumitomo Dainippon Pharma Co. Ltd., Osaka, Japan) for help with the microarray analyses; and Y. Imaizumi for expert technical assistance. We also thank all the members of the Okano laboratory for helpful comments and discussions. The research described in this study was partially supported by grants from the New Energy and Industrial Technology Development Organization and the Ministry of Education, Science, Sports and Culture (MEXT) and Ministry of Helath, Labour and Welfare (MHLW) of Japan to H.O.; by the Program for Intractable Disease Research Utilizing Disease-specific iPS Cells funded by the Japan Science and Technology Agency (JST)/Japan Agency for Medical Research and Development (A-MED) to H.O.; by Grants-in-Aid for Scientific Research (KAKENHI); and by the Ministry of Economy, Trade and Industry (METI) of Japan for “Development of Core Technologies for Innovative Drug Development Based Upon IT” in the project focused on developing key technology for discovering and manufacturing drugs for next-generation treatment and diagnosis (the biological-verifying studies) to T.K.
T.A.-N., W.A., Ta.K., and H.O. conceived and designed the experiments and wrote the manuscript. T.A.-N. performed most of the experiments, analyzed the data, and wrote the manuscript. T.A.-N., T.M., Y.O., Te.K., M.O., and H.K. generated the patient-derived hiPSCs and analyzed the culture assay results. R.Y. analyzed the microarray data. T.S. and K.N. performed the bisulfite assay and helped with data analysis. K.M. and T.K. contributed to clinical and genetic analyses of the patient-coordinated study. All the authors read and approved the final version of the manuscript.
H.O. is a scientific consultant for SanBio Co., Ltd. M.S. and A.N. are employed by Takeda Pharmaceutical Co., Ltd. The authors declare that they have no competing interest.
Below is the link to the electronic supplementary material.
Characterization of RTT-hiPSCs. (A) Isogenic RTT-hiPSCs demonstrate similar embryonic stem cell-like morphology and stain positively for the pluripotency markers, NANOG, OCT4, and TRA-1-81. Scale bar for phase contrast images, 500 μm; scale bar for NANOG, OCT4, and TRA-1-81 immunostaining, 50 μm. (B) Representative images of teratomas generated in immunodeficient mice that received an intratesticular injection of RTT-hiPSCs. The teratomas corresponded to well-defined, cystic tumors containing tissues of all three germ layers (endoderm, mesoderm, and ectoderm). Scale bar, 100 μm. (C) Images of RTT-hiPSC karyotypes. (PDF 246 kb)
Three-dimensional image of the PCA. The clustering pattern of the cells (hiPSCs vs. neural cells (N) was dependent on the MECP2 expression pattern and the presence or absence of MeCP2 protein.
Gene expression analysis of iPSCs. (A) Global gene expression/comparative microarray analyses of undifferentiated hiPSCs and differentiated neural cells. N denotes neural cells. (B) Comparison of NR3C1 gene expression in hiPSCs and differentiated neural cells by qPCR. Relative gene expression was normalized by ACTB gene expression. Primer sets are listed in Additional file 5. Data were analyzed by Student’s t-test and Welch’s t-test (*p < 0.05).
Properties of neural stem cells derived from RTT-hiPSCs. (A) Diameter (n = 3) and (B) number of neurospheres (n = 16) derived from RTT-hiPSCs at 7 days in vitro. (C) The mRNA expression levels of the neural stem cell markers, NESTIN (NES) and SOX1, were quantified in RTT hiPSC-derived neurospheres by qPCR. Each value was normalized by comparison with ACTB expression and standardized by the value for the wild-type MECP2-expressing RS1-52 M clone, which was set to 1 (n = 4). Primer sets are listed in Additional file 5. Data were analyzed by Student’s t-test and Welch’s t-test (*p < 0.05).
List of primers used in Additional file 3 and 4.