Differential Alternative Splicing Research Articles

Alternative Splicing Plays a Crucial Role in the Salt Tolerance of Foxtail Millet.

Foxtail millet is an important cereal crop that is relatively sensitive to salt stress, with its yield significantly affected by such stress. Alternative splicing (AS) widely affects plant growth, development, and adaptability to stressful environments. Through RNA-seq analysis of foxtail millet under different salt treatment periods, 2078 AS events were identified, and analyses were conducted on differential gene (DEG), differential alternative splicing gene (DASG), and overlapping gene. To investigate the regulatory mechanism of AS in response to salt stress in foxtail millet, the foxtail millet AS genes SiCYP19, with two AS variants (SiCYP19-a and SiCYP19-b), were identified and cloned. Yeast overexpression experiments indicated that SiCYP19 may be linked to the response to salt stress. Subsequently, we conducted overexpression experiments of both alternative splicing variants in foxtail millet roots to validate them experimentally. The results showed that, under salt stress, both SiCYP19-a and SiCYP19-b jointly regulated the salt tolerance of foxtail millet. Specifically, overexpression of SiCYP19-b significantly increased the proline content and reduced the accumulation of reactive oxygen species (ROS) in foxtail millet, compared to that in SiCYP19-a. This shows that SiCYP19-b plays an important role in increasing the content of proline and promoting the clearance of ROS, thus improving the salt tolerance of foxtail millet.

SRSF2 promotes glioblastoma cell proliferation by inducing alternative splicing of FSP1 and inhibiting ferroptosis

Objective: To investigate the effect of serine/arginine-rich splicing factor 2 (SRSF2) on ferroptosis and its possible mechanism in glioblastoma cells. Methods: The online database of gene expression profiling interactive analysis 2 (GEPIA 2) and Chinese Glioma Genome Atlas were used to analyze the expression of SRSF2 in glioblastoma tissue and its association with patients prognosis. To validate the findings of the online databases, the pathological sections of glioblastoma and non-tumor brain tissues from Tianjin Medical University General Hospital, Tianjin, China were collected and analyzed by using immunohistochemistry. Silencing SRSF2 gene expression in glioblastoma cells by siRNA was analyzed with Western blot. The proliferation index was detected by using CCK8 assay. The rescued experiment was conducted by using expression plasmid of pcDNA3.1(+)-SRSF2. The activity of ferroptosis was assessed by using the levels of iron ions and malondialdehyde in glioblastoma cells and the changes in the ratio of glutathione to oxidized glutathione. The changes of gene expression and differential pre-mRNA alternative splicing (PMAS) induced by SRSF2 were monitored by using the third-generation sequencing technology analysis, namely Oxford nanopore technologies (ONT) sequencing analysis. Results: SRSF2 expression was higher in glioblastoma tissues than non-tumor brain tissues. Immunohistochemistry also showed a positive rate of 88.48%±4.60% in glioblastoma tissue which was much higher than the 9.97%±4.57% in non-tumor brain tissue. The expression of SRSF2 was inversely correlated with overall and disease-free disease survivals (P<0.01). The proliferation index of glioblastoma cells was significantly reduced by silencing with SRSF2 siRNA (P<0.01) and could be reversed with transfection of exogenous SRSF2. The levels of intracellulariron ions and malondialdehyde increased (P<0.05), but the glutathione/oxidized glutathione ratio and the expression of key proteins in the glutathione pathway remained unchanged (P>0.05). ONT sequencing results showed that silencing SRSF2 in glioblastoma cells could induce a significant alternative 3' splice site change on ferroptosis suppressor protein 1 (FSP1). Conclusion: SRSF2 inhibits the ferroptosis in glioblastoma cells and promotes their proliferation, which may be achieved by regulating FSP1 PMAS.

Single-cell transcriptomic analysis of hematopoietic progenitor cells from patients with systemic lupus erythematosus reveals interferon-inducible reprogramming in early progenitors.

