Volume 11 Supplement 1
Analysis of deep sequencing microRNA expression profile from human embryonic stem cells derived mesenchymal stem cells reveals possible role of let-7 microRNA family in downstream targeting of Hepatic Nuclear Factor 4 Alpha
© Koh et al; licensee BioMed Central Ltd. 2010
Published: 10 February 2010
Recent literature has revealed that genetic exchange of microRNA between cells can be essential for cell-cell communication, tissue-specificity and developmental processes. In stem cells, as in other cells, this can be accomplished through microvesicles or exosome mediated transfer. However, molecular profiles and functions of microRNAs within the cells and in their exosomes are poorly studied. Next generation sequencing technologies could provide a broad-spectrum of microRNAs and their expression and identify possible microRNA targets. In this work, we performed deep sequencing of microRNAs to understand the profile and expression of the microRNAs in microvesicles and intracellular environment of human embryonic stem cells derived mesenchymal stem cells (hES-MSC).
We outline a workflow pertaining to visualizing, statistical analysis and interpreting deep sequencing data of known intracellular and extracellular microRNAs from hES-MSC). We utilized these results of which directed our attention towards establishing hepatic nuclear factor 4 alpha (HNF4A) as a downstream target of let-7 family of microRNAs.
In our study, significant differences in expression profile of microRNAs were found in the intracellular and extracellular environment of hES-MSC. However, a high level of let-7 family of microRNAs is predominant in both intra- and extra- cellular samples of hES-MSC. Further results derived from visualization of our alignment data and network analysis showed that let-7 family microRNAs could affect the downstream target HNF4A, which is a known endodermal differentiation marker. The elevated presence of let-7 microRNA in both intracellular and extra cellular environment further suggests a possible intercellular signalling mechanism through microvesicles transfer. We suggest that let-7 family microRNAs might play a signalling role via such a mechanism amongst populations of stem cells in maintaining self renewal property by suppressing HNF4A expression. This is in line with recent paradigm where microRNAs regulate self-renewal and differentiation pathways of embryonic stem cells by forming an integral biological network with transcription factors.
In summary, our study using a combination of alignment, statistical and network analysis tools to examine deep sequencing data of microRNAs in hES-MSC has led to a result that (i) identifies intracellular and exosome microRNA expression profiles of hES-MSCwith a possible mechanism of miRNA mediated intercellular regulation by these cells and (ii) placed HNF4A within the cross roads of regulation by the let-7 family of microRNAs.
Small RNAs play a wide range of regulatory roles from degradation to translational silencing of messenger RNA. The most studied class of small regulatory RNA is microRNAs (miRNAs). miRNAs are about 22 nucleotides long and have been identified in animals, plants and viruses. Precursor miRNAs have stem loop structures that are cleaved by Drosha and Dicer forming mature functional miRNA molecule. By forming RNA-induced silencing complexes, miRNA can either cleave messenger RNA molecules or inhibit translation. Through such mechanisms, they are involved in various cellular processes including hematopoietic differentiation and cell cycle regulation [1, 2].
The amount of genetic information regulated post-transcriptionally by miRNAs is potentially huge. Computational and indirect evidences indicate that miRNAs might regulate up to a third of all genes making direct and indirect consequences of miRNA directed regulation significant . The resulting regulatory network is very often an extensive and complex one.
Deep sequencing provides a rapid and sensitive way of obtaining miRNA profiles expressed by human embryonic derived mesenchymal stem cells (hES-MSC). These cells have the ability to differentiate into multiple mesenchymal phenotypes, such as bone, cartilage, tendon and adipose tissue . This property with a broad distribution of sources makes MSC an attractive therapeutic target.
Despite this level of interest, a clear understanding of the factors involved in regulation of MSC remains rudimentary. Global gene expression analysis has revealed that MSC differentiation into specific mature cells types is a temporally controlled and regulated process [5, 6]. miRNAs provide an attractive mechanism for temporal regulation of mRNA translation and stability. Regulation of miRNA expression pattern can then be perceived as a novel regulatory network affecting cellular function. Organization of the hierarchical order of stem cell types based on the linkage of their functional characteristics to such regulatory elements might present a novel means to understand and eventually manipulate cell fate. This miRNA mediated regulation of stem cell differentiation however will act in concert with other methods of regulation of gene expression like transcription factors, epigenetic mechanisms etc. Therefore it is important to interpret the results of this study in the context of understanding miRNA regulation in concert with other regulatory mechanisms that control cellular differentiation.
Recent literature has also revealed that genetic exchange of mRNA and miRNA between cells can be accomplished through microvesicles or exosome mediated transfer . Microvesicles can be shed from surfaces of activated cells or derived from the endosomal membrane compartment after fusion of secretory granules with the plasma membrane where they exist as intraluminal membrane-bound vesicles. Embryonic stem cells are a source of such microvesicles. It is possible that embryonic stem cell derived microvesicles contain biologically active molecules that affect growth and cell fate decision of targeted cells. Similar vesicles released from human and murine mast cell lines were shown to contain miRNAs . Employing microvesicles for transfer of genetic material would be an efficient means for these intercellular communications . It is conceivable that MSC derived from embryonic stem cells also communicate with neighbouring cells using microvesicles. In this study, we performed deep sequencing of small RNAs to understand the expression of miRNAs in both microvesicles and intracellular environment of hES-MSC.
This convergence of deep sequencing technology and the potential regulatory roles of miRNA in stem cells provide many opportunities as well as challenges in terms of data mining methodology in uncovering biological significance. Bioinformatics tools are necessary to bridge the gap between raw sequencing data and biological significance of these regulatory RNAs. A fundamental difficulty lies in identifying miRNA targets as well as deciphering the exact mechanism of translational repression of mRNAs by these snippets. Considerable work has been put into computational prediction of miRNA targets. Many challenges still remain in integrating large sequencing datasets to obtain gene-level distinctions.
In our study, we present a novel workflow using existing tools available for understanding these sequencing data sets. Our workflow extracts salient features from sequencing data for visualization and focuses on potential relationships between miRNAs with gene targets. Applying our workflow for quantitative estimation of small RNAs in hES-MSC intra and extra cellular examples, we can develop a better understanding of exosome mediated miRNA involvement in the regulation of other cells in the hES-MSC microenvironment.
Distribution of miRNAs from Deep Sequencing showed distinct phases
Considering the initial uneven distribution of transcripts with low count number, it is a possibility that this large group of transcripts constitute a form of noise since this group of transcripts is expressed sporadically with low counts. The concept of noise here is distinct from noise arising from stochastic variations in the expression level of a given gene in individual cells. It is possible that these transcripts arose from imperfect fidelity of Pol II transcriptional machinery due to stochastic components governing fundamental interactions between molecules within the cell. Alternatively, it could also be an artefact inherent to sequencing.
Transcripts of the subsequent phase with even distribution can then be inferred to be another group of miRNAs that constitute a steady state expression of miRNA genes in MSC followed by the last phase consisting of miRNAs with large count numbers that contributes significantly to the unique phenotype exhibited by such cells. From a genome wide perspective, this representation and interpretation of data presents a possibility to distinguish the relative proportions of biologically significant miRNAs from potential noise. To corroborate such a view with a quantitative aspect, a statistical test based strategy was developed.
Statistical Method employed implies transcript count of 32 as the minimum threshold
High read counts of let-7 family miRNAs transcripts are present in both intra and extra cellular samples of hES-MSC
Let 7 family of miRNAs are expressed strongly by hES-MSC within the cells. This table lists common transcripts in the intracellular samples in the high category of more than 10,000 transcripts showing a predominance of let 7 family of miRNAs.
Intracellular Transcripts Number (Sample 2)
Intracellular Transcripts Number (Sample 1)
Let 7 family of miRNAs are secreted by hES-MSC in large numbers. This table lists common transcripts in the extracellular samples that belong to the mid-range category of between 32 to 10,000 transcripts showing a predominance of let 7 family of miRNAs. It is interesting to note that no transcripts with counts greater than 10,000 were observed in the extracellular space.
Transcript Number (Sample 1)
Transcript Number (Sample 2)
Transcript Number Sample 3
Complexity Reduction using Gene interaction Networks revealed similarity in topology that suggested downstream targets for let-7 family of miRNAs
miRNAs are known to down regulate gene expression. This can be done by direct mRNA cleavage, mRNA decay by deadenylation or translational repression. The greatest challenge besetting incorporation of miRNA regulation into known gene expression mechanism is the great difficulty in predicting mRNA targets of miRNAs. Although the binding event behind the translational repression and mRNA degradation is driven primarily by complementarity between miRNA and target sites, computationally predicted target interactions generally generate a large list of targets. One such prediction algorithm, TargetScan requires perfect complementarity to the seed region of miRNA and extends these regions to account for complementarity outside the region. Conservation criteria based on the presence of seed region in an island of conservation is then incorporated by using groups of orthologous 3' UTR as input data. All of the above aims to efficiently reduce the false positive rates. However, the large amount of predictive targets for the let-7 miRNA family constitutes a complexity that can be difficult to interpret and explore. Since the accuracy of prediction of target sequences is low, it is difficult to identify bona fide gene targets regulated by miRNAs using Targetscan alone. Also such an approach will not yield genes indirectly regulated by miRNAs. Therefore we decided to focus on networks generated by Targetscan predictions rather than individual genes.
We have developed a workflow for analyzing deep sequencing data focusing on miRNA. This workflow has revealed several trends in miRNA expression profile for hES-MSCs. The majority of sequenced transcripts were identified by mapping to known small RNA libraries and a considerable proportion of these were miRNAs. Several studies have showed that certain miRNAs are associated with specific stem cell types and those that are regulated during stem cell differentiation. Characterizing miRNA expression profile for the purpose of extracting biological function from sequencing data can be a challenge.
Raw sequencing data of our samples showed distinct distribution patterns across a large range of magnitude for transcripts counts. The necessity of a minimum threshold of transcript count for biological significance was addressed by adaptive thresholding exploiting KS statistics amongst biological replicates giving a lower bound value for transcript counts. Currently KS statistics are used primarily for comparing distributions and have been used in areas such as image processing for examining pixel distribution. In our study, we employed KS statistics to estimate the lower bound where distributions of transcripts peaks begin to differ. Using a transcript count of 32 as a lower bound, we were able to define a threshold for miRNA expression based on transcript abundance. Further comparisons (as shown in Figure 4) of top ranked genes between our sequenced data with microarray data showed that there is a considerable degree of overlap between the two especially in the intracellular samples. miRNAs from the extra-cellular samples form a subset of those expressed in the intracellular samples. Among the most abundantly expressed transcripts across both intra and extra-cellular environment (Tables 1 &2) the let-7 family of miRNAs was the only overexpressed family of known miRNAs.
Amongst other over expressed miRNAs, miR199b has transcription factors like SOX4 and has been shown to be involved in liver cancer and muscular dystrophy [10, 11]. miR143 has recently been shown to cooperate with miR 145 to regulate plasticity and cell fate determination of smooth muscle cells via regulation of expression of transcription factors like Kruppel like factor 4 (Klf-4), myocardin and Elk-1 .
let-7 mutant phenotype studied in C. elegans has been associated with a lack of terminal differentiation and ongoing cell proliferation, both of which are characteristics of stem cells and cancer. let-7 targets include cell cycle regulators such as CDC25A and CDK6 ; promoters of growth including RAS and c-myc [14, 15] and a number of early embryonic genes including HMGA2, Mlin-41 and IMP-1 [16, 17]. Further, let-7 also targets Dicer [18, 19] which is the protein responsible for miRNA maturation. Therefore it is possible that the let-7 family of miRNA acts as a master regulator of miRNA function. Since the let-7 family of miRNAs were abundantly expressed in MSC and given their central role in controlling cellular differentiation and miRNA regulation, we decided to focus on this family of miRNAs for further investigations.
miRNA mediated translational repression might play an important role in stem cell self-renewal. In our study, a high level of let-7 family of miRNA transcripts was predominant in both intra and extra cellular samples for our hES-MSC. This was true for other abundantly expressed miRNAs in the extra cellular space as well. We were unable to find miRNAs that were expressed only in the extra cellular space or were expressed more abundantly in the extra cellular space compared to the intra cellular space. This suggests that the process of secretion of miRNAs is most likely a passive one unlike that of proteins, where signal peptide containing proteins are transported out of the cell as soon as they are synthesized.
Moreover, the elevated presence of miRNAs in the extra cellular environment suggests a possible intercellular signaling mechanism. Recent literature has also supported that genetic exchange of mRNA and miRNA between cells can be accomplished through microvesicles transfer [7, 8]. miRNAs like miR 21 , miR 133 and miR 30  have been shown to play a role in cardiac remodelling in infracted hearts. Since miRNA can affect cardiac remodelling and these miRNAs are transported in micro vesicles by transplanted MSC (Lim SK personal communication), it is possible that transplanted MSC regulate remodelling of the surrounding myocardium through secreted miRNAs in addition to secreted cytokines and soluble factors.
HNF4A is a transcription factor that has also been shown to be essential for morphological and functional differentiation of hepatocyte and for liver morphogenesis . HNF4A is also one of the physiological factors in the liver to activate apoB gene expression at the AF-1 site . Elevated levels of apoB have been correlated with an increased risk of atherosclerosis and coronary heart disease . HNF4A is also involved in the regulation of serum lipid levels and is linked to elevated serum cholesterol and triglyceride levels in Finnish combined familial hyperlipidemia patients . However to date there is no evidence of miRNA regulation of HNF4A. None of the target prediction algorithms predict the regulation of HNF4A by let-7 family of miRNAs. There are no predicted miRNA binding sites in the untranslated or intronic regions of the HNF4A gene. However combining the networks generated by sequence alignment of expressed miRNAs and Targetscan, we predict that HNF4A is indirectly regulated by the let-7 family of miRNAs. We were able to confirm this indirectly by measuring the levels of let-7 family of miRNAs and HNF4A in undifferentiated MSC and HEPG2 cells. In undifferentiated MSC when let-7 miRNAs are highly expressed, expression of HNF4A is very low. Conversely in HEPG2 cells where a high level of HNF4A is expressed, we find very low expression of let-7 family miRNAs. Since HNF family of transcription factors have been reported to be upregulated in hepatocytes derived from adipose tissue MSC , it is possible that let-7 regulates HNF4A levels during this process. Since genes in the HNF4A alignment network also show a similar expression profile to HNF4A in HEPG2 and hES-MSC cells, it is possible that the let-7 family miRNA regulation of HNF4A is mediated through genes in this network. It is also interesting to note that there is no predicted miRNA target site in the UTR of HNF4A. Thus miRNA based regulation of HNF4A may be indirectly achieved through regulation of genes in the network that interact with HNF4A. It will be interesting to study the effect of Lin28 overexpression on differentiation of hepatocytes from hES-MSC since Lin 28 is a transcription factor that inhibits function of let-7 family miRNAs.
let-7f miRNA is known to exert pro-angiogenic effects . Human MSC conditioned medium has been shown to reduce infarct size in patients with acute myocardial infraction . Since MSC conditioned medium contains exosomes with let-7 family miRNAs and these let-7 family miRNAs may regulate HNF4A (based on our network and expression analysis), it is highly likely that MSC conditioned medium mediated reduction of infarct size is achieved by indirect regulation of HNF4A mediated by the let-7 family of miRNAs.
Thus, our study suggests the possibility of let-7 family of miRNAs indirectly regulating this particular transcription factor to achieve physiological changes. This is in line with recent paradigm where miRNAs regulate self-renewal and differentiation pathways of embryonic stem cells by forming an integral biological network with transcription factors .
In conclusion, our study using a combination of different available tools to examine deep sequencing data by examining alignment, computer predictions, mathematical and network analysis has led us to a hypothesis that HNF4A is indirectly regulated by the let-7 family of miRNAs.
hES-MSC intracellular samples and microvesicles isolation from culture media
Human embryonic derived mesenchymal stem cells were used in this study. Two intracellular samples were obtained at passage 18 by lysing the collected cells. RNA was isolated from conditioned medium by adding three volumes of Trizol LS (Invitrogen) to one volume of conditioned medium and completing the extraction according to the manufacturer's protocol. Total RNA and small RNAs from MSC were purified using Trizol (AppliedBiosystems) and mirVana™ miRNA Isolation Kit (AppliedBiosystems), respectively. RNA was quantiated using Quanti-T™ RiboGreen® RNA Assay Kit (Invitrogen). In a separate experiment, the conditioned medium was treated with RNase, RNA extracted and compared to RNA extracted from untreated conditioned medium. No difference was observed in the RNA profile from the RNase treated and untreated samples suggesting that the RNA in the conditioned medium was contained within the exosomes (data not shown). Resulting samples were sent to Illumina Sequencing facilities for Deep Sequencing. The intracellular samples were further gel purified to only include transcripts between 18-35 nucleotides long. This is a step in the standard sample purification protocol for Deep Sequencing. This step was omitted in the samples of hES-MSC conditioned medium profiling the extra cellular small RNAs. However since all miRNAs are in the 18-35 nucleotide range this step has no effect on the data analysis described here which considered only known miRNAs.
Extracting biologically significant miRNA profiles from raw Deep Sequencing transcripts counts between intra and extra-cellular samples through data visualization tools and plots
Deep sequencing of small RNAs was performed using the Illumina/Solexa platform (sequencing depth 5.5 million reads) which produces individual sequence reads with base quality scores. Identical sequences were counted and removed. The resulting dataset consisted of unique sequences with the associated read counts. Prior to further analysis as show in the workflow (Figure 1), the adaptor sequences were trimmed computationally. Resulting sequences were mapped against known small RNA libraries to identify sequences originating from sources such as rRNA, tRNA, snRNA and snoRNA. Focusing primarily on the dataset consisting of only known human miRNAs, only small RNA sequences that align to mirBASE were retained.
These sequences were aligned to the human genome using Seqmap  with a 3 nucleotide mis-match condition. We arrived at the 3 nucleotide mismatch condition after looking at the distribution of the aligned sequences at 1, 2, 3 and 4 nucleotide mismatch. The distribution changed significantly after 3 nucleotide mismatch. Therefore we concluded that this was the maximum allowable mismatch for aligning the sequenced transcripts. The Illumina ELAND alignment tool which yields similar results as Seqmap also uses a 3 nucleotide mismatch for aligning RNA seq data. The alignment results were visualized using UCSC genome browser . Genomic graphs depict the spatial distribution across the genome and are characterized by peaks along the genome loci corresponding to transcript count that align to that region. The magnitude of each peak is a reflection of the number of transcripts. Histograms were used to describe the frequency of occurrence of theses peaks for further analysis.
Determining minimum threshold count of biological significance
The process of quantifying the number and magnitude of peaks due to each miRNA necessitates a threshold level to differentiate miRNAs of biological significance from background noise. We quantified a threshold value by deploying a strategy via applying Kolomogorov-Smirnov (KS) tests to our data iteratively.
Kolomogorov-Smirnov test is a form of minimum distance estimation. In this aspect, we used it as a nonparametric test of equality for comparing distributions from our biological replicates. The KS test statistics quantifies the distance between the cumulative distributions of our replicates and thus can be viewed as a measure of agreement between them. The threshold number of transcripts was changed with each application of the KS test and the KS statistics were recorded. The final threshold count is the one which corresponds to the first instance of a minimum point in the KS statistics value.
This condition can be thought of as the first instance when the distributions of miRNA amongst replicates began to be similar. Such an approach can be thought of as a form of adaptive thresholding. The key parameter in a thresholding process is the choice of the threshold value. Usually, this value is arbitrary chosen. Our approach computes this value automatically considering the replicates of our experiments. Threshold values can thus be uniquely generated for each experimental protocol or cell type to differentiate biological significance from background noise. Observation of the distributions of transcripts from deep sequencing suggested the presence of different distributions within each of the intracellular samples. Application of KS statistics to our data ensures that further interpretations are based on reproducible parts of the overall distributions. KS statistics also provide a lower threshold boundary which provided a quantitative justification for selecting candidates for detailed analysis in later stages of the workflow .
Having determined the threshold, we then proceeded to characterize the salient features of miRNA expression profiles by ranking the top miRNAs that accounted for the large magnitude peaks observed in the genomic graphs. The top 100 ranked miRNAs were compared to a similar list derived from a miRNA microarray experiment of the same samples as a measure of consistency.
Network Analysis for predicting indirect miRNA gene Targets
TargetScan  was used to predict gene targets for the let-7 family of miRNAs. An alternative representation of the TargetScan results was used in our study where the hypothesis is that biological pathways and networks rather than the individual genes are driving development of the range of phenotypes observed. It is such networks and interactions that we would like to understand. Using such framework to conceptualize predicted miRNA gene targets from TargetScan, the targets for the let-7 family of miRNA were subjected to pathway exploration using the Ingenuity Pathway Analysis (Ingenuity® Systems, http://www.ingenuity.com) software. Using Ingenuity Pathway Analysis and accompanying interaction database, top ranking interaction networks were generated.
Our approach in reducing complexity of gene interactions networks uses the integration of the data from different sources. This was done by comparing the network generated based on our data from mapping the respective peaks along genomic location with those derived from the TargetScan prediction. Since Targetscan over predicts miRNA targets, we compared networks generated with genes predicted by Targetscan with networks generated by genes from our alignment data. We hoped to extract biological significance from the expressed miRNAs by sieving for common nodes in these networks of predicted miRNA targets. These common nodes of gene targets strongly suggest a down stream target being regulated by the respective miRNAs.
Measurements of expression levels of let-7 family of miRNAs and HNF4A by quantitative real time reverse transcriptase polymerase chain reaction (qRT-PCR)
RNA from cultured hES-MSC and HEPG2 cells was extracted using Trizol reagent (Invitrogen, Carlsbad, CA). qRT-PCR was performed using Taqman assays according to manufacturer's instructions. For calculating relative fold change values, the CT values were normalized to U6 miRNA as internal control for miRNA and GAPDH for HNF4A. We chose U6 as the internal control because it is similar in length and structure to known miRNAs.
Further we also verified the prediction that let-7 family miRNAs regulate the network of 50 genes, by examining the expression profiles of these genes in MSC and HEPG2 cells. The mRNA expression for 12 hES-MSC samples was profiled by our lab using the Illumina Ref8 v3 Expression BeadChip (data not shown). The mRNA expression profiles for two HepG2 control samples were obtained from ArrayExpress (ID: E-MEXP-1213). We chose these samples because the mRNA profiling was also performed on Illumina microarrays (HumanWG-6 v3), allowing for more comparable expression values. Expression values for the 50 genes identified in the HNF4A alignment network were extracted for each sample. The expression values were normalized by dividing them by the level of expression of a housekeeping gene, glyceraldehyde 3-phosphate dehydrogenase (GAPDH) in each sample. Then, the normalized expression value for each gene was averaged across the 12 hES-MSC samples, and across 2 HEPG2 samples. A ratio of averaged normalized gene expression value for each gene was generated by comparing expression in HEPG2 to hES-MSC cells. The relative expression for these 50 genes was visualized using Ingenuity Pathway Analysis.
This work was supported by the Biomedical Research Council of A*STAR (Agency for Science, Technology and Research), Singapore. The authors acknowledge Drs. Igor Kurochkin, Chandra Verma and Antonis Giannakakis for useful comments and recommendations. The authors would also like to acknowledge the help and constant support provided by the Illumina team, specifically Raymond McCarley and Nobuta Kan during the initial preprocessing of the sequencing data.
This article has been published as part of BMC Genomics Volume 11 Supplement 1, 2010: International Workshop on Computational Systems Biology Approaches to Analysis of Genome Complexity and Regulatory Gene Networks. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/11?issue=S1.
- Chen C-Z, Lodish HF: MicroRNAs as regulators of mammalian hematopoiesis. Seminars in Immunology. 2005, 17 (2): 155-165. 10.1016/j.smim.2005.01.001.View ArticlePubMed
- Chivukula RR, Mendell JT: Circular reasoning: microRNAs and cell-cycle control. Trends in Biochemical Sciences. 2008, 3 (10): 474-481. 10.1016/j.tibs.2008.06.008.View Article
- Rajewsky N: microRNA target predictions in animals. Nat Genet. 2006
- Lian Q, Lye E, Suan Yeo K, Khia Way Tan E, Salto-Tellez M, Liu TM, Palanisamy N, El Oakley RM, Lee EH, Lim B: Derivation of Clinically Compliant MSCs from CD105+, CD24- Differentiated Human ESCs. Stem Cells. 2007, 25 (2): 425-436. 10.1634/stemcells.2006-0420.View ArticlePubMed
- Doi M, Nagano A, Nakamura Y: Molecular cloning and characterization of a novel gene, EMILIN-5, and its possible involvement in skeletal development. Biochemical and Biophysical Research Communications. 2004, 313 (4): 888-893. 10.1016/j.bbrc.2003.11.181.View ArticlePubMed
- Qi H, Aguiar DJ, Williams SM, La Pean A, Pan W, Verfaillie CM: Identification of genes responsible for osteoblast differentiation from human mesodermal progenitor cells. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (6): 3305-3310. 10.1073/pnas.0532693100.PubMed CentralView ArticlePubMed
- Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO: Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007, 9 (6): 654-659. 10.1038/ncb1596.View ArticlePubMed
- Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, Xiao T, Schafer J, Lee M-LT, Schmittgen TD: Detection of microRNA Expression in Human Peripheral Blood Microvesicles. PLoS ONE. 2008, 3 (11): e3694-e3694. 10.1371/journal.pone.0003694.PubMed CentralView ArticlePubMed
- Simons M, Raposo G: Exosomes--vesicular carriers for intercellular communication. Curr Opin Cell Biol. 2009, 21 (4): 575-581. 10.1016/j.ceb.2009.03.007.View ArticlePubMed
- Eisenberg I, Eran A, Nishino I, Moggio M, Lamperti C, Amato AA, Lidov HG, Kang PB, North KN, Mitrani-Rosenbaum S: Distinctive patterns of microRNA expression in primary muscular disorders. Proceedings of the National Academy of Sciences. 2007, 104 (43): 17016-17021. 10.1073/pnas.0708115104.View Article
- Gramantieri L, Ferracin M, Fornari F, Veronese A, Sabbioni S, Liu C-G, Calin GA, Giovannini C, Ferrazzi E, Grazi GL: Cyclin G1 Is a Target of miR-122a, a MicroRNA Frequently Down-regulated in Human Hepatocellular Carcinoma. Cancer Res. 2007, 67 (13): 6092-6099. 10.1158/0008-5472.CAN-06-4607.View ArticlePubMed
- Cordes KR, Sheehy NT, White MP, Berry EC, Morton SU, Muth AN, Lee T-H, Miano JM, Ivey KN, Srivastava D: miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature. 2009, 460 (7256): 705-710.PubMed CentralPubMed
- Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, Wilson M, Wang X, Shelton J, Shingara J: The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Research. 2007, 67 (16): 7713-7722. 10.1158/0008-5472.CAN-07-1083.View ArticlePubMed
- Kumar MS, Lu J, Mercer KL, Golub TR, Jacks T: Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nature Genetics. 2007, 39 (5): 673-677. 10.1038/ng2003.View ArticlePubMed
- Sampson VB, Rong NH, Han J, Yang Q, Aris V, Soteropoulos P, Petrelli NJ, Dunn SP, Krueger LJ: MicroRNA let-7a down-regulates MYC and reverts MYC-induced growth in Burkitt lymphoma cells. Cancer Research. 2007, 67 (20): 9762-9770. 10.1158/0008-5472.CAN-07-2462.View ArticlePubMed
- Lee YS, Dutta A: The tumor suppressor microRNA let-7 represses the HMGA2 oncogene. Genes & Development. 2007, 21 (9): 1025-1030. 10.1101/gad.1540407.View Article
- Mayr C, Hemann MT, Bartel DP: Disrupting the pairing between let-7 and Hmga2 enhances oncogenic transformation. Science (New York, NY). 2007, 315 (5818): 1576-1579.View Article
- Forman JJ, Legesse-Miller A, Coller HA: A search for conserved sequences in coding regions reveals that the let-7 microRNA targets Dicer within its coding sequence. Proceedings of the National Academy of Sciences. 2008, 105 (39): 14879-14884. 10.1073/pnas.0803230105.View Article
- Tokumaru S, Suzuki M, Yamada H, Nagino M, Takahashi T: let-7 regulates Dicer expression and constitutes a negative feedback loop. Carcinogenesis. 2008, 29 (11): 2073-2077. 10.1093/carcin/bgn187.View ArticlePubMed
- Thum T, Gross C, Fiedler J, Fischer T, Kissler S, Bussen M, Galuppo P, Just S, Rottbauer W, Frantz S: MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts. Nature. 2008, 456 (7224): 980-984. 10.1038/nature07511.View ArticlePubMed
- Duisters RF, Tijsen AJ, Schroen B, Leenders JJ, Lentink V, Made van der I, Herias V, van Leeuwen RE, Schellings MW, Barenbrug P: miR-133 and miR-30 regulate connective tissue growth factor: implications for a role of microRNAs in myocardial matrix remodeling. Circulation Research. 2009, 104 (2): 170-178. 10.1161/CIRCRESAHA.108.182535. 176p following 178-170-178, 176p following 178.View ArticlePubMed
- Spath GF, Weiss MC: Hepatocyte Nuclear Factor 4α Provokes Expression of Epithelial Marker Genes, Acting As a Morphogen in Dedifferentiated Hepatoma Cells. J Cell Biol. 1998, 140 (4): 935-946. 10.1083/jcb.140.4.935.PubMed CentralView ArticlePubMed
- Wang Z, Burke PA: Modulation of hepatocyte nuclear factor-4alpha function by the peroxisome-proliferator-activated receptor-gamma co-activator-1alpha in the acute-phase response. Biochem J. 2008, 415 (2): 289-296. 10.1042/BJ20080355.PubMed CentralView ArticlePubMed
- Metzger S, Halaas JL, Breslow JL, Sladek FM: Orphan receptor HNF-4 and bZip protein C/EBP alpha bind to overlapping regions of the apolipoprotein B gene promoter and synergistically activate transcription. J Biol Chem. 1993, 268 (22): 16831-16838.PubMed
- Weissglas-Volkov D, Huertas-Vazquez A, Suviolahti E, Lee J, Plaisier C, Canizales-Quinteros S, Tusie-Luna T, Aguilar-Salinas C, Taskinen M-R, Pajukanta P: Common hepatic nuclear factor-4alpha variants are associated with high serum lipid levels and the metabolic syndrome. Diabetes. 2006, 55 (7): 1970-1977. 10.2337/db06-0035.View ArticlePubMed
- Yamamoto Y, Banas A, Murata S, Ishikawa M, Lim CR, Teratani T, Hatada I, Matsubara K, Kato T, Ochiya T: A comparative analysis of the transcriptome and signal pathways in hepatic differentiation of human adipose mesenchymal stem cells. FEBS J. 2008, 275 (6): 1260-1273. 10.1111/j.1742-4658.2008.06287.x.View ArticlePubMed
- Kuehbacher A, Urbich C, Zeiher AM, Dimmeler S: Role of Dicer and Drosha for Endothelial MicroRNA Expression and Angiogenesis. Circ Res. 2007, 101 (1): 59-68. 10.1161/CIRCRESAHA.107.153916.View ArticlePubMed
- Timmers L, Lim SK, Arslan F, Armstrong JS, Hoefer IE, Doevendans PA, Piek JJ, El Oakley RM, Choo A, Lee CN: Reduction of myocardial infarct size by human mesenchymal stem cell conditioned medium. Stem Cell Research. 2007, 1 (2): 129-137. 10.1016/j.scr.2008.02.002.View ArticlePubMed
- Tsang J, Zhu J, van Oudenaarden A: MicroRNA-Mediated Feedback and Feedforward Loops Are Recurrent Network Motifs in Mammals. Molecular Cell. 2007, 26 (5): 753-767. 10.1016/j.molcel.2007.05.018.PubMed CentralView ArticlePubMed
- Jiang H, Wong WH: SeqMap: mapping massive amount of oligonucleotides to the genome. Bioinformatics. 2008, 24 (20): 2395-2396. 10.1093/bioinformatics/btn429.PubMed CentralView ArticlePubMed
- Karolchik D, Baertsch R, Diekhans M, Furey TS, Hinrichs A, Lu YT, Roskin KM, Schwartz M, Sugnet CW, Thomas DJ: The UCSC Genome Browser Database. Nucl Acids Res. 2003, 31 (1): 51-54. 10.1093/nar/gkg129.PubMed CentralView ArticlePubMed
- Kuznetsov VA, Orlov YL, Wei CL, Ruan Y: Computational analysis and modeling of genome-scale avidity distribution of transcription factor binding sites in chip-pet experiments. Genome Informatics International Conference on Genome Informatics. 2007, 19: 83-94. full_text.PubMed
- TargetScanHuman 5.1. [http://www.Targetscan.org/vert_50/]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.