Dr Chinedu Anene

Background: Long non-coding RNAs (lncRNAs) are key regulators of cellular processes that underpin cancer development and progression. DRAIC is a migration inhibitor that has been linked with lung adenocarcinoma progression; however, its mechanisms remain to be studied. Methods: Several bioinformatics tools were used to explore the role of DRAIC in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). Results: Our bioinformatics analysis illustrates that patients with low expression of DRAIC have poor overall survival outcomes. In addition, the mRNA of SH3 domain-binding kinase 1 (SBK1) was downregulated in this cohort of patients. Mechanistic analysis showed that SBK1 is under the DRAIC competing endogenous RNAs network, potentially through sponging of miRNA-92a. Conclusions: Consistent dysregulation of the DRAIC-SBK1 axis was linked to poor survival outcome in both LUAD and LUSC, suggesting a tumour inhibitor role and providing potential for new diagnostics and therapeutic approaches.

Esophageal squamous cell carcinoma (ESCC) is highly heterogeneous. Our understanding of full molecular and immune landscape of ESCC remains limited, hindering the development of personalised therapeutic strategies. To address this, we perform genomic-transcriptomic characterizations and AI-aided histopathological image analysis of 120 Chinese ESCC patients. Here we show that ESCC can be categorized into differentiated, metabolic, immunogenic and stemness subtypes based on bulk and single-cell RNA-seq, each exhibiting specific molecular and histopathological features based on an amalgamated deep-learning model. The stemness subgroup with signature genes, such as WFDC2, SFRP1, LGR6 and VWA2, has the poorest prognosis and is associated with downregulated immune activities, a high frequency of EP300 mutation/activation, functional mutation enrichment in Wnt signalling and the highest level of intratumoural heterogeneity. The immune profiling by transcriptomics and immunohistochemistry reveals ESCC cells overexpress natural killer cell markers XCL1 and CD160 as immune evasion. Strikingly, XCL1 expression also affects the sensitivity of ESCC cells to common chemotherapy drugs. This study opens avenues for ESCC treatment and provides a valuable public resource to better understand ESCC.

The nucleus is a highly organised yet dynamic environment containing distinct membraneless nuclear bodies. This spatial separation enables a subset of components to be concentrated within biomolecular condensates, allowing efficient and discrete processes to occur which regulate cellular function. One such nuclear body, paraspeckles, are comprised of multiple paraspeckle proteins (PSPs) built around the architectural RNA, NEAT1_2. Paraspeckle function is yet to be fully elucidated but has been implicated in a variety of developmental and disease scenarios. We demonstrate that Kaposi’s sarcoma-associated herpesvirus (KSHV) drives formation of structurally distinct paraspeckles with a dramatically increased size and altered protein composition that are required for productive lytic replication. We highlight these virus-modified paraspeckles form adjacent to virus replication centres, potentially functioning as RNA processing hubs for viral transcripts during infection. Notably, we reveal that PSP sequestration into virus-modified paraspeckles result in increased genome instability during both KSHV and Epstein Barr virus (EBV) infection, implicating their formation in virus-mediated tumourigenesis.

BACKGROUND: Metastasis of cutaneous squamous cell carcinoma (cSCC) is uncommon. Current staging methods are reported to have sub-optimal performances in metastasis prediction. Accurate identification of patients with tumours at high risk of metastasis would have a significant impact on management. OBJECTIVE: To develop a robust and validated gene expression profile (GEP) signature for predicting primary cSCC metastatic risk using an unbiased whole transcriptome discovery-driven approach. METHODS: Archival formalin-fixed paraffin-embedded primary cSCC with perilesional normal tissue from 237 immunocompetent patients (151 non-metastasising and 86 metastasising) were collected retrospectively from four centres. TempO-seq was used to probe the whole transcriptome and machine learning algorithms were applied to derive predictive signatures, with a 3:1 split for training and testing datasets. RESULTS: A 20-gene prognostic model was developed and validated, with an accuracy of 86.0%, sensitivity of 85.7%, specificity of 86.1%, and positive predictive value of 78.3% in the testing set, providing more stable, accurate prediction than pathological staging systems. A linear predictor was also developed, significantly correlating with metastatic risk. LIMITATIONS: This was a retrospective 4-centre study and larger prospective multicentre studies are now required. CONCLUSION: The 20-gene signature prediction is accurate, with the potential to be incorporated into clinical workflows for cSCC. CONCLUSION: The 20-gene signature prediction is accurate, with the potential to be incorporated into clinical workflows for cSCC.

