Giorgio Giurato

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

BackgroundEstrogen 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.ResultsExpression 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.ConclusionsResults 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.

Direct regulation of microRNA biogenesis and expression by estrogen receptor beta in hormone-responsive breast cancer

Estrogen effects on mammary epithelial and breast cancer (BC) cells are mediated by the nuclear receptors ERα and ERβ, transcription factors that display functional antagonism with each other, with ERβ acting as oncosuppressor and interfering with the effects of ERα on cell proliferation, tumor promotion and progression. Indeed, hormone-responsive, ERα+ BC cells often lack ERβ, which when present associates with a less aggressive clinical phenotype of the disease. Recent evidences point to a significant role of microRNAs (miRNAs) in BC, where specific miRNA expression profiles associate with distinct clinical and biological phenotypes of the lesion. Considering the possibility that ERβ might influence BC cell behavior via miRNAs, we compared miRNome expression in ERβ+ vs ERβ− hormone-responsive BC cells and found a widespread effect of this ER subtype on the expression pattern of these non-coding RNAs. More importantly, the expression pattern of 67 miRNAs, including 10 regulated by ERβ in BC cells, clearly distinguishes ERβ+, node-negative, from ERβ−, metastatic, mammary tumors. Molecular dissection of miRNA biogenesis revealed multiple mechanisms for direct regulation of this process by ERβ+ in BC cell nuclei. In particular, ERβ downregulates miR-30a by binding to two specific sites proximal to the gene and thereby inhibiting pri-miR synthesis. On the other hand, the receptor promotes miR-23b, -27b and 24-1 accumulation in the cell by binding in close proximity of the corresponding gene cluster and preventing in situ the inhibitory effects of ERα on pri-miR maturation by the p68/DDX5-Drosha microprocessor complex. These results indicate that cell autonomous regulation of miRNA expression is part of the mechanism of action of ERβ in BC cells and could contribute to establishment or maintenance of a less aggressive tumor phenotype mediated by this nuclear receptor.

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.

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

BackgroundQualitative 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.ResultsWe 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.ConclusionThe 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.

Timed regulation of P‐element‐induced wimpy testis–interacting RNA expression during rat liver regeneration

Small noncoding RNAs comprise a growing family of molecules that regulate key cellular processes, including messenger RNA (mRNA) degradation, translational repression, and transcriptional gene silencing. P‐element‐induced wimpy testis (PIWI)‐interacting RNAs (piRNAs) represent a class of small RNAs initially identified in the germline of a variety of species, where they contribute to maintenance of genome stability, and recently found expressed also in stem and somatic cells, where their role and responsiveness to physiopathological signals remain elusive. Here, we investigated piRNA expression in rat liver and its response to the stimuli exerted by regenerative proliferation of this organ. Quantitative polymerase chain reaction analysis identify in the liver the RNAs encoding PIWIL2/HILI, PIWIL4/HIWI2, and other components of the piRNA biogenesis pathways, suggesting that this is indeed functional. RNA sequencing before, during, and after the wave of cell proliferation that follows partial hepatectomy (PH) identified ∼1,400 mammalian germline piRNAs expressed in rat liver, including 72 showing timed changes in expression 24‐48 hours post‐PH, a timing that corresponds to cell transition through the S phase, returning to basal levels by 168 hours, when organ regeneration is completed and hepatocytes reach quiescence. Conclusion: The piRNA pathway is active in somatic cells of the liver and is subject to regulation during the pathophysiological process of organ regeneration, when these molecules are available to exert their regulatory functions on the cell genome and transcriptome, as demonstrated by the identification of several liver mRNAs representing candidate targets of these regulatory RNAs. (Hepatology 2014;60:798–806)

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.

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.

Identification of H ferritin-dependent and independent genes in K562 differentiating cells by targeted gene silencing and expression profiling

