The physiology of alternative splicing

DIPAN: Detecting personalized intronic polyadenylation derived neoantigens from RNA sequencing data

Intronic polyadenylation (IPA) refers to a particular type of alternative polyadenylation where a gene makes use of a polyadenylation site located within its introns. Aberrant IPA events have been observed in various types of cancer. IPA can produce noncoding transcripts or truncated protein-coding transcripts with altered coding sequences in the resulting protein product. Therefore, IPA events hold the potential to act as a reservoir of tumor neoantigens. Here, we developed a computational method termed DIPAN, which incorporates IPA detection, protein fragmentation, and MHC binding prediction to predict IPA-derived neoantigens. Utilizing RNA-seq from breast cancer cell lines and ovarian cancer clinical samples, we demonstrated the significant contribution of IPA events to the neoantigen repertoire. Through mass spectrometry immunopeptidome analysis, we further illustrated the processing and presentation of IPA-derived neoantigens on the surface of cancer cells. While most IPA-derived neoantigens are sample-specific, shared neoantigens were identified in both cancer cell lines and clinical samples. Furthermore, we demonstrated an association between IPA-derived neoantigen burden and overall survival in cancer patients.

A novel splice variant in intron 10 of PEX6 is associated with Zellweger Syndrome in a Chinese neonate

BACKGROUND

Zellweger Syndrome (ZS), or cerebrohepatorenal syndrome, is a rare disorder due to PEX gene mutations affecting peroxisome function. While PEX6 coding mutations are known to cause ZS, the impact of noncoding mutations is less clear.

METHODS

A Chinese neonate and his family were subjected to whole exome sequencing (WES) and bioinformatics to assess variant pathogenicity. A minigene assay was also performed for detailed splicing variant analysis.

RESULTS

WES identified compound heterozygous PEX6 variants: c.315G>A (p. Trp105Ter) and c.2095-3 T>G. Minigene assays indicated that the latter variant led to abnormal mRNA splicing and the loss of exon 11 in PEX6 expression, potentially causing nonsense-mediated mRNA decay (NMD) or truncated protein structure.

CONCLUSION

The study suggests that PEX6: c.2095-3 T>G might be a genetic contributor to the patient's condition, broadening the known mutation spectrum of PEX6. These insights lay groundwork for potential gene therapy for such variants.

Striving for clarity in language about gene expression

What do we mean when we say 'gene expression'? In the decades following Crick's 1958 central dogma of molecular biology, whereby genetic information flows from DNA (genes) to RNA (transcripts) to protein (products), we have learned a great deal about DNA, RNA, proteins, and the ensuing phenotypic changes. With the advent of high-throughput technologies (1990s), molecular biologists and computer scientists forged critical collaborations to understand the vast amount of data being generated, rapidly escalating gene expression research to the 'omics' level: entire sets of genes (genomes), transcribed RNAs (transcriptomes), and synthesized proteins (proteomes). However, some concessions came to be made for molecular biologists and computer scientists to understand each other-one of the most prevalent being the increasingly widespread use of 'gene' to mean 'RNAs originating from a DNA segment'. This loosening of terminology, we will argue, creates ambiguity and confusion. We propose guidelines to increase precision and clarity when communicating about gene expression, most notably to reserve 'gene' for the DNA template and 'transcript' for the RNA transcribed from that gene. Striving to use perspicuous terminology will promote rigorous gene expression science and accelerate discovery in this highly promising area of biology.

Alternative splicing and intron retention: Their profiles and roles in cutaneous fibrosis of systemic sclerosis

BACKGROUND

Alternative splicing (AS) and intron retention (IR) implicated in multiple pathophysiological processes, have rarely been reported in systemic sclerosis (SSc).

METHODS

We integrated bulk RNA-seq and 4D label-free mass spectrometry to perform a multi-omics analysis of AS and IR in SSc skin tissue and fibroblasts. RMATS and iREAD were used to identify AS and IR, which were validated by real-time PCR. Spearman correlation and the LASSO method were employed to assess correlations among clinical features, introns, splicing factors (regulators of AS) and proteins.

FINDINGS

AS profiles showed distinct alterations in SSc skin tissue, with the most pronounced changes occurring in IR. AS and IR were associated with total modified Rodnan skin score (mRSS) and local skin score. Upon TGF-β stimulation, fibroblasts exhibited significant alterations in IR profiles, affecting genes related to fibroblast proliferation and collagen fibril organization. A comprehensive integrated analysis of introns, exons, and proteome profiles revealed that IR exerted a negative impact on protein expression, with certain changes being under intronic control. RT-PCR confirmed the presence of intron and exon-derived sequences of CTTN, OGA, MED16 and PHYKPL. Additionally, notable changes were observed in the regulatory network of splicing factors in SSc skin tissues. These factors are also involved in fibrosis pathways and correlated with clinical features.

CONCLUSION

Totally, abnormal AS, IR profiles and splicing factors were identified in SSc, altered IRs and splicing factors participated in fibrosis-related pathways. IR exerted a negative impact on protein expression in TGF-β-stimulated fibroblasts. Clarification of the IR mechanisms will provide new insights into the pathophysiology of SSc.

The Exon Junction Complex component EIF4A3 plays a splicing-linked oncogenic role in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers, underscoring the urgent need for in-depth biological research. The phenomenon of alternative RNA splicing dysregulation is a common hallmark in cancer, including PDAC, presenting new avenues for understanding and developing diagnostic and therapeutic tools. Our research focuses on EIF4A3, a core component of the Exon Junction Complex intimately linked to RNA splicing, and its role in PDAC. EIF4A3 is overexpressed in PDAC tissue and associated to clinical parameters of malignancy and poorer patient survival. Mechanistically, exploration of PDAC RNA-seq data unveiled the link of EIF4A3 to diverse malignancy processes, consistent with its association to key molecular pathways. EIF4A3 targeting in vitro decreased essential functional tumor features such as proliferation, migration, colony formation and sphere formation, while its in vivo targeting reduced tumor growth. EIF4A3 silencing in PDAC cell lines severely altered its transcriptional and spliceosomic landscapes, as shown by RNA-seq analyses, suggesting a role for EIF4A3 in maintaining RNA homeostasis. Our results indicate that EIF4A3 dysregulation in PDAC has a pleiotropic regulatory role on RNA biology, influencing key cellular functions. This paves the way to explore its potential as novel biomarker and actionable target candidate for this lethal cancer.

Circular RNA in cancer

Over the past decade, circular RNA (circRNA) research has evolved into a bona fide research field shedding light on the functional consequence of this unique family of RNA molecules in cancer. Although the method of formation and the abundance of circRNAs can differ from their cognate linear mRNA, the spectrum of interacting partners and their resultant cellular functions in oncogenesis are analogous. However, with 10 times more diversity in circRNA variants compared with linear RNA variants, combined with their hyperstability in the cell, circRNAs are equipped to influence every stage of oncogenesis. This is an opportune time to address the breadth of circRNA in cancer focused on their spatiotemporal expression, mutations in biogenesis factors and contemporary functions through each stage of cancer. In this Review, we highlight examples of functional circRNAs in specific cancers, which satisfy critical criteria, including their physical co-association with the target and circRNA abundance at stoichiometrically valid quantities. These considerations are essential to develop strategies for the therapeutic exploitation of circRNAs as biomarkers and targeted anticancer agents.

