PSI-Sigma: a comprehensive splicing-detection method for short-read and long-read RNA-seq analysis

rMATS-turbo: an efficient and flexible computational tool for alternative splicing analysis of large-scale RNA-seq data

Pre-mRNA alternative splicing is a prevalent mechanism for diversifying eukaryotic transcriptomes and proteomes. Regulated alternative splicing plays a role in many biological processes, and dysregulated alternative splicing is a feature of many human diseases. Short-read RNA sequencing (RNA-seq) is now the standard approach for transcriptome-wide analysis of alternative splicing. Since 2011, our laboratory has developed and maintained Replicate Multivariate Analysis of Transcript Splicing (rMATS), a computational tool for discovering and quantifying alternative splicing events from RNA-seq data. Here we provide a protocol for the contemporary version of rMATS, rMATS-turbo, a fast and scalable re-implementation that maintains the statistical framework and user interface of the original rMATS software, while incorporating a revamped computational workflow with a substantial improvement in speed and data storage efficiency. The rMATS-turbo software scales up to massive RNA-seq datasets with tens of thousands of samples. To illustrate the utility of rMATS-turbo, we describe two representative application scenarios. First, we describe a broadly applicable two-group comparison to identify differential alternative splicing events between two sample groups, including both annotated and novel alternative splicing events. Second, we describe a quantitative analysis of alternative splicing in a large-scale RNA-seq dataset (~1,000 samples), including the discovery of alternative splicing events associated with distinct cell states. We detail the workflow and features of rMATS-turbo that enable efficient parallel processing and analysis of large-scale RNA-seq datasets on a compute cluster. We anticipate that this protocol will help the broad user base of rMATS-turbo make the best use of this software for studying alternative splicing in diverse biological systems.

Full-length transcriptome sequencing provides insights into alternative splicing under cold stress in peanut

Introduction

Peanut (Arachis hypogaea L.), also called groundnut is an important oil and cash crop grown widely in the world. The annual global production of groundnuts has increased to approximately 50 million tons, which provides a rich source of vegetable oils and proteins for humans. Low temperature (non-freezing) is one of the major factors restricting peanut growth, yield, and geographic distribution. Since the complexity of cold-resistance trait, the molecular mechanism of cold tolerance and related gene networks were largely unknown in peanut.

Methods

In this study, comparative transcriptomic analysis of two peanut cultivars (SLH vs. ZH12) with differential cold tolerance under low temperature (10°C) was performed using Oxford Nanopore Technology (ONT) platform.

Results and discussion

As a result, we identified 8,949 novel gene loci and 95,291 new/novel isoforms compared with the reference database. More differentially expressed genes (DEGs) were discovered in cold-sensitive cultivar (ZH12) than cold-tolerant cultivar (SLH), while more alternative splicing events were found in SLH compared to ZH12. Gene Ontology (GO) analyses of the common DEGs showed that the "response to stress", "chloroplast part", and "transcription factor activity" were the most enriched GO terms, indicating that photosynthesis process and transcription factors play crucial roles in cold stress response in peanut. We also detected a total of 708 differential alternative splicing genes (DASGs) under cold stress compared to normal condition. Intron retention (IR) and exon skipping (ES) were the most prevalent alternative splicing (AS) events. In total, 4,993 transcription factors and 292 splicing factors were detected, many of them had differential expression levels and/or underwent AS events in response to cold stress. Overexpression of two candidate genes (encoding trehalose-6-phosphatephosphatases, AhTPPs) in yeast improves cold tolerance. This study not only provides valuable resources for the study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

CD8+ T cell metabolic flexibility elicited by CD28-ARS2 axis-driven alternative splicing of PKM supports antitumor immunity

Metabolic flexibility has emerged as a critical determinant of CD8+ T-cell antitumor activity, yet the mechanisms driving the metabolic flexibility of T cells have not been determined. In this study, we investigated the influence of the nuclear cap-binding complex (CBC) adaptor protein ARS2 on mature T cells. In doing so, we discovered a novel signaling axis that endows activated CD8+ T cells with flexibility of glucose catabolism. ARS2 upregulation driven by CD28 signaling reinforced splicing factor recruitment to pre-mRNAs and affected approximately one-third of T-cell activation-induced alternative splicing events. Among these effects, the CD28-ARS2 axis suppressed the expression of the M1 isoform of pyruvate kinase in favor of PKM2, a key determinant of CD8+ T-cell glucose utilization, interferon gamma production, and antitumor effector function. Importantly, PKM alternative splicing occurred independently of CD28-driven PI3K pathway activation, revealing a novel means by which costimulation reprograms glucose metabolism in CD8+ T cells.