Ferritin is best known as the key molecule in intracellular iron storage, and is involved in several metabolic processes such as cell proliferation, differentiation and neoplastic transformation. We have recently demonstrated that the shRNA silencing of the ferritin heavy subunit (FHC) in a melanoma cell line is accompanied by a consistent modification of gene expression pattern leading to a reduced potential in terms of proliferation, invasiveness, and adhesion ability of the silenced cells. In this study we sought to define the repertoire of genes whose expression might be affected by FHC during the hemin-induced differentiation of the erythromyeloid cell line K562. To this aim, gene expression profiling was performed in four different sets of cells: i) wild type K562; ii) sh-RNA FHC-silenced K562; iii) hemin-treated wild-type K562; and iv) hemin-treated FHC-silenced K562. Statistical analysis of the gene expression data, performed by two-factor ANOVA, identified three distinct classes of transcripts: a) Class 1, including 657 mRNAs whose expression is modified exclusively during hemin-induced differentiation of K562 cells, independently from the FHC relative amounts; b) Class 2, containing a set of 70 mRNAs which are consistently modified by hemin and FHC-silencing; and c) Class 3, including 128 transcripts modified by FHC-silencing but not by hemin. Our data indicate that FHC may function as a modulator of gene expression during erythroid differentiation and add new findings to the knowledge of the complex gene network modulated during erythroid differentiation.

Estrogen receptor beta impacts hormone-induced alternative mRNA splicing in breast cancer cells

BackgroundEstrogens play an important role in breast cancer (BC) development and progression; when the two isoforms of the estrogen receptor (ERα and ERβ) are co-expressed each of them mediate specific effects of these hormones in BC cells. ERβ has been suggested to exert an antagonist role toward the oncogenic activities of ERα, and for this reason it is considered an oncosuppressor. As clinical evidence regarding a prognostic role for this receptor subtype in hormone-responsive BC is still limited and conflicting, more knowledge is required on the biological functions of ERβ in cancer cells. We have previously described the ERβ and ERα interactomes from BC cells, identifying specific and distinct patterns of protein interactions for the two receptors. In particular, we identified factors involved in mRNA splicing and maturation as important components of both ERα and ERβ pathways. Guided by these findings, here we performed RNA sequencing to investigate in depth the differences in the early transcriptional events and RNA splicing patterns induced by estradiol in cells expressing ERα alone or ERα and ERβ.ResultsExon skipping was the most abundant splicing event in the post-transcriptional regulation by estradiol. We identified several splicing events induced by ERα alone and by ERα + ERβ, demonstrating for the first time that ERβ significantly affects estrogen-induced splicing in BC cells, as revealed by modification of a subset of ERα-dependent splicing by ERβ, as well as by the presence of splicing isoforms only in ERβ + cells. In particular, we observed that ERβ + BC cell lines exhibited around 2-fold more splicing events than the ERβ- cells. Interestingly, we identified putative direct targets of ERβ-mediated alternative splicing by correlating the genomic locations of ERβ and ERα binding sites with estradiol-induced differential splicing in the corresponding genes.ConclusionsTaken together, these results demonstrate that ERβ significantly affects estrogen-induced early transcription and mRNA splicing in hormone-responsive BC cells, providing novel information on the biological role of ERβ in these tumors.

Specific patterns of PIWI-interacting small noncoding RNA expression in dysplastic liver nodules and hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the result of a stepwise process, often beginning with development within a cirrhotic liver of premalignant lesions, morphologically characterized by low- (LGDN) and high-grade (HGDN) dysplastic nodules. PIWI-interacting RNAs (piRNAs) are small noncoding RNAs (sncRNAs), 23–35 nucleotide-long, exerting epigenetic and post-transcriptional regulation of gene expression. Recently the PIWI-piRNA pathway, best characterized in germline cells, has been identified also in somatic tissues, including stem and cancer cells, where it influences key cellular processes. Small RNA sequencing was applied to search for liver piRNAs and to profile their expression patterns in cirrhotic nodules (CNs), LGDN, HGDN, early HCC and progressed HCC (pHCC), analyzing 55 samples (14 CN, 9 LGDN, 6 HGDN, 6 eHCC and 20 pHCC) from 17 patients, aiming at identifying possible relationships between these sncRNAs and liver carcinogenesis. We identified a 125 piRNA expression signature that characterize HCC from matched CNs, correlating also to microvascular invasion in HCC. Functional analysis of the predicted RNA targets of deregulated piRNAs indicates that these can target key signaling pathways involved in hepatocarcinogenesis and HCC progression, thereby affecting their activity. Interestingly, 24 piRNAs showed specific expression patterns in dysplastic nodules, respect to cirrhotic liver and/or pHCC. The results demonstrate that the PIWI-piRNA pathway is active in human liver, where it represents a new player in the molecular events that characterize hepatocarcinogenesis, from early stages to pHCC. Furthermore, they suggest that piRNAs might be new disease biomarkers, useful for differential diagnosis of dysplastic and neoplastic liver lesions.