Main

MicroRNAs (miRNAs) are a class of non-coding RNAs able to regulate gene expression at the post-transcriptional level, by binding to the 3′ untranslated region of target messenger RNAs (mRNAs) through partial sequence homology, and causing a block of translation and/or mRNA degradation (He and Hannon, 2004). At the time of writing of this paper, 721 human miRNA genes had been described (http://microrna.sanger.ac.uk/sequences) and it has been estimated that each of the miRNAs targets about 100 different mRNA molecules (Brennecke et al, 2005; Lewis et al, 2005; Lim et al, 2005). MiRNAs have important roles in essential processes, such as differentiation, cell growth, stress response and cell death (Miska, 2005; Zamore and Haley, 2005). Accordingly, altered miRNA expression is likely to contribute to human disease, including cancer. In fact, the spectrum of miRNAs expressed in solid cancers is very different from that of normal cells and the predicted targets for the differentially expressed miRNAs are significantly enriched for protein-coding tumour suppressors and oncogenes (Lu et al, 2005; Volinia et al, 2006).

Focussing on breast cancer, levels of specific miRNAs differ between malignant and normal breast tissue and are able to classify tumours according to clinicopathological variables, such as proliferation index, steroid hormone receptor and Her2/neu status, nodal status and tumour stage (Iorio et al, 2005; Mattie et al, 2006; Lowery et al, 2008). Furthermore, miRNAs are differentially expressed between molecular breast cancer subtypes, including luminal A, luminal B, basal-like and Her2+ (Blenkiron et al, 2007). This highlights the potential of miRNA signatures as new prognostic indicators that may contribute to the improved selection of patients for adjuvant therapy. Indeed, Foekens et al (2008) recently linked four miRNAs (miR-7, miR-128a, miR-210 and miR-516-3p) to a shorter time to distant metastasis in oestrogen receptor (ER)-positive, lymph node-negative breast cancer.

In this study, we explored whether inflammatory breast cancer (IBC) is associated with a specific miRNA signature. Although this form of breast cancer has traditionally been understudied, there are several reasons for focusing on IBC. First, IBC is arguably the deadliest form of breast cancer, with 5- and 10-year disease-free survival rates of less than 35 and 20%, respectively (Jaiyesimi et al, 1992), highlighting the need for accurate prognostic and predictive markers and the development of new molecular targeted therapies. Second, most patients with IBC have lymph node involvement at time of diagnosis and one-third of patients present with distant metastasis (Jaiyesimi et al, 1992), which makes IBC an interesting model for identifying the forces driving breast cancer aggressiveness in general. Indeed, several studies comparing IBC to non-IBC have demonstrated a specific mRNA expression signature in IBC that is significantly enriched for genes involved in cell motility, inflammation, immune response and stem cell biology (Bertucci et al, 2004, 2005; Van Laere et al, 2005, 2006, 2007). Furthermore, we recently reported that the identified IBC gene signature was associated with dismal tumour characteristics in non-IBC, such as high tumour grade, absence of ER expression and presence of Her2/neu expression, and independently predicted relapse-free survival in non-IBC (Van Laere et al, 2008). Finally, accurate and rapid diagnosis of IBC remains problematic due to the lack of uniformity in criteria for IBC diagnosis and new diagnostic markers are therefore needed.

Materials and Methods

Patients’ samples

Tumour samples were retrieved from the tissue bank of the General Hospital Sint-Augustinus (Antwerp, Belgium). Clinical and pathological data are stored in a database in accordance with hospital privacy rules. Specimens were brought to the pathologists immediately after resection and part of the tissue was placed in liquid nitrogen and subsequently stored at −180°C. A total of 20 patients with IBC and 50 patients with non-IBC were included in this study. Inflammatory breast cancer was diagnosed according to the criteria mentioned in the AJCC (American Joint Committee on Cancer)-TNM staging system (Singletary et al, 2002). All patients with IBC showed diffuse enlargement of the involved breast of sudden onset. There was erythema and oedema of the skin involving more than one-third of the breast. The presence of dermal lymphatic invasion as an isolated observation was not sufficient for the diagnosis of IBC and was not necessary for the diagnosis either. The non-IBC control group was matched for histological tumour grade, ER expression and human epidermal growth factor receptor 2 (HER2) amplification but not for tumour stage to identify true determinants of IBC. Tumour characteristics are provided in Table 1.

Table 1 Tumour characteristics
Full size table

Total RNA isolation, complementary DNA synthesis and miRNA quantification

After tissue disruption, total RNA was extracted by using the mirVana miRNA Isolation Kit (Ambion, Austin, TX, USA) according to the manufacturer's instructions for total RNA isolation. Purified total RNA was quantified by using a NanoDrop ND1000 (NanoDrop Technologies, Wilmington, DE, USA). Total RNA (100 ng) was then converted to complementary DNA by priming with a mixture of stem-looped RT primers (Megaplex RT Primers, Human Pool A, Applied Biosystems, Foster City, CA, USA) in combination with the TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems), allowing simultaneous transcription of 377 unique miRNAs and six endogenous controls. Briefly, 3 μl of total RNA was supplemented with RT primer mix (10 ×), dNTPs with dTTP (100 mM), Multiscribe Reverse Transcriptase (50 U μl−1), RT buffer (10 ×), MgCl2 (25 mM) and RNase inhibitor (20 U μl−1) in a total reaction volume of 7.5 μl. Thermal-cycling conditions were as follows: 40 cycles of 16°C for 2 min, 42°C for 1 min and 50°C for 1 s, followed by a reverse transcriptase inactivation at 85°C for 5 min. The Megaplex RT product (2.5 μl) was pre-amplified using the TaqMan PreAmp Master Mix (Applied Biosystems) and PreAmp Primers, Human Pool A (Applied Biosystems) in a 25-μl PCR reaction. Thermal-cycling conditions were as follows: 95°C for 10 min, 55°C for 2 min and 75°C for 2 min, followed by 12 cycles of 95°C for 15 s and 60°C for 4 min. MiRNA quantification was performed with the TaqMan Human MicroRNA A Array Set v2.0 (Applied Biosystems), which contains 384 TaqMan miRNA assays. The PreAmp product was diluted four-fold. Each of the eight wells was loaded with 100 μl of PCR reaction mix, containing 50 μl of TaqMan Universal PCR Master Mix, no AmpErase uracil N-glycosylase (UNG) (Applied Biosystems), 1 μl of diluted PreAmp product and 49 μl of nuclease-free water. Thermal-cycling conditions were as follows: 94.5°C for 10 min, followed by 40 cycles of 97°C for 30 s and 59.7°C for 1 min. All PCR reactions were performed on a 7900HT Fast Real-Time PCR System (Applied Biosystems). The PCR replicates were measured for four samples.

Statistics and bioinformatics

To reduce technical variation, all miRNAs with a Ct (threshold cycle) detection cut-off of less than 35 PCR cycles in at least 25% of samples were filtered, resulting in a list of 322 informative miRNAs. Next, we performed a between-sample normalisation by median-centring the distribution of miRNA expression levels for each sample (Mestdagh et al, 2009). Relative miRNA expression levels were calculated using the ΔCt method (Livak and Schmittgen, 2001).

To identify miRNA expression patterns, we performed an unsupervised hierarchical complete linkage cluster analysis using Euclidean distance as dissimilarity metric. The number of robust sample clusters was determined with the silhouette algorithm. Next, we investigated the association of the sample clustering result with clinicopathological variables using a Pearson χ2 test and a Goeman global test (Goeman et al, 2004).

To identify individual miRNAs associated with IBC, we first selected all miRNA that were significantly over- or underexpressed in IBC when compared with non-IBC using a non-parametric test. Then, each of these miRNAs was subjected to a multivariate regression analysis with N status, M status, tumour stage and HER2 amplification as covariates. In this way, we could indentify miRNAs that are specifically related to IBC and not to differences in clinicopathological variables between IBC and non-IBC known to also influence miRNA expression in breast cancer.

To investigate the biological relevance of the identified miRNAs, we adopted the strategy recently described by Wang et al (2009). First, for each of the miRNAs that were independently associated with IBC, we identified a set of putative target genes by Spearman correlation analysis, taking into account both positive and negative correlations. This analysis was performed on a subgroup of 44 samples (20 IBC and 24 non-IBC samples) for which Affymetrix HGU133 plus 2.0 gene expression profiles were available (Van Laere et al, 2007). Next, to identify the true miRNA targets, we compared the lists of putative miRNA targets with the lists of miRNA targets defined by TargetScan5.1 (http://www.targetscan.org) using a hypergeometric gene set enrichment analysis. Those miRNAs for which the miRNA targets defined by TargetScan5.1 were significantly enriched for correlation-defined miRNA targets were subjected to further analysis. In addition, all further analyses were performed only on the common genes between the two miRNA targets sets. These genes were analysed for enriched Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) categories through hypergeometric gene set enrichment analysis.

To test whether the identified miRNAs are prognostically relevant, we downloaded six publicly available gene expression data sets from the National Center for Biotechnology Information website (GSE7390, GSE9195, GSE1456, GSE11121, GSE2034 and GSE4922; http://www.ncbi.nlm.nih.gov). In addition, the data set described by van de Vijver et al (2002) was retrieved from the Rosetta website (http://www.rii.com). In total, we analysed mRNA expression data from 1504 breast cancer samples. For each sample, we calculated a score, proportional to the level of expression of a selected miRNA, by subtracting the average of the negatively correlated miRNA targets from the average of the positively correlated miRNA targets. These scores were standardised (median=0 and s.d.=1) to allow for comparisons across different data sets and analysed using a Cox proportional hazards model in both a univariate and a multivariate setting. The outcome variable was distant metastases-free survival (DMFS, N=1059), overall survival (OS, N=652) or relapse-free survival (RFS, N=1145) where possible.

Finally, within the subset of 44 samples for which Affymetrix HGU133 plus 2.0 gene expression profiles were available, we compared the expression of miRNA processing genes between IBC and non-IBC. Furthermore, we analysed the putative regulatory effect of all miRNAs from the TargetScan5.1 database on the gene expression profiles from these 44 samples. We adopted the approach recently published by Cheng et al (2009) to calculate a regulatory effect (RE) score by subtracting the average rank of the miRNA targets from the average rank from the non-miRNA targets, with high RE scores denoting a strong effect of the miRNA on the expression of the targets and vice versa. However, as this score only takes into account the negatively regulated miRNA targets, we adapted the RE score in a manner that both inhibitory and activating effects were equally weighed. Next, these RE scores were compared between IBC and non-IBC using a non-parametrical test.

All data analyses were performed using Bioconductor in R (http://www.bioconductor.org). Correction for multiple testing was performed using the Benjamini and Hochberg step-up false discovery rate controlling procedure and adjusted P-values <0.1 were considered significant.

Results

Hierarchical clustering based on miRNA expression profiles

We used quantitative RT-PCR technology in combination with a limited cycle pre-amplification to evaluate the miRNA expression profiles of 20 IBC and 50 non-IBC samples. The TaqMan Human MicroRNA A Array v2.0 contains 384 TaqMan miRNA assays enabling accurate quantification of 377 human miRNAs, six endogenous controls and one negative control. Four breast cancer samples were assayed in duplicate. Pearson correlation coefficients for these duplicates ranged from 0.98 to 0.99, indicating good assay reproducibility.

The accuracy of quantitative RT-PCR experiments is dependent on proper data normalisation. For each individual breast cancer sample, we calculated the mean expression value based on those miRNAs that were expressed according to a Ct detection cut-off of less than 35 PCR cycles in at least 25% of samples (N=322). This value was subsequently used as a normalisation factor to reduce technical variation. This normalisation strategy has been shown to be superior to the use of stable small RNA controls (Mestdagh et al, 2009).

The clustering of miRNA expression profiles derived from 70 breast cancer samples is shown in Figure 1. The dendrogram was constructed by using the 50 most varying miRNAs among all breast cancer samples (based on the s.d. value of miRNA levels across all samples). We identified three sample clusters with an average silhouette width of 0.08 (P<0.05 after permutation testing). A Pearson χ2 test demonstrated that sample clustering was associated with ER expression (P=0.011), histological tumour grade (P=0.004), N status (P=0.014) and tumour stage (P=0.010). No association with sample clustering was observed for T and M status, HER2 amplification or tumour subtype (IBC or non-IBC).

Figure 1
figure 1

Hierarchical clustering of 20 IBC and 50 non-IBC samples according to the expression pattern of the 50 most varying miRNAs. Expression values for these 50 miRNAs are represented in a matrix format, with rows indicating miRNAs and columns indicating samples. High expression values are colour-coded red and low expression values are colour-coded blue. Three robust sample clusters were identified, which were significantly associated with ER expression and histological grade. In particular, the combined first two (blue and yellow) sample clusters were enriched for ER+ breast tumours (80% of samples) when compared with the third (red) sample cluster (50% of samples) (P χ2=0.028). Notably, in the first (blue) sample cluster, 20% of samples were poorly differentiated compared with 60% of samples in the second (yellow) sample cluster (P χ2=0.004), suggestive of a subdivision of ER+ breast tumours according to the luminal A and luminal B subtype. No association of sample clustering with the difference between IBC and non-IBC was observed: 30, 25 and 45% of IBC samples grouped together in the first (blue), second (yellow) and third (red) sample cluster, respectively (P χ2=0.082). (The colour reproduction of this figure is available on the html full text version of the manuscript.)

Full size image

A Goeman global test was performed to look for associations of miRNA expression with various clinicopathological factors (Table 2), for the reduced data set consisting of the 50 most varying miRNAs. Significant associations were observed for ER expression and histological tumour grade, but not for T, N or M status, tumour stage, HER2 amplification or tumour subtype (IBC or non-IBC). Thus, the largest variation in miRNA expression in our set of breast cancer samples was attributable to differences in steroid hormone receptor expression and histological tumour grade and not to differences between IBC and non-IBC. In fact, the overall distributions of miRNA expression values in IBC and non-IBC look very similar (Figure 2).

Table 2 Association of global miRNA expression with clinicopathological factors
Full size table
Figure 2
figure 2

Distribution of miRNA expression values in (A) inflammatory breast tumours and (B) non-inflammatory breast tumours.

Full size image

Association of individual miRNAs with inflammatory breast cancer

A logistic regression analysis was performed to identify differences in individual miRNA expression levels between IBC and non-IBC samples (Table 3). Few miRNAs were observed to be independently associated with the difference between IBC and non-IBC. Increased expression of miR-335, miR-337-5p, miR-451, miR-486-3p, miR-520a-5p and miR-548d-5p was observed in the IBC subtype, whereas miR-15a, miR-24, miR-29a, miR-30b, miR-320, miR-342-5p and miR-432-3p were significantly downregulated in comparison with non-IBC.

Table 3 Association of miRNA expression with clinicopathological factors
Full size table

A similar analysis was performed to identify differences in miRNA expression between samples according to T status (low vs high), N status (positive vs negative), M status, tumour stage (low vs high), tumour grade (high vs low), ER expression and HER2 amplification. The associations of miRNAs with these clinicopathological factors are also shown in Table 3.

Prediction of miRNA targets and their involvement in biological processes

Next, we investigated the involvement of miRNA target genes in various biological processes, which may indicate the function of the miRNA, by adopting the strategy recently described by Wang et al (2009). First, we performed a correlation coefficient analysis to evaluate potential correlations between expression levels of the 13 miRNAs associated with IBC (miR-335, miR-337-5p, miR-451, miR-486-3p, miR-520a-5p, miR-548d-5p, miR-15a, miR-24, miR-29a, miR-30b, miR-320, miR-342-5p and miR-342-3p) and mRNA expression levels of 10 961 known genes in a reduced set of 44 breast cancer samples, of which 20 were IBC and 24 were non-IBC. Using a false discovery rate of <0.1, we detected significant correlations in 17 295 miRNA–mRNA pairs (Table 4). Of the 17 295 pairs, 7012 pairs and 10 283 pairs showed negative and positive correlations, respectively. For all miRNAs, with the exception of miR-520a-5p and miR-337-5p, we observed more frequently positively correlated miRNA–mRNA pairs than negatively correlated miRNA–mRNA pairs. Subsequently, we performed a gene set enrichment analysis to investigate whether these sets of miRNA-correlated genes are indeed differentially expressed in IBC and non-IBC (Table 4). Significant results were observed for all but one miRNA-correlated gene set. Expression levels of the miR-335/miR-548d-5p/miR-451-correlated gene sets were increased in IBC, whereas expression levels of the miR-15a/miR-24/miR-29a/miR-30b/miR-320/miR-342-5p/miR-342-3p/miR-337-5p/miR-520a-5p-correlated gene sets were decreased in IBC.

Table 4 Correlations between miRNAs and mRNAs in breast cancer (N=44)
Full size table

To test whether the identified miRNA-correlated genes are direct miRNA targets and not downstream genes of the miRNA targets, we downloaded the predicted miRNA targets from TargetScan5.1 and compared them with the miRNA-correlated genes. The percentages of overlapping genes between the two lists ranged from 0 to 10%. We observed four miRNAs (miR-29a, miR-30b, miR-342-3p and miR-520a-5p) of which targets were significantly enriched among their correlated genes using a false discovery rate of <0.1 (Table 4). The number of overlapping genes was 121 for miR-29a, 140 for miR-30b, 19 for miR-342-3p and 13 for miR-520a-5p (see Table 5 for all target gene symbols).

Table 5 Gene symbols for direct miRNA target genes
Full size table

To functionally classify the direct target genes of miR-29a, miR-30b, miR-342-3p and miR-520a-5p, we performed a gene set enrichment analysis through hypergeometric testing. We focussed on the GO and KEGG gene set collections. For each of the gene sets, we observed various degrees of GO and KEGG term enrichments. The top 10 terms for GO categories or KEGG pathways are listed in Table 6 for each of the four miRNA target gene sets. For miR-29a target genes, the most marked GO terms related to DNA methyltransferase activity. MiR-30b target genes mostly associated with GO terms related to the insulin receptor-signalling pathway. In the miR-342-3p target gene set, cell proliferation-related GO terms were overrepresented. For miR-520a-5p target genes, the most marked GO term was regulation of cell growth by extracellular stimulus.

Table 6 Association of direct miRNA target genes with biological processes
Full size table

Association of miRNA targets with prognosis in breast cancer

To investigate the potential prognostic relevance of miR-29a, miR-30b, miR-342-3p and miR-520a-5p in breast cancer, we performed a Cox regression analysis for seven publicly available gene expression data sets for which information regarding DMFS, RFS or OS was available using an indirect approach. First, we calculated for each sample an miRNA target gene expression score by subtracting the average relative gene expression value of the negatively correlated miRNA target genes from the average relative gene expression value of the positively correlated miRNA target genes (Figure 3). The resulting score was standardised (median=0, s.d.=1) to allow for comparison across different data sets. The miRNA target gene expression scores calculated as such were all strongly correlated with the miRNA expression values, thereby validating the score (range Spearman correlation coefficients: 0.70–0.76; P<0.0001). The association between outcome and the miRNA target gene expression score was evaluated using Kaplan–Meier and Cox regression analysis (Table 7). Overall, we observed significant associations between miRNA target gene expression and patient outcome (DMFS, RFS and OS) for miR-29a, miR-30b and miR-520a-5p, but not for miR-342-3p. In particular, high levels of miR-520a-5p target genes correlated with shorter DMFS, RFS and OS in, respectively, 3 out of 4, 3 out of 5 and 3 out of 3 data sets and performed best in a Cox regression analysis.

Figure 3
figure 3

Kaplan–Meier survival analysis of miR-520a-5p target gene expression in breast cancer with distant metastasis-free (A and C) and overall survival (B and D) as outcome. Kaplan–Meier curves are shown for the data sets of van de Vijver et al (2002) (A and B) and of Desmedt et al (2007) (C and D).

Full size image
Table 7 Association of direct miRNA target genes with prognosis in breast cancer
Full size table

Expression of miRNA processing genes in breast cancer

Recently, Cheng et al (2009) reported that the post-transcriptional regulation of miRNA expression may be important for the regulatory effect of miRNAs on their targets (Cheng et al, 2009). We therefore examined whether miRNA processing genes (trbp2, dicer, ago1, ago2 and drosha) are differentially expressed between IBC and non-IBC. We observed that among the miRNA processing genes, ago2 was significantly upregulated and dicer significantly downregulated in IBC (N=20) compared with non-IBC samples (N=24), with P-values of 0.002 and 0.004 (false discovery rate <0.1), respectively. We further examined whether this altered expression of ago2 and dicer in IBC when compared with non-IBC is also reflected in altered regulatory effects of miRNAs on their targets. We calculated an RE score for each of 153 representative miRNA families from TargetScan 5.1 in all 44 samples using a strategy adapted from Cheng et al (2009). Then, we compared RE scores between IBC and non-IBC samples. Of 153 miRNA families, 104 (68%) showed higher RE scores in IBC and 49 (32%) showed higher RE scores in non-IBC. Using the significance analysis of microarrays method, we observed 20 significant RE-changing miRNAs (miRNAs that show different regulatory effects between IBC and non-IBC).

Discussion

MiRNAs are a recently discovered class of small regulatory RNAs, which influence the stability and translational efficiency of target mRNAs (Verghese et al, 2008). In this study, we examined the expression of 384 miRNAs in 70 breast cancer samples to look for miRNAs, of which expression levels are significantly different in IBC compared with non-IBC. The spectrum of expressed miRNAs mostly varied according to steroid hormone receptor expression and histological tumour grade and not according to tumour subtype (IBC vs non-IBC). However, we did identify a set of 13 miRNAs for which expression levels were able to correctly predict the nature of the sample analysed (IBC vs non-IBC). Six of them, miR-335, miR-337-5p, miR-451, miR-486-3p, miR-520a-5p and miR-548d-5p, were upregulated in IBC and the remaining seven, miR-15a, miR-24, miR-29a, miR-30b, miR-320, miR-342-5p and miR-342-3p, were downregulated in IBC. The observation that few miRNAs are specifically associated with IBC does not support the hypothesis of a prominent role for altered miRNA expression in the phenotype of IBC. However, miRNA expression levels are not necessarily representative of the regulatory ability of an miRNA. Cheng et al (2009) recently reported that the inhibitory effects of miRNAs on their targets differs between two breast cancer subtypes (ER− and ER+) due to a differential expression of several key miRNA processing genes and not due to differences in miRNA expression levels (Cheng et al, 2009). We therefore examined the expression levels of genes in the miRNA biogenesis pathway in samples from IBC and non-IBC. We observed significantly lower expression levels for dicer, a ribonuclease, which cleaves the pre-miRNA into a single-stranded mature miRNA and significantly higher expression levels for ago2, the catalytic endonuclease of the RNA-induced silencing complex in IBC when compared with non-IBC. The latter observation may suggest higher RNA-induced silencing complex activities, and therefore may further suggest that miRNAs regulate target gene expression in IBC with higher efficiency. Indeed, miRNAs seem to have a higher regulatory effect (reflected by the RE score) in IBC samples than in non-IBC samples. It has been previously observed that low dicer expression and high ago2 expression in breast cancer is associated with the more aggressive basal-like, HER2+ and luminal B subtypes (Blenkiron et al, 2007), which are known to be overrepresented in IBC (Van Laere et al, 2006). Moreover, decreased dicer expression has been reported to be associated with poor clinical outcome in breast, ovarian and lung cancer (Karube et al, 2005; Merritt et al, 2008; Grelier et al, 2009).

At present, the lack of knowledge about bona fide miRNA target genes hampers a full understanding on the biological functions deregulated by aberrant miRNA expression. To overcome this limitation, we adopted the strategy recently described by Wang et al (2009) for the miRNAs of which expression differed between IBC and non-IBC (Wang et al, 2009). First, we used a whole genome approach to identify thousands of highly correlated miRNA–mRNA pairs. We identified a large number of both negative and positive correlations. The detection of a positive correlation between miRNA and mRNA levels suggests a positive regulatory role for miRNAs, as has been reported for some genes (Vasudevan et al, 2007). Second, we made use of a computational approach to search for predicted miRNA targets. Concordant genes are direct miRNA targets that are tightly correlated with fluctuations in miRNA expression. We observed significant concordance between miRNA-correlated genes and miRNA-predicted targets in only 4 of 13 miRNA sets (miR-29a, miR-30b, miR-342-3p and miR-520a-5p). Thus, we did not see any evidence of concordance in a majority of miRNA sets. Generally, miRNAs are believed to bind 3′ untranslated region of a target gene and regulate gene expression at protein level (Wang et al, 2009). Therefore, miRNA targets themselves may not demonstrate noticeable change at the mRNA level. However, this can also be explained by discordance in changes of expression between the key miRNA processing genes. For miR-29a, miR-30b, miR-342-3p and miR-520a-5p target genes, we were able to detect a variety of biological processes that may indicate function of those miRNAs. MiR-29a target genes were mainly related to DNA methyltransferase activity. Indeed, expression of the miR-29 family has been shown to target DNA methyltransferase 3A and 3B expression in lung cancer tissues (Fabbri et al, 2007) and in acute myeloid leukaemia (Garzon et al, 2009), inducing global hypomethylation. For miR-30b, we observed target genes related to insulin receptor signalling. Recently, it has been reported that another member of the miRNA-30 family, miR-30d, is upregulated by glucose and increases insulin gene expression (Tang et al, 2009). MiR-342-3p and miR-520a-5p target genes proved to be involved in cell proliferation.

As IBC is regarded as a model for breast cancer aggressiveness, we were interested in determining whether the IBC-specific miRNAs were associated with a dismal prognosis in non-IBC. Given the limited number of our clinical samples, we used an indirect approach and correlated miRNA target gene expression with patients’ outcome using publicly available information from microarrays studies, which could be extracted from the gene expression omnibus database. This analysis demonstrated a marked association of miR-520a-5p target gene expression with a shorter DMFS, RFS or OS in most breast cancer data sets that were investigated. This result demonstrates a possible role of miR-520a-5p as a prognostic factor in breast cancer. Although we acknowledge the limitations of using an indirect approach, these kinds of analyses may help to guide future large-scale studies on the prognostic value of miRNAs in breast cancer.

In conclusion, this study, to the best of our knowledge, represents the first integrated analysis of miRNA and mRNA expression in IBC. We identified a number of miRNAs that are differentially expressed between IBC and non-IBC. Furthermore, we reported altered expression of miRNA processing genes in IBC when compared with non-IBC. For four of the IBC-related miRNAs, we were able to detect a variety of biological processes that may indicate function of these miRNAs and to indicate their potential association with prognosis in breast cancer based on the expression levels of their target genes.