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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Alternative splicing of clock transcript mediates the response of circadian clocks to temperature changes

Circadian clocks respond to temperature changes over the calendar year, allowing organisms to adjust their daily biological rhythms to optimize health and fitness. In Drosophila, seasonal adaptations and temperature compensation are regulated by temperature-sensitive alternative splicing (AS) of period (per) and timeless (tim) genes that encode key transcriptional repressors of clock gene expression. Although clock (clk) gene encodes the critical activator of clock gene expression, AS of its transcripts and its potential role in temperature regulation of clock function have not been explored. We therefore sought to investigate whether clk exhibits AS in response to temperature and the functional changes of the differentially spliced transcripts. We observed that clk transcripts indeed undergo temperature-sensitive AS. Specifically, cold temperature leads to the production of an alternative clk transcript, hereinafter termed clk-cold, which encodes a CLK isoform with an in-frame deletion of four amino acids proximal to the DNA binding domain. Notably, serine 13 (S13), which we found to be a CK1α-dependent phosphorylation site, is among the four amino acids deleted in CLK-cold protein. Using a combination of transgenic fly, tissue culture, and in vitro experiments, we demonstrated that upon phosphorylation at CLK(S13), CLK-DNA interaction is reduced, thus decreasing CLK occupancy at clock gene promoters. This is in agreement with our findings that CLK occupancy at clock genes and transcriptional output are elevated at cold temperature, which can be explained by the higher amounts of CLK-cold isoforms that lack S13 residue. This study provides new insights into the complex collaboration between AS and phospho-regulation in shaping temperature responses of the circadian clock.

Characterization of alternative splicing events and prognostic signatures in gastric cancer

BACKGROUND

Accumulating evidences indicate that the specific alternative splicing (AS) events are linked to the occurrence and prognosis of gastric cancer (GC). Nevertheless, the impact of AS is still unclear and needed to further elucidation.

METHODS

The expression profile of GC and normal samples were downloaded from TCGA. AS events were achieved from SpliceSeq database. Cox regression together with LASSO analysis were employed to identify survival-associated AS events (SASEs) and calculate risk scores. PPI and pathway enrichment analysis were implemented to determine the function and pathways of these genes. Kaplan-Meier (K-M) analysis and Receiver Operating Characteristic Curves were used to evaluate the clinical significance of genes of SASEs. Q-PCR were applied to validate the hub genes on the survival prognosis in 47 GC samples. Drug sensitivity and immune cell infiltration analysis were conducted.

RESULTS

In total, 48 140 AS events in 10 610 genes from 361 GC and 31 normal samples were analyzed. Through univariate Cox regression, 855 SASEs in 763 genes were screened out. Further, these SASEs were analyzed by PPI and 17 hub genes were identified. Meanwhile, using Lasso and multivariate Cox regression analysis, 135 SASEs in 132 genes related to 7 AS forms were further screened and a GC prognostic model was constructed. K-M curves indicates that high-risk group has poorer prognosis. And the nomogram analysis on the basis of the multivariate Cox analysis was disclosed the interrelationships between 7 AS forms and clinical parameters in the model. Five key genes were then screened out by PPI analysis and Differential Expression Gene analysis based on TCGA and Combined-dataset, namely STAT3, RAD51B, SOCS2, POLE2 and TSR1. The expression levels of AS in STAT3, RAD51B, SOCS2, POLE2 and TSR1 were all significantly correlated with survival by qPCR verification. Nineteen drugs were sensitized to high-risk patients and eight immune cells showed significantly different infiltration between the STAD and normal groups.

CONCLUSIONS

In this research, the prognostic model constructed by SASEs can be applied to predict the prognosis of GC patients and the selected key genes are expected to become new biomarkers and therapeutical targets for GC treatment.

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