Genome-wide analyses reveal the regulatory roles of DNA methylation-regulated alternative promoter transcripts in breast cancer

A certain proportion of genes are regulated by multiple, distinct promoters, revealing a dynamic landscape of the cancer transcriptome. However, the contribution of alternative promoters (APs) in breast cancer (BRCA) remains largely unexplored. Here, we identified 3654 genes with multiple promoters in BRCA patients, and 53 of them could generate distinct AP transcripts that are dysregulated and prognosis-related in BRCA, namely prognosis-related dysregulated AP (prdeAP) transcripts. Interestingly, when we searched for the genomic signatures of these prdeAP genes, we found that the promoter regions of 92% of the prdeAP genes were enriched with abundant DNA methylation signals. Through further bioinformatic analysis and experimental validation, we showed that AP selections of TANK, UNKL, CCL28, and MAP1LC3A were regulated by DNA methylation upon their corresponding promoter regions. Functionally, by overexpressing AP variants of TANK, we found that TANK|55731 could dramatically suppress MDA-MB-231 cell proliferation and migration. Meanwhile, pan-cancer survival analyses suggested that AP variants of TANK provided more accurate prognostic predictive ability than TANK gene in a variety of tumor types, including BRCA. Together, by uncovering the DNA methylation-regulated AP transcripts with tumor prognostic features, our work revealed a novel layer of regulators in BRCA progression and provided potential targets that served as effective biomarkers for anti-BRCA treatment.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS.

Direct androgen receptor control of sexually dimorphic gene expression in the mammalian kidney

Mammalian organs exhibit distinct physiology, disease susceptibility, and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA sequencing (RNA-seq) data demonstrated that sex differences were established from 4 and 8 weeks after birth under gonadal control. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR)-mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation, whereas analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, disease, and metabolic linkage of sexually dimorphic gene activity.

Identification of Splice Variants and Isoforms in Transcriptomics and Proteomics

Alternative splicing is pivotal to the regulation of gene expression and protein diversity in eukaryotic cells. The detection of alternative splicing events requires specific omics technologies. Although short-read RNA sequencing has successfully supported a plethora of investigations on alternative splicing, the emerging technologies of long-read RNA sequencing and top-down mass spectrometry open new opportunities to identify alternative splicing and protein isoforms with less ambiguity. Here, we summarize improvements in short-read RNA sequencing for alternative splicing analysis, including percent splicing index estimation and differential analysis. We also review the computational methods used in top-down proteomics analysis regarding proteoform identification, including the construction of databases of protein isoforms and statistical analyses of search results. While many improvements in sequencing and computational methods will result from emerging technologies, there should be future endeavors to increase the effectiveness, integration, and proteome coverage of alternative splicing events.

Transcript capture and ultradeep long-read RNA sequencing (CAPLRseq) to diagnose HNPCC/Lynch syndrome

PURPOSE

Whereas most human genes encode multiple mRNA isoforms with distinct function, clinical workflows for assessing this heterogeneity are not readily available. This is a substantial shortcoming, considering that up to 25% of disease-causing gene variants are suspected of disrupting mRNA splicing or mRNA abundance. Long-read sequencing can readily portray mRNA isoform diversity, but its sensitivity is relatively low due to insufficient transcriptome penetration.

METHODS

We developed and applied capture-based target enrichment from patient RNA samples combined with Oxford Nanopore long-read sequencing for the analysis of 123 hereditary cancer transcripts (capture and ultradeep long-read RNA sequencing (CAPLRseq)).

RESULTS

Validating CAPLRseq, we confirmed 17 cases of hereditary non-polyposis colorectal cancer/Lynch syndrome based on the demonstration of splicing defects and loss of allele expression of mismatch repair genes MLH1, PMS2, MSH2 and MSH6. Using CAPLRseq, we reclassified two variants of uncertain significance in MSH6 and PMS2 as either likely pathogenic or benign.

CONCLUSION

Our data show that CAPLRseq is an automatable and adaptable workflow for effective transcriptome-based identification of disease variants in a clinical diagnostic setting.

Integrated multiomic analysis and high-throughput screening reveal potential gene targets and synergetic drug combinations for osteosarcoma therapy

Although great advances have been made over the past decades, therapeutics for osteosarcoma are quite limited. We performed long-read RNA sequencing and tandem mass tag (TMT)-based quantitative proteome on osteosarcoma and the adjacent normal tissues, next-generation sequencing (NGS) on paired osteosarcoma samples before and after neoadjuvant chemotherapy (NACT), and high-throughput drug combination screen on osteosarcoma cell lines. Single-cell RNA sequencing data were analyzed to reveal the heterogeneity of potential therapeutic target genes. Additionally, we clarified the synergistic mechanisms of doxorubicin (DOX) and HDACs inhibitors for osteosarcoma treatment. Consequently, we identified 2535 osteosarcoma-specific genes and several alternative splicing (AS) events with osteosarcoma specificity and/or patient heterogeneity. Hundreds of potential therapeutic targets were identified among them, which showed the core regulatory roles in osteosarcoma. We also identified 215 inhibitory drugs and 236 synergistic drug combinations for osteosarcoma treatment. More interestingly, the multiomic analysis pointed out the pivotal role of HDAC1 and TOP2A in osteosarcoma. HDAC inhibitors synergized with DOX to suppress osteosarcoma both in vitro and in vivo. Mechanistically, HDAC inhibitors synergized with DOX by downregulating SP1 to transcriptionally modulate TOP2A expression. This study provided a comprehensive view of molecular features, therapeutic targets, and synergistic drug combinations for osteosarcoma.

FASTdRNA: a workflow for the analysis of ONT direct RNA sequencing

Motivation

Direct RNA-seq (dRNA-seq) using Oxford Nanopore Technology (ONT) has revolutionized transcript mapping by offering enhanced precision due to its long-read length. Unlike traditional techniques, dRNA-seq eliminates the need for PCR amplification, reducing the impact of GC bias, and preserving valuable base physical information, such as RNA modification and poly(A) length estimation. However, the rapid advancement of ONT devices has set higher standards for analytical software, resulting in potential challenges of software incompatibility and reduced efficiency.

Results

We present a novel workflow, called FASTdRNA, to manipulate dRNA-seq data efficiently. This workflow comprises two modules: a data preprocessing module and a data analysis module. The preprocessing data module, dRNAmain, encompasses basecalling, mapping, and transcript counting, which are essential for subsequent analyses. The data analysis module consists of a range of downstream analyses that facilitate the estimation of poly(A) length, prediction of RNA modifications, and assessment of alternative splicing events across different conditions with duplication. The FASTdRNA workflow is designed for the Snakemake framework and can be efficiently executed locally or in the cloud. Comparative experiments have demonstrated its superior performance compared to previous methods. This innovative workflow enhances the research capabilities of dRNA-seq data analysis pipelines by optimizing existing processes and expanding the scope of analysis.

Availability and implementation

The workflow is freely available at https://github.com/Tomcxf/FASTdRNA under an MIT license. Detailed install and usage guidance can be found in the GitHub repository.

Splicing Factor SRSF1 Promotes Pancreatitis and KRASG12D-Mediated Pancreatic Cancer

Inflammation is strongly associated with pancreatic ductal adenocarcinoma (PDAC), a highly lethal malignancy. Dysregulated RNA splicing factors have been widely reported in tumorigenesis, but their involvement in pancreatitis and PDAC is not well understood. Here, we report that the splicing factor SRSF1 is highly expressed in pancreatitis, PDAC precursor lesions, and tumors. Increased SRSF1 is sufficient to induce pancreatitis and accelerate KRASG12D-mediated PDAC. Mechanistically, SRSF1 activates MAPK signaling-partly by upregulating interleukin 1 receptor type 1 (IL1R1) through alternative-splicing-regulated mRNA stability. Additionally, SRSF1 protein is destabilized through a negative feedback mechanism in phenotypically normal epithelial cells expressing KRASG12D in mouse pancreas and in pancreas organoids acutely expressing KRASG12D, buffering MAPK signaling and maintaining pancreas cell homeostasis. This negative feedback regulation of SRSF1 is overcome by hyperactive MYC, facilitating PDAC tumorigenesis. Our findings implicate SRSF1 in the etiology of pancreatitis and PDAC, and point to SRSF1-misregulated alternative splicing as a potential therapeutic target.

SIGNIFICANCE

We describe the regulation of splicing factor SRSF1 expression in the context of pancreas cell identity, plasticity, and inflammation. SRSF1 protein downregulation is involved in a negative feedback cellular response to KRASG12D expression, contributing to pancreas cell homeostasis. Conversely, upregulated SRSF1 promotes pancreatitis and accelerates KRASG12D-mediated tumorigenesis through enhanced IL1 and MAPK signaling. This article is highlighted in the In This Issue feature, p. 1501.

Direct androgen receptor regulation of sexually dimorphic gene expression in the mammalian kidney

Mammalian organs exhibit distinct physiology, disease susceptibility and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA-seq data demonstrated sex differences were established from 4 and 8 weeks after birth under gonadal control. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR) mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation while analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, and disease and metabolic linkage, of sexually dimorphic gene activity.

RBFOX2 modulates a metastatic signature of alternative splicing in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDA) is characterized by aggressive local invasion and metastatic spread, leading to high lethality. Although driver gene mutations during PDA progression are conserved, no specific mutation is correlated with the dissemination of metastases1-3. Here we analysed RNA splicing data of a large cohort of primary and metastatic PDA tumours to identify differentially spliced events that correlate with PDA progression. De novo motif analysis of these events detected enrichment of motifs with high similarity to the RBFOX2 motif. Overexpression of RBFOX2 in a patient-derived xenograft (PDX) metastatic PDA cell line drastically reduced the metastatic potential of these cells in vitro and in vivo, whereas depletion of RBFOX2 in primary pancreatic tumour cell lines increased the metastatic potential of these cells. These findings support the role of RBFOX2 as a potent metastatic suppressor in PDA. RNA-sequencing and splicing analysis of RBFOX2 target genes revealed enrichment of genes in the RHO GTPase pathways, suggesting a role of RBFOX2 splicing activity in cytoskeletal organization and focal adhesion formation. Modulation of RBFOX2-regulated splicing events, such as via myosin phosphatase RHO-interacting protein (MPRIP), is associated with PDA metastases, altered cytoskeletal organization and the induction of focal adhesion formation. Our results implicate the splicing-regulatory function of RBFOX2 as a tumour suppressor in PDA and suggest a therapeutic approach for metastatic PDA.

Identification and validation of RNA methylation-related alternative splicing gene signature for low-grade glioma to predict survival and immune landscapes

BACKGROUND

Low-grade glioma (LGG) is a crucial pathological type of glioma. The present study aimed to explore multiple RNA methylation regulator-related AS events and investigate their prognostic values in LGG.

METHODS

The prognostic model for low-grade glioma was established using the LASSO regression analysis. To validate prognostic value, we performed Kaplan-Maier survival analysis, ROC curves and nomograms. The ESTIMATE algorithm, the CIBERSORT algorithm and the ssGSEA algorithm were utilized to explore the role of the immune microenvironment in LGG. Subsequently, we then used GO, KEGG and GSEA enrichment analysis to explore the functional roles of these genes. In addition, we employed the GDSC database to screen potential chemotherapeutic agents.

RESULTS

Eight RNA methylation related AS events were involved in construct a survival and prognosis model, which had good ability of independent prediction for patients with LGG. Patients in the high-risk group had shorter life expectancy and higher mortality, while patients in the low-risk group had a better prognosis. We constructed a nomogram which showed an excellent predictive performance for individual OS. The risk score exhibited a close correlation with some immune cells and expression of immune checkpoints. Patients in high-risk group were characterized by immunosuppressive microenvironment and poor response to immunotherapy, and were sensitive to more chemotherapeutic drugs. Pathway and functional enrichment analyses further confirmed that significant differences existed in immune landscape between the two subgroups.

CONCLUSION

The prognostic RNA methylation-related alternative splicing signature constructed could constitute a promising prognostic biomarker, which could serve to optimize treatment regimens.

Comparative full-length transcriptome analysis by Oxford Nanopore Technologies reveals genes involved in anthocyanin accumulation in storage roots of sweet potatoes (Ipomoea batatas L.)

Background

Storage roots of sweet potatoes (Ipomoea batatas L.) with different colors vary in anthocyanin content, indicating different economically agronomic trait. As the newest DNA/RNA sequencing technology, Oxford Nanopore Technologies (ONT) have been applied in rapid transcriptome sequencing for investigation of genes related to nutrient metabolism. At present, few reports concern full-length transcriptome analysis based on ONT for study on the molecular mechanism of anthocyanin accumulation leading to color change of tuberous roots of sweet potato cultivars.

Results

The storage roots of purple-fleshed sweet potato (PFSP) and white-fleshed sweet potato (WFSP) at different developmental stages were subjected to anthocyanin content comparison by UV-visible spectroscopy as well as transcriptome analysis at ONT MinION platform. UV-visible spectrophotometric measurements demonstrated the anthocyanin content of PFSP was much higher than that of WFSP. ONT RNA-Seq results showed each sample generated average 2.75 GB clean data with Full-Length Percentage (FL%) over 70% and the length of N50 ranged from 1,192 to 1,395 bp, indicating reliable data for transcriptome analysis. Subsequent analysis illustrated intron retention was the most prominent splicing event present in the resulting transcripts. As compared PFSP with WFSP at the relative developmental stages with the highest (PH vs. WH) and the lowest (PL vs. WL) anthocyanin content, 282 and 216 genes were up-regulated and two and 11 genes were down-regulated respectively. The differential expression genes involved in flavonoid biosynthesis pathway include CCoAOMT, PpLDOX, DFR, Cytochrome P450, CHI, and CHS. The genes encoding oxygenase superfamily were significantly up-regulated when compared PFSP with WFSP at the relative developmental stages.