Actinic keratoses (AKs) are lesions of epidermal keratinocyte dysplasia and are precursors for invasive cutaneous squamous cell carcinoma (cSCC). Identifying the specific genomic alterations driving the progression from normal skin to skin with AK to skin with invasive cSCC is challenging because of the massive UVR-induced mutational burden characteristic at all stages of this progression. In this study, we report the largest AK whole-exome sequencing study to date and perform a mutational signature and candidate driver gene analysis on these lesions. We demonstrate in 37 AKs from both immunosuppressed and immunocompetent patients that there are significant similarities between AKs and cSCC in terms of mutational burden, copy number alterations, mutational signatures, and patterns of driver gene mutations. We identify 44 significantly mutated AK driver genes and confirm that these genes are similarly altered in cSCC. We identify azathioprine mutational signature in all AKs from patients exposed to the drug, providing further evidence for its role in keratinocyte carcinogenesis. cSCCs differ from AKs in having higher levels of intrasample heterogeneity. Alterations in signaling pathways also differ, with immune-related signaling and TGFβ signaling significantly more mutated in cSCC. Integrating our findings with independent gene expression datasets confirms that dysregulated TGFβ signaling may represent an important event in AK‒cSCC progression.

Human papillomaviruses (HPVs) cause most cervical cancers and an increasing number of anogenital and oral carcinomas, with most cases caused by HPV16 or HPV18. HPV hijacks host signalling pathways to promote carcinogenesis. Understanding these interactions could permit identification of much-needed therapeutics for HPV-driven malignancies. The Hippo signalling pathway is important in HPV+ cancers, with the downstream effector YAP playing a pro-oncogenic role. In contrast, the significance of its paralogue TAZ remains largely uncharacterised in these cancers. We demonstrate that TAZ is dysregulated in a HPV-type dependent manner by a distinct mechanism to that of YAP and controls proliferation via alternative cellular targets. Analysis of cervical cancer cell lines and patient biopsies revealed that TAZ expression was only significantly increased in HPV18+ and HPV18-like cells and TAZ knockdown reduced proliferation, migration and invasion only in HPV18+ cells. RNA-sequencing of HPV18+ cervical cells revealed that YAP and TAZ have distinct targets, suggesting they promote carcinogenesis by different mechanisms. Thus, in HPV18+ cancers, YAP and TAZ play non-redundant roles. This analysis identified TOGARAM2 as a previously uncharacterised TAZ target and demonstrates its role as a key effector of TAZ-mediated proliferation, migration and invasion in HPV18+ cancers.

Introduction Solid organ transplant recipients (SOTRs) are believed to have an increased risk of metastatic cutaneous squamous-cell carcinoma (cSCC), but reliable data are lacking regarding the precise incidence and associated risk factors. Methods In a prospective cohort study, including 19 specialist dermatology outpatient clinics in 15 countries, patient and tumor characteristics were collected using standardized questionnaires when SOTRs presented with a new cSCC. After a minimum of 2 years of follow-up, relevant data for all SOTRs were collected. Cumulative incidence of metastases was calculated by the Aalen-Johansen estimator. Fine and Gray models were used to assess multiple risk factors for metastases. Results Of 514 SOTRs who presented with 623 primary cSCCs, 37 developed metastases with a 2-year patient-based cumulative incidence of 6.2%. Risk factors for metastases included location in the head and neck area, local recurrence, size >2cm, clinical ulceration, poor differentiation grade, perineural invasion and deep invasion. A high-stage tumor that is also ulcerated showed the highest risk of metastasis, with a 2-year cumulative incidence of 46.2% (31.9% – 68.4%). Conclusions SOTRs have a high risk of cSCC metastases and well-established clinical and histological risk factors have been confirmed. High-stage, ulcerated cSCCs have the highest risk of metastasis.

Abstract

Non-coding RNA (ncRNA) regulatory networks are emerging as critical regulators of gene expression. These intricate networks of ncRNA:ncRNA interactions modulate multiple cellular pathways and impact the development and progression of multiple diseases. Herpesviruses, including Kaposi’s sarcoma-associated herpesvirus, are adept at utilising ncRNAs, encoding their own as well as dysregulating host ncRNAs to modulate virus gene expression and the host response to infection. Research has mainly focused on unidirectional ncRNA-mediated regulation of target protein-coding transcripts; however, we have identified a novel host ncRNA regulatory network essential for KSHV lytic replication in B cells. KSHV-mediated upregulation of the host cell circRNA, circHIPK3, is a key component of this network, functioning as a competing endogenous RNA of miR-30c, leading to increased levels of the miR-30c target, DLL4. Dysregulation of this network highlights a novel mechanism of cell cycle control during KSHV lytic replication in B cells. Importantly, disruption at any point within this novel ncRNA regulatory network has a detrimental effect on KSHV lytic replication, highlighting the essential nature of this network and potential for therapeutic intervention.

