Abstract
PAX6-null mice exhibit defects in multiple organs leading to neonatal lethality, but the mechanism by which this occurs has not yet fully elucidated. In this study, we generated induced pluripotent stem cells (iPSCs) from Pax6-mutant mice and investigated the effect of PAX6 on cell fate during embryoid body (EB) formation. We found that PAX6 promotes cell migration by directly downregulating miR-124, which is important for the fate transition of migratory cells during gastrulation of embryonic stem (ES) cells. Although several downstream targets of miR-124 have been reported, little is known regarding the upstream regulation of miR-124. When we observed EB formation of iPSCs from Pax6-mutant mice, we found that higher levels of miR-124 in Pax6 homozygous EBs (Homo-EBs) inhibited cell migration, whereas inhibition of miR-124 in Homo-EBs rescued the migratory phenotypes associated with PAX6 deficiency. Further, we found that PAX6 binds to the promoter regions of the miR-124-3 gene and directly represses its expression. Therefore, we propose a novel PAX6-miR-124 pathway that controls ES cell migration. Our findings may provide important information for studies on ES cell differentiation and embryonic development.
Introduction
M
The formation of embryoid bodies (EBs) is a principal step in the differentiation of ES cells in vitro [15]. The molecular and cellular morphogenic signals and events that occur within EBs recapitulate numerous aspects of embryonic development and result in cellular differentiation to the three embryonic germ layers (ie, endoderm, mesoderm, and ectoderm), similar to gastrulation of an epiblast-stage embryo in vivo [16]. The specific identity and spatial coordination of the various cell–cell interactions involved in EB formation are thought to influence the course of ES cell differentiation [17,18]. Recently, induced pluripotent stem cells (iPSCs) have been developed as a valuable tool to generate pluripotent stem cells from the somatic cells of patients. The formation of EBs from disease-specific iPSCs can be useful for exploration of disease mechanisms both in vitro and in vivo [19].
PAX6 is a transcription factor essential for the development of tissues, including those of the eyes, central nervous system (CNS), and endocrine glands of vertebrates and invertebrates. PAX6-null mutant mice exhibit defects in multiple organs, including the eye, brain, and pancreas, that lead to neonatal lethality [20 –24]. PAX6 regulates the expression of a broad range of molecules, including transcription factors, cell adhesion and short-range cell–cell signaling molecules, hormones, and structural proteins. PAX6 has also been implicated in various key biological processes, including cell proliferation, migration, adhesion, and signaling, both in normal development and in oncogenesis [25]. Our previous study demonstrated that PAX6 deficiency caused downregulation of the Pc1 gene and upregulation of the Pcsk1n gene, which in turn caused diabetes [26,27]. PAX6 appears to be important for neurogenesis in the subventricular zone (SVZ), which is abnormal in PAX6-deficient (Sey/Sey) mice [28], also for genesis and glutamatergic differentiation of late-born neurons [29]. The SVZ defect in the Sey/Sey cortex suggests that CNS development is impaired by Pax6 mutation.
In this study, we generated iPSCs from PAX6-deficient mice and then generated EBs to investigate early embryonic differentiation. Using a spontaneous EB differentiation protocol, we found that cell adhesion/migration of Pax6-homozygous EB (Homo-EB) was dramatically impaired relative to that of wild-type EB (Wt-EB). Many studies have reported that PAX6 regulates cell migration through its downstream molecules, mostly proteins [30,31], but few or no studies have investigated miRNA with respect to this issue, even though it has recently received a great deal of attention because of its important regulatory roles in various cellular processes, including cell migration [32]. The apparent importance of miRNA in cellular regulation, especially the effect of miR-124 on the fate transition of migratory cells [14], prompted us to test whether miRNA mediates the cell migration regulated by PAX6 during early development.
Materials and Methods
Mice
Animals (Pax6 R266Stop mutant mice) were used as previously described [26]. Briefly, a stop-codon mutation at a 266 aa of PAX6 was produced by injecting the chemical mutagen N-ethyl-N-nitrosourea (ENU); the truncated PAX6 could consider being null and the mutation was confirmed to be transmitted to the offspring. Pax6 heterozygous mutant mice show a small-eye phenotype, whereas homozygous mutant embryos lack eyes and die soon after birth. Pax6 heterozygous mutant mice were mated, mouse embryos collected 13.5–14.5 days post-coitus were harvested for mouse embryonic fibroblast (MEF) isolation, and the embryos were genotyped by polymerase chain reaction (PCR) analysis as previously described [26]. All animal studies were performed in accordance with Institutional Animal Care and Use Committee protocols at Peking University.
Cell culture
MEF cells of the three genotypes used for reprogramming in this study were obtained from different embryos from the same uterus. MEF isolation was performed as previously described [33]. Cells were cultured in high-glucose Dulbecco's modified Eagle's medium (DMEM) (Hyclone) containing 10% (v/v) fetal bovine serum (FBS) (Biochrom), and MEFs within three passages were used in this study. Mouse embryonic stem cells and iPSCs were maintained on mitomycin-c-treated MEFs in mouse embryonic stem cell (mESC) medium containing high-glucose DMEM (Hyclone), 20% (v/v) ES-cell-qualified FBS (Biochrom), 2 mM
Generation and culture of mouse iPSCs
The generation and structure of doxycycline-controlled Tet-on-inducible lentiviruses expressing mouse-derived Oct4, Sox2, c-Myc, and Klf4 have been described in previous studies [34]. MEF cells of three genotypes in 35-mm tissue culture plates were sequentially infected with lentiviruses two times in the presence of 8 g/mL polybrene (Sigma). After 2 days of viral infection, the culture medium was replaced by mESC medium supplemented with 2 mg/mL doxycycline (Dox) to induce reprogramming. The newly generated mouse iPS colonies were picked up after 6–8 days.
Immunostaining and alkaline phosphatase staining
iPSCs were fixed with 4% (w/v) paraformaldehyde for 30 min and permeabilized with 0.5% (v/v) Triton X-100 in phosphate-buffered saline (PBS) for 15 min at room temperature. After blocking with 5% (w/v) bovine serum albumin for 30 min, the cells were incubated with primary antibodies at 4°C overnight. The cells were washed with PBS and then incubated with secondary antibodies [DyLight™ 594-conjugated AffiniPure goat anti-rabbit immunoglobulin G (IgG) and DyLight™ 488-Conjugated AffiniPure Goat Anti-Mouse IgG; 1:200 dilution] and visualized by fluorescence microscopy (Olympus) or confocal microscopy (Olympus FV1000), after counterstaining with 4,6-diamidino-2-phenylindole. The primary antibodies used for characterization of iPSCs targeted SSEA-1 (1:100; Santa Cruz) and Nanog (1:100; Abcam). Alkaline phosphatase (AP) staining was performed according to the manufacturer's recommendations (Millipore). For EB staining, EBs were plated onto gelatin-coated cover glass (Fisher Scientific); they attached to the glass and expanded spontaneously, and then the slides were stained as described previously. The antibody used in this experiment was the anti-IQGAP1 rabbit polyclonal antibody (1:100; Bioss).
DNA microarray
Total RNA from ES cells, iPSCs, MEFs, or EBs was labeled with biotin. The samples were hybridized to a mouse oligo microarray (Affymetrix) according to the manufacturer's protocol. The arrays were scanned on Affymetrix scanners and the expression value for each gene was calculated using the GENECHIP software (Affymetrix).
Teratoma formation
iPSCs were suspended at 1×107 cells/mL in DMEM containing 10% FBS. We subcutaneously injected 100 μL of the cell suspension (1×106 cells) into the dorsal flanks of severe combined immunodeficient (SCID) mice. At 3–4 weeks after injection, the tumors were explanted, fixed in 4% paraformaldehyde, embedded in paraffin, and examined histologically using hematoxylin and eosin staining.
Spontaneous differentiation and formation of EBs
The three-step spontaneous differentiation protocol for mES cells or iPSCs was adapted from two protocols [35,36], with some modifications. Briefly, our protocol included the depletion of feeder cells, formation of EBs, and spontaneous differentiation. Details for the methods are provided in the Supplementary Data (Supplementary Data are available online at
In vitro migration assays
All migration experiments with EBs were performed in mES media without LIF on plastic substrates. Culture dishes were handled carefully to keep the EBs at the bottom, and cell migration was evaluated by measuring the distance from the edge of the EB to the nucleus of the most distant cell. Different viewing fields were randomly chosen for counting, and average values were determined. The assays were repeated at least three times.
Transwell migration assay
Transwell migration assays were performed in six-well culture plates with inserts (8-μm pore size; BD Biosciences). iPSCs were transferred to the upper-chamber inserts. In total, 2.5 mL of differentiation medium was added to the lower and upper chambers. The cells were allowed to migrate for 4 days in a humidified CO2 incubator at 37°C. Following incubation, cells on the lower side of the insert filter were quickly fixed and stained. Different viewing fields were randomly chosen for counting, and average values were determined. The assays were repeated three times.
Proliferation and apoptosis assay
Four-day-old EBs were allowed to adhere to gelatin-coated dishes for 18 or 48 h, after which they were dissociated from dishes and digested for cell counting. Cell apoptosis was determined by flow cytometry using propidium iodide and Annexin-V staining.
Western blotting
Sodium dodecyl sulfate–polyacrylamide gel electrophoresis (PAGE) and western blotting were performed according to standard procedures [26]. β-Actin was used on the same membrane as the loading control. The primary antibodies used were anti-IQGAP (1:500; Bioss), anti-PAX6 (1:1,000; C-Terminus, Millipore), and anti-β-actin (1:1,000; Santa Cruz). The experiments were repeated three times.
miRNA analysis by quantitative PCR
Total RNA containing miRNAs was extracted from cultured cells or EBs by using the miRNeasy mini kit (Qiagen) according to the manufacturer's protocol. The isolated total RNA (1 μg) containing miRNA was used as starting material for the reverse transcription reaction by using the miScript ii RT kit (Qiagen). The miScript PCR system (Qiagen), which comprises the miScript ii RT Kit, miScript SYBR Green PCR Kit, miScript primer assay, and miScript precursor assay, enables quantification of mRNA, miRNA, and precursor miRNA from the same cDNA. The Mm_miR-124_1 miScript primer assay (MS00029211; Qiagen) was used to detect the expression of mature mmu-mir-124. The Mm_mir-124-1_PR_1 miScript Precursor Assay, Mm_mir-124-1_PR_2 miScript Precursor Assay, and Mm_mir-124-1_PR_3 miScript Precursor Assay (MP00004074, MP00004081, and MP00004088; Qiagen) were used to detect the expression of 3 mmu-mir-124 precursors.
Northern blotting
Total RNA was extracted from iPSCs and EBs by using the miRNeasy mini kit (Qiagen). Northern blotting was performed using the highly sensitive miRNA Northern Blot Assay Kit (NB-1001; Signosis). Briefly, 12 μg of total RNA from each sample was separated on a 15% polyacrylamide gel, transferred to a membrane (Signosis), and UV-crosslinked. The membranes were incubated with a biotin-labeled miR-124 probe (HP-0049; Signosis) and a U6 probe (HP-1001; Signosis). Prehybridization and hybridization were conducted using a hybridization buffer, and the membranes were then exposed using a biotin-based chemiluminescence imaging system.
In situ hybridization
In situ hybridization was performed by following the Exiqon protocol with some modifications. Briefly, EBs were fixed with 4% fresh paraformaldehyde, followed by paraffin embedding. Formalin-fixed paraffin-embedded sections (6-μm thick) were cut and stored at 4°C, and paraffin was melted on the day before the in situ hybridization experiment was conducted. The slides were deparaffinized and pretreated with 10 μg/mL proteinase K for 10 min at room temperature and then dehydrated at room temperature. Subsequently, hybridization was conducted using the hybridization mix at 55°C for 2 h. The probes were labeled using digoxigenin, and in situ hybridization signals were detected using an anti-digoxigenin antibody conjugated with AP (AP-conjugated anti-DIG; Roche). AP was used to convert the soluble substrates 4-nitro-blue tetrazolium (NBT) and 5-bromo-4-chloro-3′-indolylphosphate (BCIP) into a water- and alcohol-insoluble dark-blue NBT-BCIP precipitate, and the samples were then counterstained with 0.1% Nuclear Fast Red and mounted. A scrambled-sequence probe (LNA™; Exiqon) was used as a control in the hybridization.
Recombinant adenovirus construction
Full-length Pax6 cDNA and a specific Pax6 shRNA (5′-GATCCCCAACTCCATCAGTTCTAATTCAAGAGATTA GAACTGATGGAGTTGG-3′, 5′-AGCTCCAACTCCATC AGTTCTAATCTCTTGAATTAGAACTGATGGAGTT GGG-3′) were subcloned into the adenovirus shuttle vector pDC316-EGFP and pDC316-EGFP-U6, respectively. Recombinant Ad-GFP, Ad-PAX6-GFP, Ad-Pax6-RNAi, and Ad-non-silence were generated as previously described [37,38], using the AdMax system (Microbix Biosystems). The adenovirus was purified using a ViraBind™ Adenovirus Purification Kit (Cell Biolabs). Mouse iPSCs were detached from culture dishes, grown in suspension culture for 1–2 days, and transfected with adenovirus. After 2–3 days, the cells were subjected to western blotting and quantitative PCR (qPCR).
Chromatin immunoprecipitation
Chromatin immunoprecipitation (ChIP) assays were performed as described previously [38] for EBs derived from the mouse ES cell line R1. Briefly, chemical-linked complexes of genomic DNA and PAX6 protein were immunoprecipitated with antibodies specific to PAX6 (C-Terminus; Millipore). The DNA was amplified for 21–25 cycles after ChIP using four pairs of primers: spanning −7,939 to −7,710 bp: forward, 5′-ACTGAGATCTGGCAGGGAGA-3′ and reverse, 5′-AGTGCATGTCATTGGCTCTG-3′; spanning −5,169 to −4,958 bp: forward, 5′-GAGCCTGTTGTGATGCTTGA-3′ and reverse, 5′-CAGTGTGCACGTGGAGTAGG-3′; spanning −4,049 to −3,758 bp: forward, 5′-GAAAGCACACCC TCGGAAT-3′ and reverse, 5′-TCTAACATGGGTGGTCA AGG-3′; spanning −1,754 to −1,475 bp: forward, 5′-GAGG CCTGCC TACTTCTGTG-3′ and reverse, 5′-TCCCCTTTTT CACCCTTTTT-3′.
Luciferase activity
We amplified 2 mmu-mir-124-3 DNA fragments (−1,754 bp/+28 bp, −4,394 bp/−3,432 bp), including putative PAX6-binding sites. According to the manufacturer's instructions, the −1,754 bp/+28 bp DNA fragment was cloned into the pGL3-basic vector (Promega) by using Kpn1 and Nhe1 (New England Biolabs), and the −4,394 bp/−3,432 bp fragment was cloned into the pGL3-promoter vector (Promega) by using Mlu1 and Xho1 (New England Biolabs). Correct insertion was verified by sequencing. The PCR primers were as follows: (−1,754/+28) forward, 5′-AA AGGTACCGAGGCCTGCCTACTTCTGTG-3′ and reverse, 5′-GTGGGCTAGCAATCAAGG TCCGCTGTG AAC-3′; (−4,394/−3,432) forward, 5′-ATCGACGCGTTCCCTGGG TTTTCATCTG AG-3′ and reverse, 5′-ATCGCTCGAGGGC TTGTCATTCCTCCAGAA-3′. The PAX6 and truncated PAX6 (R240Stop mutation) expression vectors (pcDNA3.1-/Pax6 and pcDNA3.1-/Pax6-tru) were constructed as previously described [37]. HEK293 cells (1×105) were seeded in a 24-well plate and transiently transfected with pGL3-basic −1,754/+28 or pGL3-promoter −4,394/−3,432, along with PAX6 (wild or truncated type) at several concentrations (0, 400, and 800 ng/well) (the empty vector was used to supply up to 800 ng/well as control) by using Lipofectamine 2000 (Invitrogen), and pRL-CMV (Promega) expressing Renilla luciferase was transfected as an internal control. Forty hours after transfection, cell lysates were extracted in passive lysis buffer (Promega), and luciferase activity was measured using the dual luciferase reporter assay system (Promega) in a Centro LB960 96-well luminometer (Berthold Technologies). The deleted promoter (−1,533/+28) was amplified using the following primers: forward, 5′-AAAGGTACCCTCCCATCC TCGTAAGCAAG-3′ and reverse, 5′-GTGGG CTAGCAAT CAAGGTCCGCTGTGAAC-3′. This sequence was then cloned into the pGL3-basic vector, and the luciferase assay was performed as described previously for the −1,754/+28 fragment.
Electrophoretic mobility shift assay
Nuclear extracts (NEs) from mES cell R1-derived EBs were prepared using the NE-PER® Nuclear and Cytoplasmic Extraction Reagent (Pierce). Electrophoretic mobility shift assay (EMSA) was conducted using the LightShift® Chemiluminescent EMSA Kit (Pierce) following the manufacturer's protocols. Custom-synthesized 5′ single-strand nucleotides labeled at the 5′-end with biotin (only the sense strands are shown) and covering −1,636 bp/−1,577 bp (F: 5′-CGCTGA TCTCTAGCCCTAAAGACCCAGATCAA
miRNA transfection
An mmu-miR-124 mimic (MSY000134; Qiagen), inhibitor (MIN000134; Qiagen), and negative controls (1027280 and 1027271) were purchased from Qiagen. Transient transfections were performed using Lipofectamine RNAi Max (Invitrogen) according to the manufacturer's instructions. In brief, iPSCs were grown in differentiation medium on low-attachment plates or the upper chambers of transwells for 1–2 days; the EBs were transfected twice on 2 consecutive days with 100 nM miR-124 mimic or inhibitor and then plated on gelatin-coated dishes or transwells for 1–2 days for migration assay.
Statistical analysis
Statistical analysis was performed using an independent sample t-test. The asterisks in each graph indicate statistically significant differences with P-value calculated by t-test and P<0.05 was considered statistically significant.
Results
Generation of iPSCs from PAX6 deficiency mice
To establish a cellular model to address the role of PAX6 in early embryonic development in mice, we generated iPSCs using MEFs from Pax6 R266Stop mutant mice [26]. These ENU-induced mutagenic mice exhibited a C-to-T mutation in the Pax6 gene that generated a polypeptide truncated at residue 266 of the PAX6 protein [26,27]. Because the mutation results in production of a truncated protein lacking the C-terminal domain, the truncated PAX6 is considered functionally inactive. Therefore, mouse homozygous mutants are widely used as Pax6-null mutants [39].
iPSC lines of three different Pax6 genotypes (designated as Wt-iPSCs, Hetz-iPSCs, and Homo-iPSCs) were established from normal siblings as well as heterozygous and homozygous Pax6 mutants; three or four independently validated iPSC clones per genotype were generated and at least two independent clones of each genotype were assayed in each experiment. Representative experimental results are shown here (Fig. 1A; Supplementary Fig. S1A). All iPSC lines exhibited morphological features similar to those of mouse ES cells—strong AP activity (Fig. 1A), expression of pluripotency markers (including Nanog and SSEA-1; Fig. 1B), and normal karyotypes (Supplementary Fig. S1B). Global gene expression was analyzed in Wt-iPSCs, Hetz-iPSCs, Homo-iPSCs, MEFs, and mouse R1 ES cells by using oligonucleotide microarrays. A hierarchical clustering analysis clearly indicated that Wt-iPSCs, Hetz-iPSCs, and Homo-iPSCs clustered closely with R1 ES cells and were distinct from MEFs (Fig. 1C). Analysis of scatter plots further demonstrated that Wt-iPSCs, Hetz-iPSCs, and Homo-iPSCs are highly similar to R1 ES cells and hence completely different from MEFs at the global transcriptional level.

Induced pluripotent stem cells (iPSCs) generated from mouse embryonic fibroblast (MEF) cells of wild-type mice and Pax6 mutants exhibit mouse embryonic stem (mES) cell characteristics.
To confirm epigenetic remodeling in the reprogrammed cells, we also evaluated the methylation status of CpG dinucleotides in the Oct4 and Nanog promoters. Bisulfite genomic sequencing analysis showed that the Oct4 and Nanog promoter regions were demethylated in Wt-iPSCs, Hetz-iPSCs, and Homo-iPSCs relative to MEFs, whereas their methylation was similar to that of R1 (Fig. 1E). The pluripotent potential of these three types of iPSCs was further indicated by teratoma formation. After iPSCs were transplanted into the hind limbs of SCID mice for 3–4 weeks, hematoxylin and eosin staining (Fig. 1F) and immunohistochemistry staining (Supplementary Fig. S1C) revealed the presence of derivatives of all three germ layers, including neural epithelium (βIII-tubulin), muscle α-smooth muscle actin, and endoderm (GATA4), within the tumors.
Migration of cells from the outer layer of EBs is impaired in Pax6 Homo-EBs
On analyzing EB formation, we detected differences between Wt-iPSCs and Homo-iPSCs. At 18 or 48 h after plating onto the dishes, the number of migrated cells and the distance from the edge of the EB to the most distant migrating cell were significantly reduced in Homo-EBs relative to those in Wt-EBs (Fig. 2A, B). Wt-EBs formed a halo around the clones, whereas Homo-EBs remained largely intact, with far fewer cells migrating. Time-lapse videos also showed that Wt-EBs were flattened and the cells adhered and then migrated quickly, while Homo-EBs remained a status of tight, domed shape and the cells showed significantly decreased migratory capacity (Supplementary Video). Moreover, when we checked the borders of teratomas derived from Wt- and Homo-iPSCs, we found relatively intact borders in homo teratomas, but no clear border was observed in Wt teratomas (Supplementary Fig. S1D).

Pax6 mutation limits migration of embryoid body (EB) cells in vitro.
This difference was quantified using a transwell migration assay (Fig. 2C, D). To exclude the possibility that the observed difference in cell number resulted from cell proliferation or apoptosis, we performed cell counting and Annexin V-fluorescence-activated cell sorting analysis and detected no difference between Wt- and Homo-iPSCs and EBs (Supplementary Fig. S2), which is consistent with a previous report; Quinn et al. showed that the proportions of BrdU-labeled proliferating cells were not different between Wt and Pax6−/− EBs [40]. All these findings suggested that cell migration of Homo-EBs is dramatically impaired.
Because PAX6 has been reported to be involved in the progression of differentiation toward neurons (ectoderm) and the pancreas (endoderm) [41], as well as in the migratory behavior of differentiating cells [42], to determine whether the migration defect at day 6 of EB formation resulted from blockade of differentiation progression by Pax6 mutation, we analyzed the spontaneous differentiation of Wt- and Homo-EBs using reverse-transcription PCR, qPCR, and immunofluorescence (Supplementary Fig. S3B–G). Concurrent with the downregulation of pluripotency genes (Oct4, Nanog, and Klf4), the expression of ectoderm (Sox1 and Nestin), mesoderm (Brachyury), and endodermal (Sox17) marker genes was dramatically increased to similar levels in both groups. In addition, cDNA microarray analysis was performed to compare the globe gene expression of the Wt- and Homo-EBs (5–6-day old) (Supplementary Fig. S3H). Although some slight differences were observed, the results confirmed that the expression profiles of ectoderm-, mesoderm-, and endoderm-specific genes were very similar between the two groups. All these observations suggested that PAX6 deficiency did not impair the progression of differentiation during 6 days of EB formation, which is consistent with the normal early embryo development of mice with defective Pax6 [43]. It indicates that at day 6 during the differentiation of EBs, the migration difference induced by Pax6 mutation resulted from the regulation of migration by PAX6 and its downstream molecules.
miR-124 expression is elevated in Homo-EBs
Cell migration from the outer regions of EBs has been reported to be regulated by IQGAP1, which is a Ras GTPase-activating-like protein preferentially located at the outer regions of EBs [14]. To determine IQGAP1 expression patterns in Wt-EBs and Homo-EBs, we plated the EBs on a dish for 1–2 days and then performed immunochemistry and western blotting analyses. IQGAP1 expression decreased in Homo-EBs at the outer region of the EBs (Fig. 3A) or in the whole cell lysate (Fig. 3B). Because IQGAP1 expression is inhibited by miR-124 during EB differentiation [14], we measured the expression of miR-124 in Wt- and Homo-EBs. Consistent with the findings of Lee et al. [14], miR-124 expression decreased significantly during the formation of Wt-EBs, but there was no significant decrease in Homo-EBs. Thus, the expression of miR-124 in Homo-EBs was significantly higher than that in Wt-EBs at day 4 (Fig. 3C). Three genes encode the mature sequence of miR-124—miR-124-1, miR-124-2, and miR-124-3—which are located, respectively, on chromosomes 14, 3, and 2 of the mouse genome [44]. Using Northern blotting analysis (Fig. 3D) and in situ hybridization (Fig. 3E), we further confirmed that the expression of both precursor and mature miR-124 was significantly higher in Homo-EBs than in Wt-EBs. Because the probe for mature miR-124 also recognized all three types of re-miR-124 molecules, which have very similar molecular weights, only one band for pre-miR-124 appears in the Northern blotting results (Fig. 3D, upper line). Based on these data, we hypothesized that the expression of miR-124 was affected by PAX6 during EB formation.

miR-124 inhibits migration of cells from Homo-EBs.
Expression of miR-124 is regulated by PAX6
To test this hypothesis, we investigated the changes in expression for mature miR-124 and its three precursor mRNAs in EBs with overexpression or knockdown of Pax6. As shown in Fig. 4, overexpression of Pax6 in Homo-EBs (Fig. 4A) resulted in a decrease in the expression of mature miR-124 (Fig. 4B). To analyze the transcriptional regulation of miR-124, qPCR was used to determine the expression of its precursors, that is, pre-miR-124-1, pre-miR-124-2, and pre-miR-124-3, from three different DNA fragments cleaved and processed to mature miR-124 by a nuclear enzymatic complex consisting of nuclear RNAase III and the adapter protein Drosha [45]. As shown in Fig. 4C–E, pre-miR-124-3, but not pre-miR-124-1 or pre-miR-124-2, was responsible for the decrease in mature miR-124. As expected, knockdown of Pax6 in Wt-EBs resulted in increased expression of pre-miR-124-3 and mature miR-124 (Fig. 5B–E). We confirmed the results using two siRNAs for Pax6 (siRNA ID: s71268 or s71269 Ambion) and two scrambled siRNA9 (Cat. No. AM4611 and AM4613), and used low concentrations (50 nM) to reduce off-target effects. The decrease in PAX6 expression was detected by western blotting analysis (Fig. 5A).

Overexpression of Pax6 results in decreased miR-124 expression in Homo-EBs.

Knockdown of Pax6 results in increased miR-124 expression in Wt-EBs.
Transcription of miR-124 is directly inhibited by PAX6
To determine whether miR-124 is the direct inhibitory target of PAX6, we scanned 10-kb sequences upstream of the miR-124-3 gene and found several putative PAX6-binding sites based on the TF SEARCH (Searching Transcription Factor Binding Sites,

PAX6 represses miR-124-3 transcription by binding to the promoter region of miR-124-3.
To investigate how mutation of Pax6 leads to higher expression of pre-miR-124-3, we performed a luciferase reporter assay using Pax6 and mutated Pax6. Pax6, but not its mutated form (Pax6-mut), inhibited the expression of pre-miR-124-3 in a concentration-dependent manner (Fig. 6D). In agreement with the ChIP and EMSA data, we found that the inhibitory effect of PAX6 on miR-124-3 promoter activity disappeared when the PAX6 binding site was deleted in the miR-124-3 promoter (miR-124-3-pro-del) (Fig. 6D). We therefore conclude that PAX6 directly represses the expression of pre-miR-124-3, and PAX6 deficiency maintains pre-miR-124-3 expression at a high level.
Regulation of EB cell migration by PAX6 is mediated by miR-124
To evaluate the causal role of PAX6-miR-124 in the regulation of EB cell migration, we knocked down Pax6 expression in Wt-EBs or alternatively overexpressed Pax6 in Homo-EBs to determine whether PAX6 expression rescued cell migration in EBs. Consistent with our speculation, the migration of Wt-EBs after Pax6 knockdown was inhibited similarly to that of Homo-EBs (Fig. 7A). In contrast, migration of Homo-EBs after the Pax6 overexpression was enhanced to the Wt-EB level (Fig. 7B). Moreover, we found that the migration of cells from Wt-EBs was blocked when miR-124 mimic was introduced as determined by EB morphology (Fig. 7C, upper lane) and the transwell migration assay (Fig. 7C, lower lane). Figure 7D shows the elevated expression of miR-124 in EBs after transfection of the miR-124 mimic and the statistical results of the transwell migration assay. In contrast, inhibition of miR-124 by miR-124 inhibitor rescued migration in Homo-EBs (Fig. 7E, F). All of these data indicate that the cell migration in EBs regulated by PAX6 is mediated by miR-124 expression.

PAX6 regulates cell migration during iPSC differentiation through miR-124.
Discussion
In this study, using qPCR, in situ hybridization, and Northern blotting analysis, we investigated the expression pattern of miR-124 during EB formation of iPSCs derived from wild-type and Pax6 mutant mice. We show for the first time that miR-124, which is an important regulator of cell migration during gastrulation and neurogenesis, is a direct downstream target of PAX6. Moreover, using gain-of-function and loss-of-function studies, we investigated the role of miR-124 expression in mediating the regulatory effect of PAX6 on cell migration. In the iPS stage, no significant variation of migration-related gene expression was observed from the gene array analysis, but in the EB stage, cell migration ability was impaired markedly in Homo-EBs and migration-related genes, such as IQGAP1, were prominently downregulated; in contrast, PAX6 expression was elevated in Wt-EBs.
Molecular mechanisms that either maintain stemness or precede differentiation of pluripotent stem cells in the early stages of embryonic development remain poorly understood. PAX6 is a key factor in regulating the development of the eyes, CNS, and pancreas [20,25,46,47]. In this study, we generated wild-type iPSCs as well as heterozygous and homozygous Pax6 mutant iPSCs to explore the effects of PAX6 deficiency. Using the formation of cell aggregates in EBs of iPSCs, we investigated cellular behavior during development under PAX6 deficiency. The observation of blocked cell migration in PAX6 deficiency is consistent with other studies, indicating that PAX6 plays important roles in cell migration [25] and neuronal differentiation during development [48,49]. The various cell morphological changes involved in EB formation are considered to recapitulate the course of ES cell differentiation [17,18]. Although previous reports showed that miR-124 is expressed in neural lineages, possibly regulating neurogenesis [10 –13], and is a key factor for migratory cell fate during EB formation [14], the mechanism by which miR-124 itself is regulated has not yet been determined clearly. This study provides a tool to characterize the expression of miR-124 regulated by PAX6, which controls cell motility, an essential feature of tissue and organ formation.
We confirmed that miR-124 is downregulated during formation of EBs from ES cells, as previously reported [14], but in Pax6 Homo-EBs, miR-124 was maintained at a higher level similar to that in ES cells (Fig. 3C, D). At the same time, the cell migration capacity of cells in Homo-EBs was obviously lower than that of cells in Wt-EBs. Several other studies have reported that miR-124 is involved in cell migration during neurogenesis [10 –13] and gastrulation [14]. We conclude that the regulation of cell migration by PAX6 is probably mediated through effects on the expression of miR-124. Most reported miR-124 targets, including Scp1 [10], Ptbp1 [12], Lamc1 [11], and Iqgap1, which is the target we evaluated here, show reciprocal protein expression with miR-124, which suggests that miR-124 keeps these genes silenced in Homo-EBs. Cells are released from tight adherent junctions through diminished expression of IQGAP1, which subsequently promotes the reorganization of the cytoskeleton for cell polarity and directional cell migration. Thus, PAX6 indirectly regulates downstream genes involved in iqgap1-mediated cell migration, as observed in this study.
Because PAX6 is important for the eye (ectoderm) and pancreas (endoderm) development and brain (ectoderm) patterning in animal studies [41], it remains to be determined whether the progression of differentiation is involved in the migration difference between Wt-EBs and Homo-EBs at day 6 of EB formation. It was reported that compared with transfecting iPSCs prior to EB formation with miRNAs, transfecting EBs can reduce the effects of transfection on the progression of differentiation, because only the cells in the outer layer of EBs were transfected by this way [50]. We transfected 2 (or 3)-day-old EBs with miRNA mimic and inhibitor to minimize the possibility that migration difference due to the progression of differentiation. Consistent with previous reports [51 –53], we found that defective PAX6 did not impair the spontaneous differentiation of EBs at day 6 as determined by the expression of markers of the three germ layers (Supplementary Fig. S3). The majority of migrating cells were nestin-positive neural stem cells, and the expression of nestin was not significantly different between Wt- and Homo-EBs (Supplementary Fig. S3D, E), which confirms the observation of a previous study that PAX6 is not necessary for mouse neuroectoderm specification at the early stage [51]. In addition, all-trans retinoic acid treatment induced expression of the mature neuronal marker βIII-tubulin in both Wt- and Homo-EBs, and the further neuronal differentiation did not increase the migration difference between these two groups (Supplementary Fig. S4). These results strongly suggest that progression of neural differentiation did not contribute to the migration difference. A full characterization of three germ layer differentiation of Pax6 mutant versus wild-type cells by microarray (Supplementary Fig. S3H) also indicated that Pax6 deficiency did not impair the progression of differentiation at early stage of EBs. Taking together these results suggest that Pax6 mutation blocked the migration of the cells at the early stage of spontaneous differentiation. We speculated that impaired migration would reflect in the late stage of differentiation, because Pax6 mutant mice have abnormal brain patterns and a disorganized pancreas. Together with PAX6 [23,54,55] and miR-124 [10 –12,44] are critical factors in neural (and pancreatic) differentiation, the Pax6-mir-124 pathway may be involved in mouse brain and pancreas development and the function needs to be studied in the future.
The expression of miR-124 is thought to be regulated by RE1-silencing transcription factor (REST) [56], which is a zinc-finger repressor that negatively regulates many neuronal genes in stem cells, progenitor cells, and non-neuronal cell types [57]. Similar to the previously observed upregulation of miR-124 in REST+/− ES cell lines [56,58], in the present study, miR-124 expression significantly increased in Pax6 Homo-EBs. In addition to miR-124, several other downstream targets have been identified to be downregulated by PAX6. VEGFA, which is a growth factor, and MMP2, which is a cellular matrix protein involved in tumor formation and metastasis, were found to be repressed by PAX6 [54,59]. Olig2, which is a transcriptional factor critical for glial cell fate determination during forebrain development, was also found to be directly repressed by PAX6 [55]. Most of these targets encode proteins, but our results indicate that genes encoding miRNAs are novel candidates for target genes downregulated by PAX6. Although the mechanism is not yet clear, some studies have focused on many downstream genes inhibited by PAX6 using cDNA microarray analyses [31,60]. Our finding that miR-124 is downregulated by PAX6 provides a potential target for improving cell migration and differentiation during abnormal development caused by PAX6 deficiency.
Several miRNA polycistrons—those encoding the most abundant miRNAs in ES cells and that are silenced during early cellular differentiation [14,61,62]—are occupied at their promoters by Oct4, Sox2, Nanog, and Tcf3 [8]. The miRNA expression programs directly downstream of Oct4/Sox2/Nanog/Tcf3 could help to prime ES cells for rapid and efficient differentiation. Oct4/Sox2/Nanog/Tcf3 likely contribute to this process through their occupancy of the Let-7g promoter [8]. Similarly, we now report that the promoter of miR-124-3 is occupied by PAX6, which suggests that transcription factors are involved in neural (or pancreatic) differentiation by regulating miRNAs. Indeed, our findings regarding the PAX6-mir-124 pathway are supported by previous studies both in vivo [30,63] and in vitro [64]. Maiorano et al. [63] reported that miR-124 displays three distinct expression levels: low in the VZ [65,66], intermediate in the SVZ, and high in more marginal layers [63]. This pattern of expression is opposite to that of PAX6 in the cortex, suggesting that miR-124 is downregulated by PAX6 during neurogenesis. Sansom et al. [30] analyzed mRNA levels in the E12.5 Pax6 sey/sey cortex relative to those in wild-type littermates by microarray; Mir-124 was one of 938 genes showing reproducible, statistically significant changes in expression. Delaloy et al. [64] observed that the expression pattern of miR-124 is opposite to that of PAX6 during neural differentiation of human ESCs in vitro. All of these reports are consistent with our finding that PAX6 downregulates miR-124 expression during embryonic development.
In conclusion, in addition to supporting early miR-124 expression in EBs and early cortical development, our data also strongly suggest that PAX6 modulates cell migration by regulating miR-24a expression during ES cell differentiation. The current discovery of the relationship between PAX6 and miR-124 is important for functional studies on miRNAs and differential transcriptional factors to uncover the molecular mechanisms that either maintain stemness or precede differentiation of pluripotent stem cells. Further, in light of their key roles in the nervous system and pancreas, our finding may also give new insights into the organogenesis of these tissues.
Footnotes
Acknowledgments
This research was supported by the 973 Program (2011CB966203 and 2011CBA01102) and the International S&T Cooperation Program (2011DFA31040 and 2013DFG30680) from the Ministry of Science and Technology of China, the Strategic Research Program of the Chinese Academy of Science (XDA01040408), grants from Shanghai City Committee of Science and Technology (12DZ1910900 and 12DZ1910800), and grants from the National Natural Science Foundation of China (31171417) and Beijing Natural Science Foundation (7142083). The authors also thank Junhua Zou of the Department of Medical Genetics for Peking University for G-banding chromosome analysis, and David Stenger of the Barbara Davis Center for Childhood Diabetes, University of Colorado, Denver, for editorial comments on a late version of the article.
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
