Sterile inflammation via TRPM8 RNA-dependent TLR3-NF-kB/IRF3 activation promotes antitumor immunity in prostate cancer.

Inflammation is a common condition of prostate tissue, whose impact on carcinogenesis is highly debated. Microbial colonization is a well-documented cause of a small percentage of prostatitis cases, but it remains unclear what underlies the majority of sterile inflammation reported. Here, androgen- independent fluctuations of PSA expression in prostate cells have lead us to identify a prominent function of the Transient Receptor Potential Cation Channel Subfamily M Member 8 (TRPM8) gene in sterile inflammation. Prostate cells secret TRPM8 RNA into extracellular vesicles (EVs), which primes TLR3/NF-kB-mediated inflammatory signaling after EV endocytosis by epithelial cancer cells. Furthermore, prostate cancer xenografts expressing a translation-defective form of TRPM8 RNA contain less collagen type I in the extracellular matrix, significantly more infiltrating NK cells, and larger necrotic areas as compared to control xenografts. These findings imply sustained, androgen-independent expression of TRPM8 constitutes as a promoter of anticancer innate immunity, which may constitute a clinically relevant condition affecting prostate cancer prognosis.

Apurinic/Apyrimidinic Endodeoxyribonuclease 1 modulates RNA G-quadruplex folding of miR-92b and controls its expression in cancer cells.

In the last decade, several novel functions of the mammalian Apurinic/Apyrimidinic Endodeoxyribonuclease 1 (APE1) have been discovered, going far beyond its canonical function as DNA repair enzyme and unveiling its potential roles in cancer development. Indeed, it was shown to be involved in DNA G-quadruplex biology and RNA metabolism, most importantly in the miRNA maturation pathway and the decay of oxidized or abasic miRNAs during oxidative stress conditions. In recent years, several noncanonical pathways of miRNA biogenesis have emerged, with a specific focus on guanosine-rich precursors that can form RNA G-quadruplex (rG4) structures. Here, we show that several miRNA precursors, dysregulated upon APE1 depletion, contain an rG4 motif and that their corresponding target genes are up-regulated after APE1 depletion. We also demonstrate, both by in vitro assays and by using different cancer cell lines, that APE1 can modulate the folding of an rG4 structure contained in pre-miR-92b, with a mechanism strictly dependent on lysine residues present in its N-terminal disordered region. Furthermore, APE1 cellular depletion alters the maturation process of miR-92b, mainly affecting the shuttling between the nucleus and cytosol. Bioinformatic analysis of APE1-regulated rG4-containing miRNAs supports the relevance of our findings in cancer biology. Specifically, these miRNAs exhibit high prognostic significance in lung, cervical, and liver tumors, as suggested by their involvement in several cancer-related pathways.

CAG repeat expansion in the Huntington's disease gene shapes linear and circular RNAs biogenesis.

Alternative splicing (AS) appears to be altered in Huntington's disease (HD), but its significance for early, pre-symptomatic disease stages has not been inspected. Here, taking advantage of Htt CAG knock-in mouse in vitro and in vivo models, we demonstrate a correlation between Htt CAG repeat length and increased aberrant linear AS, specifically affecting neural progenitors and, in vivo, the striatum prior to overt behavioral phenotypes stages. Remarkably, a significant proportion (36%) of the aberrantly spliced isoforms are not-functional and meant to non-sense mediated decay (NMD). The expanded Htt CAG repeats further reflect on a previously neglected, global impairment of back-splicing, leading to decreased circular RNAs production in neural progenitors. Integrative transcriptomic analyses unveil a network of transcriptionally altered micro-RNAs and RNA-binding proteins (Celf, hnRNPs, Ptbp, Srsf, Upf1, Ythd2) which might influence the AS machinery, primarily in neural cells. We suggest that this unbalanced expression of linear and circular RNAs might alter neural fitness, contributing to HD pathogenesis.

G-Quadruplexes in RNA Biology: Recent Advances and Future Directions.

