The OSMR Gene Is Involved in Hirschsprung Associated Enterocolitis Susceptibility through an Altered Downstream Signaling.

Hirschsprung (HSCR) Associated Enterocolitis (HAEC) is a common life-threatening complication in HSCR. HAEC is suggested to be due to a loss of gut homeostasis caused by impairment of immune system, barrier defense, and microbiome, likely related to genetic causes. No gene has been claimed to contribute to HAEC occurrence, yet. Genetic investigation of HAEC by Whole-Exome Sequencing (WES) on 24 HSCR patients affected (HAEC) or not affected (HSCR-only) by enterocolitis and replication of results on a larger panel of patients allowed the identification of the HAEC susceptibility variant p.H187Q in the Oncostatin-M receptor (OSMR) gene (14.6% in HAEC and 5.1% in HSCR-only, p = 0.0024). Proteomic analysis on the lymphoblastoid cell lines from one HAEC patient homozygote for this variant and one HAEC patient not carrying the variant revealed two well distinct clusters of proteins significantly up or downregulated upon OSM stimulation. A marked enrichment in immune response pathways (q < 0.0001) was shown in the HAEC H187 cell line, while proteins upregulated in the HAEC Q187 lymphoblasts sustained pathways likely involved in pathogen infection and inflammation. In conclusion, OSMR p.H187Q is an HAEC susceptibility variant and perturbates the downstream signaling cascade necessary for the gut immune response and homeostasis maintenance.

Post-transcriptional regulation of human breast cancer cell proteome by unliganded estrogen receptor ß via microRNAs.

Estrogen receptor ß (ERß) is a member of the nuclear receptor family of homeostatic regulators that is frequently lost in breast cancer (BC), where its presence correlates with a better prognosis and a less aggressive clinical outcome of the disease. In contrast to ERα, its closest homolog, ERß shows significant estrogen-independent activities, including the ability to inhibit cell cycle progression and regulate gene transcription in the absence of the ligand. Investigating the nature and extent of this constitutive activity of ERß in BC MCF-7 and ZR-75.1 cells by means of microRNA (miRNA) sequencing, we identified 30 miRNAs differentially expressed in ERß+ versus ERß- cells in the absence of ligand, including up-regulated oncosuppressor miRs such miR-30a. In addition, a significant fraction of >1,600 unique proteins identified in MCF-7 cells by iTRAQ quantitative proteomics were either increased or decreased by ERß, revealing regulation of multiple cell pathways by ligand-free receptors. Transcriptome analysis showed that for a large number of proteins regulated by ERß, the corresponding mRNAs are unaffected, including a large number of putative targets of ERß-regulated miRNAs, indicating a central role of miRNAs in mediating BC cell proteome regulation by ERß. Expression of a mimic of miR-30a-5p, a direct target and downstream effector of ERß in BC, led to the identification of several target transcripts of this miRNA, including 11 encoding proteins whose intracellular concentration was significantly affected by unliganded receptor. These results demonstrate a significant effect of ligand-free ERß on BC cell functions via modulation of the cell proteome and suggest that miRNA regulation might represent a key event in the control of the biological and clinical phenotype of hormone-responsive BC by this nuclear receptor.

IBTK Differently Modulates Gene Expression and RNA Splicing in HeLa and K562 Cells.

The IBTK gene encodes the major protein isoform IBTKα that was recently characterized as substrate receptor of Cul3-dependent E3 ligase, regulating ubiquitination coupled to proteasomal degradation of Pdcd4, an inhibitor of translation. Due to the presence of Ankyrin-BTB-RCC1 domains that mediate several protein-protein interactions, IBTKα could exert expanded regulatory roles, including interaction with transcription regulators. To verify the effects of IBTKα on gene expression, we analyzed HeLa and K562 cell transcriptomes by RNA-Sequencing before and after IBTK knock-down by shRNA transduction. In HeLa cells, 1285 (2.03%) of 63,128 mapped transcripts were differentially expressed in IBTK-shRNA-transduced cells, as compared to cells treated with control-shRNA, with 587 upregulated (45.7%) and 698 downregulated (54.3%) RNAs. In K562 cells, 1959 (3.1%) of 63128 mapped RNAs were differentially expressed in IBTK-shRNA-transduced cells, including 1053 upregulated (53.7%) and 906 downregulated (46.3%). Only 137 transcripts (0.22%) were commonly deregulated by IBTK silencing in both HeLa and K562 cells, indicating that most IBTKα effects on gene expression are cell type-specific. Based on gene ontology classification, the genes responsive to IBTK are involved in different biological processes, including in particular chromatin and nucleosomal organization, gene expression regulation, and cellular traffic and migration. In addition, IBTK RNA interference affected RNA maturation in both cell lines, as shown by the evidence of alternative 3'- and 5'-splicing, mutually exclusive exons, retained introns, and skipped exons. Altogether, these results indicate that IBTK differently modulates gene expression and RNA splicing in HeLa and K562 cells, demonstrating a novel biological role of this protein.

