MethReg: estimating the regulatory potential of DNA methylation in gene transcription.

Epigenome-wide association studies often detect many differentially methylated sites, and many are located in distal regulatory regions. To further prioritize these significant sites, there is a critical need to better understand the functional impact of CpG methylation. Recent studies demonstrated that CpG methylation-dependent transcriptional regulation is a widespread phenomenon. Here, we present MethReg, an R/Bioconductor package that analyzes matched DNA methylation and gene expression data, along with external transcription factor (TF) binding information, to evaluate, prioritize and annotate CpG sites with high regulatory potential. At these CpG sites, TF-target gene associations are often only present in a subset of samples with high (or low) methylation levels, so they can be missed by analyses that use all samples. Using colorectal cancer and Alzheimer's disease datasets, we show MethReg significantly enhances our understanding of the regulatory roles of DNA methylation in complex diseases.

Activation of bivalent factor DLX5 cooperates with master regulator TP63 to promote squamous cell carcinoma.

To reconstruct systematically hyperactive transcription factor (TF)-dependent transcription networks in squamous cell carcinomas (SCCs), a computational method (ELMER) was applied to 1293 pan-SCC patient samples, and 44 hyperactive SCC TFs were identified. As a top candidate, DLX5 exhibits a notable bifurcate re-configuration of its bivalent promoter in cancer. Specifically, DLX5 maintains a bivalent state in normal tissues; its promoter is hypermethylation, leading to DLX5 transcriptional silencing in esophageal adenocarcinoma (EAC). In stark contrast, DLX5 promoter gains active histone marks and becomes transcriptionally activated in ESCC, which is directly mediated by SOX2. Functionally, silencing of DLX5 substantially inhibits SCC viability both in vitro and in vivo. Mechanistically, DLX5 cooperates with TP63 in regulating ∼2000 enhancers and promoters, which converge on activating cancer-promoting pathways. Together, our data establish a novel and strong SCC-promoting factor and elucidate a new epigenomic mechanism - bifurcate chromatin re-configuration - during cancer development.

White paper on high-throughput process development for integrated continuous biomanufacturing.

Continuous manufacturing is an indicator of a maturing industry, as can be seen by the example of the petrochemical industry. Patent expiry promotes a price competition between manufacturing companies, and more efficient and cheaper processes are needed to achieve lower production costs. Over the last decade, continuous biomanufacturing has had significant breakthroughs, with regulatory agencies encouraging the industry to implement this processing mode. Process development is resource and time consuming and, although it is increasingly becoming less expensive and faster through high-throughput process development (HTPD) implementation, reliable HTPD technology for integrated and continuous biomanufacturing is still lacking and is considered to be an emerging field. Therefore, this paper aims to illustrate the major gaps in HTPD and to discuss the major needs and possible solutions to achieve an end-to-end Integrated Continuous Biomanufacturing, as discussed in the context of the 2019 Integrated Continuous Biomanufacturing conference. The current HTPD state-of-the-art for several unit operations is discussed, as well as the emerging technologies which will expedite a shift to continuous biomanufacturing.

ELMER v.2: an R/Bioconductor package to reconstruct gene regulatory networks from DNA methylation and transcriptome profiles.

MOTIVATION: DNA methylation has been used to identify functional changes at transcriptional enhancers and other cis-regulatory modules (CRMs) in tumors and other disease tissues. Our R/Bioconductor package ELMER (Enhancer Linking by Methylation/Expression Relationships) provides a systematic approach that reconstructs altered gene regulatory networks (GRNs) by combining enhancer methylation and gene expression data derived from the same sample set. RESULTS: We present a completely revised version 2 of ELMER that provides numerous new features including an optional web-based interface and a new Supervised Analysis mode to use pre-defined sample groupings. We show that Supervised mode significantly increases statistical power and identifies additional GRNs and associated Master Regulators, such as SOX11 and KLF5 in Basal-like breast cancer. AVAILABILITY AND IMPLEMENTATION: ELMER v.2 is available as an R/Bioconductor package at http://bioconductor.org/packages/ELMER/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

New functionalities in the TCGAbiolinks package for the study and integration of cancer data from GDC and GTEx.

