Alternative splicing: Human disease and quantitative analysis from high-throughput sequencing

Incidence of alternative splicing associated with sex and opioid effects in the axon guidance pathway

The alternative splicing of a gene results in distinct transcript isoforms that can result in proteins that differ in function. Alternative splicing processes are prevalent in the brain, have varying incidence across brain regions, and can present sexual dimorphism. Exposure to opiates and other substances of abuse can also alter the type and incidence of the splicing process and the relative abundance of the isoforms produced. The disruption of alternative splicing patterns associated with sex differences and morphine exposure in the prefrontal cortex of a pig model was studied. The numbers of genes presenting one or more significant (FDR-adjusted p-value < 0.05) alternative splicing events were 933 and 1,368 genes when comparing females relative to males and morphine- relative to saline-treated animals, respectively. The sex-dependent opioid effect was most extreme in the contrast between morphine- versus saline-treated males with 1,934 significantly differentially spliced genes. The most frequent and significant alternative splicing type was skipped exon (∼56 % event), followed by retained intron (∼15 % events). The pathways encompassing a significant number of differentially spliced genes included axon guidance, glutamatergic synapses, circadian rhythm, and lysine degradation. Genes in these pathways included ROBO1, SEMA6C, GRIN3A, GRM2, ARNTL, CLOCK, HYKK, and DOT1L. Transcription factors ETV7 and DMAP1 presented a significant number of differentially spliced target genes. The distribution of the genes presenting differential alternative splicing in the axon guidance and circadian rhythm pathways indicates that this regulatory mechanism impacts hubs and peripheral genes. The identification of sexual dimorphism in the effect of morphine across multiple pathways confirms the necessity to explore the effects of drugs of abuse within sex. Altogether, our findings advance the understanding of the response to factors that can impact the activity of excitatory synapses by modulating transcriptional mechanisms that support the plasticity of the prefrontal cortex.

SNORA37/CMTR1/ELAVL1 feedback loop drives gastric cancer progression via facilitating CD44 alternative splicing

BACKGROUND

Emerging evidence shows that small nucleolar RNA (snoRNA), a type of highly conserved non-coding RNA, is involved in tumorigenesis and aggressiveness. However, the roles of snoRNAs in regulating alternative splicing crucial for cancer progression remain elusive.

METHODS

High-throughput RNA sequencing and comprehensive analysis were performed to identify crucial snoRNAs and downstream alternative splicing events. Biotin-labeled RNA pull-down, mass spectrometry, cross-linking RNA immunoprecipitation, and in vitro binding assays were applied to explore interaction of snoRNAs with protein partners. Alternative splicing and gene expression was observed by real-time quantitative RT-PCR and western blot assays. In vitro and in vivo studies were performed to investigate biological effects of snoRNAs and their protein partners in gastric cancer. Survival analysis was undertaken by using Kaplan-Meier method and log-rank test.

RESULTS

SNORA37 was identified as an up-regulated snoRNA essential for tumorigenesis and aggressiveness of gastric cancer. Gain- and loss-of-function studies indicated that SNORA37 promoted the growth, invasion, and metastasis of gastric cancer cells in vitro and in vivo. Mechanistically, as an ELAV like RNA binding protein 1 (ELAVL1)-generated snoRNA, SNORA37 directly bound to cap methyltransferase 1 (CMTR1) to facilitate its interaction with ELAVL1, resulting in nuclear retention and activity of ELAVL1 in regulating alternative splicing of CD44. Rescue studies revealed that SNORA37 exerted oncogenic roles in gastric cancer progression via facilitating CMTR1-ELAVL1 interaction. In clinical gastric cancer cases, high levels of SNORA37, CMTR1, ELAVL1, or CD44 were associated with shorter survival and poor outcomes of patients.

CONCLUSIONS

These results indicated that SNORA37/CMTR1/ELAVL1 feedback loop drives gastric cancer progression via facilitating CD44 alternative splicing.

scTML: a pan-cancer single-cell landscape of multiple mutation types

