Dynamic changes in RNA m6A and 5 hmC influence gene expression programs during macrophage differentiation and polarisation.

RNA modifications are essential for the establishment of cellular identity. Although increasing evidence indicates that RNA modifications regulate the innate immune response, their role in monocyte-to-macrophage differentiation and polarisation is unclear. While m6A has been widely studied, other RNA modifications, including 5 hmC, remain poorly characterised. We profiled m6A and 5 hmC epitranscriptomes, transcriptomes, translatomes and proteomes of monocytes and macrophages at rest and pro- and anti-inflammatory states. Transcriptome-wide mapping of m6A and 5 hmC reveals enrichment of m6A and/or 5 hmC on specific categories of transcripts essential for macrophage differentiation. Our analyses indicate that m6A and 5 hmC modifications are present in transcripts with critical functions in pro- and anti-inflammatory macrophages. Notably, we also discover the co-occurrence of m6A and 5 hmC on alternatively-spliced isoforms and/or opposing ends of the untranslated regions (UTR) of mRNAs with key roles in macrophage biology. In specific examples, RNA 5 hmC controls the decay of transcripts independently of m6A. This study provides (i) a comprehensive dataset to interrogate the role of RNA modifications in a plastic system (ii) a resource for exploring different layers of gene expression regulation in the context of human monocyte-to-macrophage differentiation and polarisation, (iii) new insights into RNA modifications as central regulators of effector cells in innate immunity.

A multiomics dataset for the study of RNA modifications in human macrophage differentiation and polarisation.

RNA modifications have emerged as central regulators of gene expression programs. Amongst RNA modifications are N6-methyladenosine (m6A) and RNA 5-hydroxymethylcytosine (5hmC). While m6A is established as a versatile regulator of RNA metabolism, the functions of RNA 5hmC are unclear. Despite some evidence linking RNA modifications to immunity, their implications in gene expression control in macrophage development and functions remain unclear. Here we present a multi-omics dataset capturing different layers of the gene expression programs driving macrophage differentiation and polarisation. We obtained mRNA-Seq, m6A-IP-Seq, 5hmC-IP-Seq, Polyribo-Seq and LC-MS/MS data from monocytes and resting-, pro- and anti-inflammatory-like macrophages. We present technical validation showing high quality and correlation between samples for all datasets, and evidence of biological consistency of modelled macrophages at the transcriptomic, epitranscriptomic, translational and proteomic levels. This multi-omics dataset provides a resource for the study of RNA m6A and 5hmC in the context of macrophage biology and spans the gene expression process from transcripts to proteins.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer.

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

Methylglyoxal: a novel upstream regulator of DNA methylation.

BACKGROUND: Aerobic glycolysis, also known as the Warburg effect, is predominantly upregulated in a variety of solid tumors, including breast cancer. We have previously reported that methylglyoxal (MG), a very reactive by-product of glycolysis, unexpectedly enhanced the metastatic potential in triple negative breast cancer (TNBC) cells. MG and MG-derived glycation products have been associated with various diseases, such as diabetes, neurodegenerative disorders, and cancer. Glyoxalase 1 (GLO1) exerts an anti-glycation defense by detoxifying MG to D-lactate. METHODS: Here, we used our validated model consisting of stable GLO1 depletion to induce MG stress in TNBC cells. Using genome-scale DNA methylation analysis, we report that this condition resulted in DNA hypermethylation in TNBC cells and xenografts. RESULTS: GLO1-depleted breast cancer cells showed elevated expression of DNMT3B methyltransferase and significant loss of metastasis-related tumor suppressor genes, as assessed using integrated analysis of methylome and transcriptome data. Interestingly, MG scavengers revealed to be as potent as typical DNA demethylating agents at triggering the re-expression of representative silenced genes. Importantly, we delineated an epigenomic MG signature that effectively stratified TNBC patients based on survival. CONCLUSION: This study emphasizes the importance of MG oncometabolite, occurring downstream of the Warburg effect, as a novel epigenetic regulator and proposes MG scavengers to reverse altered patterns of gene expression in TNBC.