iMir: an integrated pipeline for high-throughput analysis of small non-coding RNA data obtained by smallRNA-Seq.

BACKGROUND: Qualitative and quantitative analysis of small non-coding RNAs by next generation sequencing (smallRNA-Seq) represents a novel technology increasingly used to investigate with high sensitivity and specificity RNA population comprising microRNAs and other regulatory small transcripts. Analysis of smallRNA-Seq data to gather biologically relevant information, i.e. detection and differential expression analysis of known and novel non-coding RNAs, target prediction, etc., requires implementation of multiple statistical and bioinformatics tools from different sources, each focusing on a specific step of the analysis pipeline. As a consequence, the analytical workflow is slowed down by the need for continuous interventions by the operator, a critical factor when large numbers of datasets need to be analyzed at once. RESULTS: We designed a novel modular pipeline (iMir) for comprehensive analysis of smallRNA-Seq data, comprising specific tools for adapter trimming, quality filtering, differential expression analysis, biological target prediction and other useful options by integrating multiple open source modules and resources in an automated workflow. As statistics is crucial in deep-sequencing data analysis, we devised and integrated in iMir tools based on different statistical approaches to allow the operator to analyze data rigorously. The pipeline created here proved to be efficient and time-saving than currently available methods and, in addition, flexible enough to allow the user to select the preferred combination of analytical steps. We present here the results obtained by applying this pipeline to analyze simultaneously 6 smallRNA-Seq datasets from either exponentially growing or growth-arrested human breast cancer MCF-7 cells, that led to the rapid and accurate identification, quantitation and differential expression analysis of ~450 miRNAs, including several novel miRNAs and isomiRs, as well as identification of the putative mRNA targets of differentially expressed miRNAs. In addition, iMir allowed also the identification of ~70 piRNAs (piwi-interacting RNAs), some of which differentially expressed in proliferating vs growth arrested cells. CONCLUSION: The integrated data analysis pipeline described here is based on a reliable, flexible and fully automated workflow, useful to rapidly and efficiently analyze high-throughput smallRNA-Seq data, such as those produced by the most recent high-performance next generation sequencers. iMir is available at http://www.labmedmolge.unisa.it/inglese/research/imir.

Molecular mechanisms of selective estrogen receptor modulator activity in human breast cancer cells: identification of novel nuclear cofactors of antiestrogen-ERα complexes by interaction proteomics.

Estrogen receptor alpha (ERα) is a ligand-activated transcription factor that controls key cellular pathways via protein-protein interactions involving multiple components of transcriptional coregulator and signal transduction complexes. Natural and synthetic ERα ligands are classified as agonists (17ß-estradiol/E(2)), selective estrogen receptor modulators (SERMs: Tamoxifen/Tam and Raloxifene/Ral), and pure antagonists (ICI 182,780-Fulvestrant/ICI), according to the response they elicit in hormone-responsive cells. Crystallographic analyses reveal ligand-dependent ERα conformations, characterized by specific surface docking sites for functional protein-protein interactions, whose identification is needed to understand antiestrogen effects on estrogen target tissues, in particular breast cancer (BC). Tandem affinity purification (TAP) coupled to mass spectrometry was applied here to map nuclear ERα interactomes dependent upon different classes of ligands in hormone-responsive BC cells. Comparative analyses of agonist (E(2))- vs antagonist (Tam, Ral or ICI)-bound ERα interacting proteins reveal significant differences among ER ligands that relate with their biological activity, identifying novel functional partners of antiestrogen-ERα complexes in human BC cell nuclei. In particular, the E(2)-dependent nuclear ERα interactome is different and more complex than those elicited by Tam, Ral, or ICI, which, in turn, are significantly divergent from each other, a result that provides clues to explain the pharmacological specificities of these compounds.

