Shigeyuki Magi

Transcriptional dynamics reveal critical roles for non-coding RNAs in the immediate-early response

The immediate-early response mediates cell fate in response to a variety of extracellular stimuli and is dysregulated in many cancers. However, the specificity of the response across stimuli and cell types, and the roles of non-coding RNAs are not well understood. Using a large collection of densely-sampled time series expression data we have examined the induction of the immediate-early response in unparalleled detail, across cell types and stimuli. We exploit cap analysis of gene expression (CAGE) time series datasets to directly measure promoter activities over time. Using a novel analysis method for time series data we identify transcripts with expression patterns that closely resemble the dynamics of known immediate-early genes (IEGs) and this enables a comprehensive comparative study of these genes and their chromatin state. Surprisingly, these data suggest that the earliest transcriptional responses often involve promoters generating non-coding RNAs, many of which are produced in advance of canonical protein-coding IEGs. IEGs are known to be capable of induction without de novo protein synthesis. Consistent with this, we find that the response of both protein-coding and non-coding RNA IEGs can be explained by their transcriptionally poised, permissive chromatin state prior to stimulation. We also explore the function of non-coding RNAs in the attenuation of the immediate early response in a small RNA sequencing dataset matched to the CAGE data: We identify a novel set of microRNAs responsible for the attenuation of the IEG response in an estrogen receptor positive cancer cell line. Our computational statistical method is well suited to meta-analyses as there is no requirement for transcripts to pass thresholds for significant differential expression between time points, and it is agnostic to the number of time points per dataset.

5‐Lipoxygenase and cysteinyl leukotriene receptor 1 regulate epidermal growth factor‐induced cell migration through Tiam1 upregulation and Rac1 activation

Cell migration is an essential step for tumor metastasis. The small GTPase Rac1 plays an important role in cell migration. Previously, we reported that epidermal growth factor (EGF) induced two waves of Rac1 activation; namely, at 5 min and 12 h after stimulation. A second wave of EGF‐induced Rac1 activation was required for EGF‐induced cell migration, however, the spatiotemporal regulation of the second wave of EGF‐induced Rac1 activation remains largely unclear. In this study, we found that 5‐lipoxygenase (5‐LOX) is activated in the process of EGF‐induced cell migration, and that leukotriene C4 (LTC4) produced by 5‐LOX mediated the second wave of Rac1 activation, as well as cell migration. Furthermore, these effects caused by LTC4 were found to be blocked in the presence of the antagonist of cysteinyl leukotriene receptor 1 (CysLT1). This blockage indicates that LTC4‐mediated CysLT1 signaling regulates the second EGF‐induced wave of Rac1 activation. We also found that 5‐LOX inhibitors, CysLT1 antagonists and the knockdown of CysLT1 inhibited EGF‐induced T cell lymphoma invasion and metastasis‐inducing protein 1 (Tiam1) expression. Tiam1 expression is required for the second wave of EGF‐induced Rac1 activation in A431 cells. Therefore, our results indicate that the 5‐LOX/LTC4/CysLT1 signaling pathway regulates EGF‐induced cell migration by increasing Tiam1 expression, leading to a second wave of Rac1 activation. Thus, CysLT1 may serve as a new molecular target for antimetastatic therapy. In addition, the CysLT1 antagonist, montelukast, which is used clinically for allergy treatment, might have great potential as a novel type of antimetastatic agent.

Application of Gene Expression Trajectories Initiated from ErbB Receptor Activation Highlights the Dynamics of Divergent Promoter Usage

