Functional roles of alternative splicing factors in human disease

Nematode histone H2A variant evolution reveals diverse histories of retention and loss and evidence for conserved core-like variant histone genes

Histone variants are paralogs that replace canonical histones in nucleosomes, often imparting novel functions. However, how histone variants arise and evolve is poorly understood. Reconstruction of histone protein evolution is challenging due to large differences in evolutionary rates across gene lineages and sites. Here we used intron position data from 108 nematode genomes in combination with amino acid sequence data to find disparate evolutionary histories of the three H2A variants found in Caenorhabditis elegans: the ancient H2A.ZHTZ-1, the sperm-specific HTAS-1, and HIS-35, which differs from the canonical S-phase H2A by a single glycine-to-alanine C-terminal change. Although the H2A.ZHTZ-1 protein sequence is highly conserved, its gene exhibits recurrent intron gain and loss. This pattern suggests that specific intron sequences or positions may not be important to H2A.Z functionality. For HTAS-1 and HIS-35, we find variant-specific intron positions that are conserved across species. Patterns of intron position conservation indicate that the sperm-specific variant HTAS-1 arose more recently in the ancestor of a subset of Caenorhabditis species, while HIS-35 arose in the ancestor of Caenorhabditis and its sister group, including the genus Diploscapter. HIS-35 exhibits gene retention in some descendent lineages but gene loss in others, suggesting that histone variant use or functionality can be highly flexible. Surprisingly, we find the single amino acid differentiating HIS-35 from core H2A is ancestral and common across canonical Caenorhabditis H2A sequences. Thus, we speculate that the role of HIS-35 lies not in encoding a functionally distinct protein, but instead in enabling H2A expression across the cell cycle or in distinct tissues. This work illustrates how genes encoding such partially-redundant functions may be advantageous yet relatively replaceable over evolutionary timescales, consistent with the patchwork pattern of retention and loss of both genes. Our study shows the utility of intron positions for reconstructing evolutionary histories of gene families, particularly those undergoing idiosyncratic sequence evolution.

Editorial: The RNA revolution and cancer

RNA biology has revolutionized cancer understanding and treatment, especially in endocrine-related malignancies. This editorial highlights RNA's crucial role in cancer progression, emphasizing its influence on tumor heterogeneity and behavior. Processes like alternative splicing and noncoding RNA regulation shape cancer biology, with microRNAs, long noncoding RNAs, and circular RNAs orchestrating gene expression dynamics. Aberrant RNA signatures hold promise as diagnostic and prognostic biomarkers in endocrine-related cancers. Recent findings, such as aberrant PI3Kδ splice isoforms and epithelial-mesenchymal transition-related lncRNA signatures, unveil potential therapeutic targets for personalized treatments. Insights into m6A-associated lncRNA prognostic models and the function of lncRNA LINC00659 in gastric cancer represents ongoing research in this field. As understanding of RNA's role in cancer expands, personalized therapies offer transformative potential in managing endocrine-related malignancies. This signifies a significant stride towards precision oncology, fostering innovation for more effective cancer care.

Functional analysis of the CTNS gene exonic variants predicted to affect splicing

Cystinosis is a severe, monogenic systemic disease caused by variants in CTNS gene. Currently, there is growing evidence that exonic variants in many diseases can affect pre-mRNA splicing. The impact of CTNS gene exonic variants on splicing regulation may be underestimated due to the lack of routine studies at the RNA level. Here, we analyzed 59 exonic variants in the CTNS gene using bioinformatics tools and identified candidate variants that may induce splicing alterations by minigene assays. We identified six exonic variants that induce splicing alterations by disrupting the ratio of exonic splicing enhancers/exonic splicing silencers (ESEs/ESSs) or by interfering with the recognition of classical splice sites, or both. Our results help in the correct molecular characterization of variants in cystinosis and inform emerging therapies. Furthermore, our work suggests that the combination of in silico and in vitro assays facilitates to assess the effects of DNA variants driving rare genetic diseases on splicing regulation and will enhance the clinical utility of variant functional annotation.

2'-19F labelling of ribose in RNAs: a tool to analyse RNA/protein interactions by NMR in physiological conditions

Protein-RNA interactions are central to numerous cellular processes. In this work, we present an easy and straightforward NMR-based approach to determine the RNA binding site of RNA binding proteins and to evaluate the binding of pairs of proteins to a single-stranded RNA (ssRNA) under physiological conditions, in this case in nuclear extracts. By incorporation of a 19F atom on the ribose of different nucleotides along the ssRNA sequence, we show that, upon addition of an RNA binding protein, the intensity of the 19F NMR signal changes when the 19F atom is located near the protein binding site. Furthermore, we show that the addition of pairs of proteins to a ssRNA containing two 19F atoms at two different locations informs on their concurrent binding or competition. We demonstrate that such studies can be done in a nuclear extract that mimics the physiological environment in which these protein-ssRNA interactions occur. Finally, we demonstrate that a trifluoromethoxy group (-OCF3) incorporated in the 2'ribose position of ssRNA sequences increases the sensitivity of the NMR signal, leading to decreased measurement times, and reduces the issue of RNA degradation in cellular extracts.

Exploring Multi-Tissue Alternative Splicing and Skeletal Muscle Metabolism Regulation in Obese- and Lean-Type Pigs

Alternative splicing (AS) is a crucial mechanism in post-transcriptional regulation, contributing significantly to the diversity of the transcriptome and proteome. In this study, we performed a comprehensive AS profile in nine tissues obtained from Duroc (lean-type) and Luchuan (obese-type) pigs. Notably, 94,990 AS events from 14,393 genes were identified. Among these AS events, it was observed that 80% belonged to the skipped exon (SE) type. Functional enrichment analysis showed that genes with more than ten AS events were closely associated with tissue-specific functions. Additionally, the analysis of overlap between differentially alternative splicing genes (DSGs) and differentially expressed genes (DEGs) revealed the highest number of overlapped genes in the heart and skeletal muscle. The novelty of our study is that it identified and validated three genes (PYGM, MAPK11 and CAMK2B) in the glucagon signaling pathway, and their alternative splicing differences were highly significant across two pig breeds. In conclusion, our study offers novel insights into the molecular regulation of diverse tissue physiologies and the phenotypic differences between obese- and lean-type pigs, which are helpful for pig breeding.

Genomic insights and advanced machine learning: characterizing autism spectrum disorder biomarkers and genetic interactions

Autism Spectrum Disorder (ASD) is a complex neurodevelopmental condition characterized by altered brain connectivity and function. In this study, we employed advanced bioinformatics and explainable AI to analyze gene expression associated with ASD, using data from five GEO datasets. Among 351 neurotypical controls and 358 individuals with autism, we identified 3,339 Differentially Expressed Genes (DEGs) with an adjusted p-value (≤ 0.05). A subsequent meta-analysis pinpointed 342 DEGs (adjusted p-value ≤ 0.001), including 19 upregulated and 10 down-regulated genes across all datasets. Shared genes, pathogenic single nucleotide polymorphisms (SNPs), chromosomal positions, and their impact on biological pathways were examined. We identified potential biomarkers (HOXB3, NR2F2, MAPK8IP3, PIGT, SEMA4D, and SSH1) through text mining, meriting further investigation. Additionally, ‎we shed light on the roles of RPS4Y1 and KDM5D genes in neurogenesis and neurodevelopment. Our analysis detected 1,286 SNPs linked to ASD-related conditions, of which 14 high-risk SNPs were located on chromosomes 10 and X. We highlighted potential missense SNPs associated with FGFR inhibitors, suggesting that it may serve as a promising biomarker for responsiveness to targeted therapies. Our explainable AI model identified the MID2 gene as a potential ASD biomarker. This research unveils vital genes and potential biomarkers, providing a foundation for novel gene discovery in complex diseases.