The advent of Next-Generation Sequencing (NGS) technologies has opened new perspectives in deciphering the genetic mechanisms underlying complex diseases. Nowadays, the amount of genomic data is massive and substantial efforts and new tools are required to unveil the information hidden in the data. The Genomic Data Commons (GDC) Data Portal is a platform that contains different genomic studies including the ones from The Cancer Genome Atlas (TCGA) and the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiatives, accounting for more than 40 tumor types originating from nearly 30000 patients. Such platforms, although very attractive, must make sure the stored data are easily accessible and adequately harmonized. Moreover, they have the primary focus on the data storage in a unique place, and they do not provide a comprehensive toolkit for analyses and interpretation of the data. To fulfill this urgent need, comprehensive but easily accessible computational methods for integrative analyses of genomic data that do not renounce a robust statistical and theoretical framework are required. In this context, the R/Bioconductor package TCGAbiolinks was developed, offering a variety of bioinformatics functionalities. Here we introduce new features and enhancements of TCGAbiolinks in terms of i) more accurate and flexible pipelines for differential expression analyses, ii) different methods for tumor purity estimation and filtering, iii) integration of normal samples from other platforms iv) support for other genomics datasets, exemplified here by the TARGET data. Evidence has shown that accounting for tumor purity is essential in the study of tumorigenesis, as these factors promote confounding behavior regarding differential expression analysis. With this in mind, we implemented these filtering procedures in TCGAbiolinks. Moreover, a limitation of some of the TCGA datasets is the unavailability or paucity of corresponding normal samples. We thus integrated into TCGAbiolinks the possibility to use normal samples from the Genotype-Tissue Expression (GTEx) project, which is another large-scale repository cataloging gene expression from healthy individuals. The new functionalities are available in the TCGAbiolinks version 2.8 and higher released in Bioconductor version 3.7.

GENAVi: a shiny web application for gene expression normalization, analysis and visualization.

BACKGROUND: The development of next generation sequencing (NGS) methods led to a rapid rise in the generation of large genomic datasets, but the development of user-friendly tools to analyze and visualize these datasets has not developed at the same pace. This presents a two-fold challenge to biologists; the expertise to select an appropriate data analysis pipeline, and the need for bioinformatics or programming skills to apply this pipeline. The development of graphical user interface (GUI) applications hosted on web-based servers such as Shiny can make complex workflows accessible across operating systems and internet browsers to those without programming knowledge. RESULTS: We have developed GENAVi (Gene Expression Normalization Analysis and Visualization) to provide a user-friendly interface for normalization and differential expression analysis (DEA) of human or mouse feature count level RNA-Seq data. GENAVi is a GUI based tool that combines Bioconductor packages in a format for scientists without bioinformatics expertise. We provide a panel of 20 cell lines commonly used for the study of breast and ovarian cancer within GENAVi as a foundation for users to bring their own data to the application. Users can visualize expression across samples, cluster samples based on gene expression or correlation, calculate and plot the results of principal components analysis, perform DEA and gene set enrichment and produce plots for each of these analyses. To allow scalability for large datasets we have provided local install via three methods. We improve on available tools by offering a range of normalization methods and a simple to use interface that provides clear and complete session reporting and for reproducible analysis. CONCLUSION: The development of tools using a GUI makes them practical and accessible to scientists without bioinformatics expertise, or access to a data analyst with relevant skills. While several GUI based tools are currently available for RNA-Seq analysis we improve on these existing tools. This user-friendly application provides a convenient platform for the normalization, analysis and visualization of gene expression data for scientists without bioinformatics expertise.

TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data.

The Cancer Genome Atlas (TCGA) research network has made public a large collection of clinical and molecular phenotypes of more than 10 000 tumor patients across 33 different tumor types. Using this cohort, TCGA has published over 20 marker papers detailing the genomic and epigenomic alterations associated with these tumor types. Although many important discoveries have been made by TCGA's research network, opportunities still exist to implement novel methods, thereby elucidating new biological pathways and diagnostic markers. However, mining the TCGA data presents several bioinformatics challenges, such as data retrieval and integration with clinical data and other molecular data types (e.g. RNA and DNA methylation). We developed an R/Bioconductor package called TCGAbiolinks to address these challenges and offer bioinformatics solutions by using a guided workflow to allow users to query, download and perform integrative analyses of TCGA data. We combined methods from computer science and statistics into the pipeline and incorporated methodologies developed in previous TCGA marker studies and in our own group. Using four different TCGA tumor types (Kidney, Brain, Breast and Colon) as examples, we provide case studies to illustrate examples of reproducibility, integrative analysis and utilization of different Bioconductor packages to advance and accelerate novel discoveries.

