Functional analysis of microRNA and transcription factor synergistic regulatory network based on identifying regulatory motifs in non-small cell lung cancer
- Kening Li†1,
- Zihui Li†1,
- Ning Zhao†1,
- Yaoqun Xu2,
- Yongjing Liu1,
- Yuanshuai Zhou1,
- Desi Shang1,
- Fujun Qiu1,
- Rui Zhang1,
- Zhiqiang Chang1 and
- Yan Xu1Email author
© Li et al.; licensee BioMed Central Ltd. 2013
Received: 10 June 2013
Accepted: 4 November 2013
Published: 7 November 2013
Lung cancer, especially non-small cell lung cancer, is a leading cause of malignant tumor death worldwide. Understanding the mechanisms employed by the main regulators, such as microRNAs (miRNAs) and transcription factors (TFs), still remains elusive. The patterns of their cooperation and biological functions in the synergistic regulatory network have rarely been studied.
Here, we describe the first miRNA-TF synergistic regulation network in human lung cancer. We identified important regulators (MYC, NFKB1, miR-590, and miR-570) and significant miRNA-TF synergistic regulatory motifs by random simulations. The two most significant motifs were the co-regulation of miRNAs and TFs, and TF-mediated cascade regulation. We also developed an algorithm to uncover the biological functions of the human lung cancer miRNA-TF synergistic regulatory network (regulation of apoptosis, cellular protein metabolic process, and cell cycle), and the specific functions of each miRNA-TF synergistic subnetwork. We found that the miR-17 family exerted important effects in the regulation of non-small cell lung cancer, such as in proliferation and cell cycle regulation by targeting the retinoblastoma protein (RB1) and forming a feed forward loop with the E2F1 TF. We proposed a model for the miR-17 family, E2F1, and RB1 to demonstrate their potential roles in the occurrence and development of non-small cell lung cancer.
This work will provide a framework for constructing miRNA-TF synergistic regulatory networks, function analysis in diseases, and identification of the main regulators and regulatory motifs, which will be useful for understanding the putative regulatory motifs involving miRNAs and TFs, and for predicting new targets for cancer studies.
Lung cancer, predominantly non-small cell lung cancer (NSCLC), is a common cause of malignant tumor death worldwide . Since the end of the 20th century, lung cancer has become the leading cause of malignant tumor death, with morbidity and mortality gradually increasing over the past 50 years. Active and passive tobacco smoking is the best-known risk factor for lung cancer development.
Recent advances in genomics, epigenomics, transcriptomics, and molecular pathology, as well as in the sequencing techniques, have led to the identification of many potential factors as biomarkers, which may provide possibilities for the early detection of lung cancer and personalized therapy . Several genes were identified as predictive biomarkers in NSCLC, such as the somatic mutation and gene copy gain of the epidermal growth factor receptor (EGFR) . L-myc is amplified and expressed in human small cell lung cancer . Although the oncogenicity of lung cancer-related genes has been studied extensively, there is limited knowledge of the process of malignant transformation and the regulatory mechanisms of multistep pathogenesis, especially the regulatory network of lung cancer-related genes, which urgently need to be studied .
MicroRNAs (miRNAs) are small non-coding RNAs (~23 nt long) that regulate gene expression at the post-transcriptional level. MiRNAs are encoded by genomic DNA, transcribed by RNA polymerase II and then incorporated into a RNA-induced silencing complex that binds to the 3′-UTR regions of its target mRNAs to repress translation or enhance degradation . In recent years, important roles for miRNAs were identified in developmental timing, tumorigenesis, cell proliferation, and cell death [6, 7]. MiRNAs function as oncogenes and tumor suppressors, and their regulatory effects in lung cancer development and progression have been demonstrated [8–10].
Hsa-let-7a acts as a protective miRNA that suppresses RAS and other transcriptional factors. Hsa-let-7a expression is generally reduced in NSCLC patients [11, 12]. High expression of hsa-miR-155 was reported to be associated with poor survival in lung cancer patients . Hsa-miR-128b directly regulates epidermal growth factor receptor (EGFR), and loss of heterozygosity of hsa-miR-128b was detected frequently in NSCLC patients . Higher tumor miR-92a-2* levels are associated with decreased survival in patients with small cell lung cancer. MiRNAs can act as biomarkers of human lung cancer, and this may have important clinical applications in prognosis prediction and in predicting the molecular pathogenesis of cancer, as well as in the development of targeted therapies [15–17].
At the transcriptional level, transcription factors (TFs) are the main regulators that control the transcription of their target genes by binding to specific DNA sequences in the promoter regions of the genes. TFs and miRNAs are the two largest families of trans-acting, gene regulatory molecules in multicellular organisms, and they share a common regulatory logic . Most genes in a genome are regulated not by a single factor, but instead by a synergistic network of trans-acting factors. At the network level, motifs comprising miRNAs, TFs, and common target genes were found to be widespread in diverse organisms from bacteria to human, suggesting that these motifs serve as basic building blocks of transcription networks. In our work, we have used the term “motif” to describe a small group that illustrates the regulation patterns of a miRNA, a TF, and their target genes. Common motifs, such as feedforward loops (FFLs) and feedback loops (FBLs), have been found to play crucial roles in gene regulation, such as the miR-17 cluster, the E2F1, and the c-Myc that modulates cellular proliferation in cancer . Several databases of TF-miRNA FFLs involved in tumors have been developed [20, 21]. Moreover, protein-protein interactions data have been included to construct regulatory networks for identifying novel regulatory motifs, such as the four or more node FFLs [22, 23].
Previous studies into the co-regulation between miRNAs and TFs found a variety of significant network motifs that were over-represented in the co-regulation network [24, 25], suggesting that the gene regulation system requires close synergistic regulation by transcriptional and post-transcriptional layers. However, the miRNA-TF synergistic effect may not be limited only to the FFLs. Non-loop forms, such as the cascaded form, which have also helped in understanding the regulatory mechanism, should be considered. Therefore, in this study, we have identified multiple types of motifs, including FFLs, miRNA- or TF-mediated FBLs, and miRNA or TF-mediated cascaded patterns.
Here, we used comprehensive data sources and algorithms to predict regulatory relationships of miRNAs and their targets in an attempt to provide the first miRNA-TF-mediated regulatory network in lung cancer. We also identified synergistic motifs of miRNAs and TFs. Several potential major factors were identified in subnetworks. We have developed an algorithm to predict the biological functions involved in the human lung cancer miRNA-TF regulatory network as well as the specific functions regulated by each synergistic motif. Our results showed that miRNAs of the same family exhibited similar regulatory modes, implying that miRNA family members tended to work together, particularly in regulating TFs. The miR-17 family (in an FFL with the E2F1and the RB1) was found to be an important family in the lung cancer regulation network.
This study provides a framework for constructing a lung cancer-related synergistic regulatory network and for analyzing the biological functions of the network. This approach could be applied easily to study other cancers, and may provide useful information for laboratory experiments and target validation.
MiRNA and TF synergistic regulatory network in lung cancer
We collected and curated 1990 human lung cancer-related genes from several disease-related gene databases and 1823 genes that were aberrantly expressed in NSCLC samples. From them, we selected a total of 1002 genes that met the requirements of lung cancer-related genes to use in this study. The 100 bootstrapping repetitions that we conducted on the microarray data showed that the overlap between gene sets calculated based on the re-sampled data and the original gene set (1002 genes) was quite significant, suggesting that our selected lung cancer-related genes were robust. The ratios of overlap genes to original genes are listed in Additional file 1: Table S1.
Summary of relationships in the lung cancer-related synergistic regulatory network
No. of pairs
No. of miRNAs
No. of genes
No. of TFs
Top 10 miRNAs and TFs with highest out-degree in lung cancer synergistic regulatory network
Top 10 genes, TFs, and miRNAs with highest in-degree in lung cancer synergistic regulatory network
Synergistic motif identification and subnetwork construction
Details of motifs in the lung cancer synergistic regulatory network
No. of motifs
No. of genes
No. of TFs
No. of miRNA
Examples of motifs or prognosis components of motifs
miR-106a &E2F1 &RAD51
miR-181a &TP53 &RUNX3
miR-16 &JUN &LPL
miR-17 &STAT1 &ALDH1A3
miR-21 &ESR1 &CXCL12
Let-7d &ATF1 &GSTP1
miR-200c &E2F3 &ALDH1A3
miR-141 &SOX2 & CXCL12
CTNNB1 & miR-21 & SMAD7
GADD45A& miR-34 & P53
Biological functions of the synergistic regulatory networks and subnetworks
First, mutations in the genes of ten subnetworks were analyzed. We downloaded somatic mutation profiles of 538 lung adenocarcinomas (LUAD) and 178 lung squamous cell carcinomas (LUSC) from The Cancer Genome Atlas level 2 data. Then, we selected genes in each subtype with mutation rates greater than 5% as a mutated gene set. Hypergeometric cumulative distribution was used to test the enrichment significance for genes in 10 kinds of motifs. The results showed that the 10 kinds of motifs were all significantly enriched in the mutated gene set.
Biology process terms regulated by the miRNA-TF synergistic regulatory network
Regulation of apoptosis
Regulation of cellular protein metabolic process
Enzyme linked receptor protein signaling pathway
Regulation of protein modification process
Regulation of phosphorylation
Regulation of phosphate metabolic process
Interface of mitotic cell cycle
G1/S transition of mitotic cell cycle
Regulation of protein phosphorylation
Transmembrane receptor protein tyrosine kinase signaling pathway
Regulation of cell cycle arrest
Regulation of protein kinase activity
Cell cycle checkpoint
DNA metabolic process
Regulation of kinase activity
Regulation of G1/S transition of mitotic cell cycle
S phase of mitotic cell cycle
Mitotic cell cycle checkpoint
Mitotic cell cycle G1/S transition checkpoint
G1/S transition checkpoint
Negative regulation of apoptosis
Negative regulation of programmed cell death
Positive regulation of cellular biosynthetic process
Regulation of protein serine/threonine kinase activity
Negative regulation of signal transduction
Nerve growth factor receptor signaling pathway
7/I, II, V
Regulation of transcription, DNA-dependent
Regulation of RNA biosynthetic process
Positive regulation of signal transduction
Regulation of transcription from RNA polymerase II promoter
Positive regulation of apoptosis
Positive regulation of transcription, DNA-dependent
Pathways regulated by miRNA-TF synergistic regulatory network
Direct p53 effectors
Regulation of Telomerase
Hypoxia and oxygen homeostasis regulation of HIF-1-alpha
Arf6 signaling events
Cell Cycle, Mitotic
Synthesis of DNA
Regulation of DNA replication
Orc1 removal from chromatin
Switching of origins to a post-reflective state
Removal of licensing factors from origins
Signaling events regulated by Ret tyrosine kinase
E-cadherin signaling events
FOXA transcription factor networks
E2F transcription factor network
Neurotrophic factor-mediated Trk receptor signaling
Canonical Wnt signaling pathway
We speculated whether the specific BP terms for each subnetwork were regulated specifically by each motif type. After removing the 36 common BP terms, the remaining terms were ranked by their enrichment frequency in the subnetworks and then categorized within each subnetwork. Details of the results are shown in Additional file 1: Table S5. The functions of the miRNAs in the lung cancer regulatory network were predicted from subnetwork VII because it comprised motifs with miRNAs that simultaneously regulated TFs and genes, while TFs that regulated genes or miRNAs were not included in these motifs. The predicted functions were regulation of fibroblast growth factor signaling pathway, inositol lipid-mediated signaling, response to insulin stimulus, MAPK cascade, receptor signaling pathway, cell migration, DNA replication and metabolism, phosphorylation, enzymatic activity, and meiosis. Similarly, the functions of the TFs were predicted from subnetwork X; they included regulation of protein metabolic process, apoptosis and programmed cell death, gene expression, phosphorylation, and regulation of enzyme activity. Each motif in subnetwork I comprised a FBL and FFL, and not surprisingly, their specific function was DNA replication, which requires precise and complex regulation because of its ubiquity in cells and the multiple enzymes involvement.
Interplay of miRNA and TF in the human lung cancer regulatory network
Of the 252 miRNAs in the regulatory network, 93 (36.9%) were regulated by TFs. Most of these miRNAs had low in-degree, and only 11.8% (11 of the 93) had in-degrees greater than 5. Of the 173 TFs in the network, 27 (15.6%) were regulated by miRNAs, and 37.0% (10 of 27) of them had in-degrees greater than 50. A total of 57 TFs regulated miRNAs and 244 miRNAs regulated TFs. On average, each TF was regulated by 9.33 miRNAs, while each miRNA was regulated by 1.86 TFs. By comparing the intensity and density of the interplay between the miRNAs and TFs in the lung cancer regulatory network, we found that only a small number of the TFs were regulated by miRNAs at a high intensity, while most miRNAs were regulated by TFs at a significantly lower intensity (Additional file 1: Figure S11).
A model of the miR-17 family, RB1, and E2F1 motif in lung cancer proliferation
In subnetwork I, we discovered the predicted interactions between the miR-17 family and E2F1 for the first time (Additional file 1: Table S6). Interestingly, six members of the miR-17 family (miR-17/20a/20b/106a/106b/93) clustered in one group. The miR-17 family and E2F1 formed a FBL, which was a clique. A similar phenomenon was reported for the miR-17-92 cluster (miR-17/20a/18a/19a/19b-1/92-1), which forms a FBL with E2F1, and plays roles in regulating cellular proliferation and apoptosis. The interactions between the miR-17-92 cluster and E2F1 have been verified experimentally [50–54]. The miR-17 family and the miR-17-92 cluster have two shared members, miR-17 and miR-20a, both of which were confirmed to interact with E2F1. For the other four members of the miR-17 family, we performed a sequence alignment to examine how likely they were to interact with E2F1. The conserved sequence of E2F1 among five species (Homo sapiens, Mus musculus, Pan troglodytes, Rattus norvegicus and Bos taurus) was aligned to the mature sequences of the six members of the miR-17 family. All miR-17 family members shared 8–15 bases with the E2F1 conserved sequence (Additional file 1), and their interaction was predicted by at least five algorithms (Additional file 1: Table S6), supporting the high possibility of an interaction between the miR-17 family and E2F1.
The RB1 tumor suppressor negatively regulates the cell cycle and is inactivated in a wide range of human tumors . In subnetwork I, RB1 was targeted by members of the miR-17 family (Additional file 1: Table S6) and E2F1, while the miR-17 family members and E2F1 targeted each other, thereby forming an FFL. By analyzing the miRNA expression profiles of NSCLC patients, we found that five miR-17 family members were significantly overexpressed, the expression of RB1 was significantly down-regulated, and E2F1 expression was not significantly different (Additional file 1: Table S7). MiR-20b was the exception to this because no probe was detected. Next, we examined the mechanism by which the miR-17 family regulates cell cycle and tumor progression in lung cancer using a hypothetical model. The interaction between the pRb proteins and the E2F TF family plays a central role in regulating cell cycle progression by controlling the expression of E2F-dependent cell cycle genes . The overexpressed miR-17 family may directly decrease the translation of RB1, thereby lowering the expression of the RB1 protein. In G0 or early G1 cells, Rb protein, which has been functionally inactivated by transcriptional suppression, releases the transactivation domains of E2F and activates the expression of genes that encode products necessary for S-phase progression . Moreover, E2F1 promotes the transcription of the miR-17 family, which causes overexpression of the miR-17 family members, thereby governing cell cycle and proliferation of lung tumors by targeting RB1 protein.
Here we constructed the first lung cancer-related miRNA-TF synergistic regulatory network of lung cancer. We identified 10 types of motifs and constructed 10 subnetworks. More than half of the putative hub nodes were verified by examining other published works, which indicated the robustness of the network. We developed an algorithm to understand the common and specific functions of these networks. Finally, we proposed a hypothetical model to explain the role of the miR-17 family in regulating cell cycle and tumor progression by targeting the RB1 protein in NSCLC.
In the human miRNA-TF synergistic regulatory network and subnetworks, hub genes and hub miRNAs were identified. Most either were known lung cancer-related factors or were reported to play important roles in lung cancer. The hubs with highest out-degrees in the regulatory network, Myc, TP53, and E2F1, are TFs that play roles in apoptosis, cell proliferation, and lung tumor development. The amplification and overexpression of Myc has been detected in lung cancer of different histologic subtypes . TP53 encodes tumor protein p53, abnormalities of which are frequently found in lung cancers . E2F1 overexpression was reported produce more aggressive tumors with a high proliferation rate during the progression of NSCLC . MiR-17/106a/20a/93/34a were the hubs of many subnetworks, and four of them belong to the miR-17 family. MiR-17 and miR-20a were reported to induce apoptosis in lung cancer cells  and miR-34 s was found to be dramatically down-regulated in NSCLC .
In this work, we proposed a model to predict the regulatory role of the miR-17 family in the cell cycle via RB1 and E2F1. In the model, five out of six miR-17 family members were significantly overexpressed in NSCLC cells where they enhanced the repression of the RB1 gene, which is responsible for the G1 checkpoint and blockage of S-phase entry and cell growth. Hesan et al.  confirmed the up-regulation of four members of the miR-17 family in colorectal carcinoma tissues and showed that they promote cell proliferation and tumor growth by targeting the RND3 tumor suppressor gene. A similar group, the miR-17-92 cluster with two members that were common with the miR-17 family, had diverse functions in the regulation of cellular differentiation, proliferation, and apoptosis. The two common members, miR-17 and miR-20a, were shown to temper an E2F1-induced G1 checkpoint to regulate cell cycle progression . Furthermore, the E2F and the miR-17-92 cluster could form FBLs , and in the cancer regulation network, FBLs involving miR-17-92, E2F and MYC have been reported . We checked the interactive relations of the miR-17 family with E2F1 and RB1 by sequence alignment and found a strong possibility of their interactions. Moreover, many regulatory relationships support our predictive model of the miR-17 family, E2F1, and RB1 motif, which demonstrates the effectiveness of our regulatory network.
After identifying the miRNA-TF synergistic motifs, we calculated their significance and Z-values, and ranked the motifs according to their Z scores. The first-ranked Motif VII was more significant than the second-ranked Motif VI, possibly because of the availability of abundant miRNA regulation data but insufficient TF regulation data. One reason that Motif VI was found to be the most significant regulatory motif in the network may be that genes are first regulated by TFs at the transcription level and then by miRNAs at the post-transcription level; thus, genes are significantly regulated by TFs and miRNAs separately at different times and in different locations in the cell. By comparing two Motifs, II and III, and by examining the regulatory directions between the miRNAs and TFs, we found that miRNAs tended to be significantly regulated by TFs rather than regulate TFs. This observation is despite the fact that data on the - targets of TFs are limited, while much more data on the targets of miRNAs are available. Therefore, we inferred that TFs play a dominant role in FFL regulation. This assumption is supported by the results of another study, which found that TFs held dominant positions in the global regulatory system (i.e. at the transcriptional level) compared with the miRNAs at the downstream positions (i.e. at the post-transcriptional level) . Between the linear motifs (VIII, IX), Motif IX was more significant than Motif VIII, which indicated that genes tended to be regulated directly by TFs, while miRNAs tended to regulate TFs while they were being formed rather than act as a mediated regulator between TFs and their target genes.
The expression level analysis of genes and miRNAs may help in understanding the regulatory mechanisms; therefore, the differentially expressed genes in our networks further investigated. In general, we found that up-regulated miRNAs down-regulated their target genes by degrading them at the transcript level or by repressing protein production at the translational level. We also observed the reverse, in which down-regulated miRNAs led to up-regulated target genes. However, down-regulated miRNAs that down-regulated their target genes and up-regulated miRNAs that up-regulated their target genes were also observed in our network. This may be because the expression levels of genes or miRNAs are determined by multiple factors, including environment, heredity, copy number variations, and epigenetics. Thus, miRNAs and TFs may influence expression to a great degree rather than being the decisive factors.
In this study, we analyzed the regulation of genes by miRNAs and TFs, but did not consider gene-to-gene relationships. Cui et al. [60, 61] studied the relationship among oncogenes in the context of activity/inhibitory motifs and compared the number of mutant genes and miRNA target genes in each type of motif. In the future, we will examine the activation, inhibition, and physical interactions among the genes in regulatory networks, and discuss the regularity role of miRNAs, TFs, and motifs. Our future studies will contribute to uncovering the principles of miRNA regulation in signal transduction networks.
In summary, our established miRNA and TF synergistic regulatory network in NSCLC has provided clues about the regulatory mechanisms of lung cancer and information that will help identify the core regulators. Nearly half of the hub regulators, as well as the proposed regulatory motifs, were confirmed by literature searches, which indicated the effectiveness and rationality of the network construction. The most significant motifs were of the co-regulation and TF-mediated cascade regulation types. While cooperating with miRNAs, TFs tended to play a dominant role in FFL regulations. We also developed an algorithm to analyze the functions of the human lung cancer miRNA-TF regulatory network and subnetworks. According to the full regulation subnetwork and expression analysis, we proposed a predictive model of the miR-17 family, E2F1 and RB1 in the regulation of cell cycle and cellular proliferation. Our study will provide valuable information for lung cancer investigators to identify critical elements and regulatory motifs for a better understanding of the regulatory mechanisms or for designing future experimental studies of lung cancer.
Lung cancer-related genes
Target prediction of miRNAs and TFs
For miRNA target prediction, we combined 10 popular databases or algorithms, namely, miRanda, TargetScan, PicTar5, PITA, DIANA-microT, mirSVR, RNA22, RNAhybrid, MirTarget2, and TargetMiner (Additional file 1: Table S8). To decrease the number of false-positive results, only miRNA–mRNA interactions predicted by at least three of the algorithms were accepted as positive. TFs were treated as genes when predicting miRNA–TF regulations.
TFs and their experimentally proven targets were retrieved from four databases: ORegAnno , Pazar , Transfac , and Tred . To ensure that the results were complete, we used a union set of all the retrieved data.
MiRNA precursor sequences were obtained from the miRBase database. We selected the 2-kb upstream regions of the pre-miRNAs as their putative promoter regions, and then searched the sequences for TF binding sites using the UCSC genome browser (Z score =2.33) . The selected sites were required to be conserved among human, rat, and mouse. We also incorporated experimentally supported TF–miRNA regulatory relations that were curated manually from large numbers of published papers  and from TransmiR  (Figure 5C).
Significant lung cancer-related miRNAs and TFs
where, M is the total number of all human genes in the human genome, k is the number of lung cancer-related genes, N is the number of target genes for a certain miRNA (or TF), and x is the number of target genes of a certain miRNA (or TF) that overlapped with lung cancer-related genes. To avoid the effect from false-positives in the miRNA target data, we conducted randomization tests to ensure the biological significance of the identified lung cancer-related miRNAs and TFs. Instead of using 1002 lung cancer-related genes to enrich each of the miRNAs and TFs, we randomly selected 1002 genes from all the coding genes in the human genome for enrichment, and this procedure was repeated 1000 times.
Motif identification and statistics test
where, N high is the number of random times that an acquired motif was more than or equal to the real network, and N random is 10000.
where, N real and N mean are the number of motifs observed in the real synergistic regulatory network and their mean occurrence in the random networks, respectively. SD is the standard deviation of the number of motifs in the random networks.
Motif-specific subnetwork and hub definition
All network motifs of the same type were merged to construct a motif-specific subnetwork (for example, all FFLs were merged to form the FFL subnetwork) and then the subnetworks were visualized by the Cytoscape software . When network sizes were greatly different from one another, we defined hub nodes as the top 5% highest-degree nodes of the miRNAs, TFs, and genes in both the subnetwork and whole network.
Expression level analysis of miRNAs, TFs, and genes
MiRNA and gene expression profiles were downloaded from the NCBI Gene Expression Omnibus (GEO), and differentially expressed genes were calculated using the significance analysis of microarrays (SAM) software  with the FDR set to less than 1%. The union set of differentially expressed genes in the GSE2088 and GSE11969 datasets was used to determine the lung cancer-related genes, while the intersection was shown in different colors in the motif-specific subnetwork. The miRNA expression profiles from the GSE27705 dataset, which included 20 NSCLC samples and 10 normal samples, were used to identify up- or down-regulated expressed miRNAs.
BP and pathway analysis of the genes in the human lung cancer synergistic regulatory network and synergistic subnetworks
For the function analysis of each subnetwork, we developed an algorithm to obtain the most representative functions among the numerous GO BP terms and KEGG pathways after BP and pathway enrichment analysis. For example, for the BP function analysis of subnetworki (Motifi) the algorithm: 1) found the target gene set shared by each miRNA-TF pair, which contained no less than three genes; 2) all target gene sets were enriched with BP terms by the hypergeometric test with P-value adjustment by FDR (the cutoff of P-value was 0.005); and 3) the frequency that each BP term was enriched was counted and represented as Frequent i (i∈[1, 10]), and all enriched terms were ranked in descending order according to their frequency value, so that the rank number of each BP term, represented as Motif i -rank, could be obtained.
where, stagenumber is the number of stages in motifi. The terms with the highest top 5% score were regarded as common terms, and were used to represent the functions of the whole regulatory network. We then removed the common terms from each subnetwork and ranked the remaining terms to identify the specific functions of each subnetwork. Pathway function analysis was conducted using the same procedure.
Hierarchical clustering and sequence alignment
In our analysis of the interplay between miRNAs and TFs, two-way clustering of the regulators was performed using Cluster 3.0 software  and the resultant heatmap was viewed using TreeView. MiRNA family information was obtained from miRBase . Multiple sequences were aligned using ClustalW2 software in the analysis tools framework at EMBL-EBI .
Kening Li, Zihui Li and Ning Zhao considered as co-first authors.
This work was supported in part by the National Natural Science Foundation of China (Grant No. 81372492), the Natural Science Foundation of Heilongjiang Province (Grant No. D201116), the Postdoctoral science-research developmental foundation of Heilongjiang Province (Grant No. LBH-Q11044), and the Innovation Research Fund for Graduate Students of Heilongjiang Province (YJSCX2012-211HLJ). We thank Shanzhen Zhang for the programming used in identifying motifs. We thank Yingying Wang for assistance with collecting the regulatory data.
- Siegel R, Naishadham D, Jemal A: Cancer statistics, 2012. CA: a cancer journal for clinicians. 2012, 62: 10-29. 10.3322/caac.20138.View ArticleGoogle Scholar
- Chen HY, Yu SL, Li KC, Yang PC: Biomarkers and transcriptome profiling of lung cancer. Respirology. 2012, 17: 620-626. 10.1111/j.1440-1843.2012.02154.x.PubMedView ArticleGoogle Scholar
- Dahabreh IJ, Linardou H, Siannis F, Kosmidis P, Bafaloukos D, Murray S: Somatic EGFR mutation and gene copy gain as predictive biomarkers for response to tyrosine kinase inhibitors in non-small cell lung cancer. Clin Cancer Res. 2010, 16: 291-303. 10.1158/1078-0432.CCR-09-1660.PubMedView ArticleGoogle Scholar
- Nau MM, Brooks BJ, Battey J, Sausville E, Gazdar AF, Kirsch IR, McBride OW, Bertness V, Hollis GF, Minna JD: L-myc, a new myc-related gene amplified and expressed in human small cell lung cancer. Nature. 1985, 318: 69-73. 10.1038/318069a0.PubMedView ArticleGoogle Scholar
- Ye H, Liu X, Lv M, Wu Y, Kuang S, Gong J, Yuan P, Zhong Z, Li Q, Jia H, et al: MicroRNA and transcription factor co-regulatory network analysis reveals miR-19 inhibits CYLD in T-cell acute lymphoblastic leukemia. Nucleic Acids Res. 2012, 40: 5201-5214. 10.1093/nar/gks175.PubMedPubMed CentralView ArticleGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355. 10.1038/nature02871.PubMedView ArticleGoogle Scholar
- Schickel R, Boyerinas B, Park SM, Peter ME: MicroRNAs: key players in the immune system, differentiation, tumorigenesis and cell death. Oncogene. 2008, 27: 5959-5974. 10.1038/onc.2008.274.PubMedView ArticleGoogle Scholar
- Zhang B, Pan X, Cobb GP, Anderson TA: microRNAs as oncogenes and tumor suppressors. Dev Biol. 2007, 302: 1-12. 10.1016/j.ydbio.2006.08.028.PubMedView ArticleGoogle Scholar
- Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.PubMedView ArticleGoogle Scholar
- Khoshnaw SM, Green AR, Powe DG, Ellis IO: MicroRNA involvement in the pathogenesis and management of breast cancer. J Clin Pathol. 2009, 62: 422-428. 10.1136/jcp.2008.060681.PubMedView ArticleGoogle Scholar
- Grosshans H, Johnson T, Reinert KL, Gerstein M, Slack FJ: The temporal patterning microRNA let-7 regulates several transcription factors at the larval to adult transition in C. elegans. Dev Cell. 2005, 8: 321-330. 10.1016/j.devcel.2004.12.019.PubMedView ArticleGoogle Scholar
- Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, et al: Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004, 64: 3753-3756. 10.1158/0008-5472.CAN-04-0637.PubMedView ArticleGoogle Scholar
- Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, et al: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer cell. 2006, 9: 189-198. 10.1016/j.ccr.2006.01.025.PubMedView ArticleGoogle Scholar
- Weiss GJ, Bemis LT, Nakajima E, Sugita M, Birks DK, Robinson WA, Varella-Garcia M, Bunn PA, Haney J, Helfrich BA, et al: EGFR regulation by microRNA in lung cancer: correlation with clinical response and survival to gefitinib and EGFR expression in cell lines. Ann Oncol. 2008, 19: 1053-1059. 10.1093/annonc/mdn006.PubMedView ArticleGoogle Scholar
- Bartels CL, Tsongalis GJ: MicroRNAs: novel biomarkers for human cancer. Clin Chem. 2009, 55: 623-631. 10.1373/clinchem.2008.112805.PubMedView ArticleGoogle Scholar
- Raponi M, Dossey L, Jatkoe T, Wu X, Chen G, Fan H, Beer DG: MicroRNA classifiers for predicting prognosis of squamous cell lung cancer. Cancer Res. 2009, 69: 5776-5783. 10.1158/0008-5472.CAN-09-0587.PubMedView ArticleGoogle Scholar
- Yu SL, Chen HY, Chang GC, Chen CY, Chen HW, Singh S, Cheng CL, Yu CJ, Lee YC, Chen HS, et al: MicroRNA signature predicts survival and relapse in lung cancer. Cancer cell. 2008, 13: 48-57. 10.1016/j.ccr.2007.12.008.PubMedView ArticleGoogle Scholar
- Hobert O: Gene regulation by transcription factors and microRNAs. Science. 2008, 319: 1785-1786. 10.1126/science.1151651.PubMedView ArticleGoogle Scholar
- O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435: 839-843. 10.1038/nature03677.PubMedView ArticleGoogle Scholar
- El Baroudi M, Cora D, Bosia C, Osella M, Caselle M: A curated database of miRNA mediated feed-forward loops involving MYC as master regulator. PloS one. 2011, 6: e14742-10.1371/journal.pone.0014742.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang J, Lu M, Qiu C, Cui Q: TransmiR: a transcription factor-microRNA regulation database. Nucleic acids Res. 2010, 38: D119-122. 10.1093/nar/gkp803.PubMedPubMed CentralView ArticleGoogle Scholar
- Sun J, Gong X, Purow B, Zhao Z: Uncovering MicroRNA and Transcription Factor Mediated Regulatory Networks in Glioblastoma. PLoS Comput Biol. 2012, 8: e1002488-10.1371/journal.pcbi.1002488.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin CC, Chen YJ, Chen CY, Oyang YJ, Juan HF, Huang HC: Crosstalk between transcription factors and microRNAs in human protein interaction network. BMC Syst Biol. 2012, 6: 18-10.1186/1752-0509-6-18.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen CY, Chen ST, Fuh CS, Juan HF, Huang HC: Coregulation of transcription factors and microRNAs in human transcriptional regulatory network. BMC bioinformatics. 2011, 12 (Suppl 1): S41-10.1186/1471-2105-12-S1-S41.PubMedPubMed CentralView ArticleGoogle Scholar
- Qiu C, Wang J, Yao P, Wang E, Cui Q: microRNA evolution in a human transcription factor and microRNA regulatory network. BMC Syst Biol. 2010, 4: 90-10.1186/1752-0509-4-90.PubMedPubMed CentralView ArticleGoogle Scholar
- Zajac-Kaye M: Myc oncogene: a key component in cell cycle regulation and its implication for lung cancer. Lung Cancer. 2001, 34 (Suppl 2): S43-46.PubMedView ArticleGoogle Scholar
- Toyooka S, Tsuda T, Gazdar AF: The TP53 gene, tobacco exposure, and lung cancer. Human mutation. 2003, 21: 229-239. 10.1002/humu.10177.PubMedView ArticleGoogle Scholar
- Duan HY, Cao JX, Qi JJ, Wu GS, Li SY, An GS, Jia HT, Cai WW, Ni JH: E2F1 enhances 8-chloro-adenosine-induced G2/M arrest and apoptosis in A549 and H1299 lung cancer cells. Biochemistry Biokhimiia. 2012, 77: 261-269. 10.1134/S0006297912030042.PubMedView ArticleGoogle Scholar
- Rauch TA, Wang Z, Wu X, Kernstine KH, Riggs AD, Pfeifer GP: DNA methylation biomarkers for lung cancer. Tumour Biol. 2012, 33: 287-296. 10.1007/s13277-011-0282-2.PubMedView ArticleGoogle Scholar
- Hsu TI, Wang MC, Chen SY, Yeh YM, Su WC, Chang WC, Hung JJ: Sp1 expression regulates lung tumor progression. Oncogene. 2012, 31: 3973-3988. 10.1038/onc.2011.568.PubMedPubMed CentralView ArticleGoogle Scholar
- Bankovic J, Stojsic J, Jovanovic D, Andjelkovic T, Milinkovic V, Ruzdijic S, Tanic N: Identification of genes associated with non-small-cell lung cancer promotion and progression. Lung Cancer. 2010, 67: 151-159. 10.1016/j.lungcan.2009.04.010.PubMedView ArticleGoogle Scholar
- Munksgaard Persson M, Johansson ME, Monsef N, Planck M, Beckman S, Seckl MJ, Ronnstrand L, Pahlman S, Pettersson HM: HIF-2alpha expression is suppressed in SCLC cells, which survive in moderate and severe hypoxia when HIF-1alpha is repressed. Am J Pathol. 2012, 180: 494-504. 10.1016/j.ajpath.2011.10.014.PubMedView ArticleGoogle Scholar
- Li J, Yu B, Song L, Eschrich S, Haura EB: Effects of IFN-gamma and Stat1 on gene expression, growth, and survival in non-small cell lung cancer cells. J Interferon Cytokine Res. 2007, 27: 209-220. 10.1089/jir.2006.0111.PubMedView ArticleGoogle Scholar
- Navarro A, Marrades RM, Vinolas N, Quera A, Agusti C, Huerta A, Ramirez J, Torres A, Monzo M: MicroRNAs expressed during lung cancer development are expressed in human pseudoglandular lung embryogenesis. Oncology. 2009, 76: 162-169. 10.1159/000201569.PubMedView ArticleGoogle Scholar
- Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer research. 2005, 65: 9628-9632. 10.1158/0008-5472.CAN-05-2352.PubMedView ArticleGoogle Scholar
- Terry J, De Luca A, Leung S, Peacock G, Wang Y, Elliot WM, Huntsman D: Immunohistochemical expression of neurotrophic tyrosine kinase receptors 1 and 2 in lung carcinoma: potential discriminators between squamous and nonsquamous subtypes. Arch Pathol Lab Med. 2011, 135: 433-439.PubMedGoogle Scholar
- Nakagawa M, Uramoto H, Oka S, Chikaishi Y, Iwanami T, Shimokawa H, So T, Hanagiri T, Tanaka F: Clinical significance of IGF1R expression in non-small-cell lung cancer. Clin Lung Cancer. 2012, 13: 136-142. 10.1016/j.cllc.2011.10.006.PubMedView ArticleGoogle Scholar
- Feng J, Zhang X, Zhu H, Wang X, Ni S, Huang J: High expression of FoxP1 is associated with improved survival in patients with non-small cell lung cancer. Am J Clin Pathol. 2012, 138: 230-235. 10.1309/AJCPDHQFNYJZ01YG.PubMedView ArticleGoogle Scholar
- Toyokawa G, Cho HS, Masuda K, Yamane Y, Yoshimatsu M, Hayami S, Takawa M, Iwai Y, Daigo Y, Tsuchiya E, et al: Histone lysine methyltransferase Wolf-Hirschhorn syndrome candidate 1 is involved in human carcinogenesis through regulation of the Wnt pathway. Neoplasia. 2011, 13: 887-898.PubMedPubMed CentralView ArticleGoogle Scholar
- Bandi N, Zbinden S, Gugger M, Arnold M, Kocher V, Hasan L, Kappeler A, Brunner T, Vassella E: miR-15a and miR-16 are implicated in cell cycle regulation in a Rb-dependent manner and are frequently deleted or down-regulated in non-small cell lung cancer. Cancer Res. 2009, 69: 5553-5559. 10.1158/0008-5472.CAN-08-4277.PubMedView ArticleGoogle Scholar
- Gibbons DL, Lin W, Creighton CJ, Rizvi ZH, Gregory PA, Goodall GJ, Thilaganathan N, Du L, Zhang Y, Pertsemlidis A, Kurie JM: Contextual extracellular cues promote tumor cell EMT and metastasis by regulating miR-200 family expression. Genes & development. 2009, 23: 2140-2151. 10.1101/gad.1820209.View ArticleGoogle Scholar
- Cooper CS, Nicholson AG, Foster C, Dodson A, Edwards S, Fletcher A, Roe T, Clark J, Joshi A, Norman A, et al: Nuclear overexpression of the E2F3 transcription factor in human lung cancer. Lung Cancer. 2006, 54: 155-162. 10.1016/j.lungcan.2006.07.005.PubMedView ArticleGoogle Scholar
- Suga Y, Miyajima K, Oikawa T, Maeda J, Usuda J, Kajiwara N, Ohira T, Uchida O, Tsuboi M, Hirano T, et al: Quantitative p16 and ESR1 methylation in the peripheral blood of patients with non-small cell lung cancer. Oncology Rep. 2008, 20: 1137-1142.Google Scholar
- Li X, Yang XX, Hu NY, Sun JZ, Li FX, Li M: A risk-associated single nucleotide polymorphism of SMAD7 is common to colorectal, gastric, and lung cancers in a Han Chinese population. Mol Biol Rep. 2011, 38: 5093-5097. 10.1007/s11033-010-0656-3.PubMedView ArticleGoogle Scholar
- Yamaguchi E, Nakayama T, Nanashima A, Matsumoto K, Yasutake T, Sekine I, Nagayasu T: Ets-1 proto-oncogene as a potential predictor for poor prognosis of lung adenocarcinoma. Tohoku J Exp Med. 2007, 213: 41-50. 10.1620/tjem.213.41.PubMedView ArticleGoogle Scholar
- Bai J, Hu S: Transcriptome network analysis reveals potential candidate genes for squamous lung cancer. Int J Mol Med. 2012, 29: 95-101.PubMedGoogle Scholar
- Yan Z, Shah PK, Amin SB, Samur MK, Huang N, Wang X, Misra V, Ji H, Gabuzda D, Li C: Integrative analysis of gene and miRNA expression profiles with transcription factor-miRNA feed-forward loops identifies regulators in human cancers. Nucleic acids Res. 2012, 40: e135-10.1093/nar/gks395.PubMedPubMed CentralView ArticleGoogle Scholar
- Chiba T, Uchi H, Yasukawa F, Furue M: Role of the arylhydrocarbon receptor in lung disease. Int Arch Allergy Immunol. 2011, 155 (Suppl 1): 129-134.PubMedView ArticleGoogle Scholar
- Singer S, Malz M, Herpel E, Warth A, Bissinger M, Keith M, Muley T, Meister M, Hoffmann H, Penzel R, et al: Coordinated expression of stathmin family members by far upstream sequence element-binding protein-1 increases motility in non-small cell lung cancer. Cancer Res. 2009, 69: 2234-2243. 10.1158/0008-5472.CAN-08-3338.PubMedView ArticleGoogle Scholar
- Pickering MT, Stadler BM, Kowalik TF: miR-17 and miR-20a temper an E2F1-induced G1 checkpoint to regulate cell cycle progression. Oncogene. 2009, 28: 140-145. 10.1038/onc.2008.372.PubMedPubMed CentralView ArticleGoogle Scholar
- Sylvestre Y, De Guire V, Querido E, Mukhopadhyay UK, Bourdeau V, Major F, Ferbeyre G, Chartrand P: An E2F/miR-20a autoregulatory feedback loop. J Biol Chem. 2007, 282: 2135-2143.PubMedView ArticleGoogle Scholar
- Aguda BD, Kim Y, Piper-Hunter MG, Friedman A, Marsh CB: MicroRNA regulation of a cancer network: consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proc Natl Acad Sci USA. 2008, 105: 19678-19683. 10.1073/pnas.0811166106.PubMedPubMed CentralView ArticleGoogle Scholar
- Conkrite K, Sundby M, Mukai S, Thomson JM, Mu D, Hammond SM, MacPherson D: miR-17 92 cooperates with RB pathway mutations to promote retinoblastoma. Genes & development. 2011, 25: 1734-1745. 10.1101/gad.17027411.View ArticleGoogle Scholar
- Olive V, Jiang I, He L: mir-17-92, a cluster of miRNAs in the midst of the cancer network. Int J Biochem Cell Biol. 2010, 42: 1348-1354. 10.1016/j.biocel.2010.03.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Giacinti C, Giordano A: RB and cell cycle progression. Oncogene. 2006, 25: 5220-5227. 10.1038/sj.onc.1209615.PubMedView ArticleGoogle Scholar
- Macaluso M, Montanari M, Giordano A: Rb family proteins as modulators of gene expression and new aspects regarding the interaction with chromatin remodeling enzymes. Oncogene. 2006, 25: 5263-5267. 10.1038/sj.onc.1209680.PubMedView ArticleGoogle Scholar
- Huang CL, Liu D, Nakano J, Yokomise H, Ueno M, Kadota K, Wada H: E2F1 overexpression correlates with thymidylate synthase and survivin gene expressions and tumor proliferation in non small-cell lung cancer. Clin Cancer Res. 2007, 13: 6938-6946. 10.1158/1078-0432.CCR-07-1539.PubMedView ArticleGoogle Scholar
- Bommer GT, Gerin I, Feng Y, Kaczorowski AJ, Kuick R, Love RE, Zhai Y, Giordano TJ, Qin ZS, Moore BB, et al: p53-mediated activation of miRNA34 candidate tumor-suppressor genes. Curr Biol. 2007, 17: 1298-1307. 10.1016/j.cub.2007.06.068.PubMedView ArticleGoogle Scholar
- Luo H, Zou J, Dong Z, Zeng Q, Wu D, Liu L: Up-regulated miR-17 promotes cell proliferation, tumour growth and cell cycle progression by targeting the RND3 tumour suppressor gene in colorectal carcinoma. Biochem J. 2012, 442: 311-321. 10.1042/BJ20111517.PubMedView ArticleGoogle Scholar
- Cui Q, Ma Y, Jaramillo M, Bari H, Awan A, Yang S, Zhang S, Liu L, Lu M, O'Connor-McCourt M, et al: A map of human cancer signaling. Mol Syst Biol. 2007, 3: 152-PubMedPubMed CentralView ArticleGoogle Scholar
- Cui Q, Yu Z, Purisima EO, Wang E: Principles of microRNA regulation of a human cellular signaling network. Mol Syst Biol. 2006, 2: 46-PubMedPubMed CentralView ArticleGoogle Scholar
- Yu W, Clyne M, Khoury MJ, Gwinn M: Phenopedia and Genopedia: disease-centered and gene-centered views of the evolving knowledge of human genetic associations. Bioinformatics. 2010, 26: 145-146. 10.1093/bioinformatics/btp618. 64PubMedPubMed CentralView ArticleGoogle Scholar
- Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, Jia M, Shepherd R, Leung K, Menzies A, et al: COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic acids Res. 2011, 39: D945-950. 10.1093/nar/gkq929.PubMedPubMed CentralView ArticleGoogle Scholar
- Becker KG, Barnes KC, Bright TJ, Wang SA: The genetic association database. Nat Genet. 2004, 36: 431-432. 10.1038/ng0504-431.PubMedView ArticleGoogle Scholar
- Hamosh A, Scott AF, Amberger JS, Bocchini CA, McKusick VA: Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic acids Res. 2005, 33: D514-517.PubMedPubMed CentralView ArticleGoogle Scholar
- Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al: NCBI GEO: archive for functional genomics data sets–update. Nucleic acids Res. 2013, 41: D991-995. 10.1093/nar/gks1193.PubMedPubMed CentralView ArticleGoogle Scholar
- Li J, Lenferink AE, Deng Y, Collins C, Cui Q, Purisima EO, O'Connor-McCourt MD, Wang E: Identification of high-quality cancer prognostic markers and metastasis network modules. Nat Commun. 2010, 1: 34-PubMedGoogle Scholar
- Montgomery SB, Griffith OL, Sleumer MC, Bergman CM, Bilenky M, Pleasance ED, Prychyna Y, Zhang X, Jones SJ: ORegAnno: an open access database and curation system for literature-derived promoters, transcription factor binding sites and regulatory variation. Bioinformatics. 2006, 22: 637-640. 10.1093/bioinformatics/btk027.PubMedView ArticleGoogle Scholar
- Portales-Casamar E, Kirov S, Lim J, Lithwick S, Swanson MI, Ticoll A, Snoddy J, Wasserman WW: PAZAR: a framework for collection and dissemination of cis-regulatory sequence annotation. Genome Biol. 2007, 8: R207-10.1186/gb-2007-8-10-r207.PubMedPubMed CentralView ArticleGoogle Scholar
- Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, et al: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic acids Res. 2003, 31: 374-378. 10.1093/nar/gkg108.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao F, Xuan Z, Liu L, Zhang MQ: TRED: a Transcriptional Regulatory Element Database and a platform for in silico gene regulation studies. Nucleic acids Res. 2005, 33: D103-107. 10.1093/nar/gni105.PubMedPubMed CentralView ArticleGoogle Scholar
- Dreszer TR, Karolchik D, Zweig AS, Hinrichs AS, Raney BJ, Kuhn RM, Meyer LR, Wong M, Sloan CA, Rosenbloom KR, et al: The UCSC Genome Browser database: extensions and updates 2011. Nucleic acids Res. 2012, 40: D918-923. 10.1093/nar/gkr1055.PubMedPubMed CentralView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMedPubMed CentralView ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.PubMedPubMed CentralView ArticleGoogle Scholar
- de Hoon MJ, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20: 1453-1454. 10.1093/bioinformatics/bth078.PubMedView ArticleGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic acids Res. 2011, 39: D152-157. 10.1093/nar/gkq1027.PubMedPubMed CentralView ArticleGoogle Scholar
- Goujon M, McWilliam H, Li W, Valentin F, Squizzato S, Paern J, Lopez R: A new bioinformatics analysis tools framework at EMBL-EBI. Nucleic acids research. 2010, 38: W695-699. 10.1093/nar/gkq313.PubMedPubMed CentralView ArticleGoogle Scholar
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.