Non-coding RNA (ncRNA) regulatory networks are emerging as critical regulators of gene expression. These intricate networks of ncRNA:ncRNA interactions modulate multiple cellular pathways and impact the development and progression of multiple diseases. Herpesviruses, including Kaposi’s sarcoma-associated herpesvirus, are adept at utilising ncRNAs, encoding their own as well as dysregulating host ncRNAs to modulate virus gene expression and the host response to infection. Research has mainly focused on unidirectional ncRNA-mediated regulation of target protein-coding transcripts; however, we identify a novel host ncRNA regulatory network essential for KSHV lytic replication in B cells. KSHV-mediated upregulation of the host cell circRNA, circHIPK3, is a key component of this network, functioning as a competing endogenous RNA of miR-30c, leading to increased levels of the miR-30c target, DLL4. Dysregulation of this network highlights a novel mechanism of cell cycle control during KSHV lytic replication in B cells. Importantly, disruption at any point within this novel ncRNA regulatory network has a detrimental effect on KSHV lytic replication, highlighting the essential nature of this network and potential for therapeutic intervention.

Highlights •ACSNI uses deep neural network to decompose transcriptomic measurements •ACSNI maps pathway components to target biological or clinical phenotypes •ACSNI robustly extracts regulatory output of transcription factors •ACSNI enables annotation of genes with unknown functions The bigger picture Methods that group genes into functional units to quantify pathway activities are critical in the analysis of biological systems. Although many components of biological pathways have been described in detail, these tend to be limited to well-studied genes. In contrast, the majority of possible components remain unexplored. Here, we present a machine-learning tool for constructing and predicting tissue-specific components of biological pathways from large biological datasets. Our algorithm, ACSNI, can tackle incomplete pathway descriptions and enhance current pathway analysis methods' performance. We anticipate that, by dissecting the complex signals in biological data in a flexible and context-specific manner, ACSNI can facilitate the full characterization of physiological systems of interest. Summary Determining the tissue- and disease-specific circuit of biological pathways remains a fundamental goal of molecular biology. Many components of these biological pathways still remain unknown, hindering the full and accurate characterization of biological processes of interest. Here we describe ACSNI, an algorithm that combines prior knowledge of biological processes with a deep neural network to effectively decompose gene expression profiles (GEPs) into multi-variable pathway activities and identify unknown pathway components. Experiments on public GEP data show that ACSNI predicts cogent components of mTOR, ATF2, and HOTAIRM1 signaling that recapitulate regulatory information from genetic perturbation and transcription factor binding datasets. Our framework provides a fast and easy-to-use method to identify components of signaling pathways as a tool for molecular mechanism discovery and to prioritize genes for designing future targeted experiments (https://github.com/caanene1/ACSNI).

Abstract

Enhancer RNAs (eRNAs) are best known for their role in transcriptional regulation, where they facilitate enhancer-promoter communication and chromatin remodelling. Yet growing evidence suggests that their function may extend beyond the nucleus. Here, we systematically characterise the decay kinetics of eRNAs across human cell types using time-resolved transcriptomics and kinetic modelling. While most eRNAs undergo canonical exponential decay, a subset displays non-linear dynamics, suggesting context-dependent degradation mechanisms. Perturbation of core decay regulators, including components of the m 6 A and CCR4-NOT pathways, reveals that eRNA stability is modulated by a patchwork of pathways governing mRNA turnover. Integrating transcriptome-wide ribosome profiling, RNA-Seq, and half-life data, we identify eRNAs associated with changes in mRNA stability and translation efficiency of their target protein-coding transcripts. Functional validation of one such eRNA, en4528, shows it regulates CDKN2C mRNA independently of transcription and impacts cell migration. These findings redefine the regulatory scope of eRNAs, positioning them as active participants in post-transcriptional gene control and cellular behaviour. The resulting decay profiles and regulatory annotations have been incorporated into the eRNAkit database, available at https://github.com/AneneLab/eRNAkit , enhancing its capacity for integrative systems-level analysis of eRNA function.

