Pancreatic cancer is a very aggressive disease with dismal prognosis; peculiar is the tumor microenvironment characterized by an extensive fibrotic stroma, which favors rapid tumor progression. We previously reported that pancreatic cancer patients have a selective Th2 skew in the anti-carcinoembryonic antigen (CEA) CD4(+) T cell immunity, which correlates with the presence of a predominant GATA-3(+) tumor lymphoid infiltrate. This has negative effects in both effective anti-tumor immunity and further favoring fibrinogenesis. Aim of this study was to evaluate whether the Th2 polarization of CEA-specific CD4(+) T cells from pancreatic cancer patients is stable or can be reverted by immunomodulating cytokines.

Protein-RNA interactions are implicated in a number of physiological roles as well as diseases, with molecular mechanisms ranging from defects in RNA splicing, localization and translation to the formation of aggregates. Currently, ∼1400 human proteins have experimental evidence of RNA-binding activity. However, only ∼250 of these proteins currently have experimental data on their target RNAs from various sequencing-based methods such as eCLIP. To bridge this gap, we used an established, computationally expensive protein-RNA interaction prediction method, catRAPID, to populate a large database, RNAct. RNAct allows easy lookup of known and predicted interactions and enables global views of the human, mouse and yeast protein-RNA interactomes, expanding them in a genome-wide manner far beyond experimental data (http://rnact.crg.eu).

Here we introduce omiXcore, a server for calculations of protein binding to large RNAs (> 500 nucleotides). Our webserver allows (i) use of both protein and RNA sequences without size restriction, (ii) pre-compiled library for exploration of human long intergenic RNAs interactions and (iii) prediction of binding sites.

Proteins and RNAs assemble in membrane-less organelles that organize intracellular spaces and regulate biochemical reactions. The ability of proteins and RNAs to form condensates is encoded in their sequences, yet it is unknown which domains drive the phase separation (PS) process and what are their specific roles. Here, we systematically investigated the human and yeast proteomes to find regions promoting condensation. Using advanced computational methods to predict the PS propensity of proteins, we designed a set of experiments to investigate the contributions of Prion-Like Domains (PrLDs) and RNA-binding domains (RBDs). We found that one PrLD is sufficient to drive PS, whereas multiple RBDs are needed to modulate the dynamics of the assemblies. In the case of stress granule protein Pub1 we show that the PrLD promotes sequestration of protein partners and the RBD confers liquid-like behaviour to the condensate. Our work sheds light on the fine interplay between RBDs and PrLD to regulate formation of membrane-less organelles, opening up the avenue for their manipulation.

Long noncoding RNAs (lncRNAs) are major regulators of physiological and disease-related gene expression, particularly in the central nervous system. Dysregulated lncRNA expression has been documented in several human cancers, and their tissue-specificity makes them attractive candidates as diagnostic/prognostic biomarkers and/or therapeutic agents. Here we show that linc-NeD125, which we previously characterized as a neuronal-induced lncRNA, is significantly overexpressed in Group 4 medulloblastomas (G4 MBs), the largest and least well characterized molecular MB subgroup. Mechanistically, linc-NeD125 is able to recruit the miRNA-induced silencing complex (miRISC) and to directly bind the microRNAs miR-19a-3p, miR-19b-3p and miR-106a-5p. Functionally, linc-NeD125 acts as a competing endogenous RNA (ceRNA) that, sequestering the three miRNAs, leads to de-repression of their targets CDK6, MYCN, SNCAIP, and KDM6A, which are major driver genes of G4 MB. Accordingly, linc-NeD125 downregulation reduces G4 cell proliferation. Moreover, we also provide evidence that linc-NeD125 ectopic expression in the aggressive Group 3 MB cells attenuates their proliferation, migration and invasion.This study unveils the first lncRNA-based ceRNA network in central nervous system tumours and provides a novel molecular circuit underlying the enigmatic Group 4 medulloblastoma.

To compare the secondary structure profiles of RNA molecules we developed the CROSSalign method. CROSSalign is based on the combination of the Computational Recognition Of Secondary Structure (CROSS) algorithm to predict the RNA secondary structure profile at single-nucleotide resolution and the Dynamic Time Warping (DTW) method to align profiles of different lengths. We applied CROSSalign to investigate the structural conservation of long non-coding RNAs such as XIST and HOTAIR as well as ssRNA viruses including HIV. CROSSalign performs pair-wise comparisons and is able to find homologs between thousands of matches identifying the exact regions of similarity between profiles of different lengths. In a pool of sequences with the same secondary structure CROSSalign accurately recognizes repeat A of XIST and domain D2 of HOTAIR and outperforms other methods based on covariance modeling. The algorithm is freely available at the webpage http://service.tartaglialab.com//new_submission/crossalign.

Exon skipping is an effective strategy for the treatment of many Duchenne Muscular Dystrophy (DMD) mutations. Natural exon skipping observed in several DMD cases can help in identifying novel therapeutic tools. Here, we show a DMD study case where the lack of a splicing factor (Celf2a), which results in exon skipping and dystrophin rescue, is due to a maternally inherited trans-generational epigenetic silencing. We found that the study case and his mother express a repressive long non-coding RNA, DUXAP8, whose presence correlates with silencing of the Celf2a coding region. We also demonstrate that DUXAP8 expression is lost upon cell reprogramming and that, upon induction of iPSCs into myoblasts, Celf2a expression is recovered leading to the loss of exon skipping and loss of dystrophin synthesis. Finally, CRISPR/Cas9 inactivation of the splicing factor Celf2a was proven to ameliorate the pathological state in other DMD backgrounds establishing Celf2a ablation or inactivation as a novel therapeutic approach for the treatment of Duchenne Muscular Dystrophy.

Biological condensates are membraneless organelles with different material properties. Proteins and RNAs are the main components, but most of their interactions are still unknown. Here, we introduce PRALINE, a database for the interrogation of proteins and RNAs contained in stress granules, processing bodies and other assemblies including droplets and amyloids. PRALINE provides information about the predicted and experimentally validated protein-protein, protein-RNA and RNA-RNA interactions. For proteins, it reports the liquid-liquid phase separation and liquid-solid phase separation propensities. For RNAs, it provides information on predicted secondary structure content. PRALINE shows detailed information on human single-nucleotide variants, their clinical significance and presence in protein and RNA binding sites, and how they can affect condensates' physical properties.

Pancreatic cancer is a very aggressive disease characterized by a marked desmoplasia with a predominant Th2 (GATA-3+) over Th1 (T-bet+) lymphoid infiltrate. We found that the ratio of GATA-3+/T-bet+ tumor-infiltrating lymphoid cells is an independent predictive marker of patient survival. Patients surgically treated for stage IB/III disease with a ratio inferior to the median value had a statistically significant prolonged overall survival, implying an active role for Th2 responses in disease progression. Thymic stromal lymphopoietin (TSLP), which favors Th2 cell polarization through myeloid dendritic cell (DC) conditioning, was secreted by cancer-associated fibroblasts (CAFs) after activation with tumor-derived tumor necrosis factor α and interleukin 1β. TSLP-containing supernatants from activated CAFs induced in vitro myeloid DCs to up-regulate the TSLP receptor (TSLPR), secrete Th2-attracting chemokines, and acquire TSLP-dependent Th2-polarizing capability in vitro. In vivo, Th2 chemoattractants were expressed in the tumor and in the stroma, and TSLPR-expressing DCs were present in the tumor stroma and in tumor-draining but not in nondraining lymph nodes. Collectively, this study identifies in pancreatic cancer a cross talk between tumor cells and CAFs, resulting in a TSLP-dependent induction of Th2-type inflammation which associates with reduced patient survival. Thus, blocking TSLP production by CAFs might help to improve prognosis in pancreatic cancer.

Specific elements of viral genomes regulate interactions within host cells. Here, we calculated the secondary structure content of >2000 coronaviruses and computed >100 000 human protein interactions with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genomic regions display different degrees of conservation. SARS-CoV-2 domain encompassing nucleotides 22 500-23 000 is conserved both at the sequence and structural level. The regions upstream and downstream, however, vary significantly. This part of the viral sequence codes for the Spike S protein that interacts with the human receptor angiotensin-converting enzyme 2 (ACE2). Thus, variability of Spike S is connected to different levels of viral entry in human cells within the population. Our predictions indicate that the 5' end of SARS-CoV-2 is highly structured and interacts with several human proteins. The binding proteins are involved in viral RNA processing, include double-stranded RNA specific editases and ATP-dependent RNA-helicases and have strong propensity to form stress granules and phase-separated assemblies. We propose that these proteins, also implicated in viral infections such as HIV, are selectively recruited by SARS-CoV-2 genome to alter transcriptional and post-transcriptional regulation of host cells and to promote viral replication.

Mutations in MECP2 cause several neurological disorders of which Rett syndrome (RTT) represents the best-defined condition. Although mainly working as a transcriptional repressor, MeCP2 is a multifunctional protein revealing several activities, the involvement of which in RTT remains obscure. Besides being mainly localized in the nucleus, MeCP2 associates with the centrosome, an organelle from which primary cilia originate. Primary cilia function as "sensory antennae" protruding from most cells, and a link between primary cilia and mental illness has recently been reported. We herein demonstrate that MeCP2 deficiency affects ciliogenesis in cultured cells, including neurons and RTT fibroblasts, and in the mouse brain. Consequently, the cilium-related Sonic Hedgehog pathway, which is essential for brain development and functioning, is impaired. Microtubule instability participates in these phenotypes that can be rescued by HDAC6 inhibition together with the recovery of RTT-related neuronal defects. Our data indicate defects of primary cilium as a novel pathogenic mechanism that by contributing to the clinical features of RTT might impact on proper cerebellum/brain development and functioning, thus providing a novel therapeutic target.

Changes in the epigenetic regulation of gene expression have a central role in evolution. Here, we extensively profiled a panel of human, chimpanzee, gorilla, orangutan, and macaque lymphoblastoid cell lines (LCLs), using ChIP-seq for five histone marks, ATAC-seq and RNA-seq, further complemented with whole genome sequencing (WGS) and whole genome bisulfite sequencing (WGBS). We annotated regulatory elements (RE) and integrated chromatin contact maps to define gene regulatory architectures, creating the largest catalog of RE in primates to date. We report that epigenetic conservation and its correlation with sequence conservation in primates depends on the activity state of the regulatory element. Our gene regulatory architectures reveal the coordination of different types of components and highlight the role of promoters and intragenic enhancers (gE) in the regulation of gene expression. We observe that most regulatory changes occur in weakly active gE. Remarkably, novel human-specific gE with weak activities are enriched in human-specific nucleotide changes. These elements appear in genes with signals of positive selection and human acceleration, tissue-specific expression, and particular functional enrichments, suggesting that the regulatory evolution of these genes may have contributed to human adaptation.

Long Interspersed Nuclear Elements-1s (L1s) are transposable elements that constitute most of the genome's transcriptional output yet have still largely unknown functions. Here we show that L1s are required for proper mouse brain corticogenesis operating as regulatory long non-coding RNAs. They contribute to the regulation of the balance between neuronal progenitors and differentiation, the migration of post-mitotic neurons and the proportions of different cell types. In cortical cultured neurons, L1 RNAs are mainly associated to chromatin and interact with the Polycomb Repressive Complex 2 (PRC2) protein subunits enhancer of Zeste homolog 2 (Ezh2) and suppressor of zeste 12 (Suz12). L1 RNA silencing influences PRC2's ability to bind a portion of its targets and the deposition of tri-methylated histone H3 (H3K27me3) marks. Our results position L1 RNAs as crucial signalling hubs for genome-wide chromatin remodelling, enabling the fine-tuning of gene expression during brain development and evolution.

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

microRNA abundance has been shown to depend on the amount of the microprocessor components or, in some cases, on specific auxiliary co-factors. In this paper, we show that the FUS/TLS (fused in sarcoma/translocated in liposarcoma) protein, associated with familial forms of Amyotrophic Lateral Sclerosis (ALS), contributes to the biogenesis of a specific subset of microRNAs. Among them, species with roles in neuronal function, differentiation and synaptogenesis were identified. We also show that FUS/TLS is recruited to chromatin at sites of their transcription and binds the corresponding pri-microRNAs. Moreover, FUS/TLS depletion leads to decreased Drosha level at the same chromatin loci. Limited FUS/TLS depletion leads to a reduced microRNA biogenesis and we suggest a possible link between FUS mutations affecting nuclear/cytoplasmic partitioning of the protein and altered neuronal microRNA biogenesis in ALS pathogenesis.

Endoribonucleases participate in almost every step of eukaryotic RNA metabolism, acting either as degradative or biosynthetic enzymes. We previously identified the founding member of the Eukaryotic EndoU ribonuclease family, whose components display unique biochemical features and are flexibly involved in important biological processes, such as ribosome biogenesis, tumorigenesis and viral replication. Here we report the discovery of the CG3303 gene product, which we named DendoU, as a novel family member in Drosophila. Functional characterisation revealed that DendoU is essential for Drosophila viability and nervous system activity. Pan-neuronal silencing of dendoU resulted in fly immature phenotypes, highly reduced lifespan and dramatic motor performance defects. Neuron-subtype selective silencing showed that DendoU is particularly important in cholinergic circuits. At the molecular level, we unveiled that DendoU is a positive regulator of the neurodegeneration-associated protein dTDP-43, whose downregulation recapitulates the ensemble of dendoU-dependent phenotypes. This interdisciplinary work, which comprehends in silico, in vitro and in vivo studies, unveils a relevant role for DendoU in Drosophila nervous system physio-pathology and highlights that DendoU-mediated neurotoxicity is, at least in part, contributed by dTDP-43 loss-of-function.

The combination of high-throughput sequencing and in vivo crosslinking approaches leads to the progressive uncovering of the complex interdependence between cellular transcriptome and proteome. Yet, the molecular determinants governing interactions in protein-RNA networks are not well understood. Here we investigated the relationship between the structure of an RNA and its ability to interact with proteins. Analysing in silico, in vitro and in vivo experiments, we find that the amount of double-stranded regions in an RNA correlates with the number of protein contacts. This relationship -which we call structure-driven protein interactivity- allows classification of RNA types, plays a role in gene regulation and could have implications for the formation of phase-separated ribonucleoprotein assemblies. We validate our hypothesis by showing that a highly structured RNA can rearrange the composition of a protein aggregate. We report that the tendency of proteins to phase-separate is reduced by interactions with specific RNAs.

Lamins (A/C and B) are major constituents of the nuclear lamina (NL). Structurally conserved lamina-associated domains (LADs) are formed by genomic regions that contact the NL. Lamins are also found in the nucleoplasm, with a yet unknown function. Here we map the genome-wide localization of lamin B1 in an euchromatin-enriched fraction of the mouse genome and follow its dynamics during the epithelial-to-mesenchymal transition (EMT). Lamin B1 associates with actively expressed and open euchromatin regions, forming dynamic euchromatin lamin B1-associated domains (eLADs) of about 0.3 Mb. Hi-C data link eLADs to the 3D organization of the mouse genome during EMT and correlate lamin B1 enrichment at topologically associating domain (TAD) borders with increased border strength. Having reduced levels of lamin B1 alters the EMT transcriptional signature and compromises the acquisition of mesenchymal traits. Thus, during EMT, the process of genome reorganization in mouse involves dynamic changes in eLADs.

Cells must coordinate the activation of thousands of replication origins dispersed throughout their genome. Active transcription is known to favor the formation of mammalian origins, although the role that RNA plays in this process remains unclear. We show that the ORC1 subunit of the human Origin Recognition Complex interacts with RNAs transcribed from genes with origins in their transcription start sites (TSSs), displaying a positive correlation between RNA binding and origin activity. RNA depletion, or the use of ORC1 RNA-binding mutant, result in inefficient activation of proximal origins, linked to impaired ORC1 chromatin release. ORC1 RNA binding activity resides in its intrinsically disordered region, involved in intra- and inter-molecular interactions, regulation by phosphorylation, and phase-separation. We show that RNA binding favors ORC1 chromatin release, by regulating its phosphorylation and subsequent degradation. Our results unveil a non-coding function of RNA as a dynamic component of the chromatin, orchestrating the activation of replication origins.

Enzymes of intermediary metabolism are often reported to have moonlighting functions as RNA-binding proteins and have regulatory roles beyond their primary activities. Human serine hydroxymethyltransferase (SHMT) is essential for the one-carbon metabolism, which sustains growth and proliferation in normal and tumour cells. Here, we characterize the RNA-binding function of cytosolic SHMT (SHMT1) in vitro and using cancer cell models. We show that SHMT1 controls the expression of its mitochondrial counterpart (SHMT2) by binding to the 5'untranslated region of the SHMT2 transcript (UTR2). Importantly, binding to RNA is modulated by metabolites in vitro and the formation of the SHMT1-UTR2 complex inhibits the serine cleavage activity of the SHMT1, without affecting the reverse reaction. Transfection of UTR2 in cancer cells controls SHMT1 activity and reduces cell viability. We propose a novel mechanism of SHMT regulation, which interconnects RNA and metabolites levels to control the cross-talk between cytosolic and mitochondrial compartments of serine metabolism.