Conclusions

Comparative full-length transcriptome analysis based on ONT serves as an effective approach to detect the differences in anthocyanin accumulation in the storage roots of different sweet potato cultivars at transcript level, with noting that some key genes can now be closely related to flavonoids biosynthesis. This study helps to improve understanding of molecular mechanism for anthocyanin accumulation in sweet potatoes and also provides a theoretical basis for high-quality sweet potato breeding.

Oncogenic ZMYND11-MBTD1 fusion protein anchors the NuA4/TIP60 histone acetyltransferase complex to the coding region of active genes

A recurrent chromosomal translocation found in acute myeloid leukemia leads to an in-frame fusion of the transcription repressor ZMYND11 to MBTD1, a subunit of the NuA4/TIP60 histone acetyltransferase complex. To understand the abnormal molecular events that ZMYND11-MBTD1 expression can create, we perform a biochemical and functional characterization comparison to each individual fusion partner. ZMYND11-MBTD1 is stably incorporated into the endogenous NuA4/TIP60 complex, leading to its mislocalization on the body of genes normally bound by ZMYND11. This can be correlated to increased chromatin acetylation and altered gene transcription, most notably on the MYC oncogene, and alternative splicing. Importantly, ZMYND11-MBTD1 expression favors Myc-driven pluripotency during embryonic stem cell differentiation and self-renewal of hematopoietic stem/progenitor cells. Altogether, these results indicate that the ZMYND11-MBTD1 fusion functions primarily by mistargeting the NuA4/TIP60 complex to the body of genes, altering normal transcription of specific genes, likely driving oncogenesis in part through the Myc regulatory network.

Identification of alternative splicing-derived cancer neoantigens for mRNA vaccine development

Messenger RNA (mRNA) vaccines have shown great potential for anti-tumor therapy due to the advantages in safety, efficacy and industrial production. However, it remains a challenge to identify suitable cancer neoantigens that can be targeted for mRNA vaccines. Abnormal alternative splicing occurs in a variety of tumors, which may result in the translation of abnormal transcripts into tumor-specific proteins. High-throughput technologies make it possible for systematic characterization of alternative splicing as a source of suitable target neoantigens for mRNA vaccine development. Here, we summarized difficulties and challenges for identifying alternative splicing-derived cancer neoantigens from RNA-seq data and proposed a conceptual framework for designing personalized mRNA vaccines based on alternative splicing-derived cancer neoantigens. In addition, several points were presented to spark further discussion toward improving the identification of alternative splicing-derived cancer neoantigens.

Alternative Splicing Regulation of Glycine-Rich Proteins via Target of Rapamycin-Reactive Oxygen Species Pathway in Arabidopsis Seedlings Upon Glucose Stress

Glucose can serve as both the source of energy and regulatory signaling molecule in plant. Due to the environmental and metabolic change, sugar levels could affect various developmental processes. High glucose environment is hardly conductive to the plant growth but cause development arrest. Increasing evidence indicate that alternative splicing (AS) plays a pivotal role in sugar signaling. However, the regulatory mechanism upon glucose stress remains unclear. The full-length transcriptomes were obtained from the samples of Arabidopsis seedlings with 3% glucose and mock treatment, using Oxford Nanopore sequencing technologies. Further analysis indicated that many genes involved in photosynthesis were significantly repressed and many genes involved in glycolysis, mitochondrial function, and the response to oxidative stress were activated. In total, 1,220 significantly differential alternative splicing (DAS) events related to 619 genes were identified, among which 75.74% belong to intron retention (IR). Notably, more than 20% of DAS events come from a large set of glycine-rich protein (GRP) family genes, such as GRP7, whose AS types mostly belong to IR. Besides the known productive GRP transcript isoforms, we identified a lot of splicing variants with diverse introns spliced in messenger RNA (mRNA) region coding the glycine-rich (GR) domain. The AS pattern of GRPs changed and particularly, the productive GRPs increased upon glucose stress. These ASs of GRP pre-mRNAs triggered by glucose stress could be abolished by AZD-8055, which is an ATP competitive inhibitor for the target of rapamycin (TOR) kinase but could be mimicked by H₂O₂. Additionally, AS pattern change of arginine/serine-rich splicing factor 31(RS31) via TOR pathway, which was previously described in response to light and sucrose signaling, was also induced in a similar manner by both glucose stress and reactive oxygen species (ROS). Here we conclude that (i) glucose stress suppresses photosynthesis and activates the glycolysis-mitochondria energy relay and ROS scavenging system; (ii) glucose stress triggers transcriptome-wide AS pattern changes including a large set of splicing factors, such as GRPs and RS31; (iii) high sugars regulate AS pattern change of both GRPs and RS31 via TOR-ROS pathway. The results from this study will deepen our understanding of the AS regulation mechanism in sugar signaling.