Background: Follicular lymphoma (FL) is an incurable indolent B cell malignancy characterized in the majority of cases by the t(14;18) translocation. While the mutational landscape of the coding genome is nearing completion, less is known about the characteristics of its noncoding genome. Our expectation is that the distribution of noncoding mutations will be non-random, determined by epigenomics features, such as chromatin modification and accessibility. Our strategy is therefore to integrate whole genome epigenomic marks and mutations in order to enrich for variants with regulatory potential and identify unique mutational processes specific to thesecis-regulatory elements (CREs).Allele-specific expression (ASE) patterns allow further refinement to resolve functional CREs and bona fide mutations associated with changes in gene expression.

Methods: The H3K27Ac consensus CREs were determined using the ChIP-Seq data of 9 FL patient primary cells (Koues et al., Immunity2015). DNase I-hypersensitive sites (DHSs) and Hi-C of the B-lymphocyte cell line GM12878 was downloaded from ENCODE. Whole genome sequencing and RNA-Seq of 70 FL patients with relative high tumour cellularity (≥30% for DNA tumour purity and ≥25% for RNA B-cell content) were obtained from the International Cancer Genome Consortium project ICGC MMML-Seq. The mutation rate within CREs and DHSs in comparison to their flanking sequences was estimated using a previous pipeline (Sabarinathan et al., Nature2016). We further developed cis-ASE, an integrated analytic pipeline for the identification of recurrent ASE genes, and significantly associated CREs and mutations (Fig1.A).

Results: In total 1.04 million noncoding mutations, corresponding to 14.8K mutations per sample with a median of 9,991 noncoding mutations were identified in our series of 70 FL samples. 62K (6.0%) mutations were located within H3K27Ac bound CREs and there was an elevated mutation rate in H3K27Ac CREs compared to the corresponding flanking regions of 1kb up and downstream (χ² test p = 1.95e-19, Fig1.B). For DHSs, we observed the opposite pattern, with a lower mutation rate in DHSs than in the flanking regions (p = 0.04, Fig1.C), most likely reflecting the higher accessibility to global genome repair machinery in DHSs in relation to flanking sequences. Dividing mutations into high (≥0.3) and low (<0.3) adjusted variant allele frequency (VAF) groups (accounting for tumour purity), we observed significantly higher mutation rate in CREs than in the flanking sequences, that was specific for mutations with higher VAFs (p= 6.95e−38), as the difference was much weaker for low VAF mutations (p= 0.02). Mutation signatures 6 and 20, linked to defective DNA mismatch repair, were highly enriched for mutations within H3K27Ac CREs (44%) versus other regions outside (27%) (p< 0.001).

ASE was assessed using 45.6k informative SNPs (1.7%) per sample wherein matched genotype data was available from WGS and RNA-Seq profiles. Cis-ASE identified on average 480 ASE genes per sample (binomial test, adjusted p<0.05), corresponding to 1,943 recurrent ASE genes with a minimum threshold of ≥5 samples. These ASE events were not significantly influenced by local copy number changes or promoter methylation. KEGG pathway analysis of recurrent ASE genes identified adherens junction, B cell receptor signaling pathway and Fc gamma R-mediated phagocytosis as the most overrepresented pathways. Cis-ASE further identified 18 ASE-CRE interactions where CRE mutations were significantly correlated with an imbalance in elevated alternative allele ratios. These ASE genes included recognized lymphoma related genes including BCL2, STAT6, MAF, and additional novel targets. The pattern of aberrant somatic hypermutation (aSHM) was assessed for these 18 CREs, and we narrowed down to 15 significant ASE-CRE interactions not strongly affected by aSHM, consisting of 92 bona fide mutation candidates present in 37 FL samples of our cohort (52.9%).

Conclusion: Our study identified unique mutation processes operating in H3K27Ac CREs. Using an integrated genomic approach of whole genome mutations, chromatin marks and RNA-Seq, we explored ASE-CRE interactions, and identified 15 H3K27Ac bound CREs enriched for cis-regulatory mutations significantly associated with ASE and total expression of targeted genes, deserving for further exploration to enrich our understanding of FL noncoding genome.

Abstract