RNA G-quadruplexes (RG4s) are four-stranded structures known to control gene expression mechanisms, from transcription to protein synthesis, and DNA-related processes. Their potential impact on RNA biology allows these structures to shape cellular processes relevant to disease development, making their targeting for therapeutic purposes an attractive option. We review here the current knowledge on RG4s, focusing on the latest breakthroughs supporting the notion of transient structures that fluctuate dynamically in cellulo, their interplay with RNA modifications, their role in cell compartmentalization, and their deregulation impacting the host immune response. We emphasize RG4-binding proteins as determinants of their transient conformation and effectors of their biological functions.

Ythdf is a N6-methyladenosine reader that modulates Fmr1 target mRNA selection and restricts axonal growth in Drosophila.

N6-methyladenosine (m6 A) regulates a variety of physiological processes through modulation of RNA metabolism. This modification is particularly enriched in the nervous system of several species, and its dysregulation has been associated with neurodevelopmental defects and neural dysfunctions. In Drosophila, loss of m6 A alters fly behavior, albeit the underlying molecular mechanism and the role of m6 A during nervous system development have remained elusive. Here we find that impairment of the m6 A pathway leads to axonal overgrowth and misguidance at larval neuromuscular junctions as well as in the adult mushroom bodies. We identify Ythdf as the main m6 A reader in the nervous system, being required to limit axonal growth. Mechanistically, we show that the m6 A reader Ythdf directly interacts with Fmr1, the fly homolog of Fragile X mental retardation RNA binding protein (FMRP), to inhibit the translation of key transcripts involved in axonal growth regulation. Altogether, this study demonstrates that the m6 A pathway controls development of the nervous system and modulates Fmr1 target transcript selection.

QUADRatlas: the RNA G-quadruplex and RG4-binding proteins database.

RNA G-quadruplexes (RG4s) are non-canonical, disease-associated post-transcriptional regulators of gene expression whose functions are driven by RNA-binding proteins (RBPs). Being able to explore transcriptome-wide RG4 formation and interaction with RBPs is thus paramount to understanding how they are regulated and exploiting them as potential therapeutic targets. Towards this goal, we present QUADRatlas (https://rg4db.cibio.unitn.it), a database of experimentally-derived and computationally predicted RG4s in the human transcriptome, enriched with biological function and disease associations. As RBPs are key to their function, we mined known interactions of RG4s with such proteins, complemented with an extensive RBP binding sites dataset. Users can thus intersect RG4s with their potential regulators and effectors, enabling the formulation of novel hypotheses on RG4 regulation, function and pathogenicity. To support this capability, we provide analysis tools for predicting whether an RBP can bind RG4s, RG4 enrichment in a gene set, and de novo RG4 prediction. Genome-browser and table views allow exploring, filtering, and downloading the data quickly for individual genes and in batch. QUADRatlas is a significant step forward in our ability to understand the biology of RG4s, offering unmatched data content and enabling the integrated analysis of RG4s and their interactions with RBPs.

Translational enhancement by base editing of the Kozak sequence rescues haploinsufficiency.

A variety of single-gene human diseases are caused by haploinsufficiency, a genetic condition by which mutational inactivation of one allele leads to reduced protein levels and functional impairment. Translational enhancement of the spare allele could exert a therapeutic effect. Here we developed BOOST, a novel gene-editing approach to rescue haploinsufficiency loci by the change of specific single nucleotides in the Kozak sequence, which controls translation by regulating start codon recognition. We evaluated for translational strength 230 Kozak sequences of annotated human haploinsufficient genes and 4621 derived variants, which can be installed by base editing, by a high-throughput reporter assay. Of these variants, 149 increased the translation of 47 Kozak sequences, demonstrating that a substantial proportion of haploinsufficient genes are controlled by suboptimal Kozak sequences. Validation of 18 variants for 8 genes produced an average enhancement in an expression window compatible with the rescue of the genetic imbalance. Base editing of the NCF1 gene, whose monoallelic loss causes chronic granulomatous disease, resulted in the desired increase of NCF1 (p47phox) protein levels in a relevant cell model. We propose BOOST as a fine-tuned approach to modulate translation, applicable to the correction of dozens of haploinsufficient monogenic disorders independently of the causing mutation.

MODOMICS: a database of RNA modification pathways. 2021 update.

The MODOMICS database has been, since 2006, a manually curated and centralized resource, storing and distributing comprehensive information about modified ribonucleosides. Originally, it only contained data on the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. Over the years, prompted by the accumulation of new knowledge and new types of data, it has been updated with new information and functionalities. In this new release, we have created a catalog of RNA modifications linked to human diseases, e.g., due to mutations in genes encoding modification enzymes. MODOMICS has been linked extensively to RCSB Protein Data Bank, and sequences of experimentally determined RNA structures with modified residues have been added. This expansion was accompanied by including nucleotide 5'-monophosphate residues. We redesigned the web interface and upgraded the database backend. In addition, a search engine for chemically similar modified residues has been included that can be queried by SMILES codes or by drawing chemical molecules. Finally, previously available datasets of modified residues, biosynthetic pathways, and RNA-modifying enzymes have been updated. Overall, we provide users with a new, enhanced, and restyled tool for research on RNA modification. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

CHD8 suppression impacts on histone H3 lysine 36 trimethylation and alters RNA alternative splicing.

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD.

Genome-wide dynamics of RNA synthesis, processing, and degradation without RNA metabolic labeling.

The quantification of the kinetic rates of RNA synthesis, processing, and degradation are largely based on the integrative analysis of total and nascent transcription, the latter being quantified through RNA metabolic labeling. We developed INSPEcT-, a computational method based on the mathematical modeling of premature and mature RNA expression that is able to quantify kinetic rates from steady-state or time course total RNA-seq data without requiring any information on nascent transcripts. Our approach outperforms available solutions, closely recapitulates the kinetic rates obtained through RNA metabolic labeling, improves the ability to detect changes in transcript half-lives, reduces the cost and complexity of the experiments, and can be adopted to study experimental conditions in which nascent transcription cannot be readily profiled. Finally, we applied INSPEcT- to the characterization of post-transcriptional regulation landscapes in dozens of physiological and disease conditions. This approach was included in the INSPEcT Bioconductor package, which can now unveil RNA dynamics from steady-state or time course data, with or without the profiling of nascent RNA.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies.

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

Regulation of HuR structure and function by dihydrotanshinone-I.

The Human antigen R protein (HuR) is an RNA-binding protein that recognizes U/AU-rich elements in diverse RNAs through two RNA-recognition motifs, RRM1 and RRM2, and post-transcriptionally regulates the fate of target RNAs. The natural product dihydrotanshinone-I (DHTS) prevents the association of HuR and target RNAs in vitro and in cultured cells by interfering with the binding of HuR to RNA. Here, we report the structural determinants of the interaction between DHTS and HuR and the impact of DHTS on HuR binding to target mRNAs transcriptome-wide. NMR titration and Molecular Dynamics simulation identified the residues within RRM1 and RRM2 responsible for the interaction between DHTS and HuR. RNA Electromobility Shifts and Alpha Screen Assays showed that DHTS interacts with HuR through the same binding regions as target RNAs, stabilizing HuR in a locked conformation that hampers RNA binding competitively. HuR ribonucleoprotein immunoprecipitation followed by microarray (RIP-chip) analysis showed that DHTS treatment of HeLa cells paradoxically enriched HuR binding to mRNAs with longer 3'UTR and with higher density of U/AU-rich elements, suggesting that DHTS inhibits the association of HuR to weaker target mRNAs. In vivo, DHTS potently inhibited xenograft tumor growth in a HuR-dependent model without systemic toxicity.

DynaMIT: the dynamic motif integration toolkit.

De-novo motif search is a frequently applied bioinformatics procedure to identify and prioritize recurrent elements in sequences sets for biological investigation, such as the ones derived from high-throughput differential expression experiments. Several algorithms have been developed to perform motif search, employing widely different approaches and often giving divergent results. In order to maximize the power of these investigations and ultimately be able to draft solid biological hypotheses, there is the need for applying multiple tools on the same sequences and merge the obtained results. However, motif reporting formats and statistical evaluation methods currently make such an integration task difficult to perform and mostly restricted to specific scenarios. We thus introduce here the Dynamic Motif Integration Toolkit (DynaMIT), an extremely flexible platform allowing to identify motifs employing multiple algorithms, integrate them by means of a user-selected strategy and visualize results in several ways; furthermore, the platform is user-extendible in all its aspects. DynaMIT is freely available at http://cibioltg.bitbucket.org.

Exploring the Role of Fallopian Ciliated Cells in the Pathogenesis of High-Grade Serous Ovarian Cancer.

High-grade serous epithelial ovarian cancer (HGSOC) is the fifth leading cause of cancer death in women and the first among gynecological malignancies. Despite an initial response to standard chemotherapy, most HGSOC patients relapse. To improve treatment options, we must continue investigating tumor biology. Tumor characteristics (e.g., risk factors and epidemiology) are valuable clues to accomplish this task. The two most frequent risk factors for HGSOC are the lifetime number of ovulations, which is associated with increased oxidative stress in the pelvic area caused by ovulation fluid, and a positive family history due to genetic factors. In the attempt to identify novel genetic factors (i.e., genes) associated with HGSOC, we observed that several genes in linkage with HGSOC are expressed in the ciliated cells of the fallopian tube. This finding made us hypothesize that ciliated cells, despite not being the cell of origin for HGSOC, may take part in HGSOC tumor initiation. Specifically, malfunction of the ciliary beat impairs the laminar fluid flow above the fallopian tube epithelia, thus likely reducing the clearance of oxidative stress caused by follicular fluid. Herein, we review the up-to-date findings dealing with HGSOC predisposition with the hypothesis that fallopian ciliated cells take part in HGSOC onset. Finally, we review the up-to-date literature concerning genes that are located in genomic loci associated with epithelial ovarian cancer (EOC) predisposition that are expressed by the fallopian ciliated cells.

The BioMart community portal: an innovative alternative to large, centralized data repositories.

The BioMart Community Portal (www.biomart.org) is a community-driven effort to provide a unified interface to biomedical databases that are distributed worldwide. The portal provides access to numerous database projects supported by 30 scientific organizations. It includes over 800 different biological datasets spanning genomics, proteomics, model organisms, cancer data, ontology information and more. All resources available through the portal are independently administered and funded by their host organizations. The BioMart data federation technology provides a unified interface to all the available data. The latest version of the portal comes with many new databases that have been created by our ever-growing community. It also comes with better support and extensibility for data analysis and visualization tools. A new addition to our toolbox, the enrichment analysis tool is now accessible through graphical and web service interface. The BioMart community portal averages over one million requests per day. Building on this level of service and the wealth of information that has become available, the BioMart Community Portal has introduced a new, more scalable and cheaper alternative to the large data stores maintained by specialized organizations.

tRanslatome: an R/Bioconductor package to portray translational control.

UNLABELLED: High-throughput technologies have led to an explosion of genomic data available for automated analysis. The consequent possibility to simultaneously sample multiple layers of variation along the gene expression flow requires computational methods integrating raw information from different '-omics'. It has been recently demonstrated that translational control is a widespread phenomenon, with profound and still underestimated regulation capabilities. Although detecting changes in the levels of total messenger RNAs (mRNAs; the transcriptome), of polysomally loaded mRNAs (the translatome) and of proteins (the proteome) is experimentally feasible in a high-throughput way, the integration of these levels is still far from being robustly approached. Here we introduce tRanslatome, a new R/Bioconductor package, which is a complete platform for the simultaneous pairwise analysis of transcriptome, translatome and proteome data. The package includes most of the available statistical methods developed for the analysis of high-throughput data, allowing the parallel comparison of differentially expressed genes and the corresponding differentially enriched biological themes. Notably, it also enables the prediction of translational regulatory elements on mRNA sequences. The utility of this tool is demonstrated with two case studies. AVAILABILITY AND IMPLEMENTATION: tRanslatome is available in Bioconductor.

Hyper conserved elements in vertebrate mRNA 3'-UTRs reveal a translational network of RNA-binding proteins controlled by HuR.

Little is known regarding the post-transcriptional networks that control gene expression in eukaryotes. Additionally, we still need to understand how these networks evolve, and the relative role played in them by their sequence-dependent regulatory factors, non-coding RNAs (ncRNAs) and RNA-binding proteins (RBPs). Here, we used an approach that relied on both phylogenetic sequence sharing and conservation in the whole mapped 3'-untranslated regions (3'-UTRs) of vertebrate species to gain knowledge on core post-transcriptional networks. The identified human hyper conserved elements (HCEs) were predicted to be preferred binding sites for RBPs and not for ncRNAs, namely microRNAs and long ncRNAs. We found that the HCE map identified a well-known network that post-transcriptionally regulates histone mRNAs. We were then able to discover and experimentally confirm a translational network composed of RNA Recognition Motif (RRM)-type RBP mRNAs that are positively controlled by HuR, another RRM-type RBP. HuR shows a preference for these RBP mRNAs bound in stem-loop motifs, confirming its role as a 'regulator of regulators'. Analysis of the transcriptome-wide HCE distribution revealed a profile of prevalently small clusters separated by unconserved intercluster RNA stretches, which predicts the formation of discrete small ribonucleoprotein complexes in the 3'-UTRs.

