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.

Zebrafish tsc1 and cxcl12a increase susceptibility to mycobacterial infection.

Regulation of host miRNA expression is a contested node that controls the host immune response to mycobacterial infection. The host must counter subversive efforts of pathogenic mycobacteria to launch a protective immune response. Here, we examine the role of miR-126 in the zebrafish-Mycobacterium marinum infection model and identify a protective role for infection-induced miR-126 through multiple effector pathways. We identified a putative link between miR-126 and the tsc1a and cxcl12a/ccl2/ccr2 signalling axes resulting in the suppression of non-tnfa expressing macrophage accumulation at early M. marinum granulomas. Mechanistically, we found a detrimental effect of tsc1a expression that renders zebrafish embryos susceptible to higher bacterial burden and increased cell death via mTOR inhibition. We found that macrophage recruitment driven by the cxcl12a/ccl2/ccr2 signalling axis was at the expense of the recruitment of classically activated tnfa-expressing macrophages and increased cell death around granulomas. Together, our results delineate putative pathways by which infection-induced miR-126 may shape an effective immune response to M. marinum infection in zebrafish embryos.

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.

Ku70 senses cytosolic DNA and assembles a tumor-suppressive signalosome.

The innate immune response contributes to the development or attenuation of acute and chronic diseases, including cancer. Microbial DNA and mislocalized DNA from damaged host cells can activate different host responses that shape disease outcomes. Here, we show that mice and humans lacking a single allele of the DNA repair protein Ku70 had increased susceptibility to the development of intestinal cancer. Mechanistically, Ku70 translocates from the nucleus into the cytoplasm where it binds to cytosolic DNA and interacts with the GTPase Ras and the kinase Raf, forming a tripartite protein complex and docking at Rab5+Rab7+ early-late endosomes. This Ku70-Ras-Raf signalosome activates the MEK-ERK pathways, leading to impaired activation of cell cycle proteins Cdc25A and CDK1, reducing cell proliferation and tumorigenesis. We also identified the domains of Ku70, Ras, and Raf involved in activating the Ku70 signaling pathway. Therapeutics targeting components of the Ku70 signalosome could improve the treatment outcomes in cancer.

CYP eicosanoid pathway mediates colon cancer-promoting effects of dietary linoleic acid.

Human and animal studies support that consuming a high level of linoleic acid (LA, 18:2ω-6), an essential fatty acid and key component of the human diet, increases the risk of colon cancer. However, results from human studies have been inconsistent, making it challenging to establish dietary recommendations for optimal LA intake. Given the importance of LA in the human diet, it is crucial to better understand the molecular mechanisms underlying its potential colon cancer-promoting effects. Using LC-MS/MS-based targeted lipidomics, we find that the cytochrome P450 (CYP) monooxygenase pathway is a major pathway for LA metabolism in vivo. Furthermore, CYP monooxygenase is required for the colon cancer-promoting effects of LA, since the LA-rich diet fails to exacerbate colon cancer in CYP monooxygenase-deficient mice. Finally, CYP monooxygenase mediates the pro-cancer effects of LA by converting LA to epoxy octadecenoic acids (EpOMEs), which have potent effects on promoting colon tumorigenesis via gut microbiota-dependent mechanisms. Overall, these results support that CYP monooxygenase-mediated conversion of LA to EpOMEs plays a crucial role in the health effects of LA, establishing a unique mechanistic link between dietary fatty acid intake and cancer risk. These results could help in developing more effective dietary guidelines for optimal LA intake and identifying subpopulations that may be especially vulnerable to LA's negative effects.

Overexpression of VIRMA confers vulnerability to breast cancers via the m6A-dependent regulation of unfolded protein response.

Virilizer-like m6A methyltransferase-associated protein (VIRMA) maintains the stability of the m6A writer complex. Although VIRMA is critical for RNA m6A deposition, the impact of aberrant VIRMA expression in human diseases remains unclear. We show that VIRMA is amplified and overexpressed in 15-20% of breast cancers. Of the two known VIRMA isoforms, the nuclear-enriched full-length but not the cytoplasmic-localised N-terminal VIRMA promotes m6A-dependent breast tumourigenesis in vitro and in vivo. Mechanistically, we reveal that VIRMA overexpression upregulates the m6A-modified long non-coding RNA, NEAT1, which contributes to breast cancer cell growth. We also show that VIRMA overexpression enriches m6A on transcripts that regulate the unfolded protein response (UPR) pathway but does not promote their translation to activate the UPR under optimal growth conditions. Under stressful conditions that are often present in tumour microenvironments, VIRMA-overexpressing cells display enhanced UPR and increased susceptibility to death. Our study identifies oncogenic VIRMA overexpression as a vulnerability that may be exploited for cancer therapy.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Tumor suppressor CEBPA interacts with and inhibits DNMT3A activity.

