Uncovering the impacts of alternative splicing on the proteome with current omics techniques

Unveiling the veil of RNA binding protein phase separation in cancer biology and therapy

RNA-binding protein (RBP) phase separation in oncology reveals a complex interplay crucial for understanding tumor biology and developing novel therapeutic strategies. Aberrant phase separation of RBPs significantly influences gene regulation, signal transduction, and metabolic reprogramming, contributing to tumorigenesis and drug resistance. Our review highlights the integral roles of RBP phase separation in stress granule dynamics, mRNA stabilization, and the modulation of transcriptional and translational processes. Furthermore, interactions between RBPs and non-coding RNAs add a layer of complexity, providing new insights into their collaborative roles in cancer progression. The intricate relationship between RBPs and phase separation poses significant challenges but also opens up novel opportunities for targeted therapeutic interventions. Advancing our understanding of the molecular mechanisms and regulatory networks governing RBP phase separation could lead to breakthroughs in cancer treatment strategies.

Relevance and regulation of alternative splicing in plant secondary metabolism: current understanding and future directions

The secondary metabolism of plants is an essential life process enabling organisms to navigate various stages of plant development and cope with ever-changing environmental stresses. Secondary metabolites, abundantly found in nature, possess significant medicinal value. Among the regulatory mechanisms governing these metabolic processes, alternative splicing stands out as a widely observed post-transcriptional mechanism present in multicellular organisms. It facilitates the generation of multiple mRNA transcripts from a single gene by selecting different splicing sites. Selective splicing events in plants are widely induced by various signals, including external environmental stress and hormone signals. These events ultimately regulate the secondary metabolic processes and the accumulation of essential secondary metabolites in plants by influencing the synthesis of primary metabolites, hormone metabolism, biomass accumulation, and capillary density. Simultaneously, alternative splicing plays a crucial role in enhancing protein diversity and the abundance of the transcriptome. This paper provides a summary of the factors inducing alternative splicing events in plants and systematically describes the progress in regulating alternative splicing with respect to different secondary metabolites, including terpenoid, phenolic compounds, and nitrogen-containing compounds. Such elucidation offers critical foundational insights for understanding the role of alternative splicing in regulating plant metabolism and presents novel avenues and perspectives for bioengineering.

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

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

IS-PRM-Based Peptide Targeting Informed by Long-Read Sequencing for Alternative Proteome Detection

Alternative splicing is a major contributor of transcriptomic complexity, but the extent to which transcript isoforms are translated into stable, functional protein isoforms is unclear. Furthermore, detection of relatively scarce isoform-specific peptides is challenging, with many protein isoforms remaining uncharted due to technical limitations. Recently, a family of advanced targeted MS strategies, termed internal standard parallel reaction monitoring (IS-PRM), have demonstrated multiplexed, sensitive detection of predefined peptides of interest. Such approaches have not yet been used to confirm existence of novel peptides. Here, we present a targeted proteogenomic approach that leverages sample-matched long-read RNA sequencing (lrRNA-seq) data to predict potential protein isoforms with prior transcript evidence. Predicted tryptic isoform-specific peptides, which are specific to individual gene product isoforms, serve as "triggers" and "targets" in the IS-PRM method, Tomahto. Using the model human stem cell line WTC11, LR RNaseq data were generated and used to inform the generation of synthetic standards for 192 isoform-specific peptides (114 isoforms from 55 genes). These synthetic "trigger" peptides were labeled with super heavy tandem mass tags (TMT) and spiked into TMT-labeled WTC11 tryptic digest, predicted to contain corresponding endogenous "target" peptides. Compared to DDA mode, Tomahto increased detectability of isoforms by 3.6-fold, resulting in the identification of five previously unannotated isoforms. Our method detected protein isoform expression for 43 out of 55 genes corresponding to 54 resolved isoforms. This lrRNA-seq-informed Tomahto targeted approach is a new modality for generating protein-level evidence of alternative isoforms─a critical first step in designing functional studies and eventually clinical assays.

RBM10 loss induces aberrant splicing of cytoskeletal and extracellular matrix mRNAs and promotes metastatic fitness

RBM10 modulates transcriptome-wide cassette exon splicing. Loss-of-function RBM10 mutations are enriched in thyroid cancers with distant metastases. Analysis of transcriptomes and genes mis-spliced by RBM10 loss showed pro-migratory and RHO/RAC signaling signatures. RBM10 loss increases cell velocity. Cytoskeletal and ECM transcripts subject to exon-inclusion events included vinculin (VCL), tenascin C (TNC) and CD44. Knockdown of the VCL exon inclusion transcript in RBM10-null cells reduced cell velocity, whereas knockdown of TNC and CD44 exon-inclusion isoforms reduced invasiveness. RAC1-GTP levels were increased in RBM10-null cells. Mouse Hras G12V /Rbm1O KO thyrocytes develop metastases that are reversed by RBM10 or by combined knockdown of VCL, CD44 and TNC inclusion isoforms. Thus, RBM10 loss generates exon inclusions in transcripts regulating ECM-cytoskeletal interactions, leading to RAC1 activation and metastatic competency. Moreover, a CRISPR-Cas9 screen for synthetic lethality with RBM10 loss identified NFkB effectors as central to viability, providing a therapeutic target for these lethal thyroid cancers.

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

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

Population genetics and phylogeographic history of the insular lizard Podarcis lilfordi (Gunther, 1874) from the Balearic Islands based on genome-wide polymorphic data

Islands provide a great system to explore the processes that maintain genetic diversity and promote local adaptation. We explored the genomic diversity of the Balearic lizard Podarcis lilfordi, an endemic species characterized by numerous small insular populations with large phenotypic diversity. Using the newly available genome for this species, we characterized more than 300,000 SNPs, merging genotyping-by-sequencing (GBS) data with previously published restriction site-associated DNA sequencing (RAD-Seq) data, providing a dataset of 16 island populations (191 individuals) across the range of species distribution (Menorca, Mallorca, and Cabrera). Results indicate that each islet hosts a well-differentiated population (F ST = 0.247 ± 0.09), with no recent immigration/translocation events. Contrary to expectations, most populations harbor a considerable genetic diversity (mean nucleotide diversity, P i = 0.144 ± 0.021), characterized by overall low inbreeding values (F IS < 0.1). While the genetic diversity significantly decreased with decreasing islet surface, maintenance of substantial genetic diversity even in tiny islets suggests variable selection or other mechanisms that buffer genetic drift. Maximum-likelihood tree based on concatenated SNP data confirmed the existence of the two major independent lineages of Menorca and Mallorca/Cabrera. Multiple lines of evidence, including admixture and root testing, robustly placed the origin of the species in the Mallorca Island, rather than in Menorca. Outlier analysis mainly retrieved a strong signature of genome differentiation between the two major archipelagos, especially in the sexual chromosome Z. A set of proteins were target of multiple outliers and primarily associated with binding and catalytic activity, providing interesting candidates for future selection studies. This study provides the framework to explore crucial aspects of the genetic basis of phenotypic divergence and insular adaptation.

Updating mRNA variants of the human RSK4 gene and their expression in different stressed situations

We determined RNA spectrum of the human RSK4 (hRSK4) gene (also called RPS6KA6) and identified 29 novel mRNA variants derived from alternative splicing, which, plus the NCBI-documented ones and the five we reported previously, totaled 50 hRSK4 RNAs that, by our bioinformatics analyses, encode 35 hRSK4 protein isoforms of 35-762 amino acids. Many of the mRNAs are bicistronic or tricistronic for hRSK4. The NCBI-normalized NM_014496.5 and the protein it encodes are designated herein as the Wt-1 mRNA and protein, respectively, whereas the NM_001330512.1 and the long protein it encodes are designated as the Wt-2 mRNA and protein, respectively. Many of the mRNA variants responded differently to different situations of stress, including serum starvation, a febrile temperature, treatment with ethanol or ethanol-extracted clove buds (an herbal medicine), whereas the same stressed situation often caused quite different alterations among different mRNA variants in different cell lines. Mosifloxacin, an antibiotics and also a functional inhibitor of hRSK4, could inhibit the expression of certain hRSK4 mRNA variants. The hRSK4 gene likely uses alternative splicing as a handy tool to adapt to different stressed situations, and the mRNA and protein multiplicities may partly explain the incongruous literature on its expression and comports.

