|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
EMBRYONIC STEM CELLS |
aStem Cell and Developmental Biology, Genome Institute of Singapore, Singapore;
bBioinformatics Institute, Singapore;
cInformation and Mathematical Sciences, Genome Institute of Singapore, Singapore;
dISIS Pharmaceuticals, Carlsbad, California, USA;
eBioinformatics and Pattern Discovery Group, IBM Thomas J. Watson Research Center, Yorktown Heights, New York, USA;
fDepartment of Chemical Engineering, Massachusetts Institute of Technology, Boston, Massachusetts, USA;
gDepartment of Medicine, Harvard Medical School, Boston, Massachusetts, USA
Key Words. Embryonic stem cells • MicroRNAs • Cell differentiation • Retinoic acid
Correspondence: Correspondence: Bing Lim, M.D., Ph.D., Stem Cell and Developmental Biology, Genome Institute of Singapore, #02-01 Genome, Singapore 138672. Telephone: 65-6478-8186; Fax: 65-6478-9005; e-mail: limb1{at}gis.a-star.edu.sg
Received on April 23, 2007;
accepted for publication on September 28, 2007.
First published online in STEM CELLS EXPRESS October 4, 2007.
| ABSTRACT |
|---|
|
|
|---|
Disclosure of potential conflicts of interest is found at the end of this article.
| INTRODUCTION |
|---|
|
|
|---|
Retinoic acid (RA) promotes differentiation of ESCs to ectodermal lineages [7–9]. RA is a morphogen that is required during development of the central nervous system, lung, and kidneys and during proximodistal patterning and limb generation [10–13]. Dysfunction in RA synthesis or RA receptor activation can lead to embryonic lethality [14]. Neuroectodermal differentiation of mESCs can also be induced under conventional serum-free culture conditions with N2 and B27 supplements (N2B27) [15].
Recent studies in invertebrate model systems have identified lsy-6, the first miRNA found to play a role in neuronal patterning [16], and microRNA (miR)-9a, which ensures the generation of the precise number of neuronal precursor cells during development [17]. Related work has also indicated key roles for miRNAs during neural differentiation in vitro [18, 19] and in vertebrate central nervous system development [20–22]. In particular, miR-134 has recently been identified as a potential regulator of dendritic spine volume and synapse formation in mature rat hippocampal neurons in vitro through the localized repression of a protein kinase, LimK1 [23]. The mouse homologue of miR-134, which demonstrates conservation across rodents and primates, was originally identified by cloning from the mouse cortex [24] and is located in a large imprinted miRNA gene cluster at the mouse Dlk1-Gtl2 domain [3].
Here, we report on our findings that implicate miR-134 in a previously uncharacterized role in enhancing differentiation of mESCs to ectodermal lineages. Beginning with microarray analysis, miRNAs whose expression changed significantly during mESC differentiation were identified. Among the miRNAs that were altered during mESC differentiation, miR-134 was upregulated by RA and N2B27 treatment but not during embryoid body differentiation. Interestingly, elevated levels of miR-134 enhanced mRNA levels of Sox1, Nestin, and Neurogenin-2, markers associated with differentiation toward ectoderm [15, 25]. This effect was blocked using an miR-134 antagonist. To investigate potential mechanisms by which miR-134 affects mESCs, the miRNA target prediction algorithm, rna22 [26], was used to identify putative mRNA targets of miR-134. Endogenous protein levels of four of these target genes tested, including Nanog and LRH1, were reduced by elevating miR-134 levels without concomitant decreases of the target mRNA levels. Together, these data provide evidence that miR-134 enhances differentiation of mESCs to ectodermal lineages and establish a functional role for miR-134 in modulating mESC differentiation through its activity against genes that include Nanog and LRH1, both of which are known to be positive regulators of Oct4 and of mESC growth.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Pre-miR miRNA precursors (Pre-miRs), Anti-miR miRNA inhibitors (Anti-miRs), and scrambled RNA oligomer (Scr; AGACUAGCGGUAUCUUUAUCCC) were purchased from Ambion. These were transfected into mESCs at a final concentration of 100 nM using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, http://www.invitrogen.com) as per the manufacturer's instructions. For transcript and protein analysis, 1.5 µg of gene-specific short hairpin RNA (shRNA) plasmids, nonsilencing shRNA plasmid, or empty vector were transfected into mESCs in 12-well plates, and drug selection was performed with 1 µg/ml of puromycin (Sigma-Aldrich) 24 hours post-transfection for a period of 2 days.
RNA Extraction and Microarray Analysis
For microarray, Northern blot, and reverse transcription (RT)-polymerase chain reaction (PCR) analyses, total RNA was extracted from cells using Trizol reagent as per the manufacturer's instructions and subsequently column-purified with RNeasy kits (Qiagen, Hilden, Germany, http://www1.qiagen.com). miRNA microarray was performed using a service provider (LC Sciences, Houston, TX, http://www.lcsciences.com). A brief description is given in the supplemental online data. The total array profile has Minimum Information About a Microarray Experiment (MIAME) accession number E-MEXP-977 (http://www.ebi.ac.uk/arrayexpress). Arrays were performed in duplicate, and 134 microRNAs were selected as significant (
twofold change). Expression profiling of coding genes was carried out using Illumina MouseRef-8 BeadArrays as per the manufacturer's instructions (Illumina Inc., San Diego, http://www.illumina.com). Total chip data have been deposited for public access with GEO repository accession number GSE4522
[NCBI GEO]
. Arrays were performed in duplicate, and 1,204 genes were selected as significant (
twofold change).
Northern Blot Analysis
miRNAs were probed using
32P uridine triphosphate (UTP)-labeled RNA probes transcribed in vitro from a T7 promoter driving DNA oligomer templates. 5S ribosomal RNA (rRNA) was probed using
32P-end-labeled DNA probe. The specific activity of each probe was routinely 107 cpm/µg RNA/DNA. Ten micrograms of total RNA was size-fractionated by 15% Tris-borate-EDTA-urea polyacrylamide gel electrophoresis (PAGE) and electroblotted onto MagnaProbe nylon membrane (GE Osmonics Labstore, Minnetonka, MN, http://www.osmolabstore.com). The membrane was probed overnight at 42°C in 6x standard saline citrate (SSC), 0.2% SDS, 5x Denhardt's solution with 32P-labeled RNA probes (1 x 106 cpm/ml) and washed three times at 30°C in 6x SSC, 0.2% SDS. Complexes were detected by Phosphoimager (GE Healthcare, Little Chalfont, Buckinghamshire, U.K., http://www.gehealthcare.com).
RT-PCR
cDNA synthesis for miRNAs was performed with 50 ng of total RNA according to the manufacturer's instructions (Applied Biosystems, Foster City, CA, http://www.appliedbiosystems.com). Endogenous miRNA levels were measured with inventoried TaqMan probes, and PCR Master Mix (Applied Biosystems). miR-16 was used as an internal control [27]. cDNA synthesis for coding genes was performed with 1 µg of total RNA according to the manufacturer's instructions (Applied Biosystems). Endogenous mRNA levels of pluripotency and differentiation markers were measured with inventoried TaqMan probes and PCR Master Mix. Endogenous mRNA levels of LRH1, FADD, and G
o were measured using SYBR Green PCR Master Mix (Applied Biosystems). Primer sequences are given in the supplemental online data. β-Actin was used as an internal control. All amplicons were analyzed using ABI Prism 7900HT Sequence Detection System 2.2 software (Applied Biosystems).
Protein Extraction and Western Blot Analysis
Cell pellets were washed in chilled phosphate-buffered saline (PBS) and incubated for 20 minutes in ice cold lysis buffer containing freshly added protease inhibitors. Lysates were cleared by centrifugation at 4°C for 10 minutes at 12,100g, and protein concentrations were determined using Bradford dye (Bio-Rad, Hercules, CA, http://www.bio-rad.com). For Western blot analysis, 10 µg of total protein was size-fractionated by SDS-PAGE on 10% bis-Tris acrylamide NuPAGE gels and transferred to Hybond-P polyvinylidene difluoride (PVDF) membrane (GE Healthcare) in 1x NuPAGE transfer buffer (Invitrogen) with 10% methanol. The membrane was probed with specific primary antibodies (Oct4, sc8628; Sox2, sc17320; bone morphogenetic protein [BMP] 4, sc6896; Gata4, sc1237; Nestin, sc21248; G
o, sc387; FADD, sc6036; Santa Cruz Biotechnology Inc., Santa Cruz, CA, http://www.scbt.com; LRH-1, F4100–44; United States Biological, Swampscott, Massachusetts, http://www.usbio.net; Nanog, AB5731; β-actin, MAB1501; Chemicon) and secondary horseradish peroxidase-conjugated antibodies (anti-goat horseradish peroxidase [HRP], sc2768; anti-rabbit HRP, sc2030; anti-mouse, sc2005; Santa Cruz Biotechnology). Antibody-protein complexes were identified by ECL-Plus (Amersham Biosciences, Piscataway, NJ, http://www.amersham.com) and film.
Colony Formation Assay
E14 cells mESCs were transfected with Pre-miRs (100 nM) and/or Anti-miRs (100 nM), scrambled RNA oligomer (100 nM), or Lipofectamine only (Mock Transfection) in 12-well culture plates at a density of 5 x 105 cells per well. These cells were subsequently trypsinized at 48 hours and resuspended in ESC medium. Various numbers of cells (200 and 400) were plated onto mouse feeder layers in six-well culture plates to form secondary ESC colonies. After 7 days, emerging colonies were stained using Wright-Giemsa stain (Sigma-Aldrich) and counted. Colony morphology and number provide an indication of the number of colony-forming undifferentiated mESCs present in a population of cells [28].
Immunostaining and Detection of Alkaline Phosphatase
Cells were fixed in 4% paraformaldehyde, washed twice with PBS, and then incubated for 5 minutes at –20°C in 95% ethanol (vol/vol in PBS). Cells were then washed three times with PBS, blocked for 1 hour in 5% normal goat serum in PBS with 0.1 x Triton X-100, and incubated overnight with anti-class III β-tubulin (Ab18660; Abcam, Cambridge, U.K., http://www.abcam.com) primary antibody at 4°C. Cells were next washed twice with PBS and incubated for 1 hour with the corresponding secondary antibody (goat anti-rabbit IgG 568; Molecular Probes, Eugene, OR, http://probes.invitrogen.com). Cells were washed and mounted. Immunofluorescence was observed using a Leica microscope (Heerbrugg, Switzerland, http://www.leica.com), and image analysis was conducted using IM50 software (Leica Microsystems, http://www.leica-microsystems.com/). The Alkaline Phosphatase Detection Kit (Chemicon) was used to determine alkaline phosphatase activity, according to the manufacturer's instructions.
Quantitative Analysis of Immunofluorescence
Image acquisition and subsequent neurite outgrowth measurements were performed using the Cellomics ArrayScan II High Content Screening platform (Swallowfield, U.K., http://www.cellomics.com). Two separate fluorescent filters at 20x objective magnification were used; channel 1 illuminated Hoechst 33342-labeled nuclei, allowing automated focusing upon the cells, and channel 2 excited objects stained with Alexa Fluor 568-labeled anti-class III β-tubulin. Both images were acquired through a 535 nm x 35 nm bandwidth dichroic emission filter. Exposure times for each wavelength were determined empirically. All images were then analyzed using Cellomics neurite outgrowth software, and data analysis for the number of neurons per total nucleus count in each field and average intensity of neuronal cell bodies and neurites are presented. Data shown are an average of values from 16 random fields from three replicate wells in two independent experiments.
Vector Construction
pOct4-Luciferase and pNanog-Luciferase were as described previously [29]. pLuc-microRNA response element (MRE) constructs were generated as described previously [26]. The pLuc-Nanog3'-untranslated region (UTR) construct was generated by cloning the entire Nanog 3'UTR into psiCHECK-2 (Promega, Madison, WI, http://www.promega.com) at the NotI restriction site. To generate the Nanog 3'UTR mutant construct, PCR-directed mutagenesis was performed with a pair of primers containing the mutant MRE sequence (GAGGACGTGTTAACTAGTTTCC). The PCR product was subsequently used as a template for full-length PCR. The Nanog 3'UTR mutant was also cloned into psiCHECK-2 at the NotI restriction site. For shRNA constructs, 19-base pair gene-specific regions for RNA interference (RNAi) were designed and cloned into pSUPER.puro (BglII and HindIII sites) as described previously [29]. All vector and oligonucleotide sequences are available on request.
Luciferase Assays
The pOct4/pNanog-Luciferase assays were performed as described previously [29]. A brief description is given in the supplemental online data. miRNA target validation assays were performed as described previously [26]. A brief description is given in the supplemental online data.
In Situ Hybridization
Whole-mount and sectioned in situ hybridization were performed as previously described [30, 31].
Statistical Analysis
Unless otherwise stated, the unpaired Student's t test was used to determine statistical significance. For microarray data, the Pearson's correlation was used for the similarity test. Dendrograms were drawn using Eisenplots (clustering and treeview tools), with arm lengths proportional to the respective correlation coefficients.
| RESULTS |
|---|
|
|
|---|
In particular, miR-134 expression increased after 4 days of RA treatment, with subsidence by day 6, but still at elevated levels above control (Fig. 1A–1C; supplemental online Fig. 1E). When an alternative neuroectodermal lineage differentiation promoting regime was used, the levels of endogenous miR-134 increased significantly within 2 days with N2B27 treatment (three times greater induction compared with RA), and again the levels of miR-134 subsided over the time course but remained elevated above control (Fig. 1B). Interestingly, miR-134 levels exhibited a small, but significant, decrease during EB formation (Fig. 1B).
|
miR-134 Modulates mESC Differentiation, Even in the Presence of LIF
In addition to the Oct4- and Nanog-promoter activity assays, colony-forming unit assays demonstrated that there was an
30% decrease in the ability of mESCs treated for 2 days with Pre-miR-134 to form mESC colonies compared with scrambled RNA oligomer-treated controls (Fig. 1F). This suggested that miR-134 alone could modulate mESC differentiation.
To examine whether miR-134 alone could induce changes in cell phenotype as determined by changes in the mRNA and protein levels of key marker genes associated with pluripotency and differentiation, quantitative PCR of total RNA samples taken from E14 and D3 mESCs 3 days post-transfection with Pre-miR-134 was performed. This showed that the Oct4 mRNA levels were significantly reduced, consistent with the reduction of promoter activity (Fig. 2A). Nanog mRNA levels were reduced significantly in D3, but not E14, mESCs transfected with Pre-miR-134. Sox2 mRNA levels were reduced significantly in both D3 and E14 mESCs after Pre-miR-134 transfection.
|
Using stringent thresholds for the rna22 algorithm [26] screening entire mRNA sequences, an MRE for miR-134 was predicted in the 3'UTR of Nanog (discussed below). No predicted MREs were found in the 5'UTR, coding region, or 3'UTR of Oct4 under the stringent thresholds used in the context of this paper. However, by relaxing these thresholds, rna22 predicts that there may be several MREs present in Oct4 mRNA. This is the subject of ongoing investigation.
In addition, elevated mRNA levels of early primitive ectoderm marker (Fgf5) and neuroectoderm marker (Nestin) were observed, with no change or a downregulation of endoderm (Gata4 and Sox17) and mesoderm (Bmp4 and Nkx2.5) marker transcripts (Fig. 2A). Conversely, transfection of other miRNAs (e.g., miR-22) did not elicit any substantial change in the level of these transcripts (data not shown). These findings indicate that elevated levels of miR-134 alone in undifferentiated mESCs promotes a transcriptional expression profile that is suggestive of differentiation toward ectoderm. Concomitant with changes in mRNA levels in E14 mESCs, Pre-miR-134 transfection led to increased protein levels of Nestin and decreased protein levels of BMP4. Gata4 protein levels were not affected (Fig. 2C).
miR-134 Induces Morphological Changes Indicative of mESC Differentiation
It was next determined whether miR-134 could elicit phenotypic effects. E14 mESCs transfected with water only (mock transfection), scrambled RNA oligomer (Scr), Pre-miR-Let-7i (a miRNA that did not promote mESC differentiation), Anti-miR-134, or Anti-miR-134+Pre-miR-134 did not display morphological changes for up to 4 days post-transfection (Fig. 2D, for Scr). The cells maintained the characteristic domed colony structures of mESCs. However, transfection with Pre-miR-134 alone induced visible morphological changes in these cells within 4 days post-transfection (Fig. 2D), with the cells acquiring a flattened epithelial-like morphology typical of differentiation. These data are consistent with reduced alkaline phosphatase activity in Pre-miR-134-treated mESCs (Fig. 2D), where alkaline phosphatase is expressed at higher levels in pluripotent ESCs compared with differentiated cells.
The mRNA Expression Patterns Between RA-Treated and miR-134-Transfected mESCs Demonstrate a High Degree of Correlation
The microarray profile of mESCs treated with Pre-miR-134 was compared with that of RA treatment alone to investigate the effects of miR-134 during RA-induced differentiation of mESCs. For RA-induced differentiation, total RNA was collected from untreated mESCs, and from mESCs after 2 and 4 days of RA treatment (in the absence of LIF), whereas for Pre-miR-induced differentiation, total RNA was collected 3 days after transfection of mESCs (in the presence of LIF) with Pre-miR-134, water only (mock transfection), scrambled oligo, or Pre-miR-Let-7i. RNA samples were labeled and hybridized to Illumina microarray chips, and heat maps were generated to compare the spectrum of transcript levels altered by RA or Pre-miR-134 treatment.
There was a high degree of correlation between the expression pattern of genes altered after 4 days of RA treatment in mESCs and all Pre-miR-134-responsive genes (Fig. 3A, 3B). Prior to Pearson correlation analysis, the intensity data were log2-transformed and subtracted from the medium intensity. The mean Pearson correlation coefficient between these two expression profiles (and replicates) was
0.57, whereas the mean Pearson correlation coefficient between Pre-miR-134 transfection/4 days of RA treatment and other treatments (ESC/mock transfection using water [MT]/Scr/Pre-Let-7i) was
–0.36. When comparing the means of the above Pearson correlation coefficients between RAd4 and miR-134 to the corresponding means between miR-134/RAd4 and ESC/MT/Scr/Pre-Let-7i, all p values were <.001. Analysis of mRNAs upregulated in response to either RA or Pre-miR-134 indicated increased levels of transcripts associated with ectodermal differentiation (Fig. 3C) [33–38], supporting the evidence presented in Figure 2. In contrast, levels of nonectoderm markers were not significantly altered by either RA or Pre-miR-134 [33–38]. Transfection of Pre-Let-7i, an miRNA that neither induces mESC differentiation nor is induced by RA, resulted in a very different pattern of gene expression compared with Pre-miR-134 treatment (Fig. 3A).
|
|
miR-134 Targets Nanog and LRH1, Among Other Genes
miRNAs play a role in all aspects of post-transcriptional regulation where they regulate proteins involved in splicing, translation inhibition, and mRNA degradation [39–41]. To begin to elucidate the potential role of miR-134, its predicted mRNA targets were ascertained. Using rna22 [26], 2,800 potential 3'UTR miR-134 targets were identified. Of these, 158 predicted targets were selected for validation in vitro with luciferase assays. For each target transcript, a single copy of the sequence segment that rna22 identified as a putative miR-134 binding site (MRE), was cloned within the 3'UTR of Renilla luciferase reporter mRNAs. As was reported previously [26], a total of 129 of 158 tested targets of miR-134 were validated using luciferase assays, thus suggesting that a high proportion (>80%) of the rna22-predicted targets may be post-transcriptionally regulated by miR-134.
Gene expression analysis revealed that 1,051 of the genes potentially targeted by miR-134 were expressed in mESCs (supplemental online Fig. 3A, 3B). Of these, approximately 50% were upregulated, 20% remained unchanged, and 30% were downregulated relative to water only (mock transfection) and scrambled oligomer control transfections. This suggests that miR-134 can have different effects, direct or indirect, on target mRNAs, where it may affect mRNA translation, degradation, and splicing. Interestingly, the expression profile of predicted target genes after Pre-miR-134 transfection was highly correlated with their corresponding expression profile after 4 days of RA treatment (supplemental online Fig. 3A, 3B). The mean Pearson correlation coefficient between these two expression profiles was 0.51, whereas the mean Pearson correlation coefficient between predicted target gene expression in Pre-miR-134 transfection or after 4 days of RA treatment and other treatments (ESC/MT/Scr/Pre-Let-7i) was
–0.39.
To begin to understand the role of miR-134 in the post-transcriptional regulation of its potential direct targets, four predicted mRNA targets were selected for further characterization. Upon Pre-miR-134 treatment, LRH1, G
o, and Nanog mRNA levels did not change, and FADD mRNA levels increased. The predicted MREs in the 3'UTRs of these four genes are depicted in supplemental online Fig. 3C. The orphan nuclear receptor liver receptor homologue 1 (LRH1, also known as NR5A2), is a transcription factor that binds to the proximal enhancer (PE) and proximal promoter (PP) regions of the Oct4 promoter (supplemental online Fig. 2A) and regulates its expression [42]; it is also known to modulate cell proliferation through its interaction with β-catenin [43]. LRH1 is an important modulator of stem cell fate and is required for maintaining Oct4 expression in the epiblast of the embryo; loss of LRH1 leads to early embryonic lethality [42]. G
o (guanine nucleotide binding protein,
o) and FADD (Fas-associated via death domain) have been shown to play key roles in neurite extension [44], apoptosis [45], and embryo development [46].
For these four targets, miR-134 elicited a significant decrease in luciferase-MRE activity (Fig. 5A). No effect was observed with Pre-miR-124a, Pre-miR-21, or Scr control, which were not predicted to target these MREs. Secreted frizzled receptor protein 2 (Sfrp2) is an example of a predicted MRE that was not affected by Pre-miR-134 (Fig. 5A). Luciferase activity of miR-134 reverse complement (RC) was suppressed by 92% with Pre-miR-134 transfection (Fig. 5A), confirming Pre-miR-134 activity in this assay. To ensure that miR-134's effects are not due to nonspecific random targeting of any MRE, the RC sequence of another miRNA, miR-21, was cloned into pLuc-MRE to act as a negative control. No effect of Pre-miR-134 against the miR-21 RC was observed. In addition, miR-134 could specifically suppress luciferase activity when the Nanog 3'UTR was cloned within the 3'UTR of the Renilla luciferase reporter (Fig. 5B). The effect of miR-134 on Nanog's 3'UTR was abolished when its predicted MRE was mutated while the integrity of the rest of the 3'UTR was maintained (Fig. 5B, 5C). Other miRNAs that were not predicted to target Nanog had no effect.
|
o, and Nanog
o (lanes 10 and 12), and Nanog (lanes 14 and 16) in both the presence and the absence of RA. RA also downregulated the endogenous protein levels of LRH1 (Fig. 5D, lanes 1 and 3), G
o (lanes 9 and 11), and Nanog (lanes 13 and 15), consistent with the induction of miR-134 by RA. Importantly, and for all three tested targets, protein downregulation was not associated with a concomitant downregulation in mRNA levels (Fig. 5F), suggesting that miR-134 acts, at least in part, by translationally inhibiting these mRNAs in E14 mESCs. An intriguing finding was that miR-134 on its own led to an increase (+50%) in FADD mRNA levels (Fig. 5F), whereas FADD protein levels decreased (–80%) (Fig. 5D, lanes 5 and 6), suggesting post-transcriptional suppression of FADD mRNA by miR-134 in the absence of RA. RA induced a significant increase in FADD mRNA (+100%) and enhanced the upregulation of FADD mRNA by Pre-miR-134 (+250%) (Fig. 5F, columns 7 and 8). However, in the presence of RA, miR-134 lost its ability to suppress FADD protein levels as both FADD mRNA and protein levels increased (Fig. 5D, lanes 7 and 8). This highlights that miR-134 on its own can affect mRNA levels and that its effect on protein levels can be context-dependent. RA- or N2B27-mediated changes on other components of the post-transcriptional machinery [47] may alter the mechanism of action by which miR-134 affects FADD translation.
Knockdown of LRH1, G
o, and Nanog Results in Differentiation of mESCs
To examine whether modulation of individual miRNA-targeted genes could directly promote mESC differentiation, FADD, LRH1, G
o, and Nanog levels were perturbed by RNAi. RNAi was effective in significantly reducing transcript levels of LRH1, FADD, G
o, and Nanog (Fig. 6A). Knockdown of Nanog, LRH1, and G
o resulted in differentiation of mESCs as determined by the reduction in Oct4 promoter activity (Fig. 6B). To compare the effects of miR-134 and target gene RNAi on mESC differentiation, transcripts associated with pluripotency (Oct4 and Sox2) and differentiation (Nestin, Otx2, Fgf5, Sox17, and Bmp4; Fig. 6C) were examined. RNAi against LRH1, G
o, and Nanog decreased Oct4 and Sox2, with concomitant elevation of ectoderm-associated Fgf5, but had variable effects on Nestin, Otx2, Sox17, and Bmp4 (Fig. 6C).
|
|
| DISCUSSION |
|---|
|
|
|---|
1,051 mESC-expressed predicted target mRNAs, Nanog, LRH1, and G
o, through heteroduplex formation with 3'UTR cis-elements. Interestingly, our data showed that miR-134 was not upregulated during differentiation of EBs. As EB formation results in a random differentiation of ESCs to endodermal, mesodermal, and ectodermal lineages, the heterogeneity of cell types present in EBs may explain why we do not see an overall increase in miR-134 levels, even though its expression may increase in the population of cells differentiating toward ectoderm. The data suggest that upregulation of miR-134 and its effects on target mRNAs are not absolute requirements for mESC differentiation per se but highlight its role in promoting ectodermal lineage differentiation that we observed (Figs. 2, 4). Indeed, Anti-miR-134 was able to attenuate the RA-induced changes in Sox1, Nestin, and Neurogenin2 mRNA levels but could not attenuate the Oct4 mRNA downregulation. Again, this suggests that miR-134 elevation by RA and N2B27 is not an absolute requirement for differentiation, but it enhances the RA- and N2B27-induced differentiation and aids the specification to ectodermal lineages.
Expression profiling demonstrated that miR-134 levels were notably higher in the adult brain and spinal cord relative to other tissues. The peak of miR-134 expression during embryogenesis between E13.5 and E17.5 correlates with the progressive expansion and stratification of the neuroepithelium (a result of neuronal proliferation, migration and differentiation) and a transition from neuro- to gliogenesis. We also observed a dichotomy between the expression of miR-124a and miR-134 during nervous system development, where miR-134 is higher relative to miR-124a in embryonic, but not adult, tissues. miR-124a expression is correlated with the transition from non-nervous system to nervous system-specific alternative splicing patterns by inhibiting the splicing suppressor PTBP1 [39]. PTBP1 controls the splicing of the pre-mRNA of its homolog PTBP2 and vice versa [39]. Interestingly, PTBP2 is a predicted target of miR-134 (supplemental online Fig. 3C). One might speculate that the relative levels of miR-124a and miR-134 expression may alter specific splicing patterns in a developmental stage-specific manner. In any case, these data indicate a developmental role for miR-134 prior to its role in dendritic spine formation in terminally differentiated neurons in vitro [23].
It has been proposed that miRNAs function in vertebrates by subtly modulating cell types [48], and most of our augmentation and reduction of miRNA expression experiments did not show any significant phenotype change in mESCs, apart from experiments with miR-134. This may be linked to the number of potential mRNAs with which an miRNA can form heteroduplexes. The algorithm for calculating miRNA-mRNA interactions, rna22, predicted 2,800 mRNA 3'UTR cis-elements in the mouse genome as targets for forming heteroduplexes with miR-134, of which we confirmed 129 of 158 as positive (>80%) [26]. More than 1,000 genes expressed in mESCs comprised miR-134 potential targets. This implies that the mechanism by which miR-134 post-transcriptionally regulates mRNA to enhance mESC differentiation is complex. Interestingly, RNAi of any single target of miR-134 was not sufficient to replicate the Pre-miR-134 induced changes in mESC transcript levels, again suggestive of a wide-ranging miR-134 target network. Together, these data indicate that RA-induced elevation of miR-134, or exogenously elevated levels of miR-134 alone, may promote mESC differentiation through the coordinate regulation of a potentially large target-gene pool, which includes the pluripotency-associated target genes, such as Nanog, LRH1, and G
o. Interestingly, the sequences of the predicted miR-134 target sites in LRH1 and G
o are conserved in mouse, rat, and human and in mouse, rat, human, dog, and chicken, respectively.
As stated previously, miR-134 suppresses the protein levels of Nanog and LRH1. Nanog and LRH1 are known regulators of mESC self-renewal, where their interaction with specific promoter regions is integral to the transcriptional regulation of many genes, including Oct4 [32]. miR-134 overexpression resulted in the reduction of three other endogenous proteins tested to date, in part, through translation repression, an effect that was selectively blocked by Anti-miR-134. These genes, Wnt15, G
o, and FADD, have all been shown to play key roles in development [49–51]. Although Nanog mRNA levels remained unchanged in E14 mESCs during Pre-miR-134 treatment over 3 days but decreased in D3 mESCs, this may reflect subtle differences across cell types or that D3 mESCs may be more sensitive to elevated levels of Pre-miR-134 than E14 ESCs. Intriguingly, we observed that miR-134 on its own led to an increase in FADD mRNA levels, whereas FADD protein levels decreased. RA induced a significant increase in FADD mRNA and enhanced the effect of Pre-miR miRNA precursor 134, with concomitant increase in FADD protein. This suggests that the suppressive effect of miR-134 on FADD mRNA translation is lost in the presence of RA, and this highlights that the effect of an miRNA on a target or nontarget mRNA may be context-dependent. Our array profile data show that there was upregulation, no change, or downregulation of target and nontarget mRNAs, which again highlights that miR-134 may act through multiple gene-regulatory mechanisms.
In conclusion, elevated levels of miR-134 alone can enhance differentiation of mESCs, where this promiscuous miRNA exerts its effects through post-transcriptional regulation of multiple mRNAs, including translation attenuation of Nanog and LRH1. Elevation of miR-134 enhances the ectodermal differentiation of mESCs treated with RA or N2B27, where it augments the RA- and N2B27-induced elevation of ectodermal markers (β-III-tubulin, Nestin, Neurogenin2, and Sox1). These findings will aid in our understanding of miRNA function in embryonic stem cells and facilitate our ability to inhibit or enhance the differentiation of ectodermal lineages.
| DISCLOSURE OF POTENTIAL CONFLICTS OF INTEREST |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
I. Rigoutsos New Tricks for Animal MicroRNAs: Targeting of Amino Acid Coding Regions at Conserved and Nonconserved Sites Cancer Res., April 15, 2009; 69(8): 3245 - 3248. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T.N. Le, C. Teh, N. Shyh-Chang, H. Xie, B. Zhou, V. Korzh, H. F. Lodish, and B. Lim MicroRNA-125b is a novel negative regulator of p53 Genes & Dev., April 1, 2009; 23(7): 862 - 876. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Chen, F. Fang, Y.-C. Liou, and H.-H. Ng Zfp143 Regulates Nanog Through Modulation of Oct4 Binding Stem Cells, November 1, 2008; 26(11): 2759 - 2767. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Song and L. Wang Transcriptional mechanism for the paired miR-433 and miR-127 genes by nuclear receptors SHP and ERR{gamma} Nucleic Acids Res., October 1, 2008; 36(18): 5727 - 5735. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Yu and N. B. Hecht The DNA/RNA-Binding Protein, Translin, Binds microRNA122a and Increases Its In Vivo Stability J Androl, September 1, 2008; 29(5): 572 - 579. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Thum, D. Catalucci, and J. Bauersachs MicroRNAs: novel regulators in cardiac development and disease Cardiovasc Res, September 1, 2008; 79(4): 562 - 570. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Chen and G. Q. Daley Molecular basis of pluripotency Hum. Mol. Genet., April 15, 2008; 17(R1): R23 - R27. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. W. Yeo, N. Coufal, S. Aigner, B. Winner, J. A. Scolnick, M. C.N. Marchetto, A. R. Muotri, C. Carson, and F. H. Gage Multiple layers of molecular controls modulate self-renewal and neuronal lineage specification of embryonic stem cells Hum. Mol. Genet., April 15, 2008; 17(R1): R67 - R75. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| STEM CELLS | THE ONCOLOGIST | CME | ALPHAMED PRESS JOURNALS |