Computational Model for Predicting the Relationship Between Micro-RNAs and Their Target Messenger RNAs in Breast and Colon Cancers

Abstract

Motivation:

Uncovering the relationship between micro-RNAs (miRNAs) and their target messenger RNAs (mRNAs) can provide critical information regarding the mechanisms underlying certain types of cancers. In this context, we have proposed a computational method, referred to as prediction analysis by optimization method (PAOM), to predict miRNA-mRNA relations using data from normal and cancer tissues, and then applying the relevant algorithms to colon and breast cancers. Specifically, we used 26 miRNAs and 26 mRNAs with 676 (= 26 × 26) relationships to be recovered as unknown parameters.

Results:

Optimization methods were used to detect 61 relationships in breast cancer and 32 relationships in colon cancer. Using sequence filtering, we detected 18 relationships in breast cancer and 15 relationships in colon cancer. Among the 18 relationships, CD24 is the target gene of let-7a and miR-98, and E2F1 is the target gene of miR-20. In addition, the frequencies of the target genes of miR-223, miR-23a, and miR-20 were significant in breast cancer, and the frequencies of the target genes of miR-17, miR-124, and miR-30a were found to be significant in colon cancer.

Availability:

The numerical code is available from the authors on request.

Keywords

miRNA- mRNA relationship computational method cancer data

Introduction

Micro-RNAs (miRNAs) are a large class of small, non-coding RNAs consisting of 22 nucleotides that are expressed from longer and endogenous hairpin-shaped transcripts generally referred to as pre-miRNAs.¹ The miRNAs regulate protein-coding gene expression post-transcriptionally via the translational repression or transcript degradation of their target messenger RNAs (mRNAs),^1,2 thereby indicating that miRNAs perform crucial roles in a variety of biological functions. The results of recent studies have demonstrated that miRNAs are deregulated in cancers.^2,3 Several computational methods have been proposed to determine how miRNAs pair with their target mRNAs,^3,4 in an effort to unravel the roles of miRNAs in the deregulated expression of their target mRNAs during cancer development^5,6 and some relationships were identified between miRNAs and their target mRNAs.^5,7–11 However, the development of a computational method for the identification of such relationships in cancer remains a difficult issue. Thus far, 2 computational methods have been developed: either the identification of miRNAs that are conserved in different species or stem loop prosecutors¹¹ or the identification of the relationship between miRNAs and their target mRNAs via the use of sequence homologues. Previously, we proposed a numerical optimization method for multi (miRNAs) -to -multi (mRNAs), which was used for the identification of 16 miRNAs-mRNA relations in colon cancer microarray profiles.⁷ The proposed method was used to identify 207 relationships successfully, out of 484. Some relationships detected in that study were verified through previous experimental evidences. For example, the relationship between miR-17 and its target E2F1 was identified in a previous colon cancer study.¹² Like previous studies, however, the difference between the expression profiles of normal and cancer tissues, which we believe to be a critical factor in determining the relationship between miRNAs and their target mRNAs, was ignored.

Here, we propose a novel method, referred to as prediction analysis by optimization method (PAOM), which is composed of a mathematical model and computational method designed to predict miRNA-mRNA relations in the context of cancer development. In this model, the predicted relations are filtered using sequence analysis resources. For mathematical modeling, we employed linear system equations to obtain the inhibiting parameters. The role of the computational method is to optimize the relations and to allow for comparisons of the differences in the parameters between normal and cancer genes. We employed 2 optimization methods—the Broyden-Fletcher-Goldfarb-Shannon (BFGS) and the Powell method,¹³ both of which are well known for the optimization of a multi-dimensional matrix problem. For filtering sequence analysis, we used PicTar, based on the scanning multiple alignment of 3′ UTR (untranslated region) sequences and a search set of miRNA sequences, and scored the overlapping position.¹¹ MiRanda is a miRNA target prediction algorithm that optimizes sequence complementarity using position-specific rules and relies on strict interspecies conservation requirements.¹⁴ In this study, we considered breast and colon cancer data with 676 relationships between 26 miRNAs and 26 mRNAs.

Methods

The mathematical formulation and computational method of the PAOM are described in this section.

Mathematical formulation

It is generally accepted that miRNAs regulate gene expression via either the transcript cleavage or translational repression of their specific target mRNAs, whereas 1 mRNA expression is regulated by multiple miRNAs. For this mechanism, we have a linear equation model, in which 1 mRNA is affected by several miRNAs $(x_{1}, x_{2}, \dots, x_{m})$ as follows:

y_{i} = a_{i 1} x_{1} + a_{i 2} x_{2} + \dots + a_{i m} x_{m}, i = 1, \dots, n

(1)

where $a_{i j}$ represents the influence of the jth miRNA on the ith mRNA, $x_{j}$ represents the expression level of the jth miRNA, and $y_{i}$ represents the expression level of the ith mRNA. To observe the relation between m miRNAs and n mRNAs simultaneously, we rewrite the system of linear equations in matrix form.

[\begin{matrix} y_{1} \\ y_{2} \\ ⋮ \\ y_{n} \end{matrix}] = [\begin{matrix} a_{11} & a_{12} & \dots & a_{1 m} \\ a_{21} & a_{22} & \dots & a_{2 m} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ a_{n 1} & a_{n 2} & \dots & a_{n m} \end{matrix}] [\begin{matrix} x_{1} \\ x_{2} \\ ⋮ \\ x_{m} \end{matrix}],

(2)

in which the measurements of the expression of m different miRNAs are denoted by $(x_{1}, \dots, x_{m})$ and of the expression of n different mRNAs are denoted by $(y_{1}, \dots, y_{n})$ . From K times experiments,

Y_{n \times K} = A_{n \times M} X_{m \times K}

(3)

We solve the equation as an inverse problem, and then obtain $a_{i j}$ as an unknown parameter, where $a_{i j}$ mainly represents the effect of $x_{i}$ and $y_{j}$ , even if $a_{i j}$ is affected by the levels of other $x_{i}^{’} s (i ’ \neq i)$ .

Computational method

In this section, we present a computational scheme for the identification of the relationships between miRNAs and their target mRNAs. The computational method is composed of 3 components—the optimization routine, the objective function, and the direct code. Among them, the optimization routine is the principal component of the computational method, the role of which is to obtain a new set of parameter estimations by solving the inverse problem shown in Equation 3. In particular, we employ the BFGS method for optimization, the so-called quasi-Newton method, which requires second derivatives of the objective function and thereby makes a quadratic convergent to the minimum of error norm, coupled with a drastic reduction of the computational burden. The second method is the Powell method known as direction set method, in which no such second derivatives of the objective function are required. Both methods are useful for multi-dimensional optimization, but they do not work successfully in all cases. After testing both methods with small nodes, we selected an appropriate one for the current cases. Next, the objective function is to provide the criterion for further processing to the next iteration on the basis of iterative error norms. Regarding the error norm, we employed the $L_{1}$ and $L_{2}$ norms, which are as follows:

\begin{array}{l} L_{1} = | y (x, a_{i j}) - y_{r e a l} | \\ L_{2} = {(y (x, a_{i j}) - y_{r e a l})}^{2} \end{array}

where $y (x, a_{i j})$ is numerically computed data and $y_{r e a l}$ is the experimentally obtained expression data. The $L_{1}$ norm worked successfully for the BFGS method and the $L_{2}$ norm worked better for Powell method. Finally, the direct solver is used to generate the computational data by solving the proposed mathematical model. Those computational data are then compared with the experimental data in an objective function. As both cancer and normal data were employed in this study, a note was that the predicted relations using cancer data are called $P_{c a n c e r}$ , and the predicted relations using normal data are called $P_{n o r m a l}$ . Then, the relationships between the cancer and normal data sets were computed and compared with each other. Therefore, the numerical scheme is composed of 3 individual algorithms, such as 2 subroutines and 1 main routine, as shown in Figure 1. In the main routines, we compared the inhibitory relations of the normal data with the inhibitory relations of the cancer data, and then calculated the comparison values (CVs) using relative error as follows:

C V = \frac{(P_{n o r m a l} - P_{c a n c e r})}{| P_{n o r m a l} |}

Figure 1.

Comparison of real data set with reconstructed data set using obtained parameters from Powell and Broyden-Fletcher-Goldfarb-Shannon methods with $L_{1}$ , and $L_{2}$ norms.

As CV becomes bigger, the relation is proportional to the strength of the relation between mRNA and miRNA. With the normal and cancer expression data sets, we calculated the unknown parameters as inhibitory relations in the subroutine. The overall numerical scheme of the proposed algorithm follows.

Data setting

We extracted experimentally known 26 miRNAs and 26 mRNAs from the RNA expression profiles of human cancers reported by Lu et al.¹⁵ In colon cancers, each gene consists of 4 normal and 7 cancer data points. In breast cancer, each gene consists of 3 normal and 6 cancer data points. As the data sets are much smaller than unknown parameters, malpositioning problem frequently occurs in the algorithm. For filtering analysis, we used PicTar (https://pictar.mdc-berlin.de/) and miRanda (http://www.microrna.org). PicTar is a computational method used to identify common targets of miRNAs, based on scanning multiple alignments of 3′ UTR sequences and a search set of miRNA sequences, followed by the scoring of the overlapping positions combining the PicTar scores of orthologous transcripts.¹¹ miRanda is a miRNA target prediction algorithm that optimizes sequence complementarity using position-specific rules and relies on strict interspecies conservation requirements.¹⁴ Neither sequence filtering has any relation with cancers.

Main routine

Main input: data sets

Call subroutine

Input : miRNA cancer, miRNA normal mRNA cancer, and mRNA normal data sets.

Output: $P_{n o r m a l} \cdot P_{c a n c e r}$

2. Compare the values of parameters

$\frac{(P_{n o r m a l} - P_{c a n c e r})}{| P_{n o r m a l} |} > 0.5$

Output represents the relationship between miRNAs and their target mRNAs

Subroutine

Set $g t o l = 10^{- 14}$

Read experimental data of the expression profiles of mRNAs and miRNAs;

Set initial guesses to zero

Construct linear model of Equation 2, and generate numerical data using miRNAs:D

Run Optimization method (BFGS)

Read expression microarray data profiles (normal:N or cancer:C)

Implement Objective function: $f (e r r) = | D - N (or C) |$

$I f f (e r r) > g t o l$ then go to 2

$I f f (e r r) < g t o l$ then

Output

$P_{n o r m a l}$ : Parameters from normal data

$P_{c a n c e r}$ : Parameters from cancer data

Results

Figure 1 represents the comparisons of the real and reconstructed data based on the parameters derived from the Powell and BFGS methods. We tested 2 optimization methods: Powell and BFGS, with 4 nodes of mRNAs (MAPK14, E2F1, HMGA2, and NOTCH1) and 4 nodes of miRNAs (miR-34a, miR-17a, miR-24, and miR-30) from data points of the normal colon tissues. We also assessed each method with the L₁ and L₂ objective functions. The BFGS method with the $L_{1}$ and $L_{2}$ norms successfully reconstructed the real data points. The Powell method worked well with the $L_{2}$ norm for the reconstruction but not with $L_{1}$ , which suggests that more numerical studies will be necessary to evaluate the performance of Powell method. When the BFGS method with the $L_{1}$ and $L_{2}$ norms was used for large nodes in the normal breast data, an error term of 0.011 was obtained with $L_{2}$ norm, thereby indicating that its accuracy was quite poor as compared with that of 5.97E-7 with the $L_{1}$ norm. In this study, therefore, we employed the BFGS method with $L_{1}$ . Figure 2 shows the frequency of negative relationships with each of the presented miRNAs in breast cancer by cutting the PAOM score 5.

Figure 2.

Distribution of the number of mRNAs of each miRNA with 5 PAOM in breast cancer. The most significant genes are miR-23a, miR-223, and miR-20. mRNA indicates messenger RNA; miRNA, micro-RNA; PAOM, prediction analysis by optimization method.

The frequency of miR-20, miR-23a, and miR-223 are relatively high, which suggests that those miRNAs do exert some effect on breast cancer.^16,17 Zhang et al⁶ presented miR-20 with the copy number lost in breast cancer. Figure 3 represents the frequency of negative relation with each of the colon cancer miRNAs by cutting PAOM score 5. The frequencies of miR-17, miR-30a, and miR-124 were highest, which implies that those miRNAs exert an effect on colon cancer. Recently, Monzo et al¹² observed that miR-17 was detected in human colon cancer development, and Silber et al¹⁸ observed upregulated miR-124 in colon (HCT-116) cancer cell lines. Tables 1 and 2 show the miRNAs-mRNAs relationships using the proposed method with filtering sequence analysis.

Figure 3.

Distribution of the number of mRNAs of each miRNA with 5.0 PAOM in colon cancer. miR-124, miR-30a, and miR-17 are the most significant. mRNA indicates messenger RNA; miRNA, micro-RNA; PAOM, prediction analysis by optimization method.

Table 1.

The relation of miRNAs and target mRNAs with sequence analysis of breast cancer, experimental analysis, and optimization analysis.

mRNA	miRNA	miRandaa	PicTar^a	PAOM^b
MAPK14	miR-124	O	o	5.5
CLOCK	miR-141	O		7.9
NF2	miR-141	O		7.6
	miR-23a	O		5.0
	miR-27a	O		6.0
	miR-15	O	o	56.4
NFIA	miR-155	O		41.5
HMGA2	miR-20	O	o	10.5
THBS1	miR-19a	O		25.5
CXCL12	miR-141		o	8.3
E2F1	miR-20	O	o	57.3
NOTCH1	miR-155	O		8.5
SERP1	miR-223	O		5.5
CD24	let-7a	O		48.2
	miR-98a	O		12.0
POLR2	miR-223	O		6.1
	miR-98a	O		24.0
	miR-16	O		5.1

Abbreviations: mRNA, messenger RNA; miRNA, micro-RNA; PAOM, prediction analysis by optimization method.

Predicted targets based on sequence analysis.

Predicted targets based on the proposed method.

Table 2.

The relation of miRNAs and target mRNAs with sequence analysis, experimental analysis, and the proposed optimization analysis based on colon cancer.

mRNA	miRNA	miRanda^a	PicTar^a	PAOM^b
MAKP14	miR-141	O		8.9
	miR-124	O	o	6.6
CLOCK	miR-30a	O	o	15.2
	miR-17	O	o	10.0
	miR-20	O	o	5.0
HMGA2	miR-30a	O	o	5.0
THBS1	miR-155	O		7.5
CXCL12	miR-124	O		6.1
HRAS	let-7a	O		8.7
SIP1	miR-21	O		7.5
NOTCH1	miR-155	O		5.1
MTPN	miR-223	O	o	5.3
E2F3	miR-124	O		11.1
SERP1	miR-223	O		12.5
DVL2	miR-125a	O		20.1

Abbreviations: mRNA, messenger RNA; miRNA, micro-RNA; PAOM, prediction analysis by optimization method.

Predicted targets based on sequence analysis.

Predicted targets based on the proposed methods.

Prediction of breast cancer

Table 1 shows the predicted relationships between miRNAs and mRNAs with sequence analysis in breast cancer. Based on our proposed method, we acquired 61 strong candidates out of a total of 676 miRNA-mRNA relations. With the integration of filtering sequence analysis, we predicted 18 miRNA-mRNA relationships in the breast cancer cells. Filtering analysis include sequence conservation-based miRanda and PicTar applications.

Those filtering resources are not generally associated with breast cancer. However, the current methods predicted some relationships using normal breast and cancer breast expression data. The results of our proposed method indicate that let-7a, miR-223, miR-98a, and miR-34a downregulate CD24 expression in cases of breast cancer. With sequence analysis, we suggest that let-7a (PAOM score 48.2) and miR-98a (PAOM score 12.0) are predicted to be strongly associated with CD24 in breast cancer. Those relationships were experimentally verified.^19,20 Kaipparettu et al¹⁹ found that CD24 expression was downregulated by estrogen in breast cancer stem cells. Verghese et al²⁰ found that let-7 family was downregulated significantly in breast tumor-initiating cells (BT-1Cs) that were enriched with CD24. Dai et al²¹ showed that NF2 was a tumor suppressor gene in human breast cancer. The findings of the current study reveal that miR-24, miR-141, miR-23a, miR-19a, miR-27a, and miR-15a are involved in the downregulation of NF2 expression in human breast cancer cells. With further sequence analysis, we found that miR-141, miR-23a, and miR-27a may be involved in the downregulation of NF2 expression in breast cancer cells. Recent findings have shown that miR-20 regulates E2F1 negatively.^8,22 Recently, Yu et al¹⁷ discovered a novel regulatory mechanism of breast cancer involving miR-20, which we also found here in PAOM 57.3. The following target genes have yet been verified by their miRNAs in breast cancer; however, the supporting evidences suggest that the relations are strongly associated with the breast cancer. CXCL12 expression is downregulated in primary breast carcinomas.^23,24 Inactivation of the product of MAPK14 via PPM1D overexpression was also previously discovered in breast tumor cell lines.²⁵ Supporting evidence for the role of NOTCH1 in breast cancer cells involves the fact that the rate of NOTCH1 expression in human breast cancer was found to be significantly higher than those of normal breast tissues at the margin of tumor sections.²⁶ Zang et al²⁷ showed that Notch signaling is overexpressed and highly activated in breast cancer. HMGA2 has been reported to be expressed in invasive and non-invasive breast cell lines.²⁸ Yang et al²⁹ experimentally verified that a core circadian CLOCK gene evidences tumor suppression properties and is downregulated in human breast cancer cells. Overexpression of THBS1 (TSP1) was detected in breast carcinoma and melanoma cells by interferon (IFN)-gamma-differentiated U937 cells in vitro via the release of reactive oxygen species.³⁰ Finally, the expression of SERP1 was suppressed in papillary thyroid cancer.³¹ However, SERP1 and POLR2 have not yet been identified in breast cancer cells.

Prediction of colon cancer

Table 2 shows the prediction of 15 relationships in colon cancer using the proposed method with filtering sequence analysis. HMGA2 (high mobility group [HMGI]) was observed to be abundantly expressed in human colorectal carcinomas.^32–34 MAPK14 that is regulated by miR-124^11,35 maintains a high level of ERbeta for E2 anti-proliferative effects in colon cancer cells³⁶ and in giloblastoma.³⁷ In addition, MAPK14 is involved in apoptosis in colorectal cancer induced by growth factors.³⁸ The activation of the Wnt signaling pathway appears to suppress the expression of the THBS1 gene in colon cancer cells.³⁹ Jung et al⁴⁰ observed that SIP1 (ZEB2), an E-cadherin transcriptional repressor, is induced by overexpressing TMPRSS2 in colon cancer cells, and affects the loss of E-cadherin-mediated cell-cell adhesion, resulting in an increase in cellular motility. Krugluger et al⁴¹ found that CLOCK is more abundantly expressed in colon cancer tissues than in normal tissue.

In addition, Kiriakidou et al⁴² experimentally reported that CLOCK is a target gene of miR-141. Recently, Zhang et al⁴³ demonstrated that the Notch1 signal transduction pathway mediates the effect of COX-2 selective inhibitors on colorectal cancer cells, and also discovered the mechanism of the Notch1 pathway which regulates the proliferation and apoptosis of colorectal cancer cells. E2F3, a member of E2F family, is downregulated in the HCT 116 and RKO colon cancer cell lines.⁴⁴ You et al⁴⁵ determined that with the differential expression of dishevelled segment polarity protein 2 (DVL2), the Wnt signaling pathway may contribute to colon carcinogenesis. You et al⁴⁶ reported that DVL2 was expressed in sporadic colon cancer tissues. Wendt et al⁴⁷ recently reported that the expression of CXCL12 in human colorectal carcinoma cells reduced orthotopic tumor formation and inhibited metastasis in severe combined immunodeficient mice. HRAS and MTPN are yet to be confirmed.

Discussion

Most previous computational studies have been conducted to predict miRNA-mRNA relations on the basis of DNA sequence data. The resultant large number of the sequence predictions makes biological validation quite difficult. On the contrary, a variety of previous studies have demonstrated that an miRNA deregulates its target mRNA in a cancer type-specific manner. For example, miR-34a deregulates E2F in cancer cell lines⁴⁴ whereas miR-17 deregulates E2F in breast cancer cells.²⁰ In this article, we suggested a PAOM consisting of a mathematical model and computational method using microarray data sets and filtering sequence analysis, such that the cancer-specific relationships between miRNAs and mRNAs can be predicted. The proposed PAOM was assessed and compared with normal and cancer microarray profiles in both breast and colon cancers. Among 676 relationships, we predicted 61 and 28 miRNA-mRNA relationships that might exert some effects on breast and colon cancer development, respectively. According to the results of sequence analysis filtering, we uncovered 18 breast putative relations with 12 mRNAs and 14 miRNAs and 15 colon putative relations with 12 mRNAs and 10 miRNAs. We confirmed that 8 genes—MAPK14, CLOCK, NF2, HMGA2, CXCL12, E2F1, NOTCH1, and CD24—are associated with breast cancer. Most importantly, we demonstrated that miR-20, a member of miR-17 cluster, is the target for E2F1 with PAOM score 57.3. Yu et al¹⁷ and Verghese et al²⁰ independently verified that miR-20 regulates the development of breast cancer. In addition, Hossain et al⁴⁸ confirmed that miR-17 down regulates E2F1 expression in breast cancer cells. Therefore, miR-20 is a strong candidate for targeting E2F1 mRNA in breast cancer. In particular, we predicted that CD24 is the target of let-7a and miR-98 with PAOM scores 48.2 and 12.0, respectively, in breast cancer,¹⁹ which was previously verified by Kaipparettu et al.¹⁹ Further studies will be necessary to verify those findings. In colon cancer cells, we predicted 15 relationships with 12 genes and verified 9 genes—MAPK14, CLOCK, HMGA2, THBS1, CXCL12, SIP1, NOTCH1, E2F3, and DVL2—influence on colon cancer. Overall, the proposed method described in this study was successful in the detection of some potential relationships, and may provide information for experimental studies targeting toward the identification of miRNA-mRNA relationships for specific cancers. However, there is no doubt that many unidentified relations continue to exist. Therefore, a novel approach using both computational methods and experimental validation is yet to be proposed for better outcomes in the prediction of miRNA-mRNA relationships in cancers.

Footnotes

Declaration of conflicting interests:

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding:

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF), the Ministry of Science, ICT & Future Planning (NRF-2017R1A2B4010684).

Author Contributions

SK performed all analyses, drafted the manuscript and designed the project.

References

Bartel

DP.

MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–297.

Michael

O’Connor

van Holst Pellekaan

Young

James

RJ.

Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol Cancer Res. 2003;1:882–891.

Wang

Correlation analysis of the mRNA and miRNA expression profiles in the nascent synthetic allotetraploid Raphanobrassica. Sci Rep. 2016;6:37416.

Watanabe

Tomita

Kanai

Computational methods for microRNA target prediction. Methods Enzymol. 2007;427:65–86.

Cummins

Leary

et al . The colorectal microRNAome. Proc Natl Acad Sci U S A. 2006;103:3687–3692.

Zhang

Huang

Yang

et al . microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci U S A. 2006;103:9136–9141.

Kim

Choi

Cho

KH.

Identifying the target mRNAs of microRNAs in colorectal cancer. Comput Biol Chem. 2009;33:94–99.

O’Donnell

Wentzel

Zeller

Dang

Mendell

JT.

c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005;435:839–843.

Joung

Shin

Seong

Zhang

BT.

Identification of regulatory modules by co-clustering latent variable models: stem cell differentiation. Bioinformatics. 2006;22:2005–2011.

10.

Betel

Wilson

Gabow

Marks

Sander

The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008;36:D149–D153.

11.

Krek

Grun

Poy

et al . Combinatorial microRNA target predictions. Nat Genet. 2005;37:495–500.

12.

Monzo

Navarro

Bandres

et al . Overlapping expression of microRNAs in human embryonic colon and colorectal cancer. Cell Res. 2008;18:823-833.

13.

Press

WH.

Numerical Recipes in C: The Art of Scientific Computing. 2nd ed. Cambridge, UK; New York, NY: Cambridge University Press; 1992.

14.

John

Enright

Aravin

Tuschl

Sander

Marks

DS.

Human microRNA targets. PLoS Biol. 2004;2:e363.

15.

Getz

Miska

et al . MicroRNA expression profiles classify human cancers. Nature. 2005;435:834-838.

16.

Volinia

Calin

Liu

et al . A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006;103:2257-2261.

17.

Wang

et al . A cyclin D1/microRNA 17/20 regulatory feedback loop in control of breast cancer cell proliferation. J Cell Biol. 2008;182:509–517.

18.

Silber

Lim

Petritsch

et al . miR-124 and miR-137 inhibit proliferation of glioblastoma multiforme cells and induce differentiation of brain tumor stem cells. BMC Med. 2008;6:14.

19.

Kaipparettu

Malik

Konduri

et al . Estrogen-mediated downregulation of CD24 in breast cancer cells. Int J Cancer. 2008;123:66–72.

20.

Verghese

Hanby

Speirs

Hughes

TA.

Small is beautiful: microRNAs and breast cancer—where are we now?

J Pathol. 2008;215:214–221.

21.

Dai

Wang

Sahin

Broemeling

Schutte

Pan

NHERF (Na+/H+ exchanger regulatory factor) gene mutations in human breast cancer. Oncogene. 2004;23:8681–8687.

22.

Sylvestre

De Guire

Querido

et al . An E2F/miR-20a autoregulatory feedback loop. J Biol Chem. 2007;282:2135–2143.

23.

Hsu

Chen

Westbrook

et al . CXCR4 and CXCL12 down-regulation: a novel mechanism for the chemoprotection of 3,3’-diindolylmethane for breast and ovarian cancers. Cancer Lett. 2008;265:113–123.

24.

Zhou

Jiang

Liu

et al . Down-regulation of CXCL12 mRNA expression by promoter hypermethylation and its association with metastatic progression in human breast carcinomas. J Cancer Res Clin Oncol. 2009;135:91–102.

25.

Bulavin

Demidov

Saito

et al . Amplification of PPM1D in human tumors abrogates p53 tumor-suppressor activity. Nat Genet. 2002;31:210–215.

26.

Dong

et al . Notch signaling in human breast cancer. Zhonghua Zhong Liu Za Zhi. 2007;29:425–428.

27.

Zang

et al . A study on Notch signaling in human breast cancer. Neoplasma. 2007;54:304–310.

28.

Fabjani

Tong

Wolf

et al . HMGA2 is associated with invasiveness but not a suitable marker for the detection of circulating tumor cells in breast cancer. Oncol Rep. 2005;14:737–741.

29.

Yang

Wood

Du-Quiton

Ansell

Hrushesky

WJ.

Down regulation of circadian clock gene period 2 accelerates breast cancer growth by altering its daily growth rhythm. Breast Cancer Res Treat. 2009;117:423–431.

30.

Martin-Manso

Galli

Ridnour

Tsokos

Wink

Roberts

DD.

Thrombospondin 1 promotes tumor macrophage recruitment and enhances tumor cell cytotoxicity of differentiated U937 cells. Cancer Res. 2008;68:7090–7099.

31.

Puskas

Juhasz

Zarva

Hackler

Farid

NR.

Gene profiling identifies genes specific for well-differentiated epithelial thyroid tumors. Cell Mol Biol (Noisy-le-grand). 2005;51:177–186.

32.

Fedele

Bandiera

Chiappetta

et al . Human colorectal carcinomas express high levels of high mobility group HMGI(Y) proteins. Cancer Res. 1996;56:1896–1901.

33.

Ohmachi

Tanaka

Mimori

Inoue

Yanaga

Mori

Clinical significance of TROP2 expression in colorectal cancer. Clin Cancer Res. 2006;12:3057–3063.

34.

Helmke

Markowski

Meyer

Bullerdiek

The expression of HMGA2 varies strongly among colon carcinomas. Anticancer Res. 2012;32:1589–1593.

35.

Lim

Lau

Garrett-Engele

et al . Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433:769–773.

36.

Caiazza

Merritt

Brothers

O’Toole

GA.

Inverse regulation of biofilm formation and swarming motility by Pseudomonas aeruginosa PA14. J Bacteriol. 2007;189:3603–3612.

37.

Mucaj

Lee

Skuli

et al . MicroRNA-124 expression counteracts pro-survival stress responses in glioblastoma. Oncogene. 2015;34:2204–2214.

38.

Fang

Richardson

BC.

The MAPK signalling pathways and colorectal cancer. Lancet Oncol. 2005;6:322–327.

39.

Mizukami

Duerr

Zukerberg

Chung

DC.

Wnt signaling can repress thrombospondin-1 expression in colonic tumorigenesis. Cancer Biol Ther. 2005;4:1361–1366.

40.

Jung

Lee

Park

et al . TMPRSS4 promotes invasion, migration and metastasis of human tumor cells by facilitating an epithelial-mesenchymal transition. Oncogene. 2008;27:2635–2647.

41.

Krugluger

Brandstaetter

Kallay

et al . Regulation of genes of the circadian clock in human colon cancer: reduced period-1 and dihydropyrimidine dehydrogenase transcription correlates in high-grade tumors. Cancer Res. 2007;67:7917–7922.

42.

Kiriakidou

Nelson

Kouranov

et al . A combined computational-experimental approach predicts human microRNA targets. Genes Dev. 2004;18:1165–1178.

43.

Zhang

Bai

Wang

The COX-2 selective inhibitor-independent COX-2 effect on colon carcinoma cells is associated with the Delta1/Notch1 pathway. Dig Dis Sci. 2008;53:2195–2203.

44.

Tazawa

Tsuchiya

Izumiya

Nakagama

Tumor-suppressive miR-34a induces senescence-like growth arrest through modulation of the E2F pathway in human colon cancer cells. Proc Natl Acad Sci U S A. 2007;104:15472–15477.

45.

You

Nguyen

Albers

Lin

Holcombe

RF.

Wnt pathway-related gene expression in inflammatory bowel disease. Dig Dis Sci. 2008;53:1013–1019.

46.

You

Bryant

Jurnak

Holcombe

RF.

Expression of Wnt pathway components frizzled and disheveled in colon cancer arising in patients with inflammatory bowel disease. Oncol Rep. 2007;18:691–694.

47.

Wendt

Drury

Vongsa

Dwinell

MB.

Constitutive CXCL12 expression induces anoikis in colorectal carcinoma cells. Gastroenterology. 2008;135:508–517.

48.

Hossain

Kuo

Saunders

GF.

Mir-17-5p regulates breast cancer cell proliferation by inhibiting translation of AIB1 mRNA. Mol Cell Biol. 2006;26:8191–8201.