NONO promotes gallbladder cancer cell proliferation by enhancing oncogenic RNA splicing of DLG1 through interaction with IGF2BP3/RBM14

Gallbladder cancer (GBC) is a highly malignant and rapidly progressing tumor of the human biliary system, and there is an urgent need to develop new therapeutic targets and modalities. Non-POU domain-containing octamer-binding protein (NONO) is an RNA-binding protein involved in the regulation of transcription, mRNA splicing, and DNA repair. NONO expression is elevated in multiple tumors and can act as an oncogene to promote tumor progression. Here, we found that NONO was highly expressed in GBC and promoted tumor cells growth. The dysregulation of RNA splicing is a molecular feature of almost all tumor types. Accordingly, mRNA-seq and RIP-seq analysis showed that NONO promoted exon6 skipping in DLG1, forming two isomers (DLG1-FL and DLG1-S). Furthermore, lower Percent-Spliced-In (PSI) values of DLG1 were detected in tumor tissue relative to the paraneoplastic tissue, and were associated with poor patient prognosis. Moreover, DLG1-S and DLG1-FL act as tumor promoters and tumor suppressors, respectively, by regulating the YAP1/JUN pathway. N6-methyladenosine (m6A) is the most common and abundant RNA modification involved in alternative splicing processes. We identified an m6A reader, IGF2BP3, which synergizes with NONO to promote exon6 skipping in DLG1 in an m6A-dependent manner. Furthermore, IP/MS results showed that RBM14 was bound to NONO and interfered with NONO-mediated exon6 skipping of DLG1. In addition, IGF2BP3 disrupted the binding of RBM14 to NONO. Overall, our data elucidate the molecular mechanism by which NONO promotes DLG1 exon skipping, providing a basis for new therapeutic targets in GBC treatment.

IS-PRM-based peptide targeting informed by long-read sequencing for alternative proteome detection

Alternative splicing is a major contributor of transcriptomic complexity, but the extent to which transcript isoforms are translated into stable, functional protein isoforms is unclear. Furthermore, detection of relatively scarce isoform-specific peptides is challenging, with many protein isoforms remaining uncharted due to technical limitations. Recently, a family of advanced targeted MS strategies, termed internal standard parallel reaction monitoring (IS-PRM), have demonstrated multiplexed, sensitive detection of pre-defined peptides of interest. Such approaches have not yet been used to confirm existence of novel peptides. Here, we present a targeted proteogenomic approach that leverages sample-matched long-read RNA sequencing (LR RNAseq) data to predict potential protein isoforms with prior transcript evidence. Predicted tryptic isoform-specific peptides, which are specific to individual gene product isoforms, serve as "triggers" and "targets" in the IS-PRM method, Tomahto. Using the model human stem cell line WTC11, LR RNAseq data were generated and used to inform the generation of synthetic standards for 192 isoform-specific peptides (114 isoforms from 55 genes). These synthetic "trigger" peptides were labeled with super heavy tandem mass tags (TMT) and spiked into TMT-labeled WTC11 tryptic digest, predicted to contain corresponding endogenous "target" peptides. Compared to DDA mode, Tomahto increased detectability of isoforms by 3.6-fold, resulting in the identification of five previously unannotated isoforms. Our method detected protein isoform expression for 43 out of 55 genes corresponding to 54 resolved isoforms. This LR RNA seq-informed Tomahto targeted approach, called LRP-IS-PRM, is a new modality for generating protein-level evidence of alternative isoforms - a critical first step in designing functional studies and eventually clinical assays.

Biosurfer for systematic tracking of regulatory mechanisms leading to protein isoform diversity