Understanding how cells use complex transcriptional programs to alter their fate in response to specific stimuli is an important question in biology. For the MCF-7 human breast cancer cell line, we applied gene expression trajectory models to identify the genes involved in driving cell fate transitions. We modified trajectory models to account for the scenario where cells were exposed to different stimuli, in this case epidermal growth factor and heregulin, to arrive at different cell fates, i.e. proliferation and differentiation respectively. Using genome-wide CAGE time series data collected from the FANTOM5 consortium, we identified the sets of promoters that were involved in the transition of MCF-7 cells to their specific fates versus those with expression changes that were generic to both stimuli. Of the 1,552 promoters identified, 1,091 had stimulus-specific expression while 461 promoters had generic expression profiles over the time course surveyed. Many of these stimulus-specific promoters mapped to key regulators of the ERK (extracellular signal-regulated kinases) signaling pathway such as FHL2 (four and a half LIM domains 2). We observed that in general, generic promoters peaked in their expression early on in the time course, while stimulus-specific promoters tended to show activation of their expression at a later stage. The genes that mapped to stimulus-specific promoters were enriched for pathways that control focal adhesion, p53 signaling and MAPK signaling while generic promoters were enriched for cell death, transcription and the cell cycle. We identified 162 genes that were controlled by an alternative promoter during the time course where a subset of 37 genes had separate promoters that were classified as stimulus-specific and generic. The results of our study highlighted the degree of complexity involved in regulating a cell fate transition where multiple promoters mapping to the same gene can demonstrate quite divergent expression profiles.

Promoter-level expression clustering identifies time development of transcriptional regulatory cascades initiated by ErbB receptors in breast cancer cells

The analysis of CAGE (Cap Analysis of Gene Expression) time-course has been proposed by the FANTOM5 Consortium to extend the understanding of the sequence of events facilitating cell state transition at the level of promoter regulation. To identify the most prominent transcriptional regulations induced by growth factors in human breast cancer, we apply here the Complexity Invariant Dynamic Time Warping motif EnRichment (CIDER) analysis approach to the CAGE time-course datasets of MCF-7 cells stimulated by epidermal growth factor (EGF) or heregulin (HRG). We identify a multi-level cascade of regulations rooted by the Serum Response Factor (SRF) transcription factor, connecting the MAPK-mediated transduction of the HRG stimulus to the negative regulation of the MAPK pathway by the members of the DUSP family phosphatases. The finding confirms the known primary role of FOS and FOSL1, members of AP-1 family, in shaping gene expression in response to HRG induction. Moreover, we identify a new potential regulation of DUSP5 and RARA (known to antagonize the transcriptional regulation induced by the estrogen receptors) by the activity of the AP-1 complex, specific to HRG response. The results indicate that a divergence in AP-1 regulation determines cellular changes of breast cancer cells stimulated by ErbB receptors.

A chemical genomic study identifying diversity in cell migration signaling in cancer cells.

The aim of this study was to analyze the diversity and consistency of regulatory signaling in cancer cell migration, using a chemical genomic approach. The effects of 34 small molecular compounds were assessed quantitatively by wound healing assay in ten types of migrating cells. Hierarchical clustering was performed on the subsequent migration inhibition profile of the compounds and cancer cell types. The result was that hierarchical clustering accurately classified the compounds according to their targets. Furthermore, the cancer cells tested in this study were classified into three clusters, and the compounds were grouped into four clusters. An inhibitor of JNK suppressed all types of cell migration; however, inhibitors of ROCK, GSK-3 and p38MAPK only inhibited the migration of a subset of cell lines. Thus, our analytical system could easily distinguish between the common and cell type-specific signals responsible for cell migration.

Feedforward regulation of mRNA stability by prolonged extracellular signal‐regulated kinase activity

Extracellular signal‐regulated kinase (ERK) plays a central role in signal transduction networks and cell fate decisions. Sustained ERK activation induces cell differentiation, whereas transient ERK results in the proliferation of several types of cells. Sustained ERK activity stabilizes the proteins of early‐response gene products. However, the effect of ERK activity duration on mRNA stability is unknown. We analyzed the quantitative relationship between the duration of four ERK activity kinetics and the mRNA expression profile in growth factor‐treated cells. Time‐course transcriptome analysis revealed that the cells with prolonged ERK activity generally showed sustained mRNA expression of late response genes but not early or mid genes. Selected late response genes decayed more rapidly in the presence of a specific ERK inhibitor than a general transcription inhibitor and the decay rate was not related to the number of AU‐rich elements. Our results suggest that sustained ERK activity plays an important role in the lifespan of the mRNA encoded by late response genes, in addition to the previously demonstrated role in protein stabilization of early‐response genes, including transcription factors regulating the transcription of mid and late genes. This double‐positive regulation of ligand‐induced genes, also termed feedforward regulation, is critical in cell fate decisions.

