RNA-Based Strategies for Cancer Therapy: In Silico Design and Evaluation of ASOs for Targeted Exon Skipping.

Precision medicine in oncology has made significant progress in recent years by approving drugs that target specific genetic mutations. However, many cancer driver genes remain challenging to pharmacologically target ("undruggable"). To tackle this issue, RNA-based methods like antisense oligonucleotides (ASOs) that induce targeted exon skipping (ES) could provide a promising alternative. In this work, a comprehensive computational procedure is presented, focused on the development of ES-based cancer treatments. The procedure aims to produce specific protein variants, including inactive oncogenes and partially restored tumor suppressors. This novel computational procedure encompasses target-exon selection, in silico prediction of ES products, and identification of the best candidate ASOs for further experimental validation. The method was effectively employed on extensively mutated cancer genes, prioritized according to their suitability for ES-based interventions. Notable genes, such as NRAS and VHL, exhibited potential for this therapeutic approach, as specific target exons were identified and optimal ASO sequences were devised to induce their skipping. To the best of our knowledge, this is the first computational procedure that encompasses all necessary steps for designing ASO sequences tailored for targeted ES, contributing with a versatile and innovative approach to addressing the challenges posed by undruggable cancer driver genes and beyond.

Interleukin-8 in Colorectal Cancer: A Systematic Review and Meta-Analysis of Its Potential Role as a Prognostic Biomarker.

Among soluble actors that have emerged as druggable factors, the chemokine interleukin-8 (IL-8) has emerged as a possible determinant of response to immunotherapy and targeted treatment in several cancer types; however, its prognostic/predictive role in colorectal cancer (CRC) remains to be established. We: (i) conducted a systematic review of published literature on IL-8 expression in CRC; (ii) searched public transcriptomics databases; (iii) investigated IL-8 expression, by tumor and infiltrating cells, in a series of CRC samples; and (iv) carried out a meta-analysis of published literature correlating IL-8 expression and CRC prognosis. IL-8 possesses an important role as a mediator of the bidirectional crosstalk between tumor/stromal cells. Transcriptomic analysis indicated that specific IL-8 transcripts were significantly overexpressed in CRC compared to normal colon mucosa. Moreover, in our series we observed a statistically significant correlation between PTEN-loss and IL-8 expression by infiltrating mononuclear and tumor cells. In total, 12 papers met our meta-analysis inclusion criteria, demonstrating that high IL-8 levels significantly correlated with shorter overall survival and progression-free survival. Sensitivity analysis demonstrated a highly significant correlation with outcome for circulating, but not for tissue-detected, IL-8. IL-8 is overexpressed in CRC tissues and differentially produced by tumor or stromal components depending on CRC genetic background. Moreover, circulating IL-8 represents a strong prognostic factor in CRC, suggesting its use in the refining of prognostic CRC assessment and potentially the tailoring of therapeutic strategies in individual CRC patients.

Mouse dendritic cells in the steady state: Hypoxia, autophagy, and stem cell factor.

Dendritic cells (DCs) are innate immune cells with a central role in immunity and tolerance. Under steady-state, DCs are scattered in tissues as resting cells. Upon infection or injury, DCs get activated and acquire the full capacity to prime antigen-specific CD4+ and CD8+ T cells, thus bridging innate and adaptive immunity. By secreting different sets of cytokines and chemokines, DCs orchestrate diverse types of immune responses, from a classical proinflammatory to an alternative pro-repair one. DCs are highly heterogeneous, and physiological differences in tissue microenvironments greatly contribute to variations in DC phenotype. Oxygen tension is normally low in some lymphoid areas, including bone marrow (BM) hematopoietic niches; nevertheless, the possible impact of tissue hypoxia on DC physiology has been poorly investigated. We assessed whether DCs are hypoxic in BM and spleen, by staining for hypoxia-inducible-factor-1α subunit (HIF-1α), the master regulator of hypoxia-induced response, and pimonidazole (PIM), a hypoxic marker, and by flow cytometric analysis. Indeed, we observed that mouse DCs have a hypoxic phenotype in spleen and BM, and showed some remarkable differences between DC subsets. Notably, DCs expressing membrane c-kit, the receptor for stem cell factor (SCF), had a higher PIM median fluorescence intensity (MFI) than c-kit- DCs, both in the spleen and in the BM. To determine whether SCF (a.k.a. kit ligand) has a role in DC hypoxia, we evaluated molecular pathways activated by SCF in c-kit+ BM-derived DCs cultured in hypoxic conditions. Gene expression microarrays and gene set enrichment analysis supported the hypothesis that SCF had an impact on hypoxia response and inhibited autophagy-related gene sets. Our results suggest that hypoxic response and autophagy, and their modulation by SCF, can play a role in DC homeostasis at the steady state, in agreement with our previous findings on SCF's role in DC survival.

