Counting pseudoalignments to novel splicing events.

MOTIVATION: Alternative splicing (AS) of introns from pre-mRNA produces diverse sets of transcripts across cell types and tissues, but is also dysregulated in many diseases. Alignment-free computational methods have greatly accelerated the quantification of mRNA transcripts from short RNA-seq reads, but they inherently rely on a catalog of known transcripts and might miss novel, disease-specific splicing events. By contrast, alignment of reads to the genome can effectively identify novel exonic segments and introns. Event-based methods then count how many reads align to predefined features. However, an alignment is more expensive to compute and constitutes a bottleneck in many AS analysis methods. RESULTS: Here, we propose fortuna, a method that guesses novel combinations of annotated splice sites to create transcript fragments. It then pseudoaligns reads to fragments using kallisto and efficiently derives counts of the most elementary splicing units from kallisto's equivalence classes. These counts can be directly used for AS analysis or summarized to larger units as used by other widely applied methods. In experiments on synthetic and real data, fortuna was around 7× faster than traditional align and count approaches, and was able to analyze almost 300 million reads in just 15 min when using four threads. It mapped reads containing mismatches more accurately across novel junctions and found more reads supporting aberrant splicing events in patients with autism spectrum disorder than existing methods. We further used fortuna to identify novel, tissue-specific splicing events in Drosophila. AVAILABILITY AND IMPLEMENTATION: fortuna source code is available at https://github.com/canzarlab/fortuna.

Contributions of alternative splicing to muscle type development and function.

Animals possess a wide variety of muscle types that support different kinds of movements. Different muscles have distinct locations, morphologies and contractile properties, raising the question of how muscle diversity is generated during development. Normal aging processes and muscle disorders differentially affect particular muscle types, thus understanding how muscles normally develop and are maintained provides insight into alterations in disease and senescence. As muscle structure and basic developmental mechanisms are highly conserved, many important insights into disease mechanisms in humans as well as into basic principles of muscle development have come from model organisms such as Drosophila, zebrafish and mouse. While transcriptional regulation has been characterized to play an important role in myogenesis, there is a growing recognition of the contributions of alternative splicing to myogenesis and the refinement of muscle function. Here we review our current understanding of muscle type specific alternative splicing, using examples of isoforms with distinct functions from both vertebrates and Drosophila. Future exploration of the vast potential of alternative splicing to fine-tune muscle development and function will likely uncover novel mechanisms of isoform-specific regulation and a more holistic understanding of muscle development, disease and aging.

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.

Dissection of Drosophila melanogaster Flight Muscles for Omics Approaches.

Drosophila flight muscle is a powerful model to study diverse processes such as transcriptional regulation, alternative splicing, metabolism, and mechanobiology, which all influence muscle development and myofibrillogenesis. Omics data, such as those generated by mass spectrometry or deep sequencing, can provide important mechanistic insights into these biological processes. For such approaches, it is beneficial to analyze tissue-specific samples to increase both selectivity and specificity of the omics fingerprints. Here we present a protocol for dissection of fluorescent-labeled flight muscle from live pupae to generate highly enriched muscle samples for omics applications. We first describe how to dissect flight muscles at early pupal stages (<48 h after puparium formation [APF]), when the muscles are discernable by green fluorescent protein (GFP) labeling. We then describe how to dissect muscles from late pupae (>48 h APF) or adults, when muscles are distinguishable under a dissecting microscope. The accompanying video protocol will make these technically demanding dissections more widely accessible to the muscle and Drosophila research communities. For RNA applications, we assay the quantity and quality of RNA that can be isolated at different time points and with different approaches. We further show that Bruno1 (Bru1) is necessary for a temporal shift in myosin heavy chain (Mhc) splicing, demonstrating that dissected muscles can be used for mRNA-Seq, mass spectrometry, and reverse transcription polymerase chain reaction (RT-PCR) applications. This dissection protocol will help promote tissue-specific omics analyses and can be generally applied to study multiple biological aspects of myogenesis.

Conserved functions of RNA-binding proteins in muscle.

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

Adenosine A(2A) receptor up-regulates retinal wave frequency via starburst amacrine cells in the developing rat retina.

BACKGROUND: Developing retinas display retinal waves, the patterned spontaneous activity essential for circuit refinement. During the first postnatal week in rodents, retinal waves are mediated by synaptic transmission between starburst amacrine cells (SACs) and retinal ganglion cells (RGCs). The neuromodulator adenosine is essential for the generation of retinal waves. However, the cellular basis underlying adenosine's regulation of retinal waves remains elusive. Here, we investigated whether and how the adenosine A(2A) receptor (A(2A)R) regulates retinal waves and whether A(2A)R regulation of retinal waves acts via presynaptic SACs. METHODOLOGY/PRINCIPAL FINDINGS: We showed that A(2A)R was expressed in the inner plexiform layer and ganglion cell layer of the developing rat retina. Knockdown of A(2A)R decreased the frequency of spontaneous Ca²âº transients, suggesting that endogenous A(2A)R may up-regulate wave frequency. To investigate whether A(2A)R acts via presynaptic SACs, we targeted gene expression to SACs by the metabotropic glutamate receptor type II promoter. Ca²âº transient frequency was increased by expressing wild-type A(2A)R (A2AR-WT) in SACs, suggesting that A(2A)R may up-regulate retinal waves via presynaptic SACs. Subsequent patch-clamp recordings on RGCs revealed that presynaptic A(2A)R-WT increased the frequency of wave-associated postsynaptic currents (PSCs) or depolarizations compared to the control, without changing the RGC's excitability, membrane potentials, or PSC charge. These findings suggest that presynaptic A(2A)R may not affect the membrane properties of postsynaptic RGCs. In contrast, by expressing the C-terminal truncated A(2A)R mutant (A(2A)R-ΔC) in SACs, the wave frequency was reduced compared to the A(2A)R-WT, but was similar to the control, suggesting that the full-length A(2A)R in SACs is required for A(2A)R up-regulation of retinal waves. CONCLUSIONS/SIGNIFICANCE: A(2A)R up-regulates the frequency of retinal waves via presynaptic SACs, requiring its full-length protein structure. Thus, by coupling with the downstream intracellular signaling, A(2A)R may have a great capacity to modulate patterned spontaneous activity during neural circuit refinement.