The advancement of Long-Read Sequencing (LRS) techniques has significantly increased the length of sequencing to several kilobases, thereby facilitating the identification of alternative splicing events and isoform expressions. Recently, numerous computational tools for isoform detection using long-read sequencing data have been developed. Nevertheless, there remains a deficiency in comparative studies that systemically evaluate the performance of these tools, which are implemented with different algorithms, under various simulations that encompass potential influencing factors. In this study, we conducted a benchmark analysis of thirteen methods implemented in nine tools capable of identifying isoform structures from long-read RNA-seq data. We evaluated their performances using simulated data, which represented diverse sequencing platforms generated by an in-house simulator, RNA sequins (sequencing spike-ins) data, as well as experimental data. Our findings demonstrate IsoQuant as a highly effective tool for isoform detection with LRS, with Bambu and StringTie2 also exhibiting strong performance. These results offer valuable guidance for future research on alternative splicing analysis and the ongoing improvement of tools for isoform detection using LRS data.

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.

Single-Cell Omics for Transcriptome CHaracterization (SCOTCH): isoform-level characterization of gene expression through long-read single-cell RNA sequencing

The advent of long-read single-cell transcriptome sequencing (lr-scRNA-Seq) represents a significant leap forward in single-cell genomics. With the recent introduction of R10 flowcells by Oxford Nanopore, we propose that previous computational methods designed to handle high sequencing error rates are no longer relevant, and that the prevailing approach using short reads to compile "barcode space" (candidate barcode list) to de-multiplex long reads are no longer necessary. Instead, computational methods should now shift focus on harnessing the unique benefits of long reads to analyze transcriptome complexity. In this context, we introduce a comprehensive suite of computational methods named Single-Cell Omics for Transcriptome CHaracterization (SCOTCH). Our method is compatible with the single-cell library preparation platform from both 10X Genomics and Parse Biosciences, facilitating the analysis of special cell populations, such as neurons, hepatocytes and developing cardiomyocytes. We specifically re-formulated the transcript mapping problem with a compatibility matrix and addressed the multiple-mapping issue using probabilistic inference, which allows the discovery of novel isoforms as well as the detection of differential isoform usage between cell populations. We evaluated SCOTCH through analysis of real data across different combinations of single-cell libraries and sequencing technologies (10X + Illumina, Parse + Illumina, 10X + Nanopore_R9, 10X + Nanopore_R10, Parse + Nanopore_R10), and showed its ability to infer novel biological insights on cell type-specific isoform expression. These datasets enhance the availability of publicly available data for continued development of computational approaches. In summary, SCOTCH allows extraction of more biological insights from the new advancements in single-cell library construction and sequencing technologies, facilitating the examination of transcriptome complexity at the single-cell level.

CDK9 inhibition activates innate immune response through viral mimicry

Cancer cells frequently exhibit hyperactivation of transcription, which can lead to increased sensitivity to compounds targeting the transcriptional kinases, in particular CDK9. However, mechanistic details of CDK9 inhibition-induced cancer cell-selective anti-proliferative effects remain largely unknown. Here, we discover that CDK9 inhibition activates the innate immune response through viral mimicry in cancer cells. In MYC over-expressing prostate cancer cells, CDK9 inhibition leads to the gross accumulation of mis-spliced RNA. Double-stranded RNA (dsRNA)-activated kinase can recognize these mis-spliced RNAs, and we show that the activity of this kinase is required for the CDK9 inhibitor-induced anti-proliferative effects. Using time-resolved transcriptional profiling (SLAM-seq), targeted proteomics, and ChIP-seq, we show that, similar to viral infection, CDK9 inhibition significantly suppresses transcription of most genes but allows selective transcription and translation of cytokines related to the innate immune response. In particular, CDK9 inhibition activates NFκB-driven cytokine signaling at the transcriptional and secretome levels. The transcriptional signature induced by CDK9 inhibition identifies prostate cancers with a high level of genome instability. We propose that it is possible to induce similar effects in patients using CDK9 inhibition, which, we show, causes DNA damage in vitro. In the future, it is important to establish whether CDK9 inhibitors can potentiate the effects of immunotherapy against late-stage prostate cancer, a currently lethal disease.

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 residue of H3 ( i.e. , glutamine) and dopamine. In our previous studies, we discovered that the dynamics of this post-translational modification (including installation, removal, and replacement) were regulated by a single enzyme, transglutaminase 2 (TGM2), through reversible transamination. 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 TGM2 in colon cancer cells. Due to the enzyme promiscuity of TGM2, non-histone proteins have also been identified as targets of dopaminylation on glutamine residues, 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 modification sites of these dopaminylated proteins were identified, 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 to block the installation of protein dopaminylation may become a promising anti-cancer strategy.