Improving Infinium MethylationEPIC data processing: re-annotation of enhancers and long noncoding RNA genes and benchmarking of normalization methods.

Illumina Infinium DNA Methylation (5mC) arrays are a popular technology for low-cost, high-throughput, genome-scale measurement of 5mC distribution, especially in cancer and other complex diseases. After the success of its HumanMethylation450 array (450k), Illumina released the MethylationEPIC array (850k) featuring increased coverage of enhancers. Despite the widespread use of 850k, analysis of the corresponding data remains suboptimal: it still relies mostly on Illumina's default annotation, which underestimates enhancerss and long noncoding RNAs. Results: We have thus developed an approach, based on the ENCODE and LNCipedia databases, which greatly improves upon Illumina's default annotation of enhancers and long noncoding transcripts. We compared the re-annotated 850k with both 450k and reduced-representation bisulphite sequencing (RRBS), another high-throughput 5mC profiling technology. We found 850k to cover at least three times as many enhancers and long noncoding RNAs as either 450k or RRBS. We further investigated the reproducibility of the three technologies, applying various normalization methods to the 850k data. Most of these methods reduced variability to a level below that of RRBS data. We then used 850k with our new annotation and normalization to profile 5mC changes in breast cancer biopsies. 850k highlighted aberrant enhancer methylation as the predominant feature, in agreement with previous reports. Our study provides an updated processing approach for 850k data, based on refined probe annotation and normalization, allowing for improved analysis of methylation at enhancers and long noncoding RNA genes. Our findings will help to further advance understanding of the DNA methylome in health and disease.

Mesp1 controls the chromatin and enhancer landscapes essential for spatiotemporal patterning of early cardiovascular progenitors.

The mammalian heart arises from various populations of Mesp1-expressing cardiovascular progenitors (CPs) that are specified during the early stages of gastrulation. Mesp1 is a transcription factor that acts as a master regulator of CP specification and differentiation. However, how Mesp1 regulates the chromatin landscape of nascent mesodermal cells to define the temporal and spatial patterning of the distinct populations of CPs remains unknown. Here, by combining ChIP-seq, RNA-seq and ATAC-seq during mouse pluripotent stem cell differentiation, we defined the dynamic remodelling of the chromatin landscape mediated by Mesp1. We identified different enhancers that are temporally regulated to erase the pluripotent state and specify the pools of CPs that mediate heart development. We identified Zic2 and Zic3 as essential cofactors that act with Mesp1 to regulate its transcription-factor activity at key mesodermal enhancers, thereby regulating the chromatin remodelling and gene expression associated with the specification of the different populations of CPs in vivo. Our study identifies the dynamics of the chromatin landscape and enhancer remodelling associated with temporal patterning of early mesodermal cells into the distinct populations of CPs that mediate heart development.

The chromatin remodelling protein LSH/HELLS regulates the amount and distribution of DNA hydroxymethylation in the genome.

Ten-Eleven Translocation (TET) proteins convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) leading to a dynamic epigenetic state of DNA that can influence transcription and chromatin organization. While TET proteins interact with complexes involved in transcriptional repression and activation, the overall understanding of the molecular mechanisms involved in TET-mediated regulation of gene expression still remains limited. Here, we show that TET proteins interact with the chromatin remodelling protein lymphoid-specific helicase (LSH/HELLS) in vivo and in vitro. In mouse embryonic fibroblasts (MEFs) and embryonic stem cells (ESCs) knock out of Lsh leads to a significant reduction of 5-hydroxymethylation amount in the DNA. Whole genome sequencing of 5hmC in wild-type versus Lsh knock-out MEFs and ESCs showed that in absence of Lsh, some regions of the genome gain 5hmC while others lose it, with mild correlation with gene expression changes. We further show that differentially hydroxymethylated regions did not completely overlap with differentially methylated regions indicating that changes in 5hmC distribution upon Lsh knock-out are not a direct consequence of 5mC decrease. Altogether, our results suggest that LSH, which interacts with TET proteins, contributes to the regulation of 5hmC levels and distribution in MEFs and ESCs.