Long-read RNA sequencing has shed light on transcriptomic complexity, but questions remain about the functionality of downstream protein products. We introduce Biosurfer, a computational approach for comparing protein isoforms, while systematically tracking the transcriptional, splicing, and translational variations that underlie differences in the sequences of the protein products. Using Biosurfer, we analyzed the differences in 32,799 pairs of GENCODE annotated protein isoforms, finding a majority (70%) of variable N-termini are due to the alternative transcription start sites, while only 9% arise from 5' UTR alternative splicing. Biosurfer's detailed tracking of nucleotide-to-residue relationships helped reveal an uncommonly tracked source of single amino acid residue changes arising from the codon splits at junctions. For 17% of internal sequence changes, such split codon patterns lead to single residue differences, termed "ragged codons". Of variable C-termini, 72% involve splice- or intron retention-induced reading frameshifts. We found an unusual pattern of reading frame changes, in which the first frameshift is closely followed by a distinct second frameshift that restores the original frame, which we term a "snapback" frameshift. We analyzed long read RNA-seq-predicted proteome of a human cell line and found similar trends as compared to our GENCODE analysis, with the exception of a higher proportion of isoforms predicted to undergo nonsense-mediated decay. Biosurfer's comprehensive characterization of long-read RNA-seq datasets should accelerate insights of the functional role of protein isoforms, providing mechanistic explanation of the origins of the proteomic diversity driven by the alternative splicing. Biosurfer is available as a Python package at https://github.com/sheynkman-lab/biosurfer.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity

BACKGROUND

Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms.

RESULTS

In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms.

CONCLUSION

This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

Mutation of two intronic nucleotides alters RNA structure and dynamics inhibiting MBNL1 and RBFOX1 regulated splicing of the Insulin Receptor

Alternative splicing (AS) of Exon 11 of the Insulin Receptor ( INSR ) is highly regulated and disrupted in several human disorders. To better understand INSR exon 11 AS regulation, splicing activity of an INSR exon 11 minigene reporter was measured across a gradient of the AS regulator muscleblind-like 1 protein (MBNL1). The RNA-binding protein Fox-1 (RBFOX1) was added to determine its impact on MBNL1-regulated splicing. The role of the RBFOX1 UGCAUG binding site within intron 11 was assessed across the MBNL1 gradient. Mutating the UGCAUG motif inhibited RBFOX1 regulation of exon 11 and had the unexpected effect of reducing MBNL1 regulation of this exon. Molecular dynamics simulations showed that exon 11 and the adjacent RNA adopts a dynamically stable conformation. Mutation of the RBFOX1 binding site altered RNA structure and dynamics, while a mutation that created an optimal MBNL1 binding site at the RBFOX1 site shifted the RNA back to wild type. An antisense oligonucleotide (ASO) was used to confirm the structure in this region of the pre-mRNA. This example of intronic mutations shifting pre-mRNA structure and dynamics to modulate splicing suggests RNA structure and dynamics should be taken into consideration for AS regulation and therapeutic interventions targeting pre-mRNA.

Global detection of human variants and isoforms by deep proteome sequencing

An average shotgun proteomics experiment detects approximately 10,000 human proteins from a single sample. However, individual proteins are typically identified by peptide sequences representing a small fraction of their total amino acids. Hence, an average shotgun experiment fails to distinguish different protein variants and isoforms. Deeper proteome sequencing is therefore required for the global discovery of protein isoforms. Using six different human cell lines, six proteases, deep fractionation and three tandem mass spectrometry fragmentation methods, we identify a million unique peptides from 17,717 protein groups, with a median sequence coverage of approximately 80%. Direct comparison with RNA expression data provides evidence for the translation of most nonsynonymous variants. We have also hypothesized that undetected variants likely arise from mutation-induced protein instability. We further observe comparable detection rates for exon-exon junction peptides representing constitutive and alternative splicing events. Our dataset represents a resource for proteoform discovery and provides direct evidence that most frame-preserving alternatively spliced isoforms are translated.

Analyzing alternative splicing in Alzheimer's disease postmortem brain: a cell-level perspective

Alzheimer's disease (AD) is a neurodegenerative disease with no effective cure that attacks the brain's cells resulting in memory loss and changes in behavior and language skills. Alternative splicing is a highly regulated process influenced by specific cell types and has been implicated in age-related disorders such as neurodegenerative diseases. A comprehensive detection of alternative splicing events (ASEs) at the cellular level in postmortem brain tissue can provide valuable insights into AD pathology. Here, we provided cell-level ASEs in postmortem brain tissue by employing bioinformatics pipelines on a bulk RNA sequencing study sorted by cell types and two single-cell RNA sequencing studies from the prefrontal cortex. This comprehensive analysis revealed previously overlooked splicing and expression changes in AD patient brains. Among the observed alterations were changed in the splicing and expression of transcripts associated with chaperones, including CLU in astrocytes and excitatory neurons, PTGDS in astrocytes and endothelial cells, and HSP90AA1 in microglia and tauopathy-afflicted neurons, which were associated with differential expression of the splicing factor DDX5. In addition, novel, unknown transcripts were altered, and structural changes were observed in lncRNAs such as MEG3 in neurons. This work provides a novel strategy to identify the notable ASEs at the cell level in neurodegeneration, which revealed cell type-specific splicing changes in AD. This finding may contribute to interpreting associations between splicing and neurodegenerative disease outcomes.

Cas13d-mediated isoform-specific RNA knockdown with a unified computational and experimental toolbox

Alternative splicing is an essential mechanism for diversifying proteins, in which mature RNA isoforms produce proteins with potentially distinct functions. Two major challenges in characterizing the cellular function of isoforms are the lack of experimental methods to specifically and efficiently modulate isoform expression and computational tools for complex experimental design. To address these gaps, we developed and methodically tested a strategy which pairs the RNA-targeting CRISPR/Cas13d system with guide RNAs that span exon-exon junctions in the mature RNA. We performed a high-throughput essentiality screen, quantitative RT-PCR assays, and PacBio long read sequencing to affirm our ability to specifically target and robustly knockdown individual RNA isoforms. In parallel, we provide computational tools for experimental design and screen analysis. Considering all possible splice junctions annotated in GENCODE for multi-isoform genes and our gRNA efficacy predictions, we estimate that our junction-centric strategy can uniquely target up to 89% of human RNA isoforms, including 50,066 protein-coding and 11,415 lncRNA isoforms. Importantly, this specificity spans all splicing and transcriptional events, including exon skipping and inclusion, alternative 5' and 3' splice sites, and alternative starts and ends.

Inhibition of splicing factors SF3A3 and SRSF5 contributes to As3+/Se4+ combination-mediated proliferation suppression and apoptosis induction in acute promyelocytic leukemia cells

The low-dose combination of Arsenite (As³⁺) and selenite (Se⁴⁺) has the advantages of lower biological toxicity and better curative effects for acute promyelocytic leukemia (APL) therapy. However, the underlying mechanisms remain unclear. Here, based on the fact that the combination of 2 μM A³⁺ plus 4 μM Se⁴⁺ possessed a stronger anti-leukemic effect on APL cell line NB4 as compared with each individual, we employed iTRAQ-based quantitative proteomics to identify a total of 58 proteins that were differentially expressed after treatment with As³⁺/Se⁴⁺ combination rather than As³⁺ or Se⁴⁺ alone, the majority of which were involved in spliceosome pathway. Among them, eight proteins stood out by virtue of their splicing function and significant changes. They were validated as being decreased in mRNA and protein levels under As³⁺/Se⁴⁺ combination treatment. Further functional studies showed that only knockdown of two splicing factors, SF3A3 and SRSF5, suppressed the growth of NB4 cells. The reduction of SF3A3 was found to cause G1/S cell cycle arrest, which resulted in proliferation inhibition. Moreover, SRSF5 downregulation induced cell apoptosis through the activation of caspase-3. Taken together, these findings indicate that SF3A3 and SRSF5 function as pro-leukemic factors and can be potential novel therapeutic targets for APL.

DELongSeq for efficient detection of differential isoform expression from long-read RNA-seq data

Conventional gene expression quantification approaches, such as microarrays or quantitative PCR, have similar variations of estimates for all genes. However, next-generation short-read or long-read sequencing use read counts to estimate expression levels with much wider dynamic ranges. In addition to the accuracy of estimated isoform expression, efficiency, which measures the degree of estimation uncertainty, is also an important factor for downstream analysis. Instead of read count, we present DELongSeq, which employs information matrix of EM algorithm to quantify uncertainty of isoform expression estimates to improve estimation efficiency. DELongSeq uses random-effect regression model for the analysis of DE isoform, in that within-study variation represents variable precision in isoform expression estimation and between-study variation represents variation in isoform expression levels across samples. More importantly, DELongSeq allows 1 case versus 1 control comparison of differential expression, which has specific application scenarios in precision medicine (such as before versus after treatment, or tumor versus stromal tissues). Through extensive simulations and analysis of several RNA-Seq datasets, we show that the uncertainty quantification approach is computationally reliable, and can improve the power of differential expression (DE) analysis of isoforms or genes. In summary, DELongSeq allows for efficient detection of differential isoform/gene expression from long-read RNA-Seq data.