Investigating mutations, including single nucleotide variations (SNVs), gene fusions, alternative splicing and copy number variations (CNVs), is fundamental to cancer study. Recent computational methods and biological research have demonstrated the reliability and biological significance of detecting mutations from single-cell transcriptomic data. However, there is a lack of a single-cell-level database containing comprehensive mutation information in all types of cancer. Establishing a single-cell mutation landscape from the huge emerging single-cell transcriptomic data can provide a critical resource for elucidating the mechanisms of tumorigenesis and evolution. Here, we developed scTML (http://sctml.xglab.tech/), the first database offering a pan-cancer single-cell landscape of multiple mutation types. It includes SNVs, insertions/deletions, gene fusions, alternative splicing and CNVs, along with gene expression, cell states and other phenotype information. The data are from 74 datasets with 2 582 633 cells, including 35 full-length (Smart-seq2) transcriptomic single-cell datasets (all publicly available data with raw sequencing files), 23 datasets from 10X technology and 16 spatial transcriptomic datasets. scTML enables users to interactively explore multiple mutation landscapes across tumors or cell types, analyze single-cell-level mutation-phenotype associations and detect cell subclusters of interest. scTML is an important resource that will significantly advance deciphering intra-tumor and inter-tumor heterogeneity, and how mutations shape cell phenotypes.

Isoform-level expression of the constitutive androstane receptor (CAR or NR1I3) transcription factor better predicts the mRNA expression of the cytochrome P450s in human liver samples

Many factors cause interperson variability in the activity and expression of the cytochrome P450 (CYP) drug-metabolizing enzymes in the liver, leading to variable drug exposure and treatment outcomes. Several liver-enriched transcription factors are associated with CYP expression, with estrogen receptor α (ESR1) and constitutive androstane receptor (CAR or NR1I3) being the 2 top factors. ESR1 and NR1I3 undergo extensive alternative splicing that results in numerous splice isoforms, but how these splice isoforms associate with CYP expression is unknown. Here, we quantified 18 NR1I3 splice isoforms and the 3 most abundant ESR1 isoforms in 260 liver samples derived from African Americans (n = 125) and European Americans (n = 135). Our results showed variable splice isoform populations in the liver for both NR1I3 and ESR1. Multiple linear regression analyses revealed that compared with gene-level NR1I3, isoform-level NR1I3 expression better predicted the mRNA expression of most CYPs and 3 UDP-glucuronosyltransferases (UGTs), whereas ESR1 isoforms improved predictive models for the UGTs and CYP2D6 but not for most CYPs. Also, different NR1I3 isoforms were associated with different CYPs, and the associations varied depending on sample ancestry. Surprisingly, noncanonical NR1I3 isoforms having retained introns (introns 2 or 6) were abundantly expressed and associated with the expression of most CYPs and UGTs, whereas the reference isoform (NR1I3-205) was only associated with CYP2D6. Moreover, NR1I3 isoform diversity increased during the differentiation of induced pluripotent stem cells to hepatocytes, paralleling increasing CYP expression. These results suggest that isoform-level transcription factor expression may help to explain variation in CYP or UGT expression between individuals. SIGNIFICANCE STATEMENT: We quantified 18 NR1I3 splice isoforms and the 3 most abundant ESR1 splice isoforms in 260 liver samples derived from African American and European American donors and found variable NR1I3 and ESR1 splice isoform expression in the liver. Multiple linear regression analysis showed that, compared with gene-level expression, isoform-level expression of NR1I3 and ESR1 better predicted the mRNA expression of some cytochrome P450s and UDP-glucuronosyltransferases, highlighting the importance of isoform-level analyses to enhance our understanding of gene transcriptional regulatory networks controlling the expression of drug-metabolizing enzymes.

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

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

Effect of Selenium Nanoparticles on Alternative Splicing of Rainbow Trout Head Kidney under Heat Stress

Alternative splicing (AS) is an important post-transcriptional regulation, which can expand the functional diversity of gene products and is a mechanism for eukaryotes to cope with abiotic stress. However, there are few studies on AS events in rainbow trout under heat stress. In this study, RNA-Seq data were used to clarify the effect of selenium nanoparticles (SeNPs) on the AS events of rainbow trout head kidney under heat stress. The results showed that a total of 45,398 AS events were identified from 9804 genes, of which Skipped Exon (SE) was the most common type of AS event. Through the analysis of the differentially expressed genes (DEGs) in each group, we learned that DEGs were enriched in the spliceosome, and the relevant genes were significantly changed, which promoted the occurrence of AS. We found that lysine degradation, ubiquitin mediated proteolysis, RNA degradation, protein processing in endoplasmic reticulum processing and other pathways were significantly enriched after addition of SeNPs. In addition, some immune related signaling pathways, such as the mTOR signaling pathway, interact with each other to enhance the resistance of rainbow trout to heat stress. These results indicated that AS in head kidney of rainbow trout changed under heat stress and SeNPs played a key role in alleviating heat stress for rainbow trout.

Profiling genetically driven alternative splicing across the Indonesian archipelago

One of the regulatory mechanisms influencing the functional capacity of genes is alternative splicing (AS). Previous studies exploring the splicing landscape of human tissues have shown that AS has contributed to human biology, especially in disease progression and the immune response. Nonetheless, this phenomenon remains poorly characterized across human populations, and it is unclear how genetic and environmental variation contribute to AS. Here, we examine a set of 115 Indonesian samples from three traditional island populations spanning the genetic ancestry cline that characterizes Island Southeast Asia. We conduct a global AS analysis between islands to ascertain the degree of functionally significant AS events and their consequences. Using an event-based statistical model, we detected over 1,500 significant differential AS events across all comparisons. Additionally, we identify over 6,000 genetic variants associated with changes in splicing (splicing quantitative trait loci [sQTLs]), some of which are driven by Papuan-like genetic ancestry, and only show partial overlap with other publicly available sQTL datasets derived from other populations. Computational predictions of RNA binding activity reveal that a fraction of these sQTLs directly modulate the binding propensity of proteins involved in the splicing regulation of immune genes. Overall, these results contribute toward elucidating the role of genetic variation in shaping gene regulation in one of the most diverse regions in the world.

Mechanistic Insights into Alternative Gene Splicing in Oxidative Stress and Tissue Injury

Significance: Oxidative stress (OS) and inflammation are inducers of tissue injury. Alternative splicing (AS) is an essential regulatory step for diversifying the eukaryotic proteome. Human diseases link AS to OS; however, the underlying mechanisms must be better understood. Recent Advances: Genome‑wide profiling studies identify new differentially expressed genes induced by OS-dependent ischemia/reperfusion injury. Overexpression of RNA-binding protein RBFOX1 protects against inflammation. Hypoxia-inducible factor-1α directs polypyrimidine tract binding protein 1 to regulate mouse carcinoembryonic antigen-related cell adhesion molecule 1 (Ceacam1) AS under OS conditions. Heterogeneous nuclear ribonucleoprotein L variant 1 contains an RGG/RG motif that coordinates with transcription factors to influence human CEACAM1 AS. Hypoxia intervention involving short interfering RNAs directed to long-noncoding RNA 260 polarizes M2 macrophages toward an anti-inflammatory phenotype and alleviates OS by inhibiting IL-28RA gene AS. Critical Issues: Protective mechanisms that eliminate reactive oxygen species (ROS) are important for resolving imbalances that lead to chronic inflammation. Defects in AS can cause ROS generation, cell death regulation, and the activation of innate and adaptive immune factors. We propose that AS pathways link redox regulation to the activation or suppression of the inflammatory response during cellular stress. Future Directions: Emergent studies using molecule-mediated RNA splicing are being conducted to exploit the immunogenicity of AS protein products. Deciphering the mechanisms that connect misspliced OS and pathologies should remain a priority. Controlled release of RNA directly into cells with clinical applications is needed as the demand for innovative nucleic acid delivery systems continues to be demonstrated.

Alternative splicing of the conserved drug-resistant orthologue FpNcb2 is associated with its nuclear accumulation of products and full virulence of Fusarium pseudograminearum

BACKGROUND

Negative cofactor 2 NC2β (Ncb2 or Dr1) is the beta subunit of a conserved heterodimeric regulator of transcription negative cofactor 2 (NC2) complex that has been identified as key regulator of drug resistance in model fungi. However, its role in plant pathogens is still unclear.

RESULTS

We identified an NC2β orthologue, FpNcb2, in Fusarium pseudograminearum, which is not only a significant regulatory function in drug resistance, but also essential for growth, conidiation and penetration. Moreover, FpNcb2 undergoes alternative splicing which creates two mRNA isoforms. As a putative CCAAT binding protein, FpNcb2 concentrates in the nuclei, contributing to the expression of two spliced mRNA of FpNcb2 in hypha, conidiophores and conidia, with exception of FpNcb2ISOA in germlings. Expression of each spliced mRNA of FpNcb2 in Δfpncb2 mutant could full complement the defects on growth, conidiation and fungicides sensitivity to that of wild type. However, FpNcb2ISOA and FpNcb2ISOB have different effects on virulence. FpNcb2 acts as a regulator for the transcription of some genes encoding drug efflux and hydrolases.

CONCLUSION

Our analysis showed the existence of alternative mRNA splicing in the NC2β orthologue, which is associated with protein subcellular localization and fungal virulence. The further elucidation of the target genes of NC2β will provide insights into the potential regulation mechanisms in the antifungal resistance and pathogenesis of F. pseudograminearum. © 2024 Society of Chemical Industry.

Cortexa: a comprehensive resource for studying gene expression and alternative splicing in the murine brain

BACKGROUND

Gene expression and alternative splicing are strictly regulated processes that shape brain development and determine the cellular identity of differentiated neural cell populations. Despite the availability of multiple valuable datasets, many functional implications, especially those related to alternative splicing, remain poorly understood. Moreover, neuroscientists working primarily experimentally often lack the bioinformatics expertise required to process alternative splicing data and produce meaningful and interpretable results. Notably, re-analyzing publicly available datasets and integrating them with in-house data can provide substantial novel insights. However, such analyses necessitate developing harmonized data handling and processing pipelines which in turn require considerable computational resources and in-depth bioinformatics expertise.

RESULTS

Here, we present Cortexa-a comprehensive web portal that incorporates RNA-sequencing datasets from the mouse cerebral cortex (longitudinal or cell-specific) and the hippocampus. Cortexa facilitates understandable visualization of the expression and alternative splicing patterns of individual genes. Our platform provides SplicePCA-a tool that allows users to integrate their alternative splicing dataset and compare it to cell-specific or developmental neocortical splicing patterns. All standardized gene expression and alternative splicing datasets can be downloaded for further in-depth downstream analysis without the need for extensive preprocessing.

CONCLUSIONS

Cortexa provides a robust and readily available resource for unraveling the complexity of gene expression and alternative splicing regulatory processes in the mouse brain. The data portal is available at https://cortexa-rna.com/.

Transcriptomic Alterations in Spliceosome Components in Advanced Heart Failure: Status of Cardiac-Specific Alternative Splicing Factors

Heart failure (HF) is associated with global changes in gene expression. Alternative mRNA splicing (AS) is a key regulatory mechanism underlying these changes. However, the whole status of molecules involved in the splicing process in human HF is unknown. Therefore, we analysed the spliceosome transcriptome in cardiac tissue (n = 36) from control subjects and HF patients (with ischaemic (ICM) and dilated (DCM) cardiomyopathies) using RNA-seq. We found greater deregulation of spliceosome machinery in ICM. Specifically, we showed widespread upregulation of the E and C complex components, highlighting an increase in SNRPD2 (FC = 1.35, p < 0.05) and DHX35 (FC = 1.34, p < 0.001) mRNA levels. In contrast, we observed generalised downregulation of the A complex and cardiac-specific AS factors, such as the multifunctional protein PCBP2 (FC = -1.29, p < 0.001) and the RNA binding proteins QKI (FC = -1.35, p < 0.01). In addition, we found a relationship between SNPRD2 (an E complex component) and the left ventricular mass index in ICM patients (r = 0.779; p < 0.01). On the other hand, we observed the specific underexpression of DDX46 (FC = -1.29), RBM17 (FC = -1.33), SDE2 (FC = -1.35) and RBFOX1 (FC = -1.33), p < 0.05, in DCM patients. Therefore, these aetiology-related alterations may indicate the differential involvement of the splicing process in the development of ICM and DCM.

Large-scale evaluation of the ability of RNA-binding proteins to activate exon inclusion

RNA-binding proteins (RBPs) modulate alternative splicing outcomes to determine isoform expression and cellular survival. To identify RBPs that directly drive alternative exon inclusion, we developed tethered function luciferase-based splicing reporters that provide rapid, scalable and robust readouts of exon inclusion changes and used these to evaluate 718 human RBPs. We performed enhanced cross-linking immunoprecipitation, RNA sequencing and affinity purification-mass spectrometry to investigate a subset of candidates with no prior association with splicing. Integrative analysis of these assays indicates surprising roles for TRNAU1AP, SCAF8 and RTCA in the modulation of hundreds of endogenous splicing events. We also leveraged our tethering assays and top candidates to identify potent and compact exon inclusion activation domains for splicing modulation applications. Using these identified domains, we engineered programmable fusion proteins that outperform current artificial splicing factors at manipulating inclusion of reporter and endogenous exons. This tethering approach characterizes the ability of RBPs to induce exon inclusion and yields new molecular parts for programmable splicing control.

junctionCounts: comprehensive alternative splicing analysis and prediction of isoform-level impacts to the coding sequence

Alternative splicing (AS) is emerging as an important regulatory process for complex biological processes. Transcriptomic studies therefore commonly involve the identification and quantification of alternative processing events, but the need for predicting the functional consequences of changes to the relative inclusion of alternative events remains largely unaddressed. Many tools exist for the former task, albeit each constrained to its own event type definitions. Few tools exist for the latter task; each with significant limitations. To address these issues we developed junctionCounts, which captures both simple and complex pairwise AS events and quantifies them with straightforward exon-exon and exon-intron junction reads in RNA-seq data, performing competitively among similar tools in terms of sensitivity, false discovery rate and quantification accuracy. Its partner utility, cdsInsertion, identifies transcript coding sequence (CDS) information via in silico translation from annotated start codons, including the presence of premature termination codons. Finally, findSwitchEvents connects AS events with CDS information to predict the impact of individual events to the isoform-level CDS. We used junctionCounts to characterize splicing dynamics and NMD regulation during neuronal differentiation across four primates, demonstrating junctionCounts' capacity to robustly characterize AS in a variety of organisms and to predict its effect on mRNA isoform fate.

LncRNA H19 promoted alcohol-associated liver disease through dysregulation of alternative splicing and methionine metabolism

BACKGROUND AND AIMS

Long noncoding RNAs constitute a significant portion of the human genome. Among these, lncRNA H19, initially identified for its high expression during fetal development followed by a decline in the liver postnatally, re-emerges in various liver diseases. However, its specific role in alcohol-associated liver disease (ALD) remains unclear.

APPROACH AND RESULTS

Elevated H19 levels were detected in peripheral blood and livers of patients with alcohol-associated cirrhosis and hepatitis, as well as in livers of ethanol-fed mice. Hepatic overexpression of H19 exacerbated ethanol-induced liver steatosis and injury. Metabolomics analysis revealed decreased methionine levels in H19-overexpressed mouse livers, attributable to H19-mediated inhibition of betaine homocysteine methyltransferase (BHMT), a crucial enzyme in methionine synthesis. H19 regulated BHMT alternative splicing through polypyrimidine tract-binding protein 1 (PTBP1), resulting in a reduced Bhmt protein-coding variant. The maternally specific knockout of H19 (H19Mat+/-) or liver-specific knockout of the H19 differentially methylated domain (H19DMDHep-/-) in ethanol-fed mice upregulated BHMT expression and ameliorated hepatic steatosis. Furthermore, BHMT restoration counteracted H19-induced ethanol-mediated hepatic steatosis.

CONCLUSIONS

This study identifies a novel mechanism whereby H19, via PTBP1-mediated BHMT regulation, influences methionine metabolism in ALD. Targeting the H19-PTBP1-BHMT pathway may offer new therapeutic avenues for ALD.

Single-cell long-read targeted sequencing reveals transcriptional variation in ovarian cancer

Single-cell RNA sequencing predominantly employs short-read sequencing to characterize cell types, states and dynamics; however, it is inadequate for comprehensive characterization of RNA isoforms. Long-read sequencing technologies enable single-cell RNA isoform detection but are hampered by lower throughput and unintended sequencing of artifacts. Here we develop Single-cell Targeted Isoform Long-Read Sequencing (scTaILoR-seq), a hybridization capture method which targets over a thousand genes of interest, improving the median number of on-target transcripts per cell by 29-fold. We use scTaILoR-seq to identify and quantify RNA isoforms from ovarian cancer cell lines and primary tumors, yielding 10,796 single-cell transcriptomes. Using long-read variant calling we reveal associations of expressed single nucleotide variants (SNVs) with alternative transcript structures. Phasing of SNVs across transcripts enables the measurement of allelic imbalance within distinct cell populations. Overall, scTaILoR-seq is a long-read targeted RNA sequencing method and analytical framework for exploring transcriptional variation at single-cell resolution.

Exposure to elevated glucocorticoid during development primes altered transcriptional responses to acute stress in adulthood

Early life stress (ELS) is a major risk factor for developing psychiatric disorders, with glucocorticoids (GCs) implicated in mediating its effects in shaping adult phenotypes. In this process, exposure to high levels of developmental GC (hdGC) is thought to induce molecular changes that prime differential adult responses. However, identities of molecules targeted by hdGC exposure are not completely known. Here, we describe lifelong molecular consequences of hdGC exposure using a newly developed zebrafish double-hit stress model, which shows altered behaviors and stress hypersensitivity in adulthood. We identify a set of primed genes displaying altered expression only upon acute stress in hdGC-exposed adult fish brains. Interestingly, this gene set is enriched in risk factors for psychiatric disorders in humans. Lastly, we identify altered epigenetic regulatory elements following hdGC exposure. Thus, our study provides comprehensive datasets delineating potential molecular targets mediating the impact of hdGC exposure on adult responses.

Targeting the conserved active site of splicing machines with specific and selective small molecule modulators

The self-splicing group II introns are bacterial and organellar ancestors of the nuclear spliceosome and retro-transposable elements of pharmacological and biotechnological importance. Integrating enzymatic, crystallographic, and simulation studies, we demonstrate how these introns recognize small molecules through their conserved active site. These RNA-binding small molecules selectively inhibit the two steps of splicing by adopting distinctive poses at different stages of catalysis, and by preventing crucial active site conformational changes that are essential for splicing progression. Our data exemplify the enormous power of RNA binders to mechanistically probe vital cellular pathways. Most importantly, by proving that the evolutionarily-conserved RNA core of splicing machines can recognize small molecules specifically, our work provides a solid basis for the rational design of splicing modulators not only against bacterial and organellar introns, but also against the human spliceosome, which is a validated drug target for the treatment of congenital diseases and cancers.

The Expansion of Sirtuin Gene Family in Gilthead Sea Bream (Sparus aurata)-Phylogenetic, Syntenic, and Functional Insights across the Vertebrate/Fish Lineage

The Sirtuin (SIRT1-7) family comprises seven evolutionary-conserved enzymes that couple cellular NAD availability with health, nutrition and welfare status in vertebrates. This study re-annotated the sirt3/5 branch in the gilthead sea bream, revealing three paralogues of sirt3 (sirt3.1a/sirt3.1b/sirt3.2) and two of sirt5 (sirt5a/sirt5b) in this Perciform fish. The phylogeny and synteny analyses unveiled that the Sirt3.1/Sirt3.2 dichotomy was retained in teleosts and aquatic-living Sarcopterygian after early vertebrate 2R whole genome duplication (WGD). Additionally, only certain percomorphaceae and gilthead sea bream showed a conserved tandem-duplicated synteny block involving the mammalian-clustered sirt3.1 gene (psmd13-sirt3.1a/b-drd4-cdhr5-ctsd). Conversely, the expansion of the Sirt5 branch was shaped by the teleost-specific 3R WGD. As extensively reviewed in the literature, human-orthologues (sirt3.1/sirt5a) showed a high, conserved expression in skeletal muscle that increased as development advanced. However, recent sirt3.2 and sirt5b suffered an overall muscle transcriptional silencing across life, as well as an enhanced expression on immune-relevant tissues and gills. These findings fill gaps in the ontogeny and differentiation of Sirt genes in the environmentally adaptable gilthead sea bream, becoming a good starting point to advance towards a full understanding of its neo-functionalization. The mechanisms originating from these new paralogs also open new perspectives in the study of cellular energy sensing processes in vertebrates.

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.

On a kneading theory for gene-splicing

Two well-known facets in protein synthesis in eukaryotic cells are transcription of DNA to pre-RNA in the nucleus and the translation of messenger-RNA (mRNA) to proteins in the cytoplasm. A critical intermediate step is the removal of segments (introns) containing ∼97% of the nucleic-acid sites in pre-RNA and sequential alignment of the retained segments (exons) to form mRNA through a process referred to as splicing. Alternative forms of splicing enrich the proteome while abnormal splicing can enhance the likelihood of a cell developing cancer or other diseases. Mechanisms for splicing and origins of splicing errors are only partially deciphered. Our goal is to determine if rules on splicing can be inferred from data analytics on nucleic-acid sequences. Toward that end, we represent a nucleic-acid site as a point in a plane defined in terms of the anterior and posterior sub-sequences of the site. The "point-set" representation expands analytical approaches, including the use of statistical tools, to characterize genome sequences. It is found that point-sets for exons and introns are visually different, and that the differences can be quantified using a family of generalized moments. We design a machine-learning algorithm that can recognize individual exons or introns with 91% accuracy. Point-set distributions and generalized moments are found to differ between organisms.

[Targeting the spliceosome: A new therapeutic strategy to counteract chemotherapy resistance in lung cancer?]

Lung cancer is the first cancer-related cause of death worldwide. This is in partially due to therapeutic resistance, which occurs in around 70% of patients, especially those receiving platinum salts, the gold-standard chemotherapy. The massive deregulation of alternative transcript splicing processes observed in many cancers has led to the development of a new class of pharmacological agents aimed at inhibiting the activity of the splicing machinery (spliceosome). The molecular mechanisms by which these inhibitors act remain largely unknown, as do the benefits of using them in combination with other therapies. In this context, our work is focused on an inhibitor of the SRPK1 kinase, a major regulator of the spliceosome.

Isoform alterations in the ubiquitination machinery impacting gastrointestinal malignancies

The advancement of RNAseq and isoform-specific expression platforms has led to the understanding that isoform changes can alter molecular signaling to promote tumorigenesis. An active area in cancer research is uncovering the roles of ubiquitination on spliceosome assembly contributing to transcript diversity and expression of alternative isoforms. However, the effects of isoform changes on functionality of ubiquitination machineries (E1, E2, E3, E4, and deubiquitinating (DUB) enzymes) influencing onco- and tumor suppressor protein stabilities is currently understudied. Characterizing these changes could be instrumental in improving cancer outcomes via the identification of novel biomarkers and targetable signaling pathways. In this review, we focus on highlighting reported examples of direct, protein-coded isoform variation of ubiquitination enzymes influencing cancer development and progression in gastrointestinal (GI) malignancies. We have used a semi-automated system for identifying relevant literature and applied established systems for isoform categorization and functional classification to help structure literature findings. The results are a comprehensive snapshot of known isoform changes that are significant to GI cancers, and a framework for readers to use to address isoform variation in their own research. One of the key findings is the potential influence that isoforms of the ubiquitination machinery have on oncoprotein stability.

Therapeutic potential of alternative splicing in cardiovascular diseases

RNA splicing is an important RNA processing step required by multiexon protein-coding mRNAs and some noncoding RNAs. Precise RNA splicing is required for maintaining gene and cell function; however, mis-spliced RNA transcripts can lead to loss- or gain-of-function effects in human diseases. Mis-spliced RNAs induced by gene mutations or the dysregulation of splicing regulators may result in frameshifts, nonsense-mediated decay (NMD), or inclusion/exclusion of exons. Genetic animal models have characterised multiple splicing factors required for cardiac development or function. Moreover, sarcomeric and ion channel genes, which are closely associated with cardiovascular function and disease, are hotspots for AS. Here, we summarise splicing factors and their targets that are associated with cardiovascular diseases, introduce some therapies potentially related to pathological AS targets, and raise outstanding questions and future directions in this field.

Systematic characterization of protein structural features of alternative splicing isoforms using AlphaFold 2

Alternative splicing is an important cellular process in eukaryotes, altering pre-mRNA to yield multiple protein isoforms from a single gene. However, our understanding of the impact of alternative splicing events on protein structures is currently constrained by a lack of sufficient protein structural data. To address this limitation, we employed AlphaFold 2, a cutting-edge protein structure prediction tool, to conduct a comprehensive analysis of alternative splicing for approximately 3,000 human genes, providing valuable insights into its impact on the protein structural. Our investigation employed state of the art high-performance computing infrastructure to systematically characterize structural features in alternatively spliced regions and identified changes in protein structure following alternative splicing events. Notably, we found that alternative splicing tends to alter the structure of residues primarily located in coils and beta-sheets. Our research highlighted a significant enrichment of loops and highly exposed residues within human alternatively spliced regions. Specifically, our examination of the Septin-9 protein revealed potential associations between loops and alternative splicing, providing insights into its evolutionary role. Furthermore, our analysis uncovered two missense mutations in the Tau protein that could influence alternative splicing, potentially contributing to the pathogenesis of Alzheimer's disease. In summary, our work, through a thorough statistical analysis of extensive protein structural data, sheds new light on the intricate relationship between alternative splicing, evolution, and human disease.

Unsupervised analysis of whole transcriptome data from human pluripotent stem cells cardiac differentiation

The main objective of the present work was to highlight differences and similarities in gene expression patterns between different pluripotent stem cell cardiac differentiation protocols, using a workflow based on unsupervised machine learning algorithms to analyse the transcriptome of cells cultured as a 2D monolayer or as 3D aggregates. This unsupervised approach effectively allowed to portray the transcriptomic changes that occurred throughout the differentiation processes, with a visual representation of the entire transcriptome. The results allowed to corroborate previously reported data and also to unveil new gene expression patterns. In particular, it was possible to identify a correlation between low cardiomyocyte differentiation efficiencies and the early expression of a set of non-mesodermal genes, which can be further explored as predictive markers of differentiation efficiency. The workflow here developed can also be applied to analyse other stem cell differentiation transcriptomic datasets, envisaging future clinical implementation of cellular therapies.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Kallikrein-Related Peptidase 4 Promotes Proliferation, Migration, Invasion, and Pro-Angiogenesis of Endometrial Stromal Cells via Regulation of Brain-Derived Neurotrophic Factor Production in Endometriosis

Endometriosis is a common benign gynecologic condition. Endometriosis lesions are associated with endometrial cell proliferation, migration, invasion, and neovascularization, while the specific molecular mechanisms are still elusive. Transcriptome sequencing has been used for the identification of diagnostic markers in endometriosis. Here, transcriptome profiling revealed that kallikrein-related peptidase 4 (KLK4) expression was up-regulated in ectopic endometrium (EC) tissues of patients with endometriosis. KLK4 mediates the degradation of extracellular matrix proteins, and its proteolytic activity activates many tumorigenic and metastatic pathways via tumor invasion and migration. Nevertheless, whether KLK4 serves as an important regulatory factor in endometriosis remains unclear. This study confirmed that KLK4 was highly expressed in ectopic endometrial stromal cells (EC-ESCs). KLK4 overexpression promoted proliferation and suppressed apoptosis of EC-ESCs, induced cell migration and invasion, and enhanced angiogenesis in vivo. Mechanistically, KLK4 overexpression mediated the protein cleavage of pro-brain-derived neurotrophic factor in EC-ESCs. Finally, brain-derived neurotrophic factor was a vital downstream substrate of KLK4 maintained the proliferation, metastasis, and pro-angiogenesis abilities and inhibited apoptosis of ESCs through a rescue study. Together, these findings demonstrate the promotive role of KLK4 in endometriosis development. In addition, the study provides a new insight that KLK4 might be a potential therapeutic target and prognostic marker for patients with endometriosis.

CCAR-1 works together with the U2AF large subunit UAF-1 to regulate alternative splicing

The Cell Division Cycle and Apoptosis Regulator (CCAR) protein family members have recently emerged as regulators of alternative splicing and transcription, as well as having other key physiological functions. For example, mammalian CCAR2/DBC1 forms a complex with the zinc factor protein ZNF326 to integrate alternative splicing with RNA polymerase II transcriptional elongation in AT-rich regions of the DNA. Additionally, Caenorhabditis elegans CCAR-1, a homolog to mammalian CCAR2, facilitates the alternative splicing of the perlecan unc-52 gene. However, much about the CCAR family's role in alternative splicing is unknown. Here, we have examined the role of CCAR-1 in genome-wide alternative splicing in Caenorhabditis elegans and have identified new alternative splicing targets of CCAR-1 using RNA sequencing. Also, we found that CCAR-1 interacts with the spliceosome factors UAF-1 and UAF-2 using mass spectrometry, and that knockdown of ccar-1 affects alternative splicing patterns, motility, and proteostasis of UAF-1 mutant worms. Collectively, we demonstrate the role of CCAR-1 in regulating global alternative splicing in C. elegans and in conjunction with UAF-1.

SplicingLore: a web resource for studying the regulation of cassette exons by human splicing factors

One challenge faced by scientists from the alternative RNA splicing field is to decode the cooperative or antagonistic effects of splicing factors (SFs) to understand and eventually predict splicing outcomes on a genome-wide scale. In this manuscript, we introduce SplicingLore, an open-access database and web resource that help to fill this gap in a straightforward manner. The database contains a collection of RNA-sequencing-derived lists of alternative exons regulated by a total of 75 different SFs. All datasets were processed in a standardized manner, ensuring valid comparisons and correlation analyses. The user can easily retrieve a factor-specific set of differentially included exons from the database or provide a list of exons and search which SF(s) control(s) their inclusion. Our simple workflow is fast and easy to run, and it ensures a reliable calculation of correlation scores between the tested datasets. As a proof of concept, we predicted and experimentally validated a novel functional cooperation between the RNA helicases DDX17 and DDX5 and the heterogeneous nuclear ribonucleoprotein C (HNRNPC) protein. SplicingLore is available at https://splicinglore.ens-lyon.fr/. Database URL:  https://splicinglore.ens-lyon.fr/.

Comparison of Alternative Splicing Landscapes Revealed by Long-Read Sequencing in Hepatocyte-Derived HepG2 and Huh7 Cultured Cells and Human Liver Tissue

The long-read RNA sequencing developed by Oxford Nanopore Technologies provides a direct quantification of transcript isoforms, thereby making it possible to present alternative splicing (AS) profiles as arrays of single splice variants with different abundances. Additionally, AS profiles can be presented as arrays of genes characterized by the degree of alternative splicing (the DAS-the number of detected splice variants per gene). Here, we successfully utilized the DAS to reveal biological pathways influenced by the alterations in AS in human liver tissue and the hepatocyte-derived malignant cell lines HepG2 and Huh7, thus employing the mathematical algorithm of gene set enrichment analysis. Furthermore, analysis of the AS profiles as abundances of single splice variants by using the graded tissue specificity index τ provided the selection of the groups of genes expressing particular splice variants specifically in liver tissue, HepG2 cells, and Huh7 cells. The majority of these splice variants were translated into proteins products and appeal to be in focus regarding further insights into the mechanisms underlying cell malignization. The used metrics are intrinsically suitable for transcriptome-wide AS profiling using long-read sequencing.

Decoding the role of aberrant RNA alternative splicing in hepatocellular carcinoma: a comprehensive review

During eukaryotic gene expression, alternative splicing of messenger RNA precursors is critical in increasing protein diversity and regulatory complexity. Multiple transcript isoforms could be produced by alternative splicing from a single gene; they could eventually be translated into protein isoforms with deleted, added, or altered domains or produce transcripts containing premature termination codons that could be targeted by nonsense-mediated mRNA decay. Alternative splicing can generate proteins with similar, different, or even opposite functions. Increasingly strong evidence indicates that abnormal RNA splicing is a prevalent and crucial occurrence in cellular differentiation, tissue advancement, and the development and progression of cancer. Aberrant alternative splicing could affect cancer cell activities such as growth, apoptosis, invasiveness, drug resistance, angiogenesis, and metabolism. This systematic review provides a comprehensive overview of the impact of abnormal RNA alternative splicing on the development and progression of hepatocellular carcinoma.

Analysis of alternative splicing in chicken macrophages transfected with overexpression/knockdown of RIP2 gene

Receptor-interacting protein 2 (RIP2) plays a critical role in the transduction of many signaling pathways and is associated with many diseases. Alternative splicing (AS) is an essential and ubiquitous regulatory mechanism of gene expression that contributes to distinct transcript variants and many different kinds of proteins. In this present study, we characterized genome-wide AS events in wild-type chicken macrophages (WT) and RIP2 overexpression/knockdown chicken macrophages (oeRIP2/shRIP2) by high-throughput RNA sequencing technology. A total of 1901, 2061, and 817 differentially expressed (DE) AS genes were identified in the comparison of oeRIP2 vs. WT, oeRIP2 vs. shRIP2, and shRIP2 vs. WT, respectively. These DE AS genes participated in many important KEGG pathways, including regulation of autophagy, Wnt signaling pathway, Ubiquitin mediated proteolysis, MAPK signaling pathway, and Focal adhesion, etc. In conclusion, this research provided a broad atlas of the genome-wide scale of the AS event landscape in RIP2 overexpression/knockdown and wild-type chicken macrophages. This research also provides the theoretical basis of the gene network related to RIP2.

Specific recognition and sensitive quantification of mRNA splice variants via one-pot ligation-dependent loop-mediated isothermal amplification

Specific recognition and sensitive quantification of mRNA alternative splice variants have been a necessity for exploring the regulatory mechanism of RNA splicing and revealing the association between pre-mRNA splicing and transcriptome function, as well as disease diagnosis. However, their wide abundance range and high sequence homology pose enormous challenges for high sensitivity and selectivity quantification of splice variants. Herein, taking advantage of the excellent specificity of ligation and the powerful nucleic acid replication feature of loop-mediated isothermal amplification (LAMP), we developed a one-pot method (termed one-pot ligation-LAMP) for specific recognition and sensitive quantification of mRNA splicing variants based on two splicing junction-specific stem-loop DNA probe ligation and the subsequently initiating LAMP. The one-pot ligation-LAMP can specifically detect as low as 100 aM mRNA splice variants without any nonspecific signals and quantify them with a wide dynamics range spanning at least six orders of magnitude. We have demonstrated that the one-pot ligation-LAMP is a versatile and practical strategy for accurately quantifying different splicing variants in complex biological samples with high sensitivity all in one tube within 90 min, thereby providing an attractive tool for mRNA splice variant-related studies.

Neoantigen identification: Technological advances and challenges

Neoantigens have emerged as promising targets for cutting-edge immunotherapies, such as cancer vaccines and adoptive cell therapy. These neoantigens are unique to tumors and arise exclusively from somatic mutations or non-genomic aberrations in tumor proteins. They encompass a wide range of alterations, including genomic mutations, post-transcriptomic variants, and viral oncoproteins. With the advancements in technology, the identification of immunogenic neoantigens has seen rapid progress, raising new opportunities for enhancing their clinical significance. Prediction of neoantigens necessitates the acquisition of high-quality samples and sequencing data, followed by mutation calling. Subsequently, the pipeline involves integrating various tools that can predict the expression, processing, binding, and recognition potential of neoantigens. However, the continuous improvement of computational tools is constrained by the availability of datasets which contain validated immunogenic neoantigens. This review article aims to provide a comprehensive summary of the current knowledge as well as limitations in neoantigen prediction and validation. Additionally, it delves into the origin and biological role of neoantigens, offering a deeper understanding of their significance in the field of cancer immunotherapy. This article thus seeks to contribute to the ongoing efforts to harness neoantigens as powerful weapons in the fight against cancer.

Phosphorylation mediated regulation of RNA splicing in plants

For the past two decades, the study of alternative splicing (AS) and its involvement in plant development and stress response has grown in popularity. Only recently however, has the focus shifted to the study of how AS regulation (or lack-thereof) affects downstream mRNA and protein landscapes and how these AS regulatory events impact plant development and stress tolerance. In humans, protein phosphorylation represents one of the predominant mechanisms by which AS is regulated and thus the protein kinases governing these phosphorylation events are of interest for further study. Large-scale phosphoproteomic studies in plants have consistently found that RNA splicing-related proteins are extensively phosphorylated, however, the signaling pathways involved in AS regulation have not been resolved. In this mini-review, we summarize our current knowledge of the three major splicing-related protein kinase families in plants that are suggested to mediate AS phospho-regulation and draw comparisons to their metazoan orthologs. We also summarize and contextualize the phosphorylation events identified as occurring on splicing-related protein families to illustrate the high degree to which splicing-related proteins are modified, placing a new focus on elucidating the impacts of AS at the protein and PTM-level.

The RNA interactome in the Hallmarks of Cancer

Ribonucleic acid (RNA) molecules are indispensable for cellular homeostasis in healthy and malignant cells. However, the functions of RNA extend well beyond that of a protein-coding template. Rather, both coding and non-coding RNA molecules function through critical interactions with a plethora of cellular molecules, including other RNAs, DNA, and proteins. Deconvoluting this RNA interactome, including the interacting partners, the nature of the interaction, and dynamic changes of these interactions in malignancies has yielded fundamental advances in knowledge and are emerging as a novel therapeutic strategy in cancer. Here, we present an RNA-centric review of recent advances in the field of RNA-RNA, RNA-protein, and RNA-DNA interactomic network analysis and their impact across the Hallmarks of Cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

ROCK2-RNA interaction map reveals multiple biological mechanisms underlying tumor progression in renal cell carcinoma

Renal cell carcinoma (RCC) is the most common form of kidney cancer in adults. Despite new therapeutic modalities, the outcomes for RCC patients remain unsatisfactory. Rho-associated coiled-coil forming protein kinase 2 (ROCK2) has previously been shown to be upregulated in RCC, and its expression was negatively correlated with patient survival. However, the precise molecular function of ROCK2 has remained unclear. Herein, using RNA-seq analysis of ROCK2 knockdown and control cells, we identified 464 differentially expressed genes, and 1287 alternative splicing events in 786-O RCC cells. Furthermore, mapping of iRIP-seq reads in 786-O cells showed a biased distribution at 5' UTR, intronic and intergenic regions. By comparing ROCK2-regulated alternative splicing and iRIP-seq data, we found 292 overlapping genes that are enriched in multiple tumorigenic pathways. Taken together, our work defined a complex ROCK2-RNA interaction map on a genomic scale in a human RCC cell line, which deepens our understanding of the molecular function of ROCK2 in cancer development.

Mapping the splicing landscape of the human immune system

Most human genes code for more than one transcript. Different ratios of transcripts of the same gene can be found in different cell types or states, indicating differential use of transcription start sites or differential splicing. Such differential transcript use (DTUs) events provide an additional layer of regulation and protein diversity. With the exceptions of PTPRC and CIITA, there are very few reported cases of DTU events in the immune system. To rigorously map DTUs between different human immune cell types, we leveraged four publicly available RNA sequencing datasets. We identified 282 DTU events between five human healthy immune cell types that appear in at least two datasets. The patterns of the DTU events were mostly cell-type-specific or lineage-specific, in the context of the five cell types tested. DTUs correlated with the expression pattern of potential regulators, namely, splicing factors and transcription factors. Of the several immune related conditions studied, only sepsis affected the splicing of more than a few genes and only in innate immune cells. Taken together, we map the DTUs landscape in human peripheral blood immune cell types, and present hundreds of genes whose transcript use changes between cell types or upon activation.

Structure prediction of novel isoforms from uveal melanoma by AlphaFold

Alternative splicing is an important mechanism that enhances protein functional diversity. To date, our understanding of alternative splicing variants has been based on mRNA transcript data, but due to the difficulty in predicting protein structures, protein tertiary structures have been largely unexplored. However, with the release of AlphaFold, which predicts three-dimensional models of proteins, this challenge is rapidly being overcome. Here, we present a dataset of 315 predicted structures of abnormal isoforms in 18 uveal melanoma patients based on second- and third-generation transcriptome-sequencing data. This information comprises a high-quality set of structural data on recurrent aberrant isoforms that can be used in multiple types of studies, from those aimed at revealing potential therapeutic targets to those aimed at recognizing of cancer neoantigens at the atomic level.

PHF5A as a new OncoTarget and therapeutic prospects

PHF5A (PHD-finger domain protein 5A) is a highly conserved protein comprised of 110 amino acids that belong to PHD zinc finger proteins and is ubiquitously expressed in entire eukaryotic nuclei from yeast to man. PHF5A is an essential component of the SF3B splicing complex regulating protein-protein or protein-DNA interactions; particularly involved in pre-mRNA splicing. Besides its basic spliceosome-associated attributes encompassing the regulation of alternative splicing of specific genes, PHF5A also plays a pivotal role in cell cycle regulation and morphological development of cells along with their differentiation into particular tissues/organs, DNA damage repair, maintenance of pluripotent embryonic stem cells (CSCs) embryogenesis and regulation of chromatin-mediated transcription. Presently identification of spliceosome and non-spliceosome-associated attributes of PHF5A needs great attention based on its key involvement in the pathogenesis of cancer malignancies including the prognosis of lung adenocarcinoma, endometrial adenocarcinoma, breast, and colorectal cancer. PHF5A is an essential splicing factor or cofactor actively participating as an oncogenic protein in tumorigenesis via activation of downstream signaling pathway attributed to its regulation of dysregulated splicing or abnormal alternative splicing of targeted genes. Further, the participation of PHF5A in regulating the growth of cancer stem cells might not be ignored. The current review briefly overviews the structural and functional attributes of PHF5A along with its hitherto described role in the propagation of cancer malignancies and its future concern as a potential therapeutic target for cancer management/treatment.

TNFα Causes a Shift in Gene Expression of TNFRSF1A and TNFRSF1B Isoforms

Alternative splicing is a part of mRNA processing that expands the diversity of proteins encoded by a single gene. Studying the full range of proteins-products of translation of alternatively spliced mRNA is extremely important for understanding the interactions between receptor proteins and ligands since different receptor protein isoforms can provide variation in the activation of signaling pathways. In this study, we investigated the expression of isoforms of TNFR1 and TNFR2 receptors before and after exposure to TNFα in two cell lines that had previously demonstrated diverse effects on cell proliferation under TNFα incubation using RT-qPCR. We found that after incubation with TNFα: (1) expression of isoform 3 of the TNFRSF1A gene was increased in both cell lines; (2) the cell line with increased proliferation, K562, had decreased expression of isoforms 1 and 4 of the TNFRSF1A gene and expression of isoform 2 of TNFRSF1B gene was absent at all; (3) the cell line with decreased proliferation-MCF-7 had significantly increased expression of isoform 2 of TNFRSF1B gene. Thus, we can conclude that TNFα exposure to the K562 and MCF-7 cell lines leads to changes in the expression of TNFα receptor isoforms, which, in turn, can appear via diverse proliferative effects.

RNA Transcript Diversity in Neuromuscular Research

Three decades since the Human Genome Project began, scientists have now identified more then 25,000 protein coding genes in the human genome. The vast majority of the protein coding genes (>  90%) are multi-exonic, with the coding DNA being interrupted by intronic sequences, which are removed from the pre-mRNA transcripts before being translated into proteins, a process called splicing maturation. Variations in this process, i.e. by exon skipping, intron retention, alternative 5' splice site (5'ss), 3' splice site (3'ss), or polyadenylation usage, lead to remarkable transcriptome and proteome diversity in human tissues. Given its critical biological importance, alternative splicing is tightly regulated in a tissue- and developmental stage-specific manner. The central nervous system and skeletal muscle are amongst the tissues with the highest number of differentially expressed alternative exons, revealing a remarkable degree of transcriptome complexity. It is therefore not surprising that splicing mis-regulation is causally associated with a myriad of neuromuscular diseases, including but not limited to amyotrophic lateral sclerosis (ALS), spinal muscular atrophy (SMA), Duchenne muscular dystrophy (DMD), and myotonic dystrophy type 1 and 2 (DM1, DM2). A gene's transcript diversity has since become an integral and an important consideration for drug design, development and therapy. In this review, we will discuss transcript diversity in the context of neuromuscular diseases and current approaches to address splicing mis-regulation.

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Identification of TSGA10 and GGNBP2 splicing variants in 5' untranslated region with distinct expression profiles in brain tumor samples