TOC

Proteomics Identifies LUC7L3 as a Prognostic Biomarker for Hepatocellular Carcinoma

Alternative splicing has been shown to participate in tumor progression, including hepatocellular carcinoma. The poor prognosis of patients with HCC calls for molecular classification and biomarker identification to facilitate precision medicine. We performed ssGSEA analysis to quantify the pathway activity of RNA splicing in three HCC cohorts. Kaplan-Meier and Cox methods were used for survival analysis. GO and GSEA were performed to analyze pathway enrichment. We confirmed that RNA splicing is significantly correlated with prognosis, and identified an alternative splicing-associated protein LUC7L3 as a potential HCC prognostic biomarker. Further bioinformatics analysis revealed that high LUC7L3 expression indicated a more progressive HCC subtype and worse clinical features. Cell proliferation-related pathways were enriched in HCC patients with high LUC7L3 expression. Consistently, we proved that LUC7L3 knockdown could significantly inhibit cell proliferation and suppress the activation of associated signaling pathways in vitro. In this research, the relevance between RNA splicing and HCC patient prognosis was outlined. Our newly identified biomarker LUC7L3 could provide stratification for patient survival and recurrence risk, facilitating early medical intervention before recurrence or disease progression.

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.

A super minigene with a short promoter and truncated introns recapitulates essential features of transcription and splicing regulation of the SMN1 and SMN2 genes

Here we report a Survival Motor Neuron 2 (SMN2) super minigene, SMN2Sup, encompassing its own promoter, all exons, their flanking intronic sequences and the entire 3'-untranslated region. We confirm that the pre-mRNA generated from SMN2Sup undergoes splicing to produce a translation-competent mRNA. We demonstrate that mRNA generated from SMN2Sup produces more SMN than an identical mRNA generated from a cDNA clone. We uncover that overexpression of SMN triggers skipping of exon 3 of SMN1/SMN2. We define the minimal promoter and regulatory elements associated with the initiation and elongation of transcription of SMN2. The shortened introns within SMN2Sup preserved the ability of camptothecin, a transcription elongation inhibitor, to induce skipping of exons 3 and 7 of SMN2. We show that intron 1-retained transcripts undergo nonsense-mediated decay. We demonstrate that splicing factor SRSF3 and DNA/RNA helicase DHX9 regulate splicing of multiple exons in the context of both SMN2Sup and endogenous SMN1/SMN2. Prevention of SMN2 exon 7 skipping has implications for the treatment of spinal muscular atrophy (SMA). We validate the utility of the super minigene in monitoring SMN levels upon splicing correction. Finally, we demonstrate how the super minigene could be employed to capture the cell type-specific effects of a pathogenic SMN1 mutation.

Genomic deletions explain the generation of alternative BRAF isoforms conferring resistance to MAPK inhibitors in melanoma

Resistance to MAPK inhibitors (MAPKi), the main cause of relapse in BRAF-mutant melanoma, is associated with the production of alternative BRAF mRNA isoforms (altBRAFs) in up to 30% of patients receiving BRAF inhibitor monotherapy. These altBRAFs have been described as being generated by alternative pre-mRNA splicing, and splicing modulation has been proposed as a therapeutic strategy to overcome resistance. In contrast, we report that altBRAFs are generated through genomic deletions. Using different in vitro models of altBRAF-mediated melanoma resistance, we demonstrate the production of altBRAFs exclusively from the BRAF V600E allele, correlating with corresponding genomic deletions. Genomic deletions are also detected in tumor samples from melanoma and breast cancer patients expressing altBRAFs. Along with the identification of altBRAFs in BRAF wild-type and in MAPKi-naive melanoma samples, our results represent a major shift in our understanding of mechanisms leading to the generation of BRAF transcripts variants associated with resistance in melanoma.