Deep learning and direct sequencing of labeled RNA captures transcriptome dynamics

In eukaryotes, genes produce a variety of distinct RNA isoforms, each with potentially unique protein products, coding potential or regulatory signals such as poly(A) tail and nucleotide modifications. Assessing the kinetics of RNA isoform metabolism, such as transcription and decay rates, is essential for unraveling gene regulation. However, it is currently impeded by lack of methods that can differentiate between individual isoforms. Here, we introduce RNAkinet, a deep convolutional and recurrent neural network, to detect nascent RNA molecules following metabolic labeling with the nucleoside analog 5-ethynyl uridine and long-read, direct RNA sequencing with nanopores. RNAkinet processes electrical signals from nanopore sequencing directly and distinguishes nascent from pre-existing RNA molecules. Our results show that RNAkinet prediction performance generalizes in various cell types and organisms and can be used to quantify RNA isoform half-lives. RNAkinet is expected to enable the identification of the kinetic parameters of RNA isoforms and to facilitate studies of RNA metabolism and the regulatory elements that influence it.

RNA splicing factor RBFOX2 is a key factor in the progression of cancer and cardiomyopathy

BACKGROUND

Alternative splicing of pre-mRNA is a fundamental regulatory process in multicellular eukaryotes, significantly contributing to the diversification of the human proteome. RNA-binding fox-1 homologue 2 (RBFOX2), a member of the evolutionarily conserved RBFOX family, has emerged as a critical splicing regulator, playing a pivotal role in the alternative splicing of pre-mRNA. This review provides a comprehensive analysis of RBFOX2, elucidating its splicing activity through direct and indirect binding mechanisms. RBFOX2 exerts substantial influence over the alternative splicing of numerous transcripts, thereby shaping essential cellular processes such as differentiation and development.

MAIN BODY OF THE ABSTRACT

Dysregulation of RBFOX2-mediated alternative splicing has been closely linked to a spectrum of cardiovascular diseases and malignant tumours, underscoring its potential as a therapeutic target. Despite significant progress, current research faces notable challenges. The complete structural characterisation of RBFOX2 remains elusive, limiting in-depth exploration beyond its RNA-recognition motif. Furthermore, the scarcity of studies focusing on RBFOX2-targeting drugs poses a hindrance to translating research findings into clinical applications.

CONCLUSION

This review critically assesses the existing body of knowledge on RBFOX2, highlighting research gaps and limitations. By delineating these areas, this analysis not only serves as a foundational reference for future studies but also provides strategic insights for bridging these gaps. Addressing these challenges will be instrumental in unlocking the full therapeutic potential of RBFOX2, paving the way for innovative and effective treatments in various diseases.

Metabolic regulation of mRNA splicing

Alternative mRNA splicing enables the diversification of the proteome from a static genome and confers plasticity and adaptiveness on cells. Although this is often explored in development, where hard-wired programs drive the differentiation and specialization, alternative mRNA splicing also offers a way for cells to react to sudden changes in outside stimuli such as small-molecule metabolites. Fluctuations in metabolite levels and availability in particular convey crucial information to which cells react and adapt. We summarize and highlight findings surrounding the metabolic regulation of mRNA splicing. We discuss the principles underlying the biochemistry and biophysical properties of mRNA splicing, and propose how these could intersect with metabolite levels. Further, we present examples in which metabolites directly influence RNA-binding proteins and splicing factors. We also discuss the interplay between alternative mRNA splicing and metabolite-responsive signaling pathways. We hope to inspire future research to obtain a holistic picture of alternative mRNA splicing in response to metabolic cues.

Postsynaptic lncRNA Sera/Pkm2 pathway orchestrates the transition from social competition to rank by remodeling the neural ensemble in mPFC

Individuals' continuous success in competitive interactions with conspecifics strongly affects their social hierarchy. Medial prefrontal cortex (mPFC) is the key brain region mediating both social competition and hierarchy. However, the molecular regulatory mechanisms underlying the neural ensemble in the mPFC remains unclear. Here, we demonstrate that in excitatory neurons of prelimbic cortex (PL), lncRNA Sera remodels the utilization of Pkm Exon9 and Exon10, resulting in a decrease in the Pkm1/2 ratio in highly competitive mice. By employing a tet-on/off system, we disrupt or rebuild the normal Pkm1/2 ratio by controlling the expression of Pkm2 in PL excitatory neurons. We find that long-term Pkm2 modulation induces timely competition alteration and hysteretic rank change, through phosphorylating the Ser845 site of GluA1. Together, this study uncovers a crucial role of lncRNA Sera/Pkm2 pathway in the transition of social competition to rank by remodeling neural ensemble in mPFC.

Enhancing novel isoform discovery: leveraging nanopore long-read sequencing and machine learning approaches

Long-read sequencing technologies can capture entire RNA transcripts in a single sequencing read, reducing the ambiguity in constructing and quantifying transcript models in comparison to more common and earlier methods, such as short-read sequencing. Recent improvements in the accuracy of long-read sequencing technologies have expanded the scope for novel splice isoform detection and have also enabled a far more accurate reconstruction of complex splicing patterns and transcriptomes. Additionally, the incorporation and advancements of machine learning and deep learning algorithms in bioinformatic software have significantly improved the reliability of long-read sequencing transcriptomic studies. However, there is a lack of consensus on what bioinformatic tools and pipelines produce the most precise and consistent results. Thus, this review aims to discuss and compare the performance of available methods for novel isoform discovery with long-read sequencing technologies, with 25 tools being presented. Furthermore, this review intends to demonstrate the need for developing standard analytical pipelines, tools, and transcript model conventions for novel isoform discovery and transcriptomic studies.

The Functional Relationship Between RNA Splicing and the Chromatin Landscape

Chromatin is a highly regulated and dynamic structure that has been shown to play an essential role in transcriptional and co-transcriptional regulation. In the context of RNA splicing, early evidence suggested a loose connection between the chromatin landscape and splicing. More recently, it has been shown that splicing occurs in a co-transcriptional manner, meaning that the splicing process occurs in the context of chromatin. Experimental and computational evidence have also shown that chromatin dynamics can influence the splicing process and vice versa. However, much of this evidence provides mainly correlative relationships between chromatin and splicing with just a few concrete examples providing defined molecular mechanisms by which these two processes are functionally related. Nevertheless, it is clear that chromatin and RNA splicing are tightly interconnected to one another. In this review, we highlight the current state of knowledge of the relationship between chromatin and splicing.

Translatome analysis in acute ischemic stroke: Astrocytes and microglia exhibit differences in poststroke alternative splicing of expressed transcripts

Astrocytes and microglia undergo dynamic and complex morphological and functional changes following ischemic stroke, which are instrumental in both inflammatory responses and neural repair. While gene expression alterations poststroke have been extensively studied, investigations into posttranscriptional regulatory mechanisms, specifically alternative splicing (AS), remain limited. Utilizing previously reported Ribo-Tag-seq data, this study analyzed AS alterations in poststroke astrocytes and microglia from young adult male and female mice. Our findings reveal that in astrocytes, compared to the sham group, 109 differential alternative splicing (DAS) events were observed at 4 h poststroke, which increased to 320 at day 3. In microglia, these numbers were 316 and 266, respectively. Interestingly, the disparity between DAS genes and differentially expressed genes is substantial, with fewer than 10 genes shared at both poststroke time points in astrocytes and microglia. Gene ontology enrichment analysis revealed the involvement of these DAS genes in diverse functions, encompassing immune response (Adam8, Ccr1), metabolism (Acsl6, Pcyt2, Myo5a), and developmental cell growth (App), among others. Selective DAS events were further validated by semiquantitative RT-PCR. Overall, this study comprehensively describes the AS alterations in astrocytes and microglia during the hyperacute and acute phases of ischemic stroke and underscores the significance of certain hub DAS events in neuroinflammatory processes.

The Role of Sirtuin-1 Isoforms in Regulating Mitochondrial Function

The sirtuin-1 (SIRT1) gene contains multiple exons that usually undergo alternative splicing. The exclusion of one or more exons causes domain loss in the alternatively spliced isoforms and may change their functions. However, it is not completely established to what extent the loss of a non-catalytic domain could affect its regulatory function. Using muscle cells and SIRT1-knockout cells, we examined the function of the constitutively spliced isoform (SIRT1-v1) versus the alternatively spliced isoforms SIRT1-v2 and SIRT1-v3 that had lost part of the N-terminal region. Our data indicate that partial loss of the N-terminal domains in SIRT1-v2 and SIRT1-v3 attenuated their function. The full-length SIRT1-v1 significantly increased the oxidative phosphorylation and ATP production rate. Furthermore, SIRT1-v1 specifically upregulated the mitochondrial respiratory complex I without affecting the activity of complexes II, III, and IV. Additionally, domain loss affected the regulation of site-specific lysine acetylation in the histone H4 protein, the gene expression of respiratory complex I subunits, and the metabolic balance of oxidative phosphorylation versus glycolysis. Since alternatively spliced isoforms tend to increase with advancing age, the impact of SIRT1 isoforms on mitochondrial respiratory complexes warrants further investigation.

Principles of clinical genetics for rheumatologists: clinical indications and interpretation of broad-based genetic testing

Advances in DNA sequencing technologies, especially next-generation sequencing (NGS), which is the basis for whole-exome sequencing (WES) and whole-genome sequencing (WGS), have profoundly transformed immune-mediated rheumatic disease diagnosis. Recently, substantial cost reductions have facilitated access to these diagnostic tools, expanded the capacity of molecular diagnostics and enabled the pursuit of precision medicine in rheumatology. Understanding the fundamental principles of genetics and diversity in genetic variant classification is a crucial milestone in rheumatology. However, despite the growing availability of DNA sequencing platforms, a significant number of autoinflammatory diseases (AIDs), neuromuscular disorders, hereditary collagen diseases, and monogenic bone diseases remain unsolved, and variants of uncertain significance (VUS) pose a formidable challenge to addressing these unmet needs in the coming decades. This article aims to provide an overview of the clinical indications and interpretation of comprehensive genetic testing in the medical field, addressing the related complexities and implications.

Co-delivery Nano System of MS-275 and V-9302 Induces Pyroptosis and Enhances Anti-Tumor Immunity Against Uveal Melanoma

In the treatment of uveal melanoma (UVM), histone deacetylase inhibitors (HDACi) have emerged as a promising epigenetic therapy. However, their clinical efficacy is hindered by the suboptimal pharmacokinetics and the strong self-rescue of tumor cells. To overcome these limitations, reactive oxygen species (ROS)-responsive nanoparticles (NPs) are designed that encapsulate HDACi MS-275 and the glutamine metabolism inhibitor V-9302. Upon reaching the tumor microenvironment, these NPs can disintegrate, thereby releasing MS-275 to increase the level of ROS and V-9302 to reduce the production of glutathione (GSH) related to self-rescue. These synergistic effects lead to a lethal ROS storm and induce cell pyroptosis. When combined with programmed cell death protein 1 monoclonal antibodies (α-PD-1), these NPs facilitate immune cell infiltration, improving anti-tumor immunity, converting "immune-cold" tumors into "immune-hot" tumors, and enhancing immune memory in mice. The findings present a nano-delivery strategy for the co-delivery of epigenetic therapeutics and metabolic inhibitors, which induces pyroptosis in tumors cells and improves the effectiveness of chemotherapy and immunotherapy.

Geneapp: A web application for visualizing alternative splicing for biomedicine

Alternative Splicing (AS) is an essential mechanism for eukaryotes. However, the consequences of deleting a single exon can be dramatic for the organism and can lead to cancer in humans. Additionally, alternative 5' and 3' splice sites, which define the boundaries of exons, also play key roles to human disorders. Therefore, Investigating AS events is crucial for understanding the molecular basis of human diseases and developing therapeutic strategies. Workflow for AS event analysis can be sampling followed by data analysis with bioinformatics to identify the different AS events in the control and case samples, data visualization for curation, and selection of relevant targets for experimental validation. The raw output of the analysis software does not favor the inspection of events by bioinformaticians requiring custom scripts for data visualization. In this work, we propose the Geneapp application with three modules: GeneappScript, GeneappServer, and GeneappExplorer. GeneappScript is a wrapper that assists in identifying AS in samples compared in two different approaches, while GeneappServer integrates data from AS analysis already performed by the user. In GeneappExplorer, the user visualizes the previous dataset by exploring AS events in genes with functional annotation. This targeted screens that Geneapp allows to perform helps in the identification of targets for experimental validation to confirm the hypotheses under study. The Geneapp is freely available for non-commercial use at https://geneapp.net to advance research on AS for bioinformatics.

Can geroscience be translated into healthcare?

As an introduction to this special issue on geroscience, the present work focuses on the complexity of disentangling biomolecular mechanisms of aging from biopsychosocial causes of accelerated aging. Due to this complexity, the biomolecular aging hallmarks of frailty and multimorbidity as predominant aging phenotypes in geriatrics reflect single aspects of the aging process. A possible approach to facilitate the integration of geroscience into healthcare might be to consider aging as the dynamic ratio between damage accumulation at the molecular and cellular level and resilience as strategies that prevent or repair such damage. There is a large body of evidence to show that geroscience has the potential to change healthcare; however, reaching a consensus and translating the best tool to measure aging needs more research on 1) the sensitivity of biomarkers to interventions and 2) the relationship between changes in biomarkers and changes in health trajectories.

Kinase activities in pancreatic ductal adenocarcinoma with prognostic and therapeutic avenues

Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with a limited number of known driver mutations but considerable cancer cell heterogeneity. Phosphoproteomics provides a direct read-out of aberrant signaling and the resultant clinically relevant phenotype. Mass spectrometry (MS)-based proteomics and phosphoproteomics were applied to 42 PDAC tumors. Data encompassed over 19 936 phosphoserine or phosphothreonine (pS/T; in 5412 phosphoproteins) and 1208 phosphotyrosine (pY; in 501 phosphoproteins) sites and a total of 3756 proteins. Proteome data identified three distinct subtypes with tumor intrinsic and stromal features. Subsequently, three phospho-subtypes were apparent: two tumor intrinsic (Phos1/2) and one stromal (Phos3), resembling known PDAC molecular subtypes. Kinase activity was analyzed by the Integrative iNferred Kinase Activity (INKA) scoring. Phospho-subtypes displayed differential phosphorylation signals and kinase activity, such as FGR and GSK3 activation in Phos1, SRC kinase family and EPHA2 in Phos2, and EGFR, INSR, MET, ABL1, HIPK1, JAK, and PRKCD in Phos3. Kinase activity analysis of an external PDAC cohort supported our findings and underscored the importance of PI3K/AKT and ERK pathways, among others. Interestingly, unfavorable patient prognosis correlated with higher RTK, PAK2, STK10, and CDK7 activity and high proliferation, whereas long survival was associated with MYLK and PTK6 activity, which was previously unknown. Subtype-associated activity profiles can guide therapeutic combination approaches in tumor and stroma-enriched tissues, and emphasize the critical role of parallel signaling pathways. In addition, kinase activity profiling identifies potential disease markers with prognostic significance.