The shaping of mRNA translation plasticity by RNA G-quadruplexes in cancer progression and therapy resistance.

Translational reprogramming in response to oncogenic signaling or microenvironmental stress factors shapes the proteome of cancer cells, enabling adaptation and phenotypic changes underlying cell plasticity, tumor progression and response to cancer therapy. Among the mechanisms regulating translation are RNA G-quadruplexes (RG4s), non-canonical four-stranded structures whose conformational modulation by small molecule ligands and RNA-binding proteins affects the expression of cancer proteins. Here, we discuss the role of RG4s in the regulation of mRNA translation by focusing on paradigmatic examples showing their contribution to adaptive mechanisms of mRNA translation in cancer.

The three YTHDF paralogs and VIRMA are strong cross-histotype tumor driver candidates among m6A core genes.

N6-Methyladenosine (m6A) is the most abundant internal modification in mRNAs. Despite accumulating evidence for the profound impact of m6A on cancer biology, there are conflicting reports that alterations in genes encoding the m6A machinery proteins can either promote or suppress cancer, even in the same tumor type. Using data from The Cancer Genome Atlas, we performed a pan-cancer investigation of 15 m6A core factors in nearly 10000 samples from 31 tumor types to reveal underlying cross-tumor patterns. Altered expression, largely driven by copy number variations at the chromosome arm level, results in the most common mode of dysregulation of these factors. YTHDF1, YTHDF2, YTHDF3 and VIRMA are the most frequently altered factors and the only ones to be uniquely altered when tumors are grouped according to the expression pattern of the m6A factors. These genes are also the only ones with coherent, pan-cancer predictive power for progression-free survival. On the contrary, METTL3, the most intensively studied m6A factor as a cancer target, shows much lower levels of alteration and no predictive power for patient survival. Therefore, we propose the non-enzymatic YTHDF and VIRMA genes as preferred subjects to dissect the role of m6A in cancer and as priority cancer targets.

