BACKGROUND
micrornas , short  conserved endogenous non-coding rnas, inhibit messenger rna targets by repressing translation or reducing mrna stability  <cit> .
mirnas critically modulate many cellular events, including the balance between proliferation and differentiation during organ development  <cit> .
in pluripotent stem cells , mirnas play important roles in regulating stem cell bioprocesses .
mirnas modulate stem cell pluripotency and differentiation .
knocking out the key mirna processing enzymes dicer  <cit>  or dgcr <dig>  <cit>  causes stem cells to lose their pluripotency.
mirna- <dig> cluster has been proposed to regulate the core pluripotency factors like pou5f <dig> .
mirna-302- <dig> cluster has also been used to induce pluripotency  <cit> .
on the other hand, mirnas like let- <dig> induce stem cell differentiation  <cit> .
however, these recent studies have mostly focused on individual gene functions in stem cells although genome-wide data might be employed, and the conclusions drawn from these current studies are unavoidably biased on genes selected by these studies.
therefore, these studies only provide partial mechanisms of mirna functions in stem cells, and the overall systems mechanisms of how mirnas regulate stem cell processes remain largely elusive.
mirnas generally do not work alone to perform their functions  <cit> .
one mirna might target more than  <dig> genes  <cit> , and one gene can be repressed by multiple mirnas in a sequence-specific fashion  <cit> .
in turn, proteins can physically bind to the promoters and enhancers of mirnas to regulate mirna activations  <cit> .
these binary interactions between mirnas and proteins would form a complex systematic network.
this complexity of mirna interaction network presents a challenge for conventional approaches like gene-knockout to unbiasedly capture the real mechanisms of mirna functions in stem cells.
this present study employed systems physical network approaches  <cit>  and constructed a comprehensive and unbiased map of genome-wide interactions between mirnas and their targets to investigate the global basis of mirna roles in pluripotent stem cells, from the pluripotent self-renewal state to the differentiation state.
results of the present study lay a conceptual framework for future studies and applications of mirnas in stem cells.
RESULTS
physical network of mirna and protein interactions in stem cells
to systematically reveal the roles of mirnas in stem cells, this study first constructed a systems network  <cit>  of interactions between mirnas and proteins by combining multiple published datasets.
these interactions contain binary interactions from two directions, from mirnas to mrnas coding for proteins, and from proteins to mirna promoters and enhancers.
the mirna-targets were inferred from clip-seq data, which simultaneously identify mirnas-mrna interactions by measuring mirna-argonaute complexes  <cit>  .
protein-targets were inferred from chip-seq, which measures protein interactions with dna  <cit>  .
the clip-seq and chip-seq provide data of physical binding interactions.
the systems interaction network constructed here includes physical interactions of both mirna-targets and protein-targets.fig.
<dig> overall view.
a, the concept of mirna and protein interactions.
b, workflow of this study.
c, overall view of the entire network constructed by this study.
d, a sample of entire network contents shows direct interactions between mirnas and pluripotent core factors .
nodes denote mirnas or proteins as labeled; red node represents the gene up-regulation in pluripotent stem cells, green node as down-regulation, and white node as insignificant expression.
edges  represent interactions; red edge denotes mirnas targeting proteins  and green edge represents the binding of proteins with the regulatory elements of mirnas .
this labeling strategy applies to all fig.s in this study
the entire network contains ~ <dig>  nodes and ~ <dig>  interactions .
both chip-seq and clip-seq measures genome-wide targets and thus this network provides a global map of mirna targets in stem cells.
for example, the genome-wide interactions between mirnas and the pluripotent core factors  could be extracted from this network .
due to its natural interactions generated from experimental data, this network would provide accurate interactions between mirna and their targets, and results generated from this network would be reliable.
the primary role of mirnas in stem cells
to understand the primary role of all mirnas activated  at the pluripotent state in stem cells, we searched for the biological functions of the network activated by mirnas in stem cells.
to avoid the dataset biases and to ensure our results general, we included different datasets of mirnas and genes coding for proteins and selected mirnas and genes that are only activated with high frequency in all datasets .
these activated mirnas and genes were used to enrich the entire network  to get the activated network using methods as previously described  <cit> .
the network activated by overexpressed mirnas and down-regulated proteins  formed a subnetwork activated by up-regulated mirnas in the pluripotent state.
to avoid bias, we run the go  functional analysis  of this subnetwork  <cit>  by separately using references of both entire go annotation database and a set of all up-regulated genes in stem cells.
different references generated different enrichment p-value but produced the similar result that the entire subnetwork primarily functions for development, with corrected p-value <  <dig> e- <dig>  .
this activated subnetwork was further enriched by the developmental go term to obtain the developmental module .
this developmental module was decomposed into functional modules based on network topology  <cit>  and it contained  <dig> sub-functional modules .
all of these  <dig> modules primarily function for development , indicating that the primary function of over-expressed mirnas at pluripotent state is to repress developmental modules.fig.
<dig> mirnas primarily repress developmental processes in pluripotent stem cells.
a, up-regulated mirnas primarily target a developmental module, which includes  <dig> primary sub-modules functioning for development as highlighted in cycles.
b, the key modules targeted by the top  <dig> important mirnas .
c, the key modules were significantly and abundantly enriched in the developmental category
the repressing role of up-regulated mirnas was further evidenced by examining the targets of three primarily represented mirna groups, the top overexpressed mirna group, a mir-302- <dig> cluster, and a single mir- <dig>  first, a total of  <dig> out of the  <dig> most important mirnas, which were selected on the basis of the variance contribution to the system , directly target a developmental sub-network .
among the  <dig> mirnas, the top  <dig> mirnas also target a module that primarily functions for development .
furthermore, the well-known mir- <dig> cluster and even a single mirna, mir- <dig>  also target modules functionally enriched in the developmental category with respectively corr p-value < <dig> e- <dig> , and corr p-value <  <dig> e- <dig> .
therefore, overexpressed mirnas at pluripotent state primarily repress development.
biologically, to maintain the self-renewal and pluripotency, stem cells have some ways to prevent development and differentiation.
this repressing function of mirnas can help mirnas claim their contributions to the stem cell properties at the pluripotent state.fig.
<dig> modules targeted by mir-302- <dig> cluster and mirna- <dig>  a, modules targeted by the mir-302- <dig> cluster and their functional enrichment was shown in b. c, a module targeted by mirna- <dig>  and d, its functional enrichment
on the other hand, the down-regulated mirnas in stem cells directly target genes that primarily function for metabolism and pluripotency, such as sox <dig> and pou5f <dig> .
these down-regulated mirnas become up-regulated during differentiation and directly slower the entire metabolic processes and inhibit pluripotency during this state.
together, mirnas primarily and directly repress development during the pluripotent state while they repress metabolism and pluripotency during cell differentiation.
paths from activated mirnas to pluripotent core factors
mirnas like mir- <dig> cluster mediate pluripotency  <cit> , and it thus was assumed that top over-expressed mirnas such as mir- <dig> and mirna-290- <dig> cluster in stem cells might directly or indirectly target pluripotent core factors  <cit> .
to investigate if these top mirnas directly target the core factors, we systematically searched the shortest paths respectively from the top important mirnas as described above  to the three core factors .
surprisingly, all these top mirnas, including mir-302b, mir- <dig>  mir- <dig>  and mir- <dig>  do not directly target any core factors .
actually, the direct basis of these mirnas regulating the pluripotent core factors seemed blurred because all these mirnas must go through at least  <dig> steps to reach any pluripotency core factor.
these two steps include a mirna and a protein, which are not consistently expressed with these mirnas.
this suggested that these top over-expressed mirnas in stem cells do not directly mediate the pluripotency.fig.
<dig> pathways from activated mirnas to pluripotent core factors.
a, b, c, d, the shortest paths from mir-302b, mir- <dig>  mir- <dig>  and mir-292-5p respectively to the pluripotent core factors.
none of these top overexpressed mirnas directly binds to any of these core factors.
e, core factors directly targeted by activated mirnas in stem cells
to circumvent the bias of selected mirna groups, we further globally and unbiasedly searched for all up-regulated mirnas that target the pluripotent core factors .
surprisingly, we only found one mirna  that barely up-regulated   <cit>   in the pluripotent state and directly binds to sox <dig> , which was also targeted by down-regulated mirna- <dig>  this indicated that activated mirnas  do not primarily and directly target the pluripotent core factors during the pluripotent state.
on the other hand, regardless of expression, only limited mirnas directly bind to pou5f <dig>  while many mirnas directly target sox <dig> although they might not over-express in stem cells , indicating that mirnas target sox <dig> more than pou5f <dig>  this suggested that the connection from mirnas to the pluripotent core factors primarily go through sox <dig>  in contrast with the current thought that mirnas should primarily target pou5f <dig>  <cit> , a key factor for reprogramming induced pluripotent stem cells.
pluripotent network targeted by activated mirnas
we next expanded the pluripotent gene list to all pluripotent genes uncovered by hu et al.
<cit> .
we still focused on the direct mirna target and searched the first neighbor of all overexpressed mirnas  and all down-regulated mirnas .
most of these targets are shared by up- and down-regulated mirnas , indicating that the primary pluripotent genes in stem cells are carefully modulated by multiple up- and down-regulated mirnas.
the function of the entire shared network is primarily for extrinsic signal pathways associated with pluripotency .
for example, the highly connected nodes apc, rad <dig> and eif4g <dig> are involved in wnt signaling and mitotic cell cycle pathways  <cit>  .
signaling pathways with similar functions were also found in the network targeted by over-expressed mirnas only  and in modules directly targeted by the represented mirnas cluster in stem cells, such as mir-302- <dig> cluster  and mir- <dig> .
this indicated that mirnas in the pluripotent state primarily function for modulating the balance of pluripotent signal pathways instead of directly targeting pluripotent core factors.
this suggests that these regulations driven by mirnas might go through multiple steps to the pluripotent core factors.fig.
<dig> pluripotent genes targeted by over-expressed mirnas in stem cells.
a. venn diagram of pluripotent genes targeted by both up-and down-regulated mirnas.
b, activated modules commonly targeted by both up- and down-regulated mirnas.
the enlarged nodes represent the highly connected targets that are targeted by both up- and down-regulated mirnas in stem cells, and they work for pluripotent signal pathways.
c, down-regulated modules targeted by over-expressed mirnas in stem cells.
d, activated modules targeted by the mir-302- <dig> cluster in stem cells.
e, an activated network targeted by mir-294
interestingly, il6st  were regulated by multiple mirnas in pluripotent state .
a protein complex constructed by il6st  and lifr  and stabilized by lif  is a regulatory complex to maintain self-renewal in stem cells  <cit> .
targeting il6st here further suggested the important role of mirnas in modulating pluripotency.
in contrast to the up-regulated mirnas in pluripotent state, down-regulated mirnas directly target the core pluripotency factors like mirna <dig> targeting sox <dig> .
these down-regulated mirnas in the pluripotent state are actually up-regulated during the differentiation state and these pluripotency core factors become down-regulated in differentiation.
this direct linkage between up-regulated mirnas and down-regulated core factors suggested that these mirnas directly inhibit the pluripotent core factors and repress pluripotency to facilitate differentiation during stem cell differentiation.
this result parallels with recent experimental observations showing that mirnas repress the pluripotent core factors and pluripotency during differentiation  <cit> .
to summarize, mirnas do not directly target pluripotent core factors during pluripotent state but mirnas directly target and repress these core factors during differentiation.
mirnas abundantly target epigenetic system
the above result indicated that the number of mirna binding  to pluripotent genes  was very limited.
the target with the highest degree, eif4g <dig>  was only attacked by ~ <dig> up-regulated mirnas .
it was expected that a certain group of nodes should be targeted by more than that.
this drove us to further search the network hubs  in the entire network.
we systematically ranked the mirna targets by degree  and obtained the top hubs.
the top hubs mostly function for rna processing, but surprisingly, dnmt3a, a dna-methyltransferase for de novo dna methylation, was among the top hubs.
dnmt3a actually holds more than  <dig> mirna binding sites in 3'-utr region based only on 8 bp seed mapping and it was ranked within the top 1 % of the up-regulated mirna targets .
dnmt3a was even ranked higher than eif4g <dig>  the highest ranked node in the pluripotent genes , indicating that dnmt3a should be a top important node in the network directly regulated by mirnas in the pluripotent state.
this also indicated that mirnas predominately target dna methylation system, rather than the pluripotent genes.
we extracted the network of dnmt3a directly targeted by mirnas and found that the top mirna clusters , such as mir-302- <dig> and mir-290- <dig> cluster, were among the mirnas that target dnmt3a .
similarly, many well-known mirna clusters  in stem cells also target dnmt <dig> , an enzyme predominately responsible for methylation in hemimethylated cpg islands.
many down-regulated mirnas also target dnmt3a and dnmt <dig> , but their attacks would lead to differentiation instead of maintaining pluripotency in stem cells as discussed above  <cit> .
these abundantly overexpressed mirnas that target the methylation system suggest that mirnas predominately repress dnmts in the pluripotent state.fig.
<dig> mirnas abundantly target dna methylation systems.
a, dnmt3a was targeted by  <dig> up-regulated mirnas and was ranked in the top 1 % of the up-regulated mirna targets.
the network nodes  were plotted against the node degree .
b, the dnmt3a network targeted by mir nas.
the most important mirnas, including the mir-302- <dig> and mir-290- <dig> cluster, are found in the dnmt3a network.
c, mirnas target another methylation enzyme dnmt1
in addition, mirnas directly and abundantly target a core histone modification complex  , including mef2c , which was targeted by the top over-expressed mirna clusters including mir-290- <dig> and mir- <dig> cluster .
up-regulating mef2c enhances stem cells differentiation  <cit> , and down-regulated mef2c inhibits differentiation  <cit> .
the down-regulation of mef2c targeted by the top over-expressed mirna clusters suggests that mirnas repress differentiation in stem cells.
this is consistent with our discussion above on the mirna repressing development and differentiation at the pluripotent state in stem cells .
together, mirnas directly and abundantly target the epigenetic systems at the pluripotent state.fig.
<dig> mirnas directly target a histone complex.
the mir-290- <dig> cluster heavily attacks the mef2c
dna methylation mediates the mirna activation in stem cells
to search the mechanism controlling the mirna activations, this study turned to the genome-wide sequencing of dna methylation in stem cells and methylation-loss-stem cells  <cit>  .
a total of 2000 bp in each upstream and downstream of start sites of all activated mirnas were examined.
while the dna methylation in the downstream of up-regulated mirnas is not different from that of down-regulated mirnas , the down-regulated mirnas hold significantly higher methylation upstream than up-regulated mirnas .
surprisingly, the biggest difference locates in ~1000 bp up-stream instead of immediate up-stream .
furthermore, these differences are overall negatively correlated to mirna expressions with correlation coefficient of − <dig>  and p-value <  <dig>  .
this suggested that the difference in dna methylation accounts for the mirna activations.
this parallels a recent observation showing that the loss of dna methylation significantly increases mirna expressions  <cit> .
therefore, mirna activations and their network are mediated by dna methylation in ~1000 bp upstream regions.fig.
<dig> dna methylation mediates mirna activations in stem cells.
a, the methylation levels upstream and downstream from the mirna start site.
down-regulated mirnas have significantly higher methylation in upstream region when compared with up-regulated mirnas.
b, detailed methylation profiling for regions that are 2000 bp upstream from the activated mirnas.
the top  <dig> down-regulated mirnas  have higher methylation around -1000 bp  than that of the top  <dig> up-regulated mirnas .
the methylation profiling of a single down-regulated mirna- <dig> as a representative example .
c. negative correlation between dna methylation and mirna expression.
red line represents regression line
discussion
this study is the first investigation, to our knowledge, to search the primary mechanism of mirna functions in stem cells at systems level on the basis of a physical map constructed by direct interactions of mirnas and proteins.
mirnas play critical roles in maintaining overall properties like pluripotency in stem cells , but the mechanisms still remain largely elusive.
here, we systematically revealed that mirna primarily target developmental modules during the pluripotent state while another set of mirnas directly target pluripotent core factors during differentiation state.
this suggests that one of mirna primary functions is to repress development at the pluripotent state to prevent stem cell differentiation and to keep stem cell pluripotency while mirnas also degrade pluripotent core factors to facilitate differentiation during differentiation state.
this is consistent with the recent observation that overexpressions of mirnas induce pluripotency  <cit>  and mirnas also facilitate stem cells differentiation  <cit> .
the mechanism of mirnas in regulating pluripotency has been widely investigated , but the systems mechanism remains elusive.
one of fundamental questions is whether mirnas directly or indirectly regulate the pluripotency at system level.
it has been consistently observed that gene expressions of the top over-expressed mirnas are positive correlative to that of pluripotent core factors.
recent evidences also show that mirna-302- <dig> cluster could induce pluripotency  <cit> .
the observations led to a speculation that mirnas might directly target the pluripotent core factors.
although mirnas could have many functions in certain conditions, the primary functions of mirnas are for degrading and inhibiting their targets.
if the overexpressed mirnas directly target the core factors, these mirnas would likely repress the core factors as previously evidenced  <cit> , leading to down-regulations of these core factors.
a negative correlation between them was expected, but in fact, a positive correlation has been consistently observed.
this suggests that the top mirnas might not directly target and degrade pluripotent core factors.
here, we utilized the power of our system network to exhaustively search the direct linkages between mirnas and the pluripotent core factors.
our results revealed that the top mirnas such as mir- <dig> and mirna- <dig> cluster do not directly target any core pluripotent factors during the pluripotent state .
most of top mirnas only directly target the extrinsic signal pathways associated with pluripotency .
a signal from a signal pathway normally travels through multiple steps to reach its targets.
therefore, the top mirnas indirectly target the core factors and indirectly regulate the pluripotency.
in contrast to the indirect mechanism of mirnas in targeting the genetic system, the top mirnas, including mir-302- <dig> and mirna- <dig> cluster, directly and abundantly target the core enzymes of dna methylation system, dnmt3a and dnmt <dig>  targeting these dnmts by the top expressed mirnas suggests inhibition of dnmts, which highly express in somatic cells but lowly express in stem cells.
this parallels the most recent observations showing that gene expressions of up-regulated mirnas are negatively correlative with that of dna methyltransferases   <cit>  and that micrornas degrade dnmts in stem cells  <cit> .
during self-renewal and pluripotency, mouse stem cell genome are in hypomethylated status, and genome-wide three methylase  are tightly depressed  <cit> .
our data showed that abundant mirna directly target dna methylase, suggesting that mirnas play a crucial role in repressing dna methylation machinery during pluripotency.
our result of mirna directly repressing dnmts contributes to the understanding of the positive relationship existing between up-regulated mirnas and overexpressed core pluripotent factors in stem cells as discussed above.
over-expressed mirnas directly target dnmts, leading to dna methylation reductions at genome-wide level, including the loci of the pluripotent core factors.
this results in over-expression of these core factors responsible for pluripotency.
this is consistent with the observation of less methylation in the loci of the core factors during induced pluripotent stem cell reprogramming, and it is also consistent with our recent finding that the demethylation level can be used as a variable for discriminating different stem cells  <cit> .
therefore, mirnas primarily and directly target the epigenetic system that further activates pluripotent core factors in stem cells.
this parallels the most current report that mirnas degrade dnmts in stem cells  <cit> .
dna methylation might regulate expressions of a certain group of genes in stem cells  <cit> .
our data further showed that dna methylation globally mediates the mirna activations in stem cells.
these mirna activations by dna methylation in turn repress the dna methylation.
less methylation activates the mirnas and pluripotent core factors again.
together, we proposed a system-wide circuit to describe a part of mirna primary roles in modulating pluripotency in pluripotent stem cells .
in this circuit, mirnas directly repress development and directly repress the dna methylation system, while mirnas indirectly regulate pluripotency genes.
this repression of dna methylation activates both pluripotent factors and mirnas.
the activations of pluripotent factors and repression of development contribute to pluripotency in stem cells, while the activations of mirnas further inhibit both dna methylation and development.
this creates an active system-wide circuit in stem cells to maintain the pluripotent state .fig.
<dig> a system view of mirnas primary mechanistic roles in maintaining pluripotency in pluripotent stem cells.
mirnas directly repress both the dna methylation system and the development, while indirectly regulate pluripotency genes.
reduced dna methylation activates the mirnas and pluripotent core factors for pluripotency.
the activated circuit between mirnas and dna methylation, as well as the development inhibition, help stem cells to maintain the pluripotent state, see text for detail.
the solid lines are directly derived from the present study and they represent direct interactions and the dash dark line represents indirect interactions, while dash green lines denote evidences from reference papers.
the network interactions of proteins targeted by multiple mirnas also provide a systems vision of quantitative regulatory mechanism of mirnas in stem cells.
for example, our study uncovered that many up- and down-regulated mirnas directly target dnmts and wnt pathway like apc, suggesting that dna methylation and wnt pathway are delicately regulated in stem cells dynamically.
these regulations are controlled by not only simply turning on/off genes but also mediated by the quantitative concentrations of many gene expressions.
the slightly quantitative alternation of each component concentration could result in change of stem cell state, self-renewal or differentiation.
this quantitative regulatory mechanism can help to revise the controversial mechanisms derived from individual gene studies.
for example, wnt signaling and pathway components interact with pluripotent core factors  to regulate stem cell pluripotency and differentiation  <cit> , but whether wnt signaling promotes pluripotency or differentiation has been intensively debated recently  <cit> , partially due to lack of system-wide vision on it.
these components targeted by mirnas revealed here with other system-network level data could provide novel clues to design biological experiments and select mathematical models to clarify the role of wnt pathway in stem cells and to finally solve this controversial issue.
CONCLUSIONS
this study mostly focused on the primary functions of mirnas in pluripotent state.
we systematically revealed that mirnas directly repress dna methylation systems and primarily inhibiting development to help stem cells to maintain pluripotency in the pluripotent state.
however, we do not exclude other functions of mirnas in stem cells.
the results made here are based on the current limited data and these mirna functions drawn here only account for a part of mirna roles in stem cells.
with data accumulating, more functions of mirnas will certainly be explored.
recently, an exciting field of noncoding rna studies have revealed that long noncoding rnas , transcripts of more than  <dig> nucleotides, participate in the regulation of the pluripotent state and differentiation.
these incrnas even regulate the complex mirna network and interact with dna methylation machinery like methylase  <cit> .
a more exciting comprehensive map of mirna-lncrna-dna methylation-mrna can be expected.
nevertheless, our key results and conclusion regarding the fundamental basis of primary mirna functions in regulating mouse stem cell pluripotency and differentiation have actually also been observed in human  <cit> , suggesting that these mirna functions might be conserved in both mouse and human.
understanding the fundamental systems roles of mirnas studied here would broadly direct the future functional studies of mirnas in both human and mouse stem cells.
materials and methods
data resources
this study analyzed multiple datasets published in public database, geo database .
these data were generated by high-throughput methods as detailed in additional file 1: table s <dig>  including clip-seq, chip-seq, rna-seq, microarray, and bisulfite sequencing dna .
data bioinformatics analysis
to be consistent and comparable, all sequencing data were mapped to mouse genome .
all fastq files  from clip-seq, chip-seq and rna-seq were aligned using bwa  <dig> .
<dig> with default parameters  <cit> .
all pcr-duplicates were removed.
bisulfite sequencing dna was aligned by bismark  <dig> .
<dig> with tolerating one non-bisulfite mismatch per read .
the clusters from clip-seq were generated by using genomicranges  <dig>  under r  <dig>  and were subjected to second noise quartile cutoff.
the filtered clusters located in 3’utrs and cdss were used to search for mirna bindings.
the mirna bindings were searched against the perfect match of mirna seed sequences  using home-made python scripts.
the seed sequences were extracted from mirna sequences downloaded from mirbase  <dig> .
peaks from chip-seq were called using sissrs  <dig>   and the peaks were annotated using chippeakanno  <dig>  under r- <dig> .
the differential expressions from rna-seq were performed using negative binomial model using edger  <dig>  under r  <dig> .
bisulfite dna methylation was counted against the genome coordinates of mirnas .
microarray data were analyzed using limma  <dig>  under r  <dig> .
activated genes and mirnas selection
to minimize the biases from individual experiments and cell lines as well as the confounders caused by high-throughput methods, we selected genes and mirnas activated in different conditions by using different sets of data .
expression data from various resources  were employed.
genes coding for proteins with activating frequency >50 % in all observations and mirnas with activating frequency of more than 25 % out of all observations were treated as activated genes and were selected .
the up-regulation and down-regulation was based on comparison to somatic cells in each experiment .
the top important mirnas  were selected on the basis of their contributions to network structure and variance by using the algorithm as we previously published  <cit> .
briefly, the top mirnas were selected on the basis of their ranking scores calculated by the eigengene-based connectivity as defined below  <cit> .
score=didmax+2×cor|xi,e)| where di denotes the ith node degree, and dmax denotes the maximum degree of a node in the entire network.
|cor| represents the absolute pearson correlation coefficient, where xi represents a vector of ith node value, and e eigengene of the network.
network construction and analysis
the network construction and analysis were performed by approaches as our previous reports  <cit> .
briefly, we built the map  by collecting the interactions of both mirnas targets from the clip-seq and protein-binding promoters of mirnas from chip-seq data.
the interactions from clip-seq and chip-seq were signed as mirna or proteins respectively .
only the direct interactions  were selected and included, and thus this map is a physical binding network.
the network was enriched by the activated genes and mirnas selected above to get the activated network.
functional modules were further enriched by the functional genes based on gene ontology enrichment .
six sub-modules in functional developmental module  were based on network topology to identify the densely connected modules.
the target node ranking was based on the degree of each mirna node .
additional file
additional file 1:
figure s <dig>  functional enrichments of all targets that are targeted by all up- and down-regulated mirnas.
a, up-regulated mirnas.
b, down-regulated mirnas.
figure s <dig>  total  <dig> out of top  <dig> important mirnas target developmental genes.
figure s <dig>  down-regulated mirnas mediate metabolism in stem cells.
figure s <dig>  all pluripotent genes targeted by up-regulated mirnas  and down-regulated mirnas  in stem cells.
figure s <dig>  an activated network targeted by down-regulated mirnas.
figure s <dig>  down-regulated mirnas targeting dnmt3a.
table s <dig>  data sources.
table s <dig>  a list of mirnas that are frequently and significantly differential expressed between stem cells and somatic cells.
table s <dig>  a list of genes coding for proteins that are frequently and significantly differential expressed between stem cells and somatic cells.
table s <dig>  a list of top  <dig> important mirnas.
competing interests
the authors declare that they have no competing interest.
author contributions
a.w, conceived and designed the experiments, analyzed the data, contributed data/materials/analysis tools, wrote the manuscript.
y.z.
contributed data/materials/analysis tools.
q.h.
contributed data/materials/analysis tools/writing.
all authors reviewed the manuscript.
we specially thank drs.
ying du and chunxiao zhou for data collection.

 micrornas, short conserved endogenous non-coding rnas, inhibit messenger rna targets by repressing translation or reducing mrna stability. mirnas critically modulate many cellular events, including the balance between proliferation and differentiation during organ development. in pluripotent stem cells, mirnas play important roles in regulating stem cell bioprocesses. mirnas modulate stem cell pluripotency and differentiation. knocking out the key mirna processing enzymes dicer or dgcr causes stem cells to lose their pluripotency. mirna- cluster has been proposed to regulate the core pluripotency factors like pou5f. mirna-302- cluster has also been used to induce pluripotency. on the other hand, mirnas like let- induce stem cell differentiation. however, these recent studies have mostly focused on individual gene functions in stem cells although genome-wide data might be employed, and the conclusions drawn from these current studies are unavoidably biased on genes selected by these studies. therefore, these studies only provide partial mechanisms of mirna functions in stem cells, and the overall systems mechanisms of how mirnas regulate stem cell processes remain largely elusive. mirnas generally do not work alone to perform their functions. one mirna might target more than genes, and one gene can be repressed by multiple mirnas in a sequence-specific fashion. in turn, proteins can physically bind to the promoters and enhancers of mirnas to regulate mirna activations. these binary interactions between mirnas and proteins would form a complex systematic network. this complexity of mirna interaction network presents a challenge for conventional approaches like gene-knockout to unbiasedly capture the real mechanisms of mirna functions in stem cells. this present study employed systems physical network approaches and constructed a comprehensive and unbiased map of genome-wide interactions between mirnas and their targets to investigate the global basis of mirna roles in pluripotent stem cells, from the pluripotent self-renewal state to the differentiation state. results of the present study lay a conceptual framework for future studies and applications of mirnas in stem cells.  physical network of mirna and protein interactions in stem cells to systematically reveal the roles of mirnas in stem cells, this study first constructed a systems network of interactions between mirnas and proteins by combining multiple published datasets. these interactions contain binary interactions from two directions, from mirnas to mrnas coding for proteins, and from proteins to mirna promoters and enhancers. the mirna-targets were inferred from clip-seq data, which simultaneously identify mirnas-mrna interactions by measuring mirna-argonaute complexes. protein-targets were inferred from chip-seq, which measures protein interactions with dna. the clip-seq and chip-seq provide data of physical binding interactions. the systems interaction network constructed here includes physical interactions of both mirna-targets and protein-targets.fig. overall view. a, the concept of mirna and protein interactions. b, workflow of this study. c, overall view of the entire network constructed by this study. d, a sample of entire network contents shows direct interactions between mirnas and pluripotent core factors. nodes denote mirnas or proteins as labeled; red node represents the gene up-regulation in pluripotent stem cells, green node as down-regulation, and white node as insignificant expression. edges represent interactions; red edge denotes mirnas targeting proteins and green edge represents the binding of proteins with the regulatory elements of mirnas. this labeling strategy applies to all fig.s in this study the entire network contains ~ nodes and ~ interactions. both chip-seq and clip-seq measures genome-wide targets and thus this network provides a global map of mirna targets in stem cells. for example, the genome-wide interactions between mirnas and the pluripotent core factors could be extracted from this network. due to its natural interactions generated from experimental data, this network would provide accurate interactions between mirna and their targets, and results generated from this network would be reliable. the primary role of mirnas in stem cells to understand the primary role of all mirnas activated at the pluripotent state in stem cells, we searched for the biological functions of the network activated by mirnas in stem cells. to avoid the dataset biases and to ensure our results general, we included different datasets of mirnas and genes coding for proteins and selected mirnas and genes that are only activated with high frequency in all datasets. these activated mirnas and genes were used to enrich the entire network to get the activated network using methods as previously described. the network activated by overexpressed mirnas and down-regulated proteins formed a subnetwork activated by up-regulated mirnas in the pluripotent state. to avoid bias, we run the go functional analysis of this subnetwork by separately using references of both entire go annotation database and a set of all up-regulated genes in stem cells. different references generated different enrichment p-value but produced the similar result that the entire subnetwork primarily functions for development, with corrected p-value <  e-. this activated subnetwork was further enriched by the developmental go term to obtain the developmental module. this developmental module was decomposed into functional modules based on network topology and it contained sub-functional modules. all of these modules primarily function for development, indicating that the primary function of over-expressed mirnas at pluripotent state is to repress developmental modules.fig. mirnas primarily repress developmental processes in pluripotent stem cells. a, up-regulated mirnas primarily target a developmental module, which includes primary sub-modules functioning for development as highlighted in cycles. b, the key modules targeted by the top important mirnas. c, the key modules were significantly and abundantly enriched in the developmental category the repressing role of up-regulated mirnas was further evidenced by examining the targets of three primarily represented mirna groups, the top overexpressed mirna group, a mir-302- cluster, and a single mir- first, a total of out of the most important mirnas, which were selected on the basis of the variance contribution to the system, directly target a developmental sub-network. among the mirnas, the top mirnas also target a module that primarily functions for development. furthermore, the well-known mir- cluster and even a single mirna, mir- also target modules functionally enriched in the developmental category with respectively corr p-value < e-, and corr p-value <  e-. therefore, overexpressed mirnas at pluripotent state primarily repress development. biologically, to maintain the self-renewal and pluripotency, stem cells have some ways to prevent development and differentiation. this repressing function of mirnas can help mirnas claim their contributions to the stem cell properties at the pluripotent state.fig. modules targeted by mir-302- cluster and mirna- a, modules targeted by the mir-302- cluster and their functional enrichment was shown in b. c, a module targeted by mirna- and d, its functional enrichment on the other hand, the down-regulated mirnas in stem cells directly target genes that primarily function for metabolism and pluripotency, such as sox and pou5f. these down-regulated mirnas become up-regulated during differentiation and directly slower the entire metabolic processes and inhibit pluripotency during this state. together, mirnas primarily and directly repress development during the pluripotent state while they repress metabolism and pluripotency during cell differentiation. paths from activated mirnas to pluripotent core factors mirnas like mir- cluster mediate pluripotency, and it thus was assumed that top over-expressed mirnas such as mir- and mirna-290- cluster in stem cells might directly or indirectly target pluripotent core factors. to investigate if these top mirnas directly target the core factors, we systematically searched the shortest paths respectively from the top important mirnas as described above to the three core factors. surprisingly, all these top mirnas, including mir-302b, mir- mir- and mir- do not directly target any core factors. actually, the direct basis of these mirnas regulating the pluripotent core factors seemed blurred because all these mirnas must go through at least steps to reach any pluripotency core factor. these two steps include a mirna and a protein, which are not consistently expressed with these mirnas. this suggested that these top over-expressed mirnas in stem cells do not directly mediate the pluripotency.fig. pathways from activated mirnas to pluripotent core factors. a, b, c, d, the shortest paths from mir-302b, mir- mir- and mir-292-5p respectively to the pluripotent core factors. none of these top overexpressed mirnas directly binds to any of these core factors. e, core factors directly targeted by activated mirnas in stem cells to circumvent the bias of selected mirna groups, we further globally and unbiasedly searched for all up-regulated mirnas that target the pluripotent core factors. surprisingly, we only found one mirna that barely up-regulated in the pluripotent state and directly binds to sox, which was also targeted by down-regulated mirna- this indicated that activated mirnas do not primarily and directly target the pluripotent core factors during the pluripotent state. on the other hand, regardless of expression, only limited mirnas directly bind to pou5f while many mirnas directly target sox although they might not over-express in stem cells, indicating that mirnas target sox more than pou5f this suggested that the connection from mirnas to the pluripotent core factors primarily go through sox in contrast with the current thought that mirnas should primarily target pou5f, a key factor for reprogramming induced pluripotent stem cells. pluripotent network targeted by activated mirnas we next expanded the pluripotent gene list to all pluripotent genes uncovered by hu et al.. we still focused on the direct mirna target and searched the first neighbor of all overexpressed mirnas and all down-regulated mirnas. most of these targets are shared by up- and down-regulated mirnas, indicating that the primary pluripotent genes in stem cells are carefully modulated by multiple up- and down-regulated mirnas. the function of the entire shared network is primarily for extrinsic signal pathways associated with pluripotency. for example, the highly connected nodes apc, rad and eif4g are involved in wnt signaling and mitotic cell cycle pathways. signaling pathways with similar functions were also found in the network targeted by over-expressed mirnas only and in modules directly targeted by the represented mirnas cluster in stem cells, such as mir-302- cluster and mir-. this indicated that mirnas in the pluripotent state primarily function for modulating the balance of pluripotent signal pathways instead of directly targeting pluripotent core factors. this suggests that these regulations driven by mirnas might go through multiple steps to the pluripotent core factors.fig. pluripotent genes targeted by over-expressed mirnas in stem cells. a. venn diagram of pluripotent genes targeted by both up-and down-regulated mirnas. b, activated modules commonly targeted by both up- and down-regulated mirnas. the enlarged nodes represent the highly connected targets that are targeted by both up- and down-regulated mirnas in stem cells, and they work for pluripotent signal pathways. c, down-regulated modules targeted by over-expressed mirnas in stem cells. d, activated modules targeted by the mir-302- cluster in stem cells. e, an activated network targeted by mir-294 interestingly, il6st were regulated by multiple mirnas in pluripotent state. a protein complex constructed by il6st and lifr and stabilized by lif is a regulatory complex to maintain self-renewal in stem cells. targeting il6st here further suggested the important role of mirnas in modulating pluripotency. in contrast to the up-regulated mirnas in pluripotent state, down-regulated mirnas directly target the core pluripotency factors like mirna targeting sox. these down-regulated mirnas in the pluripotent state are actually up-regulated during the differentiation state and these pluripotency core factors become down-regulated in differentiation. this direct linkage between up-regulated mirnas and down-regulated core factors suggested that these mirnas directly inhibit the pluripotent core factors and repress pluripotency to facilitate differentiation during stem cell differentiation. this result parallels with recent experimental observations showing that mirnas repress the pluripotent core factors and pluripotency during differentiation. to summarize, mirnas do not directly target pluripotent core factors during pluripotent state but mirnas directly target and repress these core factors during differentiation. mirnas abundantly target epigenetic system the above result indicated that the number of mirna binding to pluripotent genes was very limited. the target with the highest degree, eif4g was only attacked by ~ up-regulated mirnas. it was expected that a certain group of nodes should be targeted by more than that. this drove us to further search the network hubs in the entire network. we systematically ranked the mirna targets by degree and obtained the top hubs. the top hubs mostly function for rna processing, but surprisingly, dnmt3a, a dna-methyltransferase for de novo dna methylation, was among the top hubs. dnmt3a actually holds more than mirna binding sites in 3'-utr region based only on 8 bp seed mapping and it was ranked within the top 1 % of the up-regulated mirna targets. dnmt3a was even ranked higher than eif4g the highest ranked node in the pluripotent genes, indicating that dnmt3a should be a top important node in the network directly regulated by mirnas in the pluripotent state. this also indicated that mirnas predominately target dna methylation system, rather than the pluripotent genes. we extracted the network of dnmt3a directly targeted by mirnas and found that the top mirna clusters, such as mir-302- and mir-290- cluster, were among the mirnas that target dnmt3a. similarly, many well-known mirna clusters in stem cells also target dnmt, an enzyme predominately responsible for methylation in hemimethylated cpg islands. many down-regulated mirnas also target dnmt3a and dnmt, but their attacks would lead to differentiation instead of maintaining pluripotency in stem cells as discussed above. these abundantly overexpressed mirnas that target the methylation system suggest that mirnas predominately repress dnmts in the pluripotent state.fig. mirnas abundantly target dna methylation systems. a, dnmt3a was targeted by up-regulated mirnas and was ranked in the top 1 % of the up-regulated mirna targets. the network nodes were plotted against the node degree. b, the dnmt3a network targeted by mir nas. the most important mirnas, including the mir-302- and mir-290- cluster, are found in the dnmt3a network. c, mirnas target another methylation enzyme dnmt1 in addition, mirnas directly and abundantly target a core histone modification complex, including mef2c, which was targeted by the top over-expressed mirna clusters including mir-290- and mir- cluster. up-regulating mef2c enhances stem cells differentiation, and down-regulated mef2c inhibits differentiation. the down-regulation of mef2c targeted by the top over-expressed mirna clusters suggests that mirnas repress differentiation in stem cells. this is consistent with our discussion above on the mirna repressing development and differentiation at the pluripotent state in stem cells. together, mirnas directly and abundantly target the epigenetic systems at the pluripotent state.fig. mirnas directly target a histone complex. the mir-290- cluster heavily attacks the mef2c dna methylation mediates the mirna activation in stem cells to search the mechanism controlling the mirna activations, this study turned to the genome-wide sequencing of dna methylation in stem cells and methylation-loss-stem cells. a total of 2000 bp in each upstream and downstream of start sites of all activated mirnas were examined. while the dna methylation in the downstream of up-regulated mirnas is not different from that of down-regulated mirnas, the down-regulated mirnas hold significantly higher methylation upstream than up-regulated mirnas. surprisingly, the biggest difference locates in ~1000 bp up-stream instead of immediate up-stream. furthermore, these differences are overall negatively correlated to mirna expressions with correlation coefficient of − and p-value < . this suggested that the difference in dna methylation accounts for the mirna activations. this parallels a recent observation showing that the loss of dna methylation significantly increases mirna expressions. therefore, mirna activations and their network are mediated by dna methylation in ~1000 bp upstream regions.fig. dna methylation mediates mirna activations in stem cells. a, the methylation levels upstream and downstream from the mirna start site. down-regulated mirnas have significantly higher methylation in upstream region when compared with up-regulated mirnas. b, detailed methylation profiling for regions that are 2000 bp upstream from the activated mirnas. the top down-regulated mirnas have higher methylation around -1000 bp than that of the top up-regulated mirnas. the methylation profiling of a single down-regulated mirna- as a representative example. c. negative correlation between dna methylation and mirna expression. red line represents regression line discussion this study is the first investigation, to our knowledge, to search the primary mechanism of mirna functions in stem cells at systems level on the basis of a physical map constructed by direct interactions of mirnas and proteins. mirnas play critical roles in maintaining overall properties like pluripotency in stem cells, but the mechanisms still remain largely elusive. here, we systematically revealed that mirna primarily target developmental modules during the pluripotent state while another set of mirnas directly target pluripotent core factors during differentiation state. this suggests that one of mirna primary functions is to repress development at the pluripotent state to prevent stem cell differentiation and to keep stem cell pluripotency while mirnas also degrade pluripotent core factors to facilitate differentiation during differentiation state. this is consistent with the recent observation that overexpressions of mirnas induce pluripotency and mirnas also facilitate stem cells differentiation. the mechanism of mirnas in regulating pluripotency has been widely investigated, but the systems mechanism remains elusive. one of fundamental questions is whether mirnas directly or indirectly regulate the pluripotency at system level. it has been consistently observed that gene expressions of the top over-expressed mirnas are positive correlative to that of pluripotent core factors. recent evidences also show that mirna-302- cluster could induce pluripotency. the observations led to a speculation that mirnas might directly target the pluripotent core factors. although mirnas could have many functions in certain conditions, the primary functions of mirnas are for degrading and inhibiting their targets. if the overexpressed mirnas directly target the core factors, these mirnas would likely repress the core factors as previously evidenced, leading to down-regulations of these core factors. a negative correlation between them was expected, but in fact, a positive correlation has been consistently observed. this suggests that the top mirnas might not directly target and degrade pluripotent core factors. here, we utilized the power of our system network to exhaustively search the direct linkages between mirnas and the pluripotent core factors. our results revealed that the top mirnas such as mir- and mirna- cluster do not directly target any core pluripotent factors during the pluripotent state. most of top mirnas only directly target the extrinsic signal pathways associated with pluripotency. a signal from a signal pathway normally travels through multiple steps to reach its targets. therefore, the top mirnas indirectly target the core factors and indirectly regulate the pluripotency. in contrast to the indirect mechanism of mirnas in targeting the genetic system, the top mirnas, including mir-302- and mirna- cluster, directly and abundantly target the core enzymes of dna methylation system, dnmt3a and dnmt targeting these dnmts by the top expressed mirnas suggests inhibition of dnmts, which highly express in somatic cells but lowly express in stem cells. this parallels the most recent observations showing that gene expressions of up-regulated mirnas are negatively correlative with that of dna methyltransferases and that micrornas degrade dnmts in stem cells. during self-renewal and pluripotency, mouse stem cell genome are in hypomethylated status, and genome-wide three methylase are tightly depressed. our data showed that abundant mirna directly target dna methylase, suggesting that mirnas play a crucial role in repressing dna methylation machinery during pluripotency. our result of mirna directly repressing dnmts contributes to the understanding of the positive relationship existing between up-regulated mirnas and overexpressed core pluripotent factors in stem cells as discussed above. over-expressed mirnas directly target dnmts, leading to dna methylation reductions at genome-wide level, including the loci of the pluripotent core factors. this results in over-expression of these core factors responsible for pluripotency. this is consistent with the observation of less methylation in the loci of the core factors during induced pluripotent stem cell reprogramming, and it is also consistent with our recent finding that the demethylation level can be used as a variable for discriminating different stem cells. therefore, mirnas primarily and directly target the epigenetic system that further activates pluripotent core factors in stem cells. this parallels the most current report that mirnas degrade dnmts in stem cells. dna methylation might regulate expressions of a certain group of genes in stem cells. our data further showed that dna methylation globally mediates the mirna activations in stem cells. these mirna activations by dna methylation in turn repress the dna methylation. less methylation activates the mirnas and pluripotent core factors again. together, we proposed a system-wide circuit to describe a part of mirna primary roles in modulating pluripotency in pluripotent stem cells. in this circuit, mirnas directly repress development and directly repress the dna methylation system, while mirnas indirectly regulate pluripotency genes. this repression of dna methylation activates both pluripotent factors and mirnas. the activations of pluripotent factors and repression of development contribute to pluripotency in stem cells, while the activations of mirnas further inhibit both dna methylation and development. this creates an active system-wide circuit in stem cells to maintain the pluripotent state.fig. a system view of mirnas primary mechanistic roles in maintaining pluripotency in pluripotent stem cells. mirnas directly repress both the dna methylation system and the development, while indirectly regulate pluripotency genes. reduced dna methylation activates the mirnas and pluripotent core factors for pluripotency. the activated circuit between mirnas and dna methylation, as well as the development inhibition, help stem cells to maintain the pluripotent state, see text for detail. the solid lines are directly derived from the present study and they represent direct interactions and the dash dark line represents indirect interactions, while dash green lines denote evidences from reference papers. the network interactions of proteins targeted by multiple mirnas also provide a systems vision of quantitative regulatory mechanism of mirnas in stem cells. for example, our study uncovered that many up- and down-regulated mirnas directly target dnmts and wnt pathway like apc, suggesting that dna methylation and wnt pathway are delicately regulated in stem cells dynamically. these regulations are controlled by not only simply turning on/off genes but also mediated by the quantitative concentrations of many gene expressions. the slightly quantitative alternation of each component concentration could result in change of stem cell state, self-renewal or differentiation. this quantitative regulatory mechanism can help to revise the controversial mechanisms derived from individual gene studies. for example, wnt signaling and pathway components interact with pluripotent core factors to regulate stem cell pluripotency and differentiation, but whether wnt signaling promotes pluripotency or differentiation has been intensively debated recently, partially due to lack of system-wide vision on it. these components targeted by mirnas revealed here with other system-network level data could provide novel clues to design biological experiments and select mathematical models to clarify the role of wnt pathway in stem cells and to finally solve this controversial issue.  this study mostly focused on the primary functions of mirnas in pluripotent state. we systematically revealed that mirnas directly repress dna methylation systems and primarily inhibiting development to help stem cells to maintain pluripotency in the pluripotent state. however, we do not exclude other functions of mirnas in stem cells. the results made here are based on the current limited data and these mirna functions drawn here only account for a part of mirna roles in stem cells. with data accumulating, more functions of mirnas will certainly be explored. recently, an exciting field of noncoding rna studies have revealed that long noncoding rnas, transcripts of more than nucleotides, participate in the regulation of the pluripotent state and differentiation. these incrnas even regulate the complex mirna network and interact with dna methylation machinery like methylase. a more exciting comprehensive map of mirna-lncrna-dna methylation-mrna can be expected. nevertheless, our key results and conclusion regarding the fundamental basis of primary mirna functions in regulating mouse stem cell pluripotency and differentiation have actually also been observed in human, suggesting that these mirna functions might be conserved in both mouse and human. understanding the fundamental systems roles of mirnas studied here would broadly direct the future functional studies of mirnas in both human and mouse stem cells. materials and methods data resources this study analyzed multiple datasets published in public database, geo database. these data were generated by high-throughput methods as detailed in additional file 1: table s including clip-seq, chip-seq, rna-seq, microarray, and bisulfite sequencing dna. data bioinformatics analysis to be consistent and comparable, all sequencing data were mapped to mouse genome. all fastq files from clip-seq, chip-seq and rna-seq were aligned using bwa. with default parameters. all pcr-duplicates were removed. bisulfite sequencing dna was aligned by bismark. with tolerating one non-bisulfite mismatch per read. the clusters from clip-seq were generated by using genomicranges under r and were subjected to second noise quartile cutoff. the filtered clusters located in 3’utrs and cdss were used to search for mirna bindings. the mirna bindings were searched against the perfect match of mirna seed sequences using home-made python scripts. the seed sequences were extracted from mirna sequences downloaded from mirbase. peaks from chip-seq were called using sissrs and the peaks were annotated using chippeakanno under r-. the differential expressions from rna-seq were performed using negative binomial model using edger under r. bisulfite dna methylation was counted against the genome coordinates of mirnas. microarray data were analyzed using limma under r. activated genes and mirnas selection to minimize the biases from individual experiments and cell lines as well as the confounders caused by high-throughput methods, we selected genes and mirnas activated in different conditions by using different sets of data. expression data from various resources were employed. genes coding for proteins with activating frequency >50 % in all observations and mirnas with activating frequency of more than 25 % out of all observations were treated as activated genes and were selected. the up-regulation and down-regulation was based on comparison to somatic cells in each experiment. the top important mirnas were selected on the basis of their contributions to network structure and variance by using the algorithm as we previously published. briefly, the top mirnas were selected on the basis of their ranking scores calculated by the eigengene-based connectivity as defined below. score=didmax+2×cor|xi,e)| where di denotes the ith node degree, and dmax denotes the maximum degree of a node in the entire network. |cor| represents the absolute pearson correlation coefficient, where xi represents a vector of ith node value, and e eigengene of the network. network construction and analysis the network construction and analysis were performed by approaches as our previous reports. briefly, we built the map by collecting the interactions of both mirnas targets from the clip-seq and protein-binding promoters of mirnas from chip-seq data. the interactions from clip-seq and chip-seq were signed as mirna or proteins respectively. only the direct interactions were selected and included, and thus this map is a physical binding network. the network was enriched by the activated genes and mirnas selected above to get the activated network. functional modules were further enriched by the functional genes based on gene ontology enrichment. six sub-modules in functional developmental module were based on network topology to identify the densely connected modules. the target node ranking was based on the degree of each mirna node. additional file additional file 1: figure s functional enrichments of all targets that are targeted by all up- and down-regulated mirnas. a, up-regulated mirnas. b, down-regulated mirnas. figure s total out of top important mirnas target developmental genes. figure s down-regulated mirnas mediate metabolism in stem cells. figure s all pluripotent genes targeted by up-regulated mirnas and down-regulated mirnas in stem cells. figure s an activated network targeted by down-regulated mirnas. figure s down-regulated mirnas targeting dnmt3a. table s data sources. table s a list of mirnas that are frequently and significantly differential expressed between stem cells and somatic cells. table s a list of genes coding for proteins that are frequently and significantly differential expressed between stem cells and somatic cells. table s a list of top important mirnas. competing interests the authors declare that they have no competing interest. author contributions a.w, conceived and designed the experiments, analyzed the data, contributed data/materials/analysis tools, wrote the manuscript. y.z. contributed data/materials/analysis tools. q.h. contributed data/materials/analysis tools/writing. all authors reviewed the manuscript. we specially thank drs. ying du and chunxiao zhou for data collection.