Immune cells that contribute to the pathogenesis of systemic lupus erythematosus (SLE) derive from adult hematopoietic stem and progenitor cells (HSPCs) within the bone marrow (BM). For this reason, we reasoned that fundamental abnormalities in SLE can be traced to a BM-derived HSPC inflammatory signature. BM samples from four SLE patients, six healthy controls, and two umbilical cord blood (CB) samples were used. CD34+ cells were isolated from BM and CB samples, and single-cell RNA-sequencing was performed. A total of 426 cells and 24,473 genes were used in the analysis. Clustering analysis resulted in seven distinct clusters of cell types. Mutually exclusive markers, which were characteristic of each cell type, were identified. We identified three HSPC subpopulations, one of which consisted of proliferating cells (MKI67 expressing cells), one T-like, one B-like, and two myeloid-like progenitor subpopulations. Differential expression analysis revealed i) cell cycle-associated signatures, in healthy BM of HSPC clusters 3 and 4 when compared with CB, and ii) interferon (IFN) signatures in SLE BM of HSPC clusters 3 and 4 and myeloid-like progenitor cluster 5 when compared with healthy controls. The IFN signature in SLE appeared to be deregulated following TF regulatory network analysis and differential alternative splicing analysis between SLE and healthy controls in HSPC subpopulations. This study revealed both quantitative-as evidenced by decreased numbers of non-proliferating early progenitors-and qualitative differences-characterized by an IFN signature in SLE, which is known to drive loss of function and depletion of HSPCs. Chronic IFN exposure affects early hematopoietic progenitors in SLE, which may account for the immune aberrancies and the cytopenias in SLE.

Lactylome analyses suggest systematic lysine-lactylated substrates in oral squamous cell carcinoma under normoxia and hypoxia

Cancer cells tend to live in hypoxic environment characterized by enhanced glycolysis and accumulation of lactate. Intracellular lactate is shown to drive a novel type of post-translational modification (PTM), lysine lactylation (Kla). Kla has been confirmed to affect the malignant progression of tumors such as hepatocellular carcinoma (HCC) and colon cancer, whereas the global lactylomic profiling of oral squamous cell carcinoma (OSCC) is unclear. Here, the integrative lactylome and proteome analyses by using liquid chromatography-tandem mass spectrometry (LC-MS/MS) identified 1011 Kla sites within 532 proteins and 1197 Kla sites within 608 proteins in SCC25 cells under normoxic and hypoxic environments, respectively. Among these lactylated proteins, histones accounted for only a small fraction, suggesting the presence of Kla modification of OSCC in a large number of non-histone proteins. Notably, Kla preferred to enrich in spliceosome, ribosome and glycolysis/gluconeogenesis pathway in both normoxic and hypoxic cultures. Compared with normoxia, 589 differential proteins with 898 differentially lactylated sites were detected under hypoxia, which were mainly associated with the glycolysis/gluconeogenesis pathway by KEGG analysis. Importantly, we verified the presence of lactylation modification in the spliceosomal proteins hnRNPA1, SF3A1, hnRNPU and SLU7, as well as in glycolytic enzyme PFKP. In addition, the differential alternative splicing analysis described the divergence of pre-mRNA splicing patterns in the presence or absence of sodium lactate and at different oxygen concentrations. Finally, a negative correlation between tissue Kla levels and the prognosis of OSCC patients was revealed by immunohistochemistry. Our study is the first report to elucidate the lactylome and its biological function in OSCC, which deepens our understanding of the mechanisms underlying OSCC progression and provides a novel strategy for targeted therapy for OSCC.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera

BackgroundAlternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio’s SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera.ResultsWe identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA’s role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis.ConclusionsThis study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

Sex Differences in Cochlear Transcriptomes in Horseshoe Bats.