ASO-based PKM splice-switching therapy inhibits hepatocellular carcinoma growth

The M2 pyruvate kinase (PKM2) isoform is upregulated in most cancers and plays a crucial role in regulation of the Warburg effect, which is characterized by the preference for aerobic glycolysis over oxidative phosphorylation for energy metabolism. PKM2 is an alternative-splice isoform of the PKM gene and is a potential therapeutic target. Antisense oligonucleotides (ASO) that switch PKM splicing from the cancer-associated PKM2 to the PKM1 isoform have been shown to induce apoptosis in cultured glioblastoma cells when delivered by lipofection. Here, we explore the potential of ASO-based PKM splice switching as a targeted therapy for liver cancer. A more potent lead cEt/DNA ASO induced PKM splice-switching and inhibited the growth of cultured hepatocellular carcinoma (HCC) cells. This PKM isoform switch increased pyruvate-kinase activity and altered glucose metabolism. In an orthotopic HCC xenograft mouse model, the lead ASO and a second ASO targeting a non-overlapping site inhibited tumor growth. Finally, in a genetic HCC mouse model, a surrogate mouse-specific ASO induced Pkm splice switching and inhibited tumorigenesis, without observable toxicity. These results lay the groundwork for a potential ASO-based splicing therapy for HCC.

Genome-wide screening identifies cell cycle control as a synthetic lethal pathway with SRSF2P95H mutation

Current strategies to target RNA splicing mutant myeloid cancers proposes targeting the remaining splicing apparatus. This approach has only been modestly sensitizing and is also toxic to non-mutant bearing wild-type cells. To explore potentially exploitable genetic interactions with spliceosome mutations, we combined data mining and functional screening for synthetic lethal interactions with an Srsf2P95H/+ mutation. Analysis of mis-splicing events in a series of both human and murine SRSF2P95H mutant samples across multiple myeloid diseases (AML, MDS, CMML) was performed to identify conserved mis-splicing events. From this analysis, we identified that the cell cycle and DNA repair pathways were overrepresented within the conserved mis-spliced transcript sets. In parallel, to functionally define pathways essential for survival and proliferation of Srsf2P95H/+ cells, we performed a genome-wide CRISPR loss of function screen using Hoxb8 immortalized R26-CreERki/+ Srsf2P95H/+ and R26-CreERki/+ Srsf2+/+ cell lines. We assessed loss of sgRNA representation at three timepoints: immediately after Srsf2P95H/+ activation, and at one week and two weeks post Srsf2P95H/+ mutation. Pathway analysis demonstrated that the cell cycle and DNA damage response pathways were amongst the top synthetic lethal pathways with Srsf2P95H/+ mutation. Based on the loss of guide RNAs targeting Cdk6, we identified that Palbociclib, a CDK6 inhibitor, showed preferential sensitivity in Srsf2P95H/+ cell lines and in primary non-immortalized lin-cKIT+Sca-1+ cells compared to wild type controls. Our data strongly suggest that the cell cycle and DNA damage response pathways are required for Srsf2P95H/+ cell survival, and that Palbociclib could be an alternative therapeutic option for targeting SRSF2 mutant cancers.

Multi-omics of the esophageal microenvironment identifies signatures associated with progression of Barrett's esophagus

BACKGROUND

The enrichment of Gram-negative bacteria of oral origin in the esophageal microbiome has been associated with the development of metaplasia. However, to date, no study has comprehensively assessed the relationships between the esophageal microbiome and the host.

METHODS

Here, we examine the esophageal microenvironment in gastro-esophageal reflux disease and metaplasia using multi-omics strategies targeting the microbiome and host transcriptome, followed by targeted culture, comparative genomics, and host-microbial interaction studies of bacterial signatures of interest.