RNA splicing dysregulation and the hallmarks of cancer

Dysregulated RNA splicing is a molecular feature that characterizes almost all tumour types. Cancer-associated splicing alterations arise from both recurrent mutations and altered expression of trans-acting factors governing splicing catalysis and regulation. Cancer-associated splicing dysregulation can promote tumorigenesis via diverse mechanisms, contributing to increased cell proliferation, decreased apoptosis, enhanced migration and metastatic potential, resistance to chemotherapy and evasion of immune surveillance. Recent studies have identified specific cancer-associated isoforms that play critical roles in cancer cell transformation and growth and demonstrated the therapeutic benefits of correcting or otherwise antagonizing such cancer-associated mRNA isoforms. Clinical-grade small molecules that modulate or inhibit RNA splicing have similarly been developed as promising anticancer therapeutics. Here, we review splicing alterations characteristic of cancer cell transcriptomes, dysregulated splicing's contributions to tumour initiation and progression, and existing and emerging approaches for targeting splicing for cancer therapy. Finally, we discuss the outstanding questions and challenges that must be addressed to translate these findings into the clinic.

Aberrant splicing in human cancer: An RNA structural code point of view

Alternative splicing represents an essential process that occurs widely in eukaryotes. In humans, most genes undergo alternative splicing to ensure transcriptome and proteome diversity reflecting their functional complexity. Over the last decade, aberrantly spliced transcripts due to mutations in cis- or trans-acting splicing regulators have been tightly associated with cancer development, largely drawing scientific attention. Although a plethora of single proteins, ribonucleoproteins, complexed RNAs, and short RNA sequences have emerged as nodal contributors to the splicing cascade, the role of RNA secondary structures in warranting splicing fidelity has been underestimated. Recent studies have leveraged the establishment of novel high-throughput methodologies and bioinformatic tools to shed light on an additional layer of splicing regulation in the context of RNA structural elements. This short review focuses on the most recent available data on splicing mechanism regulation on the basis of RNA secondary structure, emphasizing the importance of the complex RNA G-quadruplex structures (rG4s), and other specific RNA motifs identified as splicing silencers or enhancers. Moreover, it intends to provide knowledge on newly established techniques that allow the identification of RNA structural elements and highlight the potential to develop new RNA-oriented therapeutic strategies against cancer.

Oxidative stress alters transcript localization of disease-associated genes in the retinal pigment epithelium

PURPOSE

Nuclear retention is a mechanism whereby excess RNA transcripts are stored in the event that a cell needs to quickly respond to a stimulus; maintaining proper nuclear-to-cytoplasmic balance is important for cellular homeostasis and cell function. There are many mechanisms that are employed to determine whether to retain a transcript or export it to the cytoplasm, although the extent to which tissue or cell type, internal and external stressors, and disease pathogenesis affect this process is not yet clear. As the most biochemically active tissue in the body, the retina must mitigate endogenous and exogenous stressors to maintain cell health and tissue function. Oxidative stress, believed to contribute to the pathogenesis or progression of age-related macular degeneration (AMD) and inherited retinal dystrophies (IRDs), is produced both internally from biochemical processes as well as externally from environmental insult. Here, we evaluate the effect of oxidative stress on transcript localization in the retinal pigment epithelium (RPE), with specific focus on transcripts related to RPE function and disease.

METHODS

We performed poly(A) RNA sequencing on nuclear and cytoplasmic fractions from human induced pluripotent stem cell-derived retinal pigment epithelium (iPSC-RPE) cells exposed to hydrogen peroxide (H₂O₂), as well as on untreated controls.

RESULTS

Under normal conditions, the number of mRNA transcripts retained in the nucleus exceeded that found in studies on other tissues. Further, the nuclear-to-cytoplasmic ratio of transcripts was altered following oxidative stress, as was the retention of genes associated with AMD and IRDs, as well as those that are important for RPE physiology.

