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.

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.

Dynamic intron retention modulates gene expression in the monocytic differentiation pathway.

Monocytes play a crucial role in maintaining homeostasis and mediating a successful innate immune response. They also act as central players in diverse pathological conditions, thus making them an attractive therapeutic target. Within the bone marrow, monocytes arise from a committed precursor termed Common Monocyte Progenitor (cMoP). However, molecular mechanisms that regulate the differentiation of cMoP to various monocytic subsets remain unclear. Herein, we purified murine myeloid precursors for deep poly-A-enriched RNA sequencing to understand the role of alternative splicing in the development and differentiation of monocytes under homeostasis. Our analyses revealed intron retention to be the major alternative splicing mechanism involved in the monocyte differentiation cascade, especially in the differentiation of Ly6Chi monocytes to Ly6Clo monocytes. Furthermore, we found that the intron retention of key genes involved in the differentiation of murine Ly6Chi to Ly6Clo monocytes was also conserved in humans. Our data highlight the unique role of intron retention in the regulation of the monocytic differentiation pathway.

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.

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.

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

[Figure: see text].

Accurate detection for a wide range of mutation and editing sites of microRNAs from small RNA high-throughput sequencing profiles.

Various types of mutation and editing (M/E) events in microRNAs (miRNAs) can change the stabilities of pre-miRNAs and/or complementarities between miRNAs and their targets. Small RNA (sRNA) high-throughput sequencing (HTS) profiles can contain many mutated and edited miRNAs. Systematic detection of miRNA mutation and editing sites from the huge volume of sRNA HTS profiles is computationally difficult, as high sensitivity and low false positive rate (FPR) are both required. We propose a novel method (named MiRME) for an accurate and fast detection of miRNA M/E sites using a progressive sequence alignment approach which refines sensitivity and improves FPR step-by-step. From 70 sRNA HTS profiles with over 1.3 billion reads, MiRME has detected thousands of statistically significant M/E sites, including 3'-editing sites, 57 A-to-I editing sites (of which 32 are novel), as well as some putative non-canonical editing sites. We demonstrated that a few non-canonical editing sites were not resulted from mutations in genome by integrating the analysis of genome HTS profiles of two human cell lines, suggesting the existence of new editing types to further diversify the functions of miRNAs. Compared with six existing studies or methods, MiRME has shown much superior performance for the identification and visualization of the M/E sites of miRNAs from the ever-increasing sRNA HTS profiles.

Inference of gene interaction networks using conserved subsequential patterns from multiple time course gene expression datasets.

MOTIVATION: Deciphering gene interaction networks (GINs) from time-course gene expression (TCGx) data is highly valuable to understand gene behaviors (e.g., activation, inhibition, time-lagged causality) at the system level. Existing methods usually use a global or local proximity measure to infer GINs from a single dataset. As the noise contained in a single data set is hardly self-resolved, the results are sometimes not reliable. Also, these proximity measurements cannot handle the co-existence of the various in vivo positive, negative and time-lagged gene interactions. METHODS AND RESULTS: We propose to infer reliable GINs from multiple TCGx datasets using a novel conserved subsequential pattern of gene expression. A subsequential pattern is a maximal subset of genes sharing positive, negative or time-lagged correlations of one expression template on their own subsets of time points. Based on these patterns, a GIN can be built from each of the datasets. It is assumed that reliable gene interactions would be detected repeatedly. We thus use conserved gene pairs from the individual GINs of the multiple TCGx datasets to construct a reliable GIN for a species. We apply our method on six TCGx datasets related to yeast cell cycle, and validate the reliable GINs using protein interaction networks, biopathways and transcription factor-gene regulations. We also compare the reliable GINs with those GINs reconstructed by a global proximity measure Pearson correlation coefficient method from single datasets. It has been demonstrated that our reliable GINs achieve much better prediction performance especially with much higher precision. The functional enrichment analysis also suggests that gene sets in a reliable GIN are more functionally significant. Our method is especially useful to decipher GINs from multiple TCGx datasets related to less studied organisms where little knowledge is available except gene expression data.

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.

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.

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.

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.

Epicardial CCM2 Promotes Cardiac Development and Repair Via its Regulation on Cytoskeletal Reorganization.

The epicardium provides epicardial-derived cells and molecular signals to support cardiac development and regeneration. Zebrafish and mouse studies have shown that ccm2, a cerebral cavernous malformation disease gene, is essential for cardiac development. Endocardial cell-specific deletion of Ccm2 in mice has previously established that Ccm2 is essential for maintenance of the cardiac jelly for cardiac development during early gestation. The current study aimed to explore the function of Ccm2 in epicardial cells for heart development and regeneration. Through genetic deletion of Ccm2 in epicardial cells, our in vivo and ex vivo experiments revealed that Ccm2 is required by epicardial cells to support heart development. Ccm2 regulates epicardial cell adhesion, cell polarity, cell spreading, and migration. Importantly, the loss of Ccm2 in epicardial cells delays cardiac function recovery and aggravates cardiac fibrosis following myocardial infarction. Molecularly, Ccm2 targets the production of cytoskeletal and matrix proteins to maintain epicardial cell function and behaviors. Epicardial Ccm2 plays a critical role in heart development and regeneration via its regulation of cytoskeleton reorganization.

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.

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.

The N6-methyladenosine RNA landscape in the aged mouse hippocampus.

The aged brain is associated with an inevitable decline in cognitive function and increased vulnerability to neurodegenerative disorders. Multiple molecular hallmarks have been associated with the aging nervous system through transcriptomics and proteomic studies. Recently, epitranscriptomic analysis has highlighted the role of RNA chemical modification in various biological processes. In particular, N6-methyladenosine (m6A), the most abundant internal modification in eukaryotic mRNAs, has been functionally linked to multiple aspects of RNA metabolism with the roles of m6A in processes such as learning and memory, leading to our current investigation of how the m6A-transcriptomic landscape is shaped during aging. Using the inbred C57BL/6 line, we compared the m6A-transcriptomic profiles from the hippocampi of young (3-month-old) and aged (20-month-old) mice. Methylated RNA immunoprecipitation (MeRIP)-sequencing analysis revealed hyper- and hypomethylation in 426 and 102 genes, respectively, in the aged hippocampus (fold change >1.5, false discovery rate <0.05). By correlating the methylation changes to their steady-state transcript levels in the RNA-Seq data, we found a significant concordance between m6A and transcript levels in both directions. Notably, the myelin regulator gene Gpr17 was downregulated in the aged hippocampus concomitant with reduced m6A levels in its 3'UTR. Using reporter constructs and mutagenesis analysis, we demonstrated that the putative m6A sites in the 3'UTR of Gpr17 are important for mRNA translation but not for regulating transcript stability. Overall, the positive correlation between m6A and the transcript expression levels indicates a co-transcriptional regulation of m6A with gene expression changes that occur in the aged mouse hippocampus.