Sexual dimorphism of calls is common in animals, whereas studies on the molecular basis underlying this phenotypic variation are still scarce. In this study, we used comparative transcriptomics of cochlea to investigate the sex-related difference in gene expression and alternative splicing in four Rhinolophus taxa. Based on 31 cochlear transcriptomes, we performed differential gene expression (DGE) and alternative splicing (AS) analyses between the sexes in each taxon. Consistent with the degree of difference in the echolocation pulse frequency between the sexes across the four taxa, we identified the largest number of differentially expressed genes (DEGs) and alternatively spliced genes (ASGs) in R. sinicus. However, we also detected multiple DEGs and ASGs in taxa without sexual differences in echolocation pulse frequency, suggesting that these genes might be related to other parameters of echolocation pulse rather than the frequency component. Some DEGs and ASGs are related to hearing loss or deafness genes in human or mice and they can be considered to be candidates associated with the sexual differences of echolocation pulse in bats. We also detected more than the expected overlap of DEGs and ASGs in two taxa. Overall, our current study supports the important roles of both DGE and AS in generating or maintaining sexual differences in animals.

Renal Endothelial Single-Cell Transcriptomics Reveals Spatiotemporal Regulation and Divergent Roles of Differential Gene Transcription and Alternative Splicing in Murine Diabetic Nephropathy.

Endothelial cell (EC) injury is a crucial contributor to the progression of diabetic kidney disease (DKD), but the specific EC populations and mechanisms involved remain elusive. Kidney ECs (n = 5464) were collected at three timepoints from diabetic BTBRob/ob mice and non-diabetic littermates. Their heterogeneity, transcriptional changes, and alternative splicing during DKD progression were mapped using SmartSeq2 single-cell RNA sequencing (scRNAseq) and elucidated through pathway, network, and gene ontology enrichment analyses. We identified 13 distinct transcriptional EC phenotypes corresponding to different kidney vessel subtypes, confirmed through in situ hybridization and immunofluorescence. EC subtypes along nephrons displayed extensive zonation related to their functions. Differential gene expression analyses in peritubular and glomerular ECs in DKD underlined the regulation of DKD-relevant pathways including EIF2 signaling, oxidative phosphorylation, and IGF1 signaling. Importantly, this revealed the differential alteration of these pathways between the two EC subtypes and changes during disease progression. Furthermore, glomerular and peritubular ECs also displayed aberrant and dynamic alterations in alternative splicing (AS), which is strongly associated with DNA repair. Strikingly, genes displaying differential transcription or alternative splicing participate in divergent biological processes. Our study reveals the spatiotemporal regulation of gene transcription and AS linked to DKD progression, providing insight into pathomechanisms and clues to novel therapeutic targets for DKD treatment.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

MRNA-Seq of testis and liver tissues reveals a testis-specific gene and alternative splicing associated with hybrid male sterility in dzo.

Alternative splicing (AS) plays an important role in the co-transcription and post-transcriptional regulation of gene expression during mammalian spermatogenesis. The dzo is the male F1 offspring of an interspecific hybrid between a domestic bull (Bos taurus ♂) and a yak (Bos grunniens ♀) which exhibits male sterility. This study aimed to identify the testis-specific genes and AS associated with hybrid male sterility in dzo. The iDEP90 program and rMATS software were used to identify the differentially expressed genes (DEG) and differential alternative splicing genes (DSG) based on RNA-seq data from the liver (n=9) and testis (n=6) tissues of domestic cattle, yak, and dzo. Splicing factors (SF) were obtained from the AmiGO2 and the NCBI databases, and Pearson correlation analysis was performed on the differentially expressed SFs and DSGs. We focused on the testis-specific DEGs and DSGs between dzo and cattle and yak. Among the top 3000 genes with the most significant variations between these 15 samples, a large number of genes showed testis-specific expression involved with spermatogenesis. Cluster analysis showed that the expression levels of these testis-specific genes were dysregulated during mitosis with a burst downregulation during the pachynema spermatocyte stage. The occurrence of AS events in the testis was about 2.5 fold greater than in the liver, with exon skipping being the major AS event (81.89~82.73%). A total of 74 DSGs were specifically expressed in the testis and were significantly enriched during meiosis I, synapsis, and in the piRNA biosynthesis pathways. Notably, STAG3 and DDX4 were of the exon skipping type, and DMC1 was a mutually exclusive exon. A total of 36 SFs were significantly different in dzo testis, compared with cattle and yak. DDX4, SUGP1, and EFTUD2 were potential SFs leading to abnormal AS of testis-specific genes in dzo. These results show that AS of testis-specific genes can affect synapsis and the piRNA biosynthetic processes in dzo, which may be important factors associated with hybrid male sterility in dzo.