furthermore, at the pluripotent state mirnas predominately repress dna methyltransferases, the core enzymes for dna methylation. globally, mirnas directly repress the pluripotent core factors during differentiation state. surprisingly, during the pluripotent state, the top important mirnas do not directly regulate the pluripotent core factors as previously thought, but they only directly target the pluripotent signal pathways and directly repress developmental processes. the decreasing methylation repressed by mirnas in turn activates the top mirnas and pluripotent core factors, creating an active circuit system to modulate pluripotency. while mirnas directly repress pluripotent core factors to facilitate differentiation during the differentiation state, they also help stem cells to maintain pluripotency by activating pluripotent cores through directly repressing dna methylation systems and primarily inhibiting development in the pluripotent state. this study systematically analyzes multiple-omics data and builds a systems physical network including genome-wide interactions between mirnas and their targets to reveal the systems mechanism of mirna functions in mouse pluripotent stem cells. mirnasglobal mechanismstem cellssystems networkpluripotencyfunctionsissue-copyright-statement© the author 2015
 micrornas critically modulate stem cell properties like pluripotency, but the fundamental mechanism remains largely unknown. method this study systematically analyzes multiple-omics data and builds a systems physical network including genome-wide interactions between mirnas and their targets to reveal the systems mechanism of mirna functions in mouse pluripotent stem cells.  globally, mirnas directly repress the pluripotent core factors during differentiation state. surprisingly, during the pluripotent state, the top important mirnas do not directly regulate the pluripotent core factors as previously thought, but they only directly target the pluripotent signal pathways and directly repress developmental processes. furthermore, at the pluripotent state mirnas predominately repress dna methyltransferases, the core enzymes for dna methylation. the decreasing methylation repressed by mirnas in turn activates the top mirnas and pluripotent core factors, creating an active circuit system to modulate pluripotency.  mirnas vary their functions with stem cell states. while mirnas directly repress pluripotent core factors to facilitate differentiation during the differentiation state, they also help stem cells to maintain pluripotency by activating pluripotent cores through directly repressing dna methylation systems and primarily inhibiting development in the pluripotent state. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. key words mirnasglobal mechanismstem cellssystems networkpluripotencyfunctionsissue-copyright-statement© the author 2015 