Next Article in Journal
Transcriptome Profiling Unravels a Vital Role of Pectin and Pectinase in Anther Dehiscence in Chrysanthemum
Next Article in Special Issue
Correction: Zhang, Z.; et al. Roles of MicroRNAs in Establishing and Modulating Stem Cell Potential. Int. J. Mol. Sci. 2019, 20, 3643
Previous Article in Journal
Genome-Wide Analysis of the YABBY Transcription Factor Family in Pineapple and Functional Identification of AcYABBY4 Involvement in Salt Stress
Previous Article in Special Issue
RNAase III-Type Enzyme Dicer Regulates Mitochondrial Fatty Acid Oxidative Metabolism in Cardiac Mesenchymal Stem Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 23
10.3390/ijms20235864
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

microRNAs Regulating Human and Mouse Naïve Pluripotency

by
Yuliang Wang
1,2,*,
Abdiasis M. Hussein
2,3,
Logeshwaran Somasundaram
2,3,
Rithika Sankar
2,3,
Damien Detraux
2,3,4,
Julie Mathieu
2,4,* and
Hannele Ruohola-Baker
2,3,*
1
Paul G. Allen School of Computer Science & Engineering, University of Washington, Seattle, WA 98195, USA
2
Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA 98109, USA
3
Department of Biochemistry, University of Washington, Seattle, WA 98195, USA
4
Department of Comparative Medicine, University of Washington, Seattle, WA 98195, USA
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(23), 5864; https://doi.org/10.3390/ijms20235864
Submission received: 10 September 2019 / Revised: 7 November 2019 / Accepted: 13 November 2019 / Published: 22 November 2019
(This article belongs to the Special Issue MicroRNAs and Stem Cells)
Browse Figures
Review Reports Versions Notes

Abstract

:
microRNAs are ~22bp nucleotide non-coding RNAs that play important roles in the post-transcriptional regulation of gene expression. Many studies have established that microRNAs are important for cell fate choices, including the naïve to primed pluripotency state transitions, and their intermediate state, the developmentally suspended diapause state in early development. However, the full extent of microRNAs associated with these stage transitions in human and mouse remain under-explored. By meta-analysis of microRNA-seq, RNA-seq, and metabolomics datasets from human and mouse, we found a set of microRNAs, and importantly, their experimentally validated target genes that show consistent changes in naïve to primed transitions (microRNA up, target genes down, or vice versa). The targets of these microRNAs regulate developmental pathways (e.g., the Hedgehog-pathway), primary cilium, and remodeling of metabolic processes (oxidative phosphorylation, fatty acid metabolism, and amino acid transport) during the transition. Importantly, we identified 115 microRNAs that significantly change in the same direction in naïve to primed transitions in both human and mouse, many of which are novel candidate regulators of pluripotency. Furthermore, we identified 38 microRNAs and 274 target genes that may be involved in diapause, where embryonic development is temporarily suspended prior to implantation to uterus. The upregulated target genes suggest that microRNAs activate stress response in the diapause stage. In conclusion, we provide a comprehensive resource of microRNAs and their target genes involved in naïve to primed transition and in the paused intermediate, the embryonic diapause stage.

1. Introduction

MicroRNAs are highly conserved small non coding RNA that post-transcriptionally regulate their target genes by binding mainly to the 3′UTR region of protein-coding mRNA, resulting in translational repression or transcript destabilization [1]. MicroRNAs have been shown to be essential for animal development [2,3], which is not surprising since they play key roles in various cellular processes such as cell cycle progression, proliferation, apoptosis, and differentiation [1,4,5]. MicroRNAs can regulate a large fraction of the transcriptome and function as rapid switches by regulating gene expression in a coordinated, fine-tuned manner. Indeed, a single microRNA can target hundreds of genes at the same time, and multiple microRNAs can target the same gene. MicroRNAs have therefore emerged as major regulators of cell fate determination, from maintenance of stem cell self-renewal to commitment toward specific lineages [6,7,8]. Disruption of enzymes involved in microRNA biogenesis (Drosha, DGCR8, Dicer) leads to the failure of embryonic stem cells (ESC) to downregulate pluripotency markers and to differentiate after exposure to a differentiation signals [9,10,11,12,13,14]. In addition, microRNAs are required for the generation of induced pluripotent stem cells (iPSC) [15] and regeneration in general, including heart regeneration and its opposing process, maturation [16,17,18,19,20].
We and others have shown that pluripotent stem cells express a very specific set of microRNAs [7,21,22,23,24], including miR-302 cluster, miR-372 (miR-290-295 cluster in mouse), the let-7 family and the primate specific imprinted CMC19 cluster containing the miR-520 family. Some of these microRNAs, named Embryonic Stem cell specific Cell Cycle regulating miRNA (ESCC miRNA) control the proliferation of ESC by repressing genes involved in cell cycle regulation [12,14,25]. In order to sustain their proliferation rate, pluripotent stem cells rely on glycolysis as their main source of energy [26,27,28]. The ESCC miRNA have also been shown to promote glycolysis in mESC by regulating the Mbd2-Myc-glycolytic enzymes axis [29].
Mouse and human ESC have been stabilized in culture at different pluripotency states corresponding to pre- and post-implantation stages in vivo and are referred as naïve and primed cells, respectively [30,31,32] (Figure 1A). Even though these cells are close in a developmental timeline, they are very different in terms of signaling requirements, gene expression, epigenetic landscape, and metabolic signature [26,30,31,32]. In the past few years it has become clear that pluripotency is a very dynamic stage and cells progress through a continuum of pluripotent states with unique properties for each state [30,33,34]. The pre-to-post-implantation transition can be suspended under certain conditions, and this stage is called diapause [35] (Figure 1A). Let-7 has been previously shown to be a potential regulator of diapause [36,37]. Additional microRNA regulators of diapause and their target genes remain under-explored.
Several microRNA-seq studies have revealed distinct microRNA signatures in naïve and primed mESC [38,39,40], and Dgcr8 KO experiments have shown that microRNA are essential for the transition from naïve mESC to primed mEpiSC [40]. In particular, the miR-302 cluster is expressed at higher levels in mEpiSC compared to mESC [22,38,40] and facilitates the exit of naive pluripotency in part by promoting the activity of MEK pathway [38]. To our knowledge, no study has compared the expression of microRNAs in naïve and primed human pluripotent stem cells. However, low concentration of the HDAC inhibitor sodium butyrate, which induces primed hESC to de-differentiate to an earlier stage in development [41], increases expression of miR-302 cluster while decreasing expression of miR-372 cluster [22], suggesting common microRNAs are involved in mouse and human naïve-to-primed transition. In this paper we compared naïve hESC grown in 2iLIF media [26,27,42,43,44] with primed H1 for their microRNA profile and analyzed it in parallel with their transcriptomic and metabolomic profiles. In addition, we combined existing datasets in mouse pluripotent cells [38,39,40] in order to find microRNAs regulating important pathways during the naïve to primed transition, and in naïve and primed states. We also identified 38 microRNAs as potential regulators of diapause by combining existing microRNA expression data [37] with our RNA-seq of diapause and post-implantation embryos [35].
We found 2184 consistent microRNA-target gene connections between 280 microRNAs and 647 target genes in human, and 435 consistent microRNA-target gene interactions between 80 microRNAs and 241 target genes in mouse. Importantly, we identified 115 microRNAs that significantly change in the same direction in naïve to primed transition in both human and mouse, many of which have not been previously reported, and serve as a resource for future studies. These microRNAs and their target genes regulate developmental (e.g., Hedgehog pathway) and metabolic pathways (e.g., fatty acid oxidation, OXPHOS) important for pluripotency. Interestingly, we found that microRNAs are likely to repress Sonic Hedgehog (shh) activity in human pluripotent cells. Indeed, microRNAs could down-regulate shh components in the naïve state. A negative regulator of shh pathway (GPR161) is upregulated in the primed state, since its regulator microRNA is reduced. These two miRNA based control systems keep shh activity low in both states despite the emergence of cilia at the post-implantation stage.

2. Results

2.1. microRNAs Regulating Human Naïve to Primed ESCs Transition

To dissect how microRNAs effect pre-to-post-implantation transition in early development (Figure 1A) we identified microRNAs and their experimentally validated target genes showing consistent expression changes between naïve and primed pluripotency using published datasets (Table 1). We first identified 357 microRNAs that are significantly differentially expressed between human naïve and primed ESCs (Table S1). We then identified 1146 protein-coding genes showing consistent differential expression in two naïve vs. primed studies [27,45]. We then used mirTarBase, a database of experimentally validated microRNA-target gene relations [46], to connect changes in microRNA and mRNA expression: a microRNA and its target gene is considered consistent if their expression levels change in the opposite direction in naïve vs. primed comparison. In total, we identified 2184 consistent microRNA-target gene connections between 280 out of the 357 differentially expressed microRNAs and 647 target genes out of the 1146 differentially expressed protein-coding genes (Table S2). Consistent microRNA-target connections are cases where a microRNA is up (down)-regulated, and its target gene is down (up)-regulated (Figure 1A,B). Of the 647 target genes, 243 are significantly higher in naïve hESC compared to primed hESC; 404 target genes are significantly lower in naïve hESCs. Although the microRNA-target gene interactions have prior experimental support from mirTarBase, and show consistent expression changes in naïve vs. primed comparison, our analysis is still based on expression correlation. Therefore, the functional relevance of any specific microRNA-target gene interactions from our comprehensive list in naïve to primed transition still requires experimental validation. However, our integrative analysis does provide a high confidence starting point for future studies.
To identify biological processes under microRNA regulation in pre-to-post-implantation development, we performed Gene Ontology enrichment for the up- and down-regulated microRNA target genes. In this analysis, we compared significantly down-regulated microRNA-target genes against all significantly down-regulated genes (as opposed to all genes, which is typically done). The purpose is to identify whether microRNAs regulate specific biological processes. Target genes up-regulated in naïve hESCs (i.e., down-regulated in primed hESCs) did not show significant enrichment in specific pathways. Interestingly, we found that target genes down-regulated in naïve hESCs compared to primed hESCs (the corresponding microRNAs up-regulated in naïve hESC) are significantly enriched for mRNA transcription and multiple developmental processes (e.g., NOTCH3, ETV5, GLI2, PTCH1; Figure 1C). In particular, GLI2, a key transcription factor in the sonic hedgehog (shh) pathway, is significantly down-regulated in naïve hESC compared to primed hESC, and its microRNA regulators, miR-654-5p, miR-92a-2-5p, and miR-541-3p are up-regulated in naïve hESC compared to primed hESC.