Stochasticity in genetics and gene regulation

Development from fertilized egg to functioning multi-cellular organism requires precision. There is no precision, and often no survival, without plasticity. Plasticity is conferred partly by stochastic variation, present inherently in all biological systems. Gene expression levels fluctuate ubiquitously through transcription, alternative splicing, translation and turnover. Small differences in gene expression are exploited to trigger early differentiation, conferring distinct function on selected individual cells and setting in motion regulatory interactions. Non-selected cells then acquire new functions along the spatio-temporal developmental trajectory. The differentiation process has many stochastic components. Meiotic segregation, mitochondrial partitioning, X-inactivation and the dynamic DNA binding of transcription factor assemblies-all exhibit randomness. Non-random X-inactivation generally signals deleterious X-linked mutations. Correct neural wiring, such as retina to brain, arises through repeated confirmatory activity of connections made randomly. In immune system development, both B-cell antibody generation and the emergence of balanced T-cell categories begin through stochastic trial and error followed by functional selection. Aberrant selection processes lead to immune dysfunction. DNA sequence variants also arise through stochastic events: some involving environmental fluctuation (radiation or presence of pollutants), or genetic repair system malfunction. The phenotypic outcome of mutations is also fluid. Mutations may be advantageous in some circumstances, deleterious in others. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

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.

Lipotoxicity as a therapeutic target in obesity and diabetic cardiomyopathy

Unhealthy sources of fats, ultra-processed foods with added sugars, and a sedentary lifestyle make humans more susceptible to developing overweight and obesity. While lipids constitute an integral component of the organism, excessive and abnormal lipid accumulation that exceeds the storage capacity of lipid droplets disrupts the intracellular composition of fatty acids and results in the release of deleterious lipid species, thereby giving rise to a pathological state termed lipotoxicity. This condition induces endoplasmic reticulum stress, mitochondrial dysfunction, inflammatory responses, and cell death. Recent advances in omics technologies and analytical methodologies and clinical research have provided novel insights into the mechanisms of lipotoxicity, including gut dysbiosis, epigenetic and epitranscriptomic modifications, dysfunction of lipid droplets, post-translational modifications, and altered membrane lipid composition. In this review, we discuss the recent knowledge on the mechanisms underlying the development of lipotoxicity and lipotoxic cardiometabolic disease in obesity, with a particular focus on lipotoxic and diabetic cardiomyopathy.

Comprehensive Long-Read Sequencing Analysis Discloses the Transcriptome Features of Papillary Thyroid Microcarcinoma

CONTEXT

Papillary thyroid microcarcinoma (PTMC) is the most common type of thyroid cancer. It has been shown that lymph node metastasis is associated with poor prognosis in patients with PTMC.

OBJECTIVE

We aim to characterize the PTMC transcriptome landscape and identify the candidate transcripts that are associated with lateral neck lymph node metastasis of PTMC.

METHODS

We performed full-length transcriptome sequencing in 64 PTMC samples. Standard bioinformatic pipelines were applied to characterize and annotate the full-length expression profiles of 2 PTMC subtypes. Functional open reading frame (ORF) annotation of the known and novel transcripts were predicted by HMMER, DeepLoc, and DeepTMHMM tools. Candidate transcripts associated with the pN1b subtype were identified after transcript quantification and differential gene expression analyses.

RESULTS

We found that skipping exons accounted for the more than 27.82% of the alternative splicing events. At least 42.56% of the discovered transcripts were novel isoforms of annotated genes. A total of 39 193 ORFs in novel transcripts and 18 596 ORFs in known transcripts were identified. Distribution patterns of the characterized transcripts in functional domain, subcellular localization, and transmembrane structure were predicted. In total, 1033 and 1204 differentially expressed genes were identified in the pN0 and pN1b groups, respectively. Moreover, novel isoforms of FRMD3, NOD1, and SHROOM4 were highlighted for their association with pN1b subtype.