Introduction

Brain tumors (BTs) are perceived as one of the most common malignancies among children. The specific regulation of each gene can play a critical role in cancer progression. The present study aimed to determine the transcripts of the TSGA10 and GGNBP2 genes, considering the alternative 5'UTR region, and investigating the expression of these different transcripts in BTs.

Material and methods

Public data on brain tumor microarray datasets in GEO were analyzed with R software to evaluate the expression levels of TSGA10 and GGNBP2 genes (the Pheatmap package in R was also used to plot DEGs in a heat map). In addition, to validate our in-silico data analysis, RT-PCR was performed to determine the splicing variants of TSGA10 and GGNBP2 genes in testis and brain tumor samples. The expression levels of splice variants of these genes were analyzed in 30 brain tumor samples and two testicular tissue samples as a positive control.

Results

In silico results show that the differential expression levels of TSGA10 and GGNBP2 were significant in the GEO datasets of BTs compared to normal samples (with adjusted p-value<0.05 and log fold change > 1). This study's experimental results showed that the TSGA10 gene produces four different transcripts with two distinct promoter regions and splicing exon 4. The relative mRNA expression of transcripts without exon 4 was higher than transcripts with exon 4 in BT samples (p-value<001). In GGNBP2, exon 2 in the 5'UTR region and exon 6 in the coding sequence were spliced. The expression analysis results showed that the relative mRNA expression of transcript variants without exon 2 was higher than other transcript variants with exon 2 in BT samples (p-value<001).