2.2. Hh Pathway in Naïve-to-Primed Transition

Interestingly, multiple components of the shh pathway are downregulated in naïve hESC compared to primed hESC, including GLI1, GLI3, PTCH1, PITCH2 and SMO (Figure 2B). The shh pathway is an important regulator of embryonic development (Figure 2A, [49]). However, despite the presence of functioning components in the primed state, the shh pathway seems to play only a minimal role in the maintenance of pluripotency [50]. Binding of Nanog to shh components has been proposed as one of the mechanisms to inhibit shh pathway in pluripotent cells [51]. Our data suggest that microRNAs could down-regulate the expression of the components of shh-signaling pathway in naïve hESC (Figure 2H). Furthermore, naïve, preimplantation mouse ESC do not contain cilia, the key mammalian cellular organelle for Hh signaling. However, this organelle is found in the post-implantation stage [52,53,54]. To test the presence of cilia in human naïve-to-primed hESC transition we stained naïve and primed hESC with ARL13B. Importantly, we observed a significant increase in cilia number between these two stages (20-fold increase between naïve and primed hESC; (Figure 2C,D). These data argue that cilia morphogenesis occurs in naïve-to-primed hESC transition, licensing primed hESC for Hh pathway activity. It is therefore perhaps even more surprising that the Hh pathway is not active in the primed hESC stage. Through a whole genome CRISPR screen we found that GPR161, a ciliary G-protein coupled receptor, is required for the naïve to primed transition in hESC [55]. GPR161 is up-regulated in primed hESC (Figure 2B) and some of the microRNA predicted to target GPR161 (miR-372, miR-520, and miR-22) are down-regulated in primed compared to naive hESC [22,27]. GPR161 can negatively regulate the shh pathway [56], hence it could explain why the Hh-pathway is not active in primed hESC despite the presence of cilia and up-regulation of Hh-pathway components (Figure 2H). To test this hypothesis, we generated a GPR161 mutant hESC line using CRISPR Cas9 (Figure 2E–G). Most of the cells from the GPR161 mutant line express a GPR161 protein with a truncated C terminus due to the insertion of one nucleotide, resulting in a frameshift at residue 392 and the introduction of a premature STOP codon in position 426 (Figure 2E). The missing region in the truncated mutant line has been shown to be essential for GPR161 function and regulation of PKA and shh [57,58]. The mutant GPR161 hESC line exhibits primary cilia (Figure 2F) and one of the shh components, SMO, is up-regulated in GPR161 mutant hESCs compared to controls (Figure 2G). These data suggest that lack of GPR161 may induce promiscuous Hh pathway activity in hESCs. We propose that GPR161 is required for repressing the Hh pathway in the primed hESC stage (Figure 2H).

2.3. Metabolism in Naïve-to-Primed Transition

Naïve and primed hESCs have distinct metabolic profiles (Table 2 and Table 3, [26]). Some of these metabolic differences have been shown to be regulated by microRNAs [27,29]. We previously found that naïve hESCs can utilize fatty acids as an energy source, while primed hESCs cannot [27] (Table 2 and Table 3). Our analysis confirmed higher expression of miR-10a in naïve compared to primed hESC. Accordingly, miR-10a target genes in fatty acid activation (ACSL1), synthesis (FASN), elongation (ELOVL7), and desaturation (SCD) are downregulated in naïve compared to primed hESCs, suggesting that fatty acid synthesis is repressed in naïve hESC by critical miRNAs (Table 3). In contrast, CPT1A, a rate limiting transporter in fatty acid beta-oxidation is highly upregulated in naïve compared to primed hESC [27]. We found 11 microRNAs targeting CPT1A were significantly downregulated in naïve compared to primed hESC (miR-106b-5p, miR-20a-5p, miR-302c-3p, miR-17-5p, miR-20b-5p, miR-106a-5p, miR-302a-3p, miR-150-5p, miR-16-5p, miR-93-5p, miR-4517; see Table S1 for more details), suggesting that CPT1 and thereby beta-oxidation is down-regulated in primed hESC by critical miRNAs (Table 3). These data support our previous observation that fatty acid oxidation is more active in naïve hESCs compared to primed hESCs (Table 2 and Table 3). Thus, microRNAs may regulate the fatty acid oxidation (higher in naïve) vs. synthesis (higher in primed) differences between naive and primed hESCs (Table 3).
Additionally, naïve hESCs have higher oxidative phosphorylation (OXPHOS) activity than primed hESCs [27,59]. COX7B, a protein that is indispensable for cytochrome c oxidase (complex IV) assembly, activity, and mitochondrial respiration [60], is expressed higher in naïve than primed hESCs. Interestingly, hsa-miR-18a-5p and miR-17-5p that target COX7B are upregulated in in primed hESCs, and might therefore downregulate COX7B at that stage. miR-18a-5p and miR-17-5p belong to the miR-17~92 microRNA Cluster, which is a global regulator of cancer metabolism, known to suppress the expression of LKB1 [61]. Another member of the cluster, miR-20a, also targets CPT1A, a key step in fatty acid oxidation, suggesting the miR-17~92 cluster’s general role in the metabolic switch between naïve and primed hESC stages.
We further integrated microRNA expression, mRNA expression, and metabolomics data, and found consistent changes across three datasets in the polyamine metabolic pathway (Figure 3). ODC1 (ornithine decarboxylase) is a rate-limiting step in polyamine metabolism, and is significantly up-regulated in naïve hESCs. miR-218-5p, which targets ODC1, is down-regulated in naïve hESCs. Additionally, the substrate of ODC1, ornithine, is significantly depleted in naïve hESCs based on our metabolomics data. Spermidine, a key product of polyamine synthesis, is more abundant in naïve hESCs. Consistent with the conversion of ornithine to spermidine, the alternative route of ornithine catabolism, citrulline, is less abundant in naïve hESCs. Other enzymes or metabolites in the polyamine pathway did not change, possibly because some enzymes are regulated at the translational level, not transcriptional level (e.g., AMD1) [67]. ODC1 has been shown to be essential for mouse ESC self-renewal [68], and high polyamine (e.g., spermidine) levels have been shown to promote mouse ESC self-renewal [69]. Our analysis suggests that ODC1 and spermidine may also promote human naïve pluripotency, and down-regulation of miR-218-5p may be important for maintaining high polyamine metabolism in naïve hESCs. miR-218-5p has been widely studied in cancer [70].

2.4. microRNAs Regulating Mouse Naïve to Primed ESCs Transition

microRNAs also play important roles in regulating pluripotency in mouse. Several studies have profiled microRNAs in mouse embryonic stem cells and epiblast stem cells, or during mouse embryonic stem cell differentiation [38,39,40] (Table 1). To generate a robust consensus set of microRNAs associated with mouse naïve to primed transition, we integrated microRNA-seq datasets from three different studies. We identified 221 microRNAs that show consistent and significant changes in at least two out of the three studies (same direction and at least 1.5-fold change for either the -3p or -5p form, Table S3 and Methods). Similar to the human microRNA analysis, we combined these consistently changing microRNAs with mirTarBase information, and mESC vs. mEpiSC mRNA-seq data, and found 836 consistent microRNA-target gene interactions between 119 microRNAs and 304 target genes (Figure 4A, only a subset of the 221 microRNAs are connected to target genes changing in the opposite direction. The full list of 836 interactions are in Table S4).
Through Gene Ontology enrichment analysis, we found that compared to all down-regulated genes, down-regulated genes that are also targets of these consistently changing microRNAs are significantly enriched for multiple developmental processes (Figure 4B). This is similar to human naïve to primed transition, where microRNAs also target developmental genes (Figure 1B).
MicroRNAs also regulate the metabolic differences between naïve and primed states in mouse, especially in amino acid metabolism based on our analysis. We identified two metabolic genes, Slc1a2 and Slc25a12, that satisfy the following three criteria: (1) significant transcriptomic change; (2) their microRNA regulators change in the opposite direction; and (3) their connected metabolite (e.g., substrate/product) show significant change between naïve and primed mESCs. Slc1a2, a glutamate transporter, is 11.7-fold higher in mEpiSC compared to mESC, while its microRNA regulator, miR-199b-3p, is consistently up-regulated in all naïve state (down-regulated in mEpiSC) in three microRNA-seq studies. Slc1a2 is important for glutamate homeostasis [71]. Consistently, glutamate level is 1.93-fold higher in mEpiSC compared to mESC based on our metabolomics data. Glutamate can be converted to α-ketoglutarate and utilized in the TCA cycle. Higher expression of the glutamate transporter and abundance of glutamate in EpiSC suggest that glutamate may be an important fuel for EpiSC. Supporting this hypothesis, uptake of glutamate and its conversion to α-ketoglutarate has been previously reported during the initial differentiation of human primed pluripotent stem cells [72].
Slc25a12 (mitochondrial aspartate/glutamate carrier, commonly called aralar) is 5.5-fold higher in mESC compared to EpiSC. Its microRNA regulators (miR-466h-5p, 466m-5p, and 669m-5p) are all lower in mESC. Slc25a12 is a critical component of the mitochondria malate-aspartate shuttle, and deficient Slc25a12 expression is known to increase lactate production from pyruvate, and reduce glucose oxidation in the brain [73]—a metabolic phenotype similar to the primed state. Additionally, deficient Slc25a12 expression results in reduced aspartate level in the brain [74], and our metabolomics data showed that aspartate is significantly lower in EpiSC (where Slc25a12 expression is lower) compared to mESC. The expression pattern of Slc25a12 and abundance of aspartate may influence the shift from glucose oxidation to glycolysis and lactate production in the naïve to primed transition.

2.5. Consistent microRNA Changes across Human and Mouse Pluripotency

Previous studies have shown that certain microRNAs such as let-7 and microRNAs in the imprinted Dlk1-Dio3 locus are important regulators of both human and mouse pluripotency [10,39]. To systematically identify more shared microRNAs regulating both human and mouse pluripotency, we focused on 115 microRNAs (Table S5) that satisfy two criteria: (1) that show significant differences between human naïve and primed ESCs; and (2) that change in the same direction in at least one mouse naïve vs. primed microRNA-seq studies (>1.5-fold). We recovered members of the let-7 family (let7-b, d, f, and i) and microRNAs from the Dlk1-Dio3 locus (miR-541-5p, miR-410-3p, miR-381-3p, and miR-495-3p) as consistently higher in human and mouse naïve ESCs. We also found miR-143-3p to be the only microRNA that is 1.5-fold higher in naïve ESCs across all four different datasets. miR-143 is significantly induced by vitamin C and vitamin C is known to promote naïve pluripotency [75]. Additionally, miR-143 over-expression reduced differentiation and up-regulates key pluripotency genes (Oct4, Klf4, and Esrrb) [76]. Thus, this list of 115 consistently changing microRNAs provide a useful resource to mine microRNAs regulating pluripotency state transitions. We further identified microRNAs showing consistent changes in three or more pluripotency studies (Table 4).
To reveal the potential functions of these 115 microRNAs, we used Co-expression Meta-analysis of miRNA Targets (CoMeTa) [77] to assign enriched Gene Ontology terms to 30 of the 115 microRNAs (Table S5). These inferred functions include small GTPase signaling (miR-582, higher in naïve), angiogenesis (miR-369, higher in naïve), cell death (miR-338-3p, higher in primed), and cell adhesion (miR-671-5p, higher in primed). Previous studies have found that low cell-matrix adhesion promotes homogeneous self-renewal of naïve mESCs [78], which supports the computationally predicted function of miR-671-5p.
Of the 115 microRNAs, 73 have at least one target gene that changes in the same direction between human and mouse. That is, both the microRNA and their target genes change in the same direction in both species (Table S6). Some of the consistently changing microRNAs regulate the metabolic transitions between naïve and primed pluripotency. For example, miR-615-3p is higher in the naïve state, and its target gene ATP13A2 is lower in both human and mouse naïve stem cells. ATP13A2 is important for mitochondria maintenance, and ATP13A2 knockdown cells show increased oxygen consumption [79]. The increased oxygen consumption in ATP13A2 knockdown cells is similar to naïve stem cells where ATP13A2 expression is lower, and oxygen consumption is higher compared to primed stem cells [27]. Similarly, miR-485-5p is higher, while its target gene HIF3A is lower in both human and mouse naïve state. HIF3A expression is known to increase in response to hypoxia in pluripotent stem cells [80], and its lower expression in the naïve state agrees with our previous results that the HIF pathway is activated in the primed state but not in naïve [27]. Additionally, VLDLR (very low density lipoprotein receptor) is higher in primed (lower in naïve), and three microRNAs targeting VLDLR, miR-493-3p, 376c-3p, and 409-5p are all lower in primed (higher in naïve). Higher VLDLR expression is associated with increased amount of lipid droplet formation, which is consistent with our previous observation [27]. Six microRNAs (miR-302c-3p, miR-20b-5p, miR-106a-5p, miR-302a-3p, miR-455-3p, and miR-338-3p) targeting ACOT9 (Acyl-CoA thioesterase 9) are all lower in naïve; and ACOT9 is higher in naïve. ACOT9 hydrolyzes both short- and long-chain acyl-CoAs, and links amino acid and fatty acid metabolism in the mitochondria [81].

2.6. microRNAs and Their Target Genes Involved in Diapause Regulation

The transition from pre-to-post-implantation can be suspended in an intermediate state called diapause. To obtain a more complete picture of microRNAs potentially involved in regulating the diapause state, we combined a publicly available microarray-based profiling of diapause and re-activated embryos [37] with our RNA-seq of diapause and post-implantation embryos [35], using the strategy outlined in Figure 1A. We identified 379 consistent connections between 38 microRNA and 274 target genes, where a microRNA is up (down)-regulated, and its target gene is down (up)-regulated (Table S7). We compared target genes up-regulated in diapause (their microRNA regulators down in diapause) against all up-regulated genes in diapause, and found that these microRNA target genes show more enrichment in stress response (e.g., EGR1) [82], hypoxia (BACH1) [83], and hormone response (e.g., ZBTB7A) [84] pathways, which is consistent with the diapause phenotype (Figure 5A). On the other hand, target genes down-regulated in diapause (their microRNAs up-regulated in diapause) showed enrichment in cellular morphogenesis, differentiation, and development (Figure 5B), which is consistent with diapause as a suspended state in developmental progression.

3. Materials and Methods

Raw mRNA-seq data from Sperber et al. [27] and Theunissen et al. [45] were aligned to hg19/GRCh37 with STAR aligner57. Transcript quantification was performed with htseq-count from HTSeq package58 using GENCODE v15. Differential expression analysis was performed with DESeq after filtering out genes whose total read count across samples are below the 40th quantile of all genes. RNA-seq of mouse ESC and EpiSC from Factor et al. [48] was analyzed similarly using the mouse genome reference. Differentially expressed microRNAs were downloaded from supplementary material from Sperber et al. [27], Jouneau et al. [38], Gu et al. [40], and Moradi et al. [39]. topGO R package [85] was used for Gene Ontology enrichment analysis. Processed metabolomics data were downloaded from Sperber et al. [27]. Single cell RNA-seq of monkey in vivo implantation (Figure 3C) was downloaded from Nakamura et al. [86]. We noticed that when previous studies established the role of microRNAs in pluripotency (e.g., miR-290), they did not distinguish miR-290-3p vs. miR-290-5p (Gu et al., Cell Research 26:350–366 (2016)). We therefore combined the -3p/-5p form when analyzing the three mouse microRNA-seq datasets. A microRNA is considered consistently changing if any one of the three criteria is satisfied: (1) the -3p form changed ≥1.5-fold in at least two out of three studies; (2) the -5p form changed ≥1.5-fold in at least two out of three studies; (3) the -3p form changed 1.5-fold in one study, the -5p form changed >1.5-fold in the same direction as the -3p form, but in a different study.

3.1. Culture of Naïve and Primed Human Pluripotent Stem Cells

Naïve hESC Elf-1 (NIH_hESC Registry #0156, [42]) and primed human iPSC (WTC-11) [87] were cultured as previously described [26,27,43,44]. Naïve cells were grown on a feeder layer of irradiated primary mouse embryonic fibroblasts (MEF) in naïve hESC 2iL-I-F media: DMEM/F-12 media supplemented with 20% knock-out serum replacer (KSR), 0.1 mM nonessential amino acids (NEAA), 1 mM sodium pyruvate, penicillin/streptomycin (all from Invitrogen, Carlsbad, CA, USA), 0.1 mM β-mercaptoethanol (Sigma-Aldrich, St. Louis, MO, USA), 1 µM GSK3 inhibitor (CHIR99021, Selleckchem, Houston, TX, USA), 1 µM of MEK inhibitor (PD0325901, Selleckchem), 10 ng/mL human LIF (Chemicon, Temecula, CA, USA), 5 ng/mL IGF1 (Peprotech, Rocky Hill, NJ, USA) and 10 ng/mL bFGF. Cells were passaged using 0.05% Trypsin-EDTA (Life Technologies, Carlsbad, CA, USA). Cells were transferred to matrigel-coated plates prior molecular analysis. Naive Elf1 cells were pushed to the primed state by culturing the cells in mTeSR1 media (StemCell Technologies, Vancouver, BC, Canada) on matrigel-coated plates. Primed pluripotent stem cells WTC11 were cultured in mTeSR1 media on Matrigel-coated plates. All cells were cultured at 37 °C, 5% CO2, and 5% O2.

3.2. Generation of GPR161 Mutant hESC

Naïve hESC (Elf1 2iL-I-F) were spin-infected with lentiCRISPR-v2 lentiviral prep expressing GPR161 sgRNA (CCCACACCTCACTGCGCTCA) in presence of 4 μg/mL polybrene and plated onto irradiated DR4 MEF. Two days later, cells were selected with puromycin (0.5 µg/mL, for 3 days). Genomic DNA was extracted using DNAzol reagent (Invitrogen) according to manufacturer’s instructions and quantified using Nanodrop ND-1000. Genomic regions flanking the CRISPR target sites were PCR amplified (primers Fwd: TCACCTTGGTGGTGGTTGAC; Rev: AAAACGCGACAGGTGAGAGG) and sent for Sanger sequencing. Analysis of the CRISPR editing data through the ICE tool software (ICE v2; Synthego, Sunnyvale, CA, USA) predicted a score of 73% KO in the pooled population.

3.3. Immunofluorescence Staining and Confocal Imaging

Naïve and primed pluripotent stem cells were cultured on Matrigel-coated glass cover slips at the bottom of the culture dishes. Cells are fixed with 4% paraformaldehyde (PFA) for 15 min, washed with PBS 1X (3 times, 5 min each) and blocked with 3% BSA (VWR 0332) and 0.1% Triton (Sigma 9002-93-1) (blocking solution) for 1 h at room temperature (RT). Cells were then incubated overnight with ARL13B- (Proteintech 17711-1-AP, 1:250), and ZO-1 (Invitrogen 33-9100, 1:250) antibodies diluted in blocking solution. After PBS 1X washes (3 times, 5 min each) cells were incubated with anti-rabbit Alexa 488 and anti-mouse 647 conjugated secondary antibodies (1:250) at RT for 1 h. The slides were then washed 4 times with PBS 1X and stained with DAPI DAPI during the 2nd wash (1:10) for 10 min. Glass cover slips were mounted on slides with vectashield Hardset (H-1400). The images were taken on Leica SPE confocal and Nikon A1 confocal and analyzed using Fiji ImageJ, v1.52g. The percentage of cells with the presence of a cilia was assessed by counting Arl-13B positive cells. At least 500 cells were counted for each condition.

3.4. RT-qPCR Analysis of GPR161 Mutant

RNA was extracted using Trizol (Life Technologies) according to manufacturer’s instructions. RNA samples were treated with Turbo DNase (ThermoFischer, Waltham, MA, USA) and quantified using Nanodrop ND-1000 (Thermo Scientific). Reverse transcription was performed using iScript reverse transcription kit (BioRAD, Hercules, CA, USA). 10 ng of cDNA was used to perform qRT-PCR using SYBR Green (Applied Biosystems, Foster City, CA, USA). Real-time RT-PCR analysis was performed on 7300 real time PCR system (Applied Biosystems) and β-actin was used as an endogenous control. Primer sequences for SMO are Fwd: 5′-ACGAGGACGTGGAGGGCTG-3′ and Rev: 5′-CGCACGGTATCGGTAGTTCT-3′; and for β-actin Fwd: TCCCTGGAGAAGAGCTACG and Rev: GTAGTTTCGTGGATGCCACA.

4. Conclusions

microRNAs play important roles in regulating cell fate choices. Through integrative analysis of mRNA-seq, microRNA-seq, and metabolomics data in human and mouse, we identified a comprehensive resource for microRNAs and their target genes involved in regulating the developmental and metabolic processes in naïve to primed transition, as well as in their paused, intermediate stage, diapause.
We identified a set of microRNAs with robust changes across multiple studies in both human and mouse by integrating transcriptomic data of mRNA and microRNA, as well as metabolomics data of naïve to primed transition in both human and mouse. The target genes of these consistently changing microRNAs regulate important metabolic and developmental pathways in naïve to primed transition. Importantly, we propose that the dynamics of the paused stage, embryonic diapause are affected by microRNA regulated stress and inflammation response genes. This is interesting in the light of recent publications that have revealed importance of stress response and inflammation in regulation of embryonic diapause and stem cell quiescence [35,88,89], and calls for further analysis on microRNA contribution in stress induced regulation of pause in development.
These studies also propose that microRNAs could play an important role in the maintenance of pluripotency. In particular, we suggest that miRNAs repress Hh pathway activity in both naïve and primed hESCs, preventing their differentiation. Interestingly we found in human ESC, as previously observed in mouse embryos, that primary cilia appear only at the primed stage of hESCs, in naïve hESCs the presence of primary cilia is minimal. Since primary cilium is essential for Hh-pathway activity, this could explain the lack of Hh signaling in the naïve pluripotent stage. Furthermore, microRNA may target components of the shh pathway at the naïve hESC stage. Our data suggest that microRNAs could target a repressor of shh, the G-protein coupled receptor GPR161, usually located in primary cilia [56]. As cells transition from the naïve to primed stage, primary cilia are formed and Hh-pathway components, for example GPR 161, PTCH, GLI, and SMO are expressed, however, the Hh-pathway is not activated. We now propose that in order to keep the Hh-pathway inactive in the primed hESC, the cells utilize GPR161. In the naïve to primed transition microRNAs targeting GPR161 are downregulated and GPR161 is expressed, resulting in the repression of the Hh pathway activity, despite the existing cilia and expression of the Hh pathway components. In summary, our study provides a valuable resource for dissecting microRNA based regulation of naïve to primed transition.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/20/23/5864/s1.

Author Contributions

Y.W., J.M. and H.R.-B. wrote the manuscript. Y.W. performed the bioinformatic analysis J.M. and H.R.-B. conceived and designed the experiments. A.M.H., L.S., R.S. and D.D. performed the experiments. All of the authors participated in the preparation of the figures.

Funding

This work was partially supported by a gift from the Hahn Family, and grants from the National Institutes of Health 1P01GM081619, R01GM097372, R01GM97372-03S1 and R01GM083867 and the NHLBI Progenitor Cell Biology Consortium (U01HL099997; UO1HL099993) for HRB, and the ISCRM Innovation Pilot Award for JM.

Acknowledgments

We thank members of the Ruohola-Baker laboratory for helpful discussions. We thank Sonia Sidhu and Christopher Cavanaugh for technical help, and Benjamin Freedman for providing the ARL13B and ZO-1 antibodies.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bartel, D.P. MicroRNAs: Target recognition and regulatory functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [PubMed]
  2. Alvarez-Garcia, I.; Miska, E.A. MicroRNA functions in animal development and human disease. Development 2005, 132, 4653–4662. [Google Scholar] [CrossRef]
  3. Gebert, L.F.R.; MacRae, I.J. Regulation of microRNA function in animals. Nat. Rev. Mol. Cell Biol. 2019, 20, 21–37. [Google Scholar] [CrossRef] [PubMed]
  4. Hatfield, S.D.; Shcherbata, H.R.; Fischer, K.A.; Nakahara, K.; Carthew, R.W.; Ruohola-Baker, H. Stem cell division is regulated by the microRNA pathway. Nature 2005, 435, 974–978. [Google Scholar] [CrossRef] [PubMed]
  5. Shcherbata, H.R.; Hatfield, S.; Ward, E.J.; Reynolds, S.; Fischer, K.A.; Ruohola-Baker, H. The MicroRNA pathway plays a regulatory role in stem cell division. Cell Cycle 2006, 5, 172–175. [Google Scholar] [CrossRef] [PubMed]
  6. Mathieu, J.; Ruohola-Baker, H. Regulation of stem cell populations by microRNAs. Adv. Exp. Med. Biol. 2013, 786, 329–351. [Google Scholar] [PubMed]
  7. Bar, M.; Wyman, S.K.; Fritz, B.R.; Qi, J.; Garg, K.S.; Parkin, R.K.; Kroh, E.M.; Bendoraite, A.; Mitchell, P.S.; Nelson, A.M.; et al. MicroRNA discovery and profiling in human embryonic stem cells by deep sequencing of small RNA libraries. Stem Cells 2008, 26, 2496–2505. [Google Scholar] [CrossRef] [PubMed]
  8. Kuppusamy, K.T.; Sperber, H.; Ruohola-Baker, H. MicroRNA regulation and role in stem cell maintenance, cardiac differentiation and hypertrophy. Curr. Mol. Med. 2013, 13, 757–764. [Google Scholar] [CrossRef]
  9. Murchison, E.P.; Partridge, J.F.; Tam, O.H.; Cheloufi, S.; Hannon, G.J. Characterization of Dicer-deficient murine embryonic stem cells. Proc. Natl. Acad. Sci. USA 2005, 102, 12135–12140. [Google Scholar] [CrossRef]
  10. Melton, C.; Judson, R.L.; Blelloch, R. Opposing microRNA families regulate self-renewal in mouse embryonic stem cells. Nature 2010, 463, 621–626. [Google Scholar] [CrossRef]
  11. Kanellopoulou, C.; Muljo, S.A.; Kung, A.L.; Ganesan, S.; Drapkin, R.; Jenuwein, T.; Livingston, D.M.; Rajewsky, K. Dicer-deficient mouse embryonic stem cells are defective in differentiation and centromeric silencing. Genes Dev. 2005, 19, 489–501. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Y.; Baskerville, S.; Shenoy, A.; Babiarz, J.E.; Baehner, L.; Blelloch, R. Embryonic stem cell-specific microRNAs regulate the G1-S transition and promote rapid proliferation. Nat. Genet. 2008, 40, 1478–1483. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, Y.; Medvid, R.; Melton, C.; Jaenisch, R.; Blelloch, R. DGCR8 is essential for microRNA biogenesis and silencing of embryonic stem cell self-renewal. Nat. Genet. 2007, 39, 380–385. [Google Scholar] [CrossRef] [PubMed]
  14. Qi, J.; Yu, J.Y.; Shcherbata, H.R.; Mathieu, J.; Wang, A.J.; Seal, S.; Zhou, W.; Stadler, B.M.; Bourgin, D.; Wang, L.; et al. microRNAs regulate human embryonic stem cell division. Cell Cycle 2009, 8, 3729–3741. [Google Scholar] [CrossRef]
  15. Subramanyam, D.; Lamouille, S.; Judson, R.L.; Liu, J.Y.; Bucay, N.; Derynck, R.; Blelloch, R. Multiple targets of miR-302 and miR-372 promote reprogramming of human fibroblasts to induced pluripotent stem cells. Nat. Biotechnol. 2011, 29, 443–448. [Google Scholar] [CrossRef]
  16. Kuppusamy, K.T.; Jones, D.C.; Sperber, H.; Madan, A.; Fischer, K.A.; Rodriguez, M.L.; Pabon, L.; Zhu, W.Z.; Tulloch, N.L.; Yang, X.; et al. Let-7 family of microRNA is required for maturation and adult-like metabolism in stem cell-derived cardiomyocytes. Proc. Natl. Acad. Sci. USA 2015, 112, E2785–E2794. [Google Scholar] [CrossRef]
  17. Lim, G.B. MicroRNA-directed cardiac repair after myocardial infarction in pigs. Nat. Rev. Cardiol. 2019, 16, 454–455. [Google Scholar] [CrossRef]
  18. Gao, F.; Kataoka, M.; Liu, N.; Liang, T.; Huang, Z.P.; Gu, F.; Ding, J.; Liu, J.; Zhang, F.; Ma, Q.; et al. Therapeutic role of miR-19a/19b in cardiac regeneration and protection from myocardial infarction. Nat. Commun. 2019, 10, 1802. [Google Scholar] [CrossRef]
  19. Gabisonia, K.; Prosdocimo, G.; Aquaro, G.D.; Carlucci, L.; Zentilin, L.; Secco, I.; Ali, H.; Braga, L.; Gorgodze, N.; Bernini, F.; et al. MicroRNA therapy stimulates uncontrolled cardiac repair after myocardial infarction in pigs. Nature 2019, 569, 418–422. [Google Scholar] [CrossRef]
  20. Miklas, J.W.; Clark, E.; Levy, S.; Detraux, D.; Leonard, A.; Beussman, K.; Showalter, M.R.; Smith, A.T.; Hofsteen, P.; Yang, X.; et al. TFPa/HADHA is required for fatty acid beta-oxidation and cardiolipin re-modeling in human cardiomyocytes. Nat. Commun. 2019, 10, 4671. [Google Scholar] [CrossRef]
  21. Stadler, B.M.; Ruohola-Baker, H. Small RNAs: Keeping stem cells in line. Cell 2008, 132, 563–566. [Google Scholar] [CrossRef] [PubMed]
  22. Stadler, B.; Ivanovska, I.; Mehta, K.; Song, S.; Nelson, A.; Tan, Y.; Mathieu, J.; Darby, C.; Blau, C.A.; Ware, C.; et al. Characterization of microRNAs involved in embryonic stem cell states. Stem Cells Dev. 2010, 19, 935–950. [Google Scholar] [CrossRef] [PubMed]
  23. Houbaviy, H.B.; Murray, M.F.; Sharp, P.A. Embryonic stem cell-specific MicroRNAs. Dev. Cell 2003, 5, 351–358. [Google Scholar] [CrossRef]
  24. Suh, M.R.; Lee, Y.; Kim, J.Y.; Kim, S.K.; Moon, S.H.; Lee, J.Y.; Cha, K.Y.; Chung, H.M.; Yoon, H.S.; Moon, S.Y.; et al. Human embryonic stem cells express a unique set of microRNAs. Dev. Biol. 2004, 270, 488–498. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, Z.; Hong, Y.; Xiang, D.; Zhu, P.; Wu, E.; Li, W.; Mosenson, J.; Wu, W.S. MicroRNA-302/367 cluster governs hESC self-renewal by dually regulating cell cycle and apoptosis pathways. Stem Cell Rep. 2015, 4, 645–657. [Google Scholar] [CrossRef] [PubMed]
  26. Mathieu, J.; Ruohola-Baker, H. Metabolic remodeling during the loss and acquisition of pluripotency. Development 2017, 144, 541–551. [Google Scholar] [CrossRef] [PubMed]
  27. Sperber, H.; Mathieu, J.; Wang, Y.; Ferreccio, A.; Hesson, J.; Xu, Z.; Fischer, K.A.; Devi, A.; Detraux, D.; Gu, H.; et al. The metabolome regulates the epigenetic landscape during naive-to-primed human embryonic stem cell transition. Nat. Cell Biol. 2015, 17, 1523–1535. [Google Scholar] [CrossRef]
  28. Gu, W.; Gaeta, X.; Sahakyan, A.; Chan, A.B.; Hong, C.S.; Kim, R.; Braas, D.; Plath, K.; Lowry, W.E.; Christofk, H.R. Glycolytic Metabolism Plays a Functional Role in Regulating Human Pluripotent Stem Cell State. Cell Stem Cell 2016, 19, 476–490. [Google Scholar] [CrossRef]
  29. Cao, Y.; Guo, W.T.; Tian, S.; He, X.; Wang, X.W.; Liu, X.; Gu, K.L.; Ma, X.; Huang, D.; Hu, L.; et al. miR-290/371-Mbd2-Myc circuit regulates glycolytic metabolism to promote pluripotency. EMBO J. 2015, 34, 609–623. [Google Scholar] [CrossRef]
  30. Manor, Y.S.; Massarwa, R.; Hanna, J.H. Establishing the human naive pluripotent state. Curr. Opin. Genet. Dev. 2015, 34, 35–45. [Google Scholar] [CrossRef]
  31. Nichols, J.; Smith, A. Naive and primed pluripotent states. Cell Stem Cell 2009, 4, 487–492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ware, C.B. Concise Review: Lessons from Naive Human Pluripotent Cells. Stem Cells 2017, 35, 35–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kinoshita, M.; Smith, A. Pluripotency Deconstructed. Dev. Growth Differ. 2018, 60, 44–52. [Google Scholar] [CrossRef] [PubMed]
  34. Morgani, S.; Nichols, J.; Hadjantonakis, A.K. The many faces of Pluripotency: In vitro adaptations of a continuum of in vivo states. BMC Dev. Biol. 2017, 17, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hussein, A.M.; Wang, Y.; Mathieu, J.; Margaretha, L.; Song, C.; Jones, D.C.; Cavanaugh, C.; Miklas, J.W.; Mahen, E.; Showalter, M.R.; et al. Metabolic control over mTOR dependent diapause-like state. Dev. Cell 2019. (accepted). [Google Scholar]
  36. Cheong, A.W.; Pang, R.T.; Liu, W.M.; Kottawatta, K.S.; Lee, K.F.; Yeung, W.S. MicroRNA Let-7a and dicer are important in the activation and implantation of delayed implanting mouse embryos. Hum. Reprod. 2014, 29, 750–762. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, W.M.; Pang, R.T.; Cheong, A.W.; Ng, E.H.; Lao, K.; Lee, K.F.; Yeung, W.S. Involvement of microRNA lethal-7a in the regulation of embryo implantation in mice. PLoS ONE 2012, 7, e37039. [Google Scholar] [CrossRef] [Green Version]
  38. Jouneau, A.; Ciaudo, C.; Sismeiro, O.; Brochard, V.; Jouneau, L.; Vandormael-Pournin, S.; Coppee, J.Y.; Zhou, Q.; Heard, E.; Antoniewski, C.; et al. Naive and primed murine pluripotent stem cells have distinct miRNA expression profiles. RNA 2012, 18, 253–264. [Google Scholar] [CrossRef] [Green Version]
  39. Moradi, S.; Sharifi-Zarchi, A.; Ahmadi, A.; Mollamohammadi, S.; Stubenvoll, A.; Gunther, S.; Salekdeh, G.H.; Asgari, S.; Braun, T.; Baharvand, H. Small RNA Sequencing Reveals Dlk1-Dio3 Locus-Embedded MicroRNAs as Major Drivers of Ground-State Pluripotency. Stem Cell Rep. 2017, 9, 2081–2096. [Google Scholar] [CrossRef] [Green Version]
  40. Gu, K.L.; Zhang, Q.; Yan, Y.; Li, T.T.; Duan, F.F.; Hao, J.; Wang, X.W.; Shi, M.; Wu, D.R.; Guo, W.T.; et al. Pluripotency-associated miR-290/302 family of microRNAs promote the dismantling of naive pluripotency. Cell Res. 2016, 26, 350–366. [Google Scholar] [CrossRef] [Green Version]
  41. Ware, C.B.; Wang, L.; Mecham, B.H.; Shen, L.; Nelson, A.M.; Bar, M.; Lamba, D.A.; Dauphin, D.S.; Buckingham, B.; Askari, B.; et al. Histone deacetylase inhibition elicits an evolutionarily conserved self-renewal program in embryonic stem cells. Cell Stem Cell 2009, 4, 359–369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Ware, C.B.; Nelson, A.M.; Mecham, B.; Hesson, J.; Zhou, W.; Jonlin, E.C.; Jimenez-Caliani, A.J.; Deng, X.; Cavanaugh, C.; Cook, S.; et al. Derivation of naive human embryonic stem cells. Proc. Natl. Acad. Sci. USA 2014, 111, 4484–4489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Moody, J.D.; Levy, S.; Mathieu, J.; Xing, Y.; Kim, W.; Dong, C.; Tempel, W.; Robitaille, A.M.; Dang, L.T.; Ferreccio, A.; et al. First critical repressive H3K27me3 marks in embryonic stem cells identified using designed protein inhibitor. Proc. Natl. Acad. Sci. USA 2017, 114, 10125–10130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ferreccio, A.; Mathieu, J.; Detraux, D.; Somasundaram, L.; Cavanaugh, C.; Sopher, B.; Fischer, K.; Bello, T.; Hussein, A.M.; Levy, S.; et al. Inducible CRISPR genome editing platform in naive human embryonic stem cells reveals JARID2 function in self-renewal. Cell Cycle 2018, 17, 535–549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Theunissen, T.W.; Powell, B.E.; Wang, H.; Mitalipova, M.; Faddah, D.A.; Reddy, J.; Fan, Z.P.; Maetzel, D.; Ganz, K.; Shi, L.; et al. Systematic Identification of Culture Conditions for Induction and Maintenance of Naive Human Pluripotency. Cell Stem Cell 2014, 15, 524–526. [Google Scholar] [CrossRef] [Green Version]
  46. Chou, C.H.; Shrestha, S.; Yang, C.D.; Chang, N.W.; Lin, Y.L.; Liao, K.W.; Huang, W.C.; Sun, T.H.; Tu, S.J.; Lee, W.H.; et al. miRTarBase update 2018: A resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018, 46, D296–D302. [Google Scholar] [CrossRef]
  47. Grow, E.J.; Flynn, R.A.; Chavez, S.L.; Bayless, N.L.; Wossidlo, M.; Wesche, D.J.; Martin, L.; Ware, C.B.; Blish, C.A.; Chang, H.Y.; et al. Intrinsic retroviral reactivation in human preimplantation embryos and pluripotent cells. Nature 2015, 522, 221–225. [Google Scholar] [CrossRef] [Green Version]
  48. Factor, D.C.; Corradin, O.; Zentner, G.E.; Saiakhova, A.; Song, L.; Chenoweth, J.G.; McKay, R.D.; Crawford, G.E.; Scacheri, P.C.; Tesar, P.J. Epigenomic comparison reveals activation of “seed” enhancers during transition from naive to primed pluripotency. Cell Stem Cell 2014, 14, 854–863. [Google Scholar] [CrossRef] [Green Version]
  49. Varjosalo, M.; Taipale, J. Hedgehog: Functions and mechanisms. Genes Dev. 2008, 22, 2454–2472. [Google Scholar] [CrossRef] [Green Version]
  50. Wu, S.M.; Choo, A.B.; Yap, M.G.; Chan, K.K. Role of Sonic hedgehog signaling and the expression of its components in human embryonic stem cells. Stem Cell Res. 2010, 4, 38–49. [Google Scholar] [CrossRef] [Green Version]
  51. Li, Q.; Lex, R.K.; Chung, H.; Giovanetti, S.M.; Ji, Z.; Ji, H.; Person, M.D.; Kim, J.; Vokes, S.A. The Pluripotency Factor NANOG Binds to GLI Proteins and Represses Hedgehog-mediated Transcription. J. Biol. Chem. 2016, 291, 7171–7182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Kiprilov, E.N.; Awan, A.; Desprat, R.; Velho, M.; Clement, C.A.; Byskov, A.G.; Andersen, C.Y.; Satir, P.; Bouhassira, E.E.; Christensen, S.T.; et al. Human embryonic stem cells in culture possess primary cilia with hedgehog signaling machinery. J. Cell Biol. 2008, 180, 897–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Bangs, F.K.; Schrode, N.; Hadjantonakis, A.K.; Anderson, K.V. Lineage specificity of primary cilia in the mouse embryo. Nat. Cell Biol. 2015, 17, 113–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Freedman, B.S.; Lam, A.Q.; Sundsbak, J.L.; Iatrino, R.; Su, X.; Koon, S.J.; Wu, M.; Daheron, L.; Harris, P.C.; Zhou, J.; et al. Reduced ciliary polycystin-2 in induced pluripotent stem cells from polycystic kidney disease patients with PKD1 mutations. J. Am. Soc. Nephrol. 2013, 24, 1571–1586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Mathieu, J.; Detraux, D.; Kuppers, D.; Wang, Y.; Cavanaugh, C.; Sidhu, S.; Levy, S.; Robitaille, A.M.; Ferreccio, A.; Bottorff, T.; et al. Folliculin regulates mTORC1/2 and WNT pathways in early human pluripotency. Nat. Commun. 2019, 10, 632. [Google Scholar] [CrossRef] [Green Version]
  56. Mukhopadhyay, S.; Wen, X.; Ratti, N.; Loktev, A.; Rangell, L.; Scales, S.J.; Jackson, P.K. The ciliary G-protein-coupled receptor Gpr161 negatively regulates the Sonic hedgehog pathway via cAMP signaling. Cell 2013, 152, 210–223. [Google Scholar] [CrossRef] [Green Version]
  57. Matteson, P.G.; Desai, J.; Korstanje, R.; Lazar, G.; Borsuk, T.E.; Rollins, J.; Kadambi, S.; Joseph, J.; Rahman, T.; Wink, J.; et al. The orphan G protein-coupled receptor, Gpr161, encodes the vacuolated lens locus and controls neurulation and lens development. Proc. Natl. Acad. Sci. USA 2008, 105, 2088–2093. [Google Scholar] [CrossRef] [Green Version]
  58. Bachmann, V.A.; Mayrhofer, J.E.; Ilouz, R.; Tschaikner, P.; Raffeiner, P.; Rock, R.; Courcelles, M.; Apelt, F.; Lu, T.W.; Baillie, G.S.; et al. Gpr161 anchoring of PKA consolidates GPCR and cAMP signaling. Proc. Natl. Acad. Sci. USA 2016, 113, 7786–7791. [Google Scholar] [CrossRef] [Green Version]
  59. Takashima, Y.; Guo, G.; Loos, R.; Nichols, J.; Ficz, G.; Krueger, F.; Oxley, D.; Santos, F.; Clarke, J.; Mansfield, W.; et al. Resetting transcription factor control circuitry toward ground-state pluripotency in human. Cell 2014, 158, 1254–1269. [Google Scholar] [CrossRef] [Green Version]
  60. Indrieri, A.; van Rahden, V.A.; Tiranti, V.; Morleo, M.; Iaconis, D.; Tammaro, R.; D’Amato, I.; Conte, I.; Maystadt, I.; Demuth, S.; et al. Mutations in COX7B cause microphthalmia with linear skin lesions, an unconventional mitochondrial disease. Am. J. Hum. Genet. 2012, 91, 942–949. [Google Scholar] [CrossRef] [Green Version]
  61. Izreig, S.; Samborska, B.; Johnson, R.M.; Sergushichev, A.; Ma, E.H.; Lussier, C.; Loginicheva, E.; Donayo, A.O.; Poffenberger, M.C.; Sagan, S.M.; et al. The miR-17 approximately 92 microRNA Cluster Is a Global Regulator of Tumor Metabolism. Cell Rep. 2016, 16, 1915–1928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zhou, W.; Choi, M.; Margineantu, D.; Margaretha, L.; Hesson, J.; Cavanaugh, C.; Blau, C.A.; Horwitz, M.S.; Hockenbery, D.; Ware, C.; et al. HIF1alpha induced switch from bivalent to exclusively glycolytic metabolism during ESC-to-EpiSC/hESC transition. EMBO J. 2012, 31, 2103–2116. [Google Scholar] [CrossRef] [PubMed]
  63. Zhang, J.; Ratanasirintrawoot, S.; Chandrasekaran, S.; Wu, Z.; Ficarro, S.B.; Yu, C.; Ross, C.A.; Cacchiarelli, D.; Xia, Q.; Seligson, M.; et al. LIN28 Regulates Stem Cell Metabolism and Conversion to Primed Pluripotency. Cell Stem Cell 2016, 19, 66–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Sun, Z.; Zhu, M.; Lv, P.; Cheng, L.; Wang, Q.; Tian, P.; Yan, Z.; Wen, B. The Long Noncoding RNA Lncenc1 Maintains Naive States of Mouse ESCs by Promoting the Glycolysis Pathway. Stem Cell Rep. 2018, 11, 741–755. [Google Scholar] [CrossRef] [Green Version]
  65. Bahat, A.; Goldman, A.; Zaltsman, Y.; Khan, D.H.; Halperin, C.; Amzallag, E.; Krupalnik, V.; Mullokandov, M.; Silberman, A.; Erez, A.; et al. MTCH2-mediated mitochondrial fusion drives exit from naive pluripotency in embryonic stem cells. Nat. Commun. 2018, 9, 5132. [Google Scholar] [CrossRef]
  66. Cornacchia, D.; Zhang, C.; Zimmer, B.; Chung, S.Y.; Fan, Y.; Soliman, M.A.; Tchieu, J.; Chambers, S.M.; Shah, H.; Paull, D.; et al. Lipid Deprivation Induces a Stable, Naive-to-Primed Intermediate State of Pluripotency in Human PSCs. Cell Stem Cell 2019, 25, 120–136 e10. [Google Scholar] [CrossRef]
  67. Zhang, D.; Zhao, T.; Ang, H.S.; Chong, P.; Saiki, R.; Igarashi, K.; Yang, H.; Vardy, L.A. AMD1 is essential for ESC self-renewal and is translationally down-regulated on differentiation to neural precursor cells. Genes Dev. 2012, 26, 461–473. [Google Scholar] [CrossRef] [Green Version]
  68. Zhao, T.; Goh, K.J.; Ng, H.H.; Vardy, L.A. A role for polyamine regulators in ESC self-renewal. Cell Cycle 2012, 11, 4517–4523. [Google Scholar] [CrossRef] [Green Version]
  69. James, C.; Zhao, T.Y.; Rahim, A.; Saxena, P.; Muthalif, N.A.; Uemura, T.; Tsuneyoshi, N.; Ong, S.; Igarashi, K.; Lim, C.Y.; et al. MINDY1 Is a Downstream Target of the Polyamines and Promotes Embryonic Stem Cell Self-Renewal. Stem Cells 2018, 36, 1170–1178. [Google Scholar] [CrossRef] [Green Version]
  70. Zhu, K.; Ding, H.; Wang, W.; Liao, Z.; Fu, Z.; Hong, Y.; Zhou, Y.; Zhang, C.Y.; Chen, X. Tumor-suppressive miR-218-5p inhibits cancer cell proliferation and migration via EGFR in non-small cell lung cancer. Oncotarget 2016, 7, 28075–28085. [Google Scholar] [CrossRef]
  71. Bjornsen, L.P.; Hadera, M.G.; Zhou, Y.; Danbolt, N.C.; Sonnewald, U. The GLT-1 (EAAT2; slc1a2) glutamate transporter is essential for glutamate homeostasis in the neocortex of the mouse. J. Neurochem. 2014, 128, 641–649. [Google Scholar] [CrossRef] [PubMed]
  72. TeSlaa, T.; Chaikovsky, A.C.; Lipchina, I.; Escobar, S.L.; Hochedlinger, K.; Huang, J.; Graeber, T.G.; Braas, D.; Teitell, M.A. alpha-Ketoglutarate Accelerates the Initial Differentiation of Primed Human Pluripotent Stem Cells. Cell Metab. 2016, 24, 485–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Contreras, L.; Ramirez, L.; Du, J.; Hurley, J.B.; Satrustegui, J.; de la Villa, P. Deficient glucose and glutamine metabolism in Aralar/AGC1/Slc25a12 knockout mice contributes to altered visual function. Mol. Vis. 2016, 22, 1198–1212. [Google Scholar] [PubMed]
  74. Jalil, M.A.; Begum, L.; Contreras, L.; Pardo, B.; Iijima, M.; Li, M.X.; Ramos, M.; Marmol, P.; Horiuchi, M.; Shimotsu, K.; et al. Reduced N-acetylaspartate levels in mice lacking aralar, a brain- and muscle-type mitochondrial aspartate-glutamate carrier. J. Biol. Chem. 2005, 280, 31333–31339. [Google Scholar] [CrossRef] [Green Version]
  75. Hore, T.A.; von Meyenn, F.; Ravichandran, M.; Bachman, M.; Ficz, G.; Oxley, D.; Santos, F.; Balasubramanian, S.; Jurkowski, T.P.; Reik, W. Retinol and ascorbate drive erasure of epigenetic memory and enhance reprogramming to naive pluripotency by complementary mechanisms. Proc. Natl. Acad. Sci. USA 2016, 113, 12202–12207. [Google Scholar] [CrossRef] [Green Version]
  76. Gao, Y.; Han, Z.; Li, Q.; Wu, Y.; Shi, X.; Ai, Z.; Du, J.; Li, W.; Guo, Z.; Zhang, Y. Vitamin C induces a pluripotent state in mouse embryonic stem cells by modulating microRNA expression. FEBS J. 2015, 282, 685–699. [Google Scholar] [CrossRef]
  77. Gennarino, V.A.; D’Angelo, G.; Dharmalingam, G.; Fernandez, S.; Russolillo, G.; Sanges, R.; Mutarelli, M.; Belcastro, V.; Ballabio, A.; Verde, P.; et al. Identification of microRNA-regulated gene networks by expression analysis of target genes. Genome. Res. 2012, 22, 1163–1172. [Google Scholar] [CrossRef] [Green Version]
  78. Chowdhury, F.; Li, Y.; Poh, Y.C.; Yokohama-Tamaki, T.; Wang, N.; Tanaka, T.S. Soft substrates promote homogeneous self-renewal of embryonic stem cells via downregulating cell-matrix tractions. PLoS ONE 2010, 5, e15655. [Google Scholar] [CrossRef]
  79. Gusdon, A.M.; Zhu, J.; Van Houten, B.; Chu, C.T. ATP13A2 regulates mitochondrial bioenergetics through macroautophagy. Neurobiol. Dis. 2012, 45, 962–972. [Google Scholar] [CrossRef] [Green Version]
  80. Forristal, C.E.; Wright, K.L.; Hanley, N.A.; Oreffo, R.O.; Houghton, F.D. Hypoxia inducible factors regulate pluripotency and proliferation in human embryonic stem cells cultured at reduced oxygen tensions. Reproduction 2010, 139, 85–97. [Google Scholar] [CrossRef] [Green Version]
  81. Tillander, V.; Arvidsson Nordstrom, E.; Reilly, J.; Strozyk, M.; Van Veldhoven, P.P.; Hunt, M.C.; Alexson, S.E. Acyl-CoA thioesterase 9 (ACOT9) in mouse may provide a novel link between fatty acid and amino acid metabolism in mitochondria. Cell Mol. Life Sci. 2014, 71, 933–948. [Google Scholar] [CrossRef] [PubMed]
  82. Lim, C.P.; Jain, N.; Cao, X. Stress-induced immediate-early gene, egr-1, involves activation of p38/JNK1. Oncogene 1998, 16, 2915–2926. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Kitamuro, T.; Takahashi, K.; Ogawa, K.; Udono-Fujimori, R.; Takeda, K.; Furuyama, K.; Nakayama, M.; Sun, J.; Fujita, H.; Hida, W.; et al. Bach1 functions as a hypoxia-inducible repressor for the heme oxygenase-1 gene in human cells. J. Biol. Chem. 2003, 278, 9125–9133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Molloy, M.E.; Lewinska, M.; Williamson, A.K.; Nguyen, T.T.; Kuser-Abali, G.; Gong, L.; Yan, J.; Little, J.B.; Pandolfi, P.P.; Yuan, Z.M. ZBTB7A governs estrogen receptor alpha expression in breast cancer. J. Mol. Cell Biol. 2018, 10, 273–284. [Google Scholar] [CrossRef]
  85. Alexa, A.; Rahnenfuhrer, J.; Lengauer, T. Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics 2006, 22, 1600–1607. [Google Scholar] [CrossRef] [Green Version]
  86. Nakamura, T.; Okamoto, I.; Sasaki, K.; Yabuta, Y.; Iwatani, C.; Tsuchiya, H.; Seita, Y.; Nakamura, S.; Yamamoto, T.; Saitou, M. A developmental coordinate of pluripotency among mice, monkeys and humans. Nature 2016, 537, 57–62. [Google Scholar] [CrossRef]
  87. Miyaoka, Y.; Chan, A.H.; Judge, L.M.; Yoo, J.; Huang, M.; Nguyen, T.D.; Lizarraga, P.P.; So, P.L.; Conklin, B.R. Isolation of single-base genome-edited human iPS cells without antibiotic selection. Nat. Methods 2014, 11, 291–293. [Google Scholar] [CrossRef] [Green Version]
  88. He, B.; Zhang, H.; Wang, J.; Liu, M.; Sun, Y.; Guo, C.; Lu, J.; Wang, H.; Kong, S. Blastocyst activation engenders transcriptome reprogram affecting X-chromosome reactivation and inflammatory trigger of implantation. Proc. Natl. Acad. Sci. USA 2019, 116, 16621–16630. [Google Scholar] [CrossRef] [Green Version]
  89. Chen, M.; Reed, R.R.; Lane, A.P. Chronic Inflammation Directs an Olfactory Stem Cell Functional Switch from Neuroregeneration to Immune Defense. Cell Stem Cell 2019, 25, 501–513 e5. [Google Scholar] [CrossRef]
Ijms 20 05864 g001 550
Figure 1. microRNAs regulating human naïve to primed ESCs transition: (A) A schematic figure of early embryonic developmental stages. (B) Analysis workflow. We first identified 357 differentially expressed microRNAs and 1146 differentially expressed protein-coding genes in two naïve-primed studies [27,45]. We then used mirTarBase to connect changes in microRNA and their experimentally validated target genes, and filtered down to 2184 miR-target gene connections where microRNA is up and its target is down (or vice versa). Green √ means the microRNA-gene connection is considered consistent; red × means the connection is not consistent. (C) Gene ontology enrichment of microRNA target genes with lower expression in human naïve ESCs (the microRNA regulators are higher in naïve). x-axis is negative log10 of enrichment p-value (larger means more significant).
Figure 1. microRNAs regulating human naïve to primed ESCs transition: (A) A schematic figure of early embryonic developmental stages. (B) Analysis workflow. We first identified 357 differentially expressed microRNAs and 1146 differentially expressed protein-coding genes in two naïve-primed studies [27,45]. We then used mirTarBase to connect changes in microRNA and their experimentally validated target genes, and filtered down to 2184 miR-target gene connections where microRNA is up and its target is down (or vice versa). Green √ means the microRNA-gene connection is considered consistent; red × means the connection is not consistent. (C) Gene ontology enrichment of microRNA target genes with lower expression in human naïve ESCs (the microRNA regulators are higher in naïve). x-axis is negative log10 of enrichment p-value (larger means more significant).
Ijms 20 05864 g001
Ijms 20 05864 g002 550
Figure 2. The Hh pathway in naïve-to-primed transition: (A) Schematic representation of the role of shh pathway during early development. Shh pathway is not active in ESC but plays an important role in various developmental processes. (B) Expression of shh pathway components in naïve hESC (Elf1 2iLIF, Sperber et al. and Grow et al.) and primed hESC (Elf1 AF, Grow et al. and H1 Sperber et al.). (C,D) Cilia are present only in primed hESC, not in naïve hESC, as visualized by Arl13B staining in naïve (Elf1 2iLIF) and primed (Elf1 TeSR) hESC. Scale bars represent 5 μm, triangles indicate the presence of primary cilia. Standard Error of Mean (s.e.m); *** p < 0.001; 2-tailed t-test. (E) Sanger sequencing results of the GPR161 gene around the gRNA in Elf1 wild type (WT) and GPR161 mutant show an insertion of the nucleotide T, resulting in a frame shift and premature STOP codon (left panel). This mutation leads to the generation of a truncated GPR161 protein missing the functional C terminus domain responsible for regulation of PKA and shh activity (right panel). (F) Cilia (ARL13B) immunofluorescence staining in Elf1 WT and Elf1 GPR161 mutant grown in primed conditions (TeSR). On the right side are high magnification images of the boxed regions shown on the left side. Scale bars represent 5 μm. (G) Shh component SMO is up-regulated in Elf1 GPR161 mutant cells compared to WT in primed conditions (qPCR analysis). S.e.m.; * p < 0.05; 2-tailed t-test. (H) Model of repression of shh activity in naïve and primed hESC through microRNAs. T-bars indicate repression, arrow upwards upregulation and downwards downregulation.
Figure 2. The Hh pathway in naïve-to-primed transition: (A) Schematic representation of the role of shh pathway during early development. Shh pathway is not active in ESC but plays an important role in various developmental processes. (B) Expression of shh pathway components in naïve hESC (Elf1 2iLIF, Sperber et al. and Grow et al.) and primed hESC (Elf1 AF, Grow et al. and H1 Sperber et al.). (C,D) Cilia are present only in primed hESC, not in naïve hESC, as visualized by Arl13B staining in naïve (Elf1 2iLIF) and primed (Elf1 TeSR) hESC. Scale bars represent 5 μm, triangles indicate the presence of primary cilia. Standard Error of Mean (s.e.m); *** p < 0.001; 2-tailed t-test. (E) Sanger sequencing results of the GPR161 gene around the gRNA in Elf1 wild type (WT) and GPR161 mutant show an insertion of the nucleotide T, resulting in a frame shift and premature STOP codon (left panel). This mutation leads to the generation of a truncated GPR161 protein missing the functional C terminus domain responsible for regulation of PKA and shh activity (right panel). (F) Cilia (ARL13B) immunofluorescence staining in Elf1 WT and Elf1 GPR161 mutant grown in primed conditions (TeSR). On the right side are high magnification images of the boxed regions shown on the left side. Scale bars represent 5 μm. (G) Shh component SMO is up-regulated in Elf1 GPR161 mutant cells compared to WT in primed conditions (qPCR analysis). S.e.m.; * p < 0.05; 2-tailed t-test. (H) Model of repression of shh activity in naïve and primed hESC through microRNAs. T-bars indicate repression, arrow upwards upregulation and downwards downregulation.
Ijms 20 05864 g002
Ijms 20 05864 g003 550
Figure 3. The polyamine pathway is dysregulated in naïve to primed transition: (A) Ornithine is converted to polyamines (spermine and spermidine) by the rate limiting enzyme ornithine decarboxylase (ODC1). ODC1 is consistently up-regulated in the naïve state, and its microRNA regulator, miR-218-bp, is lower in the naïve state. Each rectangle box is a metabolite; each arrow is a reaction. Red means higher abundance in naïve state; blue means lower abundance in naïve state. “--|” (T-bar) indicates repression of ODC1 by miR-218-5p. (B) Log2 fold change of ODC1 in three separate naïve vs. primed cell line comparisons. (C) ODC1 expression changes in monkey pre-to-post-implantation transition in vivo. ODC1 is also expressed higher in the naïve state in vivo. Blue line represents a local regression fit of the data (D) Bar plot of ornithine, citrulline, and spermidine abundance in naïve (Elf1, H12i) and primed (H1) cell lines.
Figure 3. The polyamine pathway is dysregulated in naïve to primed transition: (A) Ornithine is converted to polyamines (spermine and spermidine) by the rate limiting enzyme ornithine decarboxylase (ODC1). ODC1 is consistently up-regulated in the naïve state, and its microRNA regulator, miR-218-bp, is lower in the naïve state. Each rectangle box is a metabolite; each arrow is a reaction. Red means higher abundance in naïve state; blue means lower abundance in naïve state. “--|” (T-bar) indicates repression of ODC1 by miR-218-5p. (B) Log2 fold change of ODC1 in three separate naïve vs. primed cell line comparisons. (C) ODC1 expression changes in monkey pre-to-post-implantation transition in vivo. ODC1 is also expressed higher in the naïve state in vivo. Blue line represents a local regression fit of the data (D) Bar plot of ornithine, citrulline, and spermidine abundance in naïve (Elf1, H12i) and primed (H1) cell lines.
Ijms 20 05864 g003
Ijms 20 05864 g004 550
Figure 4. microRNAs regulating mouse naïve to primed ESCs transition: (A) Venn diagram of microRNAs changing in the same direction in naïve to primed transition across three microRNA-seq studies of mouse ESC vs. mouse EpiSC. (B) Gene ontology enrichment of microRNA target genes with lower expression in mouse ESCs (the microRNA regulators are higher in naïve). The x-axis is the negative log10 of enrichment p-value (larger means more significant).
Figure 4. microRNAs regulating mouse naïve to primed ESCs transition: (A) Venn diagram of microRNAs changing in the same direction in naïve to primed transition across three microRNA-seq studies of mouse ESC vs. mouse EpiSC. (B) Gene ontology enrichment of microRNA target genes with lower expression in mouse ESCs (the microRNA regulators are higher in naïve). The x-axis is the negative log10 of enrichment p-value (larger means more significant).
Ijms 20 05864 g004
Ijms 20 05864 g005 550
Figure 5. Consistent microRNA changes across human and mouse pluripotency: (A) Left panel: microRNAs with lower expression in the diapause state compared to reactivated state. MicroRNAs that are expressed lower only in diapause but not in the naïve state are highlighted in blue. Right panel: Gene ontology enrichment of microRNA target genes with higher expression in mouse diapause embryos compared to post-implantation embryos (the microRNA regulators are lower in diapause, shown on the left panel). The x-axis is the negative log10 of enrichment p-value (larger means more significant). (B) Left panel: Top 10 microRNAs with highest expression in the diapause state compared to reactivated state. These microRNAs are highlighted red because they are higher in diapause but not in naïve state. Right panel: Gene ontology enrichment of microRNA target genes with lower expression in mouse diapause embryos compared to post-implantation embryos (the microRNA regulators are higher in diapause, shown on the left panel).
Figure 5. Consistent microRNA changes across human and mouse pluripotency: (A) Left panel: microRNAs with lower expression in the diapause state compared to reactivated state. MicroRNAs that are expressed lower only in diapause but not in the naïve state are highlighted in blue. Right panel: Gene ontology enrichment of microRNA target genes with higher expression in mouse diapause embryos compared to post-implantation embryos (the microRNA regulators are lower in diapause, shown on the left panel). The x-axis is the negative log10 of enrichment p-value (larger means more significant). (B) Left panel: Top 10 microRNAs with highest expression in the diapause state compared to reactivated state. These microRNAs are highlighted red because they are higher in diapause but not in naïve state. Right panel: Gene ontology enrichment of microRNA target genes with lower expression in mouse diapause embryos compared to post-implantation embryos (the microRNA regulators are higher in diapause, shown on the left panel).
Ijms 20 05864 g005
Table
Table 1. Datasets used in this study.
Table 1. Datasets used in this study.
StudyTypeConditionGrowth MediaAccession Number
Sperber H. et al. [27]microRNA-seqNaïve hESC2iL-I-FGSE60995
Sperber H. et al. [27]mRNA-seqNaïve and primed hESC2iL-I-F/KSR + FGFGSE60995
Jouneau A. et al. [38]microRNA-seqMouse ESC and EpiSC15% FBS-LIF/AFSupplemental data from original paper
Gu et al. [40]microRNA-seqMouse ESC differentiation time courseN2B27-2iLIF/KSR + FGFSupplemental data from original paper
Moradi S. et al. [39]microRNA-seqMouse ground ESC and naïve ESCN2B27 2iLIF or 3iLIF/15% FBS-LIFGSE87174
ENCODE projectmicroRNA-seqprimed ESCTeSRhttps://www.encodeproject.org
Roadmap Epigenome projectmicroRNA-seqprimed ESCTeSRhttp://www.roadmapepigenomics.org
Theunissen et al. [45]mRNA-seqNaïve/primed hESC5iLA/serumGSE59435
Grow et al. [47]mRNA-seqNaïve/primed hESC2iLIF/KSR + FGFGSE63570
Factor et al. [48]mRNA-seqMouse ESC and EpiSCKSR LIF/KSR + FGFGSE57409
Hussein et al. [35]mRNA-seqMouse pre-implantation, diapause, post-implantationFrom embryos
Liu et al. [37]microRNA microarrayMouse diapause and re-activated embryosFrom embryosSupplemental data from original paper
Table
Table 2. Metabolic differences between naïve and primed pluripotent stem cells.
Table 2. Metabolic differences between naïve and primed pluripotent stem cells.
NaivePrimedAssayReferences
Glycolysis ↑Glycolysis ↑RNAseq, Metabolomics, Seahorse Flux Analyzer, TMREZhou et al. [62]
Takashima et al. [59]
Sperber et al. [27]
Gu et al. [28]
Zhang et al. [63]
Sun et al. [64]
Baha et al. [65]
Cornacchia et al. [66]
OXPHOS ↑OXPHOS ↓
Fatty acid synthesis ↓Fatty acid synthesis ↑RNAseq, Metabolomics, Seahorse Flux Analyzer, Oil RedO and Bodipy stainingSperber et al. [27]
Cornacchia et al. [66]
Fatty acid oxidation ↑Fatty acid oxidation ↓
↑ indicates up-regulation, ↓ indicates down-regulation.
Table
Table 3. Hypothetic role of microRNAs in regulation of fatty acid metabolism of naïve and primed embryonic stem cells (ESC).
Table 3. Hypothetic role of microRNAs in regulation of fatty acid metabolism of naïve and primed embryonic stem cells (ESC).
NaivePrimedAssayReferences
FA synthesis, activation, elongation, desaturation
(FASN, ACSL1, ELOVL7, SCD)
Ijms 20 05864 i001
miR-10a		FA transporter (rate limiting FAO)
CPT1A
Ijms 20 05864 i001
miR-106b-5p, miR-20a-5p, miR-302c-3p, miR-17-5p, miR-20b-5p, miR-106a-5p, miR-302a-3p, miR-150-5p, miR-16-5p, miR-93-5p, miR-4517		RNAseq, microRNAseq, RT-qPCRSperber et al. [27]
T bars indicate microRNA based repression of the target genes.
Table
Table 4. microRNAs showing consistent changes in at least three out of four datasets.
Table 4. microRNAs showing consistent changes in at least three out of four datasets.
IDlog2FC.Ref_Epigenomelog2FC.ENCODEpadj.Ref_Epigenomepadj.ENCODEJouneauGuMoradiHow Many Consistent
hsa-mir-143, hsa-miR-143-3p4.107.502.932 × 10−31.163 × 10−72.161.700.794
hsa-mir-200c, hsa-miR-200c-3p−6.60−8.978.956 × 10−965.294 × 10−936.44−1.03−0.883
hsa-mir-205, hsa-miR-205-5p3.239.284.715 × 10−121.392 × 10−882.230.93−1.283
hsa-mir-302c, hsa-miR-302c-3p−1.68−6.405.552 × 10−131.086 × 10−66−5.47−1.122.633
hsa-mir-20b, hsa-miR-20b-5p−1.93−5.371.923 × 10−162.778 × 10−510.00−1.60−0.753
hsa-mir-335,hsa-miR-335-5p−1.76−5.071.075 × 10−136.196 × 10−470.00−1.00−0.743
hsa-mir-302a, hsa-miR-302a-3p−1.50−4.461.090 × 10−108.146 × 10−39−5.15−2.081.673
hsa-mir-412, hsa-miR-412-5p7.1715.023.547 × 10−641.038 × 10−28−2.600.923.343
hsa-mir-433, hsa-miR-433-3p5.3815.026.536 × 10−397.545 × 10−210.000.624.293
hsa-mir-495, hsa-miR-495-3p4.687.472.492 × 10−341.098 × 10−190.000.604.273
hsa-mir-582, hsa-miR-582-3p2.725.917.063 × 10−124.664 × 10−141.400.960.483
hsa-mir-196a-1, hsa-miR-196a-5p8.9615.021.038 × 10−151.389 × 10−130.002.460.783
hsa-mir-196a-2, hsa-miR-196a-5p8.6915.022.845 × 10−122.660 × 10−120.002.460.783
hsa-mir-200a, hsa-miR-200a-5p−3.37−4.711.511 × 10−65.775 × 10−90.00−0.89−1.593
hsa-mir-1247, hsa-miR-1247-3p3.4015.021.677 × 10−84.527 × 10−62.850.002.173
hsa-mir-708, hsa-miR-708-3p1.131.455.997 × 10−53.027 × 10−41.870.720.083
hsa-mir-199b, hsa-miR-199b-3p3.344.274.641 × 10−26.808 × 10−32.540.301.013
hsa-let-7f-2, hsa-let-7f-5p4.794.431.363 × 10−28.710 × 10−30.001.970.783
hsa-let-7f-1, hsa-let-7f-5p4.724.311.455 × 10−29.959 × 10−30.001.970.783
hsa-mir-182, hsa-miR-182-3p3.194.107.430 × 10−42.571 × 10−20.000.640.883

Share and Cite

MDPI and ACS Style

Wang, Y.; Hussein, A.M.; Somasundaram, L.; Sankar, R.; Detraux, D.; Mathieu, J.; Ruohola-Baker, H. microRNAs Regulating Human and Mouse Naïve Pluripotency. Int. J. Mol. Sci. 2019, 20, 5864. https://doi.org/10.3390/ijms20235864

AMA Style

Wang Y, Hussein AM, Somasundaram L, Sankar R, Detraux D, Mathieu J, Ruohola-Baker H. microRNAs Regulating Human and Mouse Naïve Pluripotency. International Journal of Molecular Sciences. 2019; 20(23):5864. https://doi.org/10.3390/ijms20235864

Chicago/Turabian Style

Wang, Yuliang, Abdiasis M. Hussein, Logeshwaran Somasundaram, Rithika Sankar, Damien Detraux, Julie Mathieu, and Hannele Ruohola-Baker. 2019. "microRNAs Regulating Human and Mouse Naïve Pluripotency" International Journal of Molecular Sciences 20, no. 23: 5864. https://doi.org/10.3390/ijms20235864

APA Style

Wang, Y., Hussein, A. M., Somasundaram, L., Sankar, R., Detraux, D., Mathieu, J., & Ruohola-Baker, H. (2019). microRNAs Regulating Human and Mouse Naïve Pluripotency. International Journal of Molecular Sciences, 20(23), 5864. https://doi.org/10.3390/ijms20235864

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop