# Module Network Inference from a cancer gene expression data set identifies microRNA regulated modules

Eric Bonnet, Marianthi Tatari, Anagha Joshi, Tom Michoel, Kathleen

Marchal, Geert Berx, Yves van de Peer

► To cite this version:

Eric Bonnet, Marianthi Tatari, Anagha Joshi, Tom Michoel, Kathleen Marchal, et al.. Module Network Inference from a cancer gene expression data set identifies microRNA regulated modules. PLoS ONE, 2010, 5 (4), pp.e10162. 10.1371/journal.pone.0010162. cea-04497253

HAL Id: cea-04497253

https://cea.hal.science/cea-04497253

Submitted on 9 Mar 2024

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Module Network Inference from a cancer gene expression data set identifies microRNA regulated modules

Eric Bonnet, Marianthi Tatari, Anagha Joshi, Tom Michoel, Kathleen

Marchal, Geert Berx, Yves van de Peer

► To cite this version:

Eric Bonnet, Marianthi Tatari, Anagha Joshi, Tom Michoel, Kathleen Marchal, et al.. Module Network Inference from a cancer gene expression data set identifies microRNA regulated modules. PLoS ONE, 2010, 5 (4), pp.e10162. 10.1371/journal.pone.0010162. cea-04497253

HAL Id: cea-04497253

https://cea.hal.science/cea-04497253v1

Submitted on 9 Mar 2024

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Module Network Inference from a Cancer Gene Expression Data Set Identifies MicroRNA Regulated Modules

Eric Bonnet¹,², Marianthi Tatari³,⁴, Anagha Joshi¹,², Tom Michoel¹,², Kathleen Marchal⁵, Geert Berx³,⁴, Yves Van de Peer¹,²*

1 Department of Plant Systems Biology, VIB, Gent, Belgium, 2 Department of Molecular Genetics, Ghent University, Gent, Belgium, 3 Unit of Molecular and Cellular Oncology, Department for Molecular Biomedical Research, VIB, Gent, Belgium, 4 Department of Biomedical Molecular Biology, Ghent University, Gent, Belgium, 5 CMPG, Department Microbial and Molecular Systems, KULeuven, Leuven, Belgium

## Abstract

**Background:** MicroRNAs (miRNAs) are small RNAs that recognize and regulate mRNA target genes. Multiple lines of evidence indicate that they are key regulators of numerous critical functions in development and disease, including cancer. However, defining the place and function of miRNAs in complex regulatory networks is not straightforward. Systems approaches, like the inference of a module network from expression data, can help to achieve this goal.

**Methodology/Principal Findings:** During the last decade, much progress has been made in the development of robust and powerful module network inference algorithms. In this study, we analyze and assess experimentally a module network inferred from both miRNA and mRNA expression data, using our recently developed module network inference algorithm based on probabilistic optimization techniques. We show that several miRNAs are predicted as statistically significant regulators for various modules of tightly co-expressed genes. A detailed analysis of three of those modules demonstrates that the specific assignment of miRNAs is functionally coherent and supported by literature. We further designed a set of experiments to test the assignment of miR-200a as the top regulator of a small module of nine genes. The results strongly suggest that miR-200a is regulating the module genes via the transcription factor ZEB1. Interestingly, this module is most likely involved in epithelial homeostasis and its dysregulation might contribute to the malignant process in cancer cells.

**Conclusions/Significance:** Our results show that a robust module network analysis of expression data can provide novel insights of miRNA function in important cellular processes. Such a computational approach, starting from expression data alone, can be helpful in the process of identifying the function of miRNAs by suggesting modules of co-expressed genes in which they play a regulatory role. As shown in this study, those modules can then be tested experimentally to further investigate and refine the function of the miRNA in the regulatory network.

**Citation:** Bonnet E, Tatari M, Joshi A, Michoel T, Marchal K, et al. (2010) Module Network Inference from a Cancer Gene Expression Data Set Identifies MicroRNA Regulated Modules. PLoS ONE 5(4): e10162. doi:10.1371/journal.pone.0010162

**Editor:** Simon Rogers, University of Glasgow, United Kingdom

**Received** December 15, 2009; **Accepted** March 22, 2010; **Published** April 14, 2010

**Copyright:** © 2010 Bonnet et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work is supported by an Innovatie door Wetenschap en Technologie (IWT) grant for the Strategisch BasisOnderzoek (SBO) Bioframe project and by an Interuniversity Attraction Pole (IUAP) grant for the BioMaGNet project (Bioinformatics and Modelling: from Genomes to Networks, ref. p6/25). The funders had no role in study design, datacollection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: yves.vandepeer@psb.vib-ugent.be

## Introduction

MicroRNAs (miRNAs) are small endogenous regulatory RNAs, present in a wide variety of eukaryotic organisms. They are incorporated into an RNA induced silencing complex (RISC) that binds to sites of variable complementarity in target messenger RNAs, triggering their degradation and/or repressing their translation [1]. Evidence for the participation of miRNAs in cell growth, cell differentiation and cancer is currently piling up. Nearly half of the annotated human miRNAs map within fragile chromosomal regions, which are areas associated with various types of human cancers. Recent evidence indicates that miRNAs as well as the factors that participate in miRNA biogenesis may function as tumor suppressors and/or oncogenes [2]. According to the latest miRBase repository release [3], there are >700

human mature miRNA sequences identified with experimental support, while some computational studies expand this list to more than 1,000 [3], roughly equaling the number of transcription factors [4]. Computational and experimental studies have also predicted that between 30% and 100% of the human protein coding genes might be under the post-transcriptional regulation of miRNAs [5,6]. It is not difficult to see that even by taking the most conservative values, the regulatory network induced by such a large number of regulators and targets is potentially extremely large. Furthermore, miRNAs do not act in isolation, but are part of a complex regulatory network, involving transcription factors, signal transducers and other types of regulatory molecules [7]. Reconstructing and analyzing such regulatory networks is thus a complex but crucial challenge to tackle.

<figure><img src="image_2.png" /></figure>
# Module Network Inference from a Cancer Gene Expression Data Set Identifies MicroRNA Regulated Modules

Eric Bonnet¹,², Marianthi Tatari³,⁴, Anagha Joshi¹,², Tom Michoel¹,², Kathleen Marchal⁵, Geert Berx³,⁴, Yves Van de Peer¹,²*

1 Department of Plant Systems Biology, VIB, Gent, Belgium, 2 Department of Molecular Genetics, Ghent University, Gent, Belgium, 3 Unit of Molecular and Cellular Oncology, Department for Molecular Biomedical Research, VIB, Gent, Belgium, 4 Department of Biomedical Molecular Biology, Ghent University, Gent, Belgium, 5 CMPG, Department Microbial and Molecular Systems, KULeuven, Leuven, Belgium

## Abstract

**Background:** MicroRNAs (miRNAs) are small RNAs that recognize and regulate mRNA target genes. Multiple lines of evidence indicate that they are key regulators of numerous critical functions in development and disease, including cancer. However, defining the place and function of miRNAs in complex regulatory networks is not straightforward. Systems approaches, like the inference of a module network from expression data, can help to achieve this goal.

**Methodology/Principal Findings:** During the last decade, much progress has been made in the development of robust and powerful module network inference algorithms. In this study, we analyze and assess experimentally a module network inferred from both miRNA and mRNA expression data, using our recently developed module network inference algorithm based on probabilistic optimization techniques. We show that several miRNAs are predicted as statistically significant regulators for various modules of tightly co-expressed genes. A detailed analysis of three of those modules demonstrates that the specific assignment of miRNAs is functionally coherent and supported by literature. We further designed a set of experiments to test the assignment of miR-200a as the top regulator of a small module of nine genes. The results strongly suggest that miR-200a is regulating the module genes via the transcription factor ZEB1. Interestingly, this module is most likely involved in epithelial homeostasis and its dysregulation might contribute to the malignant process in cancer cells.

**Conclusions/Significance:** Our results show that a robust module network analysis of expression data can provide novel insights of miRNA function in important cellular processes. Such a computational approach, starting from expression data alone, can be helpful in the process of identifying the function of miRNAs by suggesting modules of co-expressed genes in which they play a regulatory role. As shown in this study, those modules can then be tested experimentally to further investigate and refine the function of the miRNA in the regulatory network.

**Citation:** Bonnet E, Tatari M, Joshi A, Michoel T, Marchal K, et al. (2010) Module Network Inference from a Cancer Gene Expression Data Set Identifies MicroRNA Regulated Modules. PLoS ONE 5(4): e10162. doi:10.1371/journal.pone.0010162

**Editor:** Simon Rogers, University of Glasgow, United Kingdom

**Received** December 15, 2009; **Accepted** March 22, 2010; **Published** April 14, 2010

**Copyright:** © 2010 Bonnet et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work is supported by an Innovatie door Wetenschap en Technologie (IWT) grant for the Strategisch BasisOnderzoek (SBO) Bioframe project and by an Interuniversity Attraction Pole (IUAP) grant for the BioMaGNet project (Bioinformatics and Modelling: from Genomes to Networks, ref. p6/25). The funders had no role in study design, datacollection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: yves.vandepeer@psb.vib-ugent.be

## Introduction

MicroRNAs (miRNAs) are small endogenous regulatory RNAs, present in a wide variety of eukaryotic organisms. They are incorporated into an RNA induced silencing complex (RISC) that binds to sites of variable complementarity in target messenger RNAs, triggering their degradation and/or repressing their translation [1]. Evidence for the participation of miRNAs in cell growth, cell differentiation and cancer is currently piling up. Nearly half of the annotated human miRNAs map within fragile chromosomal regions, which are areas associated with various types of human cancers. Recent evidence indicates that miRNAs as well as the factors that participate in miRNA biogenesis may function as tumor suppressors and/or oncogenes [2]. According to the latest miRBase repository release [3], there are >700

human mature miRNA sequences identified with experimental support, while some computational studies expand this list to more than 1,000 [3], roughly equaling the number of transcription factors [4]. Computational and experimental studies have also predicted that between 30% and 100% of the human protein coding genes might be under the post-transcriptional regulation of miRNAs [5,6]. It is not difficult to see that even by taking the most conservative values, the regulatory network induced by such a large number of regulators and targets is potentially extremely large. Furthermore, miRNAs do not act in isolation, but are part of a complex regulatory network, involving transcription factors, signal transducers and other types of regulatory molecules [7]. Reconstructing and analyzing such regulatory networks is thus a complex but crucial challenge to tackle.

<figure><img src="image_2.png" /></figure>
Various algorithms exist to infer regulatory networks from expression data [8,9,10]. One of the most powerful methods, especially for eukaryotic organisms, assumes a modular structure of the underlying regulatory network, where a group of co-expressed genes is regulated by a common set of regulators (also known as the regulatory program) [10]. The regulatory program uses the expression levels of the set of regulators to predict the condition-dependent mean expression of the co-expressed genes. Thus, modules are composed of clusters of co-expressed genes together with their associated regulators. As a regulator can be associated with more than one module, the ensemble forms a module network. We have recently developed a novel algorithm which extends the original module network concept of Segal and co-workers [10] by using probabilistic optimization techniques which enable prioritization of the statistically most significant clusters of co-expressed genes and their candidate regulators [11,12]. The main advantage of this algorithm is that it extracts more representative centroid-like solutions from an ensemble of possible statistical models, in order to avoid suboptimal solutions. By testing it on various biological datasets, we have shown that this approach generates more coherent modules, and that regulators consistently assigned to a module are more often supported by external sources of data [11,12].

In this study, we have adapted our module network algorithm to take as input a heterogeneous dataset of both miRNA and messenger RNA (mRNA) expression data measured on the same samples. Multiple miRNAs are assigned as high-scoring candidate regulators for several modules, together with well-known transcription factors or signal transducers. A detailed analysis of three modules where miRNAs are selected as high-scoring regulators shows that this assignment is highly coherent with the module function and is also supported by various external sources of data. We have also validated one of those modules experimentally, showing that over-expression or inhibition of the miRNA assigned as a regulator changes significantly the expression of a selection of module genes, thereby confirming the inference of the algorithm. Those results corroborate module network inference as a robust and useful approach to gain more precise insights into miRNA function.

# Results

## Inference of a microRNA module network from expression data

The LeMoNe algorithm, starting from an expression data matrix and a list of candidate regulators, will produce a module network, composed of modules of co-expressed genes and their associated regulators. The algorithm is also clustering the conditions (columns) for each set of co-expressed genes, creating condition clusters. The list of regulators for a given module is ordered according to their individual score. This score only takes into account the differential expression of the regulator across the different different condition clusters, and not their absolute value. This way we can simultaneously evaluate and compare mRNA and miRNA candidate regulators, using the expression levels of each class of regulators. As input for our algorithm, we used a dataset composed of expression data measured on 89 tumor and normal tissue samples (representing 11 tumor classes) both for 11,833 messenger RNAs and 124 miRNAs [13]. Unlike previous attempts [14,15,16], our approach for the integration of miRNAs in the network is not based on miRNA target prediction, or a mixture between target prediction and expression data, but relies solely on expression data. The algorithm generated a set of 76 tightly co-expressed gene clusters, corresponding to a total of 2,987

genes. We calculated the GO enrichment for all the modules [17] and found a total of 44 clusters having at least one GO category enriched ($p<0.05$), for a total of 589 enriched categories (the complete list of modules and their GO categories is available as table S2). For the assignment of regulators, we compiled a list of 1,841 candidate regulators based on their GO annotation (either transcription factors or signal transducers), plus the list of 124 miRNAs. After the assignment of regulators we took a stringent cutoff corresponding to the top 2% most significant predicted regulatory interactions (Figure 1), obtaining a final set of 294 unique regulators (the complete list of module genes and regulators is available in table S1). Within this set, ten miRNAs were selected as regulators for a total of seven modules (Table 1). In order to assess the validity and the relevance of the inferred module network, we here present a detailed analysis of three modules, with an emphasis on the typical features of miRNA mediated regulatory modules. Those modules were selected based on their intrinsic interest, functional coherence and the high number of literature references discussing their putative function.

## MiR-133 and miR-145 are assigned as regulators of a smooth muscle actomyosin module

Module 29 is a small module composed of four genes and five assigned regulators (Figure 2a). The GO over-represented categories for this module are linked to smooth muscle development and actomyosin structure (Figure 2a and table S2). MYH11 encodes a smooth muscle myosin heavy chain family member. ACTG2 is a gamma 2 actin protein found in enteric tissues. The two other genes in the module (MYLK and CNN1) are well known regulators of the actin-myosin interactions. MYLK is the myosin light chain kinase, a dedicated calcium-dependent kinase that phosphorylates a specific site on the regulatory light chain of the myosin, enhancing its activity. MYLK is ubiquitous in all adult tissues with the highest amounts found in smooth muscle tissues [18]. CNN1 (calponin) is a calcium binding protein that inhibits the ATPase activity of the myosin in smooth muscle. The top regulator (PPP1R12B) selected for this module is a myosin phosphatase subunit. The myosin phosphatase is also a well known core regulator of the actomyosin pathway, inhibiting the myosin activity [18]. The second high scoring candidate regulator is a miRNA, miR-133, while the third regulator is the TGF-beta stimulated clone-22 member 1 (TSC22D1) gene, which encodes a leucine zipper domain protein, a member of the TGF-beta1 pathway which is involved in the regulation of transcription. The last regulators are ANGPTL2, a vascular endothelial growth factor, and another miRNA, miR-145.

The two miRNAs selected as regulators for this module clearly show a tightly positively correlated expression pattern with the module genes (Figure 2a). As most miRNAs have been characterized so far as negative regulators of gene expression, this suggest an indirect regulation between the miRNAs and module 29 genes. Recent studies reveal several likely candidate genes that could act as intermediate regulators between the miRNAs and module 29 genes. MiR-133, selected as the second best regulator for this module (Figure 2a), was recently shown to be a key regulator for skeletal muscle development and cardiac muscle hypertrophy [19,20]. In those studies, miR-133 has been shown to directly regulate the SRF transcription factor. SRF is recognized as a vital factor for normal cytoskeletal and contractile cell activities and all the module 29 genes (*MYH11*, *CNN1*, *ACTG2*, *MYLK*) are known to be direct targets of SRF [21]. Those literature results support the hypothesis of an indirect regulatory link between miR-133a and module 29 genes via SRF. Most studies on SRF activity have so far characterized this factor as a

<figure><img src="image_3.png" /></figure>
Various algorithms exist to infer regulatory networks from expression data [8,9,10]. One of the most powerful methods, especially for eukaryotic organisms, assumes a modular structure of the underlying regulatory network, where a group of co-expressed genes is regulated by a common set of regulators (also known as the regulatory program) [10]. The regulatory program uses the expression levels of the set of regulators to predict the condition-dependent mean expression of the co-expressed genes. Thus, modules are composed of clusters of co-expressed genes together with their associated regulators. As a regulator can be associated with more than one module, the ensemble forms a module network. We have recently developed a novel algorithm which extends the original module network concept of Segal and co-workers [10] by using probabilistic optimization techniques which enable prioritization of the statistically most significant clusters of co-expressed genes and their candidate regulators [11,12]. The main advantage of this algorithm is that it extracts more representative centroid-like solutions from an ensemble of possible statistical models, in order to avoid suboptimal solutions. By testing it on various biological datasets, we have shown that this approach generates more coherent modules, and that regulators consistently assigned to a module are more often supported by external sources of data [11,12].

In this study, we have adapted our module network algorithm to take as input a heterogeneous dataset of both miRNA and messenger RNA (mRNA) expression data measured on the same samples. Multiple miRNAs are assigned as high-scoring candidate regulators for several modules, together with well-known transcription factors or signal transducers. A detailed analysis of three modules where miRNAs are selected as high-scoring regulators shows that this assignment is highly coherent with the module function and is also supported by various external sources of data. We have also validated one of those modules experimentally, showing that over-expression or inhibition of the miRNA assigned as a regulator changes significantly the expression of a selection of module genes, thereby confirming the inference of the algorithm. Those results corroborate module network inference as a robust and useful approach to gain more precise insights into miRNA function.

# Results

## Inference of a microRNA module network from expression data

The LeMoNe algorithm, starting from an expression data matrix and a list of candidate regulators, will produce a module network, composed of modules of co-expressed genes and their associated regulators. The algorithm is also clustering the conditions (columns) for each set of co-expressed genes, creating condition clusters. The list of regulators for a given module is ordered according to their individual score. This score only takes into account the differential expression of the regulator across the different different condition clusters, and not their absolute value. This way we can simultaneously evaluate and compare mRNA and miRNA candidate regulators, using the expression levels of each class of regulators. As input for our algorithm, we used a dataset composed of expression data measured on 89 tumor and normal tissue samples (representing 11 tumor classes) both for 11,833 messenger RNAs and 124 miRNAs [13]. Unlike previous attempts [14,15,16], our approach for the integration of miRNAs in the network is not based on miRNA target prediction, or a mixture between target prediction and expression data, but relies solely on expression data. The algorithm generated a set of 76 tightly co-expressed gene clusters, corresponding to a total of 2,987

genes. We calculated the GO enrichment for all the modules [17] and found a total of 44 clusters having at least one GO category enriched ($p<0.05$), for a total of 589 enriched categories (the complete list of modules and their GO categories is available as table S2). For the assignment of regulators, we compiled a list of 1,841 candidate regulators based on their GO annotation (either transcription factors or signal transducers), plus the list of 124 miRNAs. After the assignment of regulators we took a stringent cutoff corresponding to the top 2% most significant predicted regulatory interactions (Figure 1), obtaining a final set of 294 unique regulators (the complete list of module genes and regulators is available in table S1). Within this set, ten miRNAs were selected as regulators for a total of seven modules (Table 1). In order to assess the validity and the relevance of the inferred module network, we here present a detailed analysis of three modules, with an emphasis on the typical features of miRNA mediated regulatory modules. Those modules were selected based on their intrinsic interest, functional coherence and the high number of literature references discussing their putative function.

## MiR-133 and miR-145 are assigned as regulators of a smooth muscle actomyosin module

Module 29 is a small module composed of four genes and five assigned regulators (Figure 2a). The GO over-represented categories for this module are linked to smooth muscle development and actomyosin structure (Figure 2a and table S2). MYH11 encodes a smooth muscle myosin heavy chain family member. ACTG2 is a gamma 2 actin protein found in enteric tissues. The two other genes in the module (MYLK and CNN1) are well known regulators of the actin-myosin interactions. MYLK is the myosin light chain kinase, a dedicated calcium-dependent kinase that phosphorylates a specific site on the regulatory light chain of the myosin, enhancing its activity. MYLK is ubiquitous in all adult tissues with the highest amounts found in smooth muscle tissues [18]. CNN1 (calponin) is a calcium binding protein that inhibits the ATPase activity of the myosin in smooth muscle. The top regulator (PPP1R12B) selected for this module is a myosin phosphatase subunit. The myosin phosphatase is also a well known core regulator of the actomyosin pathway, inhibiting the myosin activity [18]. The second high scoring candidate regulator is a miRNA, miR-133, while the third regulator is the TGF-beta stimulated clone-22 member 1 (TSC22D1) gene, which encodes a leucine zipper domain protein, a member of the TGF-beta1 pathway which is involved in the regulation of transcription. The last regulators are ANGPTL2, a vascular endothelial growth factor, and another miRNA, miR-145.

The two miRNAs selected as regulators for this module clearly show a tightly positively correlated expression pattern with the module genes (Figure 2a). As most miRNAs have been characterized so far as negative regulators of gene expression, this suggest an indirect regulation between the miRNAs and module 29 genes. Recent studies reveal several likely candidate genes that could act as intermediate regulators between the miRNAs and module 29 genes. MiR-133, selected as the second best regulator for this module (Figure 2a), was recently shown to be a key regulator for skeletal muscle development and cardiac muscle hypertrophy [19,20]. In those studies, miR-133 has been shown to directly regulate the SRF transcription factor. SRF is recognized as a vital factor for normal cytoskeletal and contractile cell activities and all the module 29 genes (*MYH11*, *CNN1*, *ACTG2*, *MYLK*) are known to be direct targets of SRF [21]. Those literature results support the hypothesis of an indirect regulatory link between miR-133a and module 29 genes via SRF. Most studies on SRF activity have so far characterized this factor as a

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 1. Random and real regulators score distributions.** The histograms represent the distribution of randomly assigned (green) and true (yellow) regulators scores for the module network. The arrow for the random regulators represents the maximum score for randomly assigned regulators with the value indicated between brackets. The arrow for the true regulators represents the cutoff score value, with the raw value indicated between brackets.

doi:10.1371/journal.pone.0010162.g001</figcaption></figure>

transcriptional activator [21], but some results also suggest that SRF might act as a transcriptional repressor of its targets genes [22,23,24]. SRF mediated gene repression is not clear, but it might involve the recruitment by SRF of transcriptional silencers [23,24]. If we hypothesize that SRF is repressing the transcription of module 29 genes, then the regulatory chain miR-133 – SRF – module 29 genes can explain the positive gene expression correlation that is observed in Figure 2a. The other miRNA selected as a regulator, miR-145 was recently shown to be an important regulator for smooth muscle cell fate [25]. This study [25] also shows that miR-145 is activating one of its direct targets, the myocardin (Myocd), which is a transcription factor well known to activate smooth muscle gene expression by interacting with SRF [26]. Thus we also have a regulatory chain miR-145 – Myocd – module 29 genes that can explain the pattern of expression observed in Figure 2a. Neither SRF nor Myocd are assigned as regulators or clustered together with the other module 29 genes. Unfortunately, the myocardin gene was not present in the microarrays used to produce the datasets [13]. The profile of the SRF transcription factor appears to correlate poorly with the expression of module 29 genes in our dataset (Data file S1), explaining why this gene could not be selected as a regulator. Several reasons could explain why the profile is divergent, like post-translational modifications or the fact that miRNAs act at the post-transcriptional level, possibly preventing the regulatory effect to be detected (by repressing the translation).

## MiR-142s are assigned as regulators of an immune response module

Module 18 is composed of six genes (Figure 2b), of which five encode immunoglobulins corresponding either to the heavy chain (IGHG4, IGHA2, IGHA1) or to the light chain (IGKV1-5,

IGLV3-21), while IGLL1 is the surrogate light chain, a critical component of the pre-B cell receptor complex. Not surprisingly, we found the GO category immune response over-represented for this module (Figure 2b and table S2). All the module genes are known to be mostly expressed in developing and mature B-cells, revealing a coherent module [27]. Nine high scoring regulators were selected for this module. The top regulator is a homeobox gene, HOXC5. The HMGA1 gene is selected as the second best regulator for this module. High mobility group proteins (HMGA) regulate the activity of a wide variety of genes by changing the DNA conformation of their target genes. HMGA1 is known to co-activate transcription in B-cells and to be important for B-cells development [28]. The third and fourth candidate regulators are two miRNAs processed from the same precursor, miR-142-5p and miR-142-3p. The HLA-DRB1 gene belongs to the HLA class II beta chain paralogues. It is known to play a central role in the immune system by presenting peptides derived from extracellular proteins [29]. CCL5 is a chemotactic cytokine playing an active role in recruiting leukocytes to inflammatory sites [30]. AXL is a receptor tyrosine kinase that is transforming in fibroblast and hematopoietic cells, and is involved in mesenchymal development [31]. CXCL14 is a small cytokine belonging to the CXC chemokine family. This gene is chemotactic for monocytes and can activate these cells in the presence of an inflammatory mediator [32]. CXCL14 expression is reduced or absent from most cancer cells [33]. This module is probably linked to an immune response triggered by various tumor states. Such persistent pro-tumor immune responses are known to potentiate primary tumor development and malignant progression.

MiR-142s are preferentially expressed in hematopoietic tissues and their expression is regulated during hematopoiesis, suggesting a role in immune cells differentiation [34]. The transcription factor

<figure><img src="image_5.png" /></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 1. Random and real regulators score distributions.** The histograms represent the distribution of randomly assigned (green) and true (yellow) regulators scores for the module network. The arrow for the random regulators represents the maximum score for randomly assigned regulators with the value indicated between brackets. The arrow for the true regulators represents the cutoff score value, with the raw value indicated between brackets.

doi:10.1371/journal.pone.0010162.g001</figcaption></figure>

transcriptional activator [21], but some results also suggest that SRF might act as a transcriptional repressor of its targets genes [22,23,24]. SRF mediated gene repression is not clear, but it might involve the recruitment by SRF of transcriptional silencers [23,24]. If we hypothesize that SRF is repressing the transcription of module 29 genes, then the regulatory chain miR-133 – SRF – module 29 genes can explain the positive gene expression correlation that is observed in Figure 2a. The other miRNA selected as a regulator, miR-145 was recently shown to be an important regulator for smooth muscle cell fate [25]. This study [25] also shows that miR-145 is activating one of its direct targets, the myocardin (Myocd), which is a transcription factor well known to activate smooth muscle gene expression by interacting with SRF [26]. Thus we also have a regulatory chain miR-145 – Myocd – module 29 genes that can explain the pattern of expression observed in Figure 2a. Neither SRF nor Myocd are assigned as regulators or clustered together with the other module 29 genes. Unfortunately, the myocardin gene was not present in the microarrays used to produce the datasets [13]. The profile of the SRF transcription factor appears to correlate poorly with the expression of module 29 genes in our dataset (Data file S1), explaining why this gene could not be selected as a regulator. Several reasons could explain why the profile is divergent, like post-translational modifications or the fact that miRNAs act at the post-transcriptional level, possibly preventing the regulatory effect to be detected (by repressing the translation).

## MiR-142s are assigned as regulators of an immune response module

Module 18 is composed of six genes (Figure 2b), of which five encode immunoglobulins corresponding either to the heavy chain (IGHG4, IGHA2, IGHA1) or to the light chain (IGKV1-5,

IGLV3-21), while IGLL1 is the surrogate light chain, a critical component of the pre-B cell receptor complex. Not surprisingly, we found the GO category immune response over-represented for this module (Figure 2b and table S2). All the module genes are known to be mostly expressed in developing and mature B-cells, revealing a coherent module [27]. Nine high scoring regulators were selected for this module. The top regulator is a homeobox gene, HOXC5. The HMGA1 gene is selected as the second best regulator for this module. High mobility group proteins (HMGA) regulate the activity of a wide variety of genes by changing the DNA conformation of their target genes. HMGA1 is known to co-activate transcription in B-cells and to be important for B-cells development [28]. The third and fourth candidate regulators are two miRNAs processed from the same precursor, miR-142-5p and miR-142-3p. The HLA-DRB1 gene belongs to the HLA class II beta chain paralogues. It is known to play a central role in the immune system by presenting peptides derived from extracellular proteins [29]. CCL5 is a chemotactic cytokine playing an active role in recruiting leukocytes to inflammatory sites [30]. AXL is a receptor tyrosine kinase that is transforming in fibroblast and hematopoietic cells, and is involved in mesenchymal development [31]. CXCL14 is a small cytokine belonging to the CXC chemokine family. This gene is chemotactic for monocytes and can activate these cells in the presence of an inflammatory mediator [32]. CXCL14 expression is reduced or absent from most cancer cells [33]. This module is probably linked to an immune response triggered by various tumor states. Such persistent pro-tumor immune responses are known to potentiate primary tumor development and malignant progression.

MiR-142s are preferentially expressed in hematopoietic tissues and their expression is regulated during hematopoiesis, suggesting a role in immune cells differentiation [34]. The transcription factor

<figure><img src="image_5.png" /></figure>
<table><caption>**Table 1.** List of modules where miRNAs have been selected as high-scoring regulators.</caption><thead><tr><th>Module ID</th><th>Number of genes in the module</th><th>GO categories enriched for the module (p&lt;0.05)</th><th>Regulators</th><th>Score</th></tr></thead><tbody><tr><td rowspan="10">15</td><td rowspan="10">73</td><td rowspan="10">None</td><td>MAFB</td><td>3065</td></tr><tr><td>NFIB</td><td>2774</td></tr><tr><td>MSRB2</td><td>713</td></tr><tr><td>IGFALS</td><td>464</td></tr><tr><td>NR2C1</td><td>457</td></tr><tr><td>CCL23</td><td>392</td></tr><tr><td>MORF4L1</td><td>288</td></tr><tr><td>FOS</td><td>249</td></tr><tr><td><strong>let-7c</strong></td><td>231</td></tr><tr><td><strong>miR-125b</strong></td><td>209</td></tr><tr><td rowspan="10">17</td><td rowspan="10">17</td><td rowspan="10">Digestion<br>Proteolysis<br>Lipid catabolic process</td><td>NFIB</td><td>3785</td></tr><tr><td>HEYL</td><td>2361</td></tr><tr><td>PLA2G1B</td><td>2164</td></tr><tr><td><strong>miR-216</strong></td><td>1085</td></tr><tr><td>CCRL2</td><td>780</td></tr><tr><td>MC5R</td><td>327</td></tr><tr><td>TRAF1</td><td>285</td></tr><tr><td>CDX1</td><td>261</td></tr><tr><td>DLX2</td><td>254</td></tr><tr><td>E2F1</td><td>251</td></tr><tr><td>IFNGR1</td><td>250</td></tr><tr><td>CCR9</td><td>244</td></tr><tr><td>TACR3</td><td>215</td></tr><tr><td rowspan="10">18</td><td rowspan="10">6</td><td rowspan="10">Immune response</td><td>HOXC5</td><td>1937</td></tr><tr><td>HMGA1</td><td>1007</td></tr><tr><td><strong>miR-142-5p</strong></td><td>863</td></tr><tr><td><strong>miR-142-3p</strong></td><td>709</td></tr><tr><td>HLA-DRB1</td><td>380</td></tr><tr><td>CCL5</td><td>241</td></tr><tr><td>AXL</td><td>215</td></tr><tr><td>CXCL14</td><td>214</td></tr><tr><td><strong>miR-200a</strong></td><td>1786</td></tr><tr><td>PPP1R1B</td><td>727</td></tr><tr><td>GPR30</td><td>512</td></tr><tr><td>PTGER3</td><td>456</td></tr><tr><td>ZNF157</td><td>320</td></tr><tr><td>GNB3</td><td>262</td></tr><tr><td>GNG5</td><td>209</td></tr><tr><td rowspan="5">29</td><td rowspan="5">4</td><td rowspan="5">Actomyosin structure organization and biogenesis<br>Smooth muscle contraction<br>Smooth muscle fiber development<br>Elastic fiber assembly<br>Striated muscle development</td><td>PPP1R12B</td><td>7346</td></tr><tr><td><strong>miR-133a</strong></td><td>1216</td></tr><tr><td>TSC22D1</td><td>1209</td></tr><tr><td>ANGPTL2</td><td>353</td></tr><tr><td><strong>miR-145</strong></td><td>220</td></tr><tr><td rowspan="6">35</td><td rowspan="6">20</td><td rowspan="6">None</td><td>NFATC1</td><td>2827</td></tr><tr><td>SH3BP1</td><td>2162</td></tr><tr><td>FMNL2</td><td>699</td></tr><tr><td>RASSF4</td><td>381</td></tr><tr><td>GNAQ</td><td>290</td></tr></tbody></table>

<figure><img src="image_6.png" /></figure>
<table><caption>**Table 1.** List of modules where miRNAs have been selected as high-scoring regulators.</caption><thead><tr><th>Module ID</th><th>Number of genes in the module</th><th>GO categories enriched for the module (p&lt;0.05)</th><th>Regulators</th><th>Score</th></tr></thead><tbody><tr><td rowspan="10">15</td><td rowspan="10">73</td><td rowspan="10">None</td><td>MAFB</td><td>3065</td></tr><tr><td>NFIB</td><td>2774</td></tr><tr><td>MSRB2</td><td>713</td></tr><tr><td>IGFALS</td><td>464</td></tr><tr><td>NR2C1</td><td>457</td></tr><tr><td>CCL23</td><td>392</td></tr><tr><td>MORF4L1</td><td>288</td></tr><tr><td>FOS</td><td>249</td></tr><tr><td><strong>let-7c</strong></td><td>231</td></tr><tr><td><strong>miR-125b</strong></td><td>209</td></tr><tr><td rowspan="10">17</td><td rowspan="10">17</td><td rowspan="10">Digestion<br>Proteolysis<br>Lipid catabolic process</td><td>NFIB</td><td>3785</td></tr><tr><td>HEYL</td><td>2361</td></tr><tr><td>PLA2G1B</td><td>2164</td></tr><tr><td><strong>miR-216</strong></td><td>1085</td></tr><tr><td>CCRL2</td><td>780</td></tr><tr><td>MC5R</td><td>327</td></tr><tr><td>TRAF1</td><td>285</td></tr><tr><td>CDX1</td><td>261</td></tr><tr><td>DLX2</td><td>254</td></tr><tr><td>E2F1</td><td>251</td></tr><tr><td>IFNGR1</td><td>250</td></tr><tr><td>CCR9</td><td>244</td></tr><tr><td>TACR3</td><td>215</td></tr><tr><td rowspan="10">18</td><td rowspan="10">6</td><td rowspan="10">Immune response</td><td>HOXC5</td><td>1937</td></tr><tr><td>HMGA1</td><td>1007</td></tr><tr><td><strong>miR-142-5p</strong></td><td>863</td></tr><tr><td><strong>miR-142-3p</strong></td><td>709</td></tr><tr><td>HLA-DRB1</td><td>380</td></tr><tr><td>CCL5</td><td>241</td></tr><tr><td>AXL</td><td>215</td></tr><tr><td>CXCL14</td><td>214</td></tr><tr><td><strong>miR-200a</strong></td><td>1786</td></tr><tr><td>PPP1R1B</td><td>727</td></tr><tr><td>GPR30</td><td>512</td></tr><tr><td>PTGER3</td><td>456</td></tr><tr><td>ZNF157</td><td>320</td></tr><tr><td>GNB3</td><td>262</td></tr><tr><td>GNG5</td><td>209</td></tr><tr><td rowspan="5">29</td><td rowspan="5">4</td><td rowspan="5">Actomyosin structure organization and biogenesis<br>Smooth muscle contraction<br>Smooth muscle fiber development<br>Elastic fiber assembly<br>Striated muscle development</td><td>PPP1R12B</td><td>7346</td></tr><tr><td><strong>miR-133a</strong></td><td>1216</td></tr><tr><td>TSC22D1</td><td>1209</td></tr><tr><td>ANGPTL2</td><td>353</td></tr><tr><td><strong>miR-145</strong></td><td>220</td></tr><tr><td rowspan="6">35</td><td rowspan="6">20</td><td rowspan="6">None</td><td>NFATC1</td><td>2827</td></tr><tr><td>SH3BP1</td><td>2162</td></tr><tr><td>FMNL2</td><td>699</td></tr><tr><td>RASSF4</td><td>381</td></tr><tr><td>GNAQ</td><td>290</td></tr></tbody></table>

<figure><img src="image_6.png" /></figure>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Module ID</th><th>Number of genes in the module</th><th>GO categories enriched for the module (p&lt;0.05)</th><th>Regulators</th><th>Score</th></tr></thead><tbody><tr><td rowspan="3">94</td><td rowspan="3">3</td><td rowspan="2">Folic acid metabolic process</td><td>CTNNBIP1</td><td>288</td></tr><tr><td>miR-181*</td><td>249</td></tr><tr><td>Regulation of Notch signaling pathway</td><td>miR-10a</td><td>253</td></tr><tr><td></td><td></td><td>Folic acid transport</td><td></td><td></td></tr></tbody></table>

All regulators that are above cutoff are listed for each module, and ordered according to their decreasing score value. MicroRNA genes are highlighted in bold.
doi:10.1371/journal.pone.0010162.t001

TCF12 is predicted as a target for miR-142-3p [35]. In a previous study, a combined dosage of the factors E2A, E2-2 and TCF12 was shown to be required for normal B-cell development. More precisely, TCF12 is important for the generation of normal numbers of pro-B-cells [36]. Because module 18 genes are expressed in developing B-cells, the regulation of TCF12 by miR-142-3p might be important for this process. Furthermore, we found conserved binding motifs for TCF12 for most module 18 genes (Data file S1), indicating that this transcription factor could be important for their regulation. Like for module 29 and SRF, the expression profile of TCF12 is highly divergent from the expression profiles of module 18 genes, explaining why this gene was not selected as module gene or regulator (data not shown).

## Mir-200a is a key regulator of a module involved in epithelial homeostasis

Module 25 is composed of nine genes (Figure 2c). SCNN1A (also known as ENaC) is the subunit alpha of the amiloride sensitive epithelial sodium channel, expressed in many epithelial tissues [37]. PRSS8 (prostasin) is a trypsinogen which regulates the activity of the epithelial sodium channel [38]. FDXY3 is a small membrane protein that is highly transcribed in tissues such as uterus, stomach and colon, and may function as a Na/K channel regulator [39]. TACSTD1, the tumor-associated calcium signal transducer 1, functions as a calcium-independent cell adhesion molecule [40]. Other genes like ATAD4 or TMEM63A are trans-membrane proteins of unknown function. RAB25 is a small GTP binding protein. RAB proteins have been involved in the regulation of vesicle trafficking [41]. The module top regulator is miR-200a. Regarding the other regulators, PPP1R1B is a phosphoprotein regulated by dopamine and cAMP, and is an inhibitor of the protein phosphatase 1. Besides its well-known role in the central nervous system, it is highly expressed in a variety of epithelial tissues where it might play a role in epithelial signaling and tumorigenesis [42]. GPR30 is a trans-membrane G protein coupled estrogen receptor [43], while PTGER3 is a G-protein coupled prostaglandin E2 receptor that is involved in various physiological processes and was shown to affect intracellular concentrations of Ca++ and cAMP [43]. ZNF157 is a zinc finger protein of unknown function while GNB3 and GNG5 are G proteins subunits involved in signal transduction. From the functions of these genes, we can conclude that most of the module 25 genes and regulators are likely involved in epithelial homeostasis, although we did not find any particular GO category enriched for this module. It is also worth noting that several of those genes are related to tumor progression [40,41,44].

MiR-200a, which was selected as the best candidate regulator for module 25, is a member of a miRNA family of five closely related miRNAs (miR-200a, miR-200b, miR-200c, miR-141 and

miR-429). Recent publications show epithelial-specific expression of miR-200a and miR-200b [45,46]. We designed a set of experiments to validate the role of miR-200a as a regulator of the expression of genes in module 25. MiR-200a was introduced in a human de-differentiated epithelial breast cancer cell line MDA-MB-231, known to express aberrantly low levels of miR-200a. The expression of six genes (RAB25, IRF6, SCNN1A, PRSS8, ATAD4, TACSTD1) out of nine belonging to module 25 was monitored using RT-qPCR (Figure 3). Without exception, the six monitored genes show a clear up-regulation upon exogenous expression of miR-200a (Figure 3a). The reverse experiment, inhibition of some members of the miR-200 family (miR-200a,b,c) in the MDA-MB-231 cells using antagonists, resulted in the significant down-regulation of four out of five tested genes, (SCNN1A is not significantly down-regulated, ATAD4 is not expressed in normal conditions in this cell line)(Figure 3b). Those results clearly show that miR-200a is a core regulator of module 25, most probably with other members of the miR-200 family.

In this module, we observe again a clear positive correlation pattern between miR-200a and the module genes expression, suggesting an indirect regulatory circuit between miR-200a and module 25 genes (Figure 2c and Figure 3a and 3b). Recent experimental work showed that miR-200 family members directly target the transcription factors ZEB1 and ZEB2 [47,48,49]. These transcription factors are known as major transcriptional repressors of epithelial differentiation orchestrating epithelial mesenchymal transition (EMT) [50]. EMT is a process that drives epithelial cells from a polarized phenotype to a highly motile, non polarized mesenchymal phenotype and is known to occur in epithelial tumors giving rise to highly malignant cancer cells. The ZEB transcription factors have been functionally related to members of the miR-200 family via a double negative feedback loop, thus promoting EMT and cancer invasion [47,51,52]. We found conserved ZEB binding motifs for several module 25 genes (Data file S1), suggesting that ZEB factors could be the intermediate regulators between miR-200 and module 25 genes. To test this hypothesis, we down-regulated the ZEB1 transcription factor in MDA-MD-231 cells with a specific siRNA while monitoring the expression of eight module 25 genes (Figure 3c). All the genes show a strong up-regulation pattern, with the exception of the gene RAB25 (Figure 3c). Those results demonstrate that the ZEB1 factor is essential for the regulation of module 25 genes. Taken together, our experimental results (Figure 3) strongly suggest the existence of a regulatory chain between miR-200 and module 25 genes via the ZEB1 transcription factor (Figure 4). As both miR-200 and ZEB1 play important roles in EMT [47,48,49,51,52] we hypothesize that module 25 repression might contribute to the malignant EMT process in cancer cells.

<figure><img src="image_7.png" /></figure>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Module ID</th><th>Number of genes in the module</th><th>GO categories enriched for the module (p&lt;0.05)</th><th>Regulators</th><th>Score</th></tr></thead><tbody><tr><td rowspan="3">94</td><td rowspan="3">3</td><td rowspan="2">Folic acid metabolic process</td><td>CTNNBIP1</td><td>288</td></tr><tr><td>miR-181*</td><td>249</td></tr><tr><td>Regulation of Notch signaling pathway</td><td>miR-10a</td><td>253</td></tr><tr><td></td><td></td><td>Folic acid transport</td><td></td><td></td></tr></tbody></table>

All regulators that are above cutoff are listed for each module, and ordered according to their decreasing score value. MicroRNA genes are highlighted in bold.
doi:10.1371/journal.pone.0010162.t001

TCF12 is predicted as a target for miR-142-3p [35]. In a previous study, a combined dosage of the factors E2A, E2-2 and TCF12 was shown to be required for normal B-cell development. More precisely, TCF12 is important for the generation of normal numbers of pro-B-cells [36]. Because module 18 genes are expressed in developing B-cells, the regulation of TCF12 by miR-142-3p might be important for this process. Furthermore, we found conserved binding motifs for TCF12 for most module 18 genes (Data file S1), indicating that this transcription factor could be important for their regulation. Like for module 29 and SRF, the expression profile of TCF12 is highly divergent from the expression profiles of module 18 genes, explaining why this gene was not selected as module gene or regulator (data not shown).

## Mir-200a is a key regulator of a module involved in epithelial homeostasis

Module 25 is composed of nine genes (Figure 2c). SCNN1A (also known as ENaC) is the subunit alpha of the amiloride sensitive epithelial sodium channel, expressed in many epithelial tissues [37]. PRSS8 (prostasin) is a trypsinogen which regulates the activity of the epithelial sodium channel [38]. FDXY3 is a small membrane protein that is highly transcribed in tissues such as uterus, stomach and colon, and may function as a Na/K channel regulator [39]. TACSTD1, the tumor-associated calcium signal transducer 1, functions as a calcium-independent cell adhesion molecule [40]. Other genes like ATAD4 or TMEM63A are trans-membrane proteins of unknown function. RAB25 is a small GTP binding protein. RAB proteins have been involved in the regulation of vesicle trafficking [41]. The module top regulator is miR-200a. Regarding the other regulators, PPP1R1B is a phosphoprotein regulated by dopamine and cAMP, and is an inhibitor of the protein phosphatase 1. Besides its well-known role in the central nervous system, it is highly expressed in a variety of epithelial tissues where it might play a role in epithelial signaling and tumorigenesis [42]. GPR30 is a trans-membrane G protein coupled estrogen receptor [43], while PTGER3 is a G-protein coupled prostaglandin E2 receptor that is involved in various physiological processes and was shown to affect intracellular concentrations of Ca++ and cAMP [43]. ZNF157 is a zinc finger protein of unknown function while GNB3 and GNG5 are G proteins subunits involved in signal transduction. From the functions of these genes, we can conclude that most of the module 25 genes and regulators are likely involved in epithelial homeostasis, although we did not find any particular GO category enriched for this module. It is also worth noting that several of those genes are related to tumor progression [40,41,44].

MiR-200a, which was selected as the best candidate regulator for module 25, is a member of a miRNA family of five closely related miRNAs (miR-200a, miR-200b, miR-200c, miR-141 and

miR-429). Recent publications show epithelial-specific expression of miR-200a and miR-200b [45,46]. We designed a set of experiments to validate the role of miR-200a as a regulator of the expression of genes in module 25. MiR-200a was introduced in a human de-differentiated epithelial breast cancer cell line MDA-MB-231, known to express aberrantly low levels of miR-200a. The expression of six genes (RAB25, IRF6, SCNN1A, PRSS8, ATAD4, TACSTD1) out of nine belonging to module 25 was monitored using RT-qPCR (Figure 3). Without exception, the six monitored genes show a clear up-regulation upon exogenous expression of miR-200a (Figure 3a). The reverse experiment, inhibition of some members of the miR-200 family (miR-200a,b,c) in the MDA-MB-231 cells using antagonists, resulted in the significant down-regulation of four out of five tested genes, (SCNN1A is not significantly down-regulated, ATAD4 is not expressed in normal conditions in this cell line)(Figure 3b). Those results clearly show that miR-200a is a core regulator of module 25, most probably with other members of the miR-200 family.

In this module, we observe again a clear positive correlation pattern between miR-200a and the module genes expression, suggesting an indirect regulatory circuit between miR-200a and module 25 genes (Figure 2c and Figure 3a and 3b). Recent experimental work showed that miR-200 family members directly target the transcription factors ZEB1 and ZEB2 [47,48,49]. These transcription factors are known as major transcriptional repressors of epithelial differentiation orchestrating epithelial mesenchymal transition (EMT) [50]. EMT is a process that drives epithelial cells from a polarized phenotype to a highly motile, non polarized mesenchymal phenotype and is known to occur in epithelial tumors giving rise to highly malignant cancer cells. The ZEB transcription factors have been functionally related to members of the miR-200 family via a double negative feedback loop, thus promoting EMT and cancer invasion [47,51,52]. We found conserved ZEB binding motifs for several module 25 genes (Data file S1), suggesting that ZEB factors could be the intermediate regulators between miR-200 and module 25 genes. To test this hypothesis, we down-regulated the ZEB1 transcription factor in MDA-MD-231 cells with a specific siRNA while monitoring the expression of eight module 25 genes (Figure 3c). All the genes show a strong up-regulation pattern, with the exception of the gene RAB25 (Figure 3c). Those results demonstrate that the ZEB1 factor is essential for the regulation of module 25 genes. Taken together, our experimental results (Figure 3) strongly suggest the existence of a regulatory chain between miR-200 and module 25 genes via the ZEB1 transcription factor (Figure 4). As both miR-200 and ZEB1 play important roles in EMT [47,48,49,51,52] we hypothesize that module 25 repression might contribute to the malignant EMT process in cancer cells.

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 2.** Modules 29 (a), 18 (b) and 25 (c) genes (MG) and assigned regulators (MR). Gene expression values are color coded, ranging from dark blue (low expression levels) to bright yellow (high expression levels). In each figure, columns represent a different sample. The color-coded bar at the bottom of the graph represents the tissue origin (see the legend), while the gray squares just below indicate whether the sample tissue was normal (light gray) or tumor (dark gray). The candidate regulators are ordered by decreasing score value (from top to bottom). The samples are</figcaption></figure>

<figure><img src="image_9.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 2.** Modules 29 (a), 18 (b) and 25 (c) genes (MG) and assigned regulators (MR). Gene expression values are color coded, ranging from dark blue (low expression levels) to bright yellow (high expression levels). In each figure, columns represent a different sample. The color-coded bar at the bottom of the graph represents the tissue origin (see the legend), while the gray squares just below indicate whether the sample tissue was normal (light gray) or tumor (dark gray). The candidate regulators are ordered by decreasing score value (from top to bottom). The samples are</figcaption></figure>

<figure><img src="image_9.png" /></figure>
grouped in leaves of homogeneous expression values, according to the hierarchical trees indicated on top of each figure. The orange boxes at the right of the figure indicate overrepresented GO categories (p≤0.05).
doi:10.1371/journal.pone.0010162.g002

Discussion

MiRNAs have emerged quite recently as a new and important
layer of regulation. Most of the studies so far have focused on their
identification and on the detection of their targets. Several
experimental studies have shown that at least some of them play
key roles in various developmental and cellular pathways.
Integrating miRNAs in regulatory networks is therefore of
fundamental importance and should ideally be done taking into
account the other types of regulatory molecules. So far, a few studies
have proposed a computational strategy to infer miRNA mediated
module networks [14,15,16]. These were mainly based on miRNA
target prediction, or on a combination of target prediction and
expression data. We have applied a robust and unbiased module
network inference algorithm to a cancer-related expression data set
of both mRNAs and miRNAs. In our approach, miRNAs were
considered as candidate regulators, together with other types of
regulators, like transcription factors and signal transducers. Even
after applying a stringent cutoff, several miRNAs were retained as
high scoring, statistically significant candidate regulators for various
modules. Through an in-depth analysis of three of those modules,
we showed that the assignment of specific miRNAs as regulators is
supported by various external sources and is functionally coherent.
Furthermore, we could show experimentally that a miRNA,
assigned as the best regulator, is indeed a key regulator for the
module genes expression. The number of miRNAs assigned in this
study (10) might seem rather modest, but this number has to be
evaluated with respect to the total number of miRNAs for which
expression was measured in the samples (124). The ratio assigned
per total number of miRNAs is equal to 8%, while the same ratio
value is 15% for the ensemble transcription factors plus signal
transducers (284/1841). The two ratio values are comparable and
therefore we can reasonably expect a higher number of miRNAs
assigned when an increased coverage of the miRNome expression
landscape will be available.

Nevertheless, just as with other similar methods, care has to be
taken for the interpretation of the inferred regulatory model. In
particular, correlation of gene expression might not always indicate
a direct interaction. Indeed, for the three modules we have
investigated in detail, we have found an indirect regulation pathway
between the regulator and the module genes. Furthermore, none of
those indirect regulator genes were assigned by the algorithm in the
regulation program or even clustered together with the module
genes. As we could show for module 29 and the SRF transcription
factor, the reason is because those indirect regulators expression
profile differ significantly from those of the module genes. Various
reasons can explain this divergence, for example the regulation
might happen at the post-transcriptional level, or might be the result
of post-translational modifications. Indirect regulators might of
course complicate the interpretation of the results but they are to be
expected, especially in higher eukaryotic organisms where regula-
tory networks are expected to be more complex [53].

Taken together, our results show that novel insights can be gained
from a robust module network analysis of miRNA and mRNA
expression data and support the view that at least some miRNAs
have key regulatory roles in important cellular processes. Our
approach has also the advantage of providing a direct view of post-
transcriptional modifications through the integration of miRNAs,
where mRNA expression alone might not be enough to reveal the
existence of regulatory interactions. All three modules for which we

did a detailed analysis in this study have each a coherent set of
genes, involved in the same process and function. Furthermore, by
connecting miRNAs to coherent modules, we believe that this
approach can help to elucidate miRNA function and could
efficiently drive experimental work towards the identification of
key regulatory components in various processes. With the rapid
proliferation of various techniques to measure with a high accuracy
the levels of expression for hundreds of miRNAs, and the
concomitant availability of mRNA expression data, it will be highly
appealing to apply computational strategies like the one we describe
here to expand our knowledge on global regulatory networks.

**Materials and Methods**

Expression data sets

We used a normalized cancer expression data set previously published [13]. We performed additional filtering steps to improve the quality of the input data set. Probesets with no known ensembl gene identifiers were discarded, as well as miRNA sequences that were not annotated as human miRNAs in the most recent miRBase release [3]. The final data matrix contained 11,833 genes and 124 miRNAs, for which expression was measured across 89 samples covering 11 different tumor classes.

Module network inference

We used the LeMoNe algorithm to infer the module network
[11,12,54]. In a first step, the algorithm is searching for a partition
of genes into clusters of co-expressed genes. In a second step, the
algorithm defines a regulatory program (a set of regulator genes)
for each cluster. To avoid local optima traps in the first step, the
algorithm uses a gibbs sampling approach for two-way clustering
of both genes and conditions [54]. For a given input expression
matrix, multiple clustering solutions are generated. For this study,
we generated 30 different cluster solutions from the initial dataset.
This ensemble of partially overlapping solutions is averaged to
produce a set of tight clusters, representing subsets of genes which
consistently cluster together in all solutions. The set of tight clusters
is extracted using a graph spectral method [54]. For the second
step, regulation programs are learned using a fuzzy decision tree
model. The two-way clustering of the first step has also generated
condition clusters (set of conditions having a similar mean and
standard deviation) for each set of co-expressed genes. The
condition clusters of a given module are first linked together in a
hierarchical decision tree. Each node in the tree is defining a split
between two sets of conditions (corresponding to low and high
expression levels). Regulators are assigned to each node of the tree
using a probabilistic score reflecting how well the expression levels
of the regulator match the genes expression levels defined by the
split value (for details about the mathematical model of the
algorithm, see [11]). Just as for the gene clusters, multiple solutions
are generated for the conditions clusters. Consequently, there are
multiple decision trees and multiple regulators assigned for each
node of each hierarchical tree. We adopt an ensemble approach
again by summing the strength with which a regulator participates
in each regulatory program for a given set of co-expressed genes. A
global score is calculated, reflecting the statistical confidence of the
regulator over all the nodes of all the hierarchical trees generated
for the set of co-expressed genes [11]. For this study, we assigned
up to 100 regulators for each node of each of the 100 hierarchical
trees defined for each module. It is worth noting that by using a

<figure><img src="image_10.png" /></figure>
grouped in leaves of homogeneous expression values, according to the hierarchical trees indicated on top of each figure. The orange boxes at the right of the figure indicate overrepresented GO categories (p≤0.05).
doi:10.1371/journal.pone.0010162.g002

Discussion

MiRNAs have emerged quite recently as a new and important
layer of regulation. Most of the studies so far have focused on their
identification and on the detection of their targets. Several
experimental studies have shown that at least some of them play
key roles in various developmental and cellular pathways.
Integrating miRNAs in regulatory networks is therefore of
fundamental importance and should ideally be done taking into
account the other types of regulatory molecules. So far, a few studies
have proposed a computational strategy to infer miRNA mediated
module networks [14,15,16]. These were mainly based on miRNA
target prediction, or on a combination of target prediction and
expression data. We have applied a robust and unbiased module
network inference algorithm to a cancer-related expression data set
of both mRNAs and miRNAs. In our approach, miRNAs were
considered as candidate regulators, together with other types of
regulators, like transcription factors and signal transducers. Even
after applying a stringent cutoff, several miRNAs were retained as
high scoring, statistically significant candidate regulators for various
modules. Through an in-depth analysis of three of those modules,
we showed that the assignment of specific miRNAs as regulators is
supported by various external sources and is functionally coherent.
Furthermore, we could show experimentally that a miRNA,
assigned as the best regulator, is indeed a key regulator for the
module genes expression. The number of miRNAs assigned in this
study (10) might seem rather modest, but this number has to be
evaluated with respect to the total number of miRNAs for which
expression was measured in the samples (124). The ratio assigned
per total number of miRNAs is equal to 8%, while the same ratio
value is 15% for the ensemble transcription factors plus signal
transducers (284/1841). The two ratio values are comparable and
therefore we can reasonably expect a higher number of miRNAs
assigned when an increased coverage of the miRNome expression
landscape will be available.

Nevertheless, just as with other similar methods, care has to be
taken for the interpretation of the inferred regulatory model. In
particular, correlation of gene expression might not always indicate
a direct interaction. Indeed, for the three modules we have
investigated in detail, we have found an indirect regulation pathway
between the regulator and the module genes. Furthermore, none of
those indirect regulator genes were assigned by the algorithm in the
regulation program or even clustered together with the module
genes. As we could show for module 29 and the SRF transcription
factor, the reason is because those indirect regulators expression
profile differ significantly from those of the module genes. Various
reasons can explain this divergence, for example the regulation
might happen at the post-transcriptional level, or might be the result
of post-translational modifications. Indirect regulators might of
course complicate the interpretation of the results but they are to be
expected, especially in higher eukaryotic organisms where regula-
tory networks are expected to be more complex [53].

Taken together, our results show that novel insights can be gained
from a robust module network analysis of miRNA and mRNA
expression data and support the view that at least some miRNAs
have key regulatory roles in important cellular processes. Our
approach has also the advantage of providing a direct view of post-
transcriptional modifications through the integration of miRNAs,
where mRNA expression alone might not be enough to reveal the
existence of regulatory interactions. All three modules for which we

did a detailed analysis in this study have each a coherent set of
genes, involved in the same process and function. Furthermore, by
connecting miRNAs to coherent modules, we believe that this
approach can help to elucidate miRNA function and could
efficiently drive experimental work towards the identification of
key regulatory components in various processes. With the rapid
proliferation of various techniques to measure with a high accuracy
the levels of expression for hundreds of miRNAs, and the
concomitant availability of mRNA expression data, it will be highly
appealing to apply computational strategies like the one we describe
here to expand our knowledge on global regulatory networks.

**Materials and Methods**

Expression data sets

We used a normalized cancer expression data set previously published [13]. We performed additional filtering steps to improve the quality of the input data set. Probesets with no known ensembl gene identifiers were discarded, as well as miRNA sequences that were not annotated as human miRNAs in the most recent miRBase release [3]. The final data matrix contained 11,833 genes and 124 miRNAs, for which expression was measured across 89 samples covering 11 different tumor classes.

Module network inference

We used the LeMoNe algorithm to infer the module network
[11,12,54]. In a first step, the algorithm is searching for a partition
of genes into clusters of co-expressed genes. In a second step, the
algorithm defines a regulatory program (a set of regulator genes)
for each cluster. To avoid local optima traps in the first step, the
algorithm uses a gibbs sampling approach for two-way clustering
of both genes and conditions [54]. For a given input expression
matrix, multiple clustering solutions are generated. For this study,
we generated 30 different cluster solutions from the initial dataset.
This ensemble of partially overlapping solutions is averaged to
produce a set of tight clusters, representing subsets of genes which
consistently cluster together in all solutions. The set of tight clusters
is extracted using a graph spectral method [54]. For the second
step, regulation programs are learned using a fuzzy decision tree
model. The two-way clustering of the first step has also generated
condition clusters (set of conditions having a similar mean and
standard deviation) for each set of co-expressed genes. The
condition clusters of a given module are first linked together in a
hierarchical decision tree. Each node in the tree is defining a split
between two sets of conditions (corresponding to low and high
expression levels). Regulators are assigned to each node of the tree
using a probabilistic score reflecting how well the expression levels
of the regulator match the genes expression levels defined by the
split value (for details about the mathematical model of the
algorithm, see [11]). Just as for the gene clusters, multiple solutions
are generated for the conditions clusters. Consequently, there are
multiple decision trees and multiple regulators assigned for each
node of each hierarchical tree. We adopt an ensemble approach
again by summing the strength with which a regulator participates
in each regulatory program for a given set of co-expressed genes. A
global score is calculated, reflecting the statistical confidence of the
regulator over all the nodes of all the hierarchical trees generated
for the set of co-expressed genes [11]. For this study, we assigned
up to 100 regulators for each node of each of the 100 hierarchical
trees defined for each module. It is worth noting that by using a

<figure><img src="image_10.png" /></figure>
a

<figure><img src="image_11.png" /></figure>

b

<figure><img src="image_12.png" /></figure>

c

<figure><img src="image_13.png" /><figcaption>**Figure 3. Validation of miR-200a as a regulator of module 25 genes expression.** (a) Real time quantitative PCR (RT qPCR) analysis of the expression of module 25 genes *RAB25*, *IRF6*, *SCNN1A*, *PRSS8*, *ATAD4* and *TACSTD1* and upon over-expression of miR-200a in MDA-MD-231 cells (mean ± standard deviation). The y-axis represents the relative mRNA expression value. miR-1 was used as the control (Ctrl), as it is not known to target any of the monitored genes (b) RT qPCR analysis of the relative expression for the genes *RAB25*, *IRF6*, *SCNN1A*, *PRSS8* and *TACSTD1* in MDA-MB-231 cells infiltrated with miR-200a,b,c antagonists. The y-axis represents the relative mRNA expression value (mean ± standard deviation). miR-1 was used as the control (Ctrl) (c) RT qPCR analysis of the relative expression levels of module 25 genes *RAB25*, *IRF6*, *SCNN1A*, *PRSS8*, *ATAD4*, *TACSTD1*, *TMEM63A* and *FXYD3* in MDA-MB-231 cells where *ZEB1* is knocked down. The first barplot (black) shows effective repression of *ZEB1* levels upon transfection with the *ZEB1* specific siRNA. Par = parental cell culture, Mock = mock transfection, si-ZEB1 = transfection with the ZEB1 specific siRNA.
doi:10.1371/journal.pone.0010162.g003</figcaption></figure>

score that only takes into account the differential expression of a
regulator across the different condition clusters, we can simulta-
neously evaluate and compare mRNA and miRNA candidate
regulators. In the end, the set of regulators assigned to each cluster
of co-expressed genes can be ranked according to their global
probabilistic score and a cutoff level can be defined, keeping only
very high-scoring regulators. In order to evaluate the statistical

significance of the assigned regulators, a second set of randomly
assigned regulators is generated along the set of “true” regulators
(Figure 1). The complete list of modules together with their high-
scoring regulators for this study is available in the table S1. The
LeMoNe software package can be downloaded from our website,
is open-source and free of charge for academics (http://
bioinformatics.psb.ugent.be/software/details/lemone).

<figure><img src="image_14.png" /></figure>
a

<figure><img src="image_11.png" /></figure>

b

<figure><img src="image_12.png" /></figure>

c

<figure><img src="image_13.png" /><figcaption>**Figure 3. Validation of miR-200a as a regulator of module 25 genes expression.** (a) Real time quantitative PCR (RT qPCR) analysis of the expression of module 25 genes *RAB25*, *IRF6*, *SCNN1A*, *PRSS8*, *ATAD4* and *TACSTD1* and upon over-expression of miR-200a in MDA-MD-231 cells (mean ± standard deviation). The y-axis represents the relative mRNA expression value. miR-1 was used as the control (Ctrl), as it is not known to target any of the monitored genes (b) RT qPCR analysis of the relative expression for the genes *RAB25*, *IRF6*, *SCNN1A*, *PRSS8* and *TACSTD1* in MDA-MB-231 cells infiltrated with miR-200a,b,c antagonists. The y-axis represents the relative mRNA expression value (mean ± standard deviation). miR-1 was used as the control (Ctrl) (c) RT qPCR analysis of the relative expression levels of module 25 genes *RAB25*, *IRF6*, *SCNN1A*, *PRSS8*, *ATAD4*, *TACSTD1*, *TMEM63A* and *FXYD3* in MDA-MB-231 cells where *ZEB1* is knocked down. The first barplot (black) shows effective repression of *ZEB1* levels upon transfection with the *ZEB1* specific siRNA. Par = parental cell culture, Mock = mock transfection, si-ZEB1 = transfection with the ZEB1 specific siRNA.
doi:10.1371/journal.pone.0010162.g003</figcaption></figure>

score that only takes into account the differential expression of a
regulator across the different condition clusters, we can simulta-
neously evaluate and compare mRNA and miRNA candidate
regulators. In the end, the set of regulators assigned to each cluster
of co-expressed genes can be ranked according to their global
probabilistic score and a cutoff level can be defined, keeping only
very high-scoring regulators. In order to evaluate the statistical

significance of the assigned regulators, a second set of randomly
assigned regulators is generated along the set of “true” regulators
(Figure 1). The complete list of modules together with their high-
scoring regulators for this study is available in the table S1. The
LeMoNe software package can be downloaded from our website,
is open-source and free of charge for academics (http://
bioinformatics.psb.ugent.be/software/details/lemone).

<figure><img src="image_14.png" /></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 4.** Module 25 hypothetical regulation model. MiR-200 genes repress ZEB factors, which in turn repress the expression of module 25 genes. The light yellow indicate genes assigned as regulators, the light green indicates module genes while the light orange indicates genes not assigned as regulators, but supported by literature (indirect regulation). This regulatory model support the positive correlation of the expression patterns between mir-200a and module 25 genes.</figcaption></figure>

doi:10.1371/journal.pone.0010162.g004

### Gene Ontology over-represented categories

For each module, we calculated GO enrichment using the BiNGO tool [17]. The complete list of GO categories enrichment for all the modules is available in the table S2.

### Transcription factor binding motifs

We used the ConTra [55] software tool to look for conserved TCF12 and ZEB motif binding sites in the promoter regions of module 18 and 25 genes. A multiple alignment of nine eutherian mammal species (Bos taurus, Canis familiaris, Equus caballus, Pan troglodytes, Pongo pygmaeus, Macacca mulatta, Mus musculus and Rattus norvegicus) and a specific position weight matrix were used to determine the conservation of the motif across all species.

### Cell culture

Human cancer cell lines were originally obtained from ATCC. MDA-MB-231 cells were maintained in Leibovitz's L-15 medium supplemented with 10% fetal calf serum (FCS), 100 µg/ml penicillin, 100 µg/ml streptomycin and 0.03% L-Glutamine. These cells were grown at 37°C without CO₂ supply.

### miRNA repression and overexpression assays

MDA-MB-231 cells were seeded at 200.000 cells per well in 6-well plates in complete medium without antibiotics one day prior to transfection. The miRNA precursors and inhibitors as well as the positive and the negative control miRNAs were transfected at a final concentration of 25 nM using DharmaFECT 1 transfection reagent (ThermoSCIENTIFIC- Dharmacon) according to the manufacturer's instructions with the modification of using 8 µl of reagent instead of 6. The medium was refreshed after 18-24 hrs for the MDA-MB-231 cells and total RNA was collected 48 hrs post transfection. The negative control is Pre-miR™ miRNA Precursor-Negative Control #1, which does not target any known mRNA

within the human or mouse transcriptome. The positive control is miR-1 Pre-miR miRNA precursor which has been shown to effectively downregulate the expression of twinfilin-1, also known as PTK9, at the mRNA level [56]. Validation of the downregulation of PTK9 was performed using a TaqMan® Gene Expression Assay (Assay ID: Hs00702289_s1). The control miRNAs and the qRT-PCR assay for human PTK9 were provided in the Pre-miR™ miRNA Starter Kit (Ambion Cat #AM1540).

### ZEB1 repression assay

MDA-MB-231 cells were seeded at 200.000 cells per well in 6-well plates in complete medium without antibiotics one day prior to transfection. ZEB1 siGENOME-SMARTpool was used for the downregulation of ZEB1 (ThermoSCIENTIFIC- Dharmacon, M-006564-02-0010), which consists of four SMART-selection designed siRNAs targeting one gene. The siZEB1 was dissolved in 1× siRNA buffer (ThermoSCIENTIFIC- Dharmacon, B-002000-UB) at a final concentration of 20 µM and was transfected at a final concentration of 25 nM using DharmaFECT 1 transfection reagent (ThermoSCIENTIFIC- Dharmacon) according to the manufacturer's instructions with the modification of using 8 µl of reagent instead of 6. As a negative control 1× siRNA buffer was used (MOCK transfection).

### Quantitative reverse transcription PCR (qRT-PCR)

Total RNA was extracted using Trizol (Invitrogen) according to the manufacturer's instructions with one modification; absolute ethanol was used in place of isopropanol. For the qPCR analysis cDNA synthesis was performed on 1 µg of total RNA using the iScript synthesis kit (BIO-RAD). The qRT-PCR for every gene was performed on 20 ng of cDNA in triplicate using the SYBRGreen I Master (Roche) or Probes Master (Roche) on a LightCycler®480 Real-time PCR System (Roche). The expression levels were determined using comparative quantification to the negative control and all quantification data were normalized against 2 reference genes, HMBS and TBP. The sequences of the RT-qPCR primers that were used are given in the data file S1.

### Supporting Information

<caption>**Table S1** Complete list of module genes and regulators.</caption>

Found at: doi:10.1371/journal.pone.0010162.s001 (0.28 MB XLS)

<caption>**Table S2** Gene Ontology (GO) categories enrichment for each module.</caption>

Found at: doi:10.1371/journal.pone.0010162.s002 (0.14 MB XLS)

**Data file S1** Gene expression profile of MYH11, CNN1, ACTG2 and MYLK compared to SRF; TCF12 binding motifs for module 18 genes; ZEB binding motifs for module 25 genes; PCR primers.

Found at: doi:10.1371/journal.pone.0010162.s003 (0.49 MB DOC)

### Acknowledgments

We thank Chris Marine for fruitful discussions and suggestions. We also thank the anonymous reviewers for valuable comments.

### Author Contributions

Conceived and designed the experiments: EB MT GB. Performed the experiments: MT. Analyzed the data: EB AJ TM. Wrote the paper: EB KM GB YVdP.

<figure><img src="image_16.png" /></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 4.** Module 25 hypothetical regulation model. MiR-200 genes repress ZEB factors, which in turn repress the expression of module 25 genes. The light yellow indicate genes assigned as regulators, the light green indicates module genes while the light orange indicates genes not assigned as regulators, but supported by literature (indirect regulation). This regulatory model support the positive correlation of the expression patterns between mir-200a and module 25 genes.</figcaption></figure>

doi:10.1371/journal.pone.0010162.g004

### Gene Ontology over-represented categories

For each module, we calculated GO enrichment using the BiNGO tool [17]. The complete list of GO categories enrichment for all the modules is available in the table S2.

### Transcription factor binding motifs

We used the ConTra [55] software tool to look for conserved TCF12 and ZEB motif binding sites in the promoter regions of module 18 and 25 genes. A multiple alignment of nine eutherian mammal species (Bos taurus, Canis familiaris, Equus caballus, Pan troglodytes, Pongo pygmaeus, Macacca mulatta, Mus musculus and Rattus norvegicus) and a specific position weight matrix were used to determine the conservation of the motif across all species.

### Cell culture

Human cancer cell lines were originally obtained from ATCC. MDA-MB-231 cells were maintained in Leibovitz's L-15 medium supplemented with 10% fetal calf serum (FCS), 100 µg/ml penicillin, 100 µg/ml streptomycin and 0.03% L-Glutamine. These cells were grown at 37°C without CO₂ supply.

### miRNA repression and overexpression assays

MDA-MB-231 cells were seeded at 200.000 cells per well in 6-well plates in complete medium without antibiotics one day prior to transfection. The miRNA precursors and inhibitors as well as the positive and the negative control miRNAs were transfected at a final concentration of 25 nM using DharmaFECT 1 transfection reagent (ThermoSCIENTIFIC- Dharmacon) according to the manufacturer's instructions with the modification of using 8 µl of reagent instead of 6. The medium was refreshed after 18-24 hrs for the MDA-MB-231 cells and total RNA was collected 48 hrs post transfection. The negative control is Pre-miR™ miRNA Precursor-Negative Control #1, which does not target any known mRNA

within the human or mouse transcriptome. The positive control is miR-1 Pre-miR miRNA precursor which has been shown to effectively downregulate the expression of twinfilin-1, also known as PTK9, at the mRNA level [56]. Validation of the downregulation of PTK9 was performed using a TaqMan® Gene Expression Assay (Assay ID: Hs00702289_s1). The control miRNAs and the qRT-PCR assay for human PTK9 were provided in the Pre-miR™ miRNA Starter Kit (Ambion Cat #AM1540).

### ZEB1 repression assay

MDA-MB-231 cells were seeded at 200.000 cells per well in 6-well plates in complete medium without antibiotics one day prior to transfection. ZEB1 siGENOME-SMARTpool was used for the downregulation of ZEB1 (ThermoSCIENTIFIC- Dharmacon, M-006564-02-0010), which consists of four SMART-selection designed siRNAs targeting one gene. The siZEB1 was dissolved in 1× siRNA buffer (ThermoSCIENTIFIC- Dharmacon, B-002000-UB) at a final concentration of 20 µM and was transfected at a final concentration of 25 nM using DharmaFECT 1 transfection reagent (ThermoSCIENTIFIC- Dharmacon) according to the manufacturer's instructions with the modification of using 8 µl of reagent instead of 6. As a negative control 1× siRNA buffer was used (MOCK transfection).

### Quantitative reverse transcription PCR (qRT-PCR)

Total RNA was extracted using Trizol (Invitrogen) according to the manufacturer's instructions with one modification; absolute ethanol was used in place of isopropanol. For the qPCR analysis cDNA synthesis was performed on 1 µg of total RNA using the iScript synthesis kit (BIO-RAD). The qRT-PCR for every gene was performed on 20 ng of cDNA in triplicate using the SYBRGreen I Master (Roche) or Probes Master (Roche) on a LightCycler®480 Real-time PCR System (Roche). The expression levels were determined using comparative quantification to the negative control and all quantification data were normalized against 2 reference genes, HMBS and TBP. The sequences of the RT-qPCR primers that were used are given in the data file S1.

### Supporting Information

<caption>**Table S1** Complete list of module genes and regulators.</caption>

Found at: doi:10.1371/journal.pone.0010162.s001 (0.28 MB XLS)

<caption>**Table S2** Gene Ontology (GO) categories enrichment for each module.</caption>

Found at: doi:10.1371/journal.pone.0010162.s002 (0.14 MB XLS)

**Data file S1** Gene expression profile of MYH11, CNN1, ACTG2 and MYLK compared to SRF; TCF12 binding motifs for module 18 genes; ZEB binding motifs for module 25 genes; PCR primers.

Found at: doi:10.1371/journal.pone.0010162.s003 (0.49 MB DOC)

### Acknowledgments

We thank Chris Marine for fruitful discussions and suggestions. We also thank the anonymous reviewers for valuable comments.

### Author Contributions

Conceived and designed the experiments: EB MT GB. Performed the experiments: MT. Analyzed the data: EB AJ TM. Wrote the paper: EB KM GB YVdP.

<figure><img src="image_16.png" /></figure>
References

1. Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116: 281-297.

2. Medina PP, Slack FJ (2008) microRNAs and cancer: an overview. Cell Cycle 7: 2485-2492.

3. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ (2008) miRBase: tools for microRNA genomics. Nucleic Acids Res 36: D154-158.

4. Berezikov E, Guryev V, van de Belt J, Wienholds E, Plasterk RH, et al. (2005) Phylogenetic shadowing and computational identification of human microRNA genes. Cell 120: 21-24.

5. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, et al. (2007) MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 27: 91-105.

6. Friedman RC, Farh KK, Burge CB, Bartel D (2008) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res.

7. Walhout AJ (2006) Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping. Genome Res 16: 1445-1454.

8. Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, et al. (2005) Reverse engineering of regulatory networks in human B cells. Nat Genet 37: 382-390.

9. Faith JJ, Hayete B, Thaden JT, Mogno I, Wierzbowski J, et al. (2007) Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS Biol 5: e8.

10. Segal E, Shapira M, Regev A, Pe'er D, Botstein D, et al. (2003) Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet 34: 166-176.

11. Joshi A, De Smet R, Marchal K, Van de Peer Y, Michael T (2009) Module networks revisited: computational assessment and prioritization of model predictions. Bioinformatics 25: 490-496.

12. Michael T, De Smet R, Joshi A, Van de Peer Y, Marchal K (2009) Comparative analysis of module-based versus direct methods for reverse-engineering transcriptional regulatory networks. BMC Syst Biol 3: 49.

13. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, et al. (2005) MicroRNA expression profiles classify human cancers. Nature 435: 834-838.

14. Joung JG, Hwang KB, Nam JW, Kim SJ, Zhang BT (2007) Discovery of microRNA-mRNA modules via population-based probabilistic learning. Bioinformatics 23: 1141-1147.

15. Yoon S, De Micheli G (2005) Prediction of regulatory modules comprising microRNAs and target genes. Bioinformatics 21 Suppl 2: ii93-100.

16. Liu B, Li J, Tsykin A, Liu L, Gaur A, et al. (2009) Exploring complex miRNA-mRNA interactions with Bayesian networks by splitting-averaging strategy. BMC Bioinformatics 10: 408.

17. Maere S, Heymans K, Kuiper M (2005) BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21: 3448-3449.

18. Kamm KE, Stull JT (2001) Dedicated myosin light chain kinases with diverse cellular functions. J Biol Chem 276: 4527-4530.

19. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, et al. (2006) The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet 38: 228-233.

20. Care A, Catalucci D, Felicetti F, Bonci D, Addario A, et al. (2007) MicroRNA-133 controls cardiac hypertrophy. Nat Med 13: 613-618.

21. Miano JM, Long X, Fujiwara K (2007) Serum response factor: master regulator of the actin cytoskeleton and contractile apparatus. Am J Physiol Cell Physiol 292: G70-81.

22. Stritt C, Stern S, Harting K, Manke T, Sinske D, et al. (2009) Paracrine control of oligodendrocyte differentiation by SRF-directed neuronal gene expression. Nat Neurosci 12: 418-427.

23. Rivera VM, Sheng M, Greenberg ME (1990) The inner core of the serum response element mediates both the rapid induction and subsequent repression of c-fos transcription following serum stimulation. Genes Dev 4: 255-268.

24. Shaw PE, Frasch S, Nordheim A (1989) Repression of c-fos transcription is mediated through p67SRF bound to the SRE. EMBO J 8: 2567-2574.

25. Cordes KR, Sheehy NT, White MP, Berry EC, Morton SU, et al. (2009) miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature 460: 705-710.

26. Wang Z, Wang DZ, Pipes GC, Olson EN (2003) Myocardin is a master regulator of smooth muscle gene expression. Proc Natl Acad Sci U S A 100: 7129-7134.

27. Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D (1998) GeneCards: a novel functional genomics compendium with automated data mining and query reformulation support. Bioinformatics 14: 656-664.

28. McCarthy KM, McDevit D, Andreucci A, Reeves R, Nikolajczyk BS (2003) HMGA1 co-activates transcription in B cells through indirect association with DNA. J Biol Chem 278: 42106-42114.

29. Turesson C, Matteson EL (2006) Genetics of rheumatoid arthritis. Mayo Clin Proc 81: 94-101.

30. Maghazachi AA, Al-Aoukaty A, Schall TJ (1996) CC chemokines induce the generation of killer cells from CD56+ cells. Eur J Immunol 26: 315-319.

31. O'Bryan JP, Fridell YW, Koski R, Varnum B, Liu ET (1995) The transforming receptor tyrosine kinase, Axl, is post-translationally regulated by proteolytic cleavage. J Biol Chem 270: 551-557.

32. Kurth I, Willimann K, Scherli P, Hunziker T, Clark-Lewis I, et al. (2001) Monocyte selectivity and tissue localization suggests a role for breast and kidney-expressed chemokine (BRAK) in macrophage development. J Exp Med 194: 855-861.

33. Frederick MJ, Henderson Y, Xu X, Deavers MT, Sahin AA, et al. (2000) In vivo expression of the novel CXC chemokine BRAK in normal and cancerous human tissue. Am J Pathol 156: 1937-1950.

34. Chen CZ, Li L, Lodish HF, Bartel DP (2004) MicroRNAs modulate hematopoietic lineage differentiation. Science 303: 83-86.

35. Liao R, Sun J, Zhang L, Lou G, Chen M, et al. (2008) MicroRNAs play a role in the development of human hematopoietic stem cells. J Cell Biochem 104: 805-817.

36. Zhuang Y, Cheng P, Weintraub H (1996) B-lymphocyte development is regulated by the combined dosage of three basic helix-loop-helix genes, E2A, E2-2, and HEB. Mol Cell Biol 16: 2898-2905.

37. McDonald EJ, Snyder PM, McCray PB, Jr., Welsh MJ (1994) Cloning, expression, and tissue distribution of a human amiloride-sensitive Na+ channel. Am J Physiol 266: L728-734.

38. Chen M, Fu YY, Lin CY, Chen LM, Chai KX (2007) Prostasin induces protease-dependent and independent molecular changes in the human prostate carcinoma cell line PC-3. Biochim Biophys Acta 1773: 1133-1140.

39. Crambert G, Li C, Claeys D, Geering K (2005) FXYD3 (Mat-8), a new regulator of Na,K-ATPase. Mol Biol Cell 16: 2363-2371.

40. Gires O, Eskofier S, Lang S, Zeidler R, Munz M (2003) Cloning and characterisation of a 1.1 kb fragment of the carcinoma-associated epithelial cell adhesion molecule promoter. Anticancer Res 23: 3255-3261.

41. Cheng KW, Lahad JP, Kuo WL, Lapuk A, Yamada K, et al. (2004) The RAB25 small GTPase determines aggressiveness of ovarian and breast cancers. Nat Med 10: 1251-1256.

42. Beckler A, Moskaluk CA, Zaika A, Hampton GM, Powell SM, et al. (2003) Overexpression of the 32-kilodalton dopamine and cyclic adenosine 3',5'-monophosphate-regulated phosphoprotein in common adenocarcinomas. Cancer 98: 1547-1551.

43. Filardo EJ (2002) Epidermal growth factor receptor (EGFR) transactivation by estrogen via the G-protein-coupled receptor, GPR30: a novel signaling pathway with potential significance for breast cancer. J Steroid Biochem Mol Biol 80: 231-238.

44. Grzmil M, Voigt S, Thelen P, Hemmerlein B, Helmke K, et al. (2004) Up-regulated expression of the MAT-8 gene in prostate cancer and its siRNA-mediated inhibition of expression induces a decrease in proliferation of human prostate carcinoma cells. Int J Oncol 24: 97-105.

45. Baskerville S, Bartel DP (2005) Microarray profiling of microRNAs reveals frequent coexpression with neighboring mRNAs and host genes. Rna 11: 241-247.

46. Thomson JM, Parker J, Perou CM, Hammond SM (2004) A custom microarray platform for analysis of microRNA gene expression. Nat Methods 1: 47-53.

47. Burk U, Schubert J, Wellner U, Schmalhofer O, Vincan E, et al. (2008) A reciprocal repression between ZEB1 and members of the miR-200 family promotes EMT and invasion in cancer cells. EMBO Rep 9: 582-589.

48. Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, et al. (2008) The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol 10: 593-601.

49. Park SM, Gaur AB, Lengyel E, Peter ME (2008) The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. Genes Dev 22: 894-907.

50. Vandewalle C, Van Roy F, Berx G (2008) The role of the ZEB family of transcription factors in development and disease. Cell Mol Life Sci.

51. Adam L, Zhong M, Choi W, Qi W, Nicoloso M, et al. (2009) miR-200 expression regulates epithelial-to-mesenchymal transition in bladder cancer cells and reverses resistance to epidermal growth factor receptor therapy. Clin Cancer Res 15: 5060-5072.

52. Kong D, Li Y, Wang Z, Banerjee S, Ahmad A, et al. (2009) miR-200 regulates PDGF-D-mediated epithelial-mesenchymal transition, adhesion, and invasion of prostate cancer cells. Stem Cells 27: 1712-1721.

53. Herrgard MJ, Covert MW, Palsson BO (2003) Reconciling gene expression data with known genome-scale regulatory network structures. Genome Res 13: 2423-2434.

54. Joshi A, Van de Peer Y, Michael T (2008) Analysis of a Gibbs sampler method for model-based clustering of gene expression data. Bioinformatics 24: 176-183.

55. Hooghe B, Hulpiau P, van Roy F, De Bleser P (2008) ConTra: a promoter alignment analysis tool for identification of transcription factor binding sites across species. Nucleic Acids Res 36: W128-132.

56. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, et al. (2005) Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature 433: 769-773.

<figure><img src="image_17.png" /></figure>
References

1. Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116: 281-297.

2. Medina PP, Slack FJ (2008) microRNAs and cancer: an overview. Cell Cycle 7: 2485-2492.

3. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ (2008) miRBase: tools for microRNA genomics. Nucleic Acids Res 36: D154–158.

4. Berezikov E, Guryev V, van de Belt J, Wienholds E, Plasterk RH, et al. (2005) Phylogenetic shadowing and computational identification of human microRNA genes. Cell 120: 21–24.

5. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, et al. (2007) MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 27: 91–105.

6. Friedman RC, Farh KK, Burge CB, Bartel D (2008) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res.

7. Walhout AJ (2006) Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping. Genome Res 16: 1445-1454.

8. Basso K, Margolin AA, Stolovitzky G, Klein U, Dalla-Favera R, et al. (2005) Reverse engineering of regulatory networks in human B cells. Nat Genet 37: 382–390.

9. Faith JJ, Hayete B, Thaden JT, Mogno I, Wierzbowski J, et al. (2007) Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS Biol 5: e8.

10. Segal E, Shapira M, Regev A, Pe'er D, Botstein D, et al. (2003) Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet 34: 166–176.

11. Joshi A, De Smet R, Marchal K, Van de Peer Y, Michael T (2009) Module networks revisited: computational assessment and prioritization of model predictions. Bioinformatics 25: 490–496.

12. Michael T, De Smet R, Joshi A, Van de Peer Y, Marchal K (2009) Comparative analysis of module-based versus direct methods for reverse-engineering transcriptional regulatory networks. BMC Syst Biol 3: 49.

13. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, et al. (2005) MicroRNA expression profiles classify human cancers. Nature 435: 834–838.

14. Joung JG, Hwang KB, Nam JW, Kim SJ, Zhang BT (2007) Discovery of microRNA-mRNA modules via population-based probabilistic learning. Bioinformatics 23: 1141–1147.

15. Yoon S, De Micheli G (2005) Prediction of regulatory modules comprising microRNAs and target genes. Bioinformatics 21 Suppl 2: ii93–100.

16. Liu B, Li J, Tsykin A, Liu L, Gaur A, et al. (2009) Exploring complex miRNA-mRNA interactions with Bayesian networks by splitting-averaging strategy. BMC Bioinformatics 10: 408.

17. Maere S, Heymans K, Kuiper M (2005) BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21: 3448–3449.

18. Kamm KE, Stull JT (2001) Dedicated myosin light chain kinases with diverse cellular functions. J Biol Chem 276: 4527–4530.

19. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, et al. (2006) The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet 38: 228–233.

20. Care A, Catalucci D, Felicetti F, Bonci D, Addario A, et al. (2007) MicroRNA-133 controls cardiac hypertrophy. Nat Med 13: 613–618.

21. Miano JM, Long X, Fujiwara K (2007) Serum response factor: master regulator of the actin cytoskeleton and contractile apparatus. Am J Physiol Cell Physiol 292: G70–81.

22. Stritt C, Stern S, Harting K, Manke T, Sinske D, et al. (2009) Paracrine control of oligodendrocyte differentiation by SRF-directed neuronal gene expression. Nat Neurosci 12: 418–427.

23. Rivera VM, Sheng M, Greenberg ME (1990) The inner core of the serum response element mediates both the rapid induction and subsequent repression of c-fos transcription following serum stimulation. Genes Dev 4: 255–268.

24. Shaw PE, Frasch S, Nordheim A (1989) Repression of c-fos transcription is mediated through p67SRF bound to the SRE. EMBO J 8: 2567–2574.

25. Cordes KR, Sheehy NT, White MP, Berry EC, Morton SU, et al. (2009) miR-145 and miR-143 regulate smooth muscle cell fate and plasticity. Nature 460: 705–710.

26. Wang Z, Wang DZ, Pipes GC, Olson EN (2003) Myocardin is a master regulator of smooth muscle gene expression. Proc Natl Acad Sci U S A 100: 7129–7134.

27. Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D (1998) GeneCards: a novel functional genomics compendium with automated data mining and query reformulation support. Bioinformatics 14: 656–664.

28. McCarthy KM, McDevit D, Andreucci A, Reeves R, Nikolajczyk BS (2003) HMGA1 co-activates transcription in B cells through indirect association with DNA. J Biol Chem 278: 42106–42114.

29. Turesson C, Matteson EL (2006) Genetics of rheumatoid arthritis. Mayo Clin Proc 81: 94–101.

30. Maghazachi AA, Al-Aoukaty A, Schall TJ (1996) CC chemokines induce the generation of killer cells from CD56+ cells. Eur J Immunol 26: 315–319.

31. O'Bryan JP, Fridell YW, Koski R, Varnum B, Liu ET (1995) The transforming receptor tyrosine kinase, Axl, is post-translationally regulated by proteolytic cleavage. J Biol Chem 270: 551–557.

32. Kurth I, Willimann K, Scherli P, Hunziker T, Clark-Lewis I, et al. (2001) Monocyte selectivity and tissue localization suggests a role for breast and kidney-expressed chemokine (BRAK) in macrophage development. J Exp Med 194: 855–861.

33. Frederick MJ, Henderson Y, Xu X, Deavers MT, Sahin AA, et al. (2000) In vivo expression of the novel CXC chemokine BRAK in normal and cancerous human tissue. Am J Pathol 156: 1937–1950.

34. Chen CZ, Li L, Lodish HF, Bartel DP (2004) MicroRNAs modulate hematopoietic lineage differentiation. Science 303: 83–86.

35. Liao R, Sun J, Zhang L, Lou G, Chen M, et al. (2008) MicroRNAs play a role in the development of human hematopoietic stem cells. J Cell Biochem 104: 805–817.

36. Zhuang Y, Cheng P, Weintraub H (1996) B-lymphocyte development is regulated by the combined dosage of three basic helix-loop-helix genes, E2A, E2-2, and HEB. Mol Cell Biol 16: 2898–2905.

37. McDonald EJ, Snyder PM, McCray PB, Jr., Welsh MJ (1994) Cloning, expression, and tissue distribution of a human amiloride-sensitive Na+ channel. Am J Physiol 266: L728–734.

38. Chen M, Fu YY, Lin CY, Chen LM, Chai KX (2007) Prostasin induces protease-dependent and independent molecular changes in the human prostate carcinoma cell line PC-3. Biochim Biophys Acta 1773: 1133–1140.

39. Crambert G, Li C, Claeys D, Geering K (2005) FXYD3 (Mat-8), a new regulator of Na,K-ATPase. Mol Biol Cell 16: 2363–2371.

40. Gires O, Eskofier S, Lang S, Zeidler R, Munz M (2003) Cloning and characterisation of a 1.1 kb fragment of the carcinoma-associated epithelial cell adhesion molecule promoter. Anticancer Res 23: 3255–3261.

41. Cheng KW, Lahad JP, Kuo WL, Lapuk A, Yamada K, et al. (2004) The RAB25 small GTPase determines aggressiveness of ovarian and breast cancers. Nat Med 10: 1251–1256.

42. Beckler A, Moskaluk CA, Zaika A, Hampton GM, Powell SM, et al. (2003) Overexpression of the 32-kilodalton dopamine and cyclic adenosine 3',5'-monophosphate-regulated phosphoprotein in common adenocarcinomas. Cancer 98: 1547–1551.

43. Filardo EJ (2002) Epidermal growth factor receptor (EGFR) transactivation by estrogen via the G-protein-coupled receptor, GPR30: a novel signaling pathway with potential significance for breast cancer. J Steroid Biochem Mol Biol 80: 231–238.

44. Grzmil M, Voigt S, Thelen P, Hemmerlein B, Helmke K, et al. (2004) Up-regulated expression of the MAT-8 gene in prostate cancer and its siRNA-mediated inhibition of expression induces a decrease in proliferation of human prostate carcinoma cells. Int J Oncol 24: 97–105.

45. Baskerville S, Bartel DP (2005) Microarray profiling of microRNAs reveals frequent coexpression with neighboring mRNAs and host genes. Rna 11: 241–247.

46. Thomson JM, Parker J, Perou CM, Hammond SM (2004) A custom microarray platform for analysis of microRNA gene expression. Nat Methods 1: 47–53.

47. Burk U, Schubert J, Wellner U, Schmalhofer O, Vincan E, et al. (2008) A reciprocal repression between ZEB1 and members of the miR-200 family promotes EMT and invasion in cancer cells. EMBO Rep 9: 582–589.

48. Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, et al. (2008) The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol 10: 593–601.

49. Park SM, Gaur AB, Lengyel E, Peter ME (2008) The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. Genes Dev 22: 894–907.

50. Vandewalle C, Van Roy F, Berx G (2008) The role of the ZEB family of transcription factors in development and disease. Cell Mol Life Sci.

51. Adam L, Zhong M, Choi W, Qi W, Nicoloso M, et al. (2009) miR-200 expression regulates epithelial-to-mesenchymal transition in bladder cancer cells and reverses resistance to epidermal growth factor receptor therapy. Clin Cancer Res 15: 5060–5072.

52. Kong D, Li Y, Wang Z, Banerjee S, Ahmad A, et al. (2009) miR-200 regulates PDGF-D-mediated epithelial-mesenchymal transition, adhesion, and invasion of prostate cancer cells. Stem Cells 27: 1712–1721.

53. Herrgard MJ, Covert MW, Palsson BO (2003) Reconciling gene expression data with known genome-scale regulatory network structures. Genome Res 13: 2423–2434.

54. Joshi A, Van de Peer Y, Michael T (2008) Analysis of a Gibbs sampler method for model-based clustering of gene expression data. Bioinformatics 24: 176–183.

55. Hooghe B, Hulpiau P, van Roy F, De Bleser P (2008) ConTra: a promoter alignment analysis tool for identification of transcription factor binding sites across species. Nucleic Acids Res 36: W128–132.

56. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, et al. (2005) Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature 433: 769–773.

<figure><img src="image_17.png" /></figure>