High-throughput sensitive screening of small molecule modulators of microexon alternative splicing using dual Nano and Firefly luciferase reporters

Disruption of alternative splicing frequently causes or contributes to human diseases and disorders. Consequently, there is a need for efficient and sensitive reporter assays capable of screening chemical libraries for compounds with efficacy in modulating important splicing events. Here, we describe a screening workflow employing dual Nano and Firefly luciferase alternative splicing reporters that affords efficient, sensitive, and linear detection of small molecule responses. Applying this system to a screen of ~95,000 small molecules identified compounds that stimulate or repress the splicing of neuronal microexons, a class of alternative exons often disrupted in autism and activated in neuroendocrine cancers. One of these compounds rescues the splicing of several analyzed microexons in the cerebral cortex of an autism mouse model haploinsufficient for Srrm4, a major activator of brain microexons. We thus describe a broadly applicable high-throughput screening system for identifying candidate splicing therapeutics, and a resource of small molecule modulators of microexons with potential for further development in correcting aberrant splicing patterns linked to human disorders and disease.

KPNA2 promotes the progression of gastric cancer by regulating the alternative splicing of related genes

RNA-binding proteins (RBPs) play critical roles in genome regulation. In this study, we explored the latent function of KPNA2, which is an essential member of the RBP family, in the regulation of alternative splicing (AS) in gastric cancer (GC). We analyzed the role of KPNA2 in regulating differential expression and AS via RNA sequencing (RNA-seq) and improved RNA immunoprecipitation sequencing (iRIP-seq). Clinical specimens were used to analyze the associations between KPNA2 expression and clinicopathological characteristics. CCK8 assays, transwell assays and wound healing assays were performed to explore the effect of KPNA2/WDR62 on GC cell progression. KPNA2 was shown to be highly expressed in GC cells and tissues and associated with lymph node metastases. KPNA2 promoted the proliferation, migration and invasion of GC cells and primarily regulated exon skipping, alternative 3's splice sites (A3SSs), alternative 5' splice sites (A5SSs), and cassette exons. We further revealed that KPNA2 participated in biological processes related to cell proliferation, and the immune response in GC via the regulation of transcription. In addition, KPNA2 preferentially bound to intron regions. Notably, KPNA2 regulated the A3SS AS mode of WDR62, and upregulation of WDR62 reversed the KPNA2 downregulation-induced inhibition of GC cell proliferation, migration and invasion. Finally, we discovered that the AS of immune-related molecules could be regulated by KPNA2. Overall, our results demonstrated for the first time that KPNA2 functions as an oncogenic splicing factor in GC that regulated the AS and differential expression of GC-related genes, and KPNA2 may be a potential target for GC treatment.

RNA structure in alternative splicing regulation: from mechanism to therapy

Alternative splicing is a highly intricate process that plays a crucial role in post-transcriptional regulation and significantly expands the functional proteome of a limited number of coding genes in eukaryotes. Its regulation is multifactorial, with RNA structure exerting a significant impact. Aberrant RNA conformations lead to dysregulation of splicing patterns, which directly affects the manifestation of disease symptoms. In this review, the molecular mechanisms of RNA secondary structure-mediated splicing regulation are summarized, with a focus on the complex interplay between aberrant RNA conformations and disease phenotypes resulted from splicing defects. This study also explores additional factors that reshape structural conformations, enriching our understanding of the mechanistic network underlying structure-mediated splicing regulation. In addition, an emphasis has been placed on the clinical role of targeting aberrant splicing corrections in human diseases. The principal mechanisms of action behind this phenomenon are described, followed by a discussion of prospective development strategies and pertinent challenges.

Novel prognostic alternative splicing events in colorectal Cancer: Impact on immune infiltration and therapy response

OBJECTIVE

This study aims to comprehensively analyze alternative splicing (AS) features in colorectal cancer (CRC) using integrative multi-omics and to elucidate their relationship with the CRC immune microenvironment.

METHODS

Transcriptomic data, clinical information, and Percent Spliced In (PSI) values of AS events for CRC patients were obtained from The Cancer Genome Atlas (TCGA) and TCGA SpliceSeq databases. Differentially expressed AS events were identified. Univariate Cox analysis was used to pinpoint prognosis-related AS events. A prognostic risk model was developed and validated using multivariate Cox analysis, patient survival analysis, and the area under the receiver operating characteristic (ROC) curve (AUC). Gene Set Enrichment Analysis (GSEA), immune infiltration, immunotherapy, chemotherapy sensitivity analyses, and regulatory relationships between AS events and splicing factors (SFs) were conducted. Single-cell sequencing was used to study the distribution of key factors. siRNA and overexpression vectors were utilized to silence/overexpress BCAS1 in CRC cells and evaluate their effects on cell growth, migration, and invasion. Furthermore, the interaction between BCAS1 and ANO7 pre-mRNA was investigated using RIP-PCR.

RESULTS

82 prognosis-related AS events were identified in CRC patients. A 15-AS prognostic model was constructed, which correlated with immune cell infiltration and showed differences in immunotherapy and chemotherapy sensitivity. BCAS1 was identified as a potential regulator of the ANO7|58341|AT splicing event in CRC. Single-cell sequencing analysis revealed the distribution of BCAS1 and ANO7 in cancer stem cells. In vitro experiments demonstrated that overexpression of BCAS1 and silencing of ANO7 inhibit the proliferation, migration, and invasion of CRC cells. Moreover, BCAS1 suppresses the progression of CRC by modulating ANO7 alternative splicing.

CONCLUSION

This study provides new insights into the role of alternative splicing in colorectal cancer, particularly the potential regulatory action of BCAS1 on the ANO7|58341|AT splicing event. It also identifies the impact of alternative splicing on the tumor microenvironment and potential implications for immunotherapy, highlighting its relevance for the in-depth study and treatment of CRC.

Zinc finger BED-type containing 3 promotes hepatic steatosis by interacting with polypyrimidine tract-binding protein 1

AIMS/HYPOTHESIS

The relationship between metabolic dysfunction-associated steatotic liver disease (MASLD) and type 2 diabetes mellitus, insulin resistance and the metabolic syndrome is well established. While zinc finger BED-type containing 3 (ZBED3) has been linked to type 2 diabetes mellitus and the metabolic syndrome, its role in MASLD remains unclear. In this study, we aimed to investigate the function of ZBED3 in the context of MASLD.

METHODS