Improving the RNA velocity approach with single-cell RNA lifecycle (nascent, mature and degrading RNAs) sequencing technologies

We presented an experimental method called FLOUR-seq, which combines BD Rhapsody and nanopore sequencing to detect the RNA lifecycle (including nascent, mature, and degrading RNAs) in cells. Additionally, we updated our HIT-scISOseq V2 to discover a more accurate RNA lifecycle using 10x Chromium and Pacbio sequencing. Most importantly, to explore how single-cell full-length RNA sequencing technologies could help improve the RNA velocity approach, we introduced a new algorithm called 'Region Velocity' to more accurately configure cellular RNA velocity. We applied this algorithm to study spermiogenesis and compared the performance of FLOUR-seq with Pacbio-based HIT-scISOseq V2. Our findings demonstrated that 'Region Velocity' is more suitable for analyzing single-cell full-length RNA data than traditional RNA velocity approaches. These novel methods could be useful for researchers looking to discover full-length RNAs in single cells and comprehensively monitor RNA lifecycle in cells.

Fluid protein fold space and its implications

Fold-switching proteins, which remodel their secondary and tertiary structures in response to cellular stimuli, suggest a new view of protein fold space. For decades, experimental evidence has indicated that protein fold space is discrete: dissimilar folds are encoded by dissimilar amino acid sequences. Challenging this assumption, fold-switching proteins interconnect discrete groups of dissimilar protein folds, making protein fold space fluid. Three recent observations support the concept of fluid fold space: (1) some amino acid sequences interconvert between folds with distinct secondary structures, (2) some naturally occurring sequences have switched folds by stepwise mutation, and (3) fold switching is evolutionarily selected and likely confers advantage. These observations indicate that minor amino acid sequence modifications can transform protein structure and function. Consequently, proteomic structural and functional diversity may be expanded by alternative splicing, small nucleotide polymorphisms, post-translational modifications, and modified translation rates.

To Explore the Binding Affinity of Human γ-Secretase Activating Protein (GSAP) Isoform 4 with APP-C99 Peptides

γ-Secretase activating protein (GSAP) is known to play an important role in the β-amyloid pathway. It acts as a modulator and accentuates the truncation of the amyloid precursor protein C-99 fragment through the γ-secretase complex. GSAP has four isoforms, out of which canonical isoform 1, a 16 kDa C-terminal portion, has been extensively studied, whereas the function of other three isoforms remains unknown. Here, we explore the GSAP isoform 4 (GSAP_I4) expression and purification from inclusion bodies followed by the refolding of the protein. The secondary structure of GSAP_I4 is predicted using circular dichroism. The protein is further characterized by western blotting and mass spectroscopy analysis. Additionally, biochemical assays and in silico molecular docking and molecular simulation are performed to investigate the binding of GSAP_I4 and APP-C99 peptide fragments. The results reflect that although GSAP_I1 and GSAP_I4 share high sequence similarity, the isoform 4 does not show any affinity toward APP-C99 peptide fragments. This hints toward the fact that GSAP_I4 might have a different role in the living system that is yet unexplored.

Roles of minor spliceosome in intron recognition and the convergence with the better understood major spliceosome

Catalyzed by spliceosomes in the nucleus, RNA splicing removes intronic sequences from precursor RNAs in eukaryotes to generate mature RNA, which also significantly increases proteome complexity and fine-tunes gene expression. Most metazoans have two coexisting spliceosomes; the major spliceosome, which removes >99.5% of introns, and the minor spliceosome, which removes far fewer introns (only 770 at present have been predicted in the human genome). Both spliceosomes are large and dynamic machineries, each consisting of five small nuclear RNAs (snRNAs) and more than 100 proteins. However, the dynamic assembly, catalysis, and protein composition of the minor spliceosome are still poorly understood. With different splicing signals, minor introns are rare and usually distributed alone and flanked by major introns in genes, raising questions of how they are recognized by the minor spliceosome and how their processing deals with the splicing of neighboring major introns. Due to large numbers of introns and close similarities between the two machinery, cooperative, and competitive recognition by the two spliceosomes has been investigated. Functionally, many minor-intron-containing genes are evolutionarily conserved and essential. Mutations in the minor spliceosome exhibit a variety of developmental defects in plants and animals and are linked to numerous human diseases. Here, we review recent progress in the understanding of minor splicing, compare currently known components of the two spliceosomes, survey minor introns in a wide range of organisms, discuss cooperation and competition of the two spliceosomes in splicing of minor-intron-containing genes, and contributions of minor splicing mutations in development and diseases. This article is categorized under: RNA Processing > Processing of Small RNAs RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

The comprehensive and systematic identification of BLCA-specific SF-regulated, survival-related AS events

Bladder urothelial carcinoma (BLCA) is a complex disease with high morbidity and mortality. Changes in alternative splicing (AS) and splicing factor (SF) can affect gene expression, thus playing an essential role in tumorigenesis. This study downloaded 412 patients' clinical information and 433 samples of transcriptome profiling data from TCGA. And we collected 48 AS signatures from SpliceSeq. LASSO and Cox analyses were used for identifying survival-related AS events in BLCA. Finally, 1,645 OS-related AS events in 1,129 genes were validated by Kaplan-Meier (KM) survival analysis, ROC analysis, risk curve analysis, and independent prognostic analysis. Finally, our survey provides an AS-SF regulation network consisting of five SFs and 46 AS events. In the end, we profiled genes that AS occurred in pan-cancer and five SFs' expression in tumor and normal samples in BLCA. We selected CLIP-seq data for validation the interaction regulated by RBP. Our study paves the way for potential therapeutic targets of BLCA.

BCL2-associated athanogene 6 exon24 contributes to testosterone synthesis and male fertility in mammals

Objectives

BCL2‐associated athanogene 6 (BAG6) plays critical roles in spermatogenesis by maintaining testicular cell survival. Our previous data showed porcine BAG6 exon24‐skipped transcript is highly expressed in immature testes compared with mature testes. The objective of this study is to reveal the functional significance of BAG6 exon24 in mammalian spermatogenesis.

Materials and Methods

CRISPR/Cas9 system was used to generate Bag6 exon24 knockout mice. Testes and cauda epididymal sperm were collected from mice. TMT proteomics analysis was used to discover the protein differences induced by Bag6 exon24 deletion. Testosterone enanthate was injected into mice to generate a high‐testosterone mice model. H&E staining, qRT‐PCR, western blotting, vector/siRNA transfection, immunofluorescence, immunoprecipitation, transmission electron microscopy, TUNEL and ELISA were performed to investigate the phenotypes and molecular basis.

Results

Bag6 exon24 knockout mice show sub‐fertility along with partially impaired blood‐testis barrier, increased apoptotic testicular cell rate and abnormal sperm morphology. Endoplasmic reticulum stress occurs in Bag6 exon24‐deficient testes and sterol regulatory element‐binding transcription factor 2 is activated; as a result, cytochrome P450 family 51 subfamily A member 1 expression is up‐regulated, which causes a high serum testosterone level. Additionally, serine/arginine‐rich splicing factor 1 down‐regulates BAG6 exon24‐skipped transcripts in porcine Sertoli cells by binding to 35–51 nt on BAG6 exon24 via its N‐terminal RNA‐recognition domain.

Conclusions

Our findings reveal the critical roles of BAG6 exon24 in testosterone biosynthesis and male fertility, which provides new insights into the regulation of spermatogenesis and pathogenesis of subfertility in mammals.

Hnrnpul1 controls transcription, splicing, and modulates skeletal and limb development in vivo

Mutations in RNA-binding proteins can lead to pleiotropic phenotypes including craniofacial, skeletal, limb, and neurological symptoms. Heterogeneous nuclear ribonucleoproteins (hnRNPs) are involved in nucleic acid binding, transcription, and splicing through direct binding to DNA and RNA, or through interaction with other proteins in the spliceosome. We show a developmental role for Hnrnpul1 in zebrafish, resulting in reduced body and fin growth and missing bones. Defects in craniofacial tendon growth and adult-onset caudal scoliosis are also seen. We demonstrate a role for Hnrnpul1 in alternative splicing and transcriptional regulation using RNA-sequencing, particularly of genes involved in translation, ubiquitination, and DNA damage. Given its cross-species conservation and role in splicing, it would not be surprising if it had a role in human development. Whole-exome sequencing detected a homozygous frameshift variant in HNRNPUL1 in 2 siblings with congenital limb malformations, which is a candidate gene for their limb malformations. Zebrafish Hnrnpul1 mutants suggest an important developmental role of hnRNPUL1 and provide motivation for exploring the potential conservation of ancient regulatory circuits involving hnRNPUL1 in human development.

RNA-Scoop: interactive visualization of transcripts in single-cell transcriptomes

Recent advances in single-cell RNA sequencing technologies have made detection of transcripts in single cells possible. The level of resolution provided by these technologies can be used to study changes in transcript usage across cell populations and help investigate new biology. Here, we introduce RNA-Scoop, an interactive cell cluster and transcriptome visualization tool to analyze transcript usage across cell categories and clusters. The tool allows users to examine differential transcript expression across clusters and investigate how usage of specific transcript expression mechanisms varies across cell groups.

Identification of tumour immune microenvironment-related alternative splicing events for the prognostication of pancreatic adenocarcinoma

PURPOSE

Pancreatic adenocarcinoma (PAAD) is characterized by low antitumour immune cell infiltration in an immunosuppressive microenvironment. This study aimed to systematically explore the impact on prognostic alternative splicing events (ASs) of tumour immune microenvironment (TIME) in PAAD.

METHODS

The ESTIMATE algorithm was implemented to compute the stromal/immune-related scores of each PAAD patient, followed by Kaplan-Meier (KM) survival analysis of patients with different scores grouped by X-tile software. TIME-related differentially expressed ASs (DEASs) were determined and evaluated through functional annotation analysis. In addition, Cox analyses were implemented to construct a TIME-related signature and an AS clinical nomogram. Moreover, comprehensive analyses, including gene set enrichment analysis (GSEA), immune infiltration, immune checkpoint gene expression, and tumour mutation were performed between the two risk groups to understand the potential mechanisms. Finally, Cytoscape was implemented to illuminate the AS-splicing factor (SF) regulatory network.

RESULTS

A total of 437 TIME-related DEASs significantly related to PAAD tumorigenesis and the formation of the TIME were identified. Additionally, a robust TIME-related prognostic signature based on seven DEASs was generated, and an AS clinical nomogram combining the signature and four clinical predictors also exhibited prominent discrimination by ROC (0.762 ~ 0.804) and calibration curves. More importantly, the fractions of CD8 T cells, regulatory T cells and activated memory CD4 T cells were lower, and the expression of four immune checkpoints-PD-L1, CD47, CD276, and PVR-was obviously higher in high-risk patients. Finally, functional analysis and tumour mutations revealed that aberrant immune signatures and activated carcinogenic pathways in high-risk patients may be the cause of the poor prognosis.

CONCLUSION

We extracted a list of DEASs associated with the TIME through the ESTIMATE algorithm and constructed a prognostic signature on the basis of seven DEASs to predict the prognosis of PAAD patients, which may guide advanced decision-making for personalized precision intervention.

Analysis of Genome-Wide Alternative Splicing Profiling and Development of Potential Drugs in Lung Adenocarcinoma

Alternative splicing (AS) is significantly related to tumor development as well as a patient’s clinical characteristics. This study was designed to systematically analyze the survival-associated AS signatures in Lung adenocarcinoma (LUAD). Among 30,735 AS events in 9,635 genes, we found that there were 1,429 AS in 1,125 genes which were conspicuously related to the overall survival of LUAD patients. Then, according to the seven types of AS events, we established AS signatures and constructed a new combined prognostic model. The Kaplan-Meier curve results showed that seven types of AS signatures and the combined prognostic model could divide patients into distinct prognoses. The ROC curve shows that all eight AS signatures had powerful predictive properties with different AUCs ranging from 0.708 to 0.849. Additionally, the elevated risk scores were positively related to higher TNM stage and metastasis. Interestingly, AS events and splicing factors (SFs) network shed light on a meaningful connection between prognostic AS genes and corresponding SFs. Moreover, we found that the combined prognostic model signature has a higher predictive ability than the mRNA signature. Furthermore, tumors at high risk might evade immune recognition by decreasing the expression of antigen presentation genes. Finally, we predicted the three most significant small molecule drugs to inhibit LUAD. Among them, NVP-AUY922 had the lowest IC50 value and might become a potential drug to prolong a patient’s survival. In conclusion, our study established a potential prognostic signature for LUAD patients, revealed a splicing network between AS and SFs and possible immune escape mechanism, and provided several small-molecule drugs to inhibit tumorigenesis.

Comprehensive Combined Proteomics and Genomics Analysis Identifies Prognostic Related Transcription Factors in Breast Cancer and Explores the Role of DMAP1 in Breast Cancer

Transcription factors (TFs) are important for regulating gene transcription and are the hallmark of many cancers. The identification of breast cancer TFs will help in developing new diagnostic and individualized cancer treatment tools. In this study, we used quantitative proteomic analyses of nuclear proteins and massive transcriptome data to identify enriched potential TFs and explore the possible role of the transcription factor DMAP1 in breast cancer. We identified 13 prognostic-related TFs and constructed their regulated genes, alternative splicing (AS) events, and splicing factor (SF) regulation networks. DMAP1 was reported less in breast cancer. The expression of DMAP1 decreased in breast cancer tumors compared with normal tissues. The poor prognosis of patients with low DMAP1 expression may relate to the activated PI3K/Akt signaling pathway, as well as other cancer-relevant pathways. This may be due to the low methylation and high expression of these pathway genes and the fact that such patients show more sensitivity to some PI3K/Akt signaling pathway inhibitors. The high expression of DMAP1 was correlated with low immune cell infiltration, and the response to immune checkpoint inhibitor treatment in patients with high DMAP1 expression was low. Our study identifies some transcription factors that are significant for breast cancer progression, which can be used as potential personalized prognostic markers in the future.

A Candidate RNAi Screen Reveals Diverse RNA-Binding Protein Phenotypes in Drosophila Flight Muscle