RESULTS

Profiling of the host transcriptome from esophageal mucosal biopsies revealed profound changes during metaplasia. Importantly, five biomarkers showed consistent longitudinal changes with disease progression from reflux disease to metaplasia. We showed for the first time that the esophageal microbiome is distinct from the salivary microbiome and the enrichment of Campylobacter species as a consistent signature in disease across two independent cohorts. Shape fitting and matrix correlation identified associations between the microbiome and host transcriptome profiles, with a novel co-exclusion relationship found between Campylobacter and napsin B aspartic peptidase. Targeted culture of Campylobacter species from the same cohort revealed a subset of isolates to have a higher capacity to survive within primary human macrophages. Comparative genomic analyses showed these isolates could be differentiated by specific genomic features, one of which was validated to be associated with intracellular fitness. Screening for these Campylobacter strain-specific signatures in shotgun metagenomics data from another cohort showed an increase in prevalence with disease progression. Comparative transcriptomic analyses of primary esophageal epithelial cells exposed to the Campylobacter isolates revealed expression changes within those infected with strains with high intracellular fitness that could explain the increased likelihood of disease progression.

CONCLUSIONS

We provide a comprehensive assessment of the esophageal microenvironment, identifying bacterial strain-specific signatures with high relevance to progression of metaplasia.

The small Cajal body-specific RNA 15 (SCARNA15) directs p53 and redox homeostasis via selective splicing in cancer cells

Small Cajal body-specific RNAs (scaRNAs) guide post-transcriptional modification of spliceosomal RNA and, while commonly altered in cancer, have poorly defined roles in tumorigenesis. Here, we uncover that SCARNA15 directs alternative splicing (AS) and stress adaptation in cancer cells. Specifically, we find that SCARNA15 guides critical pseudouridylation (Ψ) of U2 spliceosomal RNA to fine-tune AS of distinct transcripts enriched for chromatin and transcriptional regulators in malignant cells. This critically impacts the expression and function of the key tumor suppressors ATRX and p53. Significantly, SCARNA15 loss impairs p53-mediated redox homeostasis and hampers cancer cell survival, motility and anchorage-independent growth. In sum, these findings highlight an unanticipated role for SCARNA15 and Ψ in directing cancer-associated splicing programs.

N6-Methyladenosine RNA Methylation Regulator-Related Alternative Splicing (AS) Gene Signature Predicts Non-Small Cell Lung Cancer Prognosis

Aberrant N6-methyladenosine (m6A) RNA methylation regulatory genes and related gene alternative splicing (AS) could be used to predict the prognosis of non-small cell lung carcinoma. This study focused on 13 m6A regulatory genes (METTL3, METTL14, WTAP, KIAA1429, RBM15, ZC3H13, YTHDC1, YTHDC2, YTHDF1, YTHDF2, HNRNPC, FTO, and ALKBH5) and expression profiles in TCGA-LUAD (n = 504) and TCGA-LUSC (n = 479) datasets from the Cancer Genome Atlas database. The data were downloaded and bioinformatically and statistically analyzed, including the gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. There were 43,948 mRNA splicing events in lung adenocarcinoma (LUAD) and 46,020 in lung squamous cell carcinoma (LUSC), and the data suggested that m6A regulators could regulate mRNA splicing. Differential HNRNPC and RBM15 expression was associated with overall survival (OS) of LUAD and HNRNPC and METTL3 expression with the OS of LUSC patients. Furthermore, the non-small cell lung cancer prognosis-related AS events signature was constructed and divided patients into high- vs. low-risk groups using seven and 14 AS genes in LUAD and LUSC, respectively. The LUAD risk signature was associated with gender and T, N, and TNM stages, but the LUSC risk signature was not associated with any clinical features. In addition, the risk signature and TNM stage were independent prognostic predictors in LUAD and the risk signature and T stage were independent prognostic predictors in LUSC after the multivariate Cox regression and receiver operating characteristic analyses. In conclusion, this study revealed the AS prognostic signature in the prediction of LUAD and LUSC prognosis.

An interactive network of alternative splicing events with prognostic value in geriatric lung adenocarcinoma via the regulation of splicing factors

