Differential splicing analysis based on isoforms expression with NBSplice

Comprehensive benchmark of differential transcript usage analysis for static and dynamic conditions

RNA sequencing offers unique insights into transcriptome diversity, and a plethora of tools have been developed to analyze alternative splicing. One important task is to detect changes in the relative transcript abundance in differential transcript usage (DTU) analysis. The choice of the right analysis tool is non-trivial and depends on experimental factors such as the availability of single- or paired-end and bulk or single-cell data. To help users select the most promising tool for their task, we performed a comprehensive benchmark of DTU detection tools. We cover a wide array of experimental settings, using simulated bulk and single-cell RNA-seq data as well as real transcriptomics datasets, including time-series data. Our results suggest that DEXSeq, edgeR, and LimmaDS are better choices for paired-end data, while DSGseq and DEXSeq can be used for single-end data. In single-cell simulation settings, we showed that satuRn performs better than DTUrtle. In addition, we showed that Spycone is optimal for time series DTU/IS analysis based on the evidence provided using GO terms enrichment analysis.

Human IL-17A protein production is controlled through a PIP5K1α-dependent translational checkpoint

The cytokine interleukin-17 (IL-17) is secreted by T helper 17 (TH17) cells and is beneficial for microbial control; however, it also causes inflammation and pathological tissue remodeling in autoimmunity. Hence, TH17 cell differentiation and IL-17 production must be tightly regulated, but, to date, this has been defined only in terms of transcriptional control. Phosphatidylinositols are second messengers produced during T cell activation that transduce signals from the T cell receptor (TCR) and costimulatory receptors at the plasma membrane. Here, we found that phosphatidylinositol 4,5-bisphosphate (PIP2) was enriched in the nuclei of human TH17 cells, which depended on the kinase PIP5K1α, and that inhibition of PIP5K1α impaired IL-17A production. In contrast, nuclear PIP2 enrichment was not observed in TH1 or TH2 cells, and these cells did not require PIP5K1α for cytokine production. In T cells from people with multiple sclerosis, IL-17 production elicited by myelin basic protein was blocked by PIP5K1α inhibition. IL-17 protein was affected without altering either the abundance or stability of IL17A mRNA in TH17 cells. Instead, analysis of PIP5K1α-associating proteins revealed that PIP5K1α interacted with ARS2, a nuclear cap-binding complex scaffold protein, to facilitate its binding to IL17A mRNA and subsequent IL-17A protein production. These findings highlight a transcription-independent, translation-dependent mechanism for regulating IL-17A protein production that might be relevant to other cytokines.

Transcriptional and Mutational Profiling of B-Other Acute Lymphoblastic Leukemia for Improved Diagnostics

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is the most common cancer in children, and significant progress has been made in diagnostics and the treatment of this disease based on the subtypes of BCP-ALL. However, in a large proportion of cases (B-other), recurrent BCP-ALL-associated genomic alterations remain unidentifiable by current diagnostic procedures. In this study, we performed RNA sequencing and analyzed gene fusions, expression profiles, and mutations in diagnostic samples of 185 children with BCP-ALL. Gene expression clustering showed that a subset of B-other samples partially clusters with some of the known subgroups, particularly DUX4-positive. Mutation analysis coupled with gene expression profiling revealed the presence of distinctive BCP-ALL subgroups, characterized by the presence of mutations in known ALL driver genes, e.g., PAX5 and IKZF1. Moreover, we identified novel fusion partners of lymphoid lineage transcriptional factors ETV6, IKZF1 and PAX5. In addition, we report on low blast count detection thresholds and show that the use of EDTA tubes for sample collection does not have adverse effects on sequencing and downstream analysis. Taken together, our findings demonstrate the applicability of whole-transcriptome sequencing for personalized diagnostics in pediatric ALL, including tentative classification of the B-other cases that are difficult to diagnose using conventional methods.

satuRn: Scalable analysis of differential transcript usage for bulk and single-cell RNA-sequencing applications

Alternative splicing produces multiple functional transcripts from a single gene. Dysregulation of splicing is known to be associated with disease and as a hallmark of cancer. Existing tools for differential transcript usage (DTU) analysis either lack in performance, cannot account for complex experimental designs or do not scale to massive single-cell transcriptome sequencing (scRNA-seq) datasets. We introduce satuRn, a fast and flexible quasi-binomial generalized linear modelling framework that is on par with the best performing DTU methods from the bulk RNA-seq realm, while providing good false discovery rate control, addressing complex experimental designs, and scaling to scRNA-seq applications.

satuRn: Scalable analysis of differential transcript usage for bulk and single-cell RNA-sequencing applications

Alternative splicing produces multiple functional transcripts from a single gene. Dysregulation of splicing is known to be associated with disease and as a hallmark of cancer. Existing tools for differential transcript usage (DTU) analysis either lack in performance, cannot account for complex experimental designs or do not scale to massive scRNA-seq data. We introduce satuRn, a fast and flexible quasi-binomial generalized linear modelling framework that is on par with the best performing DTU methods from the bulk RNA-seq realm, while providing good false discovery rate control, addressing complex experimental designs and scaling to scRNA-seq applications.

Loss of function mutations in GEMIN5 cause a neurodevelopmental disorder

GEMIN5, an RNA-binding protein is essential for assembly of the survival motor neuron (SMN) protein complex and facilitates the formation of small nuclear ribonucleoproteins (snRNPs), the building blocks of spliceosomes. Here, we have identified 30 affected individuals from 22 unrelated families presenting with developmental delay, hypotonia, and cerebellar ataxia harboring biallelic variants in the GEMIN5 gene. Mutations in GEMIN5 perturb the subcellular distribution, stability, and expression of GEMIN5 protein and its interacting partners in patient iPSC-derived neurons, suggesting a potential loss-of-function mechanism. GEMIN5 mutations result in disruption of snRNP complex assembly formation in patient iPSC neurons. Furthermore, knock down of rigor mortis, the fly homolog of human GEMIN5, leads to developmental defects, motor dysfunction, and a reduced lifespan. Interestingly, we observed that GEMIN5 variants disrupt a distinct set of transcripts and pathways as compared to SMA patient neurons, suggesting different molecular pathomechanisms. These findings collectively provide evidence that pathogenic variants in GEMIN5 perturb physiological functions and result in a neurodevelopmental delay and ataxia syndrome.