SpidermiR: An R/Bioconductor Package for Integrative Analysis with miRNA Data.

Gene Regulatory Networks (GRNs) control many biological systems, but how such network coordination is shaped is still unknown. GRNs can be subdivided into basic connections that describe how the network members interact e.g., co-expression, physical interaction, co-localization, genetic influence, pathways, and shared protein domains. The important regulatory mechanisms of these networks involve miRNAs. We developed an R/Bioconductor package, namely SpidermiR, which offers an easy access to both GRNs and miRNAs to the end user, and integrates this information with differentially expressed genes obtained from The Cancer Genome Atlas. Specifically, SpidermiR allows the users to: (i) query and download GRNs and miRNAs from validated and predicted repositories; (ii) integrate miRNAs with GRNs in order to obtain miRNA-gene-gene and miRNA-protein-protein interactions, and to analyze miRNA GRNs in order to identify miRNA-gene communities; and (iii) graphically visualize the results of the analyses. These analyses can be performed through a single interface and without the need for any downloads. The full data sets are then rapidly integrated and processed locally.

Lamina Propria Phagocyte Profiling Reveals Targetable Signaling Pathways in Refractory Inflammatory Bowel Disease.

BACKGROUND AND AIMS: Lamina propria phagocytes are key mediators of inflammatory bowel disease (IBD). We aimed to understand the transcriptomic and functional differences in these cells based on location, disease type, inflammation state, and medication use in patients with IBD. METHODS: Phagocytic immune cells in the lamina propria, as defined by the marker CD11b, were isolated from 54 unique patients (n = 111 gut mucosal biopsies). We performed flow cytometry for cell phenotyping (n = 30) and RNA sequencing with differential gene expression analysis (n = 58). We further cultured these cells in vitro and exposed them to janus kinase inhibitors to measure cytokine output (n = 27). Finally, we matched patient genomic data to our RNA sequencing data to perform candidate gene expression quantitative trait locus analysis (n = 34). RESULTS: We found distinct differences in gene expression between CD11b+ cells from the colon vs ileum, as well as in different inflammatory states and, to a lesser degree, IBD types (Crohn's disease or ulcerative colitis). These genes mapped to targetable immune pathways and metabolic and cancer pathways. We further explored the janus kinase-signal transducer and activator of transcription pathway, which was upregulated across many comparisons including in biopsies from anti-tumor necrosis factor refractory patients. We found that isolated CD11b+ cells treated with janus kinase inhibitors had decreased secretion of cytokines tumor necrosis factora and interleukin-8 (P ≤ .05). We also found 3 genetic variants acting as expression quantitative trait loci (P ≤ .1) within our CD11b+ data set. CONCLUSIONS: Lamina propria phagocytes from IBD mucosa provide pathogenetic clues on the nature of treatment refractoriness and inform new targets for therapy.

DNA methylation and transcriptomic features are preserved throughout disease recurrence and chemoresistance in high grade serous ovarian cancers.

BACKGROUND: Little is known about the role of global DNA methylation in recurrence and chemoresistance of high grade serous ovarian cancer (HGSOC). METHODS: We performed whole genome bisulfite sequencing and transcriptome sequencing in 62 primary and recurrent tumors from 28 patients with stage III/IV HGSOC, of which 11 patients carried germline, pathogenic BRCA1 and/or BRCA2 mutations. RESULTS: Landscapes of genome-wide methylation (on average 24.2 million CpGs per tumor) and transcriptomes in primary and recurrent tumors showed extensive heterogeneity between patients but were highly preserved in tumors from the same patient. We identified significant differences in the burden of differentially methylated regions (DMRs) in tumors from BRCA1/2 compared to non-BRCA1/2 carriers (mean 659 DMRs and 388 DMRs in paired comparisons respectively). We identified overexpression of immune pathways in BRCA1/2 carriers compared to non-carriers, implicating an increased immune response in improved survival (P = 0.006) in these BRCA1/2 carriers. CONCLUSION: These findings indicate methylome and gene expression programs established in the primary tumor are conserved throughout disease progression, even after extensive chemotherapy treatment, and that changes in methylation and gene expression are unlikely to serve as drivers for chemoresistance in HGSOC.

PathwayMultiomics: An R Package for Efficient Integrative Analysis of Multi-Omics Datasets With Matched or Un-matched Samples.

Recent advances in technology have made multi-omics datasets increasingly available to researchers. To leverage the wealth of information in multi-omics data, a number of integrative analysis strategies have been proposed recently. However, effectively extracting biological insights from these large, complex datasets remains challenging. In particular, matched samples with multiple types of omics data measured on each sample are often required for multi-omics analysis tools, which can significantly reduce the sample size. Another challenge is that analysis techniques such as dimension reductions, which extract association signals in high dimensional datasets by estimating a few variables that explain most of the variations in the samples, are typically applied to whole-genome data, which can be computationally demanding. Here we present pathwayMultiomics, a pathway-based approach for integrative analysis of multi-omics data with categorical, continuous, or survival outcome variables. The input of pathwayMultiomics is pathway p-values for individual omics data types, which are then integrated using a novel statistic, the MiniMax statistic, to prioritize pathways dysregulated in multiple types of omics datasets. Importantly, pathwayMultiomics is computationally efficient and does not require matched samples in multi-omics data. We performed a comprehensive simulation study to show that pathwayMultiomics significantly outperformed currently available multi-omics tools with improved power and well-controlled false-positive rates. In addition, we also analyzed real multi-omics datasets to show that pathwayMultiomics was able to recover known biology by nominating biologically meaningful pathways in complex diseases such as Alzheimer's disease.

LR Hunting: A Random Forest Based Cell-Cell Interaction Discovery Method for Single-Cell Gene Expression Data.

Cell-cell interactions (CCIs) and cell-cell communication (CCC) are critical for maintaining complex biological systems. The availability of single-cell RNA sequencing (scRNA-seq) data opens new avenues for deciphering CCIs and CCCs through identifying ligand-receptor (LR) gene interactions between cells. However, most methods were developed to examine the LR interactions of individual pairs of genes. Here, we propose a novel approach named LR hunting which first uses random forests (RFs)-based data imputation technique to link the data between different cell types. To guarantee the robustness of the data imputation procedure, we repeat the computation procedures multiple times to generate aggregated imputed minimal depth index (IMDI). Next, we identify significant LR interactions among all combinations of LR pairs simultaneously using unsupervised RFs. We demonstrated LR hunting can recover biological meaningful CCIs using a mouse cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) dataset and a triple-negative breast cancer scRNA-seq dataset.