Prostate cancer (PCa) is genomically driven by dysregulation of transcriptional networks involving the transcriptional factors (TFs) FOXA1, ERG, AR, and HOXB13. However, the role of these specific TFs in the regulation of alternative pre-mRNA splicing (AS), which is a proven therapeutic vulnerability for cancers driven by the TF MYC, is not described. Using transcriptomic datasets from PCa patients, we tested for an association between expression of FOXA1, ERG, AR, HOXB13 , and MYC , and genes involved in AS - termed splicing-related proteins (SRPs), which have pleiotropic roles in RNA metabolism. We identified FOXA1 as the strongest predictor of dysregulated SRP gene expression, which was associated with PCa disease relapse after treatment. Subsequently, we selected a subset of FOXA1-binding and actively-transcribed SRP genes that phenocopy the FOXA1 dependency of PCa cells, and confirmed in vitro via knockdown and over-expression that FOXA1 regulates SRP gene expression. Finally, we demonstrated the persistence of a FOXA1-SRP gene association in treatment-relapsed castration-resistant PCa (CRPCa) patients. Our data demonstrate, for the first time, that FOXA1 controls dysregulated SRP gene expression, which is associated with poor PCa patient outcomes. Analogous to MYC-driven cancers, our findings implicate the therapeutic targeting of SRPs and AS in FOXA1-overexpressing PCa.

Abstract

Pseudouridylation is a prevalent RNA modification shown to occur in tRNAs, rRNAs, snoRNAs and most recently mRNAs and lncRNAs. Emerging evidence suggests that this dynamic RNA modification is implicated in altering gene expression by regulating RNA stability, modulating translation elongation and modifying amino acid substitution rates. However, the role of pseudouridylation in infection is poorly understood. Here we demonstrate that Kaposi’s sarcoma-associated herpesvirus (KSHV) manipulates the pseudouridylation pathway to enhance replication. We show the pseudouridine synthases (PUS), PUS1 and PUS7 are essential for efficient KSHV lytic replication, supported by the redistribution of both PUS1 and PUS7 to viral replication and transcription complexes. We present a comprehensive analysis of KSHV RNA pseudouridylation, revealing hundreds of modified RNAs at single-nucleotide resolution. Notably, we further demonstrate that pseudouridylation of the KSHV-encoded polyadenylated nuclear RNA (PAN) plays a significant role in the stability of PAN RNA and in the association of the KSHV ORF57 protein. Our findings reveal a novel and essential role of pseudouridine modification in the KSHV replication cycle.

The assessment of the cellular heterogeneity and abundance in bulk tissue samples is essential for characterising cellular and organismal states. Computational approaches to estimate cellular abundance from bulk RNA-Seq datasets have variable performances, often requiring benchmarking matrices to select the best performing methods for individual studies. However, such benchmarking investigations are difficult to perform and assess in typical applications because of the absence of gold standard/ground-truth cellular measurements. Here we describe Decosus, an R package that integrates seven methods and signatures for deconvoluting cell types from gene expression profiles (GEP). Benchmark analysis on a range of datasets with ground-truth measurements revealed that our integrated estimates consistently exhibited stable performances across datasets than individual methods and signatures. We further applied Decosus to characterise the immune compartment of skin samples in different settings, confirming the well-established Th1 and Th2 polarisation in psoriasis and atopic dermatitis, respectively. Secondly, we revealed immune system-related UV-induced changes in sun-exposed skin. Furthermore, a significant motivation in the design of Decosus is flexibility and the ability for the user to include new gene signatures, algorithms, and integration methods at run time.

Abstract

HNRNPA2B1 is associated with prostate cancer (PC) disease aggressiveness and underlies pro-tumourigenic cellular stress responses. By analysing >500 PC transcriptomes, we reveal that HNRNPA2B1 over-expression is associated with poor patient prognosis and stress response pathways. These include the “protein processing in the endoplasmic reticulum” (ER) pathway, which incorporates the unfolded protein response (UPR). By RNA-sequencing of HNRNPA2B1-depleted cells PC cells, we identified HNRNPA2B1-mediated down-regulation of UPR genes including the master ER-stress sensor IRE1 , which induces ER proteostasis. Consistent with IRE1 down-regulation in HNRNPA2B1-depleted cells, we observed reduced splicing of the IRE1-target and key UPR effector XBP1s. Furthermore, HNRNPA2B1 depletion up-regulates expression of the IRE1-dependent decay (RIDD) target gene BLOC1S1 , which is degraded by activated IRE1. We identify a HNRNPA2B1-IRE1-XBP1-controlled four gene prognostic biomarker signature (HIX) which classifies a subgroup of primary PC patients at high risk of disease relapse. Pharmacological targeting of IRE1 attenuated HNRNAPA2-driven PC cell growth. Taken together, our data reveal a putative novel mechanism of UPR activation in PC by HNRNPA2B1, which may promote an IRE1-dependent yet potentially-targetable recurrent disease phenotype.

