Mechanisms of alternative pre-messenger RNA splicing

Epigenetics and alternative splicing in cancer: old enemies, new perspectives

In recent years, significant strides in both conceptual understanding and technological capabilities have bolstered our comprehension of the factors underpinning cancer initiation and progression. While substantial insights have unraveled the molecular mechanisms driving carcinogenesis, there has been an overshadowing of the critical contribution made by epigenetic pathways, which works in concert with genetics. Mounting evidence demonstrates cancer as a complex interplay between genetics and epigenetics. Notably, epigenetic elements play a pivotal role in governing alternative pre-mRNA splicing, a primary contributor to protein diversity. In this review, we have provided detailed insights into the bidirectional communication between epigenetic modifiers and alternative splicing, providing examples of specific genes and isoforms affected. Notably, succinct discussion on targeting epigenetic regulators and the potential of the emerging field of epigenome editing to modulate splicing patterns is also presented. In summary, this review offers valuable insights into the intricate interplay between epigenetics and alternative splicing in cancer, paving the way for novel approaches to understanding and targeting this critical process.

Alternative splicing in ovarian cancer

Ovarian cancer is the second leading cause of gynecologic cancer death worldwide, with only 20% of cases detected early due to its elusive nature, limiting successful treatment. Most deaths occur from the disease progressing to advanced stages. Despite advances in chemo- and immunotherapy, the 5-year survival remains below 50% due to high recurrence and chemoresistance. Therefore, leveraging new research perspectives to understand molecular signatures and identify novel therapeutic targets is crucial for improving the clinical outcomes of ovarian cancer. Alternative splicing, a fundamental mechanism of post-transcriptional gene regulation, significantly contributes to heightened genomic complexity and protein diversity. Increased awareness has emerged about the multifaceted roles of alternative splicing in ovarian cancer, including cell proliferation, metastasis, apoptosis, immune evasion, and chemoresistance. We begin with an overview of altered splicing machinery, highlighting increased expression of spliceosome components and associated splicing factors like BUD31, SF3B4, and CTNNBL1, and their relationships to ovarian cancer. Next, we summarize the impact of specific variants of CD44, ECM1, and KAI1 on tumorigenesis and drug resistance through diverse mechanisms. Recent genomic and bioinformatics advances have enhanced our understanding. By incorporating data from The Cancer Genome Atlas RNA-seq, along with clinical information, a series of prognostic models have been developed, which provided deeper insights into how the splicing influences prognosis, overall survival, the immune microenvironment, and drug sensitivity and resistance in ovarian cancer patients. Notably, novel splicing events, such as PIGV|1299|AP and FLT3LG|50,941|AP, have been identified in multiple prognostic models and are associated with poorer and improved prognosis, respectively. These novel splicing variants warrant further functional characterization to unlock the underlying molecular mechanisms. Additionally, experimental evidence has underscored the potential therapeutic utility of targeting alternative splicing events, exemplified by the observation that knockdown of splicing factor BUD31 or antisense oligonucleotide-induced BCL2L12 exon skipping promotes apoptosis of ovarian cancer cells. In clinical settings, bevacizumab, a humanized monoclonal antibody that specifically targets the VEGF-A isoform, has demonstrated beneficial effects in the treatment of patients with advanced epithelial ovarian cancer. In conclusion, this review constitutes the first comprehensive and detailed exposition of the intricate interplay between alternative splicing and ovarian cancer, underscoring the significance of alternative splicing events as pivotal determinants in cancer biology and as promising avenues for future diagnostic and therapeutic intervention.

On the Roles of Protein Intrinsic Disorder in the Origin of Life and Evolution

Obviously, the discussion of different factors that could have contributed to the origin of life and evolution is clear speculation, since there is no way of checking the validity of most of the related hypotheses in practice, as the corresponding events not only already happened, but took place in a very distant past. However, there are a few undisputable facts that are present at the moment, such as the existence of a wide variety of living forms and the abundant presence of intrinsically disordered proteins (IDPs) or hybrid proteins containing ordered domains and intrinsically disordered regions (IDRs) in all living forms. Since it seems that the currently existing living forms originated from a common ancestor, their variety is a result of evolution. Therefore, one could ask a logical question of what role(s) the structureless and highly dynamic but vastly abundant and multifunctional IDPs/IDRs might have in evolution. This study represents an attempt to consider various ideas pertaining to the potential roles of protein intrinsic disorder in the origin of life and evolution.

Nanopore-based full-length transcriptome sequencing of the skin in Pseudopleuronectes yokohamae identifies novel antimicrobial peptide genes

The marbled flounder (Pseudopleuronectes yokohamae) is highly esteemed for its exceptional nutritional value and delicious taste. However, this species has extremely limited transcriptome data, which can offer priceless information for disease protection. In the study, the skin transcriptomic sequencing of P. yokohamae revealed 7.72 GB of clean data using the Nanopore sequencing platform. The results revealed 30,498 transcripts of functional annotations in the P. yokohamae transcriptome. All transcripts were searched in eight functional databases. A total of 10,337 ORFs were obtained, of which 6081 complete ORFs accounted for 58.83% of all predicted CDS. Moreover, 10,195 SSRs were detected. Meanwhile, the non-pecific immunity pathways were investigated for better understanding of the immunological reaction in P. yokohamae, and seven innate immune pathways were identified. The innate-immune related genes were highly expressed in the NOD-like receptor signaling pathway, followed by the C-type lectin receptor signaling pathway, Toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway and Cytosolic DNA-sensing pathway. In this study, four families of antimicrobial peptides (AMPs) in P. yokohamae were analysed for the first time, including piscidins, hepcidins, liver-expressed antimicrobial peptide and defensins. Seven AMPs, including Pypleurocidin-like WF3, Pypleurocidin-like WFX, Pyhepcidin 1, Pyhepcidin-like 1, PyLEAP-2, Pybeta-defensin and Pybeta-defensin-like 1, were further identified. The seven AMPs showed a highly identity in their cDNA and genomic structures and an inducible expression pattern preferable to skin in response to pathogens. The transcriptomic data and investigation of AMPs from P. yokohamae promote a deeper awareness of fish mucosal immunity and provide information in the prevention of fish diseases.

From computational models of the splicing code to regulatory mechanisms and therapeutic implications

Since the discovery of RNA splicing and its role in gene expression, researchers have sought a set of rules, an algorithm or a computational model that could predict the splice isoforms, and their frequencies, produced from any transcribed gene in a specific cellular context. Over the past 30 years, these models have evolved from simple position weight matrices to deep-learning models capable of integrating sequence data across vast genomic distances. Most recently, new model architectures are moving the field closer to context-specific alternative splicing predictions, and advances in sequencing technologies are expanding the type of data that can be used to inform and interpret such models. Together, these developments are driving improved understanding of splicing regulatory mechanisms and emerging applications of the splicing code to the rational design of RNA- and splicing-based therapeutics.

LsRTDv1, a reference transcript dataset for accurate transcript-specific expression analysis in lettuce

Accurate quantification of gene and transcript-specific expression, with the underlying knowledge of precise transcript isoforms, is crucial to understanding many biological processes. Analysis of RNA sequencing data has benefited from the development of alignment-free algorithms which enhance the precision and speed of expression analysis. However, such algorithms require a reference transcriptome. Here we generate a reference transcript dataset (LsRTDv1) for lettuce (cv. Saladin), combining long- and short-read sequencing with publicly available transcriptome annotations, and filtering to keep only transcripts with high-confidence splice junctions and transcriptional start and end sites. LsRTDv1 identifies novel genes (mostly long non-coding RNAs) and increases the number of transcript isoforms per gene in the lettuce genome from 1.4 to 2.7. We show that LsRTDv1 significantly increases the mapping rate of RNA-seq data from a lettuce time-series experiment (mock- and Botrytis cinerea-inoculated) and enables detection of genes that are differentially alternatively spliced in response to infection as well as transcript-specific expression changes. LsRTDv1 is a valuable resource for investigation of transcriptional and alternative splicing regulation in lettuce.

Cell-type-specific splicing of transcription regulators and Ptbp1 by Rbfox1/2/3 in the developing neocortex