The proper regulation of RNA processing is critical for muscle development and the fine-tuning of contractile ability among muscle fiber-types. RNA binding proteins (RBPs) regulate the diverse steps in RNA processing, including alternative splicing, which generates fiber-type specific isoforms of structural proteins that confer contractile sarcomeres with distinct biomechanical properties. Alternative splicing is disrupted in muscle diseases such as myotonic dystrophy and dilated cardiomyopathy and is altered after intense exercise as well as with aging. It is therefore important to understand splicing and RBP function, but currently, only a small fraction of the hundreds of annotated RBPs expressed in muscle have been characterized. Here, we demonstrate the utility of Drosophila as a genetic model system to investigate basic developmental mechanisms of RBP function in myogenesis. We find that RBPs exhibit dynamic temporal and fiber-type specific expression patterns in mRNA-Seq data and display muscle-specific phenotypes. We performed knockdown with 105 RNAi hairpins targeting 35 RBPs and report associated lethality, flight, myofiber and sarcomere defects, including flight muscle phenotypes for Doa, Rm62, mub, mbl, sbr, and clu. Knockdown phenotypes of spliceosome components, as highlighted by phenotypes for A-complex components SF1 and Hrb87F (hnRNPA1), revealed level- and temporal-dependent myofibril defects. We further show that splicing mediated by SF1 and Hrb87F is necessary for Z-disc stability and proper myofibril development, and strong knockdown of either gene results in impaired localization of kettin to the Z-disc. Our results expand the number of RBPs with a described phenotype in muscle and underscore the diversity in myofibril and transcriptomic phenotypes associated with splicing defects. Drosophila is thus a powerful model to gain disease-relevant insight into cellular and molecular phenotypes observed when expression levels of splicing factors, spliceosome components and splicing dynamics are altered.

Identification of prognostic alternative splicing events in sarcoma

Sarcoma is a rare malignancy with unfavorable prognoses. Accumulating evidence indicates that aberrant alternative splicing (AS) events are generally involved in cancer pathogenesis. The aim of this study was to identify the prognostic value of AS-related survival genes as potential biomarkers, and highlight the functional roles of AS events in sarcoma. RNA-sequencing and AS-event datasets were downloaded from The Cancer Genome Atlas (TCGA) sarcoma cohort and TCGA SpliceSeq, respectively. Survival-related AS events were further assessed using a univariate analysis. A multivariate Cox regression analysis was also performed to establish a survival-gene signature to predict patient survival, and the area-under-the-curve method was used to evaluate prognostic reliability. KOBAS 3.0 and Cytoscape were used to functionally annotate AS-related genes and to assess their network interactions. We detected 9674 AS events in 40,184 genes from 236 sarcoma samples, and the 15 most significant genes were then used to construct a survival regression model. We further validated the involvement of ten potential survival-related genes (TUBB3, TRIM69, ZNFX1, VAV1, KCNN2, VGLL3, AK7, ARMC4, LRRC1, and CRIP1) in the occurrence and development of sarcoma. Multivariate survival model analyses were also performed, and validated that a model using these ten genes provided good classifications for predicting patient outcomes. The present study has increased our understanding of AS events in sarcoma, and the gene-based model using AS-related events may serve as a potential predictor to determine the survival of sarcoma patients.

Epithelial Splicing Regulatory Protein 1 Is Overexpressed in Breast Cancer and Predicts Poor Prognosis for Breast Cancer Patients

Background

Epithelial splicing regulatory proteins (ESRPs), including ESRP1 and ESRP2, are important proteins for alternative splicing of mRNAs and are reported to promote or inhibit the progression of some tumors. However, the effects of ESRPs in breast cancer are still unknown.

Material/Methods

In this study, we detected the transcriptional level and alterations of ESRP1 in patients with breast cancer based on the Oncomine, Gene Expression Profiling Interactive Analysis, Gene Expression-Based Outcome for Breast Cancer Online, and cBioPortal databases. Using immunohistochemistry and quantitative polymerase chain reaction, the expression pattern of ESRP1 in breast cancer was analyzed. Analysis of the clinicopathological characteristics and function of ESRP1 in breast cancer were actualized through the University of Alabama Cancer database and Database for Annotation, Visualization and Integrated Discovery. Using the Kaplan-Meier plotter, the prognostic values of ESRP1 in patients with breast cancer were analyzed. The Encyclopedia of RNA Interactomes database was used to predict miRNAs that regulated ESRP1.

Results

We found that ESRP1 was significantly overexpressed in patients with breast cancer, compared with patients without breast cancer, and had statistically significant clinicopathological characteristics. Kaplan-Meier plotter analysis indicated that the elevated expression of ESRP1 was associated with poor prognosis in patients with breast cancer. Furthermore, hsa-miR-181c-5p was identified to be potentially involved in the regulation of ESRP1.

Conclusions

These results suggest that ESRP1 is a valuable target for the precise treatment of breast cancer and a potential biomarker for the prognosis of patients with breast cancer.

The Functional Impact of Alternative Splicing and Single Nucleotide Polymorphisms in Rheumatoid Arthritis

Advances in genomics and proteomics aid the identification of genes associated with various diseases. Genome-Wide Association Studies (GWAS) have identified multiple loci as risk alleles for susceptibility to Rheumatoid Arthritis (RA). A bisection of RA risk can be attributed to genetic factors. Over 100 associated genetic loci that encompass immune regulatory factors have been found to be linked with RA. Aberrant Single Nucleotide Polymorphisms (SNPs) and alternative splicing mechanisms in such loci induce RA. These aberrations are viewed as potential therapeutic targets due to their association with a multitude of diseases. This review presents a few imperious genes whose alterations can cause severe bone deformities culminating in RA.

Two RECK Splice Variants (Long and Short) Are Differentially Expressed in Patients with Stable and Unstable Coronary Artery Disease: A Pilot Study

Primary prevention is crucial for coronary heart disease (CAD) and the identification of new reliable biomarkers might help risk stratification or predict adverse coronary events. Alternative splicing (AS) is a less investigated genetic factors implicated in CAD etiology. We performed an RNA-seq study on PBMCs from CAD patients and control subjects (CTR) and observed 113 differentially regulated AS events (24 up and 89 downregulated) in 86 genes. The RECK (Reversion-inducing-cysteine-rich protein with Kazal motifs) gene was further analyzed in a larger case study (24 CTR subjects, 72 CAD and 32 AMI patients) for its Splicing-Index FC (FC = −2.64; p = 0.0217), the AS event involving an exon (exon 18), and its role in vascular inflammation and remodeling. We observed a significant downregulation of Long RECK splice variant (containing exon 18) in PBMCs of AMI compared to CTR subjects (FC = −3.3; p < 0.005). Interestingly, the Short RECK splice variant (lacking exon 18) was under-expressed in AMI compared to both CTR (FC = −4.5; p < 0.0001) and CAD patients (FC = −4.2; p < 0.0001). A ROC curve, constructed combining Long and Short RECK expression data, shows an AUC = 0.81 (p < 0.001) to distinguish AMI from stable CAD patients. A significant negative correlation between Long RECK and triglycerides in CTR group and a positive correlation in the AMI group was found. The combined evaluation of Long and Short RECK expression levels is a potential genomic biomarker for the discrimination of AMI from CAD patients. Our results underline the relevance of deeper studies on the expression of these two splice variants to elucidate their functional role in CAD development and progression.

Correlations Between the Characteristics of Alternative Splicing Events, Prognosis, and the Immune Microenvironment in Breast Cancer

Objective

Alternative splicing (AS) is the mechanism by which a few genes encode numerous proteins, and it redefines the concept of gene expression regulation. Recent studies showed that dysregulation of AS was an important cause of tumorigenesis and microenvironment formation. Therefore, we performed a systematic analysis to examine the role of AS in breast cancer (Breast Cancer, BrCa) progression.

Methods

The present study included 993 BrCa patients from The Cancer Genome Atlas (TCGA) database in the genome-wide analysis of AS events. We used differential and prognostic analyses and found differentially expressed alternative splicing (DEAS) events and independent prognostic factors related to patients' overall survival (OS) and disease-free survival (DFS). We divided the patients into two groups based on these AS events and analyzed their clinical features, molecular subtyping and immune characteristics. We also constructed a splicing factor (SF) regulation network for key AS events and verified the existence of AS events in tissue samples using real-time quantitative PCR.