CONCLUSION

Our data provided the global transcriptome landscape of PTMC and also revealed the novel isoforms that associated with PTMC aggressiveness.

The splicing regulators RBM5 and RBM10 are subunits of the U2 snRNP engaged with intron branch sites on chromatin

Understanding the mechanisms of pre-mRNA splicing is limited by the technical challenges to examining spliceosomes in vivo. Here, we report the isolation of RNP complexes derived from precatalytic A or B-like spliceosomes solubilized from the chromatin pellet of mammalian cell nuclei. We found that these complexes contain U2 snRNP proteins and a portion of the U2 snRNA bound with protected RNA fragments that precisely map to intronic branch sites across the transcriptome. These U2 complexes also contained the splicing regulators RBM5 and RBM10. We found RBM5 and RBM10 bound to nearly all branch site complexes and not simply those at regulated exons. The deletion of a conserved RBM5/RBM10 peptide sequence, including a zinc finger motif, disrupted U2 interaction and rendered the proteins inactive for the repression of many alternative exons. We propose a model where RBM5 and RBM10 regulate splicing as components of the U2 snRNP complex following branch site base pairing.

Dual RNA sequencing of Helicobacter pylori and host cell transcriptomes reveals ontologically distinct host-pathogen interaction

Helicobacter pylori is a highly successful pathogen that poses a substantial threat to human health. However, the dynamic interaction between H. pylori and the human gastric epithelium has not been fully investigated. In this study, using dual RNA sequencing technology, we characterized a cytotoxin-associated gene A (cagA)-modulated bacterial adaption strategy by enhancing the expression of ATP-binding cassette transporter-related genes, metQ and HP_0888, upon coculturing with human gastric epithelial cells. We observed a general repression of electron transport-associated genes by cagA, leading to the activation of oxidative phosphorylation. Temporal profiling of host mRNA signatures revealed the downregulation of multiple splicing regulators due to bacterial infection, resulting in aberrant pre-mRNA splicing of functional genes involved in the cell cycle process in response to H. pylori infection. Moreover, we demonstrated a protective effect of gastric H. pylori colonization against chronic dextran sulfate sodium (DSS)-induced colitis. Mechanistically, we identified a cluster of propionic and butyric acid-producing bacteria, Muribaculaceae, selectively enriched in the colons of H. pylori-pre-colonized mice, which may contribute to the restoration of intestinal barrier function damaged by DSS treatment. Collectively, this study presents the first dual-transcriptome analysis of H. pylori during its dynamic interaction with gastric epithelial cells and provides new insights into strategies through which H. pylori promotes infection and pathogenesis in the human gastric epithelium.

IMPORTANCE

Simultaneous profiling of the dynamic interaction between Helicobacter pylori and the human gastric epithelium represents a novel strategy for identifying regulatory responses that drive pathogenesis. This study presents the first dual-transcriptome analysis of H. pylori when cocultured with gastric epithelial cells, revealing a bacterial adaptation strategy and a general repression of electron transportation-associated genes, both of which were modulated by cytotoxin-associated gene A (cagA). Temporal profiling of host mRNA signatures dissected the aberrant pre-mRNA splicing of functional genes involved in the cell cycle process in response to H. pylori infection. We demonstrated a protective effect of gastric H. pylori colonization against chronic DSS-induced colitis through both in vitro and in vivo experiments. These findings significantly enhance our understanding of how H. pylori promotes infection and pathogenesis in the human gastric epithelium and provide evidence to identify targets for antimicrobial therapies.

The Experimental and In Silico-Based Evaluation of NRF2 Modulators, Sulforaphane and Brusatol, on the Transcriptome of Immortalized Bovine Mammary Alveolar Cells