Alternative splicing (AS), an essential process for the maturation of mRNAs, is involved in tumorigenesis and tumor progression, including angiogenesis, apoptosis, and metastasis. AS changes can be frequently observed in different tumors, especially in geriatric lung adenocarcinoma (GLAD). Previous studies have reported an association between AS events and tumorigenesis but have lacked a systematic analysis of its underlying mechanisms. In the present study, we obtained splicing event information from SpliceSeq and clinical information regarding GLAD from The Cancer Genome Atlas. Survival-associated AS events were selected to construct eight prognostic index (PI) models. We also constructed a correlation network between splicing factors (SFs) and survival-related AS events to identify a potential molecular mechanism involved in regulating AS-related events in GLAD. Our study findings confirm that AS has a strong prognostic value for GLAD and sheds light on the clinical significance of targeting SFs in the treatment of GLAD.

Recurrent SRSF2 mutations in MDS affect both splicing and NMD

Oncogenic mutations in the RNA splicing factors SRSF2, SF3B1, and U2AF1 are the most frequent class of mutations in myelodysplastic syndromes and are also common in clonal hematopoiesis, acute myeloid leukemia, chronic lymphocytic leukemia, and a variety of solid tumors. They cause genome-wide splicing alterations that affect important regulators of hematopoiesis. Several mRNA isoforms promoted by the various splicing factor mutants comprise a premature termination codon (PTC) and are therefore potential targets of nonsense-mediated mRNA decay (NMD). In light of the mechanistic relationship between splicing and NMD, we sought evidence for a specific role of mutant SRSF2 in NMD. We show that SRSF2 Pro95 hot spot mutations elicit enhanced mRNA decay, which is dependent on sequence-specific RNA binding and splicing. SRSF2 mutants enhance the deposition of exon junction complexes (EJCs) downstream from the PTC through RNA-mediated molecular interactions. This architecture then favors the association of key NMD factors to elicit mRNA decay. Gene-specific blocking of EJC deposition by antisense oligonucleotides circumvents aberrant NMD promoted by mutant SRSF2, restoring the expression of PTC-containing transcript. Our study uncovered critical effects of SRSF2 mutants in hematologic malignancies, reflecting the regulation at multiple levels of RNA metabolism, from splicing to decay.

Altered RNA Processing in Cancer Pathogenesis and Therapy

Major advances in our understanding of cancer pathogenesis and therapy have come from efforts to catalog genomic alterations in cancer. A growing number of large-scale genomic studies have uncovered mutations that drive cancer by perturbing cotranscriptional and post-transcriptional regulation of gene expression. These include alterations that affect each phase of RNA processing, including splicing, transport, editing, and decay of messenger RNA. The discovery of these events illuminates a number of novel therapeutic vulnerabilities generated by aberrant RNA processing in cancer, several of which have progressed to clinical development. SIGNIFICANCE: There is increased recognition that genetic alterations affecting RNA splicing and polyadenylation are common in cancer and may generate novel therapeutic opportunities. Such mutations may occur within an individual gene or in RNA processing factors themselves, thereby influencing splicing of many downstream target genes. This review discusses the biological impact of these mutations on tumorigenesis and the therapeutic approaches targeting cells bearing these mutations.

Coordinated alterations in RNA splicing and epigenetic regulation drive leukaemogenesis

Transcription and pre-mRNA splicing are key steps in the control of gene expression and mutations in genes regulating each of these processes are common in leukaemia1,2. Despite the frequent overlap of mutations affecting epigenetic regulation and splicing in leukaemia, how these processes influence one another to promote leukaemogenesis is not understood and, to our knowledge, there is no functional evidence that mutations in RNA splicing factors initiate leukaemia. Here, through analyses of transcriptomes from 982 patients with acute myeloid leukaemia, we identified frequent overlap of mutations in IDH2 and SRSF2 that together promote leukaemogenesis through coordinated effects on the epigenome and RNA splicing. Whereas mutations in either IDH2 or SRSF2 imparted distinct splicing changes, co-expression of mutant IDH2 altered the splicing effects of mutant SRSF2 and resulted in more profound splicing changes than either mutation alone. Consistent with this, co-expression of mutant IDH2 and SRSF2 resulted in lethal myelodysplasia with proliferative features in vivo and enhanced self-renewal in a manner not observed with either mutation alone. IDH2 and SRSF2 double-mutant cells exhibited aberrant splicing and reduced expression of INTS3, a member of the integrator complex3, concordant with increased stalling of RNA polymerase II (RNAPII). Aberrant INTS3 splicing contributed to leukaemogenesis in concert with mutant IDH2 and was dependent on mutant SRSF2 binding to cis elements in INTS3 mRNA and increased DNA methylation of INTS3. These data identify a pathogenic crosstalk between altered epigenetic state and splicing in a subset of leukaemias, provide functional evidence that mutations in splicing factors drive myeloid malignancy development, and identify spliceosomal changes as a mediator of IDH2-mutant leukaemogenesis.