Mining sponge phenomena in RNA expression data.

In the last few years, the interactions among competing endogenous RNAs (ceRNAs) have been recognized as a key post-transcriptional regulatory mechanism in cell differentiation, tissue development, and disease. Notably, such sponge phenomena substracting active microRNAs from their silencing targets have been recognized as having a potential oncosuppressive, or oncogenic, role in several cancer types. Hence, the ability to predict sponges from the analysis of large expression data sets (e.g. from international cancer projects) has become an important data mining task in bioinformatics. We present a technique designed to mine sponge phenomena whose presence or absence may discriminate between healthy and unhealthy populations of samples in tumoral or normal expression data sets, thus providing lists of candidates potentially relevant in the pathology. With this aim, we search for pairs of elements acting as ceRNA for a given miRNA, namely, we aim at discovering miRNA-RNA pairs involved in phenomena which are clearly present in one population and almost absent in the other one. The results on tumoral expression data, concerning five different cancer types, confirmed the effectiveness of the approach in mining interesting knowledge. Indeed, 32 out of 33 miRNAs and 22 out of 25 protein-coding genes identified as top scoring in our analysis are corroborated by having been similarly associated with cancer processes in independent studies. In fact, the subset of miRNAs selected by the sponge analysis results in a significant enrichment of annotation for the KEGG32 pathway "microRNAs in cancer" when tested with the commonly used bioinformatic resource DAVID. Moreover, often the cancer datasets where our sponge analysis identified a miRNA as top scoring match the one reported already in the pertaining literature.

AI applications in functional genomics.

We review the current applications of artificial intelligence (AI) in functional genomics. The recent explosion of AI follows the remarkable achievements made possible by "deep learning", along with a burst of "big data" that can meet its hunger. Biology is about to overthrow astronomy as the paradigmatic representative of big data producer. This has been made possible by huge advancements in the field of high throughput technologies, applied to determine how the individual components of a biological system work together to accomplish different processes. The disciplines contributing to this bulk of data are collectively known as functional genomics. They consist in studies of: i) the information contained in the DNA (genomics); ii) the modifications that DNA can reversibly undergo (epigenomics); iii) the RNA transcripts originated by a genome (transcriptomics); iv) the ensemble of chemical modifications decorating different types of RNA transcripts (epitranscriptomics); v) the products of protein-coding transcripts (proteomics); and vi) the small molecules produced from cell metabolism (metabolomics) present in an organism or system at a given time, in physiological or pathological conditions. After reviewing main applications of AI in functional genomics, we discuss important accompanying issues, including ethical, legal and economic issues and the importance of explainability.

Transcriptional Profiling of Three Pseudomonas syringae pv. actinidiae Biovars Reveals Different Responses to Apoplast-Like Conditions Related to Strain Virulence on the Host.

