Striking circadian neuron diversity and cycling of Drosophila alternative splicing

  1. Qingqing Wang
  2. Katharine C Abruzzi  Is a corresponding author
  3. Michael Rosbash
  4. Donald C Rio  Is a corresponding author
  1. University of California, Berkeley, United States
  2. Howard Hughes Medical Institute, Brandeis University, United States
  3. Brandeis University, United States
7 figures, 1 table and 6 additional files

Figures

Figure 1 with 2 supplements
Many novel alternatively spliced junctions were identified in isolated Drosophila neurons.

Alternatively spliced junctions (AS junctions) were identified in RNA-seq libraries generated from isolated Drosophila neurons (DN1s, LNds, LNvs, and dopaminergic (TH)) as well as from a separately …

https://doi.org/10.7554/eLife.35618.003
Figure 1—source data 1

AS junctions identified in DN1, LNd, LNv, TH neurons and heads.

https://doi.org/10.7554/eLife.35618.006
Figure 1—figure supplement 1
Many AS junctions with varying lengths were identified in isolated neuron subtypes.

Alternatively spliced junctions (AS junctions) were identified in RNA-seq libraries generated from isolated Drosophila neurons (DN1s, LNds, LNvs, and dopaminergic (TH)) as well as from a separately …

https://doi.org/10.7554/eLife.35618.004
Figure 1—figure supplement 2
Transcripts with neuron subgroup specific novel ALT junctions reflect critical neuronal functions.

Gene ontology analyses of those transcripts that present neuron subgroup-specific novel ALT junctions that are not found in whole head samples for (A) LNv neurons, (B) LNd neurons, (C) DN1 neurons …

https://doi.org/10.7554/eLife.35618.005
Cry and CG10483 transcripts show neuron subtype-specific alternative splicing that is absent in the heterogeneous head sample.

Neuron subtype-specific novel alternative splicing of cry and CG10483 are shown. RNA-seq data tracks derived from the neuronal samples are shown, with arcs representing splice junctions and the …

https://doi.org/10.7554/eLife.35618.007
Figure 3 with 3 supplements
Identification of differentially spliced AS events in circadian neurons compared to non-circadian, TH neurons.

(A) AS events that were differentially spliced in the circadian neurons (DN1, LNd and LNv) compared to the dopaminergic neurons (TH) were identified using JUM, respectively. The number and type of …

https://doi.org/10.7554/eLife.35618.008
Figure 3—source data 1

Differentially spliced AS events in DN1 neurons versus the non-circadian TH neurons.

https://doi.org/10.7554/eLife.35618.012
Figure 3—source data 2

Differentially spliced AS events in LNd neurons versus the non-circadian TH neurons.

https://doi.org/10.7554/eLife.35618.013
Figure 3—source data 3

Differentially spliced AS events in LNv neurons versus the non-circadian TH neurons.

https://doi.org/10.7554/eLife.35618.014
Figure 3—source data 4

Differentially spliced AS events in DN1 neurons versus LNd neurons.

https://doi.org/10.7554/eLife.35618.015
Figure 3—source data 5

Differentially spliced AS events in LNd neurons versus LNv neurons.

https://doi.org/10.7554/eLife.35618.016
Figure 3—source data 6

Differentially spliced AS events in DN1 neurons versus LNv neurons.

https://doi.org/10.7554/eLife.35618.017
Figure 3—source data 7

Differentially expressed gene transcripts in DN1 neurons versus the non-circadian TH neurons.

https://doi.org/10.7554/eLife.35618.018
Figure 3—source data 8

Differentially expressed gene transcripts in LNd neurons versus the non-circadian TH neurons.

https://doi.org/10.7554/eLife.35618.019
Figure 3—source data 9

Differentially expressed gene transcripts in LNv neurons versus the non-circadian TH neurons.

https://doi.org/10.7554/eLife.35618.020
Figure 3—figure supplement 1
Overlap of differentially spliced AS events in circadian neurons compared to non-circadian, TH neurons in each conventionally recognized splicing pattern category.

Venn graphs showing the overlap of differentially alternatively spliced AS events among the three circadian neuron subtypes versus TH neurons, in each splicing pattern category. Most of the …

https://doi.org/10.7554/eLife.35618.009
Figure 3—figure supplement 2
Identification of differentially spliced AS events among the circadian neuron subpopulations.

AS events that were differentially spliced in each of the circadian neuron subpopulations compared to another (DN1 vs. LNd, LNd vs. LNv and DN1 vs. LNv) were identified using JUM, respectively. The …

https://doi.org/10.7554/eLife.35618.010
Figure 3—figure supplement 3
Functional comparison of differentially spliced cassette exons and non-differentially spliced cassette exons in circadian versus non-circadian TH neurons.

(A) Cross-species conservation comparison. Average PhastCon scores across 27 insect species for the set of differentially spliced cassette exons (AS) and non-differentially spliced cassette exons …

https://doi.org/10.7554/eLife.35618.011
N-syb and Shab undergo differential splicing in circadian neurons compared to non-circadian, TH neurons.

(A) In LNvs, the Neuronal synaptobrevin (N-syb) transcript is enriched for the isoform that includes an exon present in variant J. The alternatively spliced exon is marked by ‘*’. Use of this exon …