Global analysis of estrogen receptor beta binding to breast cancer cell genome reveals an extensive interplay with estrogen receptor alpha for target gene regulation.

BACKGROUND: Estrogen receptors alpha (ERα) and beta (ERß) are transcription factors (TFs) that mediate estrogen signaling and define the hormone-responsive phenotype of breast cancer (BC). The two receptors can be found co-expressed and play specific, often opposite, roles, with ERß being able to modulate the effects of ERα on gene transcription and cell proliferation. ERß is frequently lost in BC, where its presence generally correlates with a better prognosis of the disease. The identification of the genomic targets of ERß in hormone-responsive BC cells is thus a critical step to elucidate the roles of this receptor in estrogen signaling and tumor cell biology. RESULTS: Expression of full-length ERß in hormone-responsive, ERα-positive MCF-7 cells resulted in a marked reduction in cell proliferation in response to estrogen and marked effects on the cell transcriptome. By ChIP-Seq we identified 9702 ERß and 6024 ERα binding sites in estrogen-stimulated cells, comprising sites occupied by either ERß, ERα or both ER subtypes. A search for TF binding matrices revealed that the majority of the binding sites identified comprise one or more Estrogen Response Element and the remaining show binding matrixes for other TFs known to mediate ER interaction with chromatin by tethering, including AP2, E2F and SP1. Of 921 genes differentially regulated by estrogen in ERß+ vs ERß- cells, 424 showed one or more ERß site within 10 kb. These putative primary ERß target genes control cell proliferation, death, differentiation, motility and adhesion, signal transduction and transcription, key cellular processes that might explain the biological and clinical phenotype of tumors expressing this ER subtype. ERß binding in close proximity of several miRNA genes and in the mitochondrial genome, suggests the possible involvement of this receptor in small non-coding RNA biogenesis and mitochondrial genome functions. CONCLUSIONS: Results indicate that the vast majority of the genomic targets of ERß can bind also ERα, suggesting that the overall action of ERß on the genome of hormone-responsive BC cells depends mainly on the relative concentration of both ERs in the cell.

Corrigendum: Dual-mTOR Inhibitor Rapalink-1 Reduces Prostate Cancer Patient-Derived Xenograft Growth and Alters Tumor Heterogeneity.

[This corrects the article .].

Dual-mTOR Inhibitor Rapalink-1 Reduces Prostate Cancer Patient-Derived Xenograft Growth and Alters Tumor Heterogeneity.

Bone metastasis is the leading cause of prostate cancer (PCa) mortality, frequently marking the progression to castration-resistant PCa. Dysregulation of the androgen receptor pathway is a common feature of castration-resistant PCa, frequently appearing in association with mTOR pathway deregulations. Advanced PCa is also characterized by increased tumor heterogeneity and cancer stem cell (CSC) frequency. CSC-targeted therapy is currently being explored in advanced PCa, with the aim of reducing cancer clonal divergence and preventing disease progression. In this study, we compared the molecular pathways enriched in a set of bone metastasis from breast and prostate cancer from snap-frozen tissue. To further model PCa drug resistance mechanisms, we used two patient-derived xenografts (PDX) models of bone-metastatic PCa, BM18, and LAPC9. We developed in vitro organoids assay and ex vivo tumor slice drug assays to investigate the effects of mTOR- and CSC-targeting compounds. We found that both PDXs could be effectively targeted by treatment with the bivalent mTORC1/2 inhibitor Rapalink-1. Exposure of LAPC9 to Rapalink-1 but not to the CSC-targeting drug disulfiram blocked mTORC1/2 signaling, diminished expression of metabolic enzymes involved in glutamine and lipid metabolism and reduced the fraction of CD44+ and ALDEFluorhigh cells, in vitro. Mice treated with Rapalink-1 showed a significantly delayed tumor growth compared to control and cells recovered from the tumors of treated animals showed a marked decrease of CD44 expression. Taken together these results highlight the increased dependence of advanced PCa on the mTOR pathway, supporting the development of a targeted approach for advanced, bone metastatic PCa.

Deciphering the clonal relationship between glandular and squamous components in adenosquamous carcinoma of the lung using whole exome sequencing.