Pseudomonas syringae pv. actinidiae is a phytopathogen that causes devastating bacterial canker in kiwifruit. Among five biovars defined by genetic, biochemical, and virulence traits, P. syringae pv. actinidiae biovar 3 (Psa3) is the most aggressive and is responsible for the most recent reported outbreaks; however, the molecular basis of its heightened virulence is unclear. Therefore, we designed the first P. syringae multistrain whole-genome microarray, encompassing biovars Psa1, Psa2, and Psa3 and the well-established model P. syringae pv. tomato, and analyzed early bacterial responses to an apoplast-like minimal medium. Transcriptomic profiling revealed i) the strong activation in Psa3 of all hypersensitive reaction and pathogenicity (hrp) and hrp conserved (hrc) cluster genes, encoding components of the type III secretion system required for bacterial pathogenicity and involved in responses to environmental signals; ii) potential repression of the hrp/hrc cluster in Psa2; and iii) activation of flagellum-dependent cell motility and chemotaxis genes in Psa1. The detailed investigation of three gene families encoding upstream regulatory proteins (histidine kinases, their cognate response regulators, and proteins with diguanylate cyclase or phosphodiesterase domains) indicated that cyclic di-GMP may be a key regulator of virulence in P. syringae pv. actinidiae biovars. The gene expression data were supported by the quantification of biofilm formation. Our findings suggest that diverse early responses to the host apoplast, even among bacteria belonging to the same pathovar, can lead to different virulence strategies and may explain the differing outcomes of infections. Based on our detailed structural analysis of hrp operons, we also propose a revision of hrp cluster organization and operon regulation in P. syringae.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A Census and Categorization Method of Epitranscriptomic Marks.

In the past few years, thorough investigation of chemical modifications operated in the cells on ribonucleic acid (RNA) molecules is gaining momentum. This new field of research has been dubbed "epitranscriptomics", in analogy to best-known epigenomics, to stress the potential of ensembles of RNA modifications to constitute a post-transcriptional regulatory layer of gene expression orchestrated by writer, reader, and eraser RNA-binding proteins (RBPs). In fact, epitranscriptomics aims at identifying and characterizing all functionally relevant changes involving both non-substitutional chemical modifications and editing events made to the transcriptome. Indeed, several types of RNA modifications that impact gene expression have been reported so far in different species of cellular RNAs, including ribosomal RNAs, transfer RNAs, small nuclear RNAs, messenger RNAs, and long non-coding RNAs. Supporting functional relevance of this largely unknown regulatory mechanism, several human diseases have been associated directly to RNA modifications or to RBPs that may play as effectors of epitranscriptomic marks. However, an exhaustive epitranscriptome's characterization, aimed to systematically classify all RNA modifications and clarify rules, actors, and outcomes of this promising regulatory code, is currently not available, mainly hampered by lack of suitable detecting technologies. This is an unfortunate limitation that, thanks to an unprecedented pace of technological advancements especially in the sequencing technology field, is likely to be overcome soon. Here, we review the current knowledge on epitranscriptomic marks and propose a categorization method based on the reference ribonucleotide and its rounds of modifications ("stages") until reaching the given modified form. We believe that this classification scheme can be useful to coherently organize the expanding number of discovered RNA modifications.

Plant and fungus transcriptomic data from grapevine berries undergoing artificially-induced noble rot caused by Botrytis cinerea.

Noble rot is a latent infection of grape berries caused by the necrotrophic fungus Botrytis cinerea, which develops under specific climatic conditions. The infected berries undergo biochemical and metabolic changes, associated with a rapid withering, which altogether offer interesting organoleptic features to sweet white wines. In this paper, we provide RNAseq datasets (raw and normalized counts as well as differentially expressed genes lists) of the transcriptome profiles of both grapevine berries (Vitis vinifera cv. Garganega) and B. cinerea during the establishment of noble rot, artificially induced in controlled conditions. The sequencing data are available in the NCBI GEO database under accession number GSE116741. These data were exploited in a comprehensive meta-analysis of gene expression during noble rot infection, gray mold and post-harvest withering. This highlighted an important common transcriptional reprogramming in different botrytized grape berry varieties and led to the identification of key genes specifically modulated during noble rot infection, which are described in the article entitled "Specific molecular interactions between Vitis vinifera and Botrytis cinerea are required for noble rot development in grape berries" Lovato et al., 2019.