Results

A total of 678 AS events were identified as differentially expressed, of which 13 and 10 AS events were independent prognostic factors of patients' OS and DFS, respectively. Unsupervised clustering analysis based on these prognostic factors indicated that the Cluster 1 group had a better prognosis and more immune cell infiltration. SFs were significantly related to the expression of AS events, and AA-RPS21 was significantly upregulated in tumors.

Conclusion

Alternative splicing expands the mechanism of breast cancer progression from a new perspective. Notably, alternative splicing may affect the patient's prognosis by affecting the infiltration of immune cells. Our research provides important guidance for subsequent studies of AS in breast cancer.

Unraveling synonymous and deep intronic variants causing aberrant splicing in two genetically undiagnosed epilepsy families

Background

Variants identified through parent–child trio-WES yield up to 28–55% positive diagnostic rate across a variety of Mendelian disorders, there remain numerous patients who do not receive a genetic diagnosis. Studies showed that some aberrant splicing variants, which are either not readily detectable by WES or could be miss-interpreted by regular detecting pipelines, are highly relevant to human diseases.

Methods

We retrospectively investigated the negative molecular diagnostics through trio-WES for 15 genetically undiagnosed patients whose clinical manifestations were highly suspected to be genetic disorders with well-established genotype–phenotype relationships. We scrutinized the synonymous variants from WES data and Sanger sequenced the suspected intronic region for deep intronic variants. The functional consequences of variants were analyzed by in vitro minigene experiments.

Results

Here, we report two abnormal splicing events, one of which caused exon truncating due to the activation of cryptic splicing site by a synonymous variant; the other caused partial intron retention due to the generation of splicing sites by a deep intronic variant.

Conclusions

We suggest that, despite initial negative genetic test results in clinically highly suspected genetic diseases, the combination of predictive bioinformatics and functional analysis should be considered to unveil the genetic etiology of undiagnosed rare diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01008-8.

RBM10: Structure, functions, and associated diseases

RBM10 is a nuclear RNA-binding protein (RBP) that regulates the alternative splicing of primary transcripts. Recently, research on RBM10 has become increasingly active owing to its clinical importance, as indicated by studies on RBM0 mutations that cause TARP syndrome, an X-linked congenital pleiotropic developmental anomaly, and various cancers such as lung adenocarcinoma in adults. Herein, the molecular biology of RBM10 and its significance in medicine are reviewed, focusing on the gene and protein structures of RBM10, its cell biology, molecular functions and regulation, relationship with the paralogous protein RBM5, and the mutations of RBM10 and their associated diseases. Finally, the challenges in future studies of RBM10 are discussed in the concluding remarks.

Alternative Splicing of Human Telomerase Reverse Transcriptase (hTERT) and Its Implications in Physiological and Pathological Processes

Alternative splicing (AS) of human telomerase catalytic subunit (hTERT, human telomerase reverse transcriptase) pre-mRNA strongly regulates telomerase activity. Several proteins can regulate AS in a cell type-specific manner and determine the functions of cells. In addition to being involved in telomerase activity regulation, AS provides cells with different splice variants that may have alternative biological activities. The modulation of telomerase activity through the induction of hTERT AS is involved in the development of different cancer types and embryos, and the differentiation of stem cells. Regulatory T cells may suppress the proliferation of target human and murine T and B lymphocytes and NK cells in a contact-independent manner involving activation of TERT AS. This review focuses on the mechanism of regulation of hTERT pre-mRNA AS and the involvement of splice variants in physiological and pathological processes.

Targeting natural splicing plasticity of APOBEC3B restricts its expression and mutagenic activity

APOBEC3A (A3A) and APOBEC3B (A3B) enzymes drive APOBEC-mediated mutagenesis. Identification of factors affecting the activity of these enzymes could help modulate mutagenesis and associated clinical outcomes. Here, we show that canonical and alternatively spliced A3A and A3B isoforms produce corresponding mutagenic and non-mutagenic enzymes. Increased expression of the mutagenic A3B isoform predicted shorter progression-free survival in bladder cancer. We demonstrate that the production of mutagenic vs. non-mutagenic A3B protein isoforms was considerably affected by inclusion/skipping of exon 5 in A3B. Furthermore, exon 5 skipping, resulting in lower levels of mutagenic A3B enzyme, could be increased in vitro. Specifically, we showed the effects of treatment with an SF3B1 inhibitor affecting spliceosome interaction with a branch point site in intron 4, or with splice-switching oligonucleotides targeting exon 5 of A3B. Our results underscore the clinical role of A3B and implicate alternative splicing of A3B as a mechanism that could be targeted to restrict APOBEC-mediated mutagenesis.

Biological impact of mutually exclusive exon switching

Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.

Alternative AKT2 splicing produces protein lacking the hydrophobic motif regulatory region

Three AKT serine/threonine kinase isoforms (AKT1/AKT2/AKT3) mediate proliferation, metabolism, differentiation and anti-apoptotic signals. AKT isoforms are activated downstream of PI3-kinase and also by PI3-kinase independent mechanisms. Mutations in the lipid phosphatase PTEN and PI3-kinase that increase PIP3 levels increase AKT signaling in a large proportion of human cancers. AKT and other AGC kinases possess a regulatory mechanism that relies on a conserved hydrophobic motif (HM) C-terminal to the catalytic core. In AKT, the HM is contiguous to the serine 473 and two other newly discovered (serine 477 and tyrosine 479) regulatory phosphorylation sites. In AKT genes, this regulatory HM region is encoded in the final exon. We identified a splice variant of AKT2 (AKT2-13a), which contains an alternative final exon and lacks the HM regulatory site. We validated the presence of mRNA for this AKT2-13a splice variant in different tissues, and the presence of AKT2-13a protein in extracts from HEK293 cells. When overexpressed in HEK293 cells, AKT2-13a is phosphorylated at the activation loop and at the zipper/turn motif phosphorylation sites but has reduced specific activity. Analysis of the human transcriptome corresponding to other AGC kinases revealed that all three AKT isoforms express alternative transcripts lacking the HM regulatory motif, which was not the case for SGK1-3, S6K1-2, and classical, novel and atypical PKC isoforms. The transcripts of splice variants of Akt1-3 excluding the HM regulatory region could lead to expression of deregulated forms of AKT.

Identification of Survival-Associated Alternative Splicing Signatures in Lung Squamous Cell Carcinoma

Purpose: Alternative splicing (AS) is a post-transcriptional process that plays a significant role in enhancing the diversity of transcription and protein. Accumulating evidences have demonstrated that dysregulation of AS is associated with oncogenic processes. However, AS signature specifically in lung squamous cell carcinoma (LUSC) remains unknown. This study aimed to evaluate the prognostic values of AS events in LUSC patients.

Methods: The RNA-seq data, AS events data and corresponding clinical information were obtained from The Cancer Genome Atlas (TCGA) database. Univariate Cox regression analysis was performed to identify survival-related AS events and survival-related parent genes were subjected to Gene Ontology enrichment analysis and gene network analysis. The least absolute shrinkage and selection operator (LASSO) method and multivariate Cox regression analysis were used to construct prognostic prediction models, and their predictive values were assessed by Kaplan-Meier analysis and receiver operating characteristic (ROC) curves. Then a nomogram was established to predict the survival of LUSC patients. And the interaction network of splicing factors (SFs) and survival-related AS events was constructed by Spearman correlation analysis and visualized by Cytoscape.