Adenosquamous carcinoma of the lung (ASC) is a rare subtype of non-small cell lung cancer, consisting of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) components. ASC shows morphological characteristics of classic LUAD and LUSC but behaves more aggressively. Although ASC can serve as a model of lung cancer heterogeneity and transdifferentiation, its genomic background remains poorly understood. In this study, we sought to explore the genomic landscape of macrodissected LUAD and LUSC components of three ASC using whole exome sequencing (WES). Identified truncal mutations included the pan-cancer tumor-suppressor gene TP53 but also EGFR, BRAF, and MET, which are characteristic for LUAD but uncommon in LUSC. No truncal mutation of classical LUSC driver mutations were found. Both components showed unique driver mutations that did not overlap between the three ASC. Mutational signatures of truncal mutations differed from those of the branch mutations in their descendants LUAD and LUSC. Most common signatures were related to aging (1, 5) and smoking (4). Truncal chromosomal copy number aberrations shared by all three ASC included losses of 3p, 15q and 19p, and an amplified region in 5p. Furthermore, we detected loss of STK11 and SOX2 amplification in ASC, which has previously been shown to drive transdifferentiation from LUAD to LUSC in preclinical mouse models. Conclusively, this is the first study using WES to elucidate the clonal evolution of ASC. It provides strong evidence that the LUAD and LUSC components of ASC share a common origin and that the LUAD component appears to transform to LUSC.

Role of specialized composition of SWI/SNF complexes in prostate cancer lineage plasticity.

Advanced prostate cancer initially responds to hormonal treatment, but ultimately becomes resistant and requires more potent therapies. One mechanism of resistance observed in around 10-20% of these patients is lineage plasticity, which manifests in a partial or complete small cell or neuroendocrine prostate cancer (NEPC) phenotype. Here, we investigate the role of the mammalian SWI/SNF (mSWI/SNF) chromatin remodeling complex in NEPC. Using large patient datasets, patient-derived organoids and cancer cell lines, we identify mSWI/SNF subunits that are deregulated in NEPC and demonstrate that SMARCA4 (BRG1) overexpression is associated with aggressive disease. We also show that SWI/SNF complexes interact with different lineage-specific factors in NEPC compared to prostate adenocarcinoma. These data point to a role for mSWI/SNF complexes in therapy-related lineage plasticity, which may also be relevant for other solid tumors.

MTBseq: a comprehensive pipeline for whole genome sequence analysis of Mycobacterium tuberculosis complex isolates.

Analyzing whole-genome sequencing data of Mycobacterium tuberculosis complex (MTBC) isolates in a standardized workflow enables both comprehensive antibiotic resistance profiling and outbreak surveillance with highest resolution up to the identification of recent transmission chains. Here, we present MTBseq, a bioinformatics pipeline for next-generation genome sequence data analysis of MTBC isolates. Employing a reference mapping based workflow, MTBseq reports detected variant positions annotated with known association to antibiotic resistance and performs a lineage classification based on phylogenetic single nucleotide polymorphisms (SNPs). When comparing multiple datasets, MTBseq provides a joint list of variants and a FASTA alignment of SNP positions for use in phylogenomic analysis, and identifies groups of related isolates. The pipeline is customizable, expandable and can be used on a desktop computer or laptop without any internet connection, ensuring mobile usage and data security. MTBseq and accompanying documentation is available from https://github.com/ngs-fzb/MTBseq_source.

Culture and Next-generation sequencing-based drug susceptibility testing unveil high levels of drug-resistant-TB in Djibouti: results from the first national survey.

Djibouti is a small country in the Horn of Africa with a high TB incidence (378/100,000 in 2015). Multidrug-resistant TB (MDR-TB) and resistance to second-line agents have been previously identified in the country but the extent of the problem has yet to be quantified. A national survey was conducted to estimate the proportion of MDR-TB among a representative sample of TB patients. Sputum was tested using XpertMTB/RIF and samples positive for MTB and resistant to rifampicin underwent first line phenotypic susceptibility testing. The TB supranational reference laboratory in Milan, Italy, undertook external quality assurance, genotypic testing based on whole genome and targeted-deep sequencing and phylogenetic studies. 301 new and 66 previously treated TB cases were enrolled. MDR-TB was detected in 34 patients: 4.7% of new and 31% of previously treated cases. Resistance to pyrazinamide, aminoglycosides and capreomycin was detected in 68%, 18% and 29% of MDR-TB strains respectively, while resistance to fluoroquinolones was not detected. Cluster analysis identified transmission of MDR-TB as a critical factor fostering drug resistance in the country. Levels of MDR-TB in Djibouti are among the highest on the African continent. High prevalence of resistance to pyrazinamide and second-line injectable agents have important implications for treatment regimens.