A pan-cancer analysis of CpG Island gene regulation reveals extensive plasticity within Polycomb target genes.

CpG Island promoter genes make up more than half of human genes, and a subset regulated by Polycomb-Repressive Complex 2 (PRC2+-CGI) become DNA hypermethylated and silenced in cancer. Here, we perform a systematic analysis of CGI genes across TCGA cancer types, finding that PRC2+-CGI genes are frequently prone to transcriptional upregulation as well. These upregulated PRC2+-CGI genes control important pathways such as Epithelial-Mesenchymal Transition (EMT) and TNFα-associated inflammatory response, and have greater cancer-type specificity than other CGI genes. Using publicly available chromatin datasets and genetic perturbations, we show that transcription factor binding sites (TFBSs) within distal enhancers underlie transcriptional activation of PRC2+-CGI genes, coinciding with loss of the PRC2-associated mark H3K27me3 at the linked promoter. In contrast, PRC2-free CGI genes are predominantly regulated by promoter TFBSs which are common to most cancer types. Surprisingly, a large subset of PRC2+-CGI genes that are upregulated in one cancer type are also hypermethylated/silenced in at least one other cancer type, underscoring the high degree of regulatory plasticity of these genes, likely derived from their complex regulatory control during normal development.

Sex-specific DNA methylation differences in Alzheimer's disease pathology.

Sex is an important factor that contributes to the clinical and biological heterogeneities in Alzheimer's disease (AD), but the regulatory mechanisms underlying sex disparity in AD are still not well understood. DNA methylation is an important epigenetic modification that regulates gene transcription and is known to be involved in AD. We performed the first large-scale sex-specific meta-analysis of DNA methylation differences in AD neuropathology, by re-analyzing four recent epigenome-wide association studies totaling more than 1000 postmortem prefrontal cortex brain samples using a uniform analytical pipeline. For each cohort, we employed two complementary analytical strategies, a sex-stratified analysis that examined methylation-Braak stage associations in male and female samples separately, and a sex-by-Braak stage interaction analysis that compared the magnitude of these associations between different sexes. Our analysis uncovered 14 novel CpGs, mapped to genes such as TMEM39A and TNXB that are associated with the AD Braak stage in a sex-specific manner. TMEM39A is known to be involved in inflammation, dysregulated type I interferon responses, and other immune processes. TNXB encodes tenascin proteins, which are extracellular matrix glycoproteins demonstrated to modulate synaptic plasticity in the brain. Moreover, for many previously implicated genes in AD neuropathology, such as MBP and AZU1, our analysis provided the new insights that they were predominately driven by effects in only one sex. These sex-specific DNA methylation differences were enriched in divergent biological processes such as integrin activation in females and complement activation in males. Our study implicated multiple new loci and biological processes that affected AD neuropathology in a sex-specific manner.