Changes during the production cycle of dairy cattle can leave these animals susceptible to oxidative stress and reduced antioxidant health. In particular, the periparturient period, when dairy cows must rapidly adapt to the sudden metabolic demands of lactation, is a period when the production of damaging free radicals can overwhelm the natural antioxidant systems, potentially leading to tissue damage and reduced milk production. Central to the protection against free radical damage and antioxidant defense is the transcription factor NRF2, which activates an array of genes associated with antioxidant functions and cell survival. The objective of this study was to evaluate the effect that two natural NRF2 modulators, the NRF2 agonist sulforaphane (SFN) and the antagonist brusatol (BRU), have on the transcriptome of immortalized bovine mammary alveolar cells (MACT) using both the RT-qPCR of putative NRF2 target genes, as well as RNA sequencing approaches. The treatment of cells with SFN resulted in the activation of many putative NRF2 target genes and the upregulation of genes associated with pathways involved in cell survival, metabolism, and antioxidant function while suppressing the expression of genes related to cellular senescence and DNA repair. In contrast, the treatment of cells with BRU resulted in the upregulation of genes associated with inflammation, cellular stress, and apoptosis while suppressing the transcription of genes involved in various metabolic processes. The analysis also revealed several novel putative NRF2 target genes in bovine. In conclusion, these data indicate that the treatment of cells with SFN and BRU may be effective at modulating the NRF2 transcriptional network, but additional effects associated with cellular stress and metabolism may complicate the effectiveness of these compounds to improve antioxidant health in dairy cattle via nutrigenomic approaches.

Aberrant splicing of CaV1.2 calcium channel induced by decreased Rbfox1 enhances arterial constriction during diabetic hyperglycemia

Diabetic hyperglycemia induces dysfunctions of arterial smooth muscle, leading to diabetic vascular complications. The CaV1.2 calcium channel is one primary pathway for Ca2+ influx, which initiates vasoconstriction. However, the long-term regulation mechanism(s) for vascular CaV1.2 functions under hyperglycemic condition remains unknown. Here, Sprague-Dawley rats fed with high-fat diet in combination with low dose streptozotocin and Goto-Kakizaki (GK) rats were used as diabetic models. Isolated mesenteric arteries (MAs) and vascular smooth muscle cells (VSMCs) from rat models were used to assess K+-induced arterial constriction and CaV1.2 channel functions using vascular myograph and whole-cell patch clamp, respectively. K+-induced vasoconstriction is persistently enhanced in the MAs from diabetic rats, and CaV1.2 alternative spliced exon 9* is increased, while exon 33 is decreased in rat diabetic arteries. Furthermore, CaV1.2 channels exhibit hyperpolarized current-voltage and activation curve in VSMCs from diabetic rats, which facilitates the channel function. Unexpectedly, the application of glycated serum (GS), mimicking advanced glycation end-products (AGEs), but not glucose, downregulates the expression of the splicing factor Rbfox1 in VSMCs. Moreover, GS application or Rbfox1 knockdown dynamically regulates alternative exons 9* and 33, leading to facilitated functions of CaV1.2 channels in VSMCs and MAs. Notably, GS increases K+-induced intracellular calcium concentration of VSMCs and the vasoconstriction of MAs. These results reveal that AGEs, not glucose, long-termly regulates CaV1.2 alternative splicing events by decreasing Rbfox1 expression, thereby enhancing channel functions and increasing vasoconstriction under diabetic hyperglycemia. This study identifies the specific molecular mechanism for enhanced vasoconstriction under hyperglycemia, providing a potential target for managing diabetic vascular complications.

Zygotic Splicing Activation of the Transcriptome is a Crucial Aspect of Maternal-to-Zygotic Transition and Required for the Conversion from Totipotency to Pluripotency