Inverse Problems in Systems Biology: A Critical Review.

Systems Biology may be assimilated to a symbiotic cyclic interplaying between the forward and inverse problems. Computational models need to be continuously refined through experiments and in turn they help us to make limited experimental resources more efficient. Every time one does an experiment we know that there will be some noise that can disrupt our measurements. Despite the noise certainly is a problem, the inverse problems already involve the inference of missing information, even if the data is entirely reliable. So the addition of a certain limited noise does not fundamentally change the situation but can be used to solve the so-called ill-posed problem, as defined by Hadamard. It can be seen as an extra source of information. Recent studies have shown that complex systems, among others the systems biology, are poorly constrained and ill-conditioned because it is difficult to use experimental data to fully estimate their parameters. For these reasons was born the concept of sloppy models, a sequence of models of increasing complexity that become sloppy in the limit of microscopic accuracy. Furthermore the concept of sloppy models contains also the concept of un-identifiability, because the models are characterized by many parameters that are poorly constrained by experimental data. Then a strategy needs to be designed to infer, analyze, and understand biological systems. The aim of this work is to provide a critical review to the inverse problems in systems biology defining a strategy to determine the minimal set of information needed to overcome the problems arising from dynamic biological models that generally may have many unknown, non-measurable parameters.

Characterization of epithelial-mesenchymal transition intermediate/hybrid phenotypes associated to resistance to EGFR inhibitors in non-small cell lung cancer cell lines.

Increasing evidence points to a key role played by epithelial-mesenchymal transition (EMT) in cancer progression and drug resistance. In this study, we used wet and in silico approaches to investigate whether EMT phenotypes are associated to resistance to target therapy in a non-small cell lung cancer model system harboring activating mutations of the epidermal growth factor receptor. The combination of different analysis techniques allowed us to describe intermediate/hybrid and complete EMT phenotypes respectively in HCC827- and HCC4006-derived drug-resistant human cancer cell lines. Interestingly, intermediate/hybrid EMT phenotypes, a collective cell migration and increased stem-like ability associate to resistance to the epidermal growth factor receptor inhibitor, erlotinib, in HCC827 derived cell lines. Moreover, the use of three complementary approaches for gene expression analysis supported the identification of a small EMT-related gene list, which may have otherwise been overlooked by standard stand-alone methods for gene expression analysis.

The oncogenic role of circPVT1 in head and neck squamous cell carcinoma is mediated through the mutant p53/YAP/TEAD transcription-competent complex.

BACKGROUND: Circular RNAs are a class of endogenous RNAs with various functions in eukaryotic cells. Worthy of note, circular RNAs play a critical role in cancer. Currently, nothing is known about their role in head and neck squamous cell carcinoma (HNSCC). The identification of circular RNAs in HNSCC might become useful for diagnostic and therapeutic strategies in HNSCC. RESULTS: Using samples from 115 HNSCC patients, we find that circPVT1 is over-expressed in tumors compared to matched non-tumoral tissues, with particular enrichment in patients with TP53 mutations. circPVT1 up- and down-regulation determine, respectively, an increase and a reduction of the malignant phenotype in HNSCC cell lines. We show that circPVT1 expression is transcriptionally enhanced by the mut-p53/YAP/TEAD complex. circPVT1 acts as an oncogene modulating the expression of miR-497-5p and genes involved in the control of cell proliferation. CONCLUSIONS: This study shows the oncogenic role of circPVT1 in HNSCC, extending current knowledge about the role of circular RNAs in cancer.

Erratum: SWIM: a computational tool to unveiling crucial nodes in complex biological networks.

This corrects the article DOI: 10.1038/srep44797.

SWIM: a computational tool to unveiling crucial nodes in complex biological networks.

