Francesca Rizzo

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.

Two ParaHox genes, SpLox and SpCdx, interact to partition the posterior endoderm in the formation of a functional gut

We report the characterization of the ortholog of the Xenopus XlHbox8 ParaHox gene from the sea urchin Strongylocentrotus purpuratus, SpLox. It is expressed during embryogenesis, first appearing at late gastrulation in the posterior-most region of the endodermal tube, becoming progressively restricted to the constriction between the mid- and hindgut. The physiological effects of the absence of the activity of this gene have been analyzed through knockdown experiments using gene-specific morpholino antisense oligonucleotides. We show that blocking the translation of the SpLox mRNA reduces the capacity of the digestive tract to process food, as well as eliminating the morphological constriction normally present between the mid- and hindgut. Genetic interactions of the SpLox gene are revealed by the analysis of the expression of a set of genes involved in endoderm specification. Two such interactions have been analyzed in more detail: one involving the midgut marker gene Endo16, and another involving the other endodermally expressed ParaHox gene, SpCdx. We find that SpLox is able to bind Endo16 cis-regulatory DNA, suggesting direct repression of Endo16 expression in presumptive hindgut territories. More significantly, we provide the first evidence of interaction between ParaHox genes in establishing hindgut identity, and present a model of gene regulation involving a negative-feedback loop.

Selection and validation of a set of reliable reference genes for quantitative RT-PCR studies in the brain of the Cephalopod Mollusc Octopus vulgaris.

BackgroundQuantitative real-time polymerase chain reaction (RT-qPCR) is valuable for studying the molecular events underlying physiological and behavioral phenomena. Normalization of real-time PCR data is critical for a reliable mRNA quantification. Here we identify reference genes to be utilized in RT-qPCR experiments to normalize and monitor the expression of target genes in the brain of the cephalopod mollusc Octopus vulgaris, an invertebrate. Such an approach is novel for this taxon and of advantage in future experiments given the complexity of the behavioral repertoire of this species when compared with its relatively simple neural organization.ResultsWe chose 16S, and 18S rRNA, actB, EEF1A, tubA and ubi as candidate reference genes (housekeeping genes, HKG). The expression of 16S and 18S was highly variable and did not meet the requirements of candidate HKG. The expression of the other genes was almost stable and uniform among samples. We analyzed the expression of HKG into two different set of animals using tissues taken from the central nervous system (brain parts) and mantle (here considered as control tissue) by BestKeeper, geNorm and NormFinder. We found that HKG expressions differed considerably with respect to brain area and octopus samples in an HKG-specific manner. However, when the mantle is treated as control tissue and the entire central nervous system is considered, NormFinder revealed tubA and ubi as the most suitable HKG pair. These two genes were utilized to evaluate the relative expression of the genes FoxP, creb, dat and TH in O. vulgaris.ConclusionWe analyzed the expression profiles of some genes here identified for O. vulgaris by applying RT-qPCR analysis for the first time in cephalopods. We validated candidate reference genes and found the expression of ubi and tubA to be the most appropriate to evaluate the expression of target genes in the brain of different octopuses. Our results also underline the importance of choosing a proper normalization strategy when analyzing gene expression by qPCR taking into appropriate account the experimental setting and variability of the sample of animals (and tissues), thus providing a set of HGK which expression appears to be unaffected by the experimental factor(s).

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.

Characterization of two de novoKCNT1 mutations in children with malignant migrating partial seizures in infancy.

The KCNT1 gene encodes for subunits contributing to the Na(+)-activated K(+) current (KNa), expressed in many cell types. Mutations in KCNT1 have been found in patients affected with a wide spectrum of early-onset epilepsies, including Malignant Migrating Partial Seizures in Infancy (MMPSI), a severe early-onset epileptic encephalopathy characterized by pharmacoresistant focal seizures migrating from one brain region or hemisphere to another and neurodevelopment arrest or regression, resulting in profound disability. In the present study we report identification by whole exome sequencing (WES) of two de novo, heterozygous KCNT1 mutations (G288S and, not previously reported, M516V) in two unrelated MMPSI probands. Functional studies in a heterologous expression system revealed that channels formed by mutant KCNT1 subunits carried larger currents when compared to wild-type KCNT1 channels, both as homo- and heteromers with these last. Both mutations induced a marked leftward shift in homomeric channel activation gating. Interestingly, the KCNT1 blockers quinidine (3-1000μM) and bepridil (0.03-10μM) inhibited both wild-type and mutant KCNT1 currents in a concentration-dependent manner, with mutant channels showing higher sensitivity to blockade. This latter result suggests two genotype-tailored pharmacological strategies to specifically counteract the dysfunction of KCNT1 activating mutations in MMPSI patients.

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.

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.