How master splicing regulators crosstalk with each other and to what extent transcription regulators are differentially spliced remain unclear in the developing brain. Here, cell-type-specific RNA-Seq of the developing neocortex uncover that transcription regulators are enriched for differential splicing, altering protein isoforms or inducing nonsense-mediated mRNA decay. Transient expression of Rbfox proteins in radial glia progenitors induces neuronal splicing events preferentially in transcription regulators such as Meis2 and Tead1. Surprisingly, Rbfox proteins promote the inclusion of a mammal-specific alternative exon and a previously undescribed poison exon in Ptbp1. Simultaneous ablation of Rbfox1/2/3 in the neocortex downregulates neuronal isoforms and disrupts radial neuronal migration. Furthermore, the progenitor isoform of Meis2 promotes Tgfb3 transcription, while the Meis2 neuron isoform promotes neuronal differentiation. These observations indicate that transcription regulators are differentially spliced between cell types in the developing neocortex.

Reference-informed prediction of alternative splicing and splicing-altering mutations from sequences

Alternative splicing plays a crucial role in protein diversity and gene expression regulation in higher eukaryotes, and mutations causing dysregulated splicing underlie a range of genetic diseases. Computational prediction of alternative splicing from genomic sequences not only provides insight into gene-regulatory mechanisms but also helps identify disease-causing mutations and drug targets. However, the current methods for the quantitative prediction of splice site usage still have limited accuracy. Here, we present DeltaSplice, a deep neural network model optimized to learn the impact of mutations on quantitative changes in alternative splicing from the comparative analysis of homologous genes. The model architecture enables DeltaSplice to perform "reference-informed prediction" by incorporating the known splice site usage of a reference gene sequence to improve its prediction on splicing-altering mutations. We benchmarked DeltaSplice and several other state-of-the-art methods on various prediction tasks, including evolutionary sequence divergence on lineage-specific splicing and splicing-altering mutations in human populations and neurodevelopmental disorders, and demonstrated that DeltaSplice outperformed consistently. DeltaSplice predicted ∼15% of splicing quantitative trait loci (sQTLs) in the human brain as causal splicing-altering variants. It also predicted splicing-altering de novo mutations outside the splice sites in a subset of patients affected by autism and other neurodevelopmental disorders (NDDs), including 19 genes with recurrent splicing-altering mutations. Integration of splicing-altering mutations with other types of de novo mutation burdens allowed the prediction of eight novel NDD-risk genes. Our work expanded the capacity of in silico splicing models with potential applications in genetic diagnosis and the development of splicing-based precision medicine.

Interdependent regulation of alternative splicing by SR and hnRNP proteins

Alternative pre-mRNA splicing is a combinatorial process involving SR and hnRNP splicing factors. These proteins can silence or enhance splicing based on their expression levels and binding positions. To better understand their combinatorial and interdependent regulation, computational analyses were performed using HepG2 and K562 cell knockdown and binding datasets from the ENCODE Project. Analyses of diMerential splicing for 6 SR proteins and 13 hnRNP knockdowns revealed statistically significant exon overlap among most RBP combinations, albeit at diMerent levels. Neither SR proteins nor hnRNPs showed strong preferences for collaborating with specific RBP classes in mediating exon inclusion. While SRSF1, hnRNPK, and hnRNPC stand out as major influencers of alternative splicing, they do so predominantly independent of other RBPs. Meanwhile, minor influencers of alternative splicing such as hnRNPAB and hnRNPA0 predominantly regulate exon inclusion in concert with other RBPs, indicating that inclusion can be mediated by both single and multiple RBPs. Interestingly, the higher the number of RBPs that regulate the inclusion of an exon, the more variable exon inclusion preferences become. Interdependently regulated exons are more modular and have diMerent physical characteristics such as reduced exon length compared to their independent counterparts. A comparison of RBP interdependence between HepG2 and K562 cells provides the framework that explains cell-type-specific alternative splicing. Our study highlights the importance of the interdependent regulation of alternative exons and identifies characteristics of interdependently regulated exons that diMer from independently regulated exons.

A defective splicing machinery promotes senescence through MDM4 alternative splicing

Defects in the splicing machinery are implicated in various diseases, including cancer. We observed a general reduction in the expression of spliceosome components and splicing regulators in human cell lines undergoing replicative, stress-induced, and telomere uncapping-induced senescence. Supporting the view that defective splicing contributes to senescence, splicing inhibitors herboxidiene, and pladienolide B induced senescence in normal and cancer cell lines. Furthermore, depleting individual spliceosome components also promoted senescence. All senescence types were associated with an alternative splicing transition from the MDM4-FL variant to MDM4-S. The MDM4 splicing shift was reproduced when splicing was inhibited, and spliceosome components were depleted. While decreasing the level of endogenous MDM4 promoted senescence and cell survival independently of the MDM4-S expression status, cell survival was also improved by increasing MDM4-S. Overall, our work establishes that splicing defects modulate the alternative splicing of MDM4 to promote senescence and cell survival.

Two C-terminal isoforms of Aplysia tachykinin-related peptide receptors exhibit phosphorylation-dependent and phosphorylation-independent desensitization mechanisms

Diversity, a hallmark of G protein-coupled receptor (GPCR) signaling, partly stems from alternative splicing of a single gene generating more than one isoform for a receptor. Additionally, receptor responses to ligands can be attenuated by desensitization upon prolonged or repeated ligand exposure. Both phenomena have been demonstrated and exemplified by the deuterostome tachykinin signaling system, although the role of phosphorylation in desensitization remains a subject of debate. Here, we describe the signaling system for tachykinin-related peptides (TKRPs) in a protostome, mollusk Aplysia. We cloned the Aplysia TKRP precursor, which encodes three TKRPs (apTKRP-1, apTKRP-2a, and apTKRP-2b) containing the FXGXR-amide motif. In situ hybridization and immunohistochemistry showed predominant expression of TKRP mRNA and peptide in the cerebral ganglia. TKRPs and their posttranslational modifications were observed in extracts of central nervous system ganglia using mass spectrometry. We identified two Aplysia TKRP receptors (apTKRPRs), named apTKRPR-A and apTKRPR-B. These receptors are two isoforms generated through alternative splicing of the same gene and differ only in their intracellular C termini. Structure-activity relationship analysis of apTKRP-2b revealed that both C-terminal amidation and conserved residues of the ligand are critical for receptor activation. C-terminal truncates and mutants of apTKRPRs suggested that there is a C-terminal phosphorylation-independent desensitization for both receptors. Moreover, apTKRPR-B also exhibits phosphorylation-dependent desensitization through the phosphorylation of C-terminal Ser/Thr residues. This comprehensive characterization of the Aplysia TKRP signaling system underscores the evolutionary conservation of the TKRP and TK signaling systems, while highlighting the intricacies of receptor regulation through alternative splicing and differential desensitization mechanisms.

Decoding RNA Metabolism by RNA-linked CRISPR Screening in Human Cells

RNAs undergo a complex choreography of metabolic processes in human cells that are regulated by thousands of RNA-associated proteins. While the effects of individual RNA-associated proteins on RNA metabolism have been extensively characterized, the full complement of regulators for most RNA metabolic events remain unknown. Here we present a massively parallel RNA-linked CRISPR (ReLiC) screening approach to measure the responses of diverse RNA metabolic events to knockout of 2,092 human genes encoding all known RNA-associated proteins. ReLiC screens highlight modular interactions between gene networks regulating splicing, translation, and decay of mRNAs. When combined with biochemical fractionation of polysomes, ReLiC reveals striking pathway-specific coupling between growth fitness and mRNA translation. Perturbing different components of the translation and proteostasis machineries have distinct effects on ribosome occupancy, while perturbing mRNA transcription leaves ribosome occupancy largely intact. Isoform-selective ReLiC screens capture differential regulation of intron retention and exon skipping by SF3b complex subunits. Chemogenomic screens using ReLiC decipher translational regulators upstream of mRNA decay and uncover a role for the ribosome collision sensor GCN1 during treatment with the anti-leukemic drug homoharringtonine. Our work demonstrates ReLiC as a versatile platform for discovering and dissecting regulatory principles of human RNA metabolism.

Alternative splicing regulation and its therapeutic potential in bladder cancer

Bladder cancer is one of the leading causes of mortality globally. The development of bladder cancer is closely associated with alternative splicing, which regulates human gene expression and enhances the diversity of functional proteins. Alternative splicing is a distinctive feature of bladder cancer, and as such, it may hold promise as a therapeutic target. This review aims to comprehensively discuss the current knowledge of alternative splicing in the context of bladder cancer. We review the process of alternative splicing and its regulation in bladder cancer. Moreover, we emphasize the significance of abnormal alternative splicing and splicing factor irregularities during bladder cancer progression. Finally, we explore the impact of alternative splicing on bladder cancer drug resistance and the potential of alternative splicing as a therapeutic target.

Splice-switching antisense oligonucleotides for pediatric neurological disorders

Pediatric neurological disorders are frequently devastating and present unmet needs for effective medicine. The successful treatment of spinal muscular atrophy with splice-switching antisense oligonucleotides (SSO) indicates a feasible path to targeting neurological disorders by redirecting pre-mRNA splicing. One direct outcome is the development of SSOs to treat haploinsufficient disorders by targeting naturally occurring non-productive splice isoforms. The development of personalized SSO treatment further inspired the therapeutic exploration of rare diseases. This review will discuss the recent advances that utilize SSOs to treat pediatric neurological disorders.

The interferon-regulated host factor hnRNPA0 modulates HIV-1 production by interference with LTR activity, mRNA trafficking, and programmed ribosomal frameshifting

The interplay between host factors and viral components impacts viral replication efficiency profoundly. Members of the cellular heterogeneous nuclear ribonucleoprotein family (hnRNPs) have been extensively studied as HIV-1 host dependency factors, but whether they play a role in innate immunity is currently unknown. This study aimed to identify hnRNPA0 as a type I interferon (IFN)-repressed host factor in HIV-1-infected cells. Knockdown of hnRNPA0, a situation that mirrors conditions under IFN stimulation, increased LTR activity, export of unspliced HIV-1 mRNA, viral particle production, and thus, increased infectivity. Conversely, hnRNPA0 overexpression primarily reduced plasmid-driven and integrated HIV-1 long terminal repeat (LTR) activity, significantly decreasing total viral mRNA and protein levels. In addition, high levels of hnRNPA0 significantly reduced the HIV-1 programmed ribosomal frameshifting efficiency, resulting in a shift in the HIV-1 p55/p15 ratio. The HIV-1 alternative splice site usage remained largely unaffected by altered hnRNPA0 levels suggesting that the synergistic inhibition of the LTR activity and viral mRNA transcription, as well as impaired ribosomal frameshifting efficiency, are critical factors for efficient HIV-1 replication regulated by hnRNPA0. The pleiotropic dose-dependent effects under high or low hnRNPA0 levels were further confirmed in HIV-1-infected Jurkat cells. Finally, our study revealed that hnRNPA0 levels in PBMCs were lower in therapy-naive HIV-1-infected individuals compared to healthy controls. Our findings highlight a significant role for hnRNPA0 in HIV-1 replication and suggest that its IFN-I-regulated expression levels are critical for viral fitness allowing replication in an antiviral environment.IMPORTANCERNA-binding proteins, in particular, heterogeneous nuclear ribonucleoproteins (hnRNPs), have been extensively studied. Some act as host dependency factors for HIV-1 since they are involved in multiple cellular gene expression processes. Our study revealed hnRNPA0 as an IFN-regulated host factor, that is differently expressed after IFN-I treatment in HIV-1 target cells and lower expressed in therapy-naïve HIV-1-infected individuals. Our findings demonstrate the significant pleiotropic role of hnRNPA0 in viral replication: In high concentrations, hnRNPA0 limits viral replication by negatively regulating Tat-LTR transcription, retaining unspliced mRNA in the nucleus, and significantly impairing programmed ribosomal frameshifting. Low hnRNPA0 levels as observed in IFN-treated THP-1 cells, particularly facilitate HIV LTR activity and unspliced mRNA export, suggesting a role in innate immunity in favor of HIV replication. Understanding the mode of action between hnRNPA0 and HIV-1 gene expression might help to identify novel therapeutically strategies against HIV-1 and other viruses.

The role of HnrnpF/H as a driver of oligoteratozoospermia

Male subfertility or infertility is a common condition often characterized by men producing a low number of sperm with poor quality. To gain insight into this condition, we performed a quantitative proteomic analysis of semen samples obtained from infertile and fertile men. At least 6 proteins showed significant differences in regulation of alternatively spliced isoforms. To investigate this link between aberrant alternative splicing and production of poor-quality spermatozoa, we overexpressed the hnrnpH/F-orthologue Glorund (Glo) in Drosophila, which was also found to be abundant in poor quality human sperm. Transgenic animals produced low numbers of morphologically defective spermatozoa and aberrant formation of the "dense body," an organelle akin to the mammalian manchette. Furthermore, fertility trials demonstrated that transgenic flies were either completely infertile or highly subfertile. These findings suggest that dysregulation of hnrnpH/F is likely to result in the production of low-quality semen, leading to subfertility or infertility in men.

The Rbfox1/LASR complex controls alternative pre-mRNA splicing by recognition of multi-part RNA regulatory modules

The Rbfox proteins regulate alternative pre-mRNA splicing by binding to the RNA element GCAUG. In the nucleus, most of Rbfox is bound to LASR, a complex of RNA-binding proteins that recognize additional RNA motifs. However, it remains unclear how the different subunits of the Rbfox/LASR complex act together to bind RNA and regulate splicing. We used a nuclease-protection assay to map the transcriptome-wide footprints of Rbfox1/LASR on nascent cellular RNA. In addition to GCAUG, Rbfox1/LASR binds RNA containing motifs for LASR subunits hnRNPs M, H/F, C, and Matrin3. These elements are often arranged in tandem, forming multi-part modules of RNA motifs. To distinguish contact sites of Rbfox1 from the LASR subunits, we analyzed a mutant Rbfox1(F125A) that has lost RNA binding but remains associated with LASR. Rbfox1(F125A)/LASR complexes no longer interact with GCAUG but retain binding to RNA elements for LASR. Splicing analyses reveal that in addition to activating exons through adjacent GCAUG elements, Rbfox can also stimulate exons near binding sites for LASR subunits. Mini-gene experiments demonstrate that these diverse elements produce a combined regulatory effect on a target exon. These findings illuminate how a complex of RNA-binding proteins can decode combinatorial splicing regulatory signals by recognizing groups of tandem RNA elements.

Spliceosomal complex components are critical for adjusting the C:N balance during high-light acclimation

Plant acclimation to an ever-changing environment is decisive for growth, reproduction, and survival. Light availability limits biomass production on both ends of the intensity spectrum. Therefore, the adjustment of plant metabolism is central to high-light (HL) acclimation, and the accumulation of photoprotective anthocyanins is commonly observed. However, mechanisms and factors regulating the HL acclimation response are less clear. Two Arabidopsis mutants of spliceosome components exhibiting a pronounced anthocyanin overaccumulation in HL were isolated from a forward genetic screen for new factors crucial for plant acclimation. Time-resolved physiological, transcriptome, and metabolome analysis revealed a vital function of the spliceosome components for rapidly adjusting gene expression and metabolism. Deficiency of INCREASED LEVEL OF POLYPLOIDY1 (ILP1), NTC-RELATED PROTEIN1 (NTR1), and PLEIOTROPIC REGULATORY LOCUS1 (PRL1) resulted in a marked overaccumulation of carbohydrates and strongly diminished amino acid biosynthesis in HL. While not generally limited in N-assimilation, ilp1, ntr1, and prl1 showed higher glutamate levels and reduced amino acid biosynthesis in HL. The comprehensive analysis reveals a function of the spliceosome components in the conditional regulation of the carbon:nitrogen balance and the accumulation of anthocyanins during HL acclimation. The importance of gene expression, metabolic regulation, and re-direction of carbon towards anthocyanin biosynthesis for HL acclimation are discussed.

Inteins: A Swiss army knife for synthetic biology

Inteins are proteins found in nature that execute protein splicing. Among them, split inteins stand out for their versatility and adaptability, presenting creative solutions for addressing intricate challenges in various biological applications. Their exquisite attributes, including compactness, reliability, orthogonality, low toxicity, and irreversibility, make them of interest to various fields including synthetic biology, biotechnology and biomedicine. In this review, we delve into the inherent challenges of using inteins, present approaches for overcoming these challenges, and detail their reliable use for specific cellular tasks. We will discuss the use of conditional inteins in areas like cancer therapy, drug screening, patterning, infection treatment, diagnostics and biocontainment. Additionally, we will underscore the potential of inteins in executing basic logical operations with practical implications. We conclude by showcasing their potential in crafting complex genetic circuits for performing computations and feedback control that achieves robust perfect adaptation.

Disruption of mitochondrial pyruvate oxidation in dorsal root ganglia drives persistent nociceptive sensitization and causes pervasive transcriptomic alterations

ABSTRACT

Metabolism is inextricably linked to every aspect of cellular function. In addition to energy production and biosynthesis, metabolism plays a crucial role in regulating signal transduction and gene expression. Altered metabolic states have been shown to maintain aberrant signaling and transcription, contributing to diseases like cancer, cardiovascular disease, and neurodegeneration. Metabolic gene polymorphisms and defects are also associated with chronic pain conditions, as are increased levels of nerve growth factor (NGF). However, the mechanisms by which NGF may modulate sensory neuron metabolism remain unclear. This study demonstrated that intraplantar NGF injection reprograms sensory neuron metabolism. Nerve growth factor suppressed mitochondrial pyruvate oxidation and enhanced lactate extrusion, requiring 24 hours to increase lactate dehydrogenase A and pyruvate dehydrogenase kinase 1 (PDHK1) expression. Inhibiting these metabolic enzymes reversed NGF-mediated effects. Remarkably, directly disrupting mitochondrial pyruvate oxidation induced severe, persistent allodynia, implicating this metabolic dysfunction in chronic pain. Nanopore long-read sequencing of poly(A) mRNA uncovered extensive transcriptomic changes upon metabolic disruption, including altered gene expression, splicing, and poly(A) tail lengths. By linking metabolic disturbance of dorsal root ganglia to transcriptome reprogramming, this study enhances our understanding of the mechanisms underlying persistent nociceptive sensitization. These findings imply that impaired mitochondrial pyruvate oxidation may drive chronic pain, possibly by impacting transcriptomic regulation. Exploring these metabolite-driven mechanisms further might reveal novel therapeutic targets for intractable pain.

Germline mutations of breast cancer susceptibility genes through expanded genetic analysis in unselected Colombian patients

BACKGROUND

In Colombia and worldwide, breast cancer (BC) is the most frequently diagnosed neoplasia and the leading cause of death from cancer among women. Studies predominantly involve hereditary and familial cases, demonstrating a gap in the literature regarding the identification of germline mutations in unselected patients from Latin-America. Identification of pathogenic/likely pathogenic (P/LP) variants is important for shaping national genetic analysis policies, genetic counseling, and early detection strategies. The present study included 400 women with unselected breast cancer (BC), in whom we analyzed ten genes, using Whole Exome Sequencing (WES), know to confer risk for BC, with the aim of determining the genomic profile of previously unreported P/LP variants in the affected population. Additionally, Multiplex Ligation-dependent Probe Amplification (MLPA) was performed to identify Large Genomic Rearrangements (LGRs) in the BRCA1/2 genes. To ascertain the functional impact of a recurrent intronic variant (ATM c.5496 + 2_5496 + 5delTAAG), a minigene assay was conducted.

RESULTS

We ascertained the frequency of P/LP germline variants in BRCA2 (2.5%), ATM (1.25%), BRCA1 (0.75%), PALB2 (0.50%), CHEK2 (0.50%), BARD1 (0.25%), and RAD51D (0.25%) genes in the population of study. P/LP variants account for 6% of the total population analyzed. No LGRs were detected in our study. We identified 1.75% of recurrent variants in BRCA2 and ATM genes. One of them corresponds to the ATM c.5496 + 2_5496 + 5delTAAG. Functional validation of this variant demonstrated a splicing alteration probably modifying the Pincer domain and subsequent protein structure.

CONCLUSION

This study described for the first time the genomic profile of ten risk genes in Colombian women with unselected BC. Our findings underscore the significance of population-based research, advocating the consideration of molecular testing in all women with cancer.

Interferon-Regulated Expression of Cellular Splicing Factors Modulates Multiple Levels of HIV-1 Gene Expression and Replication

Type I interferons (IFN-Is) are pivotal in innate immunity against human immunodeficiency virus I (HIV-1) by eliciting the expression of IFN-stimulated genes (ISGs), which encompass potent host restriction factors. While ISGs restrict the viral replication within the host cell by targeting various stages of the viral life cycle, the lesser-known IFN-repressed genes (IRepGs), including RNA-binding proteins (RBPs), affect the viral replication by altering the expression of the host dependency factors that are essential for efficient HIV-1 gene expression. Both the host restriction and dependency factors determine the viral replication efficiency; however, the understanding of the IRepGs implicated in HIV-1 infection remains greatly limited at present. This review provides a comprehensive overview of the current understanding regarding the impact of the RNA-binding protein families, specifically the two families of splicing-associated proteins SRSF and hnRNP, on HIV-1 gene expression and viral replication. Since the recent findings show specifically that SRSF1 and hnRNP A0 are regulated by IFN-I in various cell lines and primary cells, including intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs), we particularly discuss their role in the context of the innate immunity affecting HIV-1 replication.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

LncRNA-mediated orchestrations of alternative splicing in the landscape of breast cancer

Alternative splicing (AS) is a fundamental post-transcriptional process in eukaryotes, enabling a single gene to generate diverse mRNA transcripts, thereby enhancing protein variability. This process involves the excision of introns and the joining of exons in pre-mRNA(s) to form mature mRNA. The resulting mature mRNAs exhibit various combinations of exons, contributing to functional diversity. Dysregulation of AS can substantially modulate protein functions, impacting the onset and progression of numerous diseases, including cancer. Non-coding RNAs (ncRNAs) are distinct from protein-coding RNAs and consist of short and long types. Long non-coding RNAs (lncRNAs) play an important role in regulating several cellular processes, particularly alternative splicing, according to new research. This review provides insight into the latest discoveries concerning how lncRNAs influence alternative splicing within the realm of breast cancer. Additionally, it explores potential therapeutic strategies focused on targeting lncRNAs.

Structural insights into the cross-exon to cross-intron spliceosome switch

Early spliceosome assembly can occur through an intron-defined pathway, whereby U1 and U2 small nuclear ribonucleoprotein particles (snRNPs) assemble across the intron1. Alternatively, it can occur through an exon-defined pathway2-5, whereby U2 binds the branch site located upstream of the defined exon and U1 snRNP interacts with the 5' splice site located directly downstream of it. The U4/U6.U5 tri-snRNP subsequently binds to produce a cross-intron (CI) or cross-exon (CE) pre-B complex, which is then converted to the spliceosomal B complex6,7. Exon definition promotes the splicing of upstream introns2,8,9 and plays a key part in alternative splicing regulation10-16. However, the three-dimensional structure of exon-defined spliceosomal complexes and the molecular mechanism of the conversion from a CE-organized to a CI-organized spliceosome, a pre-requisite for splicing catalysis, remain poorly understood. Here cryo-electron microscopy analyses of human CE pre-B complex and B-like complexes reveal extensive structural similarities with their CI counterparts. The results indicate that the CE and CI spliceosome assembly pathways converge already at the pre-B stage. Add-back experiments using purified CE pre-B complexes, coupled with cryo-electron microscopy, elucidate the order of the extensive remodelling events that accompany the formation of B complexes and B-like complexes. The molecular triggers and roles of B-specific proteins in these rearrangements are also identified. We show that CE pre-B complexes can productively bind in trans to a U1 snRNP-bound 5' splice site. Together, our studies provide new mechanistic insights into the CE to CI switch during spliceosome assembly and its effect on pre-mRNA splice site pairing at this stage.

Role of alternative splicing in fish immunity

Alternative splicing serves as a pivotal source of complexity in the transcriptome and proteome, selectively connecting various coding elements to generate a diverse array of mRNAs. This process encodes multiple proteins with either similar or distinct functions, contributing significantly to the intricacies of cellular processes. The role of alternative splicing in mammalian immunity has been well studied. Remarkably, the immune system of fish shares substantial similarities with that of humans, and alternative splicing also emerges as a key player in the immune processes of fish. In this review, we offer an overview of alternative splicing and its associated functions in the immune processes of fish, and summarize the research progress on alternative splicing in the fish immunity. Furthermore, we review the impact of alternative splicing on the fish immune system's response to external stimuli. Finally, we present our perspectives on future directions in this field. Our aim is to provide valuable insights for the future investigations into the role of alternative splicing in immunity.

Concordant Gene Expression and Alternative Splicing Regulation under Abiotic Stresses in Arabidopsis

The current investigation endeavors to identify differentially expressed alternatively spliced (DAS) genes that exhibit concordant expression with splicing factors (SFs) under diverse multifactorial abiotic stress combinations in Arabidopsis seedlings. SFs serve as the post-transcriptional mechanism governing the spatiotemporal dynamics of gene expression. The different stresses encompass variations in salt concentration, heat, intensive light, and their combinations. Clusters demonstrating consistent expression profiles were surveyed to pinpoint DAS/SF gene pairs exhibiting concordant expression. Through rigorous selection criteria, which incorporate alignment with documented gene functionalities and expression patterns observed in this study, four members of the serine/arginine-rich (SR) gene family were delineated as SFs concordantly expressed with six DAS genes. These regulated SF genes encompass cactin, SR1-like, SR30, and SC35-like. The identified concordantly expressed DAS genes encode diverse proteins such as the 26.5 kDa heat shock protein, chaperone protein DnaJ, potassium channel GORK, calcium-binding EF hand family protein, DEAD-box RNA helicase, and 1-aminocyclopropane-1-carboxylate synthase 6. Among the concordantly expressed DAS/SF gene pairs, SR30/DEAD-box RNA helicase, and SC35-like/1-aminocyclopropane-1-carboxylate synthase 6 emerge as promising candidates, necessitating further examinations to ascertain whether these SFs orchestrate splicing of the respective DAS genes. This study contributes to a deeper comprehension of the varied responses of the splicing machinery to abiotic stresses. Leveraging these DAS/SF associations shows promise for elucidating avenues for augmenting breeding programs aimed at fortifying cultivated plants against heat and intensive light stresses.

Knotify_V2.0: Deciphering RNA Secondary Structures with H-Type Pseudoknots and Hairpin Loops

Accurately predicting the pairing order of bases in RNA molecules is essential for anticipating RNA secondary structures. Consequently, this task holds significant importance in unveiling previously unknown biological processes. The urgent need to comprehend RNA structures has been accentuated by the unprecedented impact of the widespread COVID-19 pandemic. This paper presents a framework, Knotify_V2.0, which makes use of syntactic pattern recognition techniques in order to predict RNA structures, with a specific emphasis on tackling the demanding task of predicting H-type pseudoknots that encompass bulges and hairpins. By leveraging the expressive capabilities of a Context-Free Grammar (CFG), the suggested framework integrates the inherent benefits of CFG and makes use of minimum free energy and maximum base pairing criteria. This integration enables the effective management of this inherently ambiguous task. The main contribution of Knotify_V2.0 compared to earlier versions lies in its capacity to identify additional motifs like bulges and hairpins within the internal loops of the pseudoknot. Notably, the proposed methodology, Knotify_V2.0, demonstrates superior accuracy in predicting core stems compared to state-of-the-art frameworks. Knotify_V2.0 exhibited exceptional performance by accurately identifying both core base pairing that form the ground truth pseudoknot in 70% of the examined sequences. Furthermore, Knotify_V2.0 narrowed the performance gap with Knotty, which had demonstrated better performance than Knotify and even surpassed it in Recall and F1-score metrics. Knotify_V2.0 achieved a higher count of true positives (tp) and a significantly lower count of false negatives (fn) compared to Knotify, highlighting improvements in Prediction and Recall metrics, respectively. Consequently, Knotify_V2.0 achieved a higher F1-score than any other platform. The source code and comprehensive implementation details of Knotify_V2.0 are publicly available on GitHub.

Unearthing a novel function of SRSF1 in binding and unfolding of RNA G-quadruplexes

SRSF1 governs splicing of over 1500 mRNA transcripts. SRSF1 contains two RNA-recognition motifs (RRMs) and a C-terminal Arg/Ser-rich region (RS). It has been thought that SRSF1 RRMs exclusively recognize single-stranded exonic splicing enhancers, while RS lacks RNA-binding specificity. With our success in solving the insolubility problem of SRSF1, we can explore the unknown RNA-binding landscape of SRSF1. We find that SRSF1 RS prefers purine over pyrimidine. Moreover, SRSF1 binds to the G-quadruplex (GQ) from the ARPC2 mRNA, with both RRMs and RS being crucial. Our binding assays show that the traditional RNA-binding sites on the RRM tandem and the Arg in RS are responsible for GQ binding. Interestingly, our FRET and circular dichroism data reveal that SRSF1 unfolds the ARPC2 GQ, with RS leading unfolding and RRMs aiding. Our saturation transfer difference NMR results discover that Arg residues in SRSF1 RS interact with the guanine base but not other nucleobases, underscoring the uniqueness of the Arg/guanine interaction. Our luciferase assays confirm that SRSF1 can alleviate the inhibitory effect of GQ on gene expression in the cell. Given the prevalence of RNA GQ and SR proteins, our findings unveil unexplored SR protein functions with broad implications in RNA splicing and translation.

The diversity of splicing modifiers acting on A-1 bulged 5'-splice sites reveals rules for rational drug design

Pharmacological modulation of RNA splicing by small molecules is an emerging facet of drug discovery. In this context, the SMN2 splicing modifier SMN-C5 was used as a prototype to understand the mode of action of small molecule splicing modifiers and propose the concept of 5'-splice site bulge repair. In this study, we combined in vitro binding assays and structure determination by NMR spectroscopy to identify the binding modes of four other small molecule splicing modifiers that switch the splicing of either the SMN2 or the HTT gene. Here, we determined the solution structures of risdiplam, branaplam, SMN-CX and SMN-CY bound to the intermolecular RNA helix epitope containing an unpaired adenine within the G-2A-1G+1U+2 motif of the 5'-splice site. Despite notable differences in their scaffolds, risdiplam, SMN-CX, SMN-CY and branaplam contact the RNA epitope similarly to SMN-C5, suggesting that the 5'-splice site bulge repair mechanism can be generalised. These findings not only deepen our understanding of the chemical diversity of splicing modifiers that target A-1 bulged 5'-splice sites, but also identify common pharmacophores required for modulating 5'-splice site selection with small molecules.

CRISPR-dCas13d-based deep screening of proximal and distal splicing-regulatory elements

Pre-mRNA splicing, a key process in gene expression, can be therapeutically modulated using various drug modalities, including antisense oligonucleotides (ASOs). However, determining promising targets is hampered by the challenge of systematically mapping splicing-regulatory elements (SREs) in their native sequence context. Here, we use the catalytically inactive CRISPR-RfxCas13d RNA-targeting system (dCas13d/gRNA) as a programmable platform to bind SREs and modulate splicing by competing against endogenous splicing factors. SpliceRUSH, a high-throughput screening method, was developed to map SREs in any gene of interest using a lentivirus gRNA library that tiles the genetic region, including distal intronic sequences. When applied to SMN2, a therapeutic target for spinal muscular atrophy, SpliceRUSH robustly identifies not only known SREs but also a previously unknown distal intronic SRE, which can be targeted to alter exon 7 splicing using either dCas13d/gRNA or ASOs. This technology enables a deeper understanding of splicing regulation with applications for RNA-based drug discovery.