Results: Totally, 467 LUSC patients were included in this study and 1,991 survival-related AS events within 1,433 genes were identified. SMAD4, FOS, POLR2L, and RNPS1 were the hub genes in the gene interaction network. Eight prognostic prediction models (seven types of AS and all AS) were constructed and all exhibited high efficiency in distinguishing good or poor survival of LUSC patients. The final integrated prediction model including all types of AS events exhibited the best prognostic power with the maximum AUC values of 0.778, 0.816, 0.814 in 1, 3, 5 years ROC curves, respectively. Meanwhile, the nomogram performed well in predicting the 1-, 3-, and 5-year survival of LUSC patients. In addition, the SF-AS regulatory network uncovered a significant correlation between SFs and survival-related AS events.

Conclusion: This is the first comprehensive study to analyze the role of AS events in LUSC specifically, which improves our understanding of the prognostic value of survival-related AS events for LUSC. And these survival-related AS events might serve as novel prognostic biomarkers and drug therapeutic targets for LUSC.

Rbm24 modulates adult skeletal muscle regeneration via regulation of alternative splicing

Rationale: The adult skeletal muscle can self-repair efficiently following mechanical or pathological damage due to its remarkable regenerative capacity. However, regulatory mechanisms underlying muscle regeneration are complicated and have not been fully elucidated. Alternative splicing (AS) is a major mechanism responsible for post-transcriptional regulation. Many aberrant AS events have been identified in patients with muscular dystrophy which is accompanied by abnormal muscle regeneration. However, little is known about the correlation between AS and muscle regeneration. It has been reported that RNA binding motif protein 24 (Rbm24), a tissue-specific splicing factor, is involved in embryo myogenesis while the role of Rbm24 in adult myogenesis (also called muscle regeneration) is poorly understood. Methods: To investigate the role of Rbm24 in adult skeletal muscle, we generated Rbm24 conditional knockout mice and satellite cell-specific knockout mice. Furthermore, a cardiotoxin (CTX)-induced injury model was utilized to assess the effects of Rbm24 on skeletal muscle regeneration. Genome-wide RNA-Seq was performed to identify the changes in AS following loss of Rbm24. Results: Rbm24 knockout mice displayed abnormal regeneration 4 months after tamoxifen treatment. Using RNA-Seq, we found that Rbm24 regulated a complex network of AS events involved in multiple biological processes, including myogenesis, muscle regeneration and muscle hypertrophy. Moreover, using a CTX-induced injury model, we showed that loss of Rbm24 in skeletal muscle resulted in myogenic fusion and differentiation defects and significantly delayed muscle regeneration. Furthermore, satellite cell-specific Rbm24 knockout mice recapitulated the defects in regeneration seen in the global Rbm24 knockout mice. Importantly, we demonstrated that Rbm24 regulated AS of Mef2d, Naca, Rock2 and Lrrfip1 which are essential for myogenic differentiation and muscle regeneration. Conclusions: The present study demonstrated that Rbm24 regulates dynamic changes in AS and is essential for adult skeletal muscle regeneration.

Linking transcription, RNA polymerase II elongation and alternative splicing

Gene expression is an intricately regulated process that is at the basis of cell differentiation, the maintenance of cell identity and the cellular responses to environmental changes. Alternative splicing, the process by which multiple functionally distinct transcripts are generated from a single gene, is one of the main mechanisms that contribute to expand the coding capacity of genomes and help explain the level of complexity achieved by higher organisms. Eukaryotic transcription is subject to multiple layers of regulation both intrinsic — such as promoter structure — and dynamic, allowing the cell to respond to internal and external signals. Similarly, alternative splicing choices are affected by all of these aspects, mainly through the regulation of transcription elongation, making it a regulatory knob on a par with the regulation of gene expression levels. This review aims to recapitulate some of the history and stepping-stones that led to the paradigms held today about transcription and splicing regulation, with major focus on transcription elongation and its effect on alternative splicing.

Large-scale comparative evaluation of user-friendly tools for predicting variant-induced alterations of splicing regulatory elements

Discriminating which nucleotide variants cause disease or contribute to phenotypic traits remains a major challenge in human genetics. In theory, any intragenic variant can potentially affect RNA splicing by altering splicing regulatory elements (SREs). However, these alterations are often ignored mainly because pioneer SRE predictors have proved inefficient. Here, we report the first large-scale comparative evaluation of four user-friendly SRE-dedicated algorithms (QUEPASA, HEXplorer, SPANR, and HAL) tested both as standalone tools and in multiple combined ways based on two independent benchmark datasets adding up to >1,300 exonic variants studied at the messenger RNA level and mapping to 89 different disease-causing genes. These methods display good predictive power, based on decision thresholds derived from the receiver operating characteristics curve analyses, with QUEPASA and HAL having the best accuracies either as standalone or in combination. Still, overall there was a tight race between the four predictors, suggesting that all methods may be of use. Additionally, QUEPASA and HEXplorer may be beneficial as well for predicting variant-induced creation of pseudoexons deep within introns. Our study highlights the potential of SRE predictors as filtering tools for identifying disease-causing candidates among the plethora of variants detected by high-throughput DNA sequencing and provides guidance for their use in genomic medicine settings.

Exploring the functional impact of alternative splicing on human protein isoforms using available annotation sources

In recent years, the emphasis of scientific inquiry has shifted from whole-genome analyses to an understanding of cellular responses specific to tissue, developmental stage or environmental conditions. One of the central mechanisms underlying the diversity and adaptability of the contextual responses is alternative splicing (AS). It enables a single gene to encode multiple isoforms with distinct biological functions. However, to date, the functions of the vast majority of differentially spliced protein isoforms are not known. Integration of genomic, proteomic, functional, phenotypic and contextual information is essential for supporting isoform-based modeling and analysis. Such integrative proteogenomics approaches promise to provide insights into the functions of the alternatively spliced protein isoforms and provide high-confidence hypotheses to be validated experimentally. This manuscript provides a survey of the public databases supporting isoform-based biology. It also presents an overview of the potential global impact of AS on the human canonical gene functions, molecular interactions and cellular pathways.

RNA-Seq Analysis Reveals Localization-Associated Alternative Splicing across 13 Cell Lines

Alternative splicing, a ubiquitous phenomenon in eukaryotes, is a regulatory mechanism for the biological diversity of individual genes. Most studies have focused on the effects of alternative splicing for protein synthesis. However, the transcriptome-wide influence of alternative splicing on RNA subcellular localization has rarely been studied. By analyzing RNA-seq data obtained from subcellular fractions across 13 human cell lines, we identified 8720 switching genes between the cytoplasm and the nucleus. Consistent with previous reports, intron retention was observed to be enriched in the nuclear transcript variants. Interestingly, we found that short and structurally stable introns were positively correlated with nuclear localization. Motif analysis reveals that fourteen RNA-binding protein (RBPs) are prone to be preferentially bound with such introns. To our knowledge, this is the first transcriptome-wide study to analyze and evaluate the effect of alternative splicing on RNA subcellular localization. Our findings reveal that alternative splicing plays a promising role in regulating RNA subcellular localization.

Innovative Therapeutic and Delivery Approaches Using Nanotechnology to Correct Splicing Defects Underlying Disease