Downregulation of the FTO m6A RNA demethylase promotes EMT-mediated progression of epithelial tumors and sensitivity to Wnt inhibitors.

Post-transcriptional modifications of RNA constitute an emerging regulatory layer of gene expression. The demethylase fat mass- and obesity-associated protein (FTO), an eraser of N6-methyladenosine (m6A), has been shown to play a role in cancer, but its contribution to tumor progression and the underlying mechanisms remain unclear. Here, we report widespread FTO downregulation in epithelial cancers associated with increased invasion, metastasis and worse clinical outcome. Both in vitro and in vivo, FTO silencing promotes cancer growth, cell motility and invasion. In human-derived tumor xenografts (PDXs), FTO pharmacological inhibition favors tumorigenesis. Mechanistically, we demonstrate that FTO depletion elicits an epithelial-to-mesenchymal transition (EMT) program through increased m6A and altered 3'-end processing of key mRNAs along the Wnt signaling cascade. Accordingly, FTO knockdown acts via EMT to sensitize mouse xenografts to Wnt inhibition. We thus identify FTO as a key regulator, across epithelial cancers, of Wnt-triggered EMT and tumor progression and reveal a therapeutically exploitable vulnerability of FTO-low tumors.

Functional role of Tet-mediated RNA hydroxymethylcytosine in mouse ES cells and during differentiation.

Tet-enzyme-mediated 5-hydroxymethylation of cytosines in DNA plays a crucial role in mouse embryonic stem cells (ESCs). In RNA also, 5-hydroxymethylcytosine (5hmC) has recently been evidenced, but its physiological roles are still largely unknown. Here we show the contribution and function of this mark in mouse ESCs and differentiating embryoid bodies. Transcriptome-wide mapping in ESCs reveals hundreds of messenger RNAs marked by 5hmC at sites characterized by a defined unique consensus sequence and particular features. During differentiation a large number of transcripts, including many encoding key pluripotency-related factors (such as Eed and Jarid2), show decreased cytosine hydroxymethylation. Using Tet-knockout ESCs, we find Tet enzymes to be partly responsible for deposition of 5hmC in mRNA. A transcriptome-wide search further reveals mRNA targets to which Tet1 and Tet2 bind, at sites showing a topology similar to that of 5hmC sites. Tet-mediated RNA hydroxymethylation is found to reduce the stability of crucial pluripotency-promoting transcripts. We propose that RNA cytosine 5-hydroxymethylation by Tets is a mark of transcriptome flexibility, inextricably linked to the balance between pluripotency and lineage commitment.

Monocytes undergo multi-step differentiation in mice during oral infection by Toxoplasma gondii.

Monocytes play a major role in the defense against pathogens. They are rapidly mobilized to inflamed sites where they exert both proinflammatory and regulatory effector functions. It is still poorly understood how this dynamic and exceptionally plastic system is controlled at the molecular level. Herein, we evaluated the differentiation process that occurs in Ly6Chi monocytes during oral infection by Toxoplasma gondii. Flow cytometry and single-cell analysis revealed distinct activation status and gene expression profiles in the bone marrow, the spleen and the lamina propria of infected mice. We provide further evidence that acquisition of effector functions, such as the capacity to produce interleukin-27, is accompanied by distinct waves of epigenetic programming, highlighting a role for STAT1/IRF1 in the bone marrow and AP-1/NF-κB in the periphery. This work broadens our understanding of the molecular events that occur in vivo during monocyte differentiation in response to inflammatory cues.

EOMES interacts with RUNX3 and BRG1 to promote innate memory cell formation through epigenetic reprogramming.