Diversity of mutations in the dystrophin gene and details of muscular lesions in porcine dystrophinopathies

During meat inspections in pigs, dystrophinopathies are among the muscle lesions targeted for disposal. In this study, the authors examined the lesions and the distribution of dystrophin expression in 25 pigs with dystrophinopathy. In addition, complementary deoxyribonucleic acid (cDNA) sequencing and western blotting were performed in 6 of the 25 cases, all of which were characterized by degeneration, necrosis, and fat replacement of muscle fibers. Comparing the results of immunohistochemistry with anti-dystrophin antibodies that recognized at different sites in the protein, the authors noted that the loss of dystrophin expression was most pronounced in the C-terminus-recognizing antibody (19/25 cases). The authors detected 5 missense mutations and 3 types of shortened transcripts generated by the skipping of exons in the cDNA, which were associated with the pathogenesis. One missense mutation had been reported previously, whereas the remaining mutations identified had not been previously documented in pigs. In the cases with shortened transcripts, normal-sized transcripts were detected together with the defective transcripts, suggesting that these mutations were caused by splicing abnormalities. In addition, they were in-frame mutations, all of which have similar pathogeneses of Becker muscular dystrophy in humans. These cases were 6 months of age and exhibited macroscopic discoloration, fatty replacement, and muscle degeneration, suggesting that the effect of these mutations on skeletal muscle was significant.

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Deciphering the role of alternative splicing in neoplastic diseases for immune-oncological therapies

Alternative splicing (AS) is an important molecular biological mechanism regulated by complex mechanisms involving a plethora of cis and trans-acting elements. Furthermore, AS is tissue specific and altered in various pathologies, including infectious, inflammatory, and neoplastic diseases. Recently developed immuno-oncological therapies include monoclonal antibodies (mAbs) and chimeric antigen receptor (CAR) T cells targeting, among others, immune checkpoint (ICP) molecules. Despite therapeutic successes have been demonstrated, only a limited number of patients showed long-term benefit from these therapies with tumor entity-related differential response rates were observed. Interestingly, splice variants of common immunotherapeutic targets generated by AS are able to completely escape and/or reduce the efficacy of mAb- and/or CAR-based tumor immunotherapies. Therefore, the analyses of splicing patterns of targeted molecules in tumor specimens prior to therapy might help correct stratification, thereby increasing therapy success by antibody panel selection and antibody dosages. In addition, the expression of certain splicing factors has been linked with the patients' outcome, thereby highlighting their putative prognostic potential. Outstanding questions are addressed to translate the findings into clinical application. This review article provides an overview of the role of AS in (tumor) diseases, its molecular mechanisms, clinical relevance, and therapy response.

RAVER1 hinders lethal EMT and modulates miR/RISC activity by the control of alternative splicing

The RAVER1 protein serves as a co-factor in guiding the polypyrimidine tract-binding protein (PTBP)-dependent control of alternative splicing (AS). Whether RAVER1 solely acts in concert with PTBPs and how it affects cancer cell fate remained elusive. Here, we provide the first comprehensive investigation of RAVER1-controlled AS in cancer cell models. This reveals a pro-oncogenic role of RAVER1 in modulating tumor growth and epithelial-mesenchymal-transition (EMT). Splicing analyses and protein-association studies indicate that RAVER1 guides AS in association with other splicing regulators, including PTBPs and SRSFs. In cancer cells, one major function of RAVER1 is the stimulation of proliferation and restriction of apoptosis. This involves the modulation of AS events within the miR/RISC pathway. Disturbance of RAVER1 impairs miR/RISC activity resulting in severely deregulated gene expression, which promotes lethal TGFB-driven EMT. Among others, RAVER1-modulated splicing events affect the insertion of protein interaction modules in factors guiding miR/RISC-dependent gene silencing. Most prominently, in all three human TNRC6 proteins, RAVER1 controls AS of GW-enriched motifs, which are essential for AGO2-binding and the formation of active miR/RISC complexes. We propose, that RAVER1 is a key modulator of AS events in the miR/RISC pathway ensuring proper abundance and composition of miR/RISC effectors. This ensures balanced expression of TGFB signaling effectors and limits TGFB induced lethal EMT.

Molecular Mechanism of Resistance to Alternaria alternata Apple Pathotype in Apple by Alternative Splicing of Transcription Factor MdMYB6-like

As a fruit tree with great economic value, apple is widely cultivated in China. However, apple leaf spot disease causes significant damage to apple quality and economic value. In our study, we found that MdMYB6-like is a transcription factor without auto-activation activity and with three alternative spliced variants. Among them, MdMYB6-like-β responded positively to the pathogen infection. Overexpression of MdMYB6-like-β increased the lignin content of leaves and improved the pathogenic resistance of apple flesh callus. In addition, all three alternative spliced variants of MdMYB6-like could bind to the promoter of MdBGLU H. Therefore, we believe that MdMYB6-like plays an important role in the infection process of the pathogen and lays a solid foundation for breeding disease-resistant cultivars of apple in the future.

Reference-informed prediction of alternative splicing and splicing-altering mutations from sequences

Alternative splicing plays a crucial role in protein diversity and gene expression regulation in higher eukaryotes and mutations causing dysregulated splicing underlie a range of genetic diseases. Computational prediction of alternative splicing from genomic sequences not only provides insight into gene-regulatory mechanisms but also helps identify disease-causing mutations and drug targets. However, the current methods for the quantitative prediction of splice site usage still have limited accuracy. Here, we present DeltaSplice, a deep neural network model optimized to learn the impact of mutations on quantitative changes in alternative splicing from the comparative analysis of homologous genes. The model architecture enables DeltaSplice to perform "reference-informed prediction" by incorporating the known splice site usage of a reference gene sequence to improve its prediction on splicing-altering mutations. We benchmarked DeltaSplice and several other state-of-the-art methods on various prediction tasks, including evolutionary sequence divergence on lineage-specific splicing and splicing-altering mutations in human populations and neurodevelopmental disorders, and demonstrated that DeltaSplice outperformed consistently. DeltaSplice predicted ~15% of splicing quantitative trait loci (sQTLs) in the human brain as causal splicing-altering variants. It also predicted splicing-altering de novo mutations outside the splice sites in a subset of patients affected by autism and other neurodevelopmental disorders, including 19 genes with recurrent splicing-altering mutations. Among the new candidate disease risk genes, MFN1 is involved in mitochondria fusion, which is frequently disrupted in autism patients. Our work expanded the capacity of in silico splicing models with potential applications in genetic diagnosis and the development of splicing-based precision medicine.

Identification of splicing factors signature predicting prognosis risk and the mechanistic roles of novel oncogenes in HNSCC

Head and neck squamous cell carcinoma (HNSCC) is the most frequent subtype of head and neck cancer, generally with a poor prognosis and limited therapeutic options due to its highly heterogeneous malignancy. In this study, we screened functional splicing regulatory RNA binding proteins (RBPs) that were closely related with the prognosis of HNSCC patients and showed significant expression differences between HNSCC tumors and normal tissues. Based on this finding, we chose six candidate genes (HNRNPC, ZCRB1, RBM12B, SF3A2, SF3B3, and SRSF11) to generate a prognostic prediction model and validated the accuracy of the prognostic model for predicting patient survival outcomes. We found that the risk score predicted by our model can serve as an independent prognostic predictor. Notably, HNSCC tumors showing higher expression of SF3B3, HNRNPC, or ZCRB1 possessed higher risk scores in the discovered prediction model. The investigation of the underlying mechanism validated that knockdown of SF3B3, HNRNPC, and ZCRB1 separately induced a substantial impairment of HNSCC cell survival. Conversely, overexpression of each of the three genes promoted tumor cellular proliferation. High throughput RNA sequencing analysis revealed that changes in the expression of SF3B3 and HNRNPC remarkably affected alternative splicing of genes related to cell cycle regulation, whereas the depletion of ZCRB1 contributed to aberrant splicing events involving in DNA damage response. In addition, the prognostic prediction model's risk score was demonstrated to be related with the immune infiltration score. Particularly, SF3B3 has a negative correlation with CD8A expression. Therefore, our findings provide promising prognosis predictors and potential therapeutic targets for better treatment efficacy of HNSCC.