Conclusion

The decreased expression levels of transcripts with longer 5'UTR in BT samples than in testicular or low-grade brain tumor samples may decrease their translation efficiency. Therefore, decreased amounts of TSGA10 and GGNBP2 as potential tumor suppressor proteins, especially in high-grade brain tumors, may cause cancer development by angiogenesis and metastasis.

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Identification of Smoking-Associated Transcriptome Aberration in Blood with Machine Learning Methods

Long-term cigarette smoking causes various human diseases, including respiratory disease, cancer, and gastrointestinal (GI) disorders. Alterations in gene expression and variable splicing processes induced by smoking are associated with the development of diseases. This study applied advanced machine learning methods to identify the isoforms with important roles in distinguishing smokers from former smokers based on the expression profile of isoforms from current and former smokers collected in one previous study. These isoforms were deemed as features, which were first analyzed by the Boruta to select features highly correlated with the target variables. Then, the selected features were evaluated by four feature ranking algorithms, resulting in four feature lists. The incremental feature selection method was applied to each list for obtaining the optimal feature subsets and building high-performance classification models. Furthermore, a series of classification rules were accessed by decision tree with the highest performance. Eventually, the rationality of the mined isoforms (features) and classification rules was verified by reviewing previous research. Features such as isoforms ENST00000464835 (expressed by LRRN3), ENST00000622663 (expressed by SASH1), and ENST00000284311 (expressed by GPR15), and pathways (cytotoxicity mediated by natural killer cell and cytokine-cytokine receptor interaction) revealed by the enrichment analysis, were highly relevant to smoking response, suggesting the robustness of our analysis pipeline.

Sources of Cancer Neoantigens beyond Single-Nucleotide Variants

The success of checkpoint blockade therapy against cancer has unequivocally shown that cancer cells can be effectively recognized by the immune system and eliminated. However, the identity of the cancer antigens that elicit protective immunity remains to be fully explored. Over the last decade, most of the focus has been on somatic mutations derived from non-synonymous single-nucleotide variants (SNVs) and small insertion/deletion mutations (indels) that accumulate during cancer progression. Mutated peptides can be presented on MHC molecules and give rise to novel antigens or neoantigens, which have been shown to induce potent anti-tumor immune responses. A limitation with SNV-neoantigens is that they are patient-specific and their accurate prediction is critical for the development of effective immunotherapies. In addition, cancer types with low mutation burden may not display sufficient high-quality [SNV/small indels] neoantigens to alone stimulate effective T cell responses. Accumulating evidence suggests the existence of alternative sources of cancer neoantigens, such as gene fusions, alternative splicing variants, post-translational modifications, and transposable elements, which may be attractive novel targets for immunotherapy. In this review, we describe the recent technological advances in the identification of these novel sources of neoantigens, the experimental evidence for their presentation on MHC molecules and their immunogenicity, as well as the current clinical development stage of immunotherapy targeting these neoantigens.