CircHTT(2,3,4,5,6) - co-evolving with the HTT CAG-repeat tract - modulates Huntington's disease phenotypes.

Circular RNA (circRNA) molecules have critical functions during brain development and in brain-related disorders. Here, we identified and validated a circRNA, circHTT(2,3,4,5,6), stemming from the Huntington's disease (HD) gene locus that is most abundant in the central nervous system (CNS). We uncovered its evolutionary conservation in diverse mammalian species, and a correlation between circHTT(2,3,4,5,6) levels and the length of the CAG-repeat tract in exon-1 of HTT in human and mouse HD model systems. The mouse orthologue, circHtt(2,3,4,5,6), is expressed during embryogenesis, increases during nervous system development, and is aberrantly upregulated in the presence of the expanded CAG tract. While an IRES-like motif was predicted in circH TT (2,3,4,5,6), the circRNA does not appear to be translated in adult mouse brain tissue. Nonetheless, a modest, but consistent fraction of circHtt(2,3,4,5,6) associates with the 40S ribosomal subunit, suggesting a possible role in the regulation of protein translation. Finally, circHtt(2,3,4,5,6) overexpression experiments in HD-relevant STHdh striatal cells revealed its ability to modulate CAG expansion-driven cellular defects in cell-to-substrate adhesion, thus uncovering an unconventional modifier of HD pathology.