RNA-binding proteins potentially regulate the alternative splicing of apoptotic genes during knee osteoarthritis progression

BACKGROUND

Alternative splicing (AS) is a principal mode of genetic regulation and one of the most widely used mechanisms to generate structurally and functionally distinct mRNA and protein variants. Dysregulation of AS may result in aberrant transcription and protein products, leading to the emergence of human diseases. Although considered important for regulating gene expression, genome-wide AS dysregulation, underlying mechanisms, and clinical relevance in knee osteoarthritis (OA) remain unelucidated. Therefore, in this study, we elucidated and validated AS events and their regulatory mechanisms during OA progression.

RESULTS

In this study, we identified differentially expressed genes between human OA and healthy meniscus samples. Among them, the OA-associated genes were primarily enriched in biological pathways such as extracellular matrix organization and ossification. The predominant OA-associated regulated AS (RAS) events were found to be involved in apoptosis during OA development. The expression of the apoptosis-related gene BCL2L13, XAF1, and NF2 were significantly different between OA and healthy meniscus samples. The construction of a covariation network of RNA-binding proteins (RBPs) and RAS genes revealed that differentially expressed RBP genes LAMA2 and CUL4B may regulate the apoptotic genes XAF1 and BCL2L13 to undergo AS events during OA progression. Finally, RT-qPCR revealed that CUL4B expression was significantly higher in OA meniscus samples than in normal controls and that the AS ratio of XAF1 was significantly different between control and OA samples; these findings were consistent with their expected expression and regulatory relationships.

CONCLUSIONS

Differentially expressed RBPs may regulate the AS of apoptotic genes during knee OA progression. XAF1 and its regulator, CUL4B, may serve as novel biomarkers and potential therapeutic targets for this disease.

Variant- and vaccination-specific alternative splicing profiles in SARS-CoV-2 infections

The COVID-19 pandemic, driven by the SARS-CoV-2 virus and its variants, highlights the important role of understanding host-viral molecular interactions influencing infection outcomes. Alternative splicing post-infection can impact both host responses and viral replication. We analyzed RNA splicing patterns in immune cells across various SARS-CoV-2 variants, considering immunization status. Using a dataset of 190 RNA-seq samples from our prior studies, we observed a substantial deactivation of alternative splicing and RNA splicing-related genes in COVID-19 patients. The alterations varied significantly depending on the infecting variant and immunization history. Notably, Alpha or Beta-infected patients differed from controls, while Omicron-infected patients displayed a splicing profile closer to controls. Particularly, vaccinated Omicron-infected individuals showed a distinct dynamic in alternative splicing patterns not widely shared among other groups. Our findings underscore the intricate interplay between SARS-CoV-2 variants, vaccination-induced immunity, and alternative splicing, emphasizing the need for further investigations to deepen understanding and guide therapeutic development.

Alternative splicing of CsWRKY21 positively regulates cold response in tea plant

Alternative splicing (AS) was an important post-transcriptional mechanism that involved in plant resistance to adversity stress. WRKY transcription factors function as transcriptional activators or repressors to modulate plant growth, development and stress response. However, the role of alternate splicing of WRKY in cold tolerance is poorly understood in tea plants. In this study, we found that the CsWRKY21 transcription factor, a member of the WRKY IId subfamily, was induced by low temperature. Subcellular localization and transcriptional activity assays showed that CsWRKY21 localized to the nucleus and had no transcriptional activation activity. Y1H and dual-luciferase reporter assays showed that CsWRKY21 suppressed expression of CsABA8H and CsUGT by binding with their promoters. Transient overexpression of CsABA8H and CsUGT reduced abscisic acid (ABA) content in tobacco leaves. Furthermore, we discovered that CsWRKY21 undergoes AS in the 5'UTR region. The AS transcript CsWRKY21-b was induced at low temperature, up to 6 folds compared to the control, while the full-length CsWRKY21-a transcript did not significantly change. Western blot analysis showed that the retention of introns in the 5'UTR region of CsWRKY21-b led to higher CsWRKY21 protein content. These results revealed that alternative splicing of CsWRKY21 involved in cold tolerance of tea plant by regulating the protein expression level and then regulating the content of ABA, and provide insights into molecular mechanisms of low temperature defense mediated by AS in tea plant.

Transcriptome Responses to Different Salinity Conditions in Litoditis marina, Revealed by Long-Read Sequencing

The marine nematode Litoditis marina is widely distributed in intertidal zones around the globe, yet the mechanisms underlying its broad adaptation to salinity remain elusive. In this study, we applied ONT long-read sequencing technology to unravel the transcriptome responses to different salinity conditions in L. marina. Through ONT sequencing under 3‱, 30‱ and 60‱ salinity environments, we obtained 131.78 G clean data and 26,647 non-redundant long-read transcripts, including 6464 novel transcripts. The DEGs obtained from the current ONT lrRNA-seq were highly correlated with those identified in our previously reported Illumina short-read RNA sequencing data. When we compared the 30‱ to the 3‱ salinity condition, we found that GO terms such as oxidoreductase activity, cation transmembrane transport and ion transmembrane transport were shared between the ONT lrRNA-seq and Illumina data. Similarly, GO terms including extracellular space, structural constituents of cuticle, substrate-specific channel activity, ion transport and substrate-specific transmembrane transporter activity were shared between the ONT and Illumina data under 60‱ compared to 30‱ salinity. In addition, we found that 79 genes significantly increased, while 119 genes significantly decreased, as the salinity increased. Furthermore, through the GO enrichment analysis of 214 genes containing DAS, in 30‱ compared to 3‱ salinity, we found that GO terms such as cellular component assembly and coenzyme biosynthetic process were enriched. Additionally, we observed that GO terms such as cellular component assembly and coenzyme biosynthetic process were also enriched in 60‱ compared to 30‱ salinity. Moreover, we found that 86, 125, and 81 genes that contained DAS were also DEGs, in comparisons between 30‱ and 3‱, 60‱ and 30‱, and 60‱ and 3‱ salinity, respectively. In addition, we demonstrated the landscape of alternative polyadenylation in marine nematode under different salinity conditions This report provides several novel insights for the further study of the mechanisms by which euryhalinity formed and evolved, and it might also contribute to the investigation of salinity dynamics induced by global climate change.

Intron retention as a productive mechanism in human MAPT: RNA species generated by retention of intron 3

BACKGROUND

Tau is a microtubule-binding protein encoded by the MAPT gene. Tau is essential for several physiological functions and associated with pathological processes, including Alzheimer's disease (AD). Six tau isoforms are typically described in the central nervous system, but current research paints a more diverse landscape and a more nuanced balance between isoforms. Recent work has described tau isoforms generated by intron 11 and intron 12 retention. This work adds to that evidence, proving the existence of MAPT transcripts retaining intron 3. Our aim is to demonstrate the existence of mature MAPT RNA species that retain intron 3 in human brain samples and to study its correlation with Alzheimer's disease across different regions.

METHODS

Initial evidence of intron-3-retaining MAPT species come from in silico analysis of RNA-seq databases. We further demonstrate the existence of these mature RNA species in a human neuroepithelioma cell line and human brain samples by quantitative PCR. We also use digital droplet PCR to demonstrate the existence of RNA species that retain either intron 3, intron 12 or both introns.

FINDINGS

Intron-3-retaining species are even more prominently present that intron-12-retaining ones. We show the presence of MAPT transcripts that retain both introns 3 and 12. These intron-retaining species are diminished in brain samples of patients with Alzheimer's disease with respect to individuals without dementia. Conversely, relative abundance of intron-3- or intron-12-retaining MAPT species with respect to double-retaining species as well as their percentage of expression with respect to total MAPT are increased in patients with Alzheimer's disease, especially in hippocampal samples. Among these TIR-MAPT species, TIR3+12 double truncation allows better classification potential of Alzheimer's disease samples. Moreover, we find a significant increase in intron-3- or intron-12-retaining species and its relative abundance with respect to double-retaining MAPT species in cerebellum in contrast to frontal lateral cortex and hippocampus in individuals with no signs of dementia.