DNA methyltransferases (DNMTs) catalyze DNA methylation, and their functions in mammalian embryonic development and diseases including cancer have been extensively studied. However, regulation of DNMTs remains under study. Here, we show that CCAAT/enhancer binding protein α (CEBPA) interacts with the long splice isoform DNMT3A, but not the short isoform DNMT3A2. CEBPA, by interacting with DNMT3A N-terminus, blocks DNMT3A from accessing DNA substrate and thereby inhibits its activity. Recurrent tumor-associated CEBPA mutations, such as preleukemic CEBPAN321D mutation, which is particularly potent in causing AML with high mortality, disrupt DNMT3A association and cause aberrant DNA methylation, notably hypermethylation of PRC2 target genes. Consequently, leukemia cells with the CEBPAN321D mutation are hypersensitive to hypomethylation agents. Our results provide insights into the functional difference between DNMT3A isoforms and the regulation of de novo DNA methylation at specific loci in the genome. Our study also suggests a therapeutic strategy for the treatment of CEBPA-mutated leukemia with DNA-hypomethylating agents.

CCM2L (Cerebral Cavernous Malformation 2 Like) Deletion Aggravates Cerebral Cavernous Malformation Through Map3k3-KLF Signaling Pathway.

[Figure: see text].

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.

Macrophage development and activation involve coordinated intron retention in key inflammatory regulators.

Monocytes and macrophages are essential components of the innate immune system. Herein, we report that intron retention (IR) plays an important role in the development and function of these cells. Using Illumina mRNA sequencing, Nanopore direct cDNA sequencing and proteomics analysis, we identify IR events that affect the expression of key genes/proteins involved in macrophage development and function. We demonstrate that decreased IR in nuclear-detained mRNA is coupled with increased expression of genes encoding regulators of macrophage transcription, phagocytosis and inflammatory signalling, including ID2, IRF7, ENG and LAT. We further show that this dynamic IR program persists during the polarisation of resting macrophages into activated macrophages. In the presence of proinflammatory stimuli, intron-retaining CXCL2 and NFKBIZ transcripts are rapidly spliced, enabling timely expression of these key inflammatory regulators by macrophages. Our study provides novel insights into the molecular factors controlling vital regulators of the innate immune response.

Murine and related chapparvoviruses are nephro-tropic and produce novel accessory proteins in infected kidneys.

Mouse kidney parvovirus (MKPV) is a member of the provisional genus Chapparvovirus that causes renal disease in immune-compromised mice, with a disease course reminiscent of polyomavirus-associated nephropathy in immune-suppressed kidney transplant patients. Here we map four major MKPV transcripts, created by alternative splicing, to a common initiator region, and use mass spectrometry to identify "p10" and "p15" as novel chapparvovirus accessory proteins produced in MKPV-infected kidneys. p15 and the splicing-dependent putative accessory protein NS2 are conserved in all near-complete amniote chapparvovirus genomes currently available (from mammals, birds and a reptile). In contrast, p10 may be encoded only by viruses with >60% amino acid identity to MKPV. We show that MKPV is kidney-tropic and that the bat chapparvovirus DrPV-1 and a non-human primate chapparvovirus, CKPV, are also found in the kidneys of their hosts. We propose, therefore, that many mammal chapparvoviruses are likely to be nephrotropic.

The long noncoding RNA MALAT1 promotes tumor-driven angiogenesis by up-regulating pro-angiogenic gene expression.

Neuroblastoma is the most common solid tumor during early childhood. One of the key features of neuroblastoma is extensive tumor-driven angiogenesis due to hypoxia. However, the mechanism through which neuroblastoma cells drive angiogenesis is poorly understood. Here we show that the long noncoding RNA MALAT1 was upregulated in human neuroblastoma cell lines under hypoxic conditions. Conditioned media from neuroblastoma cells transfected with small interfering RNAs (siRNA) targeting MALAT1, compared with conditioned media from neuroblastoma cells transfected with control siRNAs, induced significantly less endothelial cell migration, invasion and vasculature formation. Microarray-based differential gene expression analysis showed that one of the genes most significantly down-regulated following MALAT1 suppression in human neuroblastoma cells under hypoxic conditions was fibroblast growth factor 2 (FGF2). RT-PCR and immunoblot analyses confirmed that MALAT1 suppression reduced FGF2 expression, and Enzyme-Linked Immunosorbent Assays revealed that transfection with MALAT1 siRNAs reduced FGF2 protein secretion from neuroblastoma cells. Importantly, addition of recombinant FGF2 protein to the cell culture media reversed the effects of MALAT1 siRNA on vasculature formation. Taken together, our data suggest that up-regulation of MALAT1 expression in human neuroblastoma cells under hypoxic conditions increases FGF2 expression and promotes vasculature formation, and therefore plays an important role in tumor-driven angiogenesis.

Identification of lung cancer miRNA-miRNA co-regulation networks through a progressive data refining approach.