Memory CD8+ T cells have the ability to provide lifelong immunity against pathogens. Although memory features generally arise after challenge with a foreign antigen, naïve CD8 single positive (SP) thymocytes may acquire phenotypic and functional characteristics of memory cells in response to cytokines such as interleukin-4. This process is associated with the induction of the T-box transcription factor Eomesodermin (EOMES). However, the underlying molecular mechanisms remain ill-defined. Using epigenomic profiling, we show that these innate memory CD8SP cells acquire only a portion of the active enhancer repertoire of conventional memory cells. This reprograming is secondary to EOMES recruitment, mostly to RUNX3-bound enhancers. Furthermore, EOMES is found within chromatin-associated complexes containing BRG1 and promotes the recruitment of this chromatin remodelling factor. Also, the in vivo acquisition of EOMES-dependent program is BRG1-dependent. In conclusion, our results support a strong epigenetic basis for the EOMES-driven establishment of CD8+ T cell innate memory program.

TET2-Dependent Hydroxymethylome Plasticity Reduces Melanoma Initiation and Progression.

Although numerous epigenetic aberrancies accumulate in melanoma, their contribution to initiation and progression remain unclear. The epigenetic mark 5-hydroxymethylcytosine (5hmC), generated through TET-mediated DNA modification, is now referred to as the sixth base of DNA and has recently been reported as a potential biomarker for multiple types of cancer. Loss of 5hmC is an epigenetic hallmark of melanoma, but whether a decrease in 5hmc levels contributes directly to pathogenesis or whether it merely results from disease progression-associated epigenetic remodeling remains to be established. Here, we show that NRAS-driven melanomagenesis in mice is accompanied by an overall decrease in 5hmC and specific 5hmC gains in selected gene bodies. Strikingly, genetic ablation of Tet2 in mice cooperated with oncogenic NRASQ61K to promote melanoma initiation while suppressing specific gains in 5hmC. We conclude that TET2 acts as a barrier to melanoma initiation and progression, partly by promoting 5hmC gains in specific gene bodies. SIGNIFICANCE: This work emphasizes the importance of epigenome plasticity in cancer development and highlights the involvement of druggable epigenetic factors in cancer.

Comprehensive identification of long noncoding RNAs in colorectal cancer.

Colorectal cancer (CRC) is one of the most common cancers in humans and a leading cause of cancer-related deaths worldwide. As in the case of other cancers, CRC heterogeneity leads to a wide range of clinical outcomes and complicates therapy. Over the years, multiple factors have emerged as markers of CRC heterogeneity, improving tumor classification and selection of therapeutic strategies. Understanding the molecular mechanisms underlying this heterogeneity remains a major challenge. A considerable research effort is therefore devoted to identifying additional features of colorectal tumors, in order to better understand CRC etiology and to multiply therapeutic avenues. Recently, long noncoding RNAs (lncRNAs) have emerged as important players in physiological and pathological processes, including CRC. Here we looked for lncRNAs that might contribute to the various colorectal tumor phenotypes. We thus monitored the expression of 4898 lncRNA genes across 566 CRC samples and identified 282 lncRNAs reflecting CRC heterogeneity. We then inferred potential functions of these lncRNAs. Our results highlight lncRNAs that may participate in the major processes altered in distinct CRC cases, such as WNT/ß-catenin and TGF-ß signaling, immunity, the epithelial-to-mesenchymal transition (EMT), and angiogenesis. For several candidates, we provide experimental evidence supporting our functional predictions that they may be involved in the cell cycle or the EMT. Overall, our work identifies lncRNAs associated with key CRC characteristics and provides insights into their respective functions. Our findings constitute a further step towards understanding the contribution of lncRNAs to CRC heterogeneity. They may open new therapeutic opportunities.

Immunity drives TET1 regulation in cancer through NF-κB.

Ten-eleven translocation enzymes (TET1, TET2, and TET3), which induce DNA demethylation and gene regulation by converting 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), are often down-regulated in cancer. We uncover, in basal-like breast cancer (BLBC), genome-wide 5hmC changes related to TET1 regulation. We further demonstrate that TET1 repression is associated with high expression of immune markers and high infiltration by immune cells. We identify in BLBC tissues an anticorrelation between TET1 expression and the major immunoregulator family nuclear factor κB (NF-κB). In vitro and in mice, TET1 is down-regulated in breast cancer cells upon NF-κB activation through binding of p65 to its consensus sequence in the TET1 promoter. We lastly show that these findings extend to other cancer types, including melanoma, lung, and thyroid cancers. Together, our data suggest a novel mode of regulation for TET1 in cancer and highlight a new paradigm in which the immune system can influence cancer cell epigenetics.