Expression levels of ZBED3 were assessed in individuals with MASLD, as well as in cellular and animal models of MASLD. In vitro and in vivo analyses were conducted using a cellular model of MASLD induced by NEFA and an animal model of MASLD induced by a high-fat diet (HFD), respectively, to investigate the role of ZBED3 in MASLD. ZBED3 expression was increased by lentiviral infection or tail-vein injection of adeno-associated virus. RNA-seq and bioinformatics analysis were employed to examine the pathways through which ZBED3 modulates lipid accumulation. Findings from these next-generation transcriptome sequencing studies indicated that ZBED3 controls SREBP1c (also known as SREBF1; a gene involved in fatty acid de novo synthesis); thus, co-immunoprecipitation and LC-MS/MS were utilised to investigate the molecular mechanisms by which ZBED3 regulates the sterol regulatory element binding protein 1c (SREBP1c).

RESULTS

In this study, we found that ZBED3 was significantly upregulated in the liver of individuals with MASLD and in MASLD animal models. ZBED3 overexpression promoted NEFA-induced triglyceride accumulation in hepatocytes in vitro. Furthermore, the hepatocyte-specific overexpression of Zbed3 promoted hepatic steatosis. Conversely, the hepatocyte-specific knockout of Zbed3 resulted in resistance of HFD-induced hepatic steatosis. Mechanistically, ZBED3 interacts directly with polypyrimidine tract-binding protein 1 (PTBP1) and affects its binding to the SREBP1c mRNA precursor to regulate SREBP1c mRNA stability and alternative splicing.

CONCLUSIONS/INTERPRETATION

This study indicates that ZBED3 promotes hepatic steatosis and serves as a critical regulator of the progression of MASLD.

DATA AVAILABILITY