Co-regulations of miRNAs have been much less studied than the research on regulations between miRNAs and their target genes, although these two problems are equally important for understanding the entire mechanisms of complex post-transcriptional regulations. The difficulty to construct a miRNA-miRNA co-regulation network lies in how to determine reliable miRNA pairs from various resources of data related to the same disease such as expression levels, gene ontology (GO) databases, and protein-protein interactions. Here we take a novel integrative approach to the discovery of miRNA-miRNA co-regulation networks. This approach can progressively refine the various types of data and the computational analysis results. Applied to three lung cancer miRNA expression data sets of different subtypes, our method has identified a miRNA-miRNA co-regulation network and co-regulating functional modules common to lung cancer. An example of these functional modules consists of genes SMAD2, ACVR1B, ACVR2A and ACVR2B. This module is synergistically regulated by let-7a/b/c/f, is enriched in the same GO category, and has a close proximity in the protein interaction network. We also find that the co-regulation network is scale free and that lung cancer related miRNAs have more synergism in the network. According to our literature survey and database validation, many of these results are biologically meaningful for understanding the mechanism of the complex post-transcriptional regulations in lung cancer.

Connecting rules from paired miRNA and mRNA expression data sets of HCV patients to detect both inverse and positive regulatory relationships.

BACKGROUND: Intensive research based on the inverse expression relationship has been undertaken to discover the miRNA-mRNA regulatory modules involved in the infection of Hepatitis C virus (HCV), the leading cause of chronic liver diseases. However, biological studies in other fields have found that inverse expression relationship is not the only regulatory relationship between miRNAs and their targets, and some miRNAs can positively regulate a mRNA by binding at the 5' UTR of the mRNA. RESULTS: This work focuses on the detection of both inverse and positive regulatory relationships from a paired miRNA and mRNA expression data set of HCV patients through a 'change-to-change' method which can derive connected discriminatory rules. Our study uncovered many novel miRNA-mRNA regulatory modules. In particular, it was revealed that GFRA2 is positively regulated by miR-557, miR-765 and miR-17-3p that probably bind at different locations of the 5' UTR of this mRNA. The expression relationship between GFRA2 and any of these three miRNAs has not been studied before, although separate research for this gene and these miRNAs have all drawn conclusions linked to hepatocellular carcinoma. This suggests that the binding of mRNA GFRA2 with miR-557, miR-765, or miR-17-3p, or their combinations, is worthy of further investigation by experimentation. We also report another mRNA QKI which has a strong inverse expression relationship with miR-129 and miR-493-3p which may bind at the 3' UTR of QKI with a perfect sequence match. Furthermore, the interaction between hsa-miR-129-5p (previous ID: hsa-miR-129) and QKI is supported with CLIP-Seq data from starBase. Our method can be easily extended for the expression data analysis of other diseases. CONCLUSION: Our rule discovery method is useful for integrating binding information and expression profile for identifying HCV miRNA-mRNA regulatory modules and can be applied to the study of the expression profiles of other complex human diseases.

Rule discovery and distance separation to detect reliable miRNA biomarkers for the diagnosis of lung squamous cell carcinoma.

BACKGROUND: Altered expression profiles of microRNAs (miRNAs) are linked to many diseases including lung cancer. miRNA expression profiling is reproducible and miRNAs are very stable. These characteristics of miRNAs make them ideal biomarker candidates. METHOD: This work is aimed to detect 2-and 3-miRNA groups, together with specific expression ranges of these miRNAs, to form simple linear discriminant rules for biomarker identification and biological interpretation. Our method is based on a novel committee of decision trees to derive 2-and 3-miRNA 100%-frequency rules. This method is applied to a data set of lung miRNA expression profiles of 61 squamous cell carcinoma (SCC) samples and 10 normal tissue samples. A distance separation technique is used to select the most reliable rules which are then evaluated on a large independent data set. RESULTS: We obtained four 2-miRNA and three 3-miRNA top-ranked rules. One important rule is that: If the expression level of miR-98 is above 7.356 and the expression level of miR-205 is below 9.601 (log2 quantile normalized MirVan miRNA Bioarray signals), then the sample is normal rather than cancerous with specificity and sensitivity both 100%. The classification performance of our best miRNA rules remarkably outperformed that by randomly selected miRNA rules. Our data analysis also showed that miR-98 and miR-205 have two common predicted target genes FZD3 and RPS6KA3, which are actually genes associated with carcinoma according to the Online Mendelian Inheritance in Man (OMIM) database. We also found that most of the chromosomal loci of these miRNAs have a high frequency of genomic alteration in lung cancer. On the independent data set (with balanced controls), the three miRNAs miR-126, miR-205 and miR-182 from our best rule can separate the two classes of samples at the accuracy of 84.49%, sensitivity of 91.40% and specificity of 77.14%. CONCLUSION: Our results indicate that rule discovery followed by distance separation is a powerful computational method to identify reliable miRNA biomarkers. The visualization of the rules and the clear separation between the normal and cancer samples by our rules will help biology experts for their analysis and biological interpretation.