DJExpress: An Integrated Application for Differential Splicing Analysis and Visualization

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.

Expression pattern analysis and characterization of the hereditary sensory and autonomic neuropathy 2 A (HSAN2A) gene with no lysine kinase (WNK1) in human dorsal root ganglion

Inherited painless neuropathies arise due to genetic insults that either block the normal signaling of or destroy the sensory afferent neurons in the dorsal root ganglion (DRG) responsible for transducing noxious stimuli. Complete loss of these neurons leads to profound insensitivity to all sensory modalities including pain. Hereditary sensory and autonomic neuropathy type 2 (HSNAII) is a rare genetic neuropathy characterized by a progressive distal early onset sensory loss. This syndrome is caused by autosomal recessive mutations in the with-no-lysine protein kinase 1 (WNK1) serine-threonine kinase gene. Of interest, disease-associated mutations are found in the large exon, termed "HSN2," which encodes a 498 amino acid domain C-terminal to the kinase domain. These mutations lead to truncation of the HSN2-containing proteins through the addition of an early stop codon (nonsense mutation) leading to loss of the C-terminal domains of this large protein. The present study evaluates the transcripts, gene structure, and protein structure of HSN2-containing WNK1 splice variants in DRG and spinal cord in order to establish the basal expression patterns of WNK1 and HSN2-containing WNK1 splice variants using multiplex fluorescent situ hybridization. We hypothesized that these transcripts would be enriched in pain-sensing DRG neurons, and, potentially, that enrichment in nociceptive neurons was responsible for the painless phenotypes observed. However, our in-depth analyses revealed that the HSN2-WNK1 splice variants were ubiquitously expressed but were not enriched in tachykinin 1-expressing C-fiber neurons, a class of neurons with a highly nociceptive character. We subsequently identified other subpopulations of DRG neurons with higher levels of HSN2-WNK1 expression, including mechanosensory large fibers. These data are inconsistent with the hypothesis that this transcript is enriched in nociceptive fibers, and instead suggest it may be related to general axon maintenance, or that nociceptive fibers are more sensitive to the genetic insult. These findings clarify the molecular and cellular expression pattern of this painless neuropathy gene in human tissue.

TLN1 contains a cancer-associated cassette exon that alters talin-1 mechanosensitivity

Talin-1 is the core mechanosensitive adapter protein linking integrins to the cytoskeleton. The TLN1 gene is comprised of 57 exons that encode the 2,541 amino acid TLN1 protein. TLN1 was previously considered to be expressed as a single isoform. However, through differential pre-mRNA splicing analysis, we discovered a cancer-enriched, non-annotated 51-nucleotide exon in TLN1 between exons 17 and 18, which we refer to as exon 17b. TLN1 is comprised of an N-terminal FERM domain, linked to 13 force-dependent switch domains, R1-R13. Inclusion of exon 17b introduces an in-frame insertion of 17 amino acids immediately after Gln665 in the region between R1 and R2 which lowers the force required to open the R1-R2 switches potentially altering downstream mechanotransduction. Biochemical analysis of this isoform revealed enhanced vinculin binding, and cells expressing this variant show altered adhesion dynamics and motility. Finally, we showed that the TGF-β/SMAD3 signaling pathway regulates this isoform switch. Future studies will need to consider the balance of these two TLN1 isoforms.