INTERPRETATION

Intron retention constitutes a potential mechanism to generate Tau isoforms whose mature RNA expression levels correlate with Alzheimer's pathology showing its potential as a biomarker associated to the disease.

FUNDING

This research was funded by the Spanish Ministry of Science, Innovation and Universities: PGC2018-096177-B-I00 (J.A.); Spanish Ministry of Science and Innovation (MCIN): PID2020-113204GB-I00 (F.H.) and PID2021-123859OB-100 from MCIN/AEI/10.13039/501100011033/FEDER, UE (J.A.). It was also supported by CSIC through an intramural grant (201920E104) (J.A.) and the Centre for Networked Biomedical Research on Neurodegenerative Diseases (J.A.). The Centro de Biología Molecular Severo Ochoa (CBMSO) is a Severo Ochoa Center of Excellence (MICIN, award CEX2021-001154-S).

RedRibbon: A new rank-rank hypergeometric overlap for gene and transcript expression signatures

High-throughput omics technologies have generated a wealth of large protein, gene, and transcript datasets that have exacerbated the need for new methods to analyse and compare big datasets. Rank-rank hypergeometric overlap is an important threshold-free method to combine and visualize two ranked lists of P-values or fold-changes, usually from differential gene expression analyses. Here, we introduce a new rank-rank hypergeometric overlap-based method aimed at gene level and alternative splicing analyses at transcript or exon level, hitherto unreachable as transcript numbers are an order of magnitude larger than gene numbers. We tested the tool on synthetic and real datasets at gene and transcript levels to detect correlation and anticorrelation patterns and found it to be fast and accurate, even on very large datasets thanks to an evolutionary algorithm-based minimal P-value search. The tool comes with a ready-to-use permutation scheme allowing the computation of adjusted P-values at low time cost. The package compatibility mode is a drop-in replacement to previous packages. RedRibbon holds the promise to accurately extricate detailed information from large comparative analyses.

MAGOH promotes gastric cancer progression via hnRNPA1 expression inhibition-mediated RONΔ160/PI3K/AKT signaling pathway activation

BACKGROUND

Gastric cancer (GC) is associated with high mortality and heterogeneity and poses a great threat to humans. Gene therapies for the receptor tyrosine kinase RON and its spliceosomes are attracting increasing amounts of attention due to their unique characteristics. However, little is known about the mechanism involved in the formation of the RON mRNA alternative spliceosome RONΔ160.

METHODS

Fourteen human GC tissue samples and six normal gastric tissue samples were subjected to label-free relative quantitative proteomics analysis, and MAGOH was identified as a candidate protein for subsequent studies. The expression of MAGOH in clinical specimens was verified by quantitative real-time PCR and western blotting. We then determined the biological function of MAGOH in GC through in vitro and in vivo experiments. RNA pulldown, RNA sequencing and RNA immunoprecipitation (RIP) were subsequently conducted to uncover the underlying mechanism by which MAGOH regulated the formation of RONΔ160.

RESULTS

Proteomic analysis revealed that MAGOH, which is located at key nodes and participates in RNA processing and mRNA splicing, was upregulated in GC tissue and GC cell lines and was associated with poor prognosis. Functional analysis showed that MAGOH promoted the proliferation, migration and invasion of GC cells in vitro and in vivo. Mechanistically, MAGOH inhibited the expression of hnRNPA1 and reduced the binding of hnRNPA1 to RON mRNA, thereby promoting the formation of RONΔ160 to activate the PI3K/AKT signaling pathway and consequently facilitating GC progression.

CONCLUSIONS

Our study revealed that MAGOH could promote the formation of RONΔ160 and activate the PI3K/AKT signaling pathway through the inhibition of hnRNPA1 expression. We elucidate a novel mechanism and potential therapeutic targets for the growth and metastasis of GC based on the MAGOH-RONΔ160 axis, and these findings have important guiding significance and clinical value for the future development of effective therapeutic strategies for GC.

Diurnal regulation of alternative splicing associated with thermotolerance in rice by two glycine-rich RNA-binding proteins

Rice (Oryza sativa L.) production is threatened by global warming associated with extreme high temperatures, and rice heat sensitivity is differed when stress occurs between daytime and nighttime. However, the underlying molecular mechanism are largely unknown. We show here that two glycine-rich RNA binding proteins, OsGRP3 and OsGRP162, are required for thermotolerance in rice, especially at nighttime. The rhythmic expression of OsGRP3/OsGRP162 peaks at midnight, and at these coincident times, is increased by heat stress. This is largely dependent on the evening complex component OsELF3-2. We next found that the double mutant of OsGRP3/OsGRP162 is strikingly more sensitive to heat stress in terms of survival rate and seed setting rate when comparing to the wild-type plants. Interestingly, the defect in thermotolerance is more evident when heat stress occurred in nighttime than that in daytime. Upon heat stress, the double mutant of OsGRP3/OsGRP162 displays globally reduced expression of heat-stress responsive genes, and increases of mRNA alternative splicing dominated by exon-skipping. This study thus reveals the important role of OsGRP3/OsGRP162 in thermotolerance in rice, and unravels the mechanism on how OsGRP3/OsGRP162 regulate thermotolerance in a diurnal manner.

Spatial and Temporal Relationship between Epithelial-Mesenchymal Transition (EMT) and Stem Cells in Cancer

BACKGROUND

Epithelial-mesenchymal transition (EMT) is often linked with carcinogenesis. However, EMT is also important for embryo development and only reactivates in cancer. Connecting how EMT occurs during embryonic development and in cancer could help us further understand the root mechanisms of cancer diseases.

CONTENT

There are key regulatory elements that contribute to EMT and the induction and maintenance of stem cell properties during embryogenesis, tissue regeneration, and carcinogenesis. Here, we explore the implications of EMT in the different stages of embryogenesis and tissue development. We especially highlight the necessity of EMT in the mesodermal formation and in neural crest cells. Through EMT, these cells gain epithelial-mesenchymal plasticity (EMP). With this transition, crucial morphological changes occur to progress through the metastatic cascade as well as tissue regeneration after an injury. Stem-like cells, including cancer stem cells, are generated from EMT and during this process upregulate factors necessary for stem cell maintenance. Hence, it is important to understand the key regulators allowing stem cell awakening in cancer, which increases plasticity and promotes treatment resistance, to develop strategies targeting this cell population and improve patient outcomes.

SUMMARY

EMT involves multifaceted regulation to allow the fluidity needed to facilitate adaptation. This regulatory mechanism, plasticity, involves many cooperating transcription factors. Additionally, posttranslational modifications, such as splicing, activate the correct isoforms for either epithelial or mesenchymal specificity. Moreover, epigenetic regulation also occurs, such as acetylation and methylation. Downstream signaling ultimately results in the EMT which promotes tissue generation/regeneration and cancer progression.

Transcriptomic analysis of differentially alternative splicing patterns in mice with inflammatory and neuropathic pain

Although the molecular mechanisms of chronic pain have been extensively studied, a global picture of alternatively spliced genes and events in the peripheral and central nervous systems of chronic pain is poorly understood. The current study analyzed the changing pattern of alternative splicing (AS) in mouse brain, dorsal root ganglion, and spinal cord tissue under inflammatory and neuropathic pain. In total, we identified 6495 differentially alternatively spliced (DAS) genes. The molecular functions of shared DAS genes between these two models are mainly enriched in calcium signaling pathways, synapse organization, axon regeneration, and neurodegeneration disease. Additionally, we identified 509 DAS in differentially expressed genes (DEGs) shared by these two models, accounting for a small proportion of total DEGs. Our findings supported the hypothesis that the AS has an independent regulation pattern different from transcriptional regulation. Taken together, these findings indicate that AS is one of the important molecular mechanisms of chronic pain in mammals. This study presents a global description of AS profile changes in the full path of neuropathic and inflammatory pain models, providing new insights into the underlying mechanisms of chronic pain and guiding genomic clinical diagnosis methods and rational medication.