SWItchMiner (SWIM) is a wizard-like software implementation of a procedure, previously described, able to extract information contained in complex networks. Specifically, SWIM allows unearthing the existence of a new class of hubs, called "fight-club hubs", characterized by a marked negative correlation with their first nearest neighbors. Among them, a special subset of genes, called "switch genes", appears to be characterized by an unusual pattern of intra- and inter-module connections that confers them a crucial topological role, interestingly mirrored by the evidence of their clinic-biological relevance. Here, we applied SWIM to a large panel of cancer datasets from The Cancer Genome Atlas, in order to highlight switch genes that could be critically associated with the drastic changes in the physiological state of cells or tissues induced by the cancer development. We discovered that switch genes are found in all cancers we studied and they encompass protein coding genes and non-coding RNAs, recovering many known key cancer players but also many new potential biomarkers not yet characterized in cancer context. Furthermore, SWIM is amenable to detect switch genes in different organisms and cell conditions, with the potential to uncover important players in biologically relevant scenarios, including but not limited to human cancer.

Role of the long non-coding RNA PVT1 in the dysregulation of the ceRNA-ceRNA network in human breast cancer.

Recent findings have identified competing endogenous RNAs (ceRNAs) as the drivers in many disease conditions, including cancers. The ceRNAs indirectly regulate each other by reducing the amount of microRNAs (miRNAs) available to target messenger RNAs (mRNAs). The ceRNA interactions mediated by miRNAs are modulated by a titration mechanism, i.e. large changes in the ceRNA expression levels either overcome, or relieve, the miRNA repression on competing RNAs; similarly, a very large miRNA overexpression may abolish competition. The ceRNAs are also called miRNA "decoys" or miRNA "sponges" and encompass different RNAs competing with each other to attract miRNAs for interactions: mRNA, long non-coding RNAs (lncRNAs), pseudogenes, or circular RNAs. Recently, we developed a computational method for identifying ceRNA-ceRNA interactions in breast invasive carcinoma. We were interested in unveiling which lncRNAs could exert the ceRNA activity. We found a drastic rewiring in the cross-talks between ceRNAs from the physiological to the pathological condition. The main actor of this dysregulated lncRNA-associated ceRNA network was the lncRNA PVT1, which revealed a net biding preference towards the miR-200 family members in normal breast tissues. Despite its up-regulation in breast cancer tissues, mimicked by the miR-200 family members, PVT1 stops working as ceRNA in the cancerous state. The specific conditions required for a ceRNA landscape to occur are still far from being determined. Here, we emphasized the importance of the relative concentration of the ceRNAs, and their related miRNAs. In particular, we focused on the withdrawal in breast cancer tissues of the PVT1 ceRNA activity and performed a gene expression and sequence analysis of its multiple isoforms. We found that the PVT1 isoform harbouring the binding site for a representative miRNA of the miR-200 family shows a drastic decrease in its relative concentration with respect to the miRNA abundance in breast cancer tissues, providing a plausibility argument to the breakdown of the sponge program orchestrated by the oncogene PVT1.

Comprehensive RNA dataset of AGO2 associated RNAs in Jurkat cells following miR-21 over-expression.

We set out to identify miR-21 targets in Jurkat cells using a high-throughput biochemical approach (10.1016/j.biochi.2014.09.021[1]). Using a specific monoclonal antibody raised against AGO2, RISC complexes were immunopurified in Jurkat cells over-expressing miR-21 following lentiviral trasduction as well as in Jurkat control cells lines. A parallel immunoprecipitation using isotype-matched rat IgG was performed as a control. AGO2 associated mRNAs were profiled by microarray (GEO: GSE37212). AGO2 bound miRNAs were profiled by RNA-seq.

TP53 regulates miRNA association with AGO2 to remodel the miRNA-mRNA interaction network.