Alternative splicing of pre-mRNA contributes strongly to the diversity of cell- and tissue-specific protein expression patterns. Global transcriptome analyses have suggested that >90% of human multiexon genes are alternatively spliced. Alterations in the splicing process cause missplicing events that lead to genetic diseases and pathologies, including various neurological disorders, cancers, and muscular dystrophies. In recent decades, research has helped to elucidate the mechanisms regulating alternative splicing and, in some cases, to reveal how dysregulation of these mechanisms leads to disease. The resulting knowledge has enabled the design of novel therapeutic strategies for correction of splicing-derived pathologies. In this review, we focus primarily on therapeutic approaches targeting splicing, and we highlight nanotechnology-based gene delivery applications that address the challenges and barriers facing nucleic acid-based therapeutics.

Prognostic alternative mRNA splicing signatures and associated splicing factors in acute myeloid leukemia

The dysregulation of alternative splicing (AS) has emerged as a mechanism of acute myeloid leukemia (AML). However, the prognostic impact of AS events remains under-explored in AML. Here we report the prognostic value of AS events and associated splicing factors based on three datasets of AML patients. We defined the landscape of AS events in AML and identified 7033 AS events associated with the survival of AML patients. Based on these events, we further developed a composite 15 AS event-based prognostic signature, which was independent of the cytogenetic risk stratification and patient age, and showed a better performance than known gene expression signatures. More importantly, our new signature markedly improved the European LeukemiaNet (ELN) risk classification, indicating a broad applicability in the clinical management of AML. Furthermore, the splicing-regulatory network established the correlations between prognostic AS events and associated splicing factors. The finding was validated by CRISPR-based data, which indicated that the increased expression of RBM39 contributed to the higher exon inclusion of SETD5 and conferred a poor outcome. Together, AS events may serve as a novel assortment of prognosticators for AML and could refine the ELN risk stratification. The splicing regulatory network provides clues regarding the splicing factor-mediated mechanisms of AML.

FMR1 locus isoforms: potential biomarker candidates in fragile X-associated tremor/ataxia syndrome (FXTAS)

Fragile X associated tremor/ataxia syndrome (FXTAS) is a late adult-onset neurodegenerative disorder that affects movement and cognition in male and female carriers of a premutation allele of 55–200 CGG repeats in the Fragile X mental retardation (FMR1) gene. It is currently unknown if and when an individual carrier of a premutation allele will develop FXTAS, as clinical assessment fails to identify carriers at risk before significant neurological symptoms are evident. The primary objective of this study was to investigate the alternative splicing landscape at the FMR1 locus in conjunction with brain measures in male individuals with a premutation allele enrolled in a very first longitudinal study, compared to age-matched healthy male controls, with the purpose of identifying biomarkers for early diagnosis, disease prediction and, a progression of FXTAS. Our findings indicate that increased expression of FMR1 mRNA isoforms, including Iso4/4b, Iso10/10b, as well as of the ASFMR1 mRNAs Iso131bp, are present in premutation carriers as compared to non-carrier healthy controls. More specifically, we observed a higher expression of Iso4/4b and Iso10/10b, which encode for truncated proteins, only in those premutation carriers who developed symptoms of FXTAS over time as compared to non-carrier healthy controls, suggesting a potential role in the development of the disorder. In addition, we found a significant association of these molecular changes with various measurements of brain morphology, including the middle cerebellar peduncle (MCP), superior cerebellar peduncle (SCP), pons, and midbrain, indicating their potential contribution to the pathogenesis of FXTAS. Interestingly, the high expression levels of Iso4/4b observed both at visit 1 and visit 2 and found to be associated with a decrease in mean MCP width only in those individuals who developed FXTAS over time, suggests their role as potential biomarkers for early diagnosis of FXTAS.

HITS-CLIP reveals sex-differential RNA binding and alterative splicing regulation of SRm160 in Drosophila

Serine/arginine (SR)-rich proteins are critical for the regulation of alternative splicing (AS), which generates multiple mRNA isoforms from one gene and provides protein diversity for cell differentiation and tissue development. Genetic evidence suggests that Drosophila genital-specific overexpression of SR-related nuclear matrix protein of 160 kDa (SRm160), an SR protein with a PWI RNA-binding motif, causes defective development only in male flies and results in abnormal male genital structures and abnormal testis. However, the molecular characterization of SRm160 is limited. Using the high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP) method in two sex-specific embryonic cell lines, S2 from the male and Kc from the female, we first identified the genome-wide RNA-binding characteristics of SRm160, which preferred binding to the exonic tri-nucleotide repeats GCA and AAC. We then validated this binding through both in vitro gel-shift assay and in vivo splicing of minigenes and found that SRm160 level affects AS of many transcripts. Furthermore, we identified 492 differential binding sites (DBS) of SRm160 varying between the two sex-specific cell lines. Among these DBS-containing genes, splicing factors were highly enriched, including transformer, a key regulator in the sex determination cascade. Analyses of fly mutants demonstrated that the SRm160 level affects AS isoforms of transformer. These findings shed crucial light on SRm160's RNA-binding specificity and regulation of AS in Drosophila sex determination and development.

tappAS: a comprehensive computational framework for the analysis of the functional impact of differential splicing

Recent advances in long-read sequencing solve inaccuracies in alternative transcript identification of full-length transcripts in short-read RNA-Seq data, which encourages the development of methods for isoform-centered functional analysis. Here, we present tappAS, the first framework to enable a comprehensive Functional Iso-Transcriptomics (FIT) analysis, which is effective at revealing the functional impact of context-specific post-transcriptional regulation. tappAS uses isoform-resolved annotation of coding and non-coding functional domains, motifs, and sites, in combination with novel analysis methods to interrogate different aspects of the functional readout of transcript variants and isoform regulation. tappAS software and documentation are available at https://app.tappas.org.

Alcohol Causes Lasting Differential Transcription in Drosophila Mushroom Body Neurons

Repeated alcohol experiences can produce long-lasting memories for sensory cues associated with intoxication. These memories can problematically trigger relapse in individuals recovering from alcohol use disorder (AUD). The molecular mechanisms by which ethanol changes memories to become long-lasting and inflexible remain unclear. New methods to analyze gene expression within precise neuronal cell types can provide further insight toward AUD prevention and treatment. Here, we used genetic tools in Drosophila melanogaster to investigate the lasting consequences of ethanol on transcription in memory-encoding neurons. Drosophila rely on mushroom body (MB) neurons to make associative memories, including memories of ethanol-associated sensory cues. Differential expression analyses revealed that distinct transcripts, but not genes, in the MB were associated with experiencing ethanol alone compared to forming a memory of an odor cue associated with ethanol. Adult MB-specific knockdown of spliceosome-associated proteins demonstrated the necessity of RNA-processing in ethanol memory formation. These findings highlight the dynamic, context-specific regulation of transcription in cue-encoding neurons, and the lasting effect of ethanol on transcript usage during memory formation.

Srsf7 Establishes the Juvenile Transcriptome through Age-Dependent Alternative Splicing in Mice

The juvenile phase is characterized by continuously progressing physiological processes such as growth and maturation, which are accompanied by transitions in gene expression. The contribution of transcriptome dynamics to the establishment of juvenile properties remains unclear. Here, we investigated alternative splicing (AS) events in postnatal growth and elucidated the landscape of age-dependent alternative splicing (ADAS) in C57BL/6 mice. Our analysis of ADAS in the cerebral cortex, cardiomyocytes, and hepatocytes revealed numerous juvenile-specific splicing isoforms that shape the juvenile transcriptome, which in turn functions as a basis for the highly anabolic status of juvenile cells. Mechanistically, the juvenile-expressed splicing factor Srsf7 mediates ADAS, as exemplified by switching from juvenile to adult forms of anabolism-associated genes Eif4a2 and Rbm7. Suppression of Srsf7 results in “fast-forwarding” of this transcriptome transition, causing impaired anabolism and growth in mice. Thus, juvenile-specific AS is indispensable for the anabolic state of juveniles and differentiates juveniles from adults.