RNA-seq data have been deposited in the NCBI Gene Expression Omnibus ( www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE231875 ). MS proteomics data have been deposited to the ProteomeXchange Consortium via the iProX partner repository ( https://proteomecentral.proteomexchange.org/cgi/GetDataset?ID=PXD041743 ).

The Effect of Alternative Splicing Sites on Mirtron Formation and Arm Selection of Precursor microRNAs

Mirtrons represent a subclass of microRNAs (miRNAs) that rely on the splicing machinery for their maturation. However, the molecular details of this Drosha-independent processing are still not fully understood; as an example, the Microprocessor complex cannot process the mirtronic pre-miRNA from the transcript even if splice site mutations are present. To investigate the influence of alternative splicing sites on mirtron formation, we generated Enhanced Green Fluorescent Protein (EGFP) reporters containing artificial introns to compare the processing of canonical miRNAs and mirtrons. Although mutations of both splice sites generated a complex pattern of alternative transcripts, mirtron formation was always severely affected as opposed to the normal processing of the canonical hsa-mir-33b miRNA. However, we also detected that while its formation was also hindered, the mirtron-derived hsa-mir-877-3p miRNA was less affected by certain mutations than the hsa-mir-877-5p species. By knocking down Drosha, we showed that this phenomenon is not dependent on Microprocessor activity but rather points toward the potential stability difference between the miRNAs from the different arms. Our results indicate that when the major splice sites are mutated, mirtron formation cannot be rescued by nearby alternative splice sites, and stability differences between 5p and 3p species should also be considered for functional studies of mirtrons.

Efficient, specific, and combinatorial control of endogenous exon splicing with dCasRx-RBM25

Efficient targeted control of splicing is a major goal of functional genomics and therapeutic applications. Guide (g)RNA-directed, deactivated (d)Cas CRISPR enzymes fused to splicing effectors represent a promising strategy due to the flexibility of these systems. However, efficient, specific, and generalizable activation of endogenous exons using this approach has not been previously reported. By screening over 300 dCasRx-splicing factor fusion proteins tethered to splicing reporters, we identify dCasRx-RBM25 as a potent activator of exons. Moreover, dCasRx-RBM25 efficiently activates the splicing of ∼90% of targeted endogenous alternative exons and displays high on-target specificity. Using gRNA arrays for combinatorial targeting, we demonstrate that dCasRx-RBM25 enables multiplexed activation and repression of exons. Using this feature, the targeting of neural-regulated exons in Ptpb1 and Puf60 in embryonic stem cells reveals combinatorial effects on downstream alternative splicing events controlled by these factors. Collectively, our results enable versatile, combinatorial exon-resolution functional assays and splicing-directed therapeutic applications.

Genome-scale exon perturbation screens uncover exons critical for cell fitness

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

Neoantigens in cancer immunotherapy: focusing on alternative splicing

Alternative splicing (AS) functions as a crucial program in transcriptional modulation, leading to proteomic diversity and functional alterations of proteins. These splicing actions induce various neoantigens that hold prognostic significance and contribute to various aspects of cancer progression, including immune responses against cancer. The advent of immunotherapy has remarkably revolutionized tumor therapy. In this regard, AS-derived neoantigens are potent targets for cancer vaccines and chimeric antigen receptor (CAR) T cell therapies. In this review, we outline that AS-derived neoantigens serve as promising immunotherapeutic targets and guide immunotherapy strategies. This evidence contributes to a deeper comprehension of the complexity of proteomic diversity and provides novel perspectives and techniques for precision medicine in immunotherapy. Moreover, we underscore the obstacles that are awaited to be addressed for this novel approach to become clinically applicable.

Development and disease-specific regulation of RNA splicing in cardiovascular system

Alternative splicing is a complex gene regulatory process that distinguishes itself from canonical splicing by rearranging the introns and exons of an immature pre-mRNA transcript. This process plays a vital role in enhancing transcriptomic and proteomic diversity from the genome. Alternative splicing has emerged as a pivotal mechanism governing complex biological processes during both heart development and the development of cardiovascular diseases. Multiple alternative splicing factors are involved in a synergistic or antagonistic manner in the regulation of important genes in relevant physiological processes. Notably, circular RNAs have only recently garnered attention for their tissue-specific expression patterns and regulatory functions. This resurgence of interest has prompted a reevaluation of the topic. Here, we provide an overview of our current understanding of alternative splicing mechanisms and the regulatory roles of alternative splicing factors in cardiovascular development and pathological process of different cardiovascular diseases, including cardiomyopathy, myocardial infarction, heart failure and atherosclerosis.

ESRP2-microRNA-122 axis directs the postnatal onset of liver polyploidization and maturation

Hepatocyte polyploidy and maturity are critical to acquiring specialized liver functions. Multiple intra- and extracellular factors influence ploidy, but how they cooperate temporally to steer liver polyploidization and maturation or how post-transcriptional mechanisms integrate into these paradigms is unknown. Here, we identified an important regulatory hierarchy in which postnatal activation of Epithelial-Splicing-Regulatory-Protein-2 (ESRP2) stimulates biogenesis of liver-specific microRNA (miR-122), thereby facilitating polyploidization, maturation, and functional competence of hepatocytes. By determining transcriptome-wide protein-RNA interactions in vivo and integrating them with single-cell and bulk hepatocyte RNA-seq datasets, we delineate an ESRP2-driven RNA processing program that drives sequential replacement of fetal-to-adult transcript isoforms. Specifically, ESRP2 binds the primary miR-122 host gene transcript to promote its processing/biogenesis. Combining constitutive and inducible ESRP2 gain- and loss-of-function mice models with miR-122 rescue experiments, we demonstrate that timed activation of ESRP2 augments miR-122-driven program of cytokinesis failure, ensuring proper onset and extent of hepatocyte polyploidization.

Epigenetic control of skeletal muscle atrophy

Skeletal muscular atrophy is a complex disease involving a large number of gene expression regulatory networks and various biological processes. Despite extensive research on this topic, its underlying mechanisms remain elusive, and effective therapeutic approaches are yet to be established. Recent studies have shown that epigenetics play an important role in regulating skeletal muscle atrophy, influencing the expression of numerous genes associated with this condition through the addition or removal of certain chemical modifications at the molecular level. This review article comprehensively summarizes the different types of modifications to DNA, histones, RNA, and their known regulators. We also discuss how epigenetic modifications change during the process of skeletal muscle atrophy, the molecular mechanisms by which epigenetic regulatory proteins control skeletal muscle atrophy, and assess their translational potential. The role of epigenetics on muscle stem cells is also highlighted. In addition, we propose that alternative splicing interacts with epigenetic mechanisms to regulate skeletal muscle mass, offering a novel perspective that enhances our understanding of epigenetic inheritance's role and the regulatory network governing skeletal muscle atrophy. Collectively, advancements in the understanding of epigenetic mechanisms provide invaluable insights into the study of skeletal muscle atrophy. Moreover, this knowledge paves the way for identifying new avenues for the development of more effective therapeutic strategies and pharmaceutical interventions.

Chemical Proteomic Profiling of Protein Dopaminylation in Colorectal Cancer Cells

Histone dopaminylation is a newly identified epigenetic mark that plays a role in the regulation of gene transcription, where an isopeptide bond is formed between the fifth amino acid of H3 (i.e., glutamine) and dopamine. Recently, we developed a chemical probe to specifically label and enrich histone dopaminylation via bioorthogonal chemistry. Given this powerful tool, we found that histone H3 glutamine 5 dopaminylation (H3Q5dop) was highly enriched in colorectal tumors, which could be attributed to the high expression level of its regulator, transglutaminase 2 (TGM2), in colon cancer cells. Due to the enzyme promiscuity of TGM2, nonhistone proteins have also been identified as dopaminylation targets; however, the dopaminylated proteome in cancer cells still remains elusive. Here, we utilized our chemical probe to enrich dopaminylated proteins from colorectal cancer cells in a bioorthogonal manner and performed the chemical proteomics analysis. Therefore, 425 dopaminylated proteins were identified, many of which are involved in nucleic acid metabolism and transcription pathways. More importantly, a number of dopaminylation sites were identified and attributed to the successful application of our chemical probe. Overall, these findings shed light on the significant association between cellular protein dopaminylation and cancer development, further suggesting that targeting these pathways may become a promising anticancer strategy.

IsoVis - a webserver for visualization and annotation of alternative RNA isoforms

Genes commonly express multiple RNA products (RNA isoforms), which differ in exonic content and can have different functions. Making sense of the plethora of known and novel RNA isoforms being identified by transcriptomic approaches requires a user-friendly way to visualize gene isoforms and how they differ in exonic content, expression levels and potential functions. Here we introduce IsoVis, a freely available webserver that accepts user-supplied transcriptomic data and visualizes the expressed isoforms in a clear, intuitive manner. IsoVis contains numerous features, including the ability to visualize all RNA isoforms of a gene and their expression levels; the annotation of known isoforms from external databases; mapping of protein domains and features to exons, allowing changes to protein sequence and function between isoforms to be established; and extensive species compatibility. Datasets visualised on IsoVis remain private to the user, allowing analysis of sensitive data. IsoVis visualisations can be downloaded to create publication-ready figures. The IsoVis webserver enables researchers to perform isoform analyses without requiring programming skills, is free to use, and available at https://isomix.org/isovis/.

Functions of methyltransferase-like 3 in breast cancer: pathogenesis, drug resistance, and therapeutic target

Breast cancer (BC) is a highly prevalent malignancy worldwide, with complex pathogenesis and treatment challenges. Research reveals that methyltransferase-like 3 (METTL3) is widely involved in the pathogenesis of several tumors through methylation of its target RNAs, and its role and mechanisms in BC are also extensively studied. In this review, we aim to provide a comprehensive interpretation of available studies and elucidate the relationship between METTL3 and BC. This review suggests that high levels of METTL3 are associated with the pathogenesis, poor prognosis, and drug resistance of BC, suggesting METTL3 as a potential diagnostic or prognostic biomarker and therapeutic target. Collectively, this review provides a comprehensive understanding of how METTL3 functions through RNA methylation, which provides a valuable reference for future fundamental studies and clinical applications.

Expression of Truncated Products at the 5'-Terminal Region of RIPK2 and Evolutive Aspects that Support Their Biological Importance

Alternative splicing is the process of generating different mRNAs from the same primary transcript, which contributes to increase the transcriptome and proteome diversity. Abnormal splicing has been associated with the development of several diseases including cancer. Given that mutations and abnormal levels of the RIPK2 transcript and RIP-2 protein are frequent in tumors, and that RIP-2 modulates immune and inflammatory responses, we investigated alternative splicing events that result in partial deletions of the kinase domain at the N-terminus of RIP-2. We also investigated the structure and expression of the RIPK2 truncated variants and isoforms in different environments. In addition, we searched data throughout Supraprimates evolution that could support the biological importance of RIPK2 alternatively spliced products. We observed that human variants and isoforms were differentially regulated following temperature stress, and that the truncated transcript was more expressed than the long transcript in tumor samples. The inverse was found for the longer protein isoform. The truncated variant was also detected in chimpanzee, gorilla, hare, pika, mouse, rat, and tree shrew. The fact that the same variant has been preserved in mammals with divergence times up to 70 million years raises the hypothesis that it may have a functional significance.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Histamine H3 Receptor Isoforms: Insights from Alternative Splicing to Functional Complexity

Alternative splicing significantly enhances the diversity of the G protein-coupled receptor (GPCR) family, including the histamine H3 receptor (H3R). This post-transcriptional modification generates multiple H3R isoforms with potentially distinct pharmacological and physiological profiles. H3R is primarily involved in the presynaptic inhibition of neurotransmitter release in the central nervous system. Despite the approval of pitolisant for narcolepsy (Wakix®) and daytime sleepiness in adults with obstructive sleep apnea (Ozawade®) and ongoing clinical trials for other H3R antagonists/inverse agonists, the functional significance of the numerous H3R isoforms remains largely enigmatic. Recent publicly available RNA sequencing data have confirmed the expression of multiple H3R isoforms in the brain, with some isoforms exhibiting unique tissue-specific distribution patterns hinting at isoform-specific functions and interactions within neural circuits. In this review, we discuss the complexity of H3R isoforms with a focus on their potential roles in central nervous system (CNS) function. Comparative analysis across species highlights evolutionary conservation and divergence in H3R splicing, suggesting species-specific regulatory mechanisms. Understanding the functionality of H3R isoforms is crucial for the development of targeted therapeutics. This knowledge will inform the design of more precise pharmacological interventions, potentially enhancing therapeutic efficacy and reducing adverse effects in the treatment of neurological and psychiatric disorders.

Long-read single-cell sequencing reveals the transcriptional landscape of spermatogenesis in obstructive azoospermia and Sertoli cell-only patients

BACKGROUND

High-throughput single-cell RNA sequencing (scRNA-seq) is widely used in spermatogenesis. However, it only reveals short reads in germ and somatic cells, limiting the discovery of novel transcripts and genes.

AIM

This study shows the long-read transcriptional landscape of spermatogenesis in obstructive azoospermia (OA) and Sertoli cell-only patients.

DESIGN

Single cells were isolated from testicular biopsies of OA and non-obstructive azoospermia (NOA) patients. Cell culture was identified by comparing PacBio long-read single-cell sequencing (OA n = 3, NOA n = 3) with short-read scRNA-seq (OA n = 6, NOA n = 6). Ten germ cell types and eight somatic cell types were classified based on known markers.

METHODS

PacBio long-read single-cell sequencing, short-read scRNA-seq, polymerase chain reaction.

RESULTS

A total of 130 426 long-read transcripts (100 517 novel transcripts and 29 909 known transcripts) and 49 508 long-read transcripts (26 002 novel transcripts and 23 506 known transcripts) have been detected in OA and NOA patients, respectively. Moreover, 36 373 and 1642 new genes are identified in OA and NOA patients, respectively. Importantly, specific expressions of long-read transcripts were detected in germ and stomatic cells during normal spermatogenesis.

CONCLUSION

We have identified total full-length transcripts in OA and NOA, and new genes were found. Furthermore, specific expressed full-length transcripts were detected, and the genomic structure of transcripts was mapped in different cell types. These findings may provide valuable information on human spermatogenesis and the treatment of male infertility.

Alternative splicing coupled to nonsense-mediated decay coordinates downregulation of non-neuronal genes in developing mouse neurons

BACKGROUND

The functional coupling between alternative pre-mRNA splicing (AS) and the mRNA quality control mechanism called nonsense-mediated decay (NMD) can modulate transcript abundance. Previous studies have identified several examples of such a regulation in developing neurons. However, the systems-level effects of AS-NMD in this context are poorly understood.

RESULTS

We developed an R package, factR2, which offers a comprehensive suite of AS-NMD analysis functions. Using this tool, we conducted a longitudinal analysis of gene expression in pluripotent stem cells undergoing induced neuronal differentiation. Our analysis uncovers hundreds of AS-NMD events with significant potential to regulate gene expression. Notably, this regulation is significantly overrepresented in specific functional groups of developmentally downregulated genes. Particularly strong association with gene downregulation is detected for alternative cassette exons stimulating NMD upon their inclusion into mature mRNA. By combining bioinformatic analyses with CRISPR/Cas9 genome editing and other experimental approaches we show that NMD-stimulating cassette exons regulated by the RNA-binding protein PTBP1 dampen the expression of their genes in developing neurons. We also provided evidence that the inclusion of NMD-stimulating cassette exons into mature mRNAs is temporally coordinated with NMD-independent gene repression mechanisms.

CONCLUSIONS

Our study provides an accessible workflow for the discovery and prioritization of AS-NMD targets. It further argues that the AS-NMD pathway plays a widespread role in developing neurons by facilitating the downregulation of functionally related non-neuronal genes.

Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

Exploring the role of m 6 A writer RBM15 in cancer: a systematic review

In the contemporary epoch, cancer stands as the predominant cause of premature global mortality, necessitating a focused exploration of molecular markers and advanced therapeutic strategies. N6-methyladenosine (m6A), the most prevalent mRNA modification, undergoes dynamic regulation by enzymes referred to as methyltransferases (writers), demethylases (erasers), and effective proteins (readers). Despite lacking methylation activity, RNA-binding motif protein 15 (RBM15), a member of the m6A writer family, assumes a crucial role in recruiting the methyltransferase complex (MTC) and binding to mRNA. Although the impact of m6A modifications on cancer has garnered widespread attention, RBM15 has been relatively overlooked. This review briefly outlines the structure and operational mechanism, and delineates the unique role of RBM15 in various cancers, shedding light on its molecular basis and providing a groundwork for potential tumor-targeted therapies.

Long-read sequencing for brain tumors

Brain tumors and genomics have a long-standing history given that glioblastoma was the first cancer studied by the cancer genome atlas. The numerous and continuous advances through the decades in sequencing technologies have aided in the advanced molecular characterization of brain tumors for diagnosis, prognosis, and treatment. Since the implementation of molecular biomarkers by the WHO CNS in 2016, the genomics of brain tumors has been integrated into diagnostic criteria. Long-read sequencing, also known as third generation sequencing, is an emerging technique that allows for the sequencing of longer DNA segments leading to improved detection of structural variants and epigenetics. These capabilities are opening a way for better characterization of brain tumors. Here, we present a comprehensive summary of the state of the art of third-generation sequencing in the application for brain tumor diagnosis, prognosis, and treatment. We discuss the advantages and potential new implementations of long-read sequencing into clinical paradigms for neuro-oncology patients.

RNA-clique: a method for computing genetic distances from RNA-seq data

BACKGROUND

Although RNA-seq data are traditionally used for quantifying gene expression levels, the same data could be useful in an integrated approach to compute genetic distances as well. Challenges to using mRNA sequences for computing genetic distances include the relatively high conservation of coding sequences and the presence of paralogous and, in some species, homeologous genes.

RESULTS

We developed a new computational method, RNA-clique, for calculating genetic distances using assembled RNA-seq data and assessed the efficacy of the method using biological and simulated data. The method employs reciprocal BLASTn followed by graph-based filtering to ensure that only orthologous genes are compared. Each vertex in the graph constructed for filtering represents a gene in a specific sample under comparison, and an edge connects a pair of vertices if the genes they represent are best matches for each other in their respective samples. The distance computation is a function of the BLAST alignment statistics and the constructed graph and incorporates only those genes that are present in some complete connected component of this graph. As a biological testbed we used RNA-seq data of tall fescue (Lolium arundinaceum), an allohexaploid plant ( 2 n = 14 Gb ), and bluehead wrasse (Thalassoma bifasciatum), a teleost fish. RNA-clique reliably distinguished individual tall fescue plants by genotype and distinguished bluehead wrasse RNA-seq samples by individual. In tests with simulated RNA-seq data, the ground truth phylogeny was accurately recovered from the computed distances. Moreover, tests of the algorithm parameters indicated that, even with stringent filtering for orthologs, sufficient sequence data were retained for the distance computations. Although comparisons with an alternative method revealed that RNA-clique has relatively high time and memory requirements, the comparisons also showed that RNA-clique's results were at least as reliable as the alternative's for tall fescue data and were much more reliable for the bluehead wrasse data.

CONCLUSION

Results of this work indicate that RNA-clique works well as a way of deriving genetic distances from RNA-seq data, thus providing a methodological integration of functional and genetic diversity studies.

Shiba: A unified computational method for robust identification of differential RNA splicing across platforms

Alternative pre-mRNA splicing (AS) is a fundamental regulatory process that generates transcript diversity and cell type variation. We developed Shiba, a robust method integrating transcript assembly, splicing event identification, read counting, and statistical analysis, to efficiently quantify exon splicing levels across various types of RNA-seq datasets. Compared to existing pipelines, Shiba excels in capturing both annotated and unannotated or cryptic differential splicing events with superior accuracy, sensitivity, and reproducibility. Furthermore, Shiba's unique consideration of junction read imbalance and exon-body read coverage reduces false positives, essential for downstream functional analyses. We have further developed scShiba for single-cell/nucleus (sc/sn) RNA-seq data, enabling the exploration of splicing variations in heterogeneous cell populations. Both simulated and real data demonstrate Shiba's robustness across multiple sample sizes, including n=1 datasets and individual cell clusters from scRNA-seq. Application of Shiba on single replicates of RNA-seq identified new AS-NMD targets, and scShiba on snRNA-seq revealed intricate temporal AS regulation in dopaminergic neurons. Both Shiba and scShiba are provided in Docker/Singularity containers and Snakemake pipeline, enhancing accessibility and reproducibility. The comprehensive capabilities of Shiba and scShiba allow systematic and robust quantification of alternative splicing events, laying a solid foundation for mechanistic exploration of functional complexity in RNA splicing.

Alternative splicing in EMT and TGF-β signaling during cancer progression

Epithelial to mesenchymal transition (EMT) is a physiological process during development where epithelial cells transform to acquire mesenchymal characteristics, which allows them to migrate and colonize secondary tissues. Many cellular signaling pathways and master transcriptional factors exert a myriad of controls to fine tune this vital process to meet various developmental and physiological needs. Adding to the complexity of this network are post-transcriptional and post-translational regulations. Among them, alternative splicing has been shown to play important roles to drive EMT-associated phenotypic changes, including actin cytoskeleton remodeling, cell-cell junction changes, cell motility and invasiveness. In advanced cancers, transforming growth factor-β (TGF-β) is a major inducer of EMT and is associated with tumor cell metastasis, cancer stem cell self-renewal, and drug resistance. This review aims to provide an overview of recent discoveries regarding alternative splicing events and the involvement of splicing factors in the EMT and TGF-β signaling. It will emphasize the importance of various splicing factors involved in EMT and explore their regulatory mechanisms.

Conserved role of hnRNPL in alternative splicing of epigenetic modifiers enables B cell activation

The multifunctional RNA-binding protein hnRNPL is implicated in antibody class switching but its broader function in B cells is unknown. Here, we show that hnRNPL is essential for B cell activation, germinal center formation, and antibody responses. Upon activation, hnRNPL-deficient B cells show proliferation defects and increased apoptosis. Comparative analysis of RNA-seq data from activated B cells and another eight hnRNPL-depleted cell types reveals common effects on MYC and E2F transcriptional programs required for proliferation. Notably, while individual gene expression changes are cell type specific, several alternative splicing events affecting histone modifiers like KDM6A and SIRT1, are conserved across cell types. Moreover, hnRNPL-deficient B cells show global changes in H3K27me3 and H3K9ac. Epigenetic dysregulation after hnRNPL loss could underlie differential gene expression and upregulation of lncRNAs, and explain common and cell type-specific phenotypes, such as dysfunctional mitochondria and ROS overproduction in mouse B cells. Thus, hnRNPL is essential for the resting-to-activated B cell transition by regulating transcriptional programs and metabolism, at least in part through the alternative splicing of several histone modifiers.

Alternative isoform expression of key thermogenic genes in human beige adipocytes

Background

The beneficial effect of thermogenic adipocytes in maintaining body weight and protecting against metabolic disorders has raised interest in understanding the regulatory mechanisms defining white and beige adipocyte identity. Although alternative splicing has been shown to propagate adipose browning signals in mice, this has yet to be thoroughly investigated in human adipocytes.

Methods

We performed parallel white and beige adipogenic differentiation using primary adipose stem cells from 6 unrelated healthy subjects and assessed differential gene and isoform expression in mature adipocytes by RNA sequencing.

Results

We find 777 exon junctions with robust differential usage between white and beige adipocytes in all 6 subjects, mapping to 562 genes. Importantly, only 10% of these differentially spliced genes are also differentially expressed, indicating that alternative splicing constitutes an additional layer of gene expression regulation during beige adipocyte differentiation. Functional classification of alternative isoforms points to a gain of function for key thermogenic transcription factors such as PPARG and CITED1, and enzymes such as PEMT, or LPIN1. We find that a large majority of the splice variants arise from differential TSS usage, with beige-specific TSSs being enriched for PPARγ and MED1 binding compared to white-specific TSSs. Finally, we validate beige specific isoform expression at the protein level for two thermogenic regulators, PPARγ and PEMT.

Discussion

These results suggest that differential isoform expression through alternative TSS usage is an important regulatory mechanism for human adipocyte thermogenic specification.

Deep Transcriptomics Reveals Cell-Specific Isoforms of Pan-Neuronal Genes

Profiling gene expression in single neurons using single-cell RNA-Seq is a powerful method for understanding the molecular diversity of the nervous system. Profiling alternative splicing in single neurons using these methods is more challenging, however, due to low capture efficiency and sensitivity. As a result, we know much less about splicing patterns and regulation across neurons than we do about gene expression. Here we leverage unique attributes of the C. elegans nervous system to investigate deep cell-specific transcriptomes complete with biological replicates generated by the CeNGEN consortium, enabling high-confidence assessment of splicing across neuron types even for lowly-expressed genes. Global splicing maps reveal several striking observations, including pan-neuronal genes that harbor cell-specific splice variants, abundant differential intron retention across neuron types, and a single neuron highly enriched for upstream alternative 3' splice sites. We develop an algorithm to identify unique cell-specific expression patterns and use it to discover both cell-specific isoforms and potential regulatory RNA binding proteins that establish these isoforms. Genetic interrogation of these RNA binding proteins in vivo identifies three distinct regulatory factors employed to establish unique splicing patterns in a single neuron. Finally, we develop a user-friendly platform for spatial transcriptomic visualization of these splicing patterns with single-neuron resolution.

High expression of SRSF1 facilitates osteosarcoma progression and unveils its potential mechanisms

BACKGROUND

SRSF1, a member of Serine/Arginine-Rich Splicing Factors (SRSFs), has been observed to significantly influence cancer progression. However, the precise role of SRSF1 in osteosarcoma (OS) remains unclear. This study aims to investigate the functions of SRSF1 and its underlying mechanism in OS.

METHODS

SRSF1 expression level in OS was evaluated on the TCGA dataset, TAGET-OS database. qRT-PCR and Western blotting were employed to assess SRSF1 expression in human OS cell lines as well as the interfered ectopic expression states. The effect of SRSF1 on cell migration, invasion, proliferation, and apoptosis of OS cells were measured by transwell assay and flow cytometry. RNA sequence and bioinformatic analyses were conducted to elucidate the targeted genes, relevant biological pathways, and alternative splicing (AS) events regulated by SRSF1.

RESULTS

SRSF1 expression was consistently upregulated in both OS samples and OS cell lines. Diminishing SRSF1 resulted in reduced proliferation, migration, and invasion and increased apoptosis in OS cells while overexpressing SRSF1 led to enhanced growth, migration, invasion, and decreased apoptosis. Mechanistically, Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and Gene Set Enrichment Analysis (GSEA) revealed that the biological functions of SRSF1 were closely associated with the dysregulation of the protein targeting processes, location of the cytosolic ribosome, extracellular matrix (ECM), and proteinaceous extracellular matrix, along with the PI3K-AKT pathway, Wnt pathway, and HIPPO pathway. Transcriptome analysis identified AS events modulated by SRSF1, especially (Skipped Exon) SE events and (Mutually exclusive Exons) MXE events, revealing potential roles of targeted molecules in mRNA surveillance, RNA degradation, and RNA transport during OS development. qRT-PCR confirmed that SRSF1 knockdown resulted in the occurrence of alternative splicing of SRRM2, DMKN, and SCAT1 in OS.

CONCLUSIONS

Our results highlight the oncogenic role of high SRSF1 expression in promoting OS progression, and further explore the potential mechanisms of action. The significant involvement of SRSF1 in OS development suggests its potential utility as a therapeutic target in OS.