https://doi.org/10.7554/eLife.35618.021
Differentially alternatively spliced transcripts in LNvs and DN1s compared to non-circadian neurons reflect specific neuronal functions.

(A) Gene ontology analysis of those transcripts that undergo differential alternative splicing in DN1s reveals enrichments in potassium ion transport, calmodulin-dependent protein kinases, and …

https://doi.org/10.7554/eLife.35618.022
Figure 6 with 3 supplements
Identification of cycling alternative splicing in Drosophila circadian neurons.

To identify cycling time-of-day changes in alternative splicing, the frequency at which each sub-AS-junction is utilized in every profiled AS structure was examined across six timepoints taken at …

https://doi.org/10.7554/eLife.35618.023
Figure 6—source data 1

Cycling AS events in DN1, LNd, LNv and TH neurons.

https://doi.org/10.7554/eLife.35618.027
Figure 6—figure supplement 1
Predominant inclusion of exon a of the Sgg transcript in LNv and DN1 neurons.

The frequency of inclusion of exon b of the Sgg transcript in LNv (blue) and DN1 (red) neurons is shown. Low levels of exon b inclusion are found throughout the day. The x-axis shows Zeitgeber time …

https://doi.org/10.7554/eLife.35618.024
Figure 6—figure supplement 2
Identification of slo transcripts that exhibit cycling alternative splicing in Drosophila circadian neurons.

(A) Slowpoke (Slo) undergoes cycling alternative splicing in DN1 neurons. Sashimi plots highlight RNA-seq signals within the region of Slo that where exons 2a (marked by ‘*’) and exon 2b (marked by …

https://doi.org/10.7554/eLife.35618.025
Figure 6—figure supplement 3
Transcripts with neuron subgroup-specific cycling AS structures reflect crucial neuronal functions.

Gene ontology analyses of those transcripts that have neuron subgroup-specific cycling AS structures for (A) LNv neurons, (B) LNd neurons, (C) DN1 neurons, respectively.

https://doi.org/10.7554/eLife.35618.026
Figure 7 with 1 supplement
Cycling alternative splicing and mRNA expression in Drosophila circadian neurons.

(A) The phase distribution of cycling AS structures in DN1 neurons is plotted as a histogram (blue) with the phase distribution of all cycling transcripts in DN1 neurons overlaid (orange). In DN1s, …

https://doi.org/10.7554/eLife.35618.028
Figure 7—figure supplement 1
Phase distribution of cycling AS structures in LNvs and LNds.

The phase distribution of cycling AS structures in LNv and LNd neurons is plotted as a histogram (blue) with the phase distribution of all cycling transcripts overlaid (orange). In LNvs, cycling …

https://doi.org/10.7554/eLife.35618.029

Tables

Key resources table
Reagent type or resourceDesignationSource or referenceIdentifiersAdditional information
RNA-seq data(Abruzzi et al., 2017)GSE77451
RNA-seq data(Wang et al., 2016)GSE79916
Software, algorithmMorpheusN/Ahttps://software.broadinstitute.org/morpheus/
Software, algorithmBoxPlotR(Spitzer et al., 2014)http://shiny.chemgrid.org/boxplotr/
Software, algorithmSTAR 2.5.3a(Dobin et al., 2013)https://github.com/alexdobin/STAR
Software, algorithmSamtools 1.4.1(Li, 2011)http://samtools.sourceforge.net/
Software, algorithmBedtools 2.26.0N/Ahttp://bedtools.readthedocs.io/en/latest/
Software, algorithmJUM 1.3.7(Wang and Rio, 2017b)https://github.com/qqwang-berkeley/JUM
(copy archived at
https://github.com/elifesciences-publications/JUM)
Software, algorithmImageJ 1.50i(Schneider et al., 2012)https://imagej.nih.gov/ij/
Software, algorithmDavid Bioinformatics
Resources 6.8
(Huang et al., 2009b)https://david.ncifcrf.gov/home.jsp
Software, algorithmIGV 2.3.91(Robinson et al., 2011;
Thorvaldsdóttir et al., 2013)
http://software.broadinstitute.org/software/igv/

Additional files

Supplementary file 1

The distribution of splice site type (including reverse strand) in the AS junctions detected in each neuron group/sample.

The distribution of splice site type in the AS junctions (including reverse strand) detected in each neuron group/sample that are within a gene.

https://doi.org/10.7554/eLife.35618.030
Supplementary file 2

Gene ontology analyses of those transcripts that present neuron subgroup specific novel ALT junctions that are not found in whole head samples for each neuron subpopulation.

https://doi.org/10.7554/eLife.35618.031
Supplementary file 3

Gene ontology analyses of those transcripts that present neuron subgroup specific cycling AS structures.

https://doi.org/10.7554/eLife.35618.032
Supplementary file 4

Mapping statistics of the neuronal RNA-seq data to the Drosophila genome (dm3).

https://doi.org/10.7554/eLife.35618.033
Supplementary file 5

Commands for STAR mapping.

https://doi.org/10.7554/eLife.35618.034
Transparent reporting form
https://doi.org/10.7554/eLife.35618.035

Download links