•

Age-dependent alternative splicing (ADAS) was determined in mice

•

Srsf7 depletion causes loss of cellular juvenescence

•

Srsf7 mutation causes a shift from juvenile to adult-type transcriptome

•

Srsf7 promotes juvenile growth and anabolism through ADAS

Identification of prognostic splicing factors and exploration of their potential regulatory mechanisms in pancreatic adenocarcinoma

Pancreatic adenocarcinoma (PAAD), the most common subtype of pancreatic cancer, is a highly lethal disease. In this study, we integrated the expression profiles of splicing factors (SFs) of PAAD from RNA-sequencing data to provide a comprehensive view of the clinical significance of SFs. A prognostic index (PI) based on SFs was developed using the least absolute shrinkage and selection operator (LASSO) COX analysis. The PI exhibited excellent performance in predicting the status of overall survival of PAAD patients. We also used the percent spliced in (PSI) value obtained from SpliceSeq software to quantify different types of alternative splicing (AS). The prognostic value of AS events was explored using univariate COX and LASSO COX analyses; AS-based PIs were also proposed. The integration of prognosis-associated SFs and AS events suggested the potential regulatory mechanisms of splicing processes in PAAD. This study defined the markedly clinical significance of SFs and provided novel insight into their potential regulatory mechanisms.

Single-nucleotide polymorphism c.474G>A in the SEPT12 gene is a predisposing factor in male infertility

SEPT12 is a testis-specific gene involved in the terminal differentiation of male germ cells. SEPT12 protein is required for sperm head-tail formation and acts as a fundamental constituent of sperm tail annulus. In this study, we screened genetic variations in exons 5, 6, 7 of the SEPT12 and assessed the annulus status in teratozoospermic, globozoospermic, and patients with immotile short tail sperm. DNA sequencing was performed for 90 teratozoospermic and 30 normozoospermic individuals. Immunocytochemistry, transmission electron microscopy and western blotting were conducted to evaluate annulus status and the expression level of SEPT12 in patients with a distinct exonic variation (c.474G>A), respectively. Five polymorphisms identified within the desired regions of the SEPT12, among them c.474G>A had the potential to induce aberrant splicing results in the expression of a truncated protein. The annulus was detected in most of the spermatozoa from teratozoospermic and normozoospermic men with c.474G>A. In contrast, in the patient with short tail sperm defect carrying c.474G>A, 99% of spermatozoa were devoid of the annulus. Based on our findings there would be no association between exons 5, 6, 7 polymorphisms of the SEPT12 gene and the occurrence of mentioned disease but c.474G>A would be a predisposing factor in male infertility.

Esrp1-Regulated Splicing of Arhgef11 Isoforms Is Required for Epithelial Tight Junction Integrity

The epithelial-specific splicing regulators Esrp1 and Esrp2 are required for mammalian development, including establishment of epidermal barrier functions. However, the mechanisms by which Esrp ablation causes defects in epithelial barriers remain undefined. We determined that the ablation of Esrp1 and Esrp2 impairs epithelial tight junction (TJ) integrity through loss of the epithelial isoform of Rho GTP exchange factor Arhgef11. Arhgef11 is required for the maintenance of TJs via RhoA activation and myosin light chain (MLC) phosphorylation. Ablation or depletion of Esrp1/2 or Arhgef11 inhibits MLC phosphorylation and only the epithelial Arhgef11 isoform rescues MLC phosphorylation in Arhgef11 KO epithelial cells. Mesenchymal Arhgef11 transcripts contain a C-terminal exon that binds to PAK4 and inhibits RhoA activation byArhgef11. Deletion of the mesenchymal-specific Arhgef11 exon in Esrp1/2 KO epithelial cells using CRISPR/Cas9 restored TJ function, illustrating how splicing alterations can be mechanistically linked to disease phenotypes that result from impaired functions of splicing regulators.

Emerging roles of histone modifications and HDACs in RNA splicing

Histone modifications and RNA splicing, two seemingly unrelated gene regulatory processes, greatly increase proteome diversity and profoundly influence normal as well as pathological eukaryotic cellular functions. Like many histone modifying enzymes, histone deacetylases (HDACs) play critical roles in governing cellular behaviors and are indispensable in numerous biological processes. While the association between RNA splicing and histone modifications is beginning to be recognized, a lack of knowledge exists regarding the role of HDACs in splicing. Recent studies however, reveal that HDACs interact with spliceosomal and ribonucleoprotein complexes, actively control the acetylation states of splicing-associated histone marks and splicing factors, and thereby unexpectedly could modulate splicing. Here, we review the role of histone/protein modifications and HDACs in RNA splicing and discuss the convergence of two parallel fields, which supports the argument that HDACs, and perhaps most histone modifying enzymes, are much more versatile and far more complicated than their initially proposed functions. Analogously, an HDAC-RNA splicing connection suggests that splicing is regulated by additional upstream factors and pathways yet to be defined or not fully characterized. Some human diseases share common underlying causes of aberrant HDACs and dysregulated RNA splicing and, thus, further support the potential link between HDACs and RNA splicing.

Deep profiling and custom databases improve detection of proteoforms generated by alternative splicing

Alternative pre-mRNA splicing has long been proposed to contribute greatly to proteome complexity. However, the extent to which mature mRNA isoforms are successfully translated into protein remains controversial. Here, we used high-throughput RNA sequencing and mass spectrometry (MS)–based proteomics to better evaluate the translation of alternatively spliced mRNAs. To increase proteome coverage and improve protein quantitation, we optimized cell fractionation and sample processing steps at both the protein and peptide level. Furthermore, we generated a custom peptide database trained on analysis of RNA-seq data with MAJIQ, an algorithm optimized to detect and quantify differential and unannotated splice junction usage. We matched tandem mass spectra acquired by data-dependent acquisition (DDA) against our custom RNA-seq based database, as well as SWISS-PROT and RefSeq databases to improve identification of splicing-derived proteoforms by 28% compared with use of the SWISS-PROT database alone. Altogether, we identified peptide evidence for 554 alternate proteoforms corresponding to 274 genes. Our increased depth and detection of proteins also allowed us to track changes in the transcriptome and proteome induced by T-cell stimulation, as well as fluctuations in protein subcellular localization. In sum, our data here confirm that use of generic databases in proteomic studies underestimates the number of spliced mRNA isoforms that are translated into protein and provides a workflow that improves isoform detection in large-scale proteomic experiments.

Dissecting splicing decisions and cell-to-cell variability with designed sequence libraries

Most human genes are alternatively spliced, allowing for a large expansion of the proteome. The multitude of regulatory inputs to splicing limits the potential to infer general principles from investigating native sequences. Here, we create a rationally designed library of >32,000 splicing events to dissect the complexity of splicing regulation through systematic sequence alterations. Measuring RNA and protein splice isoforms allows us to investigate both cause and effect of splicing decisions, quantify diverse regulatory inputs and accurately predict (R² = 0.73–0.85) isoform ratios from sequence and secondary structure. By profiling individual cells, we measure the cell-to-cell variability of splicing decisions and show that it can be encoded in the DNA and influenced by regulatory inputs, opening the door for a novel, single-cell perspective on splicing regulation.

Alternative splicing is regulated by multiple mechanisms. Here the authors employed designed splice site libraries and massively parallel reporter assays to dissect the regulatory complexity and cell-to-cell variability of splicing decisions and to build accurate predictive models.