Multi-omics analysis identifies therapeutic vulnerabilities in triple-negative breast cancer subtypes.

Triple-negative breast cancer (TNBC) is a collection of biologically diverse cancers characterized by distinct transcriptional patterns, biology, and immune composition. TNBCs subtypes include two basal-like (BL1, BL2), a mesenchymal (M) and a luminal androgen receptor (LAR) subtype. Through a comprehensive analysis of mutation, copy number, transcriptomic, epigenetic, proteomic, and phospho-proteomic patterns we describe the genomic landscape of TNBC subtypes. Mesenchymal subtype tumors display high mutation loads, genomic instability, absence of immune cells, low PD-L1 expression, decreased global DNA methylation, and transcriptional repression of antigen presentation genes. We demonstrate that major histocompatibility complex I (MHC-I) is transcriptionally suppressed by H3K27me3 modifications by the polycomb repressor complex 2 (PRC2). Pharmacological inhibition of PRC2 subunits EZH2 or EED restores MHC-I expression and enhances chemotherapy efficacy in murine tumor models, providing a rationale for using PRC2 inhibitors in PD-L1 negative mesenchymal tumors. Subtype-specific differences in immune cell composition and differential genetic/pharmacological vulnerabilities suggest additional treatment strategies for TNBC.

DNA methylation-based signatures classify sporadic pituitary tumors according to clinicopathological features.

BACKGROUND: Distinct genome-wide methylation patterns cluster pituitary neuroendocrine tumors (PitNETs) into molecular groups associated with specific clinicopathological features. Here we aim to identify, characterize, and validate methylation signatures that objectively classify PitNET into clinicopathological groups. METHODS: Combining in-house and publicly available data, we conducted an analysis of the methylome profile of a comprehensive cohort of 177 tumors (Panpit cohort) and 20 nontumor specimens from the pituitary gland. We also retrieved methylome data from an independent PitNET cohort (N = 86) to validate our findings. RESULTS: We identified three methylation clusters associated with adenohypophyseal cell lineages and functional status using an unsupervised approach. Differentially methylated probes (DMP) significantly distinguished the Panpit clusters and accurately assigned the samples of the validation cohort to their corresponding lineage and functional subtypes memberships. The DMPs were annotated in regulatory regions enriched with enhancer elements, associated with pathways and genes involved in pituitary cell identity, function, tumorigenesis, and invasiveness. Some DMPs correlated with genes with prognostic and therapeutic values in other intra- or extracranial tumors. CONCLUSIONS: We identified and validated methylation signatures, mainly annotated in enhancer regions that distinguished PitNETs by distinct adenohypophyseal cell lineages and functional status. These signatures provide the groundwork to develop an unbiased approach to classifying PitNETs according to the most recent classification recommended by the 2017 WHO and to explore their biological and clinical relevance in these tumors.

Epigenome-wide meta-analysis of DNA methylation differences in prefrontal cortex implicates the immune processes in Alzheimer's disease.

DNA methylation differences in Alzheimer's disease (AD) have been reported. Here, we conducted a meta-analysis of more than 1000 prefrontal cortex brain samples to prioritize the most consistent methylation differences in multiple cohorts. Using a uniform analysis pipeline, we identified 3751 CpGs and 119 differentially methylated regions (DMRs) significantly associated with Braak stage. Our analysis identified differentially methylated genes such as MAMSTR, AGAP2, and AZU1. The most significant DMR identified is located on the MAMSTR gene, which encodes a cofactor that stimulates MEF2C. Notably, MEF2C cooperates with another transcription factor, PU.1, a central hub in the AD gene network. Our enrichment analysis highlighted the potential roles of the immune system and polycomb repressive complex 2 in pathological AD. These results may help facilitate future mechanistic and biomarker discovery studies in AD.

