Clonal Origin and Lineage Ambiguity in Mixed Neuroendocrine Carcinoma of the Uterine Cervix.

Small-cell neuroendocrine carcinoma (SCNEC) of the cervix is a rare disease characterized by a high incidence of mixed tumors with other types of cancer. The mechanism underlying this mixed phenotype is not well understood. This study established a panel of organoid lines from patients with SCNEC of the cervix and ultimately focused on one line, which retained a mixed tumor phenotype, both in vitro and in vivo. Histologically, both organoids and xenograft tumors showed distinct differentiation into either SCNEC or adenocarcinoma in some regions and ambiguous differentiation in others. Tracking single cells indicated the existence of cells with bipotential differentiation toward SCNEC and adenocarcinomas. Single-cell transcriptional analysis identified three distinct clusters: SCNEC-like, adenocarcinoma-like, and a cluster lacking specific differentiation markers. The expression of neuroendocrine markers was enriched in the SCNEC-like cluster but not exclusively. Human papillomavirus 18 E6 was enriched in the SCNEC-like cluster, which showed higher proliferation and lower levels of the p53 pathway. After treatment with anticancer drugs, the expression of adenocarcinoma markers increased, whereas that of SCNEC decreased. Using a reporter system for keratin 19 expression, changes in the differentiation of each cell were shown to be associated with the shift in differentiation induced by drug treatment. These data suggest that mixed SCNEC/cervical tumors have a clonal origin and are characterized by an ambiguous and flexible differentiation state.

LncRNA ZNNT1 induces p53 degradation by interfering with the interaction between p53 and the SART3-USP15 complex.

Mammalian genomes encode large number of long noncoding RNAs (lncRNAs) that play key roles in various biological processes, including proliferation, differentiation, and stem cell pluripotency. Recent studies have addressed that some lncRNAs are dysregulated in human cancers and may play crucial roles in tumor development and progression. Here, we show that the lncRNA ZNNT1 is required for the proliferation and tumorigenicity of colon cancer cells with wild-type p53. ZNNT1 knockdown leads to decreased ubiquitination and stabilization of p53 protein. Moreover, we demonstrate that ZNNT1 needs to interact with SART3 to destabilize p53 and to promote the proliferation and tumorigenicity of colon cancer cells. We further show that SART3 is associated with the ubiquitin-specific peptidase USP15 and that ZNNT1 may induce p53 destabilization by inhibiting this interaction. These results suggest that ZNNT1 interferes with the SART3-USP15 complex-mediated stabilization of p53 protein and thereby plays important roles in the proliferation and tumorigenicity of colon cancer cells. Our findings suggest that ZNNT1 may be a promising molecular target for the therapy of colon cancer.

The Number of Transcription Factors at an Enhancer Determines Switch-like Gene Expression.

NF-κB is a transcription factor that activates super enhancers (SEs) and typical enhancers (TEs) and triggers threshold and graded gene expression, respectively. However, the mechanisms by which NF-κB selectively participates in these enhancers remain unclear. Here we show using mouse primary B lymphocytes that SE activity simultaneously associates with chromatin opening and enriched NF-κB binding, resulting in a higher fold change and threshold expression upon B cell receptor (BCR) activation. The higher fold change results from longer DNA, whereas the threshold response is explained by synergy in DNA-NF-κB binding and is supported by the coexistence of PU.1 and NF-κB in a SE before cell stimulation. This model indicates that the pre-existing NF-κB functions as a seed and triggers its processive binding upon BCR activation. Our mathematical modeling of the single-cell transcriptome reveals an additional role for SEs in divergent clonal responses in B cells.

Translation-dependent unwinding of stem-loops by UPF1 licenses Regnase-1 to degrade inflammatory mRNAs.

Regnase-1-mediated mRNA decay (RMD), in which inflammatory mRNAs harboring specific stem-loop structures are degraded, is a critical part of proper immune homeostasis. Prior to initial translation, Regnase-1 associates with target stem-loops but does not carry out endoribonucleolytic cleavage. Single molecule imaging revealed that UPF1 is required to first unwind the stem-loops, thus licensing Regnase-1 to proceed with RNA degradation. Following translation, Regnase-1 physically associates with UPF1 using two distinct points of interaction: The Regnase-1 RNase domain binds to SMG1-phosphorylated residue T28 in UPF1; in addition, an intrinsically disordered segment in Regnase-1 binds to the UPF1 RecA domain, enhancing the helicase activity of UPF1. The SMG1-UPF1-Regnase-1 axis targets pioneer rounds of translation and is critical for rapid resolution of inflammation through restriction of the number of proteins translated by a given mRNA. Furthermore, small-molecule inhibition of SMG1 prevents RNA unwinding in dendritic cells, allowing post-transcriptional control of innate immune responses.