Use of WGS in Mycobacterium tuberculosis routine diagnosis.

OBJECTIVE/BACKGROUND: Whole Genome Sequencing (WGS) is becoming affordable with overall costs comparable to other tests currently in use to perform the diagnosis of drug-resistant tuberculosis (TB) and cluster analysis. The WGS approach allows an "all-in-one" approach providing results on expected sensitivity of the strains, genetic background, epidemiological data, and indication of risk of laboratory cross-contamination. METHODS: Although ideal, WGS from the direct diagnostic specimen is not yet standardized, and to date the two most promising approaches are WGS from early positive liquid culture and targeted sequencing from diagnostic specimens using Next-Generation Technology. Both have advantages and disadvantages. Sequencing from early MGIT requires positive cultures, whereas targeted sequencing can be performed from a specimen positive for Mycobacterium tuberculosis with a consistent gain in time to information. The aim of this study is to evaluate the feasibility and cost of using WGS with a centralized approach to speed up diagnosis of TB in a low-incidence country. METHODS: From March 2016 to September 2016, we collected and processed by WGS 89 early positive routine MGIT960 tubes. Time to diagnosis and accuracy of this technique were compared with those of standard testing performed in a regular laboratory. A 2-mL aliquot of early positive MGIT was processed, starting with heat inactivation. DNA was then isolated by using the Maxwell 16 Cell DNA Purification Kit and Maxwell 16 MDx for automated extraction. Paired-end libraries of read-length 75-151bp were prepared using the Nextera XT DNA Sample Preparation kit, and sequenced on Illumina Miseq/Miniseq platform (based on the 1st available run). Total variant calling was performed according to the pipeline of the Phyresse web-tool. The DNA isolation step required 30min for inactivation plus 30min for extraction. The concentration obtained ranged from 0.1 to 1ng/µL, suitable for library preparation. Samples were sequenced with a turnaround time of 24-48h. The percentage of reads mapped to the H37Rv reference genome was 83% on average. The mean read coverage was 65×. The main challenge was the presence of nonmycobacterial DNA contamination in a variable amount. Lineage detection was possible for all cases, and mutations associated with drug resistance to antitubercular drugs were examined. We observed high diagnostic accuracy for species identification and detection of full drug resistance profile compared to standard DST testing performed in MGIT. RESULTS: Two events of recent transmissions including respectively three and two patients were identified, and two laboratory cross-contaminations were investigated and confirmed based on the analysis. Time to availability of report was about 72h from MGIT positivity compared to up to 6-9weeks for XDR-TB diagnosis with standard testing. In addition to speed, the main advantages were the availability of a full prediction of resistance determinants for rifampicin-resistant cases, and the fast detection of potential cross-contaminations and clusters to guide epidemiological investigation and cross-border tracing. Cost analysis showed that the cost per strain was approximately €150 inclusive of staff cost, reagents, and machine cost. CONCLUSION: WGS is a rapid, cost-effective technique that promises to integrate and replace the other tests in routine laboratories for an accurate diagnosis of DR-TB, although it is suitable nowadays for cultured samples only.

Activating stimuli induce platelet microRNA modulation and proteome reorganisation.

Platelets carry megakaryocyte-derived mRNAs whose translation efficiency before and during activation is not known, although this can greatly affect platelet functions, both under basal conditions and in response to physiological and pathological stimuli, such as those involved in acute coronary syndromes. Aim of the present study was to determine whether changes in microRNA (miRNA) expression occur in response to activating stimuli and whether this affects activity and composition of platelet transcriptome and proteome. Purified platelet-rich plasmas from healthy volunteers were collected and activated with ADP, collagen, or thrombin receptor activating peptide. Transcriptome analysis by RNA-Seq revealed that platelet transcriptome remained largely unaffected within the first 2 hours of stimulation. In contrast, quantitative proteomics showed that almost half of > 700 proteins quantified were modulated under the same conditions. Global miRNA analysis indicated that reorganisation of platelet proteome occurring during activation reflected changes in mature miRNA expression, which therefore, appears to be the main driver of the observed discrepancy between transcriptome and proteome changes. Platelet functions significantly affected by modulated miRNAs include, among others, the integrin/cytoskeletal, coagulation and inflammatory-immune response pathways. These results demonstrate a significant reprogramming of the platelet miRNome during activation, with consequent significant changes in platelet proteome and provide for the first time substantial evidence that fine-tuning of resident mRNA translation by miRNAs is a key event in platelet pathophysiology.