Abstract 5105: Aberrant CLK1 splicing is a key dependency factor in pediatric high-grade gliomas

Abstract Pediatric brain cancer is the leading cause of disease-related mortality in children. With a paucity of therapeutically targetable alterations, there is a pressing need to identify effective strategies for many of these tumors. In this study, we systematically examine aberrant alternative splicing in pediatric brain tumors from the Open Pediatric Brain Tumor Atlas, and additional tumors from the Children's Brain Tumor Network and the Pacific Pediatric Neuro-Oncology Consortium (n = 1831). Notably, high-grade gliomas (HGGs), including those with and without H3 histone mutations, were enriched in the upper quartile of all tumors with high splicing burden (SB). HGGs demonstrated the greatest SB heterogeneity of all tumor types and we identified a total of 19,229 aberrant splicing events in 8,577 genes in midline tumors compared to non-tumor brainstem controls (ΔPSI &amp;gt;= |.20|, p &amp;lt; 0.05). We next integrated protein annotations from the UniProt Knowledgebase to assess and prioritize differential alternative splicing events affecting protein function. This approach led to the discovery of 3,754 splice variants in 5,036 genes resulting in the gain or ablation of functional sites. The protein kinase CDC-like kinase 1 (CLK1), a recognized master splicing factor and cell-cycle modulator, exhibited aberrant splicing of exon 4, whereby inclusion promoted the gain of two known phosphorylation sites. Additionally, we observed a strong correlation between CLK1 splicing, its own RNA expression and proteomic activity of the Serine/Arginine-rich (SR) family of splicing factors. The distribution of exon 4 splicing was significantly different (p = 0.019) between tumors with low vs high splicing burden; tumors with high splicing burden had increased exon 4 inclusion. Further, we performed experimental modulation of CLK1 exon 4 aberrant splicing in the brain tumor cell line KNS42 and demonstrated that forced exon 4 skipping leads to reduced CLK1 protein abundance and significantly reduces cell proliferation (p &amp;lt; 0.01 vs. control). Following bulk RNA-Seq analysis of these cell lines, we identified that forced CLK1 exon 4 skipping resulted in 4,518 dysregulated genes enriched for known cancer pathways such as mitotic spindle, E2F targets, and G2M checkpoint. Our data suggest a role for CLK1 in driving transcriptional dysregulation in pediatric high-grade gliomas and uncover a potential therapeutic vulnerability. Citation Format: Ammar S. Naqvi, Priyanka Seghal, Ryan Corbett, Komal Rathi, Brian Ennis, Bo Zhang, Miguel Brown, Daniel Miller, Adam Kraya, Joseph Dybas, Katharina Hayer, Shehbeel Arif, Zhuangzhuang Geng, Antonia Chroni, Aditya Lahiri, Karina Conkrite, Madison Hollawell, Sharon Diskin, Peter Madsen, Jessica B. Foster, Mateusz P. Koptyra, Andrei Thomas-Tikhonenko, Adam Resnick, Phillip J. Storm, Jo Lynne Rokita. Aberrant CLK1 splicing is a key dependency factor in pediatric high-grade gliomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 5105.

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