DNA damage activates TP53-regulated surveillance mechanisms that are crucial in suppressing tumorigenesis. TP53 orchestrates these responses directly by transcriptionally modulating genes, including microRNAs (miRNAs), and by regulating miRNA biogenesis through interacting with the DROSHA complex. However, whether the association between miRNAs and AGO2 is regulated following DNA damage is not yet known. Here, we show that, following DNA damage, TP53 interacts with AGO2 to induce or reduce AGO2's association of a subset of miRNAs, including multiple let-7 family members. Furthermore, we show that specific mutations in TP53 decrease rather than increase the association of let-7 family miRNAs, reducing their activity without preventing TP53 from interacting with AGO2. This is consistent with the oncogenic properties of these mutants. Using AGO2 RIP-seq and PAR-CLIP-seq, we show that the DNA damage-induced increase in binding of let-7 family members to the RISC complex is functional. We unambiguously determine the global miRNA-mRNA interaction networks involved in the DNA damage response, validating them through the identification of miRNA-target chimeras formed by endogenous ligation reactions. We find that the target complementary region of the let-7 seed tends to have highly fixed positions and more variable ones. Additionally, we observe that miRNAs, whose cellular abundance or differential association with AGO2 is regulated by TP53, are involved in an intricate network of regulatory feedback and feedforward circuits. TP53-mediated regulation of AGO2-miRNA interaction represents a new mechanism of miRNA regulation in carcinogenesis.

Identification of linc-NeD125, a novel long non coding RNA that hosts miR-125b-1 and negatively controls proliferation of human neuroblastoma cells.

The human genome contains some thousands of long non coding RNAs (lncRNAs). Many of these transcripts are presently considered crucial regulators of gene expression and functionally implicated in developmental processes in Eukaryotes. Notably, despite a huge number of lncRNAs are expressed in the Central Nervous System (CNS), only a few of them have been characterized in terms of molecular structure, gene expression regulation and function. In the present study, we identify linc-NeD125 as a novel cytoplasmic, neuronal-induced long intergenic non coding RNA (lincRNA). Linc-NeD125 represents the host gene for miR-125b-1, a microRNA with an established role as negative regulator of human neuroblastoma cell proliferation. Here, we demonstrate that these two overlapping non coding RNAs are coordinately induced during in vitro neuronal differentiation, and that their expression is regulated by different mechanisms. While the production of miR-125b-1 relies on transcriptional regulation, linc-NeD125 is controlled at the post-transcriptional level, through modulation of its stability. We also demonstrate that linc-NeD125 functions independently of the hosted microRNA, by reducing cell proliferation and activating the antiapoptotic factor BCL-2.

The role of TP53 in miRNA loading onto AGO2 and in remodelling the miRNA-mRNA interaction network.

BACKGROUND: DNA damage transactivates tumour protein p53 (TP53)-regulated surveillance, crucial in suppressing tumorigenesis. TP53 mediates this process directly by transcriptionally modulating gene and microRNA (miRNA) expression and indirectly by regulating miRNA biogenesis. However, the role of TP53 in regulating miRNA-AGO2 loading and global changes in AGO2 binding to its gene targets in response to DNA damage are unknown. These processes might be novel mechanisms by which TP53 regulates miRNAs in response to DNA damage. METHODS: To show the network of miRNA-mRNA interactions that occur in response to DNA damage, we stimulated TP53 wild-type and null cell-lines with doxorubicin and performed RNA sequencing from total RNA (RNA-Seq) and AGO2-immunoprecipitated RNA (AGO2-RIP-Seq). We used a combined AGO2 RIP-seq and AGO2 PAR-CLIP-seq (photoactivatable-ribonucleoside-enhanced cross-linking and immunoprecipitation) approach to determine the exact sites of interaction between the AGO2-bound miRNAs and their mRNA targets. FINDINGS: TP53 directly associated with AGO2, and induced and reduced loading of a subset of miRNAs, including the lethal 7 (let-7) miRNA family members, onto AGO2 in response to DNA damage. Although mutated TP53 maintained its capacity to interact with AGO2, it mediated unloading instead of loading of let-7 family miRNAs, thereby reducing their activity. We determined the miRNA-mRNA interaction networks involved in the response to DNA damage both in the presence and absence of TP53. Furthermore, we showed that miRNAs whose cellular abundance or differential loading onto AGO2 was regulated by TP53 were involved in an intricate network of regulatory feedback and feedforward circuits that fine-tuned gene expression levels in response to DNA damage to permit the repair of DNA damage or initiation of programmed cell death. INTERPRETATION: Control of AGO2 loading by TP53 is a new mechanism of miRNA regulation in carcinogenesis. FUNDING: UK Medical Research Council, Action Against Cancer.