The transcription factors Runx3 and ThPOK cross-regulate acquisition of cytotoxic function by human Th1 lymphocytes.

Cytotoxic CD4 (CD4CTX) T cells are emerging as an important component of antiviral and antitumor immunity, but the molecular basis of their development remains poorly understood. In the context of human cytomegalovirus infection, a significant proportion of CD4 T cells displays cytotoxic functions. We observed that the transcriptional program of these cells was enriched in CD8 T cell lineage genes despite the absence of ThPOK downregulation. We further show that establishment of CD4CTX-specific transcriptional and epigenetic programs occurred in a stepwise fashion along the Th1-differentiation pathway. In vitro, prolonged activation of naive CD4 T cells in presence of Th1 polarizing cytokines led to the acquisition of perforin-dependent cytotoxic activity. This process was dependent on the Th1 transcription factor Runx3 and was limited by the sustained expression of ThPOK. This work elucidates the molecular program of human CD4CTX T cells and identifies potential targets for immunotherapy against viral infections and cancer.

DNA methylation-based immune response signature improves patient diagnosis in multiple cancers.

BACKGROUND: The tumor immune response is increasingly associated with better clinical outcomes in breast and other cancers. However, the evaluation of tumor-infiltrating lymphocytes (TILs) relies on histopathological measurements with limited accuracy and reproducibility. Here, we profiled DNA methylation markers to identify a methylation of TIL (MeTIL) signature that recapitulates TIL evaluations and their prognostic value for long-term outcomes in breast cancer (BC). METHODS: MeTIL signature scores were correlated with clinical endpoints reflecting overall or disease-free survival and a pathologic complete response to preoperative anthracycline therapy in 3 BC cohorts from the Jules Bordet Institute in Brussels and in other cancer types from The Cancer Genome Atlas. RESULTS: The MeTIL signature measured TIL distributions in a sensitive manner and predicted survival and response to chemotherapy in BC better than did histopathological assessment of TILs or gene expression-based immune markers, respectively. The MeTIL signature also improved the prediction of survival in other malignancies, including melanoma and lung cancer. Furthermore, the MeTIL signature predicted differences in survival for malignancies in which TILs were not known to have a prognostic value. Finally, we showed that MeTIL markers can be determined by bisulfite pyrosequencing of small amounts of DNA from formalin-fixed, paraffin-embedded tumor tissue, supporting clinical applications for this methodology. CONCLUSIONS: This study highlights the power of DNA methylation to evaluate tumor immune responses and the potential of this approach to improve the diagnosis and treatment of breast and other cancers. FUNDING: This work was funded by the Fonds National de la Recherche Scientifique (FNRS) and Télévie, the INNOVIRIS Brussels Region BRUBREAST Project, the IUAP P7/03 program, the Belgian "Foundation against Cancer," the Breast Cancer Research Foundation (BCRF), and the Fonds Gaston Ithier.

Portraying breast cancers with long noncoding RNAs.

Evidence is emerging that long noncoding RNAs (lncRNAs) may play a role in cancer development, but this role is not yet clear. We performed a genome-wide transcriptional survey to explore the lncRNA landscape across 995 breast tissue samples. We identified 215 lncRNAs whose genes are aberrantly expressed in breast tumors, as compared to normal samples. Unsupervised hierarchical clustering of breast tumors on the basis of their lncRNAs revealed four breast cancer subgroups that correlate tightly with PAM50-defined mRNA-based subtypes. Using multivariate analysis, we identified no less than 210 lncRNAs prognostic of clinical outcome. By analyzing the coexpression of lncRNA genes and protein-coding genes, we inferred potential functions of the 215 dysregulated lncRNAs. We then associated subtype-specific lncRNAs with key molecular processes involved in cancer. A correlation was observed, on the one hand, between luminal A-specific lncRNAs and the activation of phosphatidylinositol 3-kinase, fibroblast growth factor, and transforming growth factor-ß pathways and, on the other hand, between basal-like-specific lncRNAs and the activation of epidermal growth factor receptor (EGFR)-dependent pathways and of the epithelial-to-mesenchymal transition. Finally, we showed that a specific lncRNA, which we called CYTOR, plays a role in breast cancer. We confirmed its predicted functions, showing that it regulates genes involved in the EGFR/mammalian target of rapamycin pathway and is required for cell proliferation, cell migration, and cytoskeleton organization. Overall, our work provides the most comprehensive analyses for lncRNA in breast cancers. Our findings suggest a wide range of biological functions associated with lncRNAs in breast cancer and provide a foundation for functional investigations that could lead to new therapeutic approaches.

The interplay between the lysine demethylase KDM1A and DNA methyltransferases in cancer cells is cell cycle dependent.

DNA methylation and histone modifications are key epigenetic regulators of gene expression, and tight connections are known between the two. DNA methyltransferases are upregulated in several tumors and aberrant DNA methylation profiles are a cancer hallmark. On the other hand, histone demethylases are upregulated in cancer cells. Previous work on ES cells has shown that the lysine demethylase KDM1A binds to DNMT1, thereby affecting DNA methylation. In cancer cells, the occurrence of this interaction has not been explored. Here we demonstrate in several tumor cell lines an interaction between KDM1A and both DNMT1 and DNMT3B. Intriguingly and in contrast to what is observed in ES cells, KDM1A depletion in cancer cells was found not to trigger any reduction in the DNMT1 or DNMT3B protein level or any change in DNA methylation. In the S-phase, furthermore, KDM1A and DNMT1 were found, to co-localize within the heterochromatin. Using P-LISA, we revealed substantially increased binding of KDM1A to DNMT1 during the S-phase. Together, our findings propose a mechanistic link between KDM1A and DNA methyltransferases in cancer cells and suggest that the KDM1A/DNMT1 interaction may play a role during replication. Our work also strengthens the idea that DNMTs can exert functions unrelated to act on DNA methylation.

RNA biochemistry. Transcriptome-wide distribution and function of RNA hydroxymethylcytosine.

Hydroxymethylcytosine, well described in DNA, occurs also in RNA. Here, we show that hydroxymethylcytosine preferentially marks polyadenylated RNAs and is deposited by Tet in Drosophila. We map the transcriptome-wide hydroxymethylation landscape, revealing hydroxymethylcytosine in the transcripts of many genes, notably in coding sequences, and identify consensus sites for hydroxymethylation. We found that RNA hydroxymethylation can favor mRNA translation. Tet and hydroxymethylated RNA are found to be most abundant in the Drosophila brain, and Tet-deficient fruitflies suffer impaired brain development, accompanied by decreased RNA hydroxymethylation. This study highlights the distribution, localization, and function of cytosine hydroxymethylation and identifies central roles for this modification in Drosophila.

MicroRNAs regulate KDM5 histone demethylases in breast cancer cells.

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression. Alteration of miRNA levels is common in tumors and contributes to the pathogenesis of human malignancies. In the present study we examined the role played by miR-137 in breast tumorigenesis. We found miR-137 levels to be lower in breast cancer cells than in their non-tumorigenic counterparts and observed reduced proliferation and migration of breast cancer cells overexpressing miR-137. We further identified KDM5B, a histone demethylase known to be involved in breast cancer tumorigenesis, as a target of miR-137. As the involvement of histone demethylases in cancer is still poorly understood and as the role of miRNAs in controlling epigenetic mechanisms in cancer is emerging, we broadened our study to the whole KDM5 histone demethylase family to see if the genes coding for these epigenetic enzymes might be regulated by miRNAs in cancer cells. We discovered that KDM5C is overexpressed in breast cancer cells, providing evidence that miR-138 regulates its expression. We found miR-138 overexpression to affect breast cancer cell proliferation. Altogether, our findings suggest that miRNAs may regulate KDM5 histone demethylase levels in breast cancer and thereby control breast cancer cell proliferation and migration.