Peanut (Arachis hypogaea L.), also called groundnut is an important oil and cash crop grown widely in the world. The annual global production of groundnuts has increased to approximately 50 million tons, which provides a rich source of vegetable oils and proteins for humans. Low temperature (non-freezing) is one of the major factors restricting peanut growth, yield, and geographic distribution. Since the complexity of cold-resistance trait, the molecular mechanism of cold tolerance and related gene networks were largely unknown in peanut. In this study, comparative transcriptomic analysis of two peanut cultivars (SLH vs. ZH12) with differential cold tolerance under low temperature (10°C) was performed using Oxford Nanopore Technology (ONT) platform. As a result, we identified 8,949 novel gene loci and 95,291 new/novel isoforms compared with the reference database. More differentially expressed genes (DEGs) were discovered in cold-sensitive cultivar (ZH12) than cold-tolerant cultivar (SLH), while more alternative splicing events were found in SLH compared to ZH12. Gene Ontology (GO) analyses of the common DEGs showed that the "response to stress", "chloroplast part", and "transcription factor activity" were the most enriched GO terms, indicating that photosynthesis process and transcription factors play crucial roles in cold stress response in peanut. We also detected a total of 708 differential alternative splicing genes (DASGs) under cold stress compared to normal condition. Intron retention (IR) and exon skipping (ES) were the most prevalent alternative splicing (AS) events. In total, 4,993 transcription factors and 292 splicing factors were detected, many of them had differential expression levels and/or underwent AS events in response to cold stress. Overexpression of two candidate genes (encoding trehalose-6-phosphatephosphatases, AhTPPs) in yeast improves cold tolerance. This study not only provides valuable resources for the study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset.

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity.This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

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

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

Visfatin Affects the Transcriptome of Porcine Luteal Cells during Early Pregnancy.

Visfatin/NAMPT (VIS), the hormone exerting a pleiotropic effect, is also perceived as an important factor in the regulation of reproductive processes and pregnancy maintenance. Previous studies confirmed its involvement in the control of porcine pituitary and ovary function. In this study, we hypothesized that VIS may affect the global transcriptome of luteal cells and thus regulate the functioning of the ovaries. Illumina's NovaSeq 6000 RNA sequencing was performed to investigate the differentially expressed genes (DEGs) and long non-coding RNAs (DELs) as well as the occurrence of differential alternative splicing events (DASs) in the porcine luteal cells exposed to VIS (100 ng/mL) during the implantation period. The obtained results revealed 170 DEGs (99 up- and 71 downregulated) assigned to 45 functional annotations. Moreover, we revealed 40 DELs, of which 3 were known and 37 were described for the first time. We identified 169 DASs events. The obtained results confirmed a significant effect of VIS on the transcriptome and spliceosome of luteal cells, including the genes involved in the processes crucial for successful implantation and pregnancy maintenance as angiogenesis, steroidogenesis, inflammation, cell development, migration, and proliferation.

Alternative Splicing Is a Major Factor Shaping Transcriptome Diversity in Mild and Severe Chronic Obstructive Pulmonary Disease.

The role of alternative splicing in chronic obstructive pulmonary disease (COPD) is still largely unknown. We aimed to investigate the differences in alternatively splicing events between patients with mild-to-moderate and severe COPD compared with non-COPD control subjects and to identify splicing factors associated with aberrant alternative splicing in COPD. For this purpose, we performed genome-wide RNA-sequencing analysis of bronchial brushings from 23 patients with mild-to-moderate COPD, 121 with severe COPD, and 23 non-COPD control subjects. We found a significant difference in the frequency of alternative splicing events in patients with mild-to-moderate and severe COPD compared with non-COPD control subjects. There were from two to eight times (depending on event type) more differential alternative splicing events in the severe than in the mild-to-moderate stage. The severe COPD samples showed less intron retention and more exon skipping. It is interesting that the transcript levels of the top 10 differentially expressed splicing factors were significantly correlated with the percentage of many alternatively spliced transcripts in severe COPD. The aberrant alternative splicing in severe COPD was predicted to increase the overall protein-coding capacity of gene products. In conclusion, we observed large and significant differences in alternative splicing between bronchial samples of patients with COPD and control subjects, with more events observed in severe than in mild-to-moderate COPD. The changes in the expression of several splicing factors correlated with prevalence of alternative splicing in severe COPD. Alternative splicing can indirectly impact gene expression by changing the relative abundance of protein-coding isoforms potentially influencing pathophysiological changes. The results provide a better understanding of COPD-related alternative splicing changes.

Differential alternative splicing landscape identifies potentially functional RNA binding proteins in early embryonic development in mammals

