Fig. 2
Differential transcript usage in tissues and castes. A Number of splicing isoforms per gene in HSAL50 and HSAL51. B Number of genes with select alternative splicing events (SE = skipped exon; ME = mutually exclusive exons; A5 = alternative 5′ splice site; A3 = alternative 3′ splice site; RI = retained intron; AF = alternative first exon; AL = alternative last exon) in HSAL50 and HSAL51. C Number of genes with differential transcript usage (adjusted P value (padj) < 10−5, maximum difference in proportion between isoforms (maxDprop) > 0.5, and fraction of replication iterations that support a positive DTU classification (rep_dtu_freq_threshold) > 0.8) between at least two tissues out of brain, ovary, fat body, retina, optic lobe, and antenna. Genes identified as differentially spliced in both HSAL50 and HSAL51 are represented in gray, while genes identified in only one annotation are in red. D Tissue-specific isoform usage for Mlf (identified as DTU gene, padj < 10−10) according to HSAL50 (left) or HSAL51 (right) annotation. Proportions were calculated for each replicate, then averaged within condition. *, p < 0.05. E, F Genome browser (E) and sashimi plot (F) view of Mlf showing tissue-specific alternative isoforms. Splice junction line widths are scaled to the number of reads spanning the splice junction and the total number of reads mapped to Mlf for each tissue. G Caste-specific isoform usage in the brain for Ilp2 (identified as DTU gene, padj < 10−10). Proportions were calculated for each replicate, then averaged within condition. *, p < 0.05. H Sashimi plot for the Ilp2 gene. Splice junction line widths are scaled to the number of reads spanning the splice junction and the total number of reads mapped to Ilp2 for each caste