Chemical genomic-based pathway analyses for epidermal growth factor-mediated signaling in migrating cancer cells.

To explore the diversity and consistency of the signaling pathways that regulate tumor cell migration, we chose three human tumor cell lines that migrated after treatment with EGF. We then quantified the effect of fifteen inhibitors on the levels of expression or the phosphorylation levels of nine proteins that were induced by EGF stimulation in each of these cell lines. Based on the data obtained in this study and chemical-biological assumptions, we deduced cell migration pathways in each tumor cell line, and then compared them. As a result, we found that both the MEK/ERK and JNK/c-Jun pathways were activated in all three migrating cell lines. Moreover, GSK-3 and p38 were found to regulate PI3K/Akt pathway in only EC109 cells, and JNK was found to crosstalk with p38 and Fos related pathway in only TT cells. Taken together, our analytical system could easily distinguish between the common and cell type-specific pathways responsible for tumor cell migration.

Novel derivatives of aclacinomycin A block cancer cell migration through inhibition of farnesyl transferase

In the course of screening for an inhibitor of farnesyl transferase (FTase), we identified two compounds, N-benzyl-aclacinomycin A (ACM) and N-allyl-ACM, which are new derivatives of ACM. N-benzyl-ACM and N-allyl-ACM inhibited FTase activity with IC50 values of 0.86 and 2.93 μM, respectively. Not only ACM but also C-10 epimers of each ACM derivative failed to inhibit FTase. The inhibition of FTase by N-benzyl-ACM and N-allyl-ACM seems to be specific, because these two compounds did not inhibit geranylgeranyltransferase or geranylgeranyl pyrophosphate (GGPP) synthase up to 100 μM. In cultured A431 cells, N-benzyl-ACM and N-allyl-ACM also blocked both the membrane localization of H-Ras and activation of the H-Ras-dependent PI3K/Akt pathway. In addition, they inhibited epidermal growth factor (EGF)-induced migration of A431 cells. Thus, N-benzyl-ACM and N-allyl-ACM inhibited EGF-induced migration of A431 cells by inhibiting the farnesylation of H-Ras and subsequent H-Ras-dependent activation of the PI3K/Akt pathway.

Robustness analysis of the detailed kinetic model of an ErbB signaling network by using dynamic sensitivity

The ErbB receptor signaling pathway plays an important role in the regulation of cellular proliferation, survival and differentiation, and dysregulation of the pathway is linked to various types of human cancer. Mathematical models have been developed as a practical complementary approach to deciphering the complexity of ErbB receptor signaling and elucidating how the pathways discriminate between ligands to induce different cell fates. In this study, we developed a simulator to accurately calculate the dynamic sensitivity of extracellular-signal-regulated kinase (ERK) activity (ERK*) and Akt activity (Akt*), downstream of the ErbB receptors stimulated with epidermal growth factor (EGF) and heregulin (HRG). To demonstrate the feasibility of this simulator, we estimated how the reactions critically responsible for ERK* and Akt* change with time and in response to different doses of EGF and HRG, and predicted that only a small number of reactions determine ERK* and Akt*. ERK* increased steeply with increasing HRG dose until saturation, while showing a gently rising response to EGF. Akt* had a gradual wide-range response to HRG and a blunt response to EGF. Akt* was sensitive to perturbations of intracellular kinetics, while ERK* was more robust due to multiple, negative feedback loops. Overall, the simulator predicted reactions that were critically responsible for ERK* and Akt* in response to the dose of EGF and HRG, illustrated the response characteristics of ERK* and Akt*, and estimated mechanisms for generating robustness in the ErbB signaling network.

Chemical genomic-based pathway analyses for epidermal growth factor-mediated signaling in migrating cancer cells (PLoS ONE (2014) 9, 5, (e96776) DOI: 10.1371/journal.pone.0096776)

There are errors in Figure 2. The top label indicating EGF treatments, positive or negative, are inverse. Please see the corrected Figure 2 here. Copyright: ß 2014 The PLOS ONE Staff. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.