Lack of pathogenic mutations in six patients with MMPSI.

Sequencing of the KCNT1, PLCB1, SCN1A and TBC1D24 loci was performed in six children with typical features of malignant migrating partial seizures of infancy (MMPSI), to verify the presence of potential disease-causing mutations, including those already reported to be associated with the disease. Sanger sequencing failed to identify in these genes the previously reported pathogenic mutations in these patients, while a comprehensive mutational scanning analysis of these four loci by targeted re-sequencing led to detection of both intronic and exonic new variants. Based on the current knowledge, the sequence variants identified here do not allow to predict functional phenotypes that might explain, at least in part, MMPSI symptoms.

Global transcriptome profiles of Italian Mediterranean buffalo embryos with normal and retarded growth.

The transcriptome profiles were compared for buffalo embryos with normal growth and embryos with retarded growth on Day 25 after mating. Embryos with retarded growth on Day 25 after mating have a reduced likelihood of undergoing attachment to the uterine endometrium and establishing a pregnancy. Italian Mediterranean buffaloes were mated by AI and on Day 25 underwent trans-rectal ultrasonography to ascertain embryo development. Embryos with an embryonic width (EW)>2.7 mm were classed as normal embryos and embryos with an EW<2.7 mm were classed as retarded embryos. Three buffaloes with embryos of the largest EW (3.7, 3.7 and 3.9 mm) and three buffaloes with embryos of the smallest EW (1.5, 1.6 and 1.9 mm) were slaughtered on Day 27 to recover embryos for transcriptome analysis using a bovine custom designed oligo array. A total of 1,047 transcripts were differentially expressed between embryos with normal growth and embryos with retarded growth. Retarded embryos showed 773/1,047 (74%) transcripts that were down-regulated and 274/1,047 (26%) transcripts that were up-regulated relative to normal embryos; in silico analyses focused on 680/1,047 (65%) of the differentially expressed transcripts. The most altered transcripts observed in retarded embryos were associated with membrane structure and function and with metabolic and homeostasis maintenance functions. Other notable functions altered in retarded embryos were developmental processes and in particular nervous system differentiation and function. Specific biochemical pathways such as the complement cascade and coagulation were also altered in retarded embryos. It was concluded from the findings that buffalo embryos with retarded growth on Day 25 after mating show altered gene expression compared with normal embryos, and some de-regulated functions are associated with attachment to the uterine endometrium.

Effects of oestrogen on microRNA expression in hormone-responsive breast cancer cells.

Oestrogen receptor alpha (ERα) is a ligand-dependent transcription factor that mediates oestrogen effects in hormone-responsive cells. Following oestrogenic activation, ERα directly regulates the transcription of target genes via DNA binding. MicroRNAs (miRNAs) represent a class of small noncoding RNAs that function as negative regulators of protein-coding gene expression. They are found aberrantly expressed or mutated in cancer, suggesting their crucial role as either oncogenes or tumour suppressor genes. Here, we analysed changes in miRNA expression in response to oestrogen in hormone-responsive breast cancer MCF-7 and ZR-75.1 cells by microarray-mediated expression profiling. This led to the identification of 172 miRNAs up- or down-regulated by ERα in response to 17ß-oestradiol, of which 52 are similarly regulated by the hormone in the two cell models investigated. To identify mechanisms by which ERα exerts its effects on oestrogen-responsive miRNA genes, the oestrogen-dependent miRNA expression profiles were integrated with global in vivo ERα binding site mapping in the genome by ChIP-Seq. In addition, data from miRNA and messenger RNA (mRNA) expression profiles obtained under identical experimental conditions were compared to identify relevant miRNA target transcripts. Results show that miRNAs modulated by ERα represent a novel genomic pathway to impact oestrogen-dependent processes that affect hormone-responsive breast cancer cell behaviour. MiRNome analysis in tumour tissues from breast cancer patients confirmed a strong association between expression of these small RNAs and clinical outcome of the disease, although this appears to involve only marginally the oestrogen-regulated miRNAs identified in this study.