Integrated molecular analysis to investigate the role of microRNAs in pancreatic tumour growth and progression.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNAs involved in the post-transcriptional regulation of mRNAs and are aberrantly expressed in cancer with important roles in tumorigenesis. A broad analysis of the combined effects of altered activities of miRNAs in pancreatic ductal adenocarcinoma (PDAC) has not been done, and how miRNAs might affect tumour progression or patient outcomes is unclear. METHODS: We combined data from miRNA and mRNA expression profiles from PDAC and normal pancreas samples (each n=9) and used bioinformatic analyses to identify a miRNA-mRNA regulatory network in PDAC. We validated our findings in PDAC cell-lines (PANC-1, MIA PaCa-2, LPc006, and LPc167), subcutaneous PDAC xenografts in mice, and laser capture microdissected PDACs from patients (n=91). We used this information to identify miRNAs that contributed most to tumorigenesis. FINDINGS: We identified three miRNAs (miR-21, miR-23a, and miR-27a) that acted as cooperative repressors of a network of tumour suppressor genes that included PDCD4, BTG2, and NEDD4L. Inhibition of miR-21, miR-23a, and miR-27a had synergistic effects in reducing proliferation of PDAC cells in culture and the growth of xenograft tumours. The level of inhibition was greater than that of silencing oncomiR-21 alone. In PDACs from patients, high levels of the combination of miR-21, miR-23a, and miR-27a was a strong independent predictor of short overall survival after surgical resection (hazard ratio 3·21, 95% CI 1·78-5·78). High expression of this combination was also associated with a more aggressive tumour phenotype: more microscopic tumour infiltration at resection margin and increased perineural invasion. INTERPRETATION: In an integrated data analysis, we identified functional miRNA-mRNA interactions that contribute to PDAC growth. These findings indicate that miRNAs act together to promote tumour progression and that future therapeutic strategies might require inhibition of several miRNAs. Furthermore, high tumour expression of the miR-21, miR-23a, and miR-27a combination could have potential use in the future as a prognostic signature for patients with PDAC. FUNDING: Peel Medical Research Trust, Alliance Family Foundation, Action Against Cancer, National Institute for Health Research, Association for International Cancer Research, Jason Boas Fellowship, Imperial Biomedical Research Centre, Rosetrees Trust, Joseph Ettedgui Charitable Foundation.

PVT1: a rising star among oncogenic long noncoding RNAs.

It is becoming increasingly clear that short and long noncoding RNAs critically participate in the regulation of cell growth, differentiation, and (mis)function. However, while the functional characterization of short non-coding RNAs has been reaching maturity, there is still a paucity of well characterized long noncoding RNAs, even though large studies in recent years are rapidly increasing the number of annotated ones. The long noncoding RNA PVT1 is encoded by a gene that has been long known since it resides in the well-known cancer risk region 8q24. However, a couple of accidental concurrent conditions have slowed down the study of this gene, that is, a preconception on the primacy of the protein-coding over noncoding RNAs and the prevalent interest in its neighbor MYC oncogene. Recent studies have brought PVT1 under the spotlight suggesting interesting models of functioning, such as competing endogenous RNA activity and regulation of protein stability of important oncogenes, primarily of the MYC oncogene. Despite some advancements in modelling the PVT1 role in cancer, there are many questions that remain unanswered concerning the precise molecular mechanisms underlying its functioning.