Enhancer RNAs (eRNAs) are best known for their role in transcriptional regulation, where they facilitate enhancer-promoter communication and chromatin remodelling. Yet growing evidence suggests that their function may extend beyond the nucleus. Here, we systematically characterise the decay kinetics of eRNAs across human cell types using time-resolved transcriptomics and kinetic modelling. While most eRNAs undergo canonical exponential decay, a subset displays non-linear dynamics, suggesting context-dependent degradation mechanisms. Perturbation of core decay regulators, including components of the m⁶A and CCR4-NOT pathways, reveals that eRNA stability is modulated by a patchwork of pathways governing mRNA turnover. Integrating transcriptome-wide ribosome profiling, RNA-Seq, and half-life data, we identify eRNAs associated with changes in mRNA stability and translation efficiency of their target protein-coding transcripts. Functional validation of one such eRNA, en4528, shows it regulates CDKN2C mRNA independently of transcription and impacts cell migration. These findings redefine the regulatory scope of eRNAs, positioning them as active participants in post-transcriptional gene control and cellular behaviour. The resulting decay profiles and regulatory annotations have been incorporated into the eRNAkit database, available at https://github.com/AneneLab/eRNAkit, enhancing its capacity for integrative systems-level analysis of eRNA function.

Abstract

Noncoding variants and mutations outnumber their coding counterparts but remain challenging to interpret functionally. We present TranCi, a method that prioritises human genetic variations by integrating enhancer RNA (eRNA) expression with eRNA-mRNA interactome maps. By linking variant-associated changes in eRNA to downstream gene regulation, TranCi captures functional effects missed by sequence-based or chromatin-centric approaches. In esophageal squamous cell cancer, TranCi identifies noncoding mutations with roles in disease initiation and progression. A personalised mode enables analysis at single-patient resolution, uncovering potential individual-specific regulatory variants. TranCi thus provides a mechanistic framework for interpreting noncoding variations and uniquely identifies their downstream targets, where standard methods often fall short. TranCi is available as a module within the eRNAkit R package ( https://github.com/AneneLab/eRNAkit ), leveraging its database of eRNA expression and interactions for functional variant interpretation.

Motivation Recent RNA-centric experimental methods have significantly expanded our knowledge of proteins with known RNA-binding functions. However, the complete regulatory network and pathways for many of these RNA-binding proteins (RBPs) in different cellular contexts remain unknown. Although critical to understanding the role of RBPs in health and disease, experimentally mapping the RBP–RNA interactomes in every single context is an impossible task due the cost and manpower required. Additionally, identifying relevant RNAs bound by RBPs is challenging due to their diverse binding modes and function. Results To address these challenges, we developed RBP interaction mapper RBPInper an integrative framework that discovers global RBP interactome using statistical data fusion. Experiments on splicing factor proline and glutamine rich (SFPQ) datasets revealed cogent global SFPQ interactome. Several biological processes associated with this interactome were previously linked with SFPQ function. Furthermore, we conducted tests using independent dataset to assess the transferability of the SFPQ interactome to another context. The results demonstrated robust utility in generating interactomes that transfers to unseen cellular context. Overall, RBPInper is a fast and user-friendly method that enables a systems-level understanding of RBP functions by integrating multiple molecular datasets. The tool is designed with a focus on simplicity, minimal dependencies, and straightforward input requirements. This intentional design aims to empower everyday biologists, making it easy for them to incorporate the tool into their research.

Abstract

Enhancer RNAs (eRNAs) are a class of non-coding RNAs transcribed from active enhancers that regulate various aspects of transcription. Although traditionally viewed as nuclear-localised and unstable transcripts, large transcriptomic studies reveal they vary widely in localisation and biochemical properties, including detectable accumulation in cytoplasmic compartments. These observations suggest potential non-nuclear functions for eRNAs, yet existing databases remain focused on nuclear, cis-regulatory roles, limiting systematic exploration of their broader regulatory repertoire. Here we present eRNAkit , a comprehensive and accessible resource designed to address this gap by integrating subcellular localisation, RNA-RNA interaction and expression data for annotated eRNAs. Leveraging fractionation-based RNA-Seq datasets, eRNAkit profiles eRNA distribution across nuclear and cytoplasmic compartments. It incorporates gene expression data spanning major human organs and primary cell types, enabling tissue-specific analysis of eRNA function. Crucially, eRNAkit includes experimentally derived RNA-RNA interaction data from RIC-Seq, PARIS, and KARR-Seq, supporting exploration of trans-acting and cytoplasmic roles for eRNAs grounded in physical interaction evidence. eRNAkit expands current eRNA resources beyond the enhancer-promoter paradigm, offering a robust platform for dissecting non-canonical functions of eRNAs and advancing our understanding of their full regulatory potential in human biology. The eRNAkit resource is available to download at https://github.com/AneneLab/eRNAkit .

Exosomes play a crucial role in the crosstalk between cancer associated fibroblasts (CAFs) and cancer cells, contributing to carcinogenesis and the tumour microenvironment. Recent studies have revealed that CAFs, normal fibroblasts and cancer cells all secrete exosomes that contain miRNA, establishing a cell-cell communication network within the tumour microenvironment. For example, miRNA dysregulation in melanoma has been shown to promote CAF activation via induction of epithelial-mesenchymal transition (EMT), which in turn alters the secretory phenotype of CAFs in the stroma. This review assesses the roles of melanoma exosomal miRNAs in CAF formation and how CAF exosome-mediated feedback signalling to melanoma lead to tumour progression and metastasis. Moreover, efforts to exploit exosomal miRNA-mediated network communication between tumour cells and their microenvironment, and their potential as prognostic biomarkers or novel therapeutic targets in melanoma will also be considered.

Platelet microparticle (PMP)-induced angiogenesis plays a key role in tumour metastasis and has been proposed to contribute towards cardiovascular disease by enhancing atherosclerotic plaque vulnerability. However, the mechanisms underlying PMP induced angiogenesis are ill defined. Recent reports demonstrate that PMPs deliver micro-RNAs (miRNAs) to recipient cells, controlling gene expression. We therefore evaluated whether miRNA transfer was a key regulator of PMP-induced angiogenesis. Co-culturing PMPs with human umbilical vein endothelial cells (HUVEC) on extracellular matrix gel induced robust capillary like structure formation. PMP treatment altered the release of angiogenesis modulators from HUVEC, including significantly reducing production of anti-angiogenic thrombospondin-1 (THBS-1). Both functional responses were abrogated by treating PMPs with RNase, suggesting the transfer of PMP-derived RNA was a critical event. PMPs were an abundant source of miRNA Let-7a, which was transferred to HUVEC following co-incubation. Using luciferase reporter assays we have shown that Let-7a directly targets the 3’UTR of the THBS-1 mRNA. HUVEC transfection with a Let-7a anti-sense oligonucleotide reduced the ability of PMPs to inhibit THBS-1 release, and significantly decreased PMP induced in vitro angiogenesis. Antibody neutralisation of THBS-1 reversed the anti-angiogenic effect of let-7a inhibition in PMP treated HUVEC, highlighting Let-7a dependent translational repression of THBS-1 drives angiogenesis. Importantly, plasmid overexpression of Let-7a in HUVEC alone induced robust tubule formation on extracellular matrix gel. These data reveal a new role for Let-7a in promoting angiogenesis and show for the first time PMPs induced angiogenic responses occur through miRNA regulation of HUVEC.

The multifunctional protein, splicing factor, proline- and glutamine-rich (SFPQ) has been implicated in numerous cancers often due to interaction with coding and non-coding RNAs, however, its role in melanoma remains unclear. We report that knockdown of SFPQ expression in melanoma cells decelerates several cancer-associated cell phenotypes, including cell growth, migration, epithelial to mesenchymal transition, apoptosis, and glycolysis. RIP-seq analysis revealed that the SFPQ-RNA interactome is reprogrammed in melanoma cells and specifically enriched with key melanoma-associated coding and long non-coding transcripts, including SOX10, AMIGO2 and LINC00511 and in most cases SFPQ is required for the efficient expression of these genes. Functional analysis of two SFPQ-enriched lncRNA, LINC00511 and LINC01234, demonstrated that these genes independently contribute to the melanoma phenotype and a more detailed analysis of LINC00511 indicated that this occurs in part via modulation of the miR-625-5p/PKM2 axis. Importantly, analysis of a large clinical cohort revealed that elevated expression of SFPQ in primary melanoma tumours may have utility as a prognostic biomarker. Together, these data suggest that SFPQ is an important driver of melanoma, likely due to SFPQ-RNA interactions promoting the expression of numerous oncogenic transcripts.

The underlying causes of breast cancer are diverse, however, there is a striking association between type 2 diabetes and poor patient outcomes. Platelet activation is a common feature of both type 2 diabetes and breast cancer and has been implicated in tumourigenesis through a multitude of pathways. Here transcriptomic analysis of type 2 diabetes patient-derived platelet microvesicles revealed an altered miRNA signature compared with normoglycaemic control patients. Interestingly, interrogation of these data identifies a shift towards an oncogenic signature in type 2 diabetes-derived platelet microvesicles, with increased levels of miRNAs implicated in breast cancer progression and poor prognosis. Functional studies demonstrate that platelet microvesicles isolated from type 2 diabetes patient blood are internalised by triple-negative breast cancer cells in vitro, and that co-incubation with type 2 diabetes patient-derived platelet microvesicles led to significantly increased expression of epithelial to mesenchymal transition markers and triple-negative breast cancer cell invasion compared with platelet microvesicles from healthy volunteers. Together, these data suggest that circulating PMVs in type 2 diabetes patients may contribute to the progression of triple-negative breast cancer.

Abstract

Malignant melanoma has one of the lowest 5-year survival rates of any cancer, and is recognised for being particularly invasive and metastatic, with the poorest survival outcomes in brain metastases patients. A key characteristic of these tumours is crosstalk between melanoma cells and cells of the tumour microenvironment (TME), such as cancer associated fibroblasts (CAFs). The role of melanoma-derived small extracellular vesicles (sEVs) in potentiating CAFs has been studied extensively, however the role of CAF sEVs in regulation of the local TME and distal pre-metastatic niche (PMN) is less clear. Here, we demonstrate that sEVs derived from an in vitro model of melanoma CAFs alter melanoma to promote oncogenic parameters within models of the TME and target a model of the brain PMN to promote changes associated with melanoma extravasation. Cargo profiling of these sEVs found significant differential expression of proteins, and RNA associated with pre-metastatic niche remodelling and unfavourable outcomes in patients. Together these data suggest a role for CAF sEVs in local and distal PMN formation, highlighting a potential therapeutic target for metastatic melanoma and identifying prospective liquid biomarker reservoirs.

Malignant melanoma has one of the lowest 5‐year survival rates of any cancer and is characterised by its high invasiveness and metastatic potential, with especially poor outcomes in patients who develop brain metastases. Crosstalk between melanoma cells and cells of the tumour microenvironment (TME), including cancer‐associated fibroblasts (CAFs), is a central driver of disease progression. While the role of melanoma‐derived small extracellular vesicles (sEVs) in reprogramming stromal cells has been well documented, the reciprocal effects of CAF‐derived sEVs remain less clear. Here, using an in vitro model of melanoma CAFs, we show that CAF sEVs alter melanoma cells and fibroblasts to promote oncogenic traits and remodel endothelial cells, including brain microvascular cells, in ways consistent with early pre‐metastatic niche (PMN) changes. Multi‐omics cargo profiling revealed significant differential expression of proteins and RNAs linked to extracellular matrix organisation, vascular remodelling, and patient outcomes, with functional validation identifying THBS1 as an EV cargo that restrains endothelial sprouting while potentially promoting barrier destabilisation. Together, these findings suggest that CAF‐derived sEVs contribute to local and distal PMN remodelling, highlight their potential as therapeutic targets, and identify EV cargoes with promise as circulating biomarkers in melanoma.

Merkel cell carcinoma (MCC) is an aggressive skin cancer with high rates of recurrence and metastasis. Merkel cell polyomavirus (MCPyV) is associated with the majority of MCC cases. MCPyV-induced tumourigenesis is largely dependent on the expression of the small tumour antigen (ST). Recent findings implicate MCPyV ST expression in the highly metastatic nature of MCC by promoting cell motility and migration, through differential expression of cellular proteins that lead to microtubule destabilisation, filopodium formation and breakdown of cell–cell junctions. However, the molecular mechanisms which dysregulate these cellular processes are yet to be fully elucidated. Here, we demonstrate that MCPyV ST expression activates p38 MAPK signalling to drive cell migration and motility. Notably, MCPyV ST-mediated p38 MAPK signalling occurs through MKK4, as opposed to the canonical MKK3/6 signalling pathway. In addition, our results indicate that an interaction between MCPyV ST and the cellular phospatase subunit PP4C is essential for its effect on p38 MAPK signalling. These results provide novel opportunities for the treatment of metastatic MCC given the intense interest in p38 MAPK inhibitors as therapeutic agents.