CONCLUSIONS

These results provide a localization catalog of all expressed mRNA in iPSC-RPE under normal conditions and after exposure to H₂O₂, shedding light on the extent to which H₂O₂ alters transcript localization and potentially offering insight into one mechanism through which oxidative stress may contribute to the progression of visual disorders.

APPRIS principal isoforms and MANE Select transcripts define reference splice variants

MOTIVATION

Selecting the splice variant that best represents a coding gene is a crucial first step in many experimental analyses, and vital for mapping clinically relevant variants. This study compares the longest isoforms, MANE Select transcripts, APPRIS principal isoforms, and expression data, and aims to determine which method is best for selecting biological important reference splice variants for large-scale analyses.

RESULTS

Proteomics analyses and human genetic variation data suggest that most coding genes have a single main protein isoform. We show that APPRIS principal isoforms and MANE Select transcripts best describe these main cellular isoforms, and find that using the longest splice variant as the representative is a poor strategy. Exons unique to the longest splice isoforms are not under selective pressure, and so are unlikely to be functionally relevant. Expression data are also a poor means of selecting the main splice variant. APPRIS principal and MANE Select exons are under purifying selection, while exons specific to alternative transcripts are not. There are MANE and APPRIS representatives for almost 95% of genes, and where they agree they are particularly effective, coinciding with the main proteomics isoform for over 98.2% of genes.

AVAILABILITY AND IMPLEMENTATION

APPRIS principal isoforms for human, mouse and other model species can be downloaded from the APPRIS database (https://appris.bioinfo.cnio.es), GENCODE genes (https://www.gencodegenes.org/) and the Ensembl website (https://www.ensembl.org). MANE Select transcripts for the human reference set are available from the Ensembl, GENCODE and RefSeq databases (https://www.ncbi.nlm.nih.gov/refseq/). Lists of splice variants where MANE and APPRIS coincide are available from the APPRIS database.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

The evolution and structure/function of bHLH-PAS transcription factor family

Proteins that contain basic helix-loop-helix (bHLH) and Per-Arnt-Sim motifs (PAS) function as transcription factors. bHLH-PAS proteins exhibit essential and diverse functions throughout the body, from cell specification and differentiation in embryonic development to the proper function of organs like the brain and liver in adulthood. bHLH-PAS proteins are divided into two classes, which form heterodimers to regulate transcription. Class I bHLH-PAS proteins are typically activated in response to specific stimuli, while class II proteins are expressed more ubiquitously. Here, we discuss the general structure and functions of bHLH-PAS proteins throughout the animal kingdom, including family members that do not fit neatly into the class I-class II organization. We review heterodimerization between class I and class II bHLH-PAS proteins, binding partner selectivity and functional redundancy. Finally, we discuss the evolution of bHLH-PAS proteins, and why a class I protein essential for cardiovascular development in vertebrates like chicken and fish is absent from mammals.

Transcriptional Basis of Psoriasis from Large Scale Gene Expression Studies: The Importance of Moving towards a Precision Medicine Approach

Transcriptome profiling techniques, such as microarrays and RNA sequencing (RNA-seq), are valuable tools for deciphering the regulatory network underlying psoriasis and have revealed large number of differentially expressed genes in lesional and non-lesional skin. Such approaches provide a more precise measurement of transcript levels and their isoforms than any other methods. Large cohort transcriptomic analyses have greatly improved our understanding of the physiological and molecular mechanisms underlying disease pathogenesis and progression. Here, we mostly review the findings of some important large scale psoriatic transcriptomic studies, and the benefits of such studies in elucidating potential therapeutic targets and biomarkers for psoriasis treatment. We also emphasised the importance of looking into the alternatively spliced RNA isoforms/transcripts in psoriasis, rather than focussing only on the gene-level annotation. The neutrophil and blood transcriptome signature in psoriasis is also briefly reviewed, as it provides the immune status information of patients and is a less invasive platform. The application of precision medicine in current management of psoriasis, by combining transcriptomic data, improves the clinical response outcome in individual patients. Drugs tailored to individual patient’s genetic profile will greatly improve patient outcome and cost savings for the healthcare system.