<title>Introduction</title> <p>DNA methylation usually occurs in mammalian cells at CpG dinucleotides and approximately 60–90% of cytosines at these sites are methylated ##REF##6585800##[1]##. Most CpG-rich DNA fragments, or CpG islands in genes, have been thought to remain unmethylated even in cell types that do not express the genes ##REF##11782440##[2]##. However, changes in DNA methylation have been sporadically observed in CpG islands during development and are thought to play important roles in the regulation of cell type-specific gene expression ##REF##9581285##[3]##, ##REF##15728362##[4]##. DNA methylation also participates in the regulation of differentiation and embryonic development ##REF##1606615##[5]##. For example, inactivation of <italic>Oct3/4</italic> (<italic>Pou5f1</italic>) and <italic>Nanog</italic> genes by DNA methylation is important for early development ##REF##14761969##[6]##, ##REF##17352742##[7]##. These genes are essential for maintaining pluripotency of embryonic stem (ES) cells and early embryos ##REF##10742100##[8]##, ##REF##12787504##[9]##, ##REF##12787505##[10]## and are also known as two of four genes which have been shown to reprogram somatic cells to pluripotent stem (iPS) cells with the essential characteristics of embryonic stem (ES) cells ##REF##18029452##[11]##.</p> <p>Epigenetic modification is also thought to play an important role in altering the differentiation potential of neural precursor cells (NPCs). NPCs differentiate only into neurons in the midgestational brain, while they become able to generate astrocytes as gestation proceeds ##REF##11086984##[12]##. The Janus kinase (JAK) signal transducer and activator of the transcription (STAT) pathway, which is activated by cytokines, including leukemia inhibitory factor (LIF), can effectively induce astrocyte differentiation ##REF##9334309##[13]##, ##REF##9570793##[14]##. A particular cytosine residue within a STAT3-binding site in the astrocyte-specific marker glial fibrillary acidic protein (<italic>Gfap</italic>) gene promoter is highly methylated in NPCs midgestation (E11.5) when astrogenesis does not normally occur, while it becomes demethylated in late-stage (E14.5) NPCs that are prone to differentiating into astrocytes ##REF##11740937##[15]##. Although these late-stage (E14.5) NPCs with demethylated promoter do not express <italic>Gfap</italic>, the expression can be upregulated upon differentiation by leukemia inhibitory factor (LIF), which induces STAT3-activating cytokines ##REF##9334309##[13]##, ##REF##9570793##[14]##.</p> <p>Previously, we developed a genome-wide DNA methylation analysis called Microarray-based Integrated Analysis of Methylation by Isoschizomers (MIAMI) using a microarray with genome-wide probes and used for several applications ##REF##16407832##[16]##, ##REF##18200438##[17]##. With this method, we detected DNA methylation using the methylation-sensitive restriction enzyme <italic>Hpa</italic> II and its methylation-insensitive isoschizomer <italic>Msp</italic> I.</p> <p>Although the demethylation of <italic>Gfap</italic> is known in late-stage NPCs prior to astrocyte differentiation ##REF##11740937##[15]##, it is not known whether other genes confer this astrocytic property to cells. To study the role of DNA methylation in altering the differentiation potential of NPCs, here we analyzed DNA methylation changes in mouse NPCs between the mid-(E11.5) and late (E14.5) stages of development by a genome-wide DNA methylation profiling method using microarrays. We also compared the methylation status with that of postnatal day 1 (P1) astrocytes.</p>
<title>Materials and Methods</title> <title>Methylation profiling by MIAMI</title> <p>The MIAMI method was performed using one microgram of genomic DNA as previously described ##REF##16407832##[16]##, ##REF##18200438##[17]##. The complete experimental procedure can be obtained at <ext-link ext-link-type="uri" xlink:href="http://grc.dept.med.gunma-u.ac.jp/gene/image/MIAMI20Protocol20V4.pdf">http://grc.dept.med.gunma-u.ac.jp/gene/image/MIAMI20Protocol20V4.pdf</ext-link>. The microarray was the same as that used in our previous paper ##REF##18200438##[17]##. A brief description of the MIAMI method is as follows.</p> <p>Changes in methylation were judged by using the difference in methylation-sensitive <italic>Hpa</italic> II cleavage and methylation-insensitive <italic>Msp</italic> I cleavage between samples. To detect differences in methylation-sensitive <italic>Hpa</italic> II cleavage, 0.5 microgram of genomic DNA was digested with 40 units of <italic>Hpa</italic> II overnight in a 100-microliter volume containing 33 mM Tris-acetate (pH7.9), 66 mM KOAc, 10 mM MgOAc<sub>2</sub>, 0.5 mM DTT, and 0.01% BSA. The adaptor was prepared by annealing two oligonucleotides, <named-content content-type="gene">AGCACTCTCCAGCCTCTCACCGAG</named-content> and <named-content content-type="gene">CGCTCGGTGA</named-content>. After phenol extraction and ethanol precipitation, DNA was ligated to the adaptor with 60 units of E coli DNA Ligase (Takara, Japan). The first PCR was performed using 0.1 microgram of each ligation mix as a template with a primer <named-content content-type="gene">AGCACTCTCCAGCCTCTCACCGAG</named-content> using GeneTaq DNA polymerase (Nippon Gene, Japan). The reaction mixture was incubated for 5 min at 72°C and 3min at 94°C and subjected to 5 cycles of amplification consisting of 10 sec denaturation at 94°C, 30 sec annealing at 70°C, 2.5 min extension at 72°C. The final extension was lengthened to 9.5 min. Amplified DNA was digested with <italic>Msp</italic> I by adding 35 units of the enzyme directly to the solution. After 3-hour incubation, the digested DNA was subjected to further amplification by a second PCR. The condition was the same as in the first PCR except for the number of cycles. The optimal sub-saturation cycle number was adopted (usually 10 to 13 cycles). The PCR product was purified with a MiniElute PCR purification Kit (Qiagen, USA). To detect of differences in methylation-insensitive <italic>Msp</italic> I cleavage, 0.5 microgram of genomic DNA was digested with 40 units of <italic>Msp</italic> I overnight and subjected to the same procedure as for the amplification of unmethylated DNA fragments.</p> <p>Amplified <italic>Hpa</italic> II-cleaved DNA fragments from two samples to be compared were labeled with Cy3 and Cy5, respectively, as described ##REF##12181645##[37]##, and cohybridized to a microarray. All hybridization procedures except for washing were according to the manual of Agilent Technology (USA). Amplified <italic>Msp</italic> I-cleaved DNA fragments from two samples to be compared were labeled with Cy3 and Cy5, respectively, and cohybridized to another microarray with the same probes.</p> <p>We judged spots with methylation changes using two models; one in which spots with methylation changes had different values for <italic>Hpa</italic> II and <italic>Msp</italic> I cleavage, and the other in which spots with methylation changes had a large <italic>Hpa</italic> II cleavage difference (more than 5) and a small <italic>Msp</italic> I cleavage difference (less than 2). The reproducibility of the experiment was analyzed with 885 triplicated probes on each microarray and repeated experiments.</p> <title>Bisulfite genomic analysis</title> <p>Bisulfite treatment of genomic DNA was performed using a CpGenome DNA modification kit (INTERGEN). Modified DNA was amplified with the primers described in ##SUPPL##2##Table S3##. For COBRA analysis, the restriction enzymes are described in ##SUPPL##2##Table S3##.</p> <title>Expression microarray analysis</title> <p>Expression microarray analysis was performed using the Agilent mouse whole genome array and the procedure provided by Agilent Technologies. A signal ratio of more than 2 with a P-value of less than 0.01 was judged as upregulated. A signal ratio of less than 0.5 with a P-value of less than 0.01 was judged as downregulated.</p> <title>Animals and cell preparation</title> <p>Time-pregnant ICR mice were used to prepare NPCs. The experimental protocols described below were performed according to the animal experimentation guidelines of Nara Institute of Science and Technology. NPCs were prepared from telencephalons of E11.5 and E14.5 mice and cultured as described previously ##REF##10205054##[34]##. Briefly, the telencephalons were triturated in Hank's balanced salt solution (HBSS) by mild pipetting with a 1-ml pipet tip (Gilson). Dissociated cells were cultured in N2-supplemented Dulbecco's Modified Eagle's Medium with F12 (GIBCO) containing 10 ng/ml basic FGF (R&amp;D Systems) (N2/DMEM/F12/bFGF) on culture dishes (Nunc) of a chamber slide (Nunc) which had been precoated with poly-L-ornithine (Sigma) and fibronectin (Sigma). For primary astrocyte cultures, the cerebral cortices of postnatal day 1 mice were dissociated using trypsin (GIBCO), and cultured for two weeks in DMEM containing 10% fetal calf serum.</p>
<title>Results and Discussion</title> <title>Genome-wide profiling of DNA methylation</title> <p>Methylation changes in the differentiation potential switch of NPCs were analyzed by comparing mid-(E11.5) and late (E14.5)-stage NPCs using the MIAMI method ##REF##16407832##[16]##. We also compared the methylation status with that of postnatal day 1 (P1) astrocytes. Before analysis, the key genes involved in astrocyte differentiation, such as <italic>Gfap</italic> and <italic>Stat3</italic>, were analyzed for DNA methylation in the sample we used (##FIG##0##Fig. 1##). As described previously ##REF##11740937##[15]##, the STAT3-binding site of <italic>Gfap</italic> gene was demethylated in E14.5 NPCs prior to astrocyte differentiation. On the other hand, the promoter of <italic>Stat3</italic> gene was unmethylated throughout differentiation. The microarray used consisted of probes chosen from the Agilent promoter array using an eArray system (<ext-link ext-link-type="uri" xlink:href="http://earray.chem.agilent.com/earray/">http://earray.chem.agilent.com/earray/</ext-link>). The probes are located on <italic>Hpa</italic> II fragments of less than 1 kilobasepair (kb) and cover 14,543 genes. Probes which showed methylation changes at least in E14.5 NPCs or astrocytes compared to E11.5 NPCs are presented in ##FIG##1##Fig. 2A## and ##SUPPL##0##Table S1## (Name of the probes and genes were indicated in ##SUPPL##0##Table S1##). As shown in ##FIG##1##Fig. 2B, E##14.5 NPCs are hypomethylated in 85 probes (80 genes) and hypermethylated in 15 probes (15 genes). On the other hand, astrocytes are hypomethylated in 275 probes (256 genes) and hypermethylated in 170 probes (152 genes). The reliability of the analysis was confirmed by bisulfite sequencing analysis of eight genes (##FIG##2##Fig. 3##). The methylation ratio analyzed by MIAMI (##FIG##2##Fig. 3A##) had good correlation with the methylation of two <italic>Hpa</italic> II sites adjacent to the probes (##FIG##2##Fig. 3B##). Interestingly 80% of the probes hypomethylated in E14.5 NPCs are also hypomethylated in astrocytes. (##FIG##1##Fig. 2C##). If we extend the hypomethylation change criteria to the threshold level (D10&lt;0.5, usually we use D10&lt;0.2 for hypomethylation), 89% of the probes hypomethylated in E14.5 NPCs are also hypomethylated in astrocytes. In other words, 48% of the hypomethylated probes in astrocytes are also hypomethylated in E14.5 NPCs. These include a probe for an astrocyte-specific marker gene, <italic>Gfap</italic>, which was previously shown to be hypomethylated both in E14.5 NPCs and astrocytes ##REF##11740937##[15]##. However, another hypomethylated astrocyte marker, S100β ##REF##15304345##[18]##, was not detected by the MIAMI method because the DNA sequence with methylation change does not contain any <italic>Hpa</italic> II sites, leading to underestimation of the total number of genes actually demethylated in E14.5 NPCs. Nevertheless, including the genes undetectable by MIAMI, many hypomethylated sequences in astrocytes are already demethylated in E14.5 NPCs, which are competent to differentiate into astrocytes. These include important genes for astrocyte-specific function or phenotype in addition to an astrocyte-specific marker gene, <italic>Gfap</italic>. For example, Aldolase C encodes a member of the class I fructose-biphosphate aldolase gene family specific to astrocytes (<italic>Aldoc</italic>, 19). The demethylated region of <italic>Aldoc</italic> is located on exon 1. We further confirmed this result by bisulfite sequencing and found that <italic>Aldoc</italic> was methylated in E11.5 NPCs and demethylated in E14.5 NPCs and astrocytes (##FIG##2##Fig. 3B##). Another example, <italic>Kcnj10</italic> (<italic>Kir4.1</italic>), is widely expressed in astrocytes throughout the brain ##REF##17876807##[20]##. The product of this gene is absent in immature proliferating cells, and progressive expression of the genes correlates with astrocyte differentiation, which is characterized by the establishment of a negative membrane potential and exit from the cell cycle. <italic>Kcnj10</italic> (<italic>Kir4.1</italic>) encodes a member of the inward rectifier-type potassium channel family, characterized by having a greater tendency to allow potassium to flow into, rather than out of, a cell, resulting in negative membrane potential. The encoded protein is responsible for potassium buffering action, which is a major function of astrocytes and is also responsible for promoting differentiation and inhibiting cell growth ##REF##17876807##[20]##, ##REF##16306174##[21]##. Mutations in this gene have been associated with seizure susceptibility of common idiopathic generalized epilepsy syndromes ##REF##15120748##[22]##, ##REF##15725393##[23]##. <italic>Sparcl1</italic> is known as an astrocyte marker colocalized with <italic>Gfap</italic>
##REF##8869553##[24]## and is also known as a candidate gene for multiple sclerosis ##REF##17825989##[25]##. <italic>Cbs</italic> and <italic>BC055107</italic> are also known as astrocyte markers ##REF##16160063##[26]##, ##REF##18171944##[27]##. Thus, we found that demethylation occurred not only in <italic>Gfap</italic> but also in the other genes involved in astrogenesis prior to astrocyte differentiation.</p> <p>When comparing methylation status between E14.5 NPCs and astrocytes, it is important to know what percentage of the E14.5 NPCs actually gives rise to astrocytes. After 4 day-culture in the presence of LIF, 30–40% of E14.5 NPCs become GFAP-positive, yet this does not necessary mean that <italic>Gfap</italic> promoter is methylated in all GFP-negative cells. We have previously shown that STAT site within the <italic>Gfap</italic> promoter was virtually unmethylated in E14.5 NPCs during 4-day culture in the proliferating condition. They, however, could not give rise to 100% of GFAP-positive cells ##REF##11740937##[15]##. These findings suggest the existence of other factors that inhibit astrocyte differentiation. Thus, it is conceivable that demethylation is necessary but not sufficient for NPCs to differentiate into GFAP-positive astrocytes.</p> <p>Absolute levels of DNA methylation were analyzed by MIAMI and were presented as average methylation ratio for hypomethylated and hypermethylated genes compared to E11.5 NPCs (##FIG##1##Fig. 2D##). On average, hypermethylated genes were nearly fully methylated as expected. However, hypomethylated genes had some methylation on average. In addition, some hypomethylated genes, specially, <italic>Aldoc</italic> had highly mosaic pattern, suggesting a possibility that equilibrium state may have been reached rather than DNA demethylation (##FIG##2##Fig. 3B##).</p> <p>A hallmark for Polycomb-mediated repression is methylation of lysine 27 histone H3 (H3K27), which is set up by the Polycomb repressive complex2 ##REF##12408863##[28]##, ##REF##12408864##[29]##. Many polycomb targets in embryonic stem cells (ESCs) have been shown to reside in a chromatin state characterized not only by the dual presence of repressive H3K27 methylation but also by active H3K4 methylation ##REF##16630819##[30]##. In neural differentiation, the promoters marked by H3K27 methylation in ESCs are more likely to become de novo DNA methylated ##REF##18514006##[31]##. Therefore, we analyzed the histone methylation status in ESCs for our hypermethylated genes in astrocytes compared to E11.5 NPCs using the published ESCs data ##REF##17603471##[32]##. We found that promoter that were marked by H3K27 in ESCs were significantly enriched in hypermethylated genes in astrocytes (P&lt;1×10<sup>−14</sup>).</p> <title>Expression changes in NPCs and astrocytes</title> <p>Expression changes during the differentiation potential switch of NPCs were analyzed by comparing mid-(E11.5) and late (E14.5)-stage NPCs, and astrocytes by expression microarray analysis (##FIG##3##Fig. 4##, ##SUPPL##1##Table S2##). We found that 49% of the probes upregulated in E14.5 NPCs compared to E11.5 NPCs were also upregulated in astroctyes (##FIG##2##Fig. 3B##). This is smaller than the ratio of hypomethylation where 80% of the probes hypomethylated in E14.5 NPCs were also hypothylated in astrocytes (##FIG##1##Fig. 2C##). This implies that astrocyte-specific genes, including <italic>Gfap, Aldoc,</italic> and <italic>Kcnj10</italic> (<italic>Kir4.1</italic>), are demethylated at E14.5 NPCs, while they are not expressed until cells differentiate into astrocytes. To examine this hypothesis, we analyzed the expression of genes hypomethylated both in E14.5 NPCs and astrocytes. Among these genes, those that showed upregulation in E14.5 NPCs or astrocytes compared to E11.5 NPCs are shown in ##FIG##4##Fig. 5##. Of these, 64% of the genes, including astrocyte-specific marker <italic>Gfap</italic>, showed more than twice the expression level in astrocytes compared to E14.5 NPCs. Thus, genes already demethylated, such as <italic>Gfap, Aldoc,</italic> and <italic>Kcnj10</italic> (<italic>Kir4.1</italic>), were competent but not highly expressed in E14.5 NPCs, suggesting that their expression appears after differentiation into astrocytes (##FIG##5##Fig. 6##). Astrocyte-specific genes probably require the expression and/or activation of astrocyte-inducing transcription factors, such as STAT3, to be effectively expressed (##FIG##5##Fig. 6##).</p> <p>To analyze the relation between the DNA methylation and the expression of genes, we plotted the methylation ratio calculated by MIAMI and the expression level calculated by expression microarray analysis for all genes examined (##FIG##6##Fig. 7A##). As expected, the methylation ratio is inversely correlated to the expression level. Next we analyzed the average expression levels of hypermethylated and hypomethylated genes in astrocytes compared to E11.5 NPCs (##FIG##6##Fig. 7B##). Although there are some exceptions, we found the average expression levels of hypermethylated genes were low and those of hypomethylated genes were high in astrocytes. One of the reasons which could explain the existence of exceptions is the wide distribution of the probes (between −8000 to +2000 bases from the transcription start sites). Interestingly, the expression levels of hypomethylated genes were low in E14.5 NPCs compared to ones in astrocytes. This again indicates activation of expression appears after differentiation into astrocytes.</p> <title>STAT and SMAD binding sites in demethylated and upregulated genes</title> <p>A CpG sequence in a conserved STAT-binding site in the <italic>Gfap</italic> gene promoter is methylated in NPCs midgestation when astrogliogenesis does not normally commence; however, it becomes demethylated in late gestational NPCs that have gained the potential to differentiate into astrocytes ##REF##11740937##[15]##. Thus, the <italic>Gfap</italic> gene with the hypomethylated promoter can be expressed upon astrocyte differentiation induced by LIF, which activates the JAK/STAT signaling pathway ##REF##9334309##[13]##. Therefore, we searched for conserved STAT-binding sites in the genes hypomethylated both in E14.5 NPCs and astrocytes compared to E11.5 NPCs. To avoid omitting threshold level genes, the threshold used in E14.5 NPCs was extended (D10&lt;0.5). Among these genes, those that showed upregulation in E14.5 NPCs or astrocytes compared to E11.5 NPCs were used for analysis (##FIG##7##Fig. 8##). Sequences including 500-(basepair) bp upstream and 500-bp downstream of the probes were searched for STAT3 binding sites with the CpG sequence (TTN4-6AA, CG in any N position). Such methylation-sensitive STAT binding sites were found among these 18 genes; however, none were conserved between mice and humans except for <italic>Gfap</italic>, which shows hypomethylation both in E14.5 NPCs and astrocytes (##FIG##0##Fig. 1##). Although we cannot rule out the possibility that some demethylated STAT binding sites were not detected by the MIAMI method because of its target limitation, these results suggest that not all astrocyte-specific genes are epigenetically controlled only via STAT binding site methylation.</p> <p>Bone morphogenetic proteins (BMPs) also use STATs to activate the expression of astrocyte marker genes via the association between STATs and the transactivating complex composed of the MBP-activated signaling factors Smad1 and P300/CBP ##REF##8893018##[33]##, ##REF##10205054##[34]##. Therefore, we also searched for conserved SMADs-binding sites in the genes hypomethylated both in E14.5NPCs and astrocytes compared to E11.5 NPCs. We found methylation sensitive SMAD-binding sites in 6 genes; however, only <italic>Aldoc</italic> was conserved between mouse and human (##FIG##7##Fig. 8##) Bisulfite sequencing analysis of this SMAD sites reveals demethylation in E14.5 NPCs and astrocytes (##FIG##2##Fig. 3B##).</p> <p>We could not find methylation-sensitive binding sites for transcription factor except for <italic>Gfap</italic> and <italic>Aldoc</italic>. Effective activation of demethylated genes in E14.5 NPCs probably requires cooperation with other transcription factors, which are activated in astrocyte differentiation.</p> <title>Overrepresented genes with methylation changes</title> <p>Overrepresented and underrepresented categories of genes were searched in hypomethylated and hypermethylated genes. These genes were first classified into lists using PANTHER classification categories (<ext-link ext-link-type="uri" xlink:href="http://www.pantherdb.org/tools/compareToRefListForm.jsp">http://www.pantherdb.org/tools/compareToRefListForm.jsp</ext-link>). Each list was then compared to the reference list using the binomial test ##REF##10973070##[35]## for each molecular function, biological process, or pathway term in PANTHER. Unexpectedly, there was no significant overrepresentation and underrepresentation in hypomethylated genes in E14.5 NPCs and in astrocytes (P = 0.01). On the other hand, biological process categories such as “mRNA transcription regulation”, “developmental processes”, “mRNA transcription”, “neurogenesis”, “ectoderm development”, “nucleoside, nucleotide and nucleic acids metabolism” were overrepresented among hypermethylated genes in astrocytes (P = 0.01, ##TAB##0##Table 1##). A molecular function category, “transcription factor”, was also overrepresented among hypermethylated genes in astrocytes (P = 0.01, ##TAB##0##Table 1##). It is reasonable that developmental genes are shut off by DNA methylation in terminally differentiated cells. A large part of developmental genes consists of the transcription factor; therefore, it is also reasonable that a category related to transcription factor such as “mRNA transcription regulation” is overrepresented in hypermethylated genes. It is interesting to know that the pathway category, “Notch signaling pathway”, is overrepresented 9-fold more than expected among hypermethylated genes in astrocytes, although it is not statistically significant (P = 0.03). These genes include <italic>Dll1</italic>, <italic>Hes1</italic>, <italic>Hes5</italic>, <italic>Hey2</italic>, and <italic>Ncor2</italic> (##TAB##1##Table 2##). The Notch signaling pathway has been shown to be important in astrocyte differentiation ##REF##11182080##[36]##. Inactivation of developmental and transcriptional factor genes by DNA methylation is one of the important aspects of differentiation. For example, inactivation of transcription factor genes, <italic>Oct3/4</italic> (<italic>Pou5f1</italic>) and <italic>Nanog</italic>, by DNA methylation is important for early development ##REF##14761969##[6]##, ##REF##17352742##[7]##. We also found that transcription factor genes were overrepresented among methylated genes in astrocytes. Regulation by DNA methylation of transcription factor genes is an attractive model to explain development and differentiation. Usually differentiation is caused by transcription factors, which are also regulated by upstream transcription factors. So what is the most upstream regulator? One answer is that the most upstream regulator is epigenetic change.</p> <p>In conclusion, we suggested in this study that many astrocyte-specific genes are demethylated in common in late-stage NPCs, enabling cells to become competent before astrocyte differentiation. This means DNA demethylation rather than methylation is a critical regulatory event in astrocyte differentiation. We also indicated that group of genes categorized as developmental and transcription factor genes are shut off by DNA methylation in terminally differentiated astrocytes.</p>
<title>Results and Discussion</title> <title>Genome-wide profiling of DNA methylation</title> <p>Methylation changes in the differentiation potential switch of NPCs were analyzed by comparing mid-(E11.5) and late (E14.5)-stage NPCs using the MIAMI method ##REF##16407832##[16]##. We also compared the methylation status with that of postnatal day 1 (P1) astrocytes. Before analysis, the key genes involved in astrocyte differentiation, such as <italic>Gfap</italic> and <italic>Stat3</italic>, were analyzed for DNA methylation in the sample we used (##FIG##0##Fig. 1##). As described previously ##REF##11740937##[15]##, the STAT3-binding site of <italic>Gfap</italic> gene was demethylated in E14.5 NPCs prior to astrocyte differentiation. On the other hand, the promoter of <italic>Stat3</italic> gene was unmethylated throughout differentiation. The microarray used consisted of probes chosen from the Agilent promoter array using an eArray system (<ext-link ext-link-type="uri" xlink:href="http://earray.chem.agilent.com/earray/">http://earray.chem.agilent.com/earray/</ext-link>). The probes are located on <italic>Hpa</italic> II fragments of less than 1 kilobasepair (kb) and cover 14,543 genes. Probes which showed methylation changes at least in E14.5 NPCs or astrocytes compared to E11.5 NPCs are presented in ##FIG##1##Fig. 2A## and ##SUPPL##0##Table S1## (Name of the probes and genes were indicated in ##SUPPL##0##Table S1##). As shown in ##FIG##1##Fig. 2B, E##14.5 NPCs are hypomethylated in 85 probes (80 genes) and hypermethylated in 15 probes (15 genes). On the other hand, astrocytes are hypomethylated in 275 probes (256 genes) and hypermethylated in 170 probes (152 genes). The reliability of the analysis was confirmed by bisulfite sequencing analysis of eight genes (##FIG##2##Fig. 3##). The methylation ratio analyzed by MIAMI (##FIG##2##Fig. 3A##) had good correlation with the methylation of two <italic>Hpa</italic> II sites adjacent to the probes (##FIG##2##Fig. 3B##). Interestingly 80% of the probes hypomethylated in E14.5 NPCs are also hypomethylated in astrocytes. (##FIG##1##Fig. 2C##). If we extend the hypomethylation change criteria to the threshold level (D10&lt;0.5, usually we use D10&lt;0.2 for hypomethylation), 89% of the probes hypomethylated in E14.5 NPCs are also hypomethylated in astrocytes. In other words, 48% of the hypomethylated probes in astrocytes are also hypomethylated in E14.5 NPCs. These include a probe for an astrocyte-specific marker gene, <italic>Gfap</italic>, which was previously shown to be hypomethylated both in E14.5 NPCs and astrocytes ##REF##11740937##[15]##. However, another hypomethylated astrocyte marker, S100β ##REF##15304345##[18]##, was not detected by the MIAMI method because the DNA sequence with methylation change does not contain any <italic>Hpa</italic> II sites, leading to underestimation of the total number of genes actually demethylated in E14.5 NPCs. Nevertheless, including the genes undetectable by MIAMI, many hypomethylated sequences in astrocytes are already demethylated in E14.5 NPCs, which are competent to differentiate into astrocytes. These include important genes for astrocyte-specific function or phenotype in addition to an astrocyte-specific marker gene, <italic>Gfap</italic>. For example, Aldolase C encodes a member of the class I fructose-biphosphate aldolase gene family specific to astrocytes (<italic>Aldoc</italic>, 19). The demethylated region of <italic>Aldoc</italic> is located on exon 1. We further confirmed this result by bisulfite sequencing and found that <italic>Aldoc</italic> was methylated in E11.5 NPCs and demethylated in E14.5 NPCs and astrocytes (##FIG##2##Fig. 3B##). Another example, <italic>Kcnj10</italic> (<italic>Kir4.1</italic>), is widely expressed in astrocytes throughout the brain ##REF##17876807##[20]##. The product of this gene is absent in immature proliferating cells, and progressive expression of the genes correlates with astrocyte differentiation, which is characterized by the establishment of a negative membrane potential and exit from the cell cycle. <italic>Kcnj10</italic> (<italic>Kir4.1</italic>) encodes a member of the inward rectifier-type potassium channel family, characterized by having a greater tendency to allow potassium to flow into, rather than out of, a cell, resulting in negative membrane potential. The encoded protein is responsible for potassium buffering action, which is a major function of astrocytes and is also responsible for promoting differentiation and inhibiting cell growth ##REF##17876807##[20]##, ##REF##16306174##[21]##. Mutations in this gene have been associated with seizure susceptibility of common idiopathic generalized epilepsy syndromes ##REF##15120748##[22]##, ##REF##15725393##[23]##. <italic>Sparcl1</italic> is known as an astrocyte marker colocalized with <italic>Gfap</italic>
##REF##8869553##[24]## and is also known as a candidate gene for multiple sclerosis ##REF##17825989##[25]##. <italic>Cbs</italic> and <italic>BC055107</italic> are also known as astrocyte markers ##REF##16160063##[26]##, ##REF##18171944##[27]##. Thus, we found that demethylation occurred not only in <italic>Gfap</italic> but also in the other genes involved in astrogenesis prior to astrocyte differentiation.</p> <p>When comparing methylation status between E14.5 NPCs and astrocytes, it is important to know what percentage of the E14.5 NPCs actually gives rise to astrocytes. After 4 day-culture in the presence of LIF, 30–40% of E14.5 NPCs become GFAP-positive, yet this does not necessary mean that <italic>Gfap</italic> promoter is methylated in all GFP-negative cells. We have previously shown that STAT site within the <italic>Gfap</italic> promoter was virtually unmethylated in E14.5 NPCs during 4-day culture in the proliferating condition. They, however, could not give rise to 100% of GFAP-positive cells ##REF##11740937##[15]##. These findings suggest the existence of other factors that inhibit astrocyte differentiation. Thus, it is conceivable that demethylation is necessary but not sufficient for NPCs to differentiate into GFAP-positive astrocytes.</p> <p>Absolute levels of DNA methylation were analyzed by MIAMI and were presented as average methylation ratio for hypomethylated and hypermethylated genes compared to E11.5 NPCs (##FIG##1##Fig. 2D##). On average, hypermethylated genes were nearly fully methylated as expected. However, hypomethylated genes had some methylation on average. In addition, some hypomethylated genes, specially, <italic>Aldoc</italic> had highly mosaic pattern, suggesting a possibility that equilibrium state may have been reached rather than DNA demethylation (##FIG##2##Fig. 3B##).</p> <p>A hallmark for Polycomb-mediated repression is methylation of lysine 27 histone H3 (H3K27), which is set up by the Polycomb repressive complex2 ##REF##12408863##[28]##, ##REF##12408864##[29]##. Many polycomb targets in embryonic stem cells (ESCs) have been shown to reside in a chromatin state characterized not only by the dual presence of repressive H3K27 methylation but also by active H3K4 methylation ##REF##16630819##[30]##. In neural differentiation, the promoters marked by H3K27 methylation in ESCs are more likely to become de novo DNA methylated ##REF##18514006##[31]##. Therefore, we analyzed the histone methylation status in ESCs for our hypermethylated genes in astrocytes compared to E11.5 NPCs using the published ESCs data ##REF##17603471##[32]##. We found that promoter that were marked by H3K27 in ESCs were significantly enriched in hypermethylated genes in astrocytes (P&lt;1×10<sup>−14</sup>).</p> <title>Expression changes in NPCs and astrocytes</title> <p>Expression changes during the differentiation potential switch of NPCs were analyzed by comparing mid-(E11.5) and late (E14.5)-stage NPCs, and astrocytes by expression microarray analysis (##FIG##3##Fig. 4##, ##SUPPL##1##Table S2##). We found that 49% of the probes upregulated in E14.5 NPCs compared to E11.5 NPCs were also upregulated in astroctyes (##FIG##2##Fig. 3B##). This is smaller than the ratio of hypomethylation where 80% of the probes hypomethylated in E14.5 NPCs were also hypothylated in astrocytes (##FIG##1##Fig. 2C##). This implies that astrocyte-specific genes, including <italic>Gfap, Aldoc,</italic> and <italic>Kcnj10</italic> (<italic>Kir4.1</italic>), are demethylated at E14.5 NPCs, while they are not expressed until cells differentiate into astrocytes. To examine this hypothesis, we analyzed the expression of genes hypomethylated both in E14.5 NPCs and astrocytes. Among these genes, those that showed upregulation in E14.5 NPCs or astrocytes compared to E11.5 NPCs are shown in ##FIG##4##Fig. 5##. Of these, 64% of the genes, including astrocyte-specific marker <italic>Gfap</italic>, showed more than twice the expression level in astrocytes compared to E14.5 NPCs. Thus, genes already demethylated, such as <italic>Gfap, Aldoc,</italic> and <italic>Kcnj10</italic> (<italic>Kir4.1</italic>), were competent but not highly expressed in E14.5 NPCs, suggesting that their expression appears after differentiation into astrocytes (##FIG##5##Fig. 6##). Astrocyte-specific genes probably require the expression and/or activation of astrocyte-inducing transcription factors, such as STAT3, to be effectively expressed (##FIG##5##Fig. 6##).</p> <p>To analyze the relation between the DNA methylation and the expression of genes, we plotted the methylation ratio calculated by MIAMI and the expression level calculated by expression microarray analysis for all genes examined (##FIG##6##Fig. 7A##). As expected, the methylation ratio is inversely correlated to the expression level. Next we analyzed the average expression levels of hypermethylated and hypomethylated genes in astrocytes compared to E11.5 NPCs (##FIG##6##Fig. 7B##). Although there are some exceptions, we found the average expression levels of hypermethylated genes were low and those of hypomethylated genes were high in astrocytes. One of the reasons which could explain the existence of exceptions is the wide distribution of the probes (between −8000 to +2000 bases from the transcription start sites). Interestingly, the expression levels of hypomethylated genes were low in E14.5 NPCs compared to ones in astrocytes. This again indicates activation of expression appears after differentiation into astrocytes.</p> <title>STAT and SMAD binding sites in demethylated and upregulated genes</title> <p>A CpG sequence in a conserved STAT-binding site in the <italic>Gfap</italic> gene promoter is methylated in NPCs midgestation when astrogliogenesis does not normally commence; however, it becomes demethylated in late gestational NPCs that have gained the potential to differentiate into astrocytes ##REF##11740937##[15]##. Thus, the <italic>Gfap</italic> gene with the hypomethylated promoter can be expressed upon astrocyte differentiation induced by LIF, which activates the JAK/STAT signaling pathway ##REF##9334309##[13]##. Therefore, we searched for conserved STAT-binding sites in the genes hypomethylated both in E14.5 NPCs and astrocytes compared to E11.5 NPCs. To avoid omitting threshold level genes, the threshold used in E14.5 NPCs was extended (D10&lt;0.5). Among these genes, those that showed upregulation in E14.5 NPCs or astrocytes compared to E11.5 NPCs were used for analysis (##FIG##7##Fig. 8##). Sequences including 500-(basepair) bp upstream and 500-bp downstream of the probes were searched for STAT3 binding sites with the CpG sequence (TTN4-6AA, CG in any N position). Such methylation-sensitive STAT binding sites were found among these 18 genes; however, none were conserved between mice and humans except for <italic>Gfap</italic>, which shows hypomethylation both in E14.5 NPCs and astrocytes (##FIG##0##Fig. 1##). Although we cannot rule out the possibility that some demethylated STAT binding sites were not detected by the MIAMI method because of its target limitation, these results suggest that not all astrocyte-specific genes are epigenetically controlled only via STAT binding site methylation.</p> <p>Bone morphogenetic proteins (BMPs) also use STATs to activate the expression of astrocyte marker genes via the association between STATs and the transactivating complex composed of the MBP-activated signaling factors Smad1 and P300/CBP ##REF##8893018##[33]##, ##REF##10205054##[34]##. Therefore, we also searched for conserved SMADs-binding sites in the genes hypomethylated both in E14.5NPCs and astrocytes compared to E11.5 NPCs. We found methylation sensitive SMAD-binding sites in 6 genes; however, only <italic>Aldoc</italic> was conserved between mouse and human (##FIG##7##Fig. 8##) Bisulfite sequencing analysis of this SMAD sites reveals demethylation in E14.5 NPCs and astrocytes (##FIG##2##Fig. 3B##).</p> <p>We could not find methylation-sensitive binding sites for transcription factor except for <italic>Gfap</italic> and <italic>Aldoc</italic>. Effective activation of demethylated genes in E14.5 NPCs probably requires cooperation with other transcription factors, which are activated in astrocyte differentiation.</p> <title>Overrepresented genes with methylation changes</title> <p>Overrepresented and underrepresented categories of genes were searched in hypomethylated and hypermethylated genes. These genes were first classified into lists using PANTHER classification categories (<ext-link ext-link-type="uri" xlink:href="http://www.pantherdb.org/tools/compareToRefListForm.jsp">http://www.pantherdb.org/tools/compareToRefListForm.jsp</ext-link>). Each list was then compared to the reference list using the binomial test ##REF##10973070##[35]## for each molecular function, biological process, or pathway term in PANTHER. Unexpectedly, there was no significant overrepresentation and underrepresentation in hypomethylated genes in E14.5 NPCs and in astrocytes (P = 0.01). On the other hand, biological process categories such as “mRNA transcription regulation”, “developmental processes”, “mRNA transcription”, “neurogenesis”, “ectoderm development”, “nucleoside, nucleotide and nucleic acids metabolism” were overrepresented among hypermethylated genes in astrocytes (P = 0.01, ##TAB##0##Table 1##). A molecular function category, “transcription factor”, was also overrepresented among hypermethylated genes in astrocytes (P = 0.01, ##TAB##0##Table 1##). It is reasonable that developmental genes are shut off by DNA methylation in terminally differentiated cells. A large part of developmental genes consists of the transcription factor; therefore, it is also reasonable that a category related to transcription factor such as “mRNA transcription regulation” is overrepresented in hypermethylated genes. It is interesting to know that the pathway category, “Notch signaling pathway”, is overrepresented 9-fold more than expected among hypermethylated genes in astrocytes, although it is not statistically significant (P = 0.03). These genes include <italic>Dll1</italic>, <italic>Hes1</italic>, <italic>Hes5</italic>, <italic>Hey2</italic>, and <italic>Ncor2</italic> (##TAB##1##Table 2##). The Notch signaling pathway has been shown to be important in astrocyte differentiation ##REF##11182080##[36]##. Inactivation of developmental and transcriptional factor genes by DNA methylation is one of the important aspects of differentiation. For example, inactivation of transcription factor genes, <italic>Oct3/4</italic> (<italic>Pou5f1</italic>) and <italic>Nanog</italic>, by DNA methylation is important for early development ##REF##14761969##[6]##, ##REF##17352742##[7]##. We also found that transcription factor genes were overrepresented among methylated genes in astrocytes. Regulation by DNA methylation of transcription factor genes is an attractive model to explain development and differentiation. Usually differentiation is caused by transcription factors, which are also regulated by upstream transcription factors. So what is the most upstream regulator? One answer is that the most upstream regulator is epigenetic change.</p> <p>In conclusion, we suggested in this study that many astrocyte-specific genes are demethylated in common in late-stage NPCs, enabling cells to become competent before astrocyte differentiation. This means DNA demethylation rather than methylation is a critical regulatory event in astrocyte differentiation. We also indicated that group of genes categorized as developmental and transcription factor genes are shut off by DNA methylation in terminally differentiated astrocytes.</p>
<p>Conceived and designed the experiments: IH. Performed the experiments: IH. Analyzed the data: IH. Contributed reagents/materials/analysis tools: MN SM MK TH KN. Wrote the paper: IH KN.</p> <p>Epigenetic changes are thought to lead to alterations in the property of cells, such as differentiation potential. Neural precursor cells (NPCs) differentiate only into neurons in the midgestational brain, yet they become able to generate astrocytes in the late stage of development. This differentiation-potential switch could be explained by epigenetic changes, since the promoters of astrocyte-specific marker genes, glial fibrillary acidic protein (<italic>Gfap</italic>) and <italic>S100β</italic>, have been shown to become demethylated in late-stage NPCs prior to the onset of astrocyte differentiation; however, whether demethylation occurs generally in other astrocyctic genes remains unknown. Here we analyzed DNA methylation changes in mouse NPCs between the mid-(E11.5) and late (E14.5) stage of development by a genome-wide DNA methylation profiling method using microarrays and found that many astrocytic genes are demethylated in late-stage NPCs, enabling the cell to become competent to express these genes. Although these genes are already demethylated in late-stage NPCs, they are not expressed until cells differentiate into astrocytes. Thus, late-stage NPCs have epigenetic potential which can be realized in their expression after astrocyte differentiation.</p>
<title>Supporting Information</title>
<fig id="pone-0003189-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g001</object-id><label>Figure 1</label><caption><title>DNA methylation of <italic>Gfap</italic> and <italic>Stat3</italic> gene.</title><p>DNA methylation was analyzed by bisulfite sequencing. Closed circles indicate methylated CpG sites and open circles indicate unmethylated CpG sites. STAT3 binding sites are indicated by arrows.</p></caption></fig> <fig id="pone-0003189-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g002</object-id><label>Figure 2</label><caption><title>Genome-wide profiling of DNA methylation.</title><p>(A) Hypomethylated and hypermethylated genes in E14.5 NPCs and astrocytes compared to E11.5 NPCs. Each row indicates each probe and each column indicates E14.5 NPCs or astrocytes. Probes with methylation changes are indicated in color. D10 is a value indicating the difference between methylation-sensitive <italic>Hpa</italic> II cleavage and methylation-insensitive <italic>Msp</italic> I cleavage. (B) Number of methylation changes in E14.5 NPCs and astrocytes compared to E11.5 NPCs. Numbers of methylation changes indicated in ##FIG##0##Fig. 1A## were counted and indicated by bars. (C) Venn diagram of the number of probes with methylation changes in E14.5 NPCs and astrocytes compared to E11.5 NPCs. (D) Absolute levels of DNA methylation analyzed by MIAMI presented as methylation ratio. Methylation ratios of hypermethylated and hypomethylated probes are presented as average±standard deviation.</p></caption></fig> <fig id="pone-0003189-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g003</object-id><label>Figure 3</label><caption><title>Confirmation of genome-wide methylation profiling by MIAMI.</title><p>(A) Methylation ratios of seven hypomethylated and one hypermethylated gene compared to E11.5 NPCs. Methylation ratios are calculated by MIAMI data. (B) Confirmation of seven hypomethylated genes and one hypermethylated gene by bisulfite sequencing analysis. Closed circles indicate methylated CpG sites and open circles indicate unmethylated CpG sites. Positions of <italic>Hpa</italic> II sites are indicated as arrows.</p></caption></fig> <fig id="pone-0003189-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g004</object-id><label>Figure 4</label><caption><title>Expression analysis by microarray.</title><p>(A) Number of upregulated genes in E14.5 NPCs and astrocytes compared to E11.5 NPCs. (B) Venn diagram of the number of upregulated genes in E14.5 NPCs and astrocytes compared to E11.5 NPCs.</p></caption></fig> <fig id="pone-0003189-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g005</object-id><label>Figure 5</label><caption><title>Expression of the genes hypomethylated both in E14.5 NPCs and astrocytes compared to E11.5 NPCs.</title><p>Among hypomethylated genes, genes which showed upregulation in E14.5 NPCs or astrocytes compared to E11.5 NPCs are shown. Expression is calculated from data obtained by expression microarray analysis and presented as fold change compared to E11.5 NPCs. <italic>Gfap</italic> is indicated by an arrow.</p></caption></fig> <fig id="pone-0003189-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g006</object-id><label>Figure 6</label><caption><title>Model proposed for astrocyte-specific expression of the genes demethylated in E14.5 NPCs.</title><p>Several astrocyte-specific genes are already demethylated in E14.5 NPCs; however, they are induced after diffentiation to astrocytes. This could be explained by astrocyte-specific transcription factors.</p></caption></fig> <fig id="pone-0003189-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g007</object-id><label>Figure 7</label><caption><title>DNA methylation and expression.</title><p>(A) All probes are grouped by their average methylation ratios and their expression levels are indicated as average±standard error of mean. Expression levels are described as normalized values in log scale. (B) Expression levels of hypermethylated (red) and hypomethylated (blue) genes in astrocytes. Expression levels are indicated as average±standard error of mean.</p></caption></fig> <fig id="pone-0003189-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.g008</object-id><label>Figure 8</label><caption><title>Screening of conserved STAT and SMAD binding sites in demethylated and upregulated genes.</title><p>Schematic representation of screening of conserved STAT and SMAD binding sites in hypomethylated genes both in E14.5 NPCs and astrocytes. Numbers of genes are indicated in parentheses.</p></caption></fig>
<table-wrap id="pone-0003189-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.t001</object-id><label>Table 1</label><caption><title>Overrepresented category of hypermethylated genes in astrocytes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Biological Process</td><td align="left" rowspan="1" colspan="1">Obs/Exp</td><td align="left" rowspan="1" colspan="1">P-value</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Mrna transcription regulation</td><td align="left" rowspan="1" colspan="1">3.7</td><td align="left" rowspan="1" colspan="1">1.30E-10</td></tr><tr><td align="left" rowspan="1" colspan="1">Developmental processes</td><td align="left" rowspan="1" colspan="1">2.8</td><td align="left" rowspan="1" colspan="1">1.35E-10</td></tr><tr><td align="left" rowspan="1" colspan="1">Mrna transcription</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">4.32E-08</td></tr><tr><td align="left" rowspan="1" colspan="1">Neurogenesis</td><td align="left" rowspan="1" colspan="1">3.5</td><td align="left" rowspan="1" colspan="1">4.67E-04</td></tr><tr><td align="left" rowspan="1" colspan="1">Endoderm development</td><td align="left" rowspan="1" colspan="1">3.2</td><td align="left" rowspan="1" colspan="1">5.93E-04</td></tr><tr><td align="left" rowspan="1" colspan="1">Nucleoside, nucleotide and nucleic acid metabolism</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">3.19E-03</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Molecular Function</td><td align="left" rowspan="1" colspan="1">Obs/Exp</td><td align="left" rowspan="1" colspan="1">P-value</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Transcription factor</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">4.29E-09</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003189-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003189.t002</object-id><label>Table 2</label><caption><title>Hypermethylated genes in astrocytes belong to Notch signaling pathway.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene Symbol</td><td align="left" rowspan="1" colspan="1">Gene Name</td><td align="left" rowspan="1" colspan="1">Probe Pos</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Dll1</td><td align="left" rowspan="1" colspan="1">delta-like 1</td><td align="left" rowspan="1" colspan="1">1725</td></tr><tr><td align="left" rowspan="1" colspan="1">Hes1</td><td align="left" rowspan="1" colspan="1">hairy and enhancer of split 1</td><td align="left" rowspan="1" colspan="1">−2503</td></tr><tr><td align="left" rowspan="1" colspan="1">Hes5</td><td align="left" rowspan="1" colspan="1">hairy and enhancer of split 5</td><td align="left" rowspan="1" colspan="1">−7260</td></tr><tr><td align="left" rowspan="1" colspan="1">Hey2</td><td align="left" rowspan="1" colspan="1">hairy/enhancer-of-split related with YRPW motif 2</td><td align="left" rowspan="1" colspan="1">−1761</td></tr><tr><td align="left" rowspan="1" colspan="1">Ncor2</td><td align="left" rowspan="1" colspan="1">nuclear receptor co-repressor 2</td><td align="left" rowspan="1" colspan="1">−5790</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0003189.s001"><label>Table S1</label><caption><p>Table S1 Genome-wide profiling of DNA methylation. Hypomethylated and hypermethylated genes in E14.5 NPCs and astrocytes compared to E11.5 NPCs. Each row indicates each probe. Probes with methylation changes are indicated by color as in ##FIG##0##Fig. 1A##.</p><p>(0.06 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003189.s002"><label>Table S2</label><caption><p>Table S2 List of upregulated genes compared to E11.5 NPCs. Upregulated genes in E14.5 NPCs and astrocytes compared to E11.5 NPCs.</p><p>(0.55 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003189.s003"><label>Table S3</label><caption><p>Table S3 List of primers and restriction enzymes used for COBRA and bisulfite sequencing analysis.</p><p>(0.02 MB XLS)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This study was supported in part by grants from the Japanese Science and Technology Agency, the Ministry of Education, Culture, Sports, Science and Technology of Japan, and the Ministry of Health, Labor, and Welfare of Japan. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0003189.g001"/> <graphic xlink:href="pone.0003189.g002"/> <graphic xlink:href="pone.0003189.g003"/> <graphic xlink:href="pone.0003189.g004"/> <graphic xlink:href="pone.0003189.g005"/> <graphic xlink:href="pone.0003189.g006"/> <graphic xlink:href="pone.0003189.g007"/> <graphic xlink:href="pone.0003189.g008"/> <graphic id="pone-0003189-t001-1" xlink:href="pone.0003189.t001"/> <graphic id="pone-0003189-t002-2" xlink:href="pone.0003189.t002"/>
<media xlink:href="pone.0003189.s001.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003189.s002.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003189.s003.xls"><caption><p>Click here for additional data file.</p></caption></media>
[]
PLoS One. 2008 Sep 11; 3(9):e3189