A transcriptomics resource reveals a transcriptional transition during ordered sarcomere morphogenesis in flight muscle
- Max Planck Institute of Biochemistry, Germany
- Ludwig-Maximilians-Universität München, Germany
- Aix Marseille Univ, CNRS, IBDM, France
- Foshan University, China
- Aix Marseille Univ, INSERM, TAGC, France
- University of Zurich, Switzerland
Figures
Figure 1
Development of the dorsal longitudinal indirect flight muscles.
(A–H) Time-course of DLM development. (A) Myoblasts adhering to the hinge of the larval wing disc were visualised with Him-GFP-Gma (green), F-actin was stained with phalloidin (red) and nuclei with DAPI (blue). (B–H) Time-course of DLM and myotendinous junction development at 16 hr (B), 24 hr (C), 32 hr (D), 48 hr (E), 56 hr (F), 72 hr (G) and 90 hr APF (H). F-actin was stained with phalloidin (red), Shot (green) and βPS-Integrin (blue). DLMs are highlighted in dashed yellow lines in B an C. (A’’’–C’’’) Scheme of the respective developmental stages with myoblasts and muscles in red, tendon cells in blue and wing disc or pupal thorax outline in black. The length of the muscle fibers in indicated in red. For details see text. Scale bar represents 100 µm. (I) Temporal summary of known events during myogenesis (red). Samples for mRNA-Seq were collected at time points noted in black.
-
Figure 1—source data 1
This table includes the length measurements of the indirect flight muscle fibers at the indicated stage.
- https://doi.org/10.7554/eLife.34058.004
Figure 2 with 2 supplements
Verification of mRNA-Seq time-course by ‘indicator’ gene expression.
(A,D,G,J,M,P,S) Temporal expression dynamics were evaluated by clustering standard-normal mRNA-seq counts using Mfuzz. Temporal expression profiles are plotted with high membership values in warm colours (red, pink), and lower membership values in cool colours (blue, green). (B,E,H,K,N,Q,T) The profile of one ‘indicator’ gene from each cluster is shown and coloured based on the Mfuzz membership value α. (C,F,I,L,O,R,U). Protein expression and localisation dynamics were visualised by antibody staining against Twi (C) and Aret (F) or against GFP for GFP tagged fosmid reporters for Kon (I), β-Tub60D (L), βPS-Integrin (mys) (O), Unc-89 (Obscurin) (R) and Strn-Mlck (U). Images for the same protein were acquired using the same settings, and pseudo-coloured based on intensity. Note the close correlation between mRNA and protein expression dynamics. Time points are indicated by blue dots on the mRNA expression profile. Scale bars represent 20 µm.
Figure 2—figure supplement 1
Forty distinct temporal mRNA-Seq expression profiles.
Temporal expression profiles were identified by clustering standard-normal mRNA-Seq counts using Mfuzz to group genes with similar temporal expression dynamics. Expression dynamic profiles were labelled 1–40. Each plot shows the profile for each gene in the cluster, with profiles of genes with high membership values in warm colours (red, pink) and lower membership values in cool colours (blue, green).
Figure 2—figure supplement 2
Additional examples of ‘indicator’ gene expression.
(A–P) Cluster profiles and mRNA dynamics for additional indicator genes. salm from cluster 26 (A,B), rhea (Talin) from cluster 25 (D,E), fln from cluster 3 (G,H) and Mhc, Act88F, Mf from cluster 22 (J,K,M,O) are shown. The respective protein dynamics were visualised with antibodies against the respective protein or GFP fusion protein (C,F,I,L,N,P). Note that although Mhc and Act88F show significant mRNA and protein expression already at 30 hr APF, this expression dramatically increases after 30 hr APF to reach even higher protein levels at 72 hr APF. Fln and Mf mRNA and protein are only induced to high levels after 30 hr APF. Scale bars represent 20 µm.
Figure 3
Expression dynamics reveal a temporal ordering of biological processes during muscle morphogenesis.
(Top) Heat map of Mfuzz cluster core expression profiles. Standard-normal count values for all genes with a membership value α >0.8 were averaged to generate the core expression profile for each cluster. Mfuzz expression profiles fall into eight groups (coloured dendrogram leaves) based on hierarchical clustering of their temporal expression dynamics. Time points and Mfuzz clusters are labelled. Colour scale of standard-normal count values ranges from blue (stable/no expression) to red (high expression). (Bottom) GO Biological Process and user-defined gene set (marked with *) enrichments calculated with GO-Elite. Note that proliferation, development and differentiation terms are enriched at early time points, while mitochondrial and sarcomere terms are enriched at late time points. A coloured box indicates a significant enrichment of a given term in the specified cluster (see Supplementary file 2 for details).
Figure 4 with 1 supplement
A major transition in gene expression after 30 hr APF.
(A) DESeq2 principle components analysis (PCA) of all mRNA-Seq libraries. Each triangle represents a different biological replicate coloured by time point. Note that individual replicates for a given time point cluster together. PC1 divides early (≤48 hr APF) from late (≥72 hr) stages. (B) Stacked box plot of the number of significantly differentially expressed genes up-regulated (reds) or down-regulated (blues) between sequential time points with a p-value<0.05 and a log2FC (fold change) of >1 (dark), >2 (medium) or >5 (light). The large differences between myoblast to 16 hr APF reflect muscle specification. A second large shift in expression is evident between 30 and 72 hr APF. (C) Volcano plot illustrating the strong up-regulation of sarcomeric proteins (red) from 30 hr to 72 hr APF. Significantly up- or down-regulated genes are in blue (p-value<0.05 and abs(log2FC)>2). (D) Hierarchical clustering of log2 transformed DESeq2 normalized counts for all genes that are significantly up-regulated between 30 and 72 hr APF. Note that they are either strongly induced at 48 or 72 hr APF (from yellow to red), or only turned on at 48 or 72 hr APF (blue to yellow/red), suggesting a major transition in gene expression after 30 hr APF. Colour scale of log2 count values ranges from blue (not expressed) to red (highly expressed). (E) GO-Elite and user-defined gene set (marked with *) enrichments in up- (red) and down- (blue) regulated genes from 30 hr to 72 hr APF. Note the strong enrichment of mitochondrial and sarcomere terms in the up- regulated genes. (F) Pie charts showing the proportion of genes belonging to an enriched Mfuzz cluster in the sets of genes either up- or down-regulated from 30 hr to 72 hr APF. Note that a large proportion of genes up-regulated 30 hr to 72 hr belong to cluster 22, as well as Mfuzz clusters enriched for sarcomere (yellow) or mitochondrial (green) terms.
Figure 4—figure supplement 1
Additional evidence supporting a transition in gene expression between 30 and 72 hr APF.
(A) Principle component analysis (PCA) of Mfuzz clusters. PC1 largely separates clusters based on 30 to 72 hr dynamics, with significantly ‘up-regulated’ clusters on the left (red) and ‘down regulated’ clusters on the right (green) (see Figure 4F). Sarcomeric and mitochondrial gene clusters shown in yellow script. (B) Number of differentially expressed (DE) genes that are the same between time points. Total number of DE genes on bottom. Note that sequential time points share the greatest overlap. (C,D) Venn diagrams showing the strong overlap between genes coding for sarcomeric proteins, salm core genes, Mfuzz cluster 22 and genes up-regulated (C) but not down-regulated (D) between 30 and 72 hr APF.
Figure 5 with 2 supplements
Sequential phases of sarcomere morphogenesis in flight muscles.
(A) Scheme of hemi-thoraces at 34, 48, 60, 72, 80, 90 hr APF and 1 day adults (muscle in red, tendon in blue, sarcomeres in yellow) with indicated muscle fiber length. Myofibrils and sarcomeres at these time points were stained for phalloidin (F-actin, red) and Kettin (Z-disc, green). Scale bar represents 5 µm. (B,C) Tukey box and whisker plot of sarcomere length and myofibril width. Box extends from 25 to 75%, line marks median, whiskers extend to 25/75% plus 1.5 times the interquartile range. (D) Histogram of sarcomere number per myofibril. Error bars represent SD. Note the sarcomere assembly at 34 hr, followed by sarcomere addition until ~48 hr and sarcomere maturation after ~48 hr APF. (E) Cross-sections of the DLMs at 30, 48, 60, 72, 80, 90 hr APF and 1 day adult. Scale bar represents 5 µm. (F,G) Tukey box and whisker plot of myofibril diameter and myofiber cross-sectional area. Note the lack of growth in diameter or area from 30 hr to 48 hr. (H) Histogram of the number of myofibrils per myofiber. Error bars represent SD. Note that all myofibrils are already present at 30 hr APF. Tukey’s multiple comparison p-value<0.05*,. 01**,. 001***, n.s. = not significant. N > 10 for each individual time point. (I) Stills of live movies of DLMs at 30, 36, 42, 48, 60 and 72 hr APF. Scale bar represents 50 µm. For live movies see Figure 5—video 1. Stills are a time 0.0 s image (magenta) overlaid with a time +0.65 s image (green), where a perfect overlap (white) shows no movement. (J) Quantification of spontaneous contraction events per fiber per 5 min, with single twitches in blue and double twitches in green. Fibers are first contractile at 30 hr APF, reach peak contractility at 48 hr and stop all spontaneous contraction shortly after 60 hr APF.
-
Figure 5—source data 1
This table includes calculation of the number of sarcomeres per myofibril in wild type flight muscles at the indicated stage.
- https://doi.org/10.7554/eLife.34058.013
-
Figure 5—source data 2
This table includes the sarcomere length in wild type flight muscles measured at the indicated stage.
Note that image analysis is limited by pixel size and thus can result in the same number.
- https://doi.org/10.7554/eLife.34058.014
-
Figure 5—source data 3
This table includes the sarcomere/myofibril width in wild type flight muscles measured at the indicated stage.
Note that image analysis is limited by pixel size and thus can result in the same number.
- https://doi.org/10.7554/eLife.34058.015
-
Figure 5—source data 4
This table includes the myofiber cross-sectional area measured in wild type flight muscle fibers or only fibers 3 and four at the indicated stage.
- https://doi.org/10.7554/eLife.34058.016
-
Figure 5—source data 5
This table includes the number of myofibrils per myofiber calculated for either all wild type flight muscle fibers or only fibers 3 and four at the indicated stage.
- https://doi.org/10.7554/eLife.34058.017
-
Figure 5—source data 6
This table includes the myofibril diameter measured in wild type flight muscles at the indicated stage.
- https://doi.org/10.7554/eLife.34058.018
-
Figure 5—source data 7
This table includes the muscle twitch events per wild-type fiber per 5 min recorded at the indicated stage.
- https://doi.org/10.7554/eLife.34058.019
Figure 5—figure supplement 1
Expression and localisation of thin- and thick-filament structural proteins.
(A–E) Developing flight muscle myofibrils are stained with phalloidin (red) and the respective sarcomeric proteins (green). Flightin (A) cannot be detected at 30 hr, but decorates the thick filament from 48 hr APF. Unc-89/Obscurin-GFP (B) labels the M-line from 30 hr, and is markedly refined to 72 hr APF. Myosin heavy chain (C) is visible in a regular pattern at 30 hr APF. Myofilin/Mf (D) cannot be detected at 30 hr and decorates the thick filament at 48 hr APF. Likewise, IFM-specific Strn-Mlck Isoform R (E) is not expressed at 30 hr, but is detected on the thick filament at low levels at 48 hr APF. Scale bars represent 5 µm.
Figure 5—video 1
Twitching time-course in developing DLMs.
Confocal movies of spontaneous contraction (twitching) in developing DLMs at 30, 36, 42, 48, 60 and 72 hr APF. Fibers are the same as those shown in Figure 5I. Muscle attachments are visualized using a Talin-YPet fusion, which is enriched at the muscle end. Note that fibers weakly contract at 30 hr and increase in contraction intensity and frequency until 48 hr APF, but then completely stop all contraction by 72 hr APF. Scale bar represents 20 µm. Individual movie duration (in seconds) as noted.
Figure 6 with 3 supplements
salm contributes to the transition in gene expression after 30 hr APF.
(A) Volcano plot of mRNA-Seq comparison of wild-type (WT) versus salmIR IFMs at 72 hr APF. Note the significant down-regulation of genes in salmIR, especially sarcomeric protein coding genes (red). Significantly differentially expressed (DE) genes (abs(log2 FC)>2, p<0.05) are in blue. (B) Violin plot showing down-regulation of sarcomeric proteins (in red), mitochondrial proteins (in green) and members of Mfuzz cluster 22 (in yellow) in salmIR compared to wild type at 72 hr APF. Box plots indicate the median of the data as well as the first and third quartile in the box, outliers are the dots. (C) WT mRNA-Seq fold change values of all genes significantly DE from 30 hr to 72 hr APF are ordered from lowest to highest (black). The corresponding salmIR fold change is shown in yellow. Both sarcomeric (red) and mitochondrial protein coding genes (blue) are less strongly up- or down-regulated in salmIR across the 30 hr to 72 hr APF transition. Note that many genes in salmIR (yellow dots) are not as strongly induced or even repressed compared to WT (below and above the black WT line, respectively). (D) Violin plot comparing the log2FC over the 30 hr to 72 hr transition in WT (in black) and salmIR IFM (in magenta) for members of Mfuzz Cluster22. Box plots indicate the median of the data as well as the first and third quartile in the box, outliers are the dots. Note the significant decrease in induction to 72 hr APF in the salmIR sample. ***Student’s t-test p-value<0.0005. (E–H) salm is required for the induction of some but not all sarcomeric proteins. Act88F >> salmIR in the background of GFP-tagged Strn-Mlck-IsoR (E), Fln (F), Mhc (G) and Unc-89 (H). Quantitative changes in live GFP fluorescence at 90 hr APF were measured by quantitative confocal microscopy relative to standard fluorescent beads, revealing significant decreases in induction for Strn-Mlck, Fln and Mhc between wild type control (shown in yellow, Act88F-GAL4 crossed to w1118) and Act88F >> salmIR (shown in purple). Scale bar represents 5 µm. Error bars represent SEM, Student’s t-test p-value<0.05*, 0.001***, n.s. = not significant. N > 10 for each individual sample. (E’–H’’) Intensity-coded GFP fluorescence at 90 hr APF in confocal images of fixed myofibrils.
-
Figure 6—source data 1
This table includes the fiber divided by bead fluorescence intensity measurements for quantification of the indicated fosmid-GFP or UAS-GFP-Gma expression levels at the indicated stage.
- https://doi.org/10.7554/eLife.34058.025
-
Figure 6—source data 2
This table includes the fiber divided by bead fluorescence intensity measurements for various fosmid-GFPs in the Act88F >> salmIR or control Act88F-Gal4 x w- background at 90 hr APF.
- https://doi.org/10.7554/eLife.34058.026
Figure 6—figure supplement 1
salm regulates gene expression during flight muscle development.
(A–C) Volcano plots of mRNA-Seq data comparing wild-type (WT) versus salmIR IFMs at 24 hr (A), 30 hr (B) and 1d adult. Significantly differentially expressed (DE) genes (abs(log2FC)>2, p<0.05) are in blue. Note the significant down regulation of genes at all time points, particularly sarcomeric protein coding genes (red). (D) Violin plot showing down-regulation of sarcomeric proteins (in red), mitochondrial proteins (in green) and members of Mfuzz cluster 22 (in yellow) in salmIR compared to wild type at 30 hr APF. Box plots indicate the median of the data as well as the first and third quartile in the box, outliers are the dots. Note that all three gene groups are down-regulated in in salmIR. (E) Mfuzz cluster enrichment in genes that are DE between WT and salmIR. Note that the sarcomeric and mitochondrial clusters (indicated by the green and yellow dots, respectively) are all down-regulated in salmIR. Colour scale represents enrichment Z-score. (F) GO and gene set (marked with *) enrichment in genes that are DE between WT and salmIR. Note that terms down-regulated in salmIR are enriched in genes up-regulated in WT from 30 hr to 72 hr. (G) Venn diagram showing the strong overlap between sarcomeric protein coding genes, salm core genes, Mfuzz cluster 22 and genes down-regulated in salmIR at 72 hr APF.
Figure 6—figure supplement 2
Act88F-GAL4 driven knock-down of salm is efficient.
(A–D) At 24 hr APF, both Aret and Salm expression is observed in the nuclei of wild type (A) and Act88F >> salmIR (B) DLMs. At 72 hr APF, Salm protein is present in wild type (C) but cannot be detected in the nuclei of Act88F >> salmIR (D) DLMs. Scale bars represent 10 µm.
Figure 6—figure supplement 3
Expression of sarcomere proteins strongly increases from 30 hr to 72 hr APF.
(A–E) Live assay to quantify GFP fusion protein expression by normalizing GFP expression intensity in muscle to fluorescent beads (relative fluorescence on Y-axis). Gma-GFP expression driven by Mef2-GAL4 (A) shows a decrease in expression with developmental time, likely due to a decrease in expression from the Mef2 enhancer. Mhc-GFP (B) and Unc-89-GFP (C) show detectable GFP expression at 30 hr APF, but a massive increase in expression to 72 hr APF. Fln-GFP (D) is not detectable at 30 hr APF and increases expression from 48 hr APF to 72 hr APF. Strn-Mlck-GFP (E) is weakly detectable at 48 hr APF and is massively induced to 72 hr APF. Relative fluorescence was calculated as muscle intensity/bead intensity. Error bars represent SEM, Tukey’s multiple comparison p-value<0.05*, 0.001***, n.s. = not significant. N > 10 flies for each timepoint.
Figure 7 with 3 supplements
salm is required for normal sarcomere maturation and function.
(A–H) Myofibrils of Act88F-GAL4 / + (A–D) or Act88F >> salmIR (E–H). Note that salmIR DLM remains fibrillar and appears normal at 48 hr APF (E), but at 72 hr APF (F) 90 hr APF (G) and 2 day adult (H) Z-discs widen and show actin accumulations (arrowheads). (I,J) Tukey box and whisker plot of myofibril width (I) and sarcomere length (J) in Act88F-GAL4 /+ and Act88F >> salmIR (red). Tukey’s multiple comparison p-value<0.001***. N > 10 for each individual time point. Scale bars represent 5 µm. (K–O) Stills from live movies of developing DLMs at 48 and 72 hr APF in Act88F-GAL4 / + (K, M) and Act88F >> salmIR (L, N). Scale bar represents 50 µm. Coloured as in Figure 5. (O) Quantification of spontaneous contraction events per fiber per 5 min, with single twitches in blue and double twitches in green. Error bars represent SEM. salmIR fibers continue spontaneously contracting at 72 hr APF.
-
Figure 7—source data 1
This table includes the muscle twitch events per Act88F >> salmIR fiber per 5 min and the respective controls recorded at the indicated stage.
- https://doi.org/10.7554/eLife.34058.029
-
Figure 7—source data 2
This table includes the sarcomere length measured in Act88F >> salmIR flight muscles and the respective controls at the indicated stage.
Note that image analysis is limited by pixel size and thus can result in the same number.
- https://doi.org/10.7554/eLife.34058.030
-
Figure 7—source data 3
This tables includes the sarcomere/myofibril width measured in Act88F >> salmIR flight muscles and the respective controls at the indicated stage.
Note that image analysis is limited by pixel size and thus can result in the same number.
- https://doi.org/10.7554/eLife.34058.031
Figure 7—figure supplement 1
Act88F >> salmIR and Strn-Mlck mutant flies are flightless and IFM fibers rupture in adult flies.
(A) Adult males from Act88F-GAL4 crossed to two independent salmIR hairpins as well as StrnM102893 are flightless. (B–J) At 90 hr APF, Act88F-GAL4 / + (B–D), Act88F >> salmIR (E–G) and StrnM102893 (H–J) all have six intact fibers. At 1 day and 2 days after eclosion, both Act88F >> salmIR (F,G) and StrnM102893 (I,J) show muscle tearing and atrophy. Scale bar represents 200 µm.
Figure 7—video 1
Twitching in developing Act88F-GAL4 /+ and Act88F >> salmIR DLMs at 48 hr APF.
Confocal movies of spontaneous contraction in Act88F-GAL4 / + (control) and Act88F >> salmIR DLMs at 48 hr APF. Muscle attachments are visualized using a Talin-YPet fusion. salmIR fibers show only a single twitch at 48 hr APF, while both single and double twitches are observed in the control. Scale bar represents 20 µm.
Figure 7—video 2
Twitching in developing Act88F-GAL4 /+ and Act88F >> salmIR DLMs at 72 hr APF.
72 hr APF muscle ends are labelled with a Talin-YPet fusion in Act88F-GAL4 / + (control) and Act88F-GAL4/salmIR DLMs. Control fibers show no contractions at 72 hr APF, but salmIR fibers continue to contract. Scale bar represents 20 µm.
Figure 8
Strn-Mlck contributes to sarcomere length regulation during sarcomere maturation.
(A–J) Wild type (A–E) and StrnM102893 mutant (F–J) sarcomere development at 48, 72, 80, 90 hr APF and 1 day adult. (K–N) Tukey box and whisker plot of myofibril width (K) and sarcomere length (L) in wild type and StrnM102893 mutant (red). Tukey’s multiple comparison p-value<0.001***. N > 10 for each individual time point. Histogram of fiber length (M) and number of sarcomeres per myofibril (N). Error bars represent SEM. Tukey’s multiple comparison p-value<0.01**, N > 10 for each individual time point. Note that a normal number of sarcomeres are formed in StrnM102893 mutants, but they grow too long at 80 hr APF and hyper-contract in 1 day adult.
-
Figure 8—source data 1
This table includes the numbers of sarcomeres calculated per strn-mlck mutant myofibril and the respective control at the indicated stage.
- https://doi.org/10.7554/eLife.34058.035
-
Figure 8—source data 2
This table includes the length measurements of the indirect flight muscle fibers at the indicated stage for strn-mlck mutants and wild-type controls.
- https://doi.org/10.7554/eLife.34058.036
-
Figure 8—source data 3
This table includes the calculated numbers of sarcomeres per myofibril at the indicated stage for strn-mlck mutants.
Note that image analysis is limited by pixel size and thus can result in the same number.
- https://doi.org/10.7554/eLife.34058.037
-
Figure 8—source data 4
This table includes the sarcomere length measured in strn-mlck mutant flight muscles at the indicated stage.
Note that image analysis is limited by pixel size and thus can result in the same number.
- https://doi.org/10.7554/eLife.34058.038
Figure 9
Model for ordered sarcomere morphogenesis.
Phase 1 - Sarcomere assembly: Sarcomeric proteins are expressed at low levels, enabling the self-assembly of immature sarcomeres and myofibrils at around 30 hr APF. Each of the 2000 myofibrils contains about 80 immature sarcomeres. Phase 2 - Sarcomere addition: Sarcomere protein expression increases and each immature myofibril incorporates many new sarcomeres until about 48 hr APF. These immature sarcomeres contract spontaneously. Phase 3 - Sarcomere maturation: After 48 hr APF, sarcomere protein expression increases even more, including flight muscle specific proteins/isoforms, and all sarcomeres strongly grow in width and length. This enables the flight muscle to gain stretch-activation. Muscles are shown in red, tendons in blue. Structural proteins are illustrated as cartoons and are not drawn to scale. Induction of Actin, Mhc, Myofilin, Flightin, Strn-Mlck and Unc-89/Obscurin is shown here, expression of αActinin, Projectin and Kettin/Sallimus was shown by others (Bullard et al., 2006; Burkart et al., 2007; Katzemich et al., 2013; Orfanos et al., 2015; Weitkunat et al., 2014).
Tables
Key resources table
| Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
|---|---|---|---|---|
| Gene (Drosophila melanogaster) | spalt major; salm | NA | FLYB:FBgn0261648 | |
| Gene (D. melanogaster) | Stretchin-Mlck; Strn-Mlck | NA | FLYB:FBgn0265045 | |
| Genetic reagent (D. melanogaster) | w[1118] | Bloomington | BDSC:3605; FLYB:FBst0003605; RRID:BDSC_3605 | |
| Genetic reagent (D. melanogaster) | salmIR | PMID: 22094701 | VDRC:13302; FLYB:FBst0450930 | Flybase symbol: VDRC:v13302 |
| Genetic reagent (D. melanogaster) | KK101052 | PMID: 17625558 | VDRC:101052 | |
| Genetic reagent (D. melanogaster) | Act88F-GAL4 | PMID: 22008792 | Source: Richard Cripps | |
| Genetic reagent (D. melanogaster) | Strn-Mlck-MiMIC | Bloomington | FLYB:FBal0264439 | Flybase symbol:Strn-MlckMI02893 |
| Genetic reagent (D. melanogaster) | Strn-Mlck-IR | PMID: 21460824 | BDSC:31891; FLYB:FBti0130299; RRID:BDSC_31891 | Flybase symbol:P{TRiP.JF02170}attP2 |
| Genetic reagent (D. melanogaster) | Strn-Mlck-GFP, Isoform R | PMID: 25532219 | Symbol: Strn4; | |
| Genetic reagent (D. melanogaster) | Mhc-GFP | PMID: 26896675 | VDRC:318471 | Symbol: fTRG500; |
| Genetic reagent (D. melanogaster) | Mf-GFP | PMID: 26896675 | VDRC:318132 | Symbol: fTRG501; |
| Genetic reagent (D. melanogaster) | rhea-GFP | PMID: 26896675 | VDRC:318486 | Symbol: fTRG587; |
| Genetic reagent (D. melanogaster) | Fln-GFP | PMID: 26896675 | VDRC:318238 | Symbol: fTRG876; |
| Genetic reagent (D. melanogaster) | mys-GFP | PMID: 26896675 | VDRC:318285 | Symbol: fTRG932; |
| Genetic reagent (D. melanogaster) | βTub60D-GFP | PMID: 26896675 | VDRC:318309 | Symbol: fTRG958; |
| Genetic reagent (D. melanogaster) | unc-89-GFP | PMID: 26896675 | VDRC:318326 | Symbol: fTRG1046; |
| Genetic reagent (D. melanogaster) | Act88F-GFP | PMID: 26896675 | VDRC:318362 | Symbol:fTRG10028; |
| Genetic reagent (D. melanogaster) | Him-nuc-eGFP | PMID: 19324085 | Source: Michael V. Taylor | |
| Genetic reagent (D. melanogaster) | Him-GAL4 | this paper | ||
| Genetic reagent (D. melanogaster) | UAS-BBM | PMID: 22446736 | ||
| Genetic reagent (D. melanogaster) | Him-GFP-Gma | this paper | ||
| Genetic reagent (D. melanogaster) | rhea-YPet | this paper | ||
| Genetic reagent (D. melanogaster) | kon-GFP | this paper | ||
| Genetic reagent (D. melanogaster) | Mef2-GAL4 | Bloomington | BDSC:27390; RRID:BDSC_27390 | |
| Genetic reagent (D. melanogaster) | UAS-GFP-Gma | PMID: 12324971 | Source:Don Kiehart; Description | |
| Antibody | guinea pig anti-Shot | PMID: 9832554 | (1:500); Source: Talila Volk | |
| Antibody | rat anti-Kettin (MAC155/Klg16) | Babraham Bioscience Technologies | Babraham: MAC_155(P6689) | (1:50) |
| Antibody | rabbit anti-GFP (ab290) | Abcam | Abcam:ab290 | (1:1000) |
| Antibody | rat anti-Bruno | PMID: 12591598 | (1:500); Source: Anne Ephrussi | |
| Antibody | rabbit anti-Salm | PMID: 7905822 | (1:50); Source: Reinhard Schuh | |
| Antibody | mouse anti-βPS- integrin (CF.6G11) | Developmental Studies Hybridoma Bank | DSHB:CF.6G11 | (1:500) |
| Antibody | rabbit anti-Twi | PMID: 2688897 | (1:1000); Source: Siegfried Roth | |
| Antibody | rabbit anti-Fln | PMID: 11134077 | (1:50); Source: Jim Vigoreaux | |
| Commercial assay or kit | fluorescent beads | ThermoFisher (Molecular Probes) | OrderID: InSpeckTM Green Kit I-7219 | |
| Commercial assay or kit | Dynabeads | Invitrogen | OrderID: #610.06 | |
| Commercial assay or kit | Superscript III First-Strand Synthesis System | Invitrogen | OrderID: #18080–051 | |
| Chemical compound, drug | Fluoroshield with DAPI | Sigma | OrderID: #F6057 | |
| Chemical compound, drug | Vectashield with DAPI | Biozol | OrderID: VEC-H-1200 | |
| Chemical compound, drug | Tissue-Tek O.C.T. | Weckert Labotechnik | OrderID: 4583; Sakura Finetek | |
| Chemical compound, drug | TriPure reagent | Roche | OrderID: #11667157001 | |
| Software, algorithm | Fiji (Image J) | PMID: 22743772 | ||
| Software, algorithm | MyofibrilJ | this paper | 1dbb0d4 | Source: https://imagej.net/MyofibrilJ |
| Software, algorithm | STAR | PMID: 23104886 | ||
| Software, algorithm | SAMtools | PMID: 19505943 | ||
| Software, algorithm | featureCounts | PMID: 24227677 | ||
| Software, algorithm | DESeq2 | PMID: 25516281 | ||
| Software, algorithm | R | R Project for Statistical Computing | RRID:SCR_001905 | |
| Software, algorithm | ComplexHeatmap | PMID: 27207943 | ||
| Software, algorithm | corrplot | GitHub | Source: https://github.com/taiyun/corrplot | |
| Software, algorithm | VennDiagram | PMID: 21269502 | ||
| Software, algorithm | plyr | DOI: 10.18637/jss.v040.i01 | ||
| Software, algorithm | reshape2 | DOI: 10.18637/jss.v021.i12 | ||
| Software, algorithm | ggplot2 | ISBN:978-0-387-98140-6 | ||
| Software, algorithm | RColorBrewer | Author: Erich Neuwirth; | Source: https://cran.r-project.org/web/packages/RColorBrewer/index.html | |
| Software, algorithm | Mfuzz | PMID: 16078370 | ||
| Software, algorithm | GO-Elite | PMID: 22743224 |
Additional files
-
Supplementary file 1
mRNA-Seq raw data
The file includes multiple tabs containing the raw or input counts data from bioinformatics analysis, as well as a key to all original data provided in the supplementary tables. This table includes mRNA-Seq counts data, DESeq2 normalized counts data and standard normal counts data used for Mfuzz clustering for wild-type and salmIR IFM time points. The averaged core expression profiles for each Mfuzz cluster are also listed.
- https://doi.org/10.7554/eLife.34058.040
-
Supplementary file 2
GO-Elite analysis data.
This table includes multiple tabs containing the GO-Elite analysis of enrichments in Mfuzz clusters as well as genes up- or down-regulated from 30–72 hr APF and between wild-type and salmIR IFM. It also contains a complete list of all genes included in the ‘User Defined’ gene sets.
- https://doi.org/10.7554/eLife.34058.041
-
Supplementary file 3
Summary of sarcomere and myofibril quantifications
This table includes a numerical summary of quantification values reported graphically in Figures 5, 7 and 8. Quantifications of sarcomere length, myofibril width and myofibril diameter were performed with the MyofibrilJ script (see Materials and methods). Fiber length and cross-sectional area measurements were performed in Fiji/Image J.
- https://doi.org/10.7554/eLife.34058.042
-
Supplementary file 4
DESeq2 pairwise differential expression analysis
This table contains multiple tabs containing the output data from DESeq2 differential expression analysis between sequential IFM development time points, from 30 to 72 hr APF as well as between WT and salmIR IFM.
- https://doi.org/10.7554/eLife.34058.043
-
Supplementary file 5
Statistical Data
- https://doi.org/10.7554/eLife.34058.044
-
Transparent reporting form
- https://doi.org/10.7554/eLife.34058.045
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
A transcriptomics resource reveals a transcriptional transition during ordered sarcomere morphogenesis in flight muscle
eLife 7:e34058.
https://doi.org/10.7554/eLife.34058