Interpreting pathways to discover cancer driver genes with Moonlight.

Cancer driver gene alterations influence cancer development, occurring in oncogenes, tumor suppressors, and dual role genes. Discovering dual role cancer genes is difficult because of their elusive context-dependent behavior. We define oncogenic mediators as genes controlling biological processes. With them, we classify cancer driver genes, unveiling their roles in cancer mechanisms. To this end, we present Moonlight, a tool that incorporates multiple -omics data to identify critical cancer driver genes. With Moonlight, we analyze 8000+ tumor samples from 18 cancer types, discovering 3310 oncogenic mediators, 151 having dual roles. By incorporating additional data (amplification, mutation, DNA methylation, chromatin accessibility), we reveal 1000+ cancer driver genes, corroborating known molecular mechanisms. Additionally, we confirm critical cancer driver genes by analysing cell-line datasets. We discover inactivation of tumor suppressors in intron regions and that tissue type and subtype indicate dual role status. These findings help explain tumor heterogeneity and could guide therapeutic decisions.

Lineage-Specific Epigenomic and Genomic Activation of Oncogene HNF4A Promotes Gastrointestinal Adenocarcinomas.

Gastrointestinal adenocarcinomas (GIAC) of the tubular gastrointestinal (GI) tract including esophagus, stomach, colon, and rectum comprise most GI cancers and share a spectrum of genomic features. However, the unified epigenomic changes specific to GIAC are poorly characterized. Using 907 GIAC samples from The Cancer Genome Atlas, we applied mathematical algorithms to large-scale DNA methylome and transcriptome profiles to reconstruct transcription factor (TF) networks and identify a list of functionally hyperactive master regulator (MR) TF shared across different GIAC. The top candidate HNF4A exhibited prominent genomic and epigenomic activation in a GIAC-specific manner. A complex interplay between the HNF4A promoter and three distal enhancer elements was coordinated by GIAC-specific MRTF including ELF3, GATA4, GATA6, and KLF5. HNF4A also self-regulated its own promoter and enhancers. Functionally, HNF4A promoted cancer proliferation and survival by transcriptional activation of many downstream targets, including HNF1A and factors of interleukin signaling, in a lineage-specific manner. Overall, our study provides new insights into the GIAC-specific gene regulatory networks and identifies potential therapeutic strategies against these common cancers. SIGNIFICANCE: These findings show that GIAC-specific master regulatory transcription factors control HNF4A via three distal enhancers to promote GIAC cell proliferation and survival. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/13/2722/F1.large.jpg.

RNA-Binding Protein ZFP36L1 Suppresses Hypoxia and Cell-Cycle Signaling.

ZFP36L1 is a tandem zinc-finger RNA-binding protein that recognizes conserved adenylate-uridylate-rich elements (ARE) located in 3'untranslated regions (UTR) to mediate mRNA decay. We hypothesized that ZFP36L1 is a negative regulator of a posttranscriptional hub involved in mRNA half-life regulation of cancer-related transcripts. Analysis of in silico data revealed that ZFP36L1 was significantly mutated, epigenetically silenced, and downregulated in a variety of cancers. Forced expression of ZFP36L1 in cancer cells markedly reduced cell proliferation in vitro and in vivo, whereas silencing of ZFP36L1 enhanced tumor cell growth. To identify direct downstream targets of ZFP36L1, systematic screening using RNA pull-down of wild-type and mutant ZFP36L1 as well as whole transcriptome sequencing of bladder cancer cells {plus minus} tet-on ZFP36L1 was performed. A network of 1,410 genes was identified as potential direct targets of ZFP36L1. These targets included a number of key oncogenic transcripts such as HIF1A, CCND1, and E2F1. ZFP36L1 specifically bound to the 3'UTRs of these targets for mRNA degradation, thus suppressing their expression. Dual luciferase reporter assays and RNA electrophoretic mobility shift assays showed that wild-type, but not zinc-finger mutant ZFP36L1, bound to HIF1A 3'UTR and mediated HIF1A mRNA degradation, leading to reduced expression of HIF1A and its downstream targets. Collectively, our findings reveal an indispensable role of ZFP36L1 as a posttranscriptional safeguard against aberrant hypoxic signaling and abnormal cell-cycle progression. SIGNIFICANCE: RNA-binding protein ZFP36L1 functions as a tumor suppressor by regulating the mRNA stability of a number of mRNAs involved in hypoxia and cell-cycle signaling.