During maternal-to-zygotic transition (MZT) in the embryo, mRNA undergoes complex post-transcriptional regulatory processes. However, it is unclear whether and how alternative splicing plays a functional role in MZT. By analyzing transcriptome changes in mouse and human early embryos, dynamic changes in alternative splicing during MZT are observed and a previously unnoticed process of zygotic splicing activation (ZSA) following embryonic transcriptional activation is described. As the underlying mechanism of RNA splicing, splicing factors undergo dramatic maternal-to-zygotic conversion. This conversion relies on the key maternal factors BTG4 and PABPN1L and is zygotic-transcription-dependent. CDK11-dependent phosphorylation of the key splicing factor, SF3B1, and its aggregation with SRSF2 in the subnuclear domains of 2-cell embryos are prerequisites for ZSA. Isoforms generated by erroneous splicing, such as full-length Dppa4, hinder normal embryonic development. Moreover, alternative splicing regulates the conversion of early embryonic blastomeres from totipotency to pluripotency, thereby affecting embryonic lineage differentiation. ZSA is an essential post-transcriptional process of MZT and has physiological significance in generating new life. In addition to transcriptional activation, appropriate expression of transcript isoforms is also necessary for preimplantation embryonic development.

Understanding age-related pathologic changes in TDP-43 functions and the consequence on RNA splicing and signalling in health and disease

TAR DNA binding protein-43 (TDP-43) is a key component in RNA splicing which plays a crucial role in the aging process. In neurodegenerative diseases such as amyotrophic lateral sclerosis, frontotemporal dementia and limbic-predominant age-related TDP-43 encephalopathy, TDP-43 can be mutated, mislocalised out of the nucleus of neurons and glial cells and form cytoplasmic inclusions. These TDP-43 alterations can lead to its RNA splicing dysregulation and contribute to mis-splicing of various types of RNA, such as mRNA, microRNA, and circular RNA. These changes can result in the generation of an altered transcriptome and proteome within cells, ultimately changing the diversity and quantity of gene products. In this review, we summarise the findings of novel atypical RNAs resulting from TDP-43 dysfunction and their potential as biomarkers or targets for therapeutic development.

Loss-of-function mutation in PRMT9 causes abnormal synapse development by dysregulation of RNA alternative splicing

Protein arginine methyltransferase 9 (PRMT9) is a recently identified member of the PRMT family, yet its biological function remains largely unknown. Here, by characterizing an intellectual disability associated PRMT9 mutation (G189R) and establishing a Prmt9 conditional knockout (cKO) mouse model, we uncover an important function of PRMT9 in neuronal development. The G189R mutation abolishes PRMT9 methyltransferase activity and reduces its protein stability. Knockout of Prmt9 in hippocampal neurons causes alternative splicing of ~1900 genes, which likely accounts for the aberrant synapse development and impaired learning and memory in the Prmt9 cKO mice. Mechanistically, we discover a methylation-sensitive protein-RNA interaction between the arginine 508 (R508) of the splicing factor 3B subunit 2 (SF3B2), the site that is exclusively methylated by PRMT9, and the pre-mRNA anchoring site, a cis-regulatory element that is critical for RNA splicing. Additionally, using human and mouse cell lines, as well as an SF3B2 arginine methylation-deficient mouse model, we provide strong evidence that SF3B2 is the primary methylation substrate of PRMT9, thus highlighting the conserved function of the PRMT9/SF3B2 axis in regulating pre-mRNA splicing.

Alternative mRNA splicing events and regulators in epidermal differentiation

Alternative splicing (AS) of messenger RNAs occurs in ∼95% of multi-exon human genes and generates diverse RNA and protein isoforms. We investigated AS events associated with human epidermal differentiation, a process crucial for skin function. We identified 6,413 AS events, primarily involving cassette exons. We also predicted 34 RNA-binding proteins (RBPs) regulating epidermal AS, including 19 previously undescribed candidate regulators. From these results, we identified FUS as an RBP that regulates the balance between keratinocyte proliferation and differentiation. Additionally, we characterized the function of a cassette exon AS event in MAP3K7, which encodes a kinase involved in cell signaling. We found that a switch from the short to long isoform of MAP3K7, triggered during differentiation, enforces the demarcation between proliferating basal progenitors and overlying differentiated strata. Our findings indicate that AS occurs extensively in the human epidermis and has critical roles in skin homeostasis.