PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2014; 4: 5150.
Published online 2014 June 3. doi:  10.1038/srep05150
PMCID: PMC5381490

Genome-wide survey of tissue-specific microRNA and transcription factor regulatory networks in 12 tissues

Abstract

Tissue-specific miRNAs (TS miRNA) specifically expressed in particular tissues play an important role in tissue identity, differentiation and function. However, transcription factor (TF) and TS miRNA regulatory networks across multiple tissues have not been systematically studied. Here, we manually extracted 116 TS miRNAs and systematically investigated the regulatory network of TF-TS miRNA in 12 human tissues. We identified 2,347 TF-TS miRNA regulatory relations and revealed that most TF binding sites tend to enrich close to the transcription start site of TS miRNAs. Furthermore, we found TS miRNAs were regulated widely by non-tissue specific TFs and the tissue-specific expression level of TF have a close relationship with TF-genes regulation. Finally, we describe TSmiR (http://bioeng.swjtu.edu.cn/TSmiR), a novel and web-searchable database that houses interaction maps of TF-TS miRNA in 12 tissues. Taken together, these observations provide a new suggestion to better understand the regulatory network and mechanisms of TF-TS miRNAs underlying different tissues.

Gene expression in metazoans is largely controlled by various trans-regulatory factors at various levels. At the transcriptional level, transcription factor (TF) has been considered as the primary regulator to control gene expression. By binding to specific sequences usually located in promoter regions (also known as TF binding sites, TFBS), TF can activate or repress transcription of their target genes, and form transcriptional regulatory networks1,2. In recent years, the emergence of miRNAs as another crucial suppressive regulator which share the similar regulatory logic with TFs has occurred3. MicroRNAs (miRNAs) are a class of short non-coding RNAs of 18 to 24 nucleotides in length that post-transcriptionally regulate various genes through direct degradation of the target mRNA and/or translational repression4,5,6. Abundant evidence demonstrates that miRNA control a variety of biological processes such as cell cycle, differentiation, cell proliferation, and apoptosis7,8.

Tissue-specific patterns of gene expression play fundamental roles in tissue development, distinctive features of cell types, function, and, transcriptional regulation9. Among the identified miRNAs, some of them exhibited tissue-specific or developmental-stage-specific expression pattern and contributed potential roles in maintaining tissue identity and function10,11. Tissue-specific miRNAs (referred to as TS miRNA) have been reported to be associated with various human diseases such as cardiovascular disease, diabetes and cancer12,13,14. Moreover, it has been proposed that tissue-specific gene expression patterns are controlled by combinations of TF transcriptional regulatory networks15,16. Therefore, the study of regulatory networks composed of tissue-specific miRNAs and TFs is necessary to understand tissue specificity regulation and function.

In recent years, genome-wide identification of TF-miRNA regulatory networks have been extensively studied based on the fact that TFs can regulate miRNA transcription by binding to the promoter regions of miRNA17,18,19,20. Most of these studies focus on a single tissue, or consider various tissues as a whole. However, different tissues possess a different regulatory network to perform particular functions in corresponding tissue. Therefore, systematically mapping combinatorial regulatory networks among TFs and miRNAs, especially tissue-specific miRNA and TFs across different tissues would represent a significant leap forward in disclosing molecular basis of tissue-specific gene expression, development, function and how tissue specificity is determined. In the present study, 116 experimentally validated tissue-specific miRNAs (TS miRNAs) were extracted from literatures and qRT-PCR data. It was found that half of TS miRNAs were clustered miRNAs. Here, 2,347 TF-TS miRNA regulatory relationships were identified using the TF ChIP-seq data from the ENCODE (The Encyclopedia of DNA Elements) project21 which provided high-resolution TFBS in multiple cell lines. Also, it was found that most TF-TS miRNA regulation occurred across multiple tissues, and most TFBS tend to enrich close to the transcription start site (TSS). Through the integration of TF expression data, we found that tissue-specific miRNAs were regulated widely by non-tissue specific TFs. In addition, 90 TF-TS miRNA regulatory relations were found which TF and TS miRNA specifically expressed in the same tissue; these regulatory relations suggest they perform their specific effect in a particular tissue. Furthermore, a series of TF-TS miRNA regulatory networks presented here, revealing that TF-gene regulatory relationships in network displayed two distinct types: type I : The highly tissue-specific or widely expressed TFs make less intensive interactions than other TFs; and type II: tissue-specific TFs participate in more interactions than other non-specific expression TFs. Finally, the TSmiR database was presented here (http://bioeng.swjtu.edu.cn/TSmiR) to provide interaction maps and expression data of transcription factor and tissue-specific miRNAs in 12 tissues. To our knowledge, this is the first systematic attempt to construct a regulatory network of tissue-specific miRNAs across multiple tissues, which can help to elucidate the molecular mechanisms of tissue-specific miRNAs network in tissue development and function.

Results and discussion

Identification of TS miRNAs

Using hand-curated screenings of qRT-PCR data from miRNAMap22 databases and literatures, 116 experimentally validated tissue-specific miRNAs from heart, skeletal muscle, lung, bone, kidney, liver, placenta, testis, brain, spleen, thymus and pancreas were retrieved. Twelve of 116 TS miRNAs (miR-1, miR-126, miR-208, miR-128a, miR-133a, miR-133b, miR-134, miR-146a, miR-377, miR-483, miR-92a and miR-95) were specifically expressed in two tissues; whereas, the remaining TS miRNAs were specifically expressed in only one tissue. After carrying out the clustering analysis for all tissue-specific miRNAs according to their genome context, half of them were found to be clustered miRNAs. Especially for the placenta, except miR-136, miR-184 and miR-381, 40 out of 43 (93.02%) placenta-specific miRNAs were clustered miRNAs, suggesting that these clustered miRNAs have similar expression patterns and play related or particular functions in this tissue. Other tissues also showed a high proportion of clustering miRNAs, for example, 7 out of 10 (70%) clustered miRNAs in brain, and 9 out of 14 (64.29%) clustered miRNAs in testis (see Table 1).

Table 1
Identification and cluster analysis of tissue-specific miRNAs in 12 tissues

Genome-wide identification of TF-TS miRNA regulation

Previous studies show that transcription factors, considered to regulate the transcription of microRNA, have similar mechanisms to those of protein-coding genes23. Therefore, the traditional approach of researching TF-gene regulatory is used to determine TF-miRNA relations. We first identified the TS miRNA TSS from high-throughput experimental data in four literature sources11,24,25,26. Next, 5 kb upstream and 1 kb downstream of each TS miRNA TSS were defined as the putative transcription factor binding region (promoter region) for each TS miRNAs based on preceding studies27,28. Previously, most studies usually produced high false positive results because position weight matrix (PWM) was used to predict the TFBS in miRNA putative prompter region29. Recently, ChIP-seq has emerged as a valuable and high-resolution approach to profile TFBSs in vivo30. Particularly, the Encyclopedia of DNA Elements (ENCODE) project has recently provided high-resolution TFBS in multiple cell lines to help researchers identify TF-miRNA regulation more effectively. Therefore, 2,347 TF-TS miRNA regulatory relationships were explored using highly conservative TFBS from ENCODE TF ChIP-seq data to scan the putative TS miRNA promoter region (see Supplementary Table S1).

To determine the distribution of 2,347 TF-TS miRNA regulatory relationships in 12 tissues, statistical analysis of the number of tissues in which TF-TS miRNA interactions can occur have been explored (see Supplementary Table S2).The result shows that most TF-TS miRNA regulatory relationships occurred in 5–9 tissues, in other words, most TS miRNAs regulated by a TF appeared in 5–9 tissues. The number of TF-TS miRNA relations decreased, accompanied by less or greater number of tissues in which TF-TS miRNA interactions can occur (Figure 1a). It suggests that TF-TS miRNA interactions which occur only in one tissue perform a specific biological function such as, tissue development and specific regulation in a particular tissue. Conversely, those TF-TS miRNA interactions which occur in the majority of tissues may possibly prefer performing a wider range of biological functions. Furthermore, the proportionate number of TF-TS miRNA interactions that occur in 1–12 tissues compared with the total TF-TS miRNA interactions of each tissue, was calculated. As shown in Figure 1b,it is apparent that non-specific distribution of regulatory relations of TF-TS miRNA are observed more in lung and pancreas tissue than other tissues, suggesting that TF-TS miRNA regulatory relationships in lung and pancreas are non-specific across 12 tissues. The largest proportion for specific distribution of TF-TS miRNA interactions is presented in placenta, this phenomenon may be related with the existence of a large number of clustered miRNAs which TF co-regulate through binding their shared promoter region.

Figure 1
The proportion of TF-TS miRNA regulation in 12 tissues.

TF binding profiles around TSS of TS miRNA

To determine TF binding profile around TSS (−5 kb ~ 1 kb) of TS miRNA, the pattern of TF occupancy was analyzed in each tissue. The number of TFs which regulate pancreas-specific miRNAs were too few, thus the TF binding profile of 11 tissues was examined (see Fig. 2). It was found that most TFBSs tend to enrich close to the TSS of TS miRNA (approximately −1 kb ~ 0.4 kb), which is consistent with previous work31. Significantly, for testis, heart, placenta and lung tissue-specific miRNA, there are some another enriched TFBS regions, which is located in the upstream regions of TSS of TS miRNA (−3.4 ~ −3 kb in testis, −4.8 ~ −4.4 kb in heart, −4.4 ~ −4 kb in placenta and −4 ~ −3.2 kb in lung, respectively). These TF-binding loci further away from the TSS may represent distal cis-regulatory elements for precisely regulated TS miRNA and represent involvement in tissue-specific gene expression32.

Figure 2
Heat map of TF binding profiles occupancy in −5 kb ~ 1 kb around the TSS of TS miRNA.

The distribution of tissue-specific TFs

To determine the regulatory relationship between tissue-specific miRNA and tissue-specific TFs, the following studies were done: Firstly, the tissue-specific value (TSPV) of TFs which involved in 2,347 TF-TS miRNA regulatory relationships was calculated (see method). The lower the TSPV represents the stronger the tissue specificity. For a particular tissue, the tissue-specific value in a tissue (TSVT) determines the specific expression level of a TF; the greater the value of TSVT suggests a TF is more specific in a tissue. We set TSVT > −2.5 as threshold value to indicate a TF is tissue-specific in a tissue according to the distribution of TSVT (see Supplementary Table S1 and Fig. S1).

Four tissue-specific expression levels of TFs were defined according to its TSPV: high tissue specific (tissue-specific value > = −160 and <−120, medium tissue specific (tissue-specific value > = −120 and <−80), low tissue specific (tissue-specific value > = −80 and <−60) and non-tissue-specific (tissue-specific value > = −60). For example, transcription factor CTCFL, a well-known testis-specific TF33, the TSPV of CTCFL is −107.2826303, is conservative according to our definition. CTCFL is a medium tissue specific TF (TSPV > = −120 and <−80), which with equivalence to the copy number of CTCFL in testis (442.6120899) is 288-fold greater than the mean copy number (1.53577) of TF in the other 11 tissues. Next, the proportionate number of TFs, classified by four tissue-specific levels, compared to the total number of TFs in each tissue was carried out. Results revealed that the majority of TFs involved in 2,347 TF-TS miRNA regulation relationships, represent non-tissue specific expression in a particular tissue, and the TSPV is mainly distributed from the range −60 to −43.02 (Fig. 3). Subsequently, the proportion of four tissue-specific expression levels of TFs across 12 tissues was explored. As shown in Figure 4, it is apparent that tissue-specific miRNAs were regulated widely by non-tissue specific TFs. MiRNAs are transcribed by RNAP II, which suggests that miRNA are regulated in a similar fashion as protein-coding genes34. This result was consistent with previous observations that most TF involved in tissue-specific TF-TF regulatory networks were expressed non-specifically in the corresponding tissue35. The finding indicates that TS miRNA perform specific functions in a tissue, such as tissue development and identity, mainly through the regulation of multiple signalling pathways instead of tissue-specific pathways. On the contrary, high or medium tissue-specific TFs suggest these TFs play a role in special functional regulation together with miRNAs in corresponding tissue.

Figure 3
Scatterplots of TSPV of TFs.
Figure 4
The proportion of four tissue-specific expression levels of TFs.

Tissue-specific TF-miRNA regulation

The exploration of the regulatory relationships between TFs and miRNAs both specifically expressed in the same tissue, could offer useful information to elucidate how TF-miRNA regulation plays a particular role in tissue specification or cell differentiation36. Therefore, the TF-miRNA regulation that is: TF and miRNA both specifically expressed in the same tissue, was screened for (Table 2). Finally, it was found that 38 TFs were involved in 90 TF-TS miRNA regulation in bone, brain, kidney, liver, lung, placenta, skeletal muscle, spleen, testis and thymus. For example, skeletal muscle-specific expression TF, serum response factor (SRF), regulates miR-1 cluster (miR-1 and miR-133) which is well known to express specifically in skeletal muscle; this signaling pathway has been certified to play a critical role in modulating skeletal muscle proliferation and differentiation37. The exploring of TF and TS miRNA which are specifically expressed in the same tissue will help further experimental validation studies to clarify these consistent tissue-specific TF-TS miRNA regulation and how they perform their effects in tissue specific manners.

Table 2
TF-TS miRNA regulation which TF and miRNA both specifically expressed in the same tissuea

The identification and expression analysis of TS miRNAs target genes in 12 tissues

The function of TS miRNA is achieved mainly through the miRNA target genes at the post-transcriptional level; therefore, exploring TS miRNA-target gene pairs will help researchers to further study the regulations and functions of TS miRNA in tissue specification, physiologies, differentiation, development, etc. In order to obtain highly reliable TS miRNA target genes, the experimentally verified target genes were downloaded from miRTarBase and miRecords database. If there were no experimentally verified TS miRNA target genes, TargetScan Human was used to predict the target genes and a strict threshold filter was set to ensure the reliability of prediction (see method). Finally, 3,299 TS miRNA target genes in 12 tissues were obtained: 1,419 experimentally verified target genes and 1,880 predicted target genes. Furthermore, the functional annotation of experimentally verified TS miRNA target genes were provided using annotation tools (Supplementary Table S3).The results show that these target genes were involved in various biological processes and pathways to perform a wide range of biological functions, rather than limited to the tissue-specific functions. To determine whether TS miRNA target genes were significantly expressed specifically or not in corresponding tissue, the TS miRNA target genes were injected into the TiGER (Tissue-specific Gene Expression and Regulation) database, a comprehensive human tissue-specific gene expression database. Furthermore, the ratio of TS miRNA target genes expressed specifically in a tissue compared to the sum total number of TS miRNA target genes in a corresponding tissue, was calculated.

The results show that most of the TS miRNA target genes specifically expressed in multiple tissues (Fig. 5). However, some TS miRNAs target genes specifically expressed in the same tissue with TS miRNA (Table 2). These TS miRNA-target gene pairs suggest that they form various networks to regulate the physiologies, differentiation, development and specification in particular tissues38,39. Significantly, bone-specific miRNA target genes only exist in bone-specific expression genes, not in the presence of non-bone tissue-specific genes. To determine if TS miRNA target genes were expressed significantly specifically in corresponding tissue, the Fisher's exact test was performed (Supplementary Table S4). It was found that kidney and testis specific miRNA target genes significantly specifically expressed in corresponding tissues (P-value < 0.05). However, TS miRNA target genes in other tissues did not have significant specific expression in corresponding tissues. Results suggest that most TS miRNAs are involved in different biological functions by regulating a variety of target genes in different tissues and cell types. On the contrary, the TS miRNA and their target genes which are specifically expressed in the same tissue may play an essential role in the maintenance of a specific function, such as tissue identity and differentiation in a particular tissue.

Figure 5
The proportion of TS miRNA target genes which were specifically expressed in 12 tissues.

TF-TS miRNA regulatory network

TF and miRNA, as crucial trans-regulatory factors, have been considered to play an important role in controlling gene regulation at the transcription and posttranscriptional level. Recently, the transcriptional regulatory networks of TF-miRNA have been extensively explored; however, previous studies about TF-miRNA network did not consider the properties of tissue specificity of miRNA or expression level of transcription factors. Here, we presented a series of TF-miRNA regulatory networks that integrate verified or predicted interactions and expression data in 12 tissues and counted the number of regulatory relationships between TF, TS miRNA and target genes (see Fig. 6 and Table 3). The integrated network of 12 tissues contains 5,700 protein-protein interactions, 4,203 TF-target genes, 3,299 TS miRNAs-target genes and 2,347 TF-TS miRNAs (Table 3). In addition, to make user view network more clear, the high resolution original file (Cytoscape format, cys file) of microRNA and TF regulatory networks are provided for full exploration on our website (http://bioeng.swjtu.edu.cn/TSmiR/download.asp). Users can zoom in/out and pan for browsing the network. The networks reveal that TF-genes regulation in particular tissues showed two distinct types: type I, the highly tissue-specific or widely expressed TFs make less intensive interactions than other expression levels of TFs. For example, networks in the spleen, kidney, brain, heart, placenta, skeletal muscle. Type II: widely expressed TFs have more interactions than other TFs in the bone and liver. The regulatory network shows the different cell type will form significant molecular interactions preference and network characteristics based on the expression level of each TF. These results should prove a clue to clarify how TFs, miRNAs and target genes are coordinated to perform specific and common functions in different cell types.

Figure 6
Regulatory network among TF, TS miRNA and target genes was constructed using Cytoscape.
Table 3
Statistics of the number of interactions between TF, TS miRNA and target genes

TS miRNA database

Finally, here TSmiR (http://bioeng.swjtu.edu.cn/TSmiR), a free, web-accessible database, which provides information on interaction maps of transcription factor and TS miRNAs from experimentally validated and predicted data, was presented. It currently covers 116 TS miRNAs, 101 transcription factors and 2,347 TF-miRNAs regulatory relations in 12 tissues. Furthermore, experimentally validated expression data of TF and TS miRNA was also collected. The user can use the “search-by keyword” or “search-by category” function to retrieve the TF-TS miRNA regulatory relations. In addition to browsing TSmiR, there is a “browse” button at the top of the web page which allows users to explore TSmiR by clicking 12 different tissues (Fig. 7).

Figure 7
An overview of the query interface of TSmiR, which consists of “search by keyword” and “search by category” items.

Methods

Identification of tissue-specific miRNA

To screen the tissue-specific miRNA, the following work was done: 1) the experimentally validated tissue-specific miRNAs were collected from publications11,36,39,40,41,42,43,44,45,46,47,48,49,50,51; 2) the miRNA expression data was downloaded from miRNAMap22 and tissue-specific miRNA was screened according to the copy number of miRNA in a specific tissue is 80-fold greater than the mean copy number of miRNA in other tissues. Our screening criteria is more stringent than the definition of tissue-specific miRNA (microRNA expression in a tissue is 20-fold or higher compared with the mean of microRNA expression in other tissues)40. According to a previous study19, a group of miRNAs that are consecutively located within 10 kb of distance on the same genomic strand were defined as a miRNA cluster.

Identifying transcription start site (TSS) and promoter of TS miRNA

The TS miRNA TSS from high-throughput experimental data from four literature11,24,25,26 sources was identified: If miRNA did not have available TSS experimental information, the miRNA putative TSS was identified according to start site of each miRNA cluster. Next, the 5 kb upstream and 1 kb downstream of each TS miRNA TSS was identified as the putative transcription factor binding region (promoter region) for each TS miRNA based on previous studies27,28. At last, UCSC liftOver tool was used to convert old assembly to the current genome build (GRCh37/hg19).

Genome-wide identification of TF-TS miRNA regulation

The highly conservative TFBS generated by ChIP-seq of ENCODE project from UCSC database was downloaded. This includes all the TFBS data extracted from the “Txn Factor Chip” track which combines TFBS from many various cell lines via the UCSC Table Browser. Thus, the TFBS were used to scan the putative miRNA promoter region to identify TF-TS miRNA regulation.

TF binding profiles around TSS of TS miRNA

Genomic regions from 5 kb upstream to 1 kb downstream of the TS miRNA TSS were binned into 200-bp segments and the number of TFBS were calculated for each bin according to the overlap with the bin, respectively. Then, the ratio of the number of TFBS overlapped with each bin compared to total number of TFBS of each tissue was calculated. Heat maps of TF binding occupancy pattern around TSSs were generated with Cluster and TreeView software using the data produced above.

The distribution of tissue-specific TF

The human TF quantitative RT-PCR data was downloaded from the Ravasi et al. study52. Tissue-specific values (TSPV) of 101 TFs in 12 tissues were calculated according to follow equation (1) (a simplified formula from Ravasi study52):

An external file that holds a picture, illustration, etc.
Object name is srep05150-m1.jpg

where An external file that holds a picture, illustration, etc.
Object name is srep05150-m2.jpg represents the ratio of expression level of TF a in tissue b to sum total expression value across 12 tissues. For 12 tissues, the smaller TSPV means that TF expresses more specifically in particular tissues; whereas, TSPV approximately equal to −43.02 (the maximum value) means TF is expressed uniformly across 12 tissues. For a particular tissue, we used the An external file that holds a picture, illustration, etc.
Object name is srep05150-m3.jpg (tissue-specific value in a tissue, TSVT) to indicate which TFs are specifically expressed in this tissue. The greater the value of TSVT (maximal TSVT close to 0) suggests that TF is more specifically expressed in this tissue.

The identification of TS miRNA target genes

Also, the experimentally verified 2,854 TS miRNA target genes from miRTarBase53 (release version 2.5) were downloaded, along with 5,227 validated TS miRNA target genes from miRecords54 (release version 3), respectively. If there was no experimentally verified TS miRNA target genes, TargetScan Human (release 6.2)55 were used to predict the target genes of tissue-specific miRNAs. The conserved targets were downloaded from the TargetScan and were filtered for total context score <−0.3 before further analyses. The GO term (biological process; cellular component; molecular function), Entrez gene and KEGG pathway annotation were performed using the DAVID functional annotation table tool.

The expression analysis of TS miRNA target genes in corresponding tissue

In order to determine if TS miRNA target genes are specifically expressed in corresponding tissue, the TS miRNA target genes were injected into the TiGER56, a human tissue-specific gene expression database (http://bioinfo.wilmer.jhu.edu/tiger/). Firstly, for each tissue, the numbers of TS miRNA target genes specific expression in each of the 12 tissues respectively were counted, and then each of these enrichment numbers were divided by the number of TS miRNA target genes in the corresponding tissue for percentage. In order to find significant tissue-specific enriched target genes in the corresponding tissue, fisher's exact test between 12 tissues was carried out.

TF-TS miRNA regulatory network

Finally, the 622,751 human protein-protein interaction data from BIOGRIDE57 database was downloaded (release version 3.2.96). Following which, the GREAT28 was used to predict the target genes of 101 TFs in 12 tissues, settings used as follows: Species Assembly: Human GRCh37; Gene regulatory domain: 5 kb upstream and 1 kb downstream of TSS. The network of TF, TS miRNA and target genes was constructed by cytoscape58 software (version 2.8.3).

Supplementary Material

Supplementary Information:

Supplementary Figure S1 and Table S4

Supplementary Information:

Supplementary Table S1

Supplementary Information:

Supplementary Table S2

Supplementary Information:

Supplementary Table S3

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No.31200999), and the Fundamental Research Funds for the Central Universities (No. 2682014CX049).

Footnotes

The authors declare no competing financial interests.

Author Contributions Z.Y.G. and L.L.X. designed the experiments. Z.Y.G., M.M. and Y.L.Y. prepared and analyzed the data. R.F.D. and B.Z. constructed the database. G.Z.Y., M.M. and L.L.X. wrote the manuscript. All authors have read and approved the final manuscript.

References

  • Lee T. I. et al. Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 298, 799–804 (2002). [PubMed]
  • Yu H. & Gerstein M. Genomic analysis of the hierarchical structure of regulatory networks. Proc Natl Acad Sci U S A. 103, 14724–14731 (2006). [PubMed]
  • Hobert O. Gene regulation by transcription factors and microRNAs. Science. 319, 1785–1786 (2008). [PubMed]
  • Saito Y., Suzuki H. & Hibi T. The role of microRNAs in gastrointestinal cancers. J Gastroenterol. 44 Suppl 19, 18–22 (2009). [PubMed]
  • Baek D. et al. The impact of microRNAs on protein output. Nature. 455, 64–71 (2008). [PMC free article] [PubMed]
  • Farh K. K. et al. The widespread impact of mammalian MicroRNAs on mRNA repression and evolution. Science. 310, 1817–1821 (2005). [PubMed]
  • Bartel D. P. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 116, 281–297 (2004). [PubMed]
  • Filipowicz W., Bhattacharyya S. N. & Sonenberg N. Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight? Nat Rev Genet. 9, 102–114 (2008). [PubMed]
  • Niehrs C. & Pollet N. Synexpression groups in eukaryotes. Nature. 402, 483–487 (1999). [PubMed]
  • Ambros V. The functions of animal microRNAs. Nature. 431, 350–355 (2004). [PubMed]
  • Landgraf P. et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 129, 1401–1414 (2007). [PMC free article] [PubMed]
  • Calin G. A. et al. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A. 101, 2999–3004 (2004). [PubMed]
  • van Rooij E. et al. A signature pattern of stress-responsive microRNAs that can evoke cardiac hypertrophy and heart failure. Proc Natl Acad Sci U S A. 103, 18255–18260 (2006). [PubMed]
  • Esquela-Kerscher A. & Slack F. J. Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 6, 259–269 (2006). [PubMed]
  • Odom D. T. et al. Core transcriptional regulatory circuitry in human hepatocytes. Mol Syst Biol. 2, 2006 0017 (2006). [PMC free article] [PubMed]
  • Tan K., Tegner J. & Ravasi T. Integrated approaches to uncovering transcription regulatory networks in mammalian cells. Genomics. 91, 219–231 (2008). [PubMed]
  • Liao Y. L. et al. Transcriptional regulation of miR-196b by ETS2 in gastric cancer cells. Carcinogenesis. 33, 760–769 (2012). [PMC free article] [PubMed]
  • Re A., Cora D., Taverna D. & Caselle M. Genome-wide survey of microRNA-transcription factor feed-forward regulatory circuits in human. Mol Biosyst. 5, 854–867 (2009). [PMC free article] [PubMed]
  • Shalgi R., Lieber D., Oren M. & Pilpel Y. Global and local architecture of the mammalian microRNA-transcription factor regulatory network. PLoS Comput Biol. 3, e131 (2007). [PubMed]
  • Saini H. K., Griffiths-Jones S. & Enright A. J. Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci U S A. 104, 17719–17724 (2007). [PubMed]
  • Gerstein M. B. et al. Architecture of the human regulatory network derived from ENCODE data. Nature. 489, 91–100 (2012). [PMC free article] [PubMed]
  • Hsu S. D. et al. miRNAMap 2.0: genomic maps of microRNAs in metazoan genomes. Nucleic Acids Res. 36, D165–169 (2008). [PMC free article] [PubMed]
  • Megraw M. et al. MicroRNA promoter element discovery in Arabidopsis. RNA. 12, 1612–1619 (2006). [PubMed]
  • Chien C. H. et al. Identifying transcriptional start sites of human microRNAs based on high-throughput sequencing data. Nucleic Acids Res. 39, 9345–9356 (2011). [PMC free article] [PubMed]
  • Corcoran D. L. et al. Features of mammalian microRNA promoters emerge from polymerase II chromatin immunoprecipitation data. PLoS One. 4, e5279 (2009). [PMC free article] [PubMed]
  • Ozsolak F. et al. Chromatin structure analyses identify miRNA promoters. Genes Dev. 22, 3172–3183 (2008). [PubMed]
  • Alexiou P. et al. miRGen 2.0: a database of microRNA genomic information and regulation. Nucleic Acids Res. 38, D137–141 (2010). [PMC free article] [PubMed]
  • McLean C. Y. et al. GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 28, 495–501 (2010). [PMC free article] [PubMed]
  • Bi Y., Kim H., Gupta R. & Davuluri R. V. Tree-based position weight matrix approach to model transcription factor binding site profiles. PLoS One. 6, e24210 (2011). [PMC free article] [PubMed]
  • Johnson D. S., Mortazavi A., Myers R. M. & Wold B. Genome-wide mapping of in vivo protein-DNA interactions. Science. 316, 1497–1502 (2007). [PubMed]
  • Tabach Y. et al. Wide-scale analysis of human functional transcription factor binding reveals a strong bias towards the transcription start site. PLoS One. 2, e807 (2007). [PMC free article] [PubMed]
  • Franco L. O., de O. M. C. L., Hamdi S., Sachetto-Martins G. & de Oliveira D. E. Distal regulatory regions restrict the expression of cis-linked genes to the tapetal cells. FEBS Lett. 517, 13–18 (2002). [PubMed]
  • Jelinic P., Stehle J. C. & Shaw P. The testis-specific factor CTCFL cooperates with the protein methyltransferase PRMT7 in H19 imprinting control region methylation. PLoS Biol. 4, e355 (2006). [PubMed]
  • Virtue A., Wang H. & Yang X. F. MicroRNAs and toll-like receptor/interleukin-1 receptor signaling. J Hematol Oncol. 5, 66 (2012). [PMC free article] [PubMed]
  • Yu X., Lin J., Zack D. J. & Qian J. Computational analysis of tissue-specific combinatorial gene regulation: predicting interaction between transcription factors in human tissues. Nucleic Acids Res. 34, 4925–4936 (2006). [PubMed]
  • Lagos-Quintana M. et al. Identification of tissue-specific microRNAs from mouse. Curr Biol. 12, 735–739 (2002). [PubMed]
  • Chen J. F. et al. The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet. 38, 228–233 (2006). [PMC free article] [PubMed]
  • Wienholds E. et al. MicroRNA expression in zebrafish embryonic development. Science. 309, 310–311 (2005). [PubMed]
  • Liang Y., Ridzon D., Wong L. & Chen C. Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 8, 166 (2007). [PMC free article] [PubMed]
  • Lee E. J. et al. Systematic evaluation of microRNA processing patterns in tissues, cell lines, and tumors. RNA. 14, 35–42 (2008). [PubMed]
  • Beuvink I. et al. A novel microarray approach reveals new tissue-specific signatures of known and predicted mammalian microRNAs. Nucleic Acids Res. 35, e52 (2007). [PMC free article] [PubMed]
  • Bostjancic E., Zidar N., Stajer D. & Glavac D. MicroRNAs miR-1, miR-133a, miR-133b and miR-208 are dysregulated in human myocardial infarction. Cardiology. 115, 163–169 (2010). [PubMed]
  • Fu H. et al. Identification of human fetal liver miRNAs by a novel method. FEBS Lett. 579, 3849–3854 (2005). [PubMed]
  • Bargaje R., Hariharan M., Scaria V. & Pillai B. Consensus miRNA expression profiles derived from interplatform normalization of microarray data. RNA. 16, 16–25 (2010). [PubMed]
  • Townley-Tilson W. H., Callis T. E. & Wang D. MicroRNAs 1, 133, and 206: critical factors of skeletal and cardiac muscle development, function, and disease. Int J Biochem Cell Biol. 42, 1252–1255 (2010). [PMC free article] [PubMed]
  • Liu C. G. et al. An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues. Proc Natl Acad Sci U S A. 101, 9740–9744 (2004). [PubMed]
  • Barad O. et al. MicroRNA expression detected by oligonucleotide microarrays: system establishment and expression profiling in human tissues. Genome Res. 14, 2486–2494 (2004). [PubMed]
  • Ramkissoon S. H. et al. Hematopoietic-specific microRNA expression in human cells. Leuk Res. 30, 643–647 (2006). [PubMed]
  • Mestdagh P. et al. The microRNA body map: dissecting microRNA function through integrative genomics. Nucleic Acids Res. 39, e136 (2011). [PMC free article] [PubMed]
  • Baskerville S. & Bartel D. P. Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA. 11, 241–247 (2005). [PubMed]
  • Sempere L. F. et al. Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol. 5, R13 (2004). [PMC free article] [PubMed]
  • Ravasi T. et al. An atlas of combinatorial transcriptional regulation in mouse and man. Cell. 140, 744–752 (2010). [PMC free article] [PubMed]
  • Hsu S. D. et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 39, D163–169 (2011). [PMC free article] [PubMed]
  • Xiao F. et al. miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 37, D105–110 (2009). [PMC free article] [PubMed]
  • Grimson A. et al. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 27, 91–105 (2007). [PMC free article] [PubMed]
  • Liu X., Yu X., Zack D. J., Zhu H. & Qian J. TiGER: a database for tissue-specific gene expression and regulation. BMC Bioinformatics. 9, 271 (2008). [PMC free article] [PubMed]
  • Stark C. et al. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 34, D535–539 (2006). [PMC free article] [PubMed]
  • Shannon P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003). [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group