A transcriptomics resource reveals a transcriptional transition during ordered sarcomere morphogenesis in flight muscle
Abstract
Muscles organise pseudo-crystalline arrays of actin, myosin and titin filaments to build force-producing sarcomeres. To study sarcomerogenesis, we have generated a transcriptomics resource of developing Drosophila flight muscles and identified 40 distinct expression profile clusters. Strikingly, most sarcomeric components group in two clusters, which are strongly induced after all myofibrils have been assembled, indicating a transcriptional transition during myofibrillogenesis. Following myofibril assembly, many short sarcomeres are added to each myofibril. Subsequently, all sarcomeres mature, reaching 1.5 µm diameter and 3.2 µm length and acquiring stretch-sensitivity. The efficient induction of the transcriptional transition during myofibrillogenesis, including the transcriptional boost of sarcomeric components, requires in part the transcriptional regulator Spalt major. As a consequence of Spalt knock-down, sarcomere maturation is defective and fibers fail to gain stretch-sensitivity. Together, this defines an ordered sarcomere morphogenesis process under precise transcriptional control – a concept that may also apply to vertebrate muscle or heart development.
https://doi.org/10.7554/eLife.34058.001eLife digest
Animals may have different types of muscles but they all have one thing in common: molecular machines called sarcomeres that produce a pulling force. Conserved from fruit flies to humans, these structures line up end-to-end inside muscle cells, forming long cables called myofibrils. Some of the myofibrils in a human can reach several centimetres in length, which is much longer than those in a fruit fly. However, individual sarcomeres are the same length in both humans and flies.
To build the parts of the sarcomere, an animal cell first copies the relevant genes into intermediate molecules known as mRNAs, which are then translated to build new sarcomere proteins. Developing muscle cells later tune their sarcomeres to make them sensitive to stretching. This tweaks the power and force of the mature muscle, but the details of this developmental process are not fully understood.
Now, Spletter et al. have counted all the mRNAs in the developing flight muscles of fruit flies, with the aim of generating a resource that catalogues the changes in gene activity, or expression, that occur as muscles develop. This revealed that sarcomeres form in three phases. First, the cells assembled all their myofibrils. Then, they added short sarcomeres to the ends of their myofibrils. Finally, the sarcomeres matured to their full length and diameter, and became sensitive to stretching.
Fruit fly muscles had 40 patterns of gene expression, with most of the sarcomere components having one of two specific patterns. The expression of these genes dramatically rose after the young muscle cells had finished assembling all their myofibrils, suggesting muscles express different genes when their sarcomeres mature. A protein called spalt-major helped the cell to know when to make the transition, allowing the sarcomeres to grow in length and width.
Losing spalt-major late in muscle development stopped sarcomere growth and prevented the tuning process. The sarcomeres failed to become sensitive to stretching, a crucial feature of mature muscle. Muscles without spalt-major contracted too much and without coordination, like a muscle spasm.
The similarities between fruit fly and human sarcomeres suggest this developmental sequence may also occur in human muscles too. Understanding these steps may help to improve repair after injury or muscle growth during exercise. The next step is to test whether regenerating or growing muscles develop in the same way.
https://doi.org/10.7554/eLife.34058.002Introduction
Sarcomeres are the stereotyped force producing mini-machines present in all striated muscles of bilaterians. They are built of three filament types arrayed in a pseudo-crystalline order: actin filaments are cross-linked with their plus ends at the sarcomeric Z-disc and face with their minus ends towards the sarcomere center. In the center, symmetric bipolar muscle myosin filaments, anchored at the M-line, can interact with the actin filaments. Myosin movement towards the actin plus ends thus produces force during sarcomere shortening. Both filament types are permanently linked by a third filament type, the connecting filaments, formed of titin molecules (Gautel and Djinović-Carugo, 2016; Lange et al., 2006). A remarkable feature of sarcomeres is their stereotyped size, ranging from 3.0 to 3.4 µm in relaxed human skeletal muscle fibers (Ehler and Gautel, 2008; Llewellyn et al., 2008; Regev et al., 2011). Even more remarkable, the length of each bipolar myosin filament is 1.6 µm in all mature sarcomeres of vertebrate muscles, requiring about 300 myosin hexamers to assemble per filament (Gokhin and Fowler, 2013; Tskhovrebova and Trinick, 2003).
Human muscle fibers can be several centimetres in length and both ends of each fiber need to be stably connected to tendons to achieve body movements. As sarcomeres are only a few micrometres in length, many hundreds need to assemble into long linear myofibrils that span from one muscle end to the other and thus enable force transmission from the sarcomeric series to the skeleton (Lemke and Schnorrer, 2017). Thus far, we have a very limited understanding of how sarcomeres initially assemble into long immature myofibrils during muscle development to exactly match the length of the mature muscle fiber (Sparrow and Schöck, 2009). In particular, we would like to understand how such sarcomeres mature to the very precise stereotyped machines present in mature muscle fibers.
Across evolution, both the pseudo-crystalline regularity of sarcomeres as well as their molecular components are well conserved (Ehler and Gautel, 2008; Vigoreaux, 2006). Thus, Drosophila is a valid model to investigate the biogenesis of sarcomeres as well as their maturation. In particular, the large indirect flight muscles (IFMs) that span the entire fly thorax are an ideal model system to investigate mechanisms of myofibrillogenesis. They contain thousands of myofibrils consisting of 3.2 µm long sarcomeres (Schönbauer et al., 2011; Spletter et al., 2015).
Like all Drosophila adult muscles, IFMs are formed during pupal development from a pool of undifferentiated myoblasts called adult muscle precursors (AMPs) (Bate et al., 1991). From 8 hr after puparium formation (APF), these AMPs either fuse with themselves (for the dorso-ventral flight muscles, DVMs) or with remodelled larval template muscles (for the dorso-longitudinal flight muscles, DLMs) to form myotubes (Dutta et al., 2004; Fernandes et al., 1991). These myotubes develop dynamic leading edges at both ends and initiate attachment to their respective tendon cells at 12 to 16 hr APF (Weitkunat et al., 2014). These attachments mature and mechanical tension is built up in the myotubes, followed by the formation of the first immature periodic myofibrils at 30 hr APF when the muscle fibers are about 150 µm in length. These immature myofibrils contain the earliest sarcomeres, which are about 1.8 µm in length (Weitkunat et al., 2014). During the remaining 3 days of pupal development, the muscle fibers grow to about 1 mm to fill the entire thorax and sarcomere length increases to a final length of about 3.2 µm in adult flies (Orfanos et al., 2015; Reedy and Beall, 1993).
After myoblasts have fused to myotubes, the flight muscle specific selector gene spalt major (spalt, salm) is turned on in the developing flight muscle myotubes. Spalt major is responsible for the correct fate determination and development of the flight muscles, which includes the fibrillar flight muscle morphology and the stretch-activated muscle contraction mode (Schönbauer et al., 2011; Syme and Josephson, 2002). It does so by controlling the expression of more than 700 flight muscle specific genes or gene isoforms during development (Spletter and Schnorrer, 2014; Spletter et al., 2015). However, how the interplay between all these isoforms instructs the formation of highly regular, pseudo-crystalline sarcomeres in the flight muscle is not understood.
Here, we studied the transcriptional dynamics of flight muscle development in detail. We performed a systematic mRNA-Seq time-course of isolated muscle tissue at eight time points from the myoblast stage until the mature adult muscle stage, generating a transcriptomics resource of developing flight muscle. Bioinformatic analysis of expression dynamics identified two gene clusters that are strongly enriched for sarcomeric genes. The temporal dynamics of these clusters identified a transcriptional transition that is required for sarcomere morphogenesis. We define sarcomere morphogenesis in three sequential although overlapping phases. First, immature myofibrils assemble simultaneously; second, short sarcomeres are added to each myofibril; and third, all sarcomeres mature to their final length and diameter and acquire stretch-sensitivity. Interestingly, the number of myofibrils remains constant, suggesting that every sarcomere progresses through sarcomere maturation. We show that the flight muscle selector gene spalt major contributes to the observed transcriptional transition, suggesting that muscle fiber type-specific transcription is continuously required during sarcomere formation and maturation. Together, these findings indicate that precise transcriptional control of the sarcomeric components enables their ordered assembly into sarcomeres and their maturation to pseudo-crystalline regularity.
Results
A time-course of indirect flight muscle development
To better understand muscle morphogenesis in general and myofibrillogenesis in particular, we focused on the Drosophila indirect flight muscles (IFMs). We hypothesised that major morphological transitions during IFM development may be induced by transcriptional changes, thus we aimed to generate a detailed developmental mRNA-Seq dataset from IFMs. IFMs are built from AMPs that adhere to the hinge region of the wing disc epithelium and are labelled with GFP-Gma under Him control (Figure 1A) (Soler and Taylor, 2009). At 16 hr APF, many of these myoblasts have fused to larval template muscles to build the dorsal-longitudinal flight muscle (DLM) myotubes, which initiate attachment to their tendons. At this stage, the DLM myotubes of fibers 3 and 4 have a length of about 300 µm (Figure 1B). Fusion ceases at about 24 hr APF (Figure 1C) and attachment matures until 32 hr APF, coinciding with the strong recruitment of βPS-Integrin and the spectraplakin homolog Shortstop (Shot) to the attachment sites. At this stage the myofibers have built up mechanical tension and compacted to a length of about 150 µm, coinciding with the appearance of long Shot-positive tendon extensions that anchor the muscles within the thorax. This important developmental transition is highlighted by the assembly of immature myofibrils visualised by strong F-actin staining throughout the entire muscle fiber (Figure 1D) (Weitkunat et al., 2014).
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.34058.004
After 32 hr APF, the myofibers undergo another developmental transition and begin to grow dramatically. They elongate 3-fold to reach a length of about 480 µm by 48 hr APF (Figure 1E) and about 590 µm by 56 hr APF (Figure 1F). Concomitantly, the tendon extensions shrink with the myofibers being directly connected to the basal side of the tendon cell epithelium by 72 hr APF (Figure 1G). At the end of pupal development (90 hr APF), wavy muscle fibers with a length of about 780 µm containing mature myofibrils (Figure 1H) are present within the thorax.
A transcriptomics resource of indirect flight muscle development
To quantify transcriptional dynamics across the entire developmental time-course, we focused on the major developmental transitions and isolated mRNA from dissociated myoblasts of dissected or mass-isolated third instar wing discs and from hand-dissected IFMs at 16, 24, 30, 48, 72 and 90 hr APF pupae, and adult flies 1 day after eclosion (Figure 1I). We performed mRNA-Seq using at least two biological replicates for each time point (see Materials and methods). To identify genes with similar temporal expression profiles, we used Mfuzz (Kumar and Futschik, 2007) to cluster standard normalized read counts from all genes expressed above background (12,495 of 13,322 genes). This allowed us to confidently identify 40 distinct genome-wide clusters (Figure 2—figure supplement 1), each of which contains a unique gene set ranging from 155 to 703 members (Supplementary file 1). These clusters represent various temporal expression dynamics, with high expression at early (myoblast proliferation and fusion), mid (myotube attachment and myofibril assembly) or late (myofiber growth) myogenesis stages or a combination thereof (Figure 1I, Figure 2—figure supplement 1). These distinct patterns suggest a precise temporal transcriptional regulation corresponding to observed morphological transition points.
To verify the mRNA-Seq and cluster analysis, we selected a number of ‘indicator’ genes with available antibodies or GFP fusion proteins whose expression correlates with important developmental transitions. Twist (Twi) is a myoblast nuclear marker at larval stages and its expression needs to be down-regulated after myoblast fusion in pupae (Anant et al., 1998). We find twi mRNA in Mfuzz cluster 27, with high expression in myoblasts until 16 hr APF and a significant down-regulation from 24 hr APF, which we were able to verify with antibody stainings (Figure 2A–C). The flight muscle fate selector gene spalt major (salm) (Schönbauer et al., 2011) and its target, the IFM splicing regulator arrest (aret, bruno) (Spletter et al., 2015) are members of cluster 26 and 14, respectively. Expression of both clusters is up-regulated after myoblast fusion at 16 hr APF, which we were able to verify with antibody stainings (Figure 2D–F, Figure 2—figure supplement 2A–C). The apparent salm mRNA peak at 72 hr APF, which does not appear to cause a further protein increase, represents a mere 1.3 fold increase in expression and would need further confirmation as Salm, like many transcription factors, is expressed at low levels. For the initiation of muscle attachment, we selected Kon-tiki (Kon) (Schnorrer et al., 2007; Weitkunat et al., 2014), member of cluster 15, which is transiently up-regulated after myoblast fusion, before it is down-regulated again after 30 hr APF. Consistently, we found Kon-GFP present at muscle attachment sites at 30 hr APF but not at 72 hr APF (Figure 2G–I). A similar expression peak shifted to slightly later time points is found in cluster 34, which contains β-tubulin 60D (βTub60D) (Leiss et al., 1988). Consistently, we find β-Tub60D-GFP (Sarov et al., 2016) expression in IFMs at 30 and 48 hr but not 72 hr APF (Figure 2J–L). After attachment is initiated, the attachments need to mature and be maintained. As expected, we found the essential attachment components βPS-Integrin (mys) and Talin (rhea) in clusters that are up-regulated after myoblast fusion until adulthood (clusters 7 and 25, respectively). This is consistent with continuous high protein expression of βPS-Integrin-GFP and Talin-GFP at muscle attachment sites (Figure 2M–O, Figure 2—figure supplement 2D–F). Taken together, these semi-quantitative protein localisation data nicely validate the temporal mRNA dynamics found in the mRNA-Seq data, confirming our methodology.
A transcriptional transition after 30 hr APF
Hierarchical clustering of the core expression profiles from the 40 identified Mfuzz clusters defines eight temporally ordered groups (Figure 3) that show progressive expression dynamics as muscle development proceeds. A time-dependent shift in gene ontology (GO) term enrichments is apparent between the eight groups, reflecting the different stages of IFM development. GO-Elite analysis (Zambon et al., 2012) for gene set enrichment identified GO terms related to cell proliferation and development as enriched in the early clusters (such as the twi cluster 27 or the kon cluster 15), whereas terms related to actin filament dynamics are more enriched in the middle clusters (such as the βPS-Integrin cluster seven and the Talin cluster 25), reassuring that the clustering approach is valid (Figure 3, Supplementary file 2). Strikingly, the only two clusters that display a strong enrichment for genes important for sarcomere organisation are clusters 13 and 22, both of which are late up-regulated clusters (Figure 3). Members of both clusters just become detectable at 30 hr APF (Unc-89/Obscurin-GFP, Act88F-GFP, Mhc-GFP) or even later at 48 hr APF (Strn-Mlck-GFP, Mf-GFP). In all cases, we could confirm the strong up-regulation from 30 hr to 72 hr APF in the mRNA-Seq data at the protein level using GFP fusion proteins under endogenous control (Figure 2P–U, Figure 2—figure supplement 2J–P) (Sarov et al., 2016). We additionally verified the late up-regulation of Flightin (Fln), a member of cluster 3, which is detectable at 72 hr but not at 30 hr APF (Figure 2—figure supplement 2G–I).
At late stages of flight muscle development, mitochondrial density strongly increases (Clark et al., 2006). Using GO-Elite, we found a strong enrichment for mitochondrial related pathways in four late up-regulated clusters, namely 3, 28, 39 as well as the sarcomere cluster 22 (Figure 3). By comparing the clusters to systematic functional data acquired at all stages of Drosophila muscle development (Schnorrer et al., 2010), we find enrichments in clusters throughout the time-course. Interestingly, genes highly expressed in flight muscle compared to other muscle types, identified as ‘salm-core genes’ (Spletter et al., 2015), are also enriched in the late clusters, including the mitochondrial enriched clusters 3, 28, 39 and the sarcomere enriched cluster 22 (Figure 3). These data highlight the changes in biological process enrichments that parallel expression dynamics, with a particular transition happening during later stages of muscle development after 30 hr APF. This corresponds to a time period after immature myofibrils have been assembled, which thus far remained largely unexplored.
To examine the temporal expression dynamics in more detail, we performed a principle component analysis (PCA) of the mRNA-Seq time points and Mfuzz clusters and found that the major variance is developmental time, with a notable change after 30 hr APF (Figure 4A, Figure 4—figure supplement 1A). There are a large number of genes being up-regulated as well as down-regulated between 30 and 48 hr and between 48 and 72 hr APF, with major differences between the sets of genes expressed at early (16–30 hr APF) versus late (72–90 hr APF) stages of development (Figure 4B, Figure 4—figure supplement 1B). Thus, we focused our attention on the transcriptional transition between 30 and 72 hr APF, which correlates with major growth of the flight muscle fibers (Figure 1).
A large number of genes are significantly up- or down-regulated from 30 hr to 72 hr APF, as visualized on a volcano plot displaying log2fold changes (FC) (Figure 4C), suggesting a major change in gene expression. In particular, many sarcomeric genes are strongly up-regulated. To identify fine details in expression dynamics, we took all genes significantly up-regulated between 30 and 72 hr APF and performed hierarchical clustering of their DESeq2 normalized counts values (Figure 4D). We noted that many genes are turned on from 30 hr to 72 hr APF, whereas others are already expressed at 30 hr and strongly increase their expression until 72 hr (Figure 4D), suggesting a transcriptional transition after 30 hr APF. Consistently, we found GO-terms of cell proliferation, cell cycle and Notch signalling down-regulated, whereas actin cytoskeleton, sarcomere, muscle function and mitochondrial related gene sets are strongly up-regulated from 30 hr to 72 hr APF (Figure 4E). Finally, the genes up-regulated from 30 hr to 72 hr APF are enriched for sarcomeric proteins and the ‘salm core genes’ (Figure 4—figure supplement 1C,D), and the Mfuzz gene clusters with the most up-regulated members are the mitochondrial and both sarcomeric gene containing clusters 13 and 22 (Figure 4F). Together, these data suggest that in particular expression of the sarcomeric and mitochondrial genes is strongly induced after 30 hr APF.
Ordered phases during sarcomere morphogenesis
The strong up-regulation of sarcomeric gene expression after immature myofibrils have been assembled (Figure 1) (Weitkunat et al., 2014) caught our interest and prompted us to more closely investigate the later stages of myofibrillogenesis during which myofibers grow dramatically (Figure 1). We stained the myofibers with phalloidin to reveal myofibril morphology and with the titin isoform Kettin (an isoform of the sallimus gene) to label the developing Z-discs and systematically quantified sarcomere length and myofibril width (see Materials and methods) (Figure 5A, Supplementary file 3). By measuring the total muscle fiber length, we calculated the total number of sarcomeres per myofibril at a given stage. We found that the sarcomere length and width remain relatively constant at about 2.0 and 0.5 µm, respectively until 48 hr APF (Figure 5B,C, Supplementary file 3). However, the number of sarcomeres per myofibril dramatically increases from about 100 at 34 hr to about 230 at 48 hr APF. After 48 hr only a few more sarcomeres are added, resulting in about 270 sarcomeres per myofibril at 60 hr APF. This number remains constant until the fly ecloses (Figure 5D, Supplementary file 3). Moreover, by analysing fiber cross-sections we found that the growth of the individual myofibril diameter correlates with growth of the entire muscle fiber. Both fiber diameter and myofibril diameter remain constant from 30 hr to 48 hr APF. After 48 hr the myofibril diameter grows nearly 3-fold from 0.46 µm to 1.43 µm in adult flies (Figure 5E–G), while fiber cross-sectional area grows nearly 4-fold from 1,759 µm2 to 6,970 µm2. Strikingly, during the entire time period from 30 hr APF to adults, the total number of myofibrils per muscle fiber remains constant (about 2000 per muscle fiber, Figure 5H). Taken together, these quantitative data lead us to propose ordered but somewhat overlapping phases of sarcomere morphogenesis: (1) During the sarcomere assembly phase, about 100 immature sarcomeres self-assemble within each immature myofibril. (2) Many short sarcomeres are added to each myofibril increasing its length. (3) The short sarcomeres grow in length and thickness to reach the mature pseudo-crystalline pattern. No new myofibrils are built after the initial myofibril assembly phase.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.34058.013
-
Figure 5—source data 2
- https://doi.org/10.7554/eLife.34058.014
-
Figure 5—source data 3
- https://doi.org/10.7554/eLife.34058.015
-
Figure 5—source data 4
- https://doi.org/10.7554/eLife.34058.016
-
Figure 5—source data 5
- https://doi.org/10.7554/eLife.34058.017
-
Figure 5—source data 6
- https://doi.org/10.7554/eLife.34058.018
-
Figure 5—source data 7
- https://doi.org/10.7554/eLife.34058.019
We gained additional evidence to support this ordered myofibrillogenesis model on both the molecular and functional levels. First, the initial assembly versus the later sarcomere maturation complement the transition in gene expression we observe from 30 hr to 72 hr APF. Using members of the late induced Mfuzz clusters that contain sarcomeric components, we found that indeed a subset of sarcomeric proteins, such as Unc-89/Obscurin and Mhc, are already detectable in a periodic pattern on immature myofibrils at 30 hr. By contrast, other important components, including Fln, Myofilin (Mf) and Strn-Mlck, are only incorporated into myofibrils from 48 hr APF, showing high levels by 72 hr APF (Figure 5—figure supplement 1). Second, to investigate muscle function, we used Talin-YPet as a muscle attachment marker and quantified the number of spontaneous muscle contractions in intact pupae (see Materials and methods). Interestingly, we found that immature myofibrils already start to spontaneously contract at 30 hr APF. These spontaneous contractions increase in strength and frequency until 48 hr APF, but then cease, producing no detectable spontaneous contractions at 72 hr APF (Figure 5I,J, Figure 5—video 1). This demonstrates that during the sarcomere assembly phase, immature contractile myofibrils are generated, which then likely acquire stretch-sensitivity as the immature myofibrils mature and thus cease contracting.
Salm contributes to the transcriptional transition after 30 hr APF
How does the transcriptional transition of the various sarcomeric components instruct myofibrillogenesis? As the identified sarcomeric clusters 13 and 22 are enriched for ‘salm-core genes’ (Spletter et al., 2015) (Figure 3), we chose to investigate gene expression in developing spalt-major knock-down (salmIR) flight muscles compared to wild type (Supplementary file 4). salmIR IFM shows a strong down-regulation in gene expression, notably of mRNAs coding for sarcomeric and mitochondrial components at 24 and 30 hr APF, and in particular at 72 hr APF (Figure 6A–C, Figure 6—figure supplement 1A–D). The genes down-regulated in salmIR IFM are enriched for GO terms associated with sarcomere assembly, flight behaviour and mitochondrial genes, as well as for the mitochondrial Mfuzz clusters 3, 28, 39 and the sarcomeric Mfuzz clusters 13 and 22 (Figure 6—figure supplement 1E–G). Interestingly, members of cluster 22, which is strongly enriched for sarcomeric and mitochondrial genes, are not only down-regulated at 72 hr APF in salmIR (Figure 6B) but are also less strongly induced from 30 hr to 72 hr APF in salmIR muscle compared to wild type (Figure 6D), suggesting that salm in addition to other factors is indeed required for the strong induction of sarcomeric protein expression after 30 hr APF.
-
Figure 6—source data 1
- https://doi.org/10.7554/eLife.34058.025
-
Figure 6—source data 2
- https://doi.org/10.7554/eLife.34058.026
Salm is expressed shortly after myoblast fusion and constitutive knock-down of salm with Mef2-GAL4 results in a major shift of muscle fiber fate (Schönbauer et al., 2011), which may indirectly influence transcription after 30 hr APF. Hence, we aimed to reduce Salm levels only later in development, to directly address its role during sarcomere maturation. To this end, we knocked-down salm with the flight muscle specific driver Act88F-GAL4, which is expressed from about 18 hr APF and requires salm activity for its expression (Bryantsev et al., 2012; Spletter et al., 2015). This strategy enabled us to reduce Salm protein levels at 24 hr APF resulting in undetectable Salm levels at 72 hr APF (Figure 6—figure supplement 2). To test if Salm contributes to the transcriptional boost of sarcomeric components after 30 hr APF, we performed quantitative imaging using unfixed living flight muscles expressing GFP fusion proteins under endogenous control. We used green fluorescent beads to normalise the GFP intensity between different samples, and could verify the induction of Mhc, Unc-89, Fln and Strn-Mlck on the protein level from 30 hr to 72 hr APF (Figure 6—figure supplement 3). While overall sarcomere morphology is not strongly affected in Act88F >> salmIR muscles, we found that the levels of Strn-Mlck, Fln and Mhc proteins are strongly reduced at 90 hr as compared to wild-type controls (Figure 6D–H). This suggests that salm indeed contributes to the transcriptional transition that boosts the expression of a number of sarcomeric proteins after 30 hr APF.
To investigate the consequences of late salm knock-down, we quantified the myofibril and sarcomere morphology from 48 hr onwards. The myofibrils display a fibrillar morphology, confirming that the early function of Salm to determine IFM fate was unaffected by our late knock-down. At 72 hr APF and more prominently at 90 hr APF and in adults, Act88F >> salmIR myofibrils showed actin accumulations at broadened Z-discs (Figure 7A–H), which are often a landmark of nemaline myopathies (Sevdali et al., 2013; Wallgren-Pettersson et al., 2011). The myofibril width was not significantly different in these myofibrils (Figure 7I). However, the sarcomeres of Act88F >> salmIR muscles displayed a strong defect in sarcomere length growth after 48 hr APF (Figure 7J, Supplemental File 3), with sarcomeres only obtaining a length of 2.8 µm in adult flies, demonstrating that Salm in addition to other factors is required for normal sarcomere maturation.
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.34058.029
-
Figure 7—source data 2
- https://doi.org/10.7554/eLife.34058.030
-
Figure 7—source data 3
- https://doi.org/10.7554/eLife.34058.031
Salm function contributes to gain of stretch-activation during sarcomere maturation
Given the defects in sarcomere length and sarcomere gene expression in Act88F >> salmIR muscles, we explored the function of these abnormal muscle fibers. As expected, Act88F >> salmIR flies are flightless (Figure 7—figure supplement 1A) and we observed rupturing of the adult muscle fibers within 1d after eclosion (Figure 7—figure supplement 1B–G), demonstrating the importance of proper sarcomere maturation to prevent muscle atrophy. Based on our finding that spontaneous flight muscle contractions stop by 72 hr APF, we hypothesized that if Salm truly regulates sarcomere maturation, we may see spontaneous contraction defects during development. At 48 hr APF, Act88F >> salmIR fibers twitch, but less often than and without the double twitches observed in control fibers (Figure 7K,L,O; Figure 7—video 1). Strikingly, at 72 hr APF salmIR fibers fail to stop contracting and moreover show frequent and uncoordinated spontaneous contractions in which different myofibril bundles of the same fiber twitch at different times (Figure 7M–O, Figure 7—video 2), demonstrating that sarcomere maturation is indeed disrupted, with the likely consequence that myofibrils fail to acquire normal stretch-activation sensitivity.
To directly test the function of a sarcomeric component during the sarcomere maturation phase, we investigated the role of the prominently induced Salm target Strn-Mlck. Strn-Mlck is only expressed after 30 hr APF and is largely incorporated during sarcomere maturation (Figure 5—figure supplement 1E, Figure 6—figure supplement 3E), and thus is also a bone-fide example of a gene regulated during the transcriptional transition. In Strn-Mlck mutants, sarcomere and myofibril morphology, including myofibril width, is initially normal. However, at 80 hr APF the sarcomeres overgrow, consistently reaching lengths of more than 3.5 µm and resulting in slightly longer muscle fibers at 80 hr APF (Figure 8). After overgrowing, sarcomeres appear to hyper-contract resulting in short, thick sarcomeres in 1-day-old adults (Figure 8E,J,K,L). Like Act88F >> salmIR flies, Strn-Mlck mutant adults are flightless (Figure 7—figure supplement 1A) and display ruptured fibers during the first days of life (Figure 7—figure supplement 1K–M) (Spletter et al., 2015). Together, these data demonstrate that sarcomere maturation must be precisely controlled at the transcriptional level to enable the precise growth of sarcomeres to their final mature size. This ensures the lifelong function of the contractile apparatus of muscle fibers.
-
Figure 8—source data 1
- https://doi.org/10.7554/eLife.34058.035
-
Figure 8—source data 2
- https://doi.org/10.7554/eLife.34058.036
-
Figure 8—source data 3
- https://doi.org/10.7554/eLife.34058.037
-
Figure 8—source data 4
- https://doi.org/10.7554/eLife.34058.038
Discussion
A developmental muscle transcriptomics resource
In this study, we generated a systematic developmental transcriptomics resource from Drosophila flight muscle. The resource quantifies the transcriptional dynamics across all the major stages of muscle development over five days, starting with stem cell-like myoblasts and attaching myotubes to fully differentiated, stretch-activatable muscle fibers. We have specifically focused on the transcriptional regulation of sarcomere and myofibril morphogenesis; however, the data we provide cover all other expected dynamics, such as mitochondrial biogenesis, T-tubule morphogenesis, neuromuscular junction formation, tracheal invagination, etc. Thus, together with the available systematic functional data of Drosophila muscle development (Schnorrer et al., 2010), our data should be a versatile resource for the muscle community. It nicely complements existing systematic data from vertebrate muscle, which thus far are largely restricted to postnatal stages (Brinegar et al., 2017; Drexler et al., 2012; Lang et al., 2017; Zheng et al., 2009). Furthermore, Drosophila flight muscle contains a single muscle fiber type, in contrast to the mixed fiber types found in mammals (Schiaffino and Reggiani, 2011; Spletter and Schnorrer, 2014). Hence, in this model the transcriptional dynamics of a single fiber type muscle can be followed with unprecedented precision.
A transcriptional transition correlating with ordered sarcomere morphogenesis
Earlier work has shown that the flight muscle myotubes first attach to tendon cells and then build-up mechanical tension. This tension triggers the simultaneous assembly of immature myofibrils, converting the myotube to an early myofiber (Weitkunat et al., 2014). This suggested a tension-driven self-organisation mechanism of myofibrillogenesis (Lemke and Schnorrer, 2017). Here we discovered that myofibrillogenesis is not only regulated mechanically, but to a large extent also transcriptionally. This enabled us to extend our model for ordered myofibrillogenesis also to later developmental stages and to define three sequential although somewhat overlapping phases (Figure 9). During the sarcomere self-assembly phase at about 30 hr APF, a large number of genes coding for sarcomeric proteins, including Mhc, Act88F and Unc-89/Obscurin, become up-regulated to enable the self-organization of short, immature sarcomeres within thin, immature myofibrils. Strikingly, all of the about 2000 myofibrils assemble during this phase.
This is followed by a sarcomere addition phase during which a transcriptional transition is initiated and the expression of the sarcomeric proteins increases. Concomitantly, the muscle fibers grow in length by addition of new sarcomeres to all the immature myofibrils, increasing the sarcomere number from about 80 to 230 per fibril at 48 hr APF. These sarcomeres are contractile, but remain short and thin (Figure 9).
After the transcriptional transition, myofibrillogenesis enters the final sarcomere maturation phase. Proteins present in immature myofibrils like Mhc, Act88F and Unc-89/Obscurin are expressed to even higher levels, and additional, often flight-muscle specific proteins like Mf, Fln and the titin-related isoform Strn-Mlck, begin to be expressed at high levels and are incorporated into the maturing sarcomeres. This facilitates a dramatic growth of all immature sarcomeres in length and particularly in diameter with all 2000 myofibrils reaching a pseudo-crystalline regularity within about two days of development (Figure 9). Importantly, these matured sarcomeres no longer contract spontaneously, likely because they acquired the stretch-activated mechanism of contraction described for mature Drosophila flight muscles (Bullard and Pastore, 2011; Josephson, 2006).
Our ordered sarcomere morphogenesis model is strongly supported by the observation that the number of myofibrils remains largely constant during the entire sarcomere morphogenesis period, suggesting that in flight muscles no new myofibrils are added after the initial assembly of immature myofibrils at about 30 hr APF. As all sarcomeres are present shortly after 48 hr APF, they all need to mature simultaneously to achieve their final pseudo-crystalline regularity.
Our model is further supported by previous studies. We and others found that immature myofibrils have a width of about 0.5 µm (Weitkunat et al., 2014), which corresponds to about four thick filaments across each myofibril at the EM level at 42 hr APF (at 22°C) (Reedy and Beall, 1993). This ‘core’ myofibril structure built until 48 hr APF, when most sarcomeres have formed, is expanded dramatically after 48 hr APF, reaching a mature width of 1.5 µm, corresponding to about 35 thick filaments across each myofibril at the EM-level (Reedy and Beall, 1993). Recent data showed that the ‘core’ myofibril structure built until 48 hr APF contains already highly ordered actin filaments, which gain even higher order to reach their pseudo-crystalline regularity at 90 hr APF (Loison et al., 2018). In total, each adult myofibril contains around 800 thick filaments per sarcomere (Gajewski and Schulz, 2010). The ‘core’ myofibril structure was also revealed by the preferential recruitment of over-expressed actin isoforms (Röper et al., 2005) and more importantly, by selective incorporation of a particular Mhc isoform that is only expressed at mid-stages of flight muscle development (Orfanos and Sparrow, 2013). This Mhc isoform expression switch coincides with the global transition in sarcomeric gene expression between the sarcomere assembly and the sarcomere maturation phases that we defined here. It also fits with the recent discovery that the formin family member Fhos is selectively required for actin filament elongation and recruitment of new actin filaments and thus myofibril diameter growth after 48 hr APF (Shwartz et al., 2016). Expression of Fhos is also induced after 30 hr APF. Fhos is a member of Mfuzz Cluster 28, another strongly induced cluster, underscoring the general relevance of the transcriptional transition for ordered sarcomerogenesis.
Regulated active sarcomere contractions
Mature indirect flight muscles employ a stretch-activated mechanism of muscle contraction, thus Ca2+ is not sufficient to trigger muscle contractions without additional mechanical stretch (Bullard and Pastore, 2011; Josephson, 2006). This is different to cross-striated body muscles of flies or mammals that contract synchronously with Ca2+ influx. Hence, it is intriguing that immature flight muscle myofibrils do in fact contract spontaneously, with the contraction frequencies and intensities increasing until 48 hr APF. It was recently proposed in Drosophila cross-striated abdominal muscles and in developing cross-striated zebrafish muscles that spontaneous contractions are important for the proper formation of the cross-striated pattern (Mazelet et al., 2016; Weitkunat et al., 2017). A similar role for contractions was found in C2C12 cells by stimulating the contractions optogenetically (Asano et al., 2015). This shows that spontaneous contractions are a necessary general feature for the assembly of cross-striated muscle fibers across species.
However, flight muscles are not cross-striated in the classical sense, but have a fibrillar organisation in which each myofibril remains isolated and is not aligned with its neighbouring myofibrils (Figure 5) (Josephson, 2006; Schönbauer et al., 2011). We can only speculate about the mechanism that prevents alignment of the myofibrils in the flight muscles, but it is likely related to their stretch-activated contraction mechanism. This mechanism prevents spontaneous twitching due to increased Ca2+ levels, because it additionally requires mechanical activation that can only occur during flight in the adult. Thus, flight muscle sarcomeres not only grow and mature their sarcomere structure, they also gain their stretch-activatability.
Continuous maintenance of muscle type-specific fate
We identified an important transition in gene expression between the early sarcomere assembly and the late sarcomere maturation phases. Similar large scale transcriptome changes have also been observed during postnatal stages of mouse (Brinegar et al., 2017), chicken (Zheng et al., 2009) and pig (Zhao et al., 2015) skeletal muscle development or during regeneration after injury in fish (Montfort et al., 2016) and mouse muscles (Warren et al., 2007), indicating that muscle maturation generally correlates with large scale transcriptional changes.
It is well established that general myogenic transcription factors, in particular Mef2, are continuously required in muscles for their normal differentiation (Sandmann et al., 2006; Soler et al., 2012). Mef2 regulates a suit of sarcomeric proteins in fly, fish and mouse muscle important for correct sarcomere assembly and maturation (Hinits and Hughes, 2007; Kelly et al., 2002; Potthoff et al., 2007; Stronach et al., 1999). In Drosophila, Mef2 cooperates with tissue-specific factors, such as CF2, to induce and fine-tune expression of structural genes (Gajewski and Schulz, 2010; García-Zaragoza et al., 2008; Tanaka et al., 2008). General transcriptional regulators, such as E2F, further contribute to high levels of muscle gene expression during myofibrillogenesis, in part through regulation of Mef2 itself (Zappia and Frolov, 2016). However, it is less clear if muscle type-specific identity genes are continuously required to execute muscle type-specific fate. Spalt major (Salm) is expressed after myoblast fusion in flight muscle myotubes and is required for all flight muscle type-specific gene expression: in its absence the fibrillar flight muscle is converted to tubular cross-striated muscle (Schönbauer et al., 2011; Spletter et al., 2015). Here we demonstrated that Salm is continuously required for correct sarcomere morphogenesis, as late salm knock-down leads to defects in sarcomere growth during the late sarcomere maturation phase, causing severe muscle atrophy in adults. It might do so by modifying the cooperation between Mef2 and E2F or by changing chromatin states, as vertebrate spalt homologs are recently discussed as epigenetic regulators (Yang, 2018; Zhang et al., 2015). However, Salm cannot be solely responsible for the transcriptional transition after 30 hr APF, as the transition still partially occurs in its absence.
Ordered sarcomere morphogenesis – a general mechanism?
Here we defined ordered phases of sarcomere morphogenesis in Drosophila flight muscles. Is this a general concept for sarcomere morphogenesis? Reviewing the literature, one finds that in other Drosophila muscle types which display a tubular cross-striated myofibril organisation, such as the fly abdominal muscles, the striated sarcomeres also first assemble and then grow in length (PerezPérez-Moreno et al., 2014; Weitkunat et al., 2017), suggesting a conserved mechanism. In developing zebrafish skeletal muscles, young myofibers present in younger somites show a short sarcomere length of about 1.2 µm, which increases to about 2.3 µm when somites and muscle fibers mature (Sanger et al., 2017; Sanger et al., 2009). Interestingly, sarcomere length as well as thick filament length increase simultaneously during fish muscle maturation, indicating that as in flight muscles, the length of all sarcomeres in one large muscle fiber is homogenous at a given time (Sanger et al., 2009). Similar results were obtained in mouse cardiomyocytes measuring myosin filament length at young (two somite) and older (13 somite) stages (Du et al., 2008) and even in human cardiomyocytes, in which myofibrils increase nearly threefold in width and become notably more organized and contractile from 52 to 127 days of gestation (Racca et al., 2016). These data suggest that sarcomeres generally may go through a series of ordered but overlapping developmental phases.
Interestingly, these changes in sarcomere morphology correlate with a switch in myosin heavy chain isoform expression changing from embryonic to neonatal to adult during skeletal muscle development (Schiaffino et al., 2015). Importantly, mutations in embryonic myosin (MYH3 in humans) result in severe congenital disorders characterised by multiple facial and limb contractures (Toydemir et al., 2006). As a similar ordered expression of myosin isoforms is also found during muscle regeneration after injury in adults (Ciciliot and Schiaffino, 2010; Schiaffino et al., 2015), we hypothesize that our sequential sarcomere morphogenesis model may also be applicable to vertebrate skeletal and possibly heart muscles. It will be a future challenge to identify possible feedback mechanisms that indicate the successful end of the sarcomere assembly phase or a possible re-entry into the sarcomere assembly phase during muscle regeneration or exercise induced muscle fiber growth. It is enticing to speculate that the assembling cytoskeleton itself would measure its assembly status and provide a mechanical feedback signal to modify the activity of muscle-specific transcription factors.
Materials and methods
Fly strains
Request a detailed protocolFly stocks were maintained using standard culture conditions. Characterization of normal IFM sarcomere and fiber growth was performed in w1118 grown at 27°C. salm RNAi was performed with previously characterized GD3029 (referred to as salmIR) and KK181052 (Schönbauer et al., 2011) from VDRC (http://stockcenter.vdrc.at) (Dietzl et al., 2007) at 25°C using Act88F-GAL4 to induce knock-down after 24 hr APF. Act88F-GAL4 x w1118 served as control. The Strn-Mlck-MiMIC insertion MI02893 into IFM-specific IsoR (Bloomington stock 37038) and TRiP hairpin JF02170 were obtained from Bloomington (Ni et al., 2011). The salm-EGFP line was used to sort wing discs (Marty et al., 2014). Tagged genomic fosmid reporter fly lines include strn4 (Strn-Mlck-GFP, Isoform R) (Spletter et al., 2015), fTRG500 (Mhc-GFP, Isoforms K, L, M), fTRG501 (Mf-GFP, Isoforms A, G, N), fTRG587 (Rhea-GFP, Isoforms B, E, F, G), fTRG876 (Fln-GFP), fTRG932 (mys-GFP), fTRG958 (βTub60D-GFP), fTRG1046 (unc-89-GFP), and fTRG10028 (Act88F-GFP) (Sarov et al., 2016).
To label myoblasts, we utilized the enhancer for Holes-in-muscle (Him), which is expressed in dividing myoblasts and promotes the progenitor fate. Him-nuc-eGFP flies were a gift of M. Taylor (Soler and Taylor, 2009). Him-GAL4 flies were created by cloning an EcoRI to SacII fragment of the Him enhancer (Liotta et al., 2007) upstream of GAL4 into pStinger. UAS-BBM (UAS-palmCherry) (Förster and Luschnig, 2012) was driven with Him-GAL4 to label myoblasts. Him-GFP-Gma flies were created by PCR amplifying GFP-Gma with AscI and PacI overhangs and then cloning downstream of the Him enhancer in pStinger to generate a gypsy insulator-Himenh-Gma-GFP-SV40-gypsy insulator cassette.
The rhea-YPet line used to label muscle ends for live imaging of twitch events was generated by CRISPR-mediated gene editing at the endogenous locus (S.B.L and F.S., details will be published elsewhere). The kon-GFP line was generated by inserting GFP into the kon locus after its transmembrane domain using the genomic fosmid FlyFos021621, which was integrated using Φ-C31 into VK00033 (I. Ferreira and F.S., details will be published elsewhere).
Flight tests
Request a detailed protocolFlight tests were performed as previously described (Schnorrer et al., 2010). Act88F-GAL4 crosses were kept at 25°C, as higher temperatures negatively impacted flight ability, because of the very high GAL4 expression levels in this strain. Adult males were collected on CO2 and recovered at least 24 hr at 25°C before testing. Flies were introduced into the top of a 1 m long cylinder divided into five zones. Those that landed in the top two zones were considered ‘normal fliers’, those in the next two zones ‘weak fliers’ and those that fell to the bottom of the cylinder ‘flightless’.
Immuno-staining
Request a detailed protocolWing-discs were dissected from 3rd instar wandering larvae in 1x PBS and fixed in 4% PFA in PBS-T. Discs were stained as described below for anti-GFP. Adult and pupal flight muscles were dissected and stained as previously described (Weitkunat and Schnorrer, 2014). Briefly, early pupae (16–60 hr APF) were freed from the pupal case, fixed for 20 min. in 4% PFA in relaxing solution and washed in 0.5% PBS-Triton-X100 (PBS-T). 72 hr APF and older samples were cut sagittally with a microtome blade. All samples were blocked for at least 1 hr at RT in 5% normal goat serum in PBS-T and stained with primary antibodies overnight at 4°C. Primary antibodies include: guinea pig anti-Shot 1:500 (gift of T. Volk), rat anti-Kettin 1:50 (MAC155/Klg16, Babraham Institute), rabbit anti-GFP 1:1000 (ab290, Abcam), rat anti-Bruno 1:500 (Filardo and Ephrussi, 2003), rabbit anti-Salm 1:50 (Kühnlein et al., 1994), mouse anti-βPS-integrin 1:500 (CF.6G11, DSHB), rabbit anti-Twi 1:1000 (gift of Siegfried Roth) and rabbit anti-Fln 1:50 (Reedy et al., 2000) (gift of Jim Vigoreaux). Samples were washed three times in 0.5% PBS-T and incubated overnight at 4°C with secondary conjugated antibodies (1:500) from Invitrogen (Molecular Probes) including: Alexa488 goat anti-guinea pig IgG, Alexa488 donkey anti-rat IgG, Alexa488 goat anti-mouse IgG, Alexa488 goat anti-rabbit IgG, rhodamine-phalloidin, Alexa568 goat anti-rabbit IgG and Alexa633 goat anti-mouse IgG. Samples were washed three times in 0.5% PBS-T and mounted in Vectashield containing DAPI.
Cryosections
Request a detailed protocolHead, wings and abdomen were removed from one day old w1118 flies and thoraxes were fixed overnight at 4°C in 4% PFA. For 30–90 hr APF samples, pupae were freed from the pupal case, poked 3–5 times with an insect pin in the abdomen and fixed overnight at 4°C in 4% PFA. Thoraxes or pupae were then sunk in 30% sucrose in 0.5% PBS-T overnight at 4°C on a nutator. Thoraxes or pupae were embedded in Tissue-Tek O.C.T. (Sakura Finetek) in plastic moulds (#4566, Sakura Finetek) and frozen on dry ice. Blocks were sectioned at 30 µm on a cryostat (Microm vacutome). Sections were collected on glass slides coated with 1% gelatin +0.44 μM chromium potassium sulfate dodecahydrate to facilitate tissue adherence. Slides were post-fixed for 1 min. in 4% PFA in 0.5% PBS-T at RT, washed in 0.5% PBS-T, incubated with rhodamine-phalloidin for 2 hr at RT, washed three times in 0.5% PBS-T and mounted in Fluoroshield with DAPI (#F6057, Sigma).
Microscopy and image analysis
Request a detailed protocolImages were acquired with a Zeiss LSM 780 confocal microscope equipped with an α Plan-APOCHROMAT 100x oil immersion objective lens (NA 1.46). To compare if indicator protein expression replicates the mRNA-Seq expression dynamics, we imaged three time points from each expression profile with the same confocal settings. Laser gain and pinhole settings were set on the brightest sample and reused on remaining time points in the same imaging session. All samples were additionally stained with the same antibody mix on the same day and if possible in the same tube. Images were processed with Fiji (Schindelin et al., 2012) and Photoshop, and displayed using the ‘Fire’ look-up table.
Fiber length and fiber cross-sectional area were measured with freehand drawing tools in Fiji based on rhodamine-phalloidin staining. Sarcomere length, myofibril width, and myofibril diameter were measured automatically using a custom Fiji plug-in, MyofibrilJ, available from https://imagej.net/MyofibrilJ and code from https://github.com/giocard/MyofibrilJ (Cardone, 2018; copy archived at https://github.com/elifesciences-publications/MyofibrilJ). All measurements are based on rhodamine-phalloidin staining, except 34 hr APF sarcomere lengths, which are based on both rhodamine-phalloidin and Unc-89-GFP staining. ‘Sarcomeres per fibril’ was calculated as average individual fiber length divided by sarcomere length for fiber 3 or 4. ‘Fibrils per fiber’ was calculated as average number of fibrils per unit area multiplied by individual fiber cross sectional area.
For determining myofibril diameter, samples were imaged using a 3x optical zoom (50 nm pixel size). At least 20 cross-section images from different fibers for >10 flies were acquired for each time point. The number of fibrils per section and fibril diameter were determined with the tool ‘analyze myofibrils crosswise’ from MyofibrilJ (https://imagej.net/MyofibrilJ). In this tool, an initial estimate of the diameter is obtained by finding the first minimum in the radial average profile of the autocorrelation (Goodman, 1996) of the image. This estimate is used to calibrate the optimal crop area around all the cross-sections in the image, their position previously detected by finding the local intensity peaks. All of the detected cross-sections are then combined to obtain a noise-free average representation of the fibril section. Finally, the diameter is calculated by examining the radial profile of the average and measuring the full width where the intensity is 26% of the maximum range.
For determining sarcomere length and myofibril width, for each experiment between 10 and 25 images were acquired from more than 10 individual flies. From each image, nine non-overlapping regions of interest were selected, which were rotated to orient fibrils horizontally, when necessary. The tool ‘analyze myofibrils lengthwise’ from MyofibrilJ reports the sarcomere length (indicated as repeat) and myofibril width (indicated as thickness). Because of the periodic nature of sarcomere organization, their length is estimated by means of Fourier analysis, identifying the position of the peaks on the horizontal axis of the Fourier transformed imaged. Quality of the estimate was evaluated by visual inspection of the Fourier transformed image, overlaid with the peaks detected, as generated by the plug-in. Myofibrils width is estimated from the position of the first minimum in the vertical intensity profile of the autocorrelation of the image.
Live imaging of developmental spontaneous contractions was performed on a Leica SP5 confocal microscope. Prior to imaging, a window was cut in the pupal case, and pupae were mounted in slotted slides as previously described (Lemke and Schnorrer, 2018; Weitkunat and Schnorrer, 2014). At the specified developmental time point, IFMs were recorded every 0.65 s for 5 min. General movement within the thorax was distinguished from IFM-specific contraction, and each sample was scored for the number of single or double contractions observed per 5 min time window. Data were recorded in Excel and ANOVA was performed in GraphPad Prism to determine significant differences. Movies were assembled in Fiji (Image J), cropped and edited for length to highlight a selected twitch event.
Quantitative imaging of fosmid reporter intensity was performed at 30 hr APF, 48 hr APF, 72 hr APF, 90 hr APF and in 1 day adult in live IFM by normalizing to fluorescent beads (ThermoFisher (Molecular Probes), InSpeck Green Kit I-7219). IFMs were dissected from five flies, mounted with fluorescent microspheres (0.3% or 1% relative intensity, depending on the reporter intensity) in the supplied mounting medium and immediately imaged (within 20 min). Intensity measurements were obtained at 40x for at least 10 flies in regions where both IFM and at least three beads were visible. Control Act88F-GAL4;; fosmid-GFP x w1118 and Act88F-GAL4;; fosmid-GFP x salmIR (fosmids used included Strn-Mlck-GFP, Mhc-GFP, Fln-GFP, and Unc-89-GFP) were imaged in the same imaging session. Relative fluorescence fiber to beads was calculated for each image in Fiji by averaging intensity for three fiber ROIs and three bead ROIs. Data were recorded in Excel and Student’s t-test for significance and plotting were performed in GraphPad Prism.
mRNA-Seq
Request a detailed protocolWe previously published mRNA-Seq analysis of dissected IFMs from Mef2-GAL4, UAS-GFP-Gma x w1118 at 30 hr APF, 72 hr APF and 1d adult, and Mef2-GAL4, UAS-GFP-Gma x salmIR in 1d adult (Spletter et al., 2015). We expanded this analysis in the present study to include myoblasts from third instar larval wing discs (see below) and dissected IFMs from Mef2-GAL4, UAS-GFP-Gma x w1118 at 16, 24, 30, 48, 72, 90 hr APF and from 1 day adults as well as IFMs from Mef2-GAL4, UAS-GFP-Gma x salmIR flies at 24, 30, 72 hr APF and from 1 day adults. IFMs were dissected from groups of 15 flies in 30 min to minimize changes to the transcriptome, spun down in PBS for 5 min at 7500 rpm and immediately frozen in 100 µl TriPure reagent (#11667157001, Roche) on dry ice. RNA was isolated after combining IFMs from 150 to 200 flies, with biological duplicates or triplicates for each time point.
Poly(A)+mRNA was purified using Dynabeads (#610.06, Invitrogen) and integrity was verified on a Bioanalyzer. mRNA was then fragmented by heating to 94°C for 210 s in fragmentation buffer (40 mM TrisOAc, 100 mM KOAc, 30 mM MgOAc2). First-strand cDNA synthesis was performed with the Superscript III First-Strand Synthesis System (#18080–051, Invitrogen) using random hexamers. The second strand was synthesized with dUTP and submitted to the Vienna Biocenter Core Facilities (VBCF, http://www.vbcf.ac.at) for stranded library preparation according to standard Illumina protocols and sequenced as SR100 on an Illumina HiSeq2500. Libraries were multiplexed two to four per lane using TrueSeq adaptors.
Wing disc sorting and myoblast isolation
Request a detailed protocolTo perform mRNA-Seq on fusion competent myoblasts that will form the IFMs, we first dissected wing discs from wandering third instar larvae and manually cut the hinge away from the wing pouch. mRNA was isolated in TriPure reagent and sequenced as described above. We estimate this sample (Myo1) is ~50% myoblast, as the myoblasts form a nearly uniform layer over the underlying epithelial monolayer. To obtain a purer myoblast sample, we performed large-scale imaginal disc sorting followed by dissociation. We used particle sorting to isolate imaginal discs from Him-GAL4, UAS-BBM (UAS-palmCherry); salm-EGFP flies based on the green fluorescent signal. 10–12 ml of larvae in PBS were disrupted using a GentleMACS mixer (Miltenyl Biotec) and discs were collected through a mesh sieve (#0278 in, 25 opening, 710 µm). Fat was removed by centrifugation for 10 min. at 1000 rpm at 4°C, discs were rinsed in PBS and then re-suspended in HBSS. Discs were further purified on a Ficoll gradient (25%:16%). Discs were then sorted on a Large Particle Flow Cytometer (BioSorter (FOCA1000), Union Biometrica, Inc.), obtaining 600–1000 discs per sample. Discs were spun for 5 min at 600 rcf in a Teflon Eppendorf tube and then re-suspended in the dissociation mixture (200 µl of 10x Trypsin, 200 µl HBSS, 50 µl collagenase (10 mg/mL), 50 µl dispase (10 mg/ml)). The tube was incubated for 10 min. at RT and then transferred to a thermal shaker for 30 min. at 25°C at 650 rpm. Myoblasts were filtered through a 35 µm tube-cap filter and spun at 600 rcf for 5 min. to pellet the cells. Cells were re-suspended in HBSS for evaluation or frozen in TriPure reagent for RNA extraction. We obtained samples with ~90% purity based on counting the number of red fluorescent cells/non fluorescent+green fluorescent cells in three slide regions. mRNA was isolated in TriPure reagent and sequenced as described above, generating the Myo2 and Myo3 samples.
Analysis of RNA-Seq data
Request a detailed protocolFASTA files were de-multiplexed and base called using Illumina software. Reads were trimmed using the FASTX-toolkit. Sequences were mapped using STAR (Dobin et al., 2013) to the Drosophila genome (BDGP6.80 from ENSEMBL). Mapped reads were sorted and indexed using SAMtools (Li et al., 2009), and then bam files were converted to bigwig files. Libraries were normalized based on library size and read-counts uploaded to the UCSC Browser for visualization (http://genome.ucsc.edu/cgi-bin/hgTracks?hgS_doOtherUser=submit&hgS_otherUserName=Ayeroslaviz&hgS_otherUserSessionName=IFMTP.leg.TCpaperHub1; http://genome.ucsc.edu/cgi-bin/hgTracks?hgS_doOtherUser=submit&hgS_otherUserName=Ayeroslaviz&hgS_otherUserSessionName=AretSalmIFMTP.TCpaperHub2).
Mapped sequences were run through featureCounts (Liao et al., 2014) and differential expression analysis was performed on the raw counts using DESeq2 (Love et al., 2014). We performed all pairwise comparisons across the time-course as well as between wild-type and salmIR samples (Supplementary file 1 and 4). All original and processed data can be found as supplemental data or in the Gene Expression Omnibus submission (accession number GSE107247). R packages employed in the analysis include ComplexHeatmap (Gu et al., 2016), CorrPlot (https://github.com/taiyun/corrplot), VennDiagram (Chen and Boutros, 2011), plyr (Wickham, 2011), reshape2 (Wickham, 2007), ggplot2 (Wickham, 2009) and RColorBrewer (Neuwirth, 2015).
Genome-wide soft clustering was performed in R with Mfuzz (Futschik and Carlisle, 2005), using the DESeq2 normalized count values. We filtered the dataset to include all genes expressed at one time point or more, defining expression as >100 counts after normalization. We then set all count values < 100 to 0, to remove noise below the expression threshold. DESeq2 normalized data was standardized in Mfuzz to have a mean value of zero and a standard deviation of one, to remove the influence of expression magnitude and focus on the expression dynamics. We tested ‘k’ ranging from 10 to 256. We then performed consecutive rounds of clustering to obtain three independent replicates with similar numbers of iterations, ultimately selecting a final k = 40 clusters with iterations equal to 975, 1064 and 1118. We calculated a ‘stability score’ for each cluster by calculating how many genes are found in the same cluster in each run (Supplementary file 1). Figures are from the 1064 iterations dataset. Mfuzz cluster core expression profiles were calculated as the average standard-normal expression of all genes with a membership value greater than or equal to 0.8, and then core profiles were clustered in R using Euclidean distance and complete linkage.
Enrichment analysis was performed with GO-Elite (Zambon et al., 2012) using available Gene Ontology terms for Drosophila. We additionally defined user provided gene lists for transcription factors, RNA binding proteins, microtubule associated proteins, sarcomeric proteins, genes with an RNAi phenotype in muscle (Schnorrer et al., 2010), mitochondrial genes (http://mitoXplorer.biochem.mpg.de) and salm core fibrillar genes (Spletter et al., 2015). Full results and gene lists are available in Supplementary file 2. These user-supplied lists allowed us to define more complete gene sets relevant to a particular process or with a specific localization than available in existing GO terms. Analysis was performed with 5000 iterations to generate reliable significance values.
Data availability
Request a detailed protocolProcessed data from DESeq2, Mfuzz and GO-Elite are available in Supplementary file 1, 2, 4. mRNA-Seq data are publicly available from NCBI’s Gene Expression Omnibus (GEO) under accession number GSE107247. mRNA-Seq read counts are further publicly accessible as track hubs in the UCSC Genome Browser at the following links: [1] (wild-type IFM time course) http://genome.ucsc.edu/cgi-bin/hgTracks?hgS_doOtherUser=submit&hgS_otherUserName=Ayeroslaviz&hgS_otherUserSessionName=IFMTP.leg.TCpaperHub1 and [2] (salm timecourse) http://genome.ucsc.edu/cgi-bin/hgTracks?hgS_doOtherUser=submit&hgS_otherUserName=Ayeroslaviz&hgS_otherUserSessionName=AretSalmIFMTP.TCpaperHub2. Fiji scripts for analysis of sarcomere length, myofibril width and myofibril diameter are available from https://imagej.net/MyofibrilJ. Raw data used to generate all plots presented in figure panels are available in the source data files for Figures 1, 5, 6, 7 and 8. Data on statistical test results are presented in Supplementary file 5.
Data availability
Processed data from DESeq2, Mfuzz and GO-Elite are available in Supplementary Files 1, 2, 4. mRNA-Seq data are publicly available from NCBI's Gene Expression Omnibus (GEO) under accession number GSE107247. Fiji scripts for analysis of sarcomere length, myofibril width and myofibril diameter are available from https://imagej.net/MyofibrilJ. Raw data used to generate all plots presented in figure panels are available in the source data files for Figures 1, 5, 6, 7 and 8. Data on statistical test results are presented in Supplementary File 5.
-
Systematic transcriptomics reveals a biphasic mode of sarcomere morphogenesis in flight muscles regulated by SpaltPublicly available at the NCBI Gene Expression Omnibus (accession no: GSE107247).
-
The RNA binding protein Arrest (Aret) regulates myofibril maturation in Drosophila flight musclePublicly available at the NCBI Gene Expression Omnibus (accession no: GSE63707).
References
-
Twist and notch negatively regulate adult muscle differentiation in DrosophilaDevelopment 125:1361–1369.
-
Cells with persistent twist expression are the embryonic precursors of adult muscles in DrosophilaDevelopment 113:79–89.
-
The function of elastic proteins in the oscillatory contraction of insect flight muscleJournal of Muscle Research and Cell Motility 26:479–485.https://doi.org/10.1007/s10974-005-9032-7
-
Regulating the contraction of insect flight muscleJournal of Muscle Research and Cell Motility 32:303–313.https://doi.org/10.1007/s10974-011-9278-1
-
Modular proteins from the Drosophila sallimus (sls) gene and their expression in muscles with different extensibilityJournal of Molecular Biology 367:953–969.https://doi.org/10.1016/j.jmb.2007.01.059
-
Regeneration of mammalian skeletal muscle. basic mechanisms and clinical implicationsCurrent Pharmaceutical Design 16:906–914.
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
On marathons and sprints: an integrated quantitative proteomics and transcriptomics analysis of differences between slow and fast muscle fibersMolecular & Cellular Proteomics 11:M111.010801.https://doi.org/10.1074/mcp.M111.010801
-
Cardiac myofibrillogenesis inside intact embryonic heartsDevelopmental Biology 318:236–246.https://doi.org/10.1016/j.ydbio.2008.03.011
-
The sarcomere and sarcomerogenesisAdvances in Experimental Medicine and Biology 642:1–14.
-
Development of the indirect flight muscles of DrosophilaDevelopment 113:67–77.
-
Bruno regulates gurken during Drosophila oogenesisMechanisms of Development 120:289–297.https://doi.org/10.1016/S0925-4773(02)00454-9
-
Noise-robust soft clustering of gene expression time-course dataJournal of Bioinformatics and Computational Biology 03:965–988.https://doi.org/10.1142/S0219720005001375
-
CF2 activity and enhancer integration are required for proper muscle gene expression in DrosophilaMechanisms of Development 125:617–630.https://doi.org/10.1016/j.mod.2008.03.003
-
The sarcomeric cytoskeleton: from molecules to motionJournal of Experimental Biology 219:135–145.https://doi.org/10.1242/jeb.124941
-
A two-segment model for thin filament architecture in skeletal muscleNature Reviews Molecular Cell Biology 14:113–119.https://doi.org/10.1038/nrm3510
-
BookComparative Physiology of Insect Flight MuscleIn: Josephson R, editors. Nature's Versatile Engine: Insect Flight Muscle Inside and Out. Georgetown, TX: Landes Bioscience. pp. 35–43.
-
Spalt encodes an evolutionarily conserved zinc finger protein of novel structure which provides homeotic gene function in the head and tail region of the Drosophila embryoThe EMBO Journal 13:168–179.
-
Dynamic changes in the mouse skeletal muscle proteome during denervation-induced atrophyDisease Models & Mechanisms 10:881–896.https://doi.org/10.1242/dmm.028910
-
From A to Z and back? multicompartment proteins in the sarcomereTrends in Cell Biology 16:11–18.https://doi.org/10.1016/j.tcb.2005.11.007
-
Beta 3 tubulin expression characterizes the differentiating mesodermal germ layer during Drosophila embryogenesisDevelopment 104:525–531.
-
Mechanical forces during muscle developmentMechanisms of Development 144:92–101.https://doi.org/10.1016/j.mod.2016.11.003
-
In vivo imaging of Muscle-tendon morphogenesis in Drosophila pupaeJournal of Visualized Experiments: JoVE e57312.https://doi.org/10.3791/57312
-
The Sequence Alignment/Map format and SAMtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
Sallimus and the dynamics of sarcomere assembly in Drosophila flight musclesJournal of Molecular Biology 427:2151–2158.https://doi.org/10.1016/j.jmb.2015.04.003
-
Myosin isoform switching during assembly of the Drosophila flight muscle thick filament latticeJournal of Cell Science 126:139–148.https://doi.org/10.1242/jcs.110361
-
Regulation of skeletal muscle sarcomere integrity and postnatal muscle function by Mef2cMolecular and Cellular Biology 27:8143–8151.https://doi.org/10.1128/MCB.01187-07
-
Contractile properties of developing human fetal cardiac muscleThe Journal of Physiology 594:437–452.https://doi.org/10.1113/JP271290
-
Ultrastructure of developing flight muscle in Drosophila. I. assembly of myofibrilsDevelopmental Biology 160:443–465.https://doi.org/10.1006/dbio.1993.1320
-
Flightin is essential for thick filament assembly and sarcomere stability in Drosophila flight musclesThe Journal of Cell Biology 151:1483.https://doi.org/10.1083/jcb.151.7.1483
-
Contribution of sequence variation in Drosophila actins to their incorporation into actin-based structures in vivoJournal of Cell Science 118:3937–3948.https://doi.org/10.1242/jcs.02517
-
Assembly and maintenance of myofibrils in striated muscleHandbook of Experimental Pharmacology 235:39–75.https://doi.org/10.1007/164_2016_53
-
Myofibrillogenesis in skeletal muscle cells in zebrafishCell Motility and the Cytoskeleton 66:556–566.https://doi.org/10.1002/cm.20365
-
Fiber types in mammalian skeletal musclesPhysiological Reviews 91:1447–1531.https://doi.org/10.1152/physrev.00031.2010
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
The him gene inhibits the development of Drosophila flight muscles during metamorphosisMechanisms of Development 126:595–603.https://doi.org/10.1016/j.mod.2009.03.003
-
The initial steps of myofibril assembly: integrins pave the wayNature Reviews Molecular Cell Biology 10:293–298.https://doi.org/10.1038/nrm2634
-
Muscle LIM proteins are associated with muscle sarcomeres and require dMEF2 for their expression during Drosophila myogenesisMolecular Biology of the Cell 10:2329–2342.https://doi.org/10.1091/mbc.10.7.2329
-
How to build fast muscles: synchronous and asynchronous designsIntegrative and Comparative Biology 42:762–770.https://doi.org/10.1093/icb/42.4.762
-
Myocyte enhancer factor 2 and chorion factor 2 collaborate in activation of the myogenic program in DrosophilaMolecular and Cellular Biology 28:1616–1629.https://doi.org/10.1128/MCB.01169-07
-
Titin: properties and family relationshipsNature Reviews Molecular Cell Biology 4:679–689.https://doi.org/10.1038/nrm1198
-
Muscle Development in Drosophila143–156, Molecular Basis of Muscle Structure, Muscle Development in Drosophila, New York, NY, Springer New York.
-
Nemaline myopathiesSeminars in Pediatric Neurology 18:230–238.https://doi.org/10.1016/j.spen.2011.10.004
-
Mechanisms of skeletal muscle injury and repair revealed by gene expression studies in mouse modelsThe Journal of Physiology 582:825–841.https://doi.org/10.1113/jphysiol.2007.132373
-
The Split-Apply-Combine strategy for data analysisJournal of Statistical Software 40:1–29.https://doi.org/10.18637/jss.v040.i01
-
SALL4: an emerging cancer biomarker and targetCancer Letters 357:55–62.https://doi.org/10.1016/j.canlet.2014.11.037
Article and author information
Author details
Funding
Max-Planck-Gesellschaft
- Maria L Spletter
- Christiane Barz
- Assa Yeroslaviz
- Xu Zhang
- Sandra B Lemke
- Giovanni Cardone
- Bianca H Habermann
- Frank Schnorrer
European Molecular Biology Organization (EMBO-LTR 688-2011)
- Maria L Spletter
Alexander von Humboldt-Stiftung
- Maria L Spletter
National Institute for Health Research (5F32AR062477)
- Maria L Spletter
Centre National de la Recherche Scientifique
- Xu Zhang
- Adrien Bonnard
- Bianca H Habermann
- Frank Schnorrer
Boehringer Ingelheim Fonds
- Sandra B Lemke
Center for Integrated Protein Science München (the Boehringer Ingelheim Fonds)
- Sandra B Lemke
Aix-Marseille Université (ANR-11-IDEX-0001-02)
- Bianca H Habermann
- Frank Schnorrer
H2020 European Research Council (ERC Grant 310939)
- Frank Schnorrer
Agence Nationale de la Recherche (ANR-11- LABX-0054)
- Frank Schnorrer
European Molecular Biology Organization (EMBO-YIP)
- Frank Schnorrer
Agence Nationale de la Recherche (ANR-10-INBS-04- 01)
- Frank Schnorrer
Agence Nationale de la Recherche (ANR ACHN)
- Frank Schnorrer
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Irene Ferreira for constructing the tagged kon-tiki allele and the Bloomington and VDRC stock centers for fly stocks. We are grateful to Reinhard Fässler and Andreas Ladurner for generous support and to Bettina Stender for excellent technical assistance. We thank Elena Nikonova for initial analysis of the Act88F-GAL4 x salmIR phenotype and Sandra Esser for testing the salm KK181052 hairpin. We acknowledge the VBCF (Vienna, AT) for mRNA-sequencing and the Core Facility Bioimaging at the Max Planck Institute for Biochemistry and the LMU Biomedical Center (Martinsried, DE) for help with confocal image analysis. We thank Florian Marty for help with the BioSorter, and Alexander Stark for help setting-up the mRNA-sequencing. We thank Aynur Kaya-Copur and Wouter Koolhaas for helpful discussions, and Nuno Luis and Vincent Loreau for insightful comments on the manuscript. Our work was supported by the Max Planck Society and the CNRS; postdoctoral Humboldt, EMBO long-term (688–2011), and NIH-NRSA (5F32AR062477) fellowships (MLS), the Frederich-Bauer Stiftung (MLS), the Center for Integrated Protein Science München (MLS), the Boehringer Ingelheim Fonds (SBL), a Career Development Award from the Human Frontier Science Program (FS), the EMBO Young Investigator Program (FS) and the European Research Council under the European Union's Seventh Framework Programme (FP/2007–2013)/ERC Grant 310939 (FS), the excellence initiative Aix-Marseille University AMIDEX (ANR-11-IDEX-0001–02) (FS), the ANR-ACHN (FS) and the LabEX-INFORM (ANR-11- LABX-0054) (FS) and France Bioimaging (ANR-10-INBS-04–01) (FS).
Copyright
© 2018, Spletter et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,411
- views
-
- 536
- downloads
-
- 78
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Cell Biology
Expansion microscopy (ExM) enables nanoscale imaging using a standard confocal microscope through the physical, isotropic expansion of fixed immunolabeled specimens. ExM is widely employed to image proteins, nucleic acids, and lipid membranes in single cells; however, current methods limit the number of samples that can be processed simultaneously. We developed High-throughput Expansion Microscopy (HiExM), a robust platform that enables expansion microscopy of cells cultured in a standard 96-well plate. Our method enables ~4.2 x expansion of cells within individual wells, across multiple wells, and between plates. We also demonstrate that HiExM can be combined with high-throughput confocal imaging platforms to greatly improve the ease and scalability of image acquisition. As an example, we analyzed the effects of doxorubicin, a known cardiotoxic agent, on human cardiomyocytes (CMs) as measured by the Hoechst signal across the nucleus. We show a dose-dependent effect on nuclear DNA that is not observed in unexpanded CMs, suggesting that HiExM improves the detection of cellular phenotypes in response to drug treatment. Our method broadens the application of ExM as a tool for scalable super-resolution imaging in biological research applications.
-
- Cell Biology
- Developmental Biology
Eukaryotic cells depend on exocytosis to direct intracellularly synthesized material toward the extracellular space or the plasma membrane, so exocytosis constitutes a basic function for cellular homeostasis and communication between cells. The secretory pathway includes biogenesis of secretory granules (SGs), their maturation and fusion with the plasma membrane (exocytosis), resulting in release of SG content to the extracellular space. The larval salivary gland of Drosophila melanogaster is an excellent model for studying exocytosis. This gland synthesizes mucins that are packaged in SGs that sprout from the trans-Golgi network and then undergo a maturation process that involves homotypic fusion, condensation, and acidification. Finally, mature SGs are directed to the apical domain of the plasma membrane with which they fuse, releasing their content into the gland lumen. The exocyst is a hetero-octameric complex that participates in tethering of vesicles to the plasma membrane during constitutive exocytosis. By precise temperature-dependent gradual activation of the Gal4-UAS expression system, we have induced different levels of silencing of exocyst complex subunits, and identified three temporarily distinctive steps of the regulated exocytic pathway where the exocyst is critically required: SG biogenesis, SG maturation, and SG exocytosis. Our results shed light on previously unidentified functions of the exocyst along the exocytic pathway. We propose that the exocyst acts as a general tethering factor in various steps of this cellular process.