Tamoxifen mechanically reprograms the tumor microenvironment via HIF-1A and reduces cancer cell survival.

The tumor microenvironment is fundamental to cancer progression, and the influence of its mechanical properties is increasingly being appreciated. Tamoxifen has been used for many years to treat estrogen-positive breast cancer. Here we report that tamoxifen regulates the level and activity of collagen cross-linking and degradative enzymes, and hence the organization of the extracellular matrix, via a mechanism involving both the G protein-coupled estrogen receptor (GPER) and hypoxia-inducible factor-1 alpha (HIF-1A). We show that tamoxifen reduces HIF-1A levels by suppressing myosin-dependent contractility and matrix stiffness mechanosensing. Tamoxifen also downregulates hypoxia-regulated genes and increases vascularization in PDAC tissues. Our findings implicate the GPER/HIF-1A axis as a master regulator of peri-tumoral stromal remodeling and the fibrovascular tumor microenvironment and offer a paradigm shift for tamoxifen from a well-established drug in breast cancer hormonal therapy to an alternative candidate for stromal targeting strategies in PDAC and possibly other cancers.

Hunt for the tipping point during endocrine resistance process in breast cancer by dynamic network biomarkers.

Acquired drug resistance is the major reason why patients fail to respond to cancer therapies. It is a challenging task to determine the tipping point of endocrine resistance and detect the associated molecules. Derived from new systems biology theory, the dynamic network biomarker (DNB) method is designed to quantitatively identify the tipping point of a drastic system transition and can theoretically identify DNB genes that play key roles in acquiring drug resistance. We analyzed time-course mRNA sequence data generated from the tamoxifen-treated estrogen receptor (ER)-positive MCF-7 cell line, and identified the tipping point of endocrine resistance with its leading molecules. The results show that there is interplay between gene mutations and DNB genes, in which the accumulated mutations eventually affect the DNB genes that subsequently cause the change of transcriptional landscape, enabling full-blown drug resistance. Survival analyses based on clinical datasets validated that the DNB genes were associated with the poor survival of breast cancer patients. The results provided the detection for the pre-resistance state or early signs of endocrine resistance. Our predictive method may greatly benefit the scheduling of treatments for complex diseases in which patients are exposed to considerably different drugs and may become drug resistant.

High-Quality Overlapping Paired-End Reads for the Detection of A-to-I Editing on Small RNA.

Paired-end RNA sequencing (RNA-seq) is usually applied to the quantification of long transcripts such as messenger or long non-coding RNAs, in which case overlapping pairs are discarded. In contrast, RNA-seq on short RNAs (≤ 200 nt) is typically carried out in single-end mode, as the additional cost associated with paired-end would only translate into redundant sequence information. Here, we exploit paired-end sequencing of short RNAs as a strategy to filter out sequencing errors and apply this method to the identification of adenosine-to-inosine (A-to-I) RNA editing events on human precursor microRNA (pre-miRNA) and mature miRNA. Combined with RNA immunoprecipitation sequencing (RIP-seq) of A-to-I RNA editing enzymes, this method takes full advantage of deep sequencing technology to identify RNA editing sites with unprecedented resolution in terms of editing efficiency.

Base-pairing probability in the microRNA stem region affects the binding and editing specificity of human A-to-I editing enzymes ADAR1-p110 and ADAR2.

Adenosine deaminases acting on RNA (ADARs) catalyze the deamination of adenosine (A) to inosine (I). A-to-I RNA editing targets double-stranded RNA (dsRNA), and increases the complexity of gene regulation by modulating base pairing-dependent processes such as splicing, translation, and microRNA (miRNA)-mediated gene silencing. This study investigates the genome-wide binding preferences of the nuclear constitutive isoforms ADAR1-p110 and ADAR2 on human miRNA species by RNA immunoprecipitation of ADAR-bound small RNAs (RIP-seq). Our results suggest that secondary structure predicted by base-pairing probability in the mainly double-stranded region of a pre-miRNA or mature miRNA duplex may determine ADAR isoform preference for binding distinct subpopulations of miRNAs. Furthermore, we identify 31 unique editing sites with statistical significance, 19 sites of which are novel editing sites. Editing sites are enriched in the seed region responsible for target recognition by miRNAs, and isoform-specific nucleotide motifs in the immediate vicinity and opposite of editing sites are consistent with previous studies, and further reveal that ADAR2 may edit A/C bulges more frequently than ADAR1-p110 in the context of miRNA.

Transcriptionally inducible Pleckstrin homology-like domain, family A, member 1, attenuates ErbB receptor activity by inhibiting receptor oligomerization.

Feedback control is a key mechanism in signal transduction, intimately involved in regulating the outcome of the cellular response. Here, we report a novel mechanism by which PHLDA1, Pleckstrin homology-like domain, family A, member 1, negatively regulates ErbB receptor signaling by inhibition of receptor oligomerization. We have found that the ErbB3 ligand, heregulin, induces PHILDA1 expression in MCF-7 cells. Transcriptionally-induced PHLDA1 protein directly binds to ErbB3, whereas knockdown of PHLDA1 increases complex formation between ErbB3 and ErbB2. To provide insight into the mechanism for our time-course and single-cell experimental observations, we performed a systematic computational search of network topologies of the mathematical models based on receptor dimer-tetramer formation in the ErbB activation processes. Our results indicate that only a model in which PHLDA1 inhibits formation of both dimers and tetramer can explain the experimental data. Predictions made from this model were further validated by single-molecule imaging experiments. Our studies suggest a unique regulatory feature of PHLDA1 to inhibit the ErbB receptor oligomerization process and thereby control the activity of receptor signaling network.

Genome-scale regression analysis reveals a linear relationship for promoters and enhancers after combinatorial drug treatment.

MOTIVATION: Drug combination therapy for treatment of cancers and other multifactorial diseases has the potential of increasing the therapeutic effect, while reducing the likelihood of drug resistance. In order to reduce time and cost spent in comprehensive screens, methods are needed which can model additive effects of possible drug combinations. RESULTS: We here show that the transcriptional response to combinatorial drug treatment at promoters, as measured by single molecule CAGE technology, is accurately described by a linear combination of the responses of the individual drugs at a genome wide scale. We also find that the same linear relationship holds for transcription at enhancer elements. We conclude that the described approach is promising for eliciting the transcriptional response to multidrug treatment at promoters and enhancers in an unbiased genome wide way, which may minimize the need for exhaustive combinatorial screens. AVAILABILITY AND IMPLEMENTATION: The CAGE sequence data used in this study is available in the DDBJ Sequence Read Archive (http://trace.ddbj.nig.ac.jp/index_e.html), accession number DRP001113. CONTACT: xin.gao@kaust.edu.sa or erik.arner@riken.jp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

ASBEL-TCF3 complex is required for the tumorigenicity of colorectal cancer cells.

Wnt/ß-catenin signaling plays a key role in the tumorigenicity of colon cancer. Furthermore, it has been reported that lncRNAs are dysregulated in several steps of cancer development. Here we show that ß-catenin directly activates the transcription of the long noncoding RNA (lncRNA) ASBEL [antisense ncRNA in the ANA (Abundant in neuroepithelium area)/BTG3 (B-cell translocation gene 3) locus] and transcription factor 3 (TCF3), both of which are required for the survival and tumorigenicity of colorectal cancer cells. ASBEL interacts with and recruits TCF3 to the activating transcription factor 3 (ATF3) locus, where it represses the expression of ATF3. Furthermore, we demonstrate that ASBEL-TCF3-mediated down-regulation of ATF3 expression is required for the proliferation and tumorigenicity of colon tumor cells. ATF3, in turn, represses the expression of ASBEL Our results reveal a pathway involving an lncRNA and two transcription factors that plays a key role in Wnt/ß-catenin-mediated tumorigenesis. These results may provide insights into the variety of biological and pathological processes regulated by Wnt/ß-catenin signaling.

MYU, a Target lncRNA for Wnt/c-Myc Signaling, Mediates Induction of CDK6 to Promote Cell Cycle Progression.

Aberrant activation of Wnt/ß-catenin signaling is a major driving force in colon cancer. Wnt/ß-catenin signaling induces the expression of the transcription factor c-Myc, leading to cell proliferation and tumorigenesis. c-Myc regulates multiple biological processes through its ability to directly modulate gene expression. Here, we identify a direct target of c-Myc, termed MYU, and show that MYU is upregulated in most colon cancers and required for the tumorigenicity of colon cancer cells. Furthermore, we demonstrate that MYU associates with the RNA binding protein hnRNP-K to stabilize CDK6 expression and thereby promotes the G1-S transition of the cell cycle. These results suggest that the MYU/hnRNP-K/CDK6 pathway functions downstream of Wnt/c-Myc signaling and plays a critical role in the proliferation and tumorigenicity of colon cancer cells.

Stimulus-Dependent Inhibitor of Apoptosis Protein Expression Prolongs the Duration of B Cell Signalling.

Different dynamic behaviours of signalling activity can induce distinct biological responses in a variety of cells. However, the molecular mechanisms that determine the dynamics of kinase activities in immune cells are not well understood. In this study, we showed that the duration of both IκB kinase (IKK) and extracellular signal-regulated kinase (ERK) activities in B cell receptor (BCR)- and CD40-signalling pathways in B cells were regulated by transcriptional feedback loops. We conducted a time-course transcriptome analysis after BCR or CD40 stimulation and identified the following four candidate genes as feedback regulators for IKK and ERK: inhibitor of apoptosis protein (IAP), TNF alpha-induced protein 3, dual-specificity phosphatase 5, and sprouty homolog 2. Quantitative experiments and mathematical modelling suggested that IAP inhibition shortened the duration of IKK and ERK activity following both BCR and CD40 pathway stimulation, indicating a positive role for IAP in B cell signalling. Furthermore, transient kinase activities induced by IAP blockage reduced the levels of delayed expression genes. Together, our findings suggest that IKK and ERK activity durations can be fine-tuned by the coordinated regulation of positive and negative transcriptional feedback and that these network properties determine the biological output of B cells.

Long noncoding RNA UPAT promotes colon tumorigenesis by inhibiting degradation of UHRF1.

Many long noncoding RNAs (lncRNAs) are reported to be dysregulated in human cancers and play critical roles in tumor development and progression. Furthermore, it has been reported that many lncRNAs regulate gene expression by recruiting chromatin remodeling complexes to specific genomic loci or by controlling transcriptional or posttranscriptional processes. Here we show that an lncRNA termed UPAT [ubiquitin-like plant homeodomain (PHD) and really interesting new gene (RING) finger domain-containing protein 1 (UHRF1) Protein Associated Transcript] is required for the survival and tumorigenicity of colorectal cancer cells. UPAT interacts with and stabilizes the epigenetic factor UHRF1 by interfering with its ß-transducin repeat-containing protein (TrCP)-mediated ubiquitination. Furthermore, we demonstrate that UHRF1 up-regulates Stearoyl-CoA desaturase 1 and Sprouty 4, which are required for the survival of colon tumor cells. Our study provides evidence for an lncRNA that regulates protein ubiquitination and degradation and thereby plays a critical role in the survival and tumorigenicity of tumor cells. Our results suggest that UPAT and UHRF1 may be promising molecular targets for the therapy of colon cancer.

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.

The transcription factor ATF7 mediates lipopolysaccharide-induced epigenetic changes in macrophages involved in innate immunological memory.

Immunological memory is thought to be mediated exclusively by lymphocytes. However, enhanced innate immune responses caused by a previous infection increase protection against reinfection, which suggests the presence of innate immunological memory. Here we identified an important role for the stress-response transcription factor ATF7 in innate immunological memory. ATF7 suppressed a group of genes encoding factors involved in innate immunity in macrophages by recruiting the histone H3K9 dimethyltransferase G9a. Treatment with lipopolysaccharide, which mimics bacterial infection, induced phosphorylation of ATF7 via the kinase p38, which led to the release of ATF7 from chromatin and a decrease in repressive histone H3K9me2 marks. A partially disrupted chromatin structure and increased basal expression of target genes were maintained for long periods, which enhanced resistance to pathogens. ATF7 might therefore be important in controlling memory in cells of the innate immune system.

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.

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.

Signalling mechanisms regulating phenotypic changes in breast cancer cells.

In MCF-7 breast cancer cells epidermal growth factor (EGF) induces cell proliferation, whereas heregulin (HRG)/neuregulin (NRG) induces irreversible phenotypic changes accompanied by lipid accumulation. Although these changes in breast cancer cells resemble processes that take place in the tissue, there is no understanding of signalling mechanisms regulating it. To identify molecular mechanisms mediating this cell-fate decision process, we applied different perturbations to pathways activated by these growth factors. The results demonstrate that phosphoinositide 3 (PI3) kinase (PI3K) and mammalian target of rapamycin (mTOR) complex (mTORC)1 activation is necessary for lipid accumulation that can also be induced by insulin, whereas stimulation of the extracellular-signal-regulated kinase (ERK) pathway is surprisingly dispensable. Interestingly, insulin exposure, as short as 4 h, was sufficient for triggering the lipid accumulation, whereas much longer treatment with HRG was required for achieving similar cellular response. Further, activation patterns of ATP citrate lyase (ACLY), an enzyme playing a central role in linking glycolytic and lipogenic pathways, suggest that lipids accumulated within cells are produced de novo rather than absorbed from the environment. In the present study, we demonstrate that PI3K pathway regulates phenotypic changes in breast cancer cells, whereas signal intensity and duration is crucial for cell fate decisions and commitment. Our findings reveal that MCF-7 cell fate decisions are controlled by a network of positive and negative regulators of both signalling and metabolic pathways.