Alternative splicing (AS) as one of the important post-transcriptional regulatory mechanisms has been poorly studied during embryogenesis. In this study, we comprehensively collected and analyzed the transcriptome data of early embryos from human and mouse. We found that AS plays an important role in this process and predicted candidate RNA binding protein (RBP) regulators that are associated with reproductive development. The predicted RBPs such as EIF4A3, MAK16, SRSF2, and UTP23 were found to be associated with reproductive disorders. By Smart-seq2 sequencing analysis, we identified 5445 aberrant alternative splicing events in Eif4a3-knockdown embryos. These events were preferentially associated with RNA processing. In conclusion, our work on the landscape and potential function of alternative splicing events will boost further investigation of detailed mechanisms and key factors regulating mammalian early embryo development and promote the inspiration of pharmaceutical approaches for disorders in this crucial biology process.

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

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

Alternative splicing patterns of hnrnp genes in gill tissues of rainbow trout (Oncorhynchus mykiss) during salinity changes

Alternative splicing (AS) plays an important role in various physiological processes in eukaryotes, such as the stress response. However, patterns of AS events remain largely unexplored during salinity acclimation in fishes. In this study, we conducted AS analysis using RNA-seq datasets to explore splicing patterns in the gill tissues of rainbow trout exposed to altered salinity environments, ranging from 0 ‰ (T0) to 30 ‰ (T30). The results revealed 1441, 351, 483, 1051 and 1049 differentially alternatively spliced (DAS) events in 5 pairwise comparisons, including T6 vs. T0, T12 vs. T0, T18 vs. T0, T24 vs. T0, and T30 vs. T0, respectively. These DAS events were derived from 1290, 328, 444, 963 and 948 genes. Enrichment analysis indicated that these DAS genes were related to RNA splicing and processing. Among these, 14 DAS genes were identified as members of the large heterogeneous nuclear RNP (hnRNP) gene family. Alternative 3′ splice site (A3SS), exon skipping (SE) and intron retention (RI) events resulted in the fragmentation or even loss of the functional RNA recognition motif (RRM) domains in hnrnpa0, hnrnp1a, hnrnp1b and hnrnpc genes. The incomplete RRM domains would hinder the interactions between hnRNP genes and pre-mRNAs. It would in turn influence the splicing patterns and mRNA stability of downstream target genes in response to salinity changes. The study provides insights into salinity acclimation in gill tissues of rainbow trout and serves as a significant reference on the osmoregulation mechanisms at post-transcription regulation levels in fish.

BetAS: intuitive analysis and visualization of differential alternative splicing using beta distributions.

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

Differential expression and alternative splicing analyses of multiple tissues reveal albinism-associated genes in the Wels catfish (Silurus glanis)

Albinism is a widespread departure from a typical body colouration due to altered melanin production. The Wels catfish (Silurus glanis) is among the largest freshwater fish species in the world, and albino individuals occur both in the wild and in aquaculture. Here, we performed transcriptome-wide analysis of albino and normally pigmented S. glanis using four tissues (skin, dorsal fin, whole eye and liver) to identify genes associated with albinism by exploring patterns of differential expression (DE) and differential alternative splicing (DAS). Multi-tissue analyses revealed a large number of genes in skin (n = 1355) and fin (n = 614) tissue associated with the albino phenotype in S. glanis, while the number of DE genes in eye and liver tissues was lower (n = 188, n = 189, respectively). Several DE genes across multiple tissues were detected as the most promising candidates (e.g., hsp4, hsp90b1, raph1, uqcrfs1, adcy-family and wnt-family) potentially causally linked to the albino phenotype in Wels catfish. Moreover, our findings supported earlier observations of physiological differences between albino and normally pigmented individuals, particularly in energy metabolism and immune response. In contrast, there were only a few pigmentation-related genes observed among DAS genes (4 in skin, 2 in fin), the overlap between DAS and DE genes was low (n = 25) and did not include known pigmentation-related genes. This suggests that DAS and DE in Wels catfish are, to a large extent, independent processes, and the observed alternative splicing cases are probably not causally linked with albinism in S. glanis. This work provides the first transcriptome-wide multi-tissue insights into the albinism of Wels catfish and serves as a valuable resource for further understanding the genetic mechanisms of pigmentation in fish.