Myofibril diameter is set by a finely tuned mechanism of protein oligomerization in Drosophila
Abstract
Myofibrils are huge cytoskeletal assemblies embedded in the cytosol of muscle cells. They consist of arrays of sarcomeres, the smallest contractile unit of muscles. Within a muscle type, myofibril diameter is highly invariant and contributes to its physiological properties, yet little is known about the underlying mechanisms setting myofibril diameter. Here we show that the PDZ and LIM domain protein Zasp, a structural component of Z-discs, mediates Z-disc and thereby myofibril growth through protein oligomerization. Oligomerization is induced by an interaction of its ZM domain with LIM domains. Oligomerization is terminated upon upregulation of shorter Zasp isoforms which lack LIM domains at later developmental stages. The balance between these two isoforms, which we call growing and blocking isoforms sets the stereotyped diameter of myofibrils. If blocking isoforms dominate, myofibrils become smaller. If growing isoforms dominate, myofibrils and Z-discs enlarge, eventually resulting in large pathological aggregates that disrupt muscle function.
eLife digest
Muscles are made up of many long muscle fibers, each containing thousands of cylindrical segments called sarcomeres. When animals move, proteins in the sarcomere move past each other, shortening the muscles. Inside each muscle, all sarcomeres have the same length and diameter. The protein titin controls the length of each sarcomere, but it was unknown what controls the diameter.
At the end of each sarcomere is a structure called the Z-disc that is composed of many muscle proteins. Mutations in Z-disc proteins are often involved in diseases called myopathies, where muscle structure breaks down. As the size of the Z-disc determines sarcomere diameter, improper regulation of sarcomere diameter could contribute to myopathies. One Z-disc protein called Zasp is a candidate for controlling diameter and can have many different forms in the same cells. Zasp has a similar role in most animals including humans, mice and flies.
González-Morales et al. investigated Zasp in the muscles of the fruit fly, Drosophila melanogaster. Gene editing was used to vary the amounts of different forms of Zasp inside the muscles. The results revealed two types of Zasp, those that make sarcomeres wider, and those that limit growth. Reducing the second type of Zasp resulted in bigger Z-discs and in muscle aggregates similar to the ones seen in patients with certain myopathies.
This study reveals a mechanism for coordinating the development of muscle. It also reveals the likely cause of certain myopathies and suggests a possible target for future treatment through regulation of Zasp proteins.
Introduction
Inside cells, proteins are assembled into complex functional units. The correct assembly of these units is crucial for their function (Marsh and Teichmann, 2015). Myofibrils are highly organized assemblies of cytoskeletal proteins forming an array of sarcomeres that are embedded in the cytosol of myotubes and mediate contractility (Huxley and Niedergerke, 1954b; Huxley and Hanson, 1954a; Huxley, 2004; Lemke and Schnorrer, 2017). Sarcomeres are composed of actin-containing thin filaments and myosin-containing thick filaments arranged into antiparallel cables. Thin filaments are anchored to a large multiprotein complex called the Z-disc (Luther, 2009), and thick filaments are anchored to another large multiprotein complex called the M-line (Agarkova and Perriard, 2005). Anchoring of myofibrils to the exoskeleton provides the mechanical tension that aligns sarcomeres into myofibrils and coordinates their development (Weitkunat et al., 2017; Weitkunat et al., 2014). Once aligned, sarcomeres grow to their final size (Lemke and Schnorrer, 2017). Electron and confocal microscopy studies showed that sarcomeres form initially from small structures called Z-bodies that grow eventually into mature Z-discs to which thin filaments are anchored (Loison et al., 2018; Orfanos et al., 2015; Reedy and Beall, 1993; Shafiq, 1963; Sparrow and Schöck, 2009). The size of the Z-disc therefore sets the diameter of the myofibril (Agarkova and Perriard, 2005; Luther, 2009). While mechanisms have been proposed that set the length of sarcomeres (Fernandes and Schöck, 2014; Gokhin and Fowler, 2013; Tskhovrebova and Trinick, 2017), Z-disc growth and growth termination is poorly understood.
A hallmark of genetically caused myopathies is the appearance of large aggregates composed mainly of Z-disc proteins (Kley et al., 2016; Maerkens et al., 2016). Interestingly, many myopathy-associated mutations encode Z-disc proteins. Mutations in any of the four α-actinin genes or in Zasp and other Alp/Enigma family genes in humans cause myopathies characterized by the presence of large aggregates (Murphy and Young, 2015; Selcen and Engel, 2005). The aggregation phenotype is conserved among animals: fruit flies with mutations in myopathy-related genes also develop Z-disc aggregates (González-Morales et al., 2017).
α-Actinin and Z-disc Alternatively Spliced Protein (Zasp) are conserved proteins that coordinate Z-disc formation (Faulkner et al., 1999; Katzemich et al., 2013; Murphy and Young, 2015). α-Actinin forms a rod-shaped antiparallel homodimer at the Z-disc, where it crosslinks and serves as an attachment point for actin filaments (Djinović-Carugo et al., 1999; Luther, 2009; Murphy and Young, 2015; Ribeiro et al., 2014; Rusu et al., 2017; Takahashi and Hattori, 1989). Zasp and other members of the Alp/Enigma family of proteins are scaffolding proteins with an α-actinin-binding PDZ domain (InterPro: IPR001478), an uncharacterized Zasp Motif (ZM; InterPro: IPR031847 and IPR006643) domain, and one to four protein-protein interaction LIM domains (InterPro: IPR001781) (Finn et al., 2017; Klaavuniemi et al., 2004; Liao et al., 2016). Zasp and α-actinin proteins are present at the earliest stages of Z-disc formation, and are required for Z-disc assembly (Dabiri et al., 1997; Katzemich et al., 2013).
Vertebrates have seven Alp/Enigma genes, each encoding several isoforms. The Drosophila genome has three Zasp genes, Zasp52, Zasp66, and Zasp67, which encode 21, 12, and 4 isoforms, respectively (Gramates et al., 2017). Zasp66 and Zasp67 are duplications of Zasp52 and resemble the smallest isoforms of Zasp52 (González-Morales et al., 2019). The number of isoform variants and genes adds an additional layer of complexity and regulation to sarcomere formation.
We used Drosophila indirect flight muscles and asked if the mechanism that controls Z-disc size relates to the pathological aggregation behavior known for Z-disc-related myopathies. We show that accumulation of multivalent Zasp growing isoforms (with multiple LIM domains) causes Z-disc growth, whereas upregulation of monovalent blocking isoforms at later developmental time points terminates Z-disc growth. An imbalance of growing and blocking Zasp isoforms results either in aggregate formation, enlarged Z-disc size or reduced Z-disc size. We propose that this mechanism has wide implications for diseases caused by aggregate formation.
Results
Aggregate formation upon Zasp overexpression
Dominant Zasp mutations cause aggregate formation in human myopathies (Selcen and Engel, 2005). To examine if this holds true in Drosophila, we overexpressed a full-length Zasp52-PR transgene consisting of a PDZ, a ZM, and 4 LIM domains in the indirect flight muscles (IFM) of adult flies, which causes formation of large aggregates (Figure 1A,B,I, Figure 1—figure supplement 1). To determine which domain is responsible for aggregation, we deleted them individually: the absence of LIM domains (Zasp52-PP) and ZM domain (Zasp52-PR∆ZM) abolished aggregate formation, whereas removal of the PDZ domain only slightly reduced the number of aggregates (Zasp52-PR∆PDZ, Figure 1C–E,I). The PDZ domain is required for binding α-actinin, a crucial crosslinker of actin filaments at the Z-disc (Katzemich et al., 2013; Liao et al., 2016). Reducing α-actinin levels by RNAi in Zasp52-PR overexpression reduces aggregate formation to a similar level (Figure 1—figure supplement 2), indicating that α-actinin increases aggregate number, but is not required for their initial formation. Expression of the paralogous genes Zasp66 and Zasp67 (González-Morales et al., 2019), which lack LIM domains like Zasp52-PP, also does not result in aggregate formation (Figure 1F,G and I). Finally, we tested if multiple LIM domains are required for aggregate formation by overexpressing Zasp52-PK consisting of PDZ, ZM and 1 LIM domain. No aggregates form (Figure 1H and I) indicating that the ZM domain and multiple LIM domains are required for aggregate formation.
To test if aggregates formed as a result of Z-disc overgrowth, we tested the overexpression of Zasp52-PR using UH3-Gal4, a driver line that has the same temporal and spatial expression pattern as Act88-Gal4 but is expressed at a much lower level. To better capture the size variation, we measured the Z-disc height that corresponds to the disc diameter and categorized them into specific size categories. In this condition, Z-discs are bigger than the control, but aggregates are not present (Figure 1J–L), indicating that Z-disc growth and pathological aggregation correlate with the amount of Zasp52 protein.
Domains required for Zasp self-interaction
Next, we employed the yeast two-hybrid system (Y2H) to determine if Zasp forms oligomers by interacting with itself. Zasp52-PK and Zasp52-PE, another full-length isoform, can interact with each other, whereas controls do not interact (Figure 2A). Full-length isoforms of Zasp66 and Zasp67 (Zasp66-PK and Zasp67-PD) can also interact with Zasp52-PK and Zasp52-PE (Figure 2A). These data suggest that Zasp proteins can homo- and heterodimerize.
We confirmed these results by co-immunoprecipitation of Flag-tagged Zasp52 with GFP-tagged Zasp52 from thorax muscle extracts (Figure 2—figure supplement 1A). Furthermore, bacterially purified His-Zasp52-PK-Flag dimerized in vitro in a chemical crosslinking assay (Figure 2—figure supplement 1B).
To identify the domains involved in homo- and heterodimerization, we tested the interaction of Zasp52-PK, Zasp66 and Zasp67 with each individual domain of Zasp52. Zasp52-PK interacts with the ZM domain and LIM domains, whereas Zasp66 and Zasp67 only interact with LIM domains (Figure 2B). This suggests that dimerization is mediated by a ZM-LIM domain interaction, and because Zasp66 and Zasp67 lack LIM domains, they cannot interact with the ZM domain (Figure 2C). To confirm this hypothesis, we next tested the interaction of LIM2A, LIM2B, LIM3 and Zasp52-PK with a series of Zasp66 deletion constructs. As soon as the ZM domain of Zasp66 is deleted, the interaction with LIM domains is abolished (Figure 2D). Finally, we tested Zasp66-PH, a ZM-only isoform of Zasp66 for interaction with the individual domains of Zasp52. Zasp66-PH can only interact with LIM domains of Zasp52 (Figure 2E). Thus, LIM-ZM binding mediates homo- and heterodimerization between Zasp proteins.
Bimolecular fluorescence complementation reveals Zasp interaction at the Z-disc
Next, we investigated if LIM-ZM interaction observed in a heterologous system also occurs in vivo at the Z-disc. We hypothesized that an endogenous GFP-tagged Zasp66 will be recruited to the Z-discs by Zasp52. We analyzed two Zasp52 mutants that differentially affect sarcomere structure: Zasp52MI02988 disrupts N-terminal isoforms and only partially affects the last three LIM domains, whereas Zasp52MI00979 introduces a stop codon before the last three LIM domains (Figure 3—figure supplement 1A) (Liao et al., 2016). Zasp66-GFP fluorescence is mildly reduced in Zasp52MI02988 mutants (Figure 3A,B and D). In contrast, Zasp66-GFP fluorescence is strongly reduced in Zasp52MI00979 mutants (Figure 3C and D), suggesting that Zasp52 recruits other Zasp proteins through its LIM domains, consistent with the Y2H results. α-Actinin localization at the Z-disc was not decreased in any of the Zasp mutants (Figure 3—figure supplement 1B,C).
To determine if LIM and ZM domains directly interact in flies at the Z-disc, we employed bimolecular fluorescence complementation (BiFC) assays (Figure 3—figure supplement 1D) (Ciruela, 2008; Gohl et al., 2010). We first confirmed that the ZM-only isoform Zasp66-PH fused to a Venus tag localized to Z-discs (Figure 3—figure supplement 1E). We then generated transgenic flies expressing Zasp52-PK and Zasp66-PH tagged with either the C-terminus or the N-terminus of yellow fluorescent protein (NYFP or CYFP) and quantified the fluorescence of reconstituted YFP at the Z-disc. Zasp52-PK-NYFP and Zasp52-PK-CYFP show specific fluorescence at the Z-disc, whereas controls do not (Figure 3E–G). More importantly, the ZM-only isoform Zasp66-PH-NYFP interacts specifically with Zasp52-PK-CYFP at the Z-disc (Figure 3H). Thus, ZM-LIM interaction of Zasp proteins occurs in vivo at the Z-disc.
ZM domains and canonical Z-discs co-appeared in bilateral animals
LIM domains are well-known protein-protein interaction domains (Kadrmas and Beckerle, 2004). In contrast, functional information on the ZM domain is scarce. It is an uncharacterized short motif of 26 amino acids found only in Alp/Enigma proteins (Klaavuniemi et al., 2004; Letunic et al., 2015). Alignment of Drosophila Zasp-encoding genes to the ZM consensus sequence shows weak conservation (Figure 4A). To determine the phylogenetic distribution of the ZM domain, we compiled the presence of ZM domains in all branches of metazoans using data from the PFAM database (Bateman et al., 2004; El-Gebali et al., 2019). In contrast to the universal presence of LIM domains in eukaryotes, the ZM domain is only present in a subset of animal genomes and absent from plants, fungi and unicellular metazoans (Figure 4B). Muscle striation evolved independently in cnidarians and bilaterians, with only the latter showing a canonical Z-disc structure (Steinmetz et al., 2012). Intriguingly, the ZM domain is restricted to bilateral animals (Figure 4B). Thus, our correlative data suggest that canonical Z-discs and ZM domains arose together during evolution.
Differential temporal expression and localization of Zasp growing and blocking isoforms
Zasp isoforms can be divided into isoforms with multiple LIM domains (e.g. Zasp52-PR) and isoforms with just one or no LIM domain (e.g. Zasp52-PK, Zasp52-PP, Zasp66, Zasp67). We named the former growing isoforms and the latter blocking isoforms. We then considered what mechanism might coordinate Z-disc growth: we hypothesize that multiple LIM domain-containing Zasp proteins recruit other Zasp proteins by interacting with their ZM domains. If the recruited proteins are also multivalent (contain two or more LIM domains), more Zasp proteins will be recruited leading to Z-disc growth. However, if recruited proteins lack LIM domains, they block further recruitment of Zasp proteins, and Z-disc growth terminates (Figure 5A).
To regulate Z-disc growth, growing isoforms should be overrepresented at earlier stages of Z-disc formation, and blocking isoforms at later stages, to stop Z-disc growth. We used an RNAseq dataset from developing IFM that covers the whole Z-disc formation process, from 16 hr after puparium formation (APF) to newly eclosed flies (Spletter et al., 2018). We observed earlier expression of the growing isoforms compared to the blocking isoforms (Figure 5B). For example, multivalent Zasp52 growing isoforms are already strongly expressed at 24 hr APF, whereas some blocking isoforms without any LIM domains are expressed only after 60 hr APF (Figure 5B).
Z-discs grow at the periphery of the disc, starting as a small Z-body (Orfanos et al., 2015; Shwartz et al., 2016). Our model therefore predicts growing isoforms to be enriched at the centre and blocking isoforms to be enriched at the periphery of the final-sized Z-disc. To analyze the distribution of Zasp proteins within Z-discs, we made cross sections of myofibrils and evaluated the localization pattern of blocking isoforms (using Zasp66-GFP and Zasp67-GFP), and growing isoforms (using Zasp52-GFP; Zasp52ZCL423) at the level of the Z-disc. To label the entire disc, we used actin staining as counterstain (Figure 5C). To compensate for the low signal-to-noise ratio, we used a smoothening algorithm that takes advantage of the geometrical properties of discs (Figure 5—figure supplement 1). Actin is evenly distributed in the Z-disc (Figure 5D,E). The growing isoforms are present throughout the disc, but most concentrated at the centre of the Z-disc (Figure 5D,E). The blocking isoforms Zasp66 and Zasp67 are partially excluded from the centre of the disc and form a ring-like pattern (Figure 5D,E). Zasp66 localizes more at the periphery compared to Zasp67, which correlates with its later expression peak compared to Zasp67 (middle panels of Figure 5B).
The balance between growing and blocking isoforms controls Z-disc growth and aggregation
We then asked whether decreasing the levels of growing isoforms would result in smaller Z-discs. We measured individual Z-discs in control and Zasp52 mutants (Figure 6A–C). Interestingly, only the Zasp52MI00979 mutant, which deletes multiple LIM domains, had smaller Z-discs compared to the control (Figure 6B,C). In Zasp52MI00979, the smaller size categories were significantly elevated, whereas Zasp52MI02988 was comparable to the control (Figure 6C).
We next determined if the opposite approach, increasing the levels of blocking isoforms, also results in smaller Z-discs. Overexpressing the Zasp52-PP blocking isoform led to smaller Z-discs than the control (Figure 6D). This effect requires the ZM domain, because overexpression of Zasp52-Stop143 containing only a functional extended PDZ domain, causes a milder phenotype than Zasp52-PP (Figure 6D, Figure 1—figure supplement 1). Like Zasp52-PP, overexpression of Zasp66 or Zasp67 resulted in smaller Z-discs (Figure 6E). We also confirmed these phenotypes by measuring the relative fluorescence of Zasp52-mCherry at the Z-disc. Upon overexpression of Zasp52-PP, but not of Zasp52-Stop143, Zasp52-mCherry recruitment to the Z-disc is reduced (Figure 6F, Figure 6—figure supplement 1A). Likewise, overexpression of Zasp66 and Zasp67 reduces Zasp52-mCherry recruitment to the Z-disc (Figure 6G, Figure 6—figure supplement 1A). In contrast, α-actinin levels were unchanged in Zasp overexpression conditions (Figure 6—figure supplement 1B,C). These data indicate that increasing the ZM/LIM ratio in IFM results in smaller Z-discs.
Finally, we wondered if depleting the blocking isoforms encoded by Zasp66 and Zasp67 can generate bigger Z-discs or aggregates using CRISPR Zasp66 and Zasp67 null mutants (González-Morales et al., 2019). We observed rare aggregates in Zasp66 and Zasp67 single mutants (Figure 7A,B,D). As these two proteins likely have some redundant roles, we also analyzed the Zasp66 Zasp67 double mutant, where we observed frequent aggregates (Figure 7C,D). Additionally, the single mutants had normal-sized Z-discs, whereas the double mutant had enlarged Z-discs (Figure 7E). Thus, blocking isoforms are required to prevent Z-disc overgrowth.
As many diseases are believed to be caused by the formation of aggregates (Baba et al., 1998; Selcen, 2008), we asked if we can suppress aggregate formation in our model of overexpressed GFP-Zasp52-PR (Figure 1B). Co-overexpression of both growing and the Zasp52-PP blocking isoform suppressed aggregate formation compared to a LacZ control (Figure 7F–I). The mutated blocking isoform Zasp52-Stop143 without the ZM domain did not suppress aggregation (Figure 7H,I). Co-overexpression of Zasp66 or Zasp67 also suppressed aggregation (Figure 7I).
Discussion
How myofibrils assemble and grow until they reach their final, highly invariant size has remained poorly understood. Here, we reveal a conserved mechanism by which the Z-disc scaffold protein Zasp oligomerizes to induce Z-disc and thereby myofibril growth. Furthermore, we show that Z-disc and sarcomere growth is terminated when shorter isoforms of Zasp are expressed that block multivalency-driven oligomerization.
Z-disc formation and growth is driven by multivalent oligomerization of Zasp proteins
Our findings indicate that Z-disc formation and growth is driven by multivalent oligomerization of Zasp proteins with multiple LIM domains, and eventually terminated at the proper Z-disc size by the upregulation of blocking isoforms without LIM domains (see model in Figure 5A). The process of Z body and Z-disc formation is reminiscent of membraneless organelles with compositions distinct from the surrounding cytosol, which form through a mechanism of phase separation (e.g. Cajal bodies or P bodies) (Boeynaems et al., 2018; Weber, 2017). Both have sharp boundaries between themselves and the cytoplasm, they form and organize as discrete puncta in the cytosol, and multivalent protein domains are often involved in their formation (Li et al., 2012; Weber, 2017).
Sarcomere size is stereotyped in a given muscle type but distinct among different muscles (Schönbauer et al., 2011). How can our model explain differences in sarcomere sizes? The sarcomere grows, while the Zasp growing isoforms dominate the Zasp isoform pool. Different sarcomere sizes can be achieved in two ways. First, the sarcomere growth period – the window of time in which Zasp growing isoforms dominate, might vary among muscle types. Second, the speed at which new Zasp molecules are recruited to the Z-disc, might be different among muscle types, while the growth period remains constant. Given the diversity of muscle types and therefore sarcomere sizes that exist, it is likely that a combination of these two strategies occur simultaneously.
Finally, apart from the ZM-LIM mechanism described here, additional redundant mechanisms to control Z-disc growth might exist, as evidenced by the observation that Zasp52-PR overexpression makes big Z-discs and aggregates, while the mutant removing Zasp52 LIM domains reduces Z-disc size to a comparatively small degree. Redundant mechanisms might operate through other LIM domain proteins, or the coordination of Z-disc and M-line growth, all of which may provide important buffering functions to ensure proper myofibril size, which is crucial for fully functional muscles.
The role of the ZM domain
The ZM domain is a conserved domain without a clearly defined function. On its own, the ZM domain from two mouse Zasp proteins localizes to the Z-disc (Klaavuniemi et al., 2004; Klaavuniemi and Ylänne, 2006). Our data suggest that Z-disc localization is a conserved feature of ZM domains from vertebrates to insects. ZM-containing proteins are tethered to the Z-disc by the physical interaction with the LIM domains of other Zasp proteins. In sum, the LIM domain serves as a recruitment signal for Zasp proteins and potentially other unidentified ZM-containing proteins to join the Z-disc. In addition, given the appearance of the ZM domain in bilateral animals with canonical Z-discs, we postulate that a conserved mechanism involving LIM-ZM binding underlies Z-disc growth and growth termination.
Blocking and growing isoforms in vertebrates
In vertebrates, the Zasp proteins are very diverse and are better known as Alp/Enigma family: ZASP/Cypher/Oracle/LDB3/PDLIM6, ENH/PDLIM5, PDLIM7/ENIGMA/LMP-1, CLP36/PDLIM1/Elfin/hCLIM1, PDLIM2/Mystique/SLIM, ALP/PDLIM3, and RIL/PDLIM4. The ZM/DUF4749 motif occurs in ZASP, CLP36, PDLIM2, ALP and RIL (Cheng et al., 2010; D'Cruz et al., 2016; Faulkner et al., 1999; Vallenius et al., 2004; Zheng et al., 2010; Zhou et al., 2001). The LIM domain occurs in all Zasp proteins, either as one domain in Alp family members or as three domains in Enigma family members (Zheng et al., 2010). We identified two Zasp genes that encode only blocking isoforms in fruit flies: Zasp66 and Zasp67, and one gene, Zasp52, that encodes blocking and growing isoforms. Although Zasp66 and Zasp67 genes are insect-specific (González-Morales et al., 2019), vertebrate Alp/Enigma genes also express isoforms without LIM domains that could fulfill a blocking isoform function (Cheng et al., 2011; Zheng et al., 2010). In addition, because Zasp52-PK, which only contains one LIM domain, behaves as a blocking isoform, the Alp members with only one LIM domain might also behave as blocking isoforms.
The function of the growing isoforms of Zasp requires multiple functional LIM domains. As the Enigma family members contain three C-terminal LIM domains, they are the ideal candidates to fulfill the growing role in vertebrates. Three Enigma proteins exist in vertebrates: PDLIM7/Enigma/LMP-1, ENH/PDLIM5 and ZASP/Cypher/Oracle/LDB3/PDLIM6. Functional redundancy between them at the Z-disc is likely common and demonstrated in one case (Mu et al., 2015). In Cypher knockout mice sarcomere assembly occurs normally during development, followed by immediate sarcomere failure after postnatal onset of contractility (Zhou et al., 2001). ENH mutants exhibit cardiac dilation and abnormal Z-disc structure in the heart (Cheng et al., 2010). Intriguingly, in both Cypher and ENH single mutants, as well as Cypher ENH double mutants, sarcomeres look considerably smaller in diameter in electron microscopy images (Cheng et al., 2010; Mu et al., 2015; Zhou et al., 2001). Thus, a similar role for Enigma proteins in setting sarcomere diameter in vertebrates appears likely.
Sarcomere size control in human myopathies
Is the mechanism that controls Z-disc size related to the protein aggregation defects in human myopathies? Our Z-disc oligomerization hypothesis agrees well with the observation that many myopathies present aggregates, and several human ZASP mutations have been linked to aggregate-forming myopathies (Murphy and Young, 2015; Selcen and Engel, 2005). Many ZASP mutations linked to disease lie within the ZM domain or one of the LIM domains (Selcen and Engel, 2005; Theis et al., 2006; Vatta et al., 2003). Protein aggregation in myopathy patients might be a consequence of an imbalance in the mechanism that controls sarcomere size, favoring the growing over the blocking isoforms. If this were the case, our data points to a potential therapeutic avenue: blocking the growing isoforms with short peptides containing a ZM domain.
In conclusion, we propose that a conserved mechanism involving LIM-ZM binding underlies Z-disc growth and therefore myofibril diameter.
Materials and methods
Experimental model
Request a detailed protocolWe used Drosophila melanogaster as a model organism for most experiments. Fly stocks and crosses were raised at 25°C. A comprehensive list of all strains used and generated can be found in the Key Resources Table. We used Saccharomyces cerevisiae for the yeast two-hybrid assays. Escherichia coli BL-21 strain was used to express recombinant proteins.
Transgenic flies were generated by site-directed integration using the PhiC31 integrase method into either M[3xP3-RFP.attP]ZH-58A- or M[3xP3-RFP.attP]ZH-86Fb-bearing flies to ensure comparable expression levels (Bischof et al., 2007; Markstein et al., 2008). Unless stated otherwise, Act88F-Gal4 was used to drive strong transgene expression in the IFM (Bryantsev et al., 2012).
Yeast two hybrid and protein crosslinking assays
Request a detailed protocolPartial or complete coding sequences for Zasp52, Zasp66 and Zasp67 were cloned into Y2H vectors either by PCR-ligase cloning or through Gateway cloning. Constructs and cloning details are listed in the Key Resources Table. All constructs were verified by sequencing. Selected constructs were transformed into the Matchmaker Y2H Gold strain using the lithium acetate method (Gietz and Schiestl, 2007). Double transformant colonies were selected and amplified in media lacking leucine and tryptophan (-Leu/-Trp). Serial dilutions of the selected double transformants were grown in plates lacking leucine, tryptophan, histidine and adenine (-Ade/-His/- Leu/-Trp) to test for protein-protein interactions. The experiments were done at least three times.
Recombinant 6xHis-Zasp52-PK-FLAG was expressed in Escherichia coli BL21 bacteria. Then, 6xHis-Zasp52-PK-FLAG was purified from the protein extracts using Ni-NTA agarose beads (Qiagen) for 3 hr at 4°C. The purified protein was then dialyzed overnight at 4°C. Finally, purified 6xHis-Zasp52-PK-FLAG was either incubated with ethylene glycol bis-sulfosuccinimidyl succinate (EGS) or alone. Then, the protein samples were analyzed by denaturing SDS-PAGE followed by western blotting with anti-Flag antibody (1:5000).
Muscle staining and microscopy
Request a detailed protocolThe IFM were dissected as previously described (González-Morales et al., 2017; Xiao et al., 2017). All images were acquired with comparable parameters. In agreement with previous studies, at least 10 thoraces were dissected for each condition. Samples were allocated into experimental groups according to their genotype. A big number of flies were collected for each experimental group and a subsample was randomly selected for dissection.
Bimolecular fluorescence complementation assay
View detailed protocolZasp52-PK, Zasp52-PE and Zasp66-PH Gateway pENTRY clones were obtained from the DGRC (GEO02280, GEO12859 and GEO14752). The BIFC pDEST vectors pUAS-RfB-HA-CYFP-attB and pUAS-RfB-myc-NYFP-attB vectors were a kind gift from Sven Bogdan (Gohl et al., 2010). Constructs were cloned using Gateway technology, and the final vectors were verified by sequencing. The Key Resources Table contains a list of transgenic flies used in this study.
The Act88F-Gal4 driver line was used to express CYFP- and NYFP-tagged proteins in the IFM. Then, the YFP fluorescence of individual Z-discs was quantified using the ImageJ plot profile tool (Schindelin et al., 2012). At least 20 samples were used for each condition. The data was normalized to the basal noise levels and plotted in R software.
Image analyses
Request a detailed protocolTo estimate the size of individual Z-discs and the fluorescence intensity levels of Z-disc proteins we first used a segmentation method that allows individual Z-discs to be measured (Xiao et al., 2017). Briefly, images from fluorescent Z-disc proteins were passed through a threshold filter, and then individual Z-discs were separated and measured using the analyze particles tool. Both tools can be found in any recent ImageJ/Fiji distribution (Schindelin et al., 2012). To estimate the Z-disc size, we measured the Z-disc height, that corresponds to the disc diameter. Then, the data was analyzed and plotted in R software. To estimate the relative fluorescence intensities of Z-disc proteins, we first obtained the mean intensity raw values for each Z-disc. Then, we normalized the values to a control genotype. The resulting data was then analyzed in R software. Rare segmentation mistakes result from Z-discs that are in close proximity. We could filter these mistakes out by size exclusion in the control images, but not in mutant conditions. We decided to leave the segmentation mistakes for consistency. If anything, this method underestimates the size differences between control and mutants.
Batch Macro processing code for measuring individual Z-disc side-views in ImageJ
setAutoThreshold("Default dark");
run("Analyze Particles...", "size = 0.2 Infinity show = Outlines display exclude");
To estimate the density of aggregates, present in a given IFM sample, we used confocal images of GFP- or mCherry-tagged Zasp. Each 36 × 36 µm image was divided into 256 2 × 2 µm images. An automatic threshold filter based on the original image was then used on the small images, and the area above the given threshold was measured. Then, we used R software to count the number of images with aggregates.
Batch Macro processing code for aggregation estimates in ImageJ
run("Montage to Stack...", "images_per_row = 16 images_per_column = 16 border = 0");
macro "Measure Stack" {
saveSettings;
setOption("Stack position", true);
for (n = 1; n <= nSlices; n++) {
setSlice(n);
setAutoThreshold("Default dark stack");
//run("Threshold...");
run("Create Selection");;
run("Measure");
}
restoreSettings;
}
Cross-section images of Z-discs
Request a detailed protocolThoraces of Zasp52-GFP (Zasp52ZCL423), Zasp66-GFP (Zasp66ZCL0663), and Zasp67-GFP (Zasp67 fTRG) flies were fixed in paraformaldehyde and embedded into low-melting agarose blocks. Then, we sectioned the thoraces into 100-µm-thick sections using a standard vibratome. The sections were stained with TRITC-phalloidin and imaged using a LEICA SP8 confocal microscope and a 63x/1.4 oil objective. Images were then processed to isolate individual properly oriented Z-discs. To avoid artifacts from misoriented Z-discs, we selected only those in which the GFP signal coming from the Z-disc would cover the whole sarcomere area stained with phalloidin. To remove the background noise from cross-section Z-disc images, we took advantage of the geometrical properties of the Z-disc. We drew a selection circle around individual Z-discs and rotated the image by 2 degrees 180 times. Then, we calculated the average from all the rotated images.
We compared the Z-disc size categorical data using Fisher's exact test for count data, p-values above 0.001 were noted as not significant. We used Welch’s two-sample t-test to compare the relative fluorescence intensities between groups. We used Welch’s two-sample t-test to compare the aggregation frequency estimates.
Phylogenetic distribution of ZM domain proteins
Request a detailed protocolAll species that contain at least one protein with either the ZM (PF15936) or the LIM (PF00412) domain were downloaded from the PFAM web (Bateman et al., 2004; El-Gebali et al., 2019). A general phylogenetic tree was plotted based on the NCBI taxonomy data (Federhen, 2012) using the ape:plot.phyllo tool (Paradis, 2012).
RNAseq analysis
Request a detailed protocolWe uploaded all individual RNA-seq SRA reads from the project: PRJNA419412 (GEO: GSE107247) (Spletter et al., 2018). Each SRA read corresponds to a different time point in IFM development. We extracted the single-end reads using the fastq-dump tool. We then used Salmon v0.8.2 (Patro et al., 2017) to quantify the expression of all transcripts as Relative Transcripts Per Million (TPM) values. We used all the cDNAs from the D. melanogaster r6 reference transcriptome in fasta format (Adams et al., 2000). Then, the individual transcripts corresponding to the three Zasp genes were grouped according to their encoded protein domain architectures. Finally, we calculated the mean TPM values between similar isoforms and plotted the results over developmental time using R software.
The accession number for all reads used are: SRR1665023, SRR1665024, SRR1665025, SRR1665026, SRR1665027, SRR1665028, SRR6314253, SRR6314254, SRR6314255, SRR6314256, SRR6314257, SRR6314258, SRR6314259, SRR6314260, SRR6314261, SRR6314262, SRR6314263, SRR6314273, SRR6314274, SRR6314275, SRR6314276, SRR6314277, SRR6314278.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
-
NCBI Gene Expression OmnibusID GSE107247. Systematic transcriptomics reveals a biphasic mode of sarcomere morphogenesis in flight muscles regulated by Spalt.
References
-
The M-band: an elastic web that crosslinks thick filaments in the center of the sarcomereTrends in Cell Biology 15:477–485.https://doi.org/10.1016/j.tcb.2005.07.001
-
Aggregation of alpha-synuclein in lewy bodies of sporadic Parkinson's disease and dementia with Lewy bodiesThe American Journal of Pathology 152:879–884.
-
The pfam protein families databaseNucleic Acids Research 32:138D–141.https://doi.org/10.1093/nar/gkh121
-
Protein phase separation: a new phase in cell biologyTrends in Cell Biology 28:420–435.https://doi.org/10.1016/j.tcb.2018.02.004
-
Loss of enigma homolog protein results in dilated cardiomyopathyCirculation Research 107:348–356.https://doi.org/10.1161/CIRCRESAHA.110.218735
-
Selective deletion of long but not short cypher isoforms leads to late-onset dilated cardiomyopathyHuman Molecular Genetics 20:1751–1762.https://doi.org/10.1093/hmg/ddr050
-
Fluorescence-based methods in the study of protein-protein interactions in living cellsCurrent Opinion in Biotechnology 19:338–343.https://doi.org/10.1016/j.copbio.2008.06.003
-
PDLIM7 is a novel target of the ubiquitin ligase Nedd4-1 in skeletal muscleBiochemical Journal 473:267–276.https://doi.org/10.1042/BJ20150222
-
The pfam protein families database in 2019Nucleic Acids Research 47:D427–D432.https://doi.org/10.1093/nar/gky995
-
ZASP: a new Z-band alternatively spliced PDZ-motif proteinThe Journal of Cell Biology 146:465–475.https://doi.org/10.1083/jcb.146.2.465
-
The nebulin repeat protein lasp regulates I-band architecture and filament spacing in myofibrilsThe Journal of Cell Biology 206:559–572.https://doi.org/10.1083/jcb.201401094
-
InterPro in 2017-beyond protein family and domain annotationsNucleic Acids Research 45:D190–D199.https://doi.org/10.1093/nar/gkw1107
-
WAVE forms hetero- and homo-oligomeric complexes at integrin junctions in Drosophila visualized by bimolecular fluorescence complementationJournal of Biological Chemistry 285:40171–40179.https://doi.org/10.1074/jbc.M110.139337
-
A two-segment model for thin filament architecture in skeletal muscleNature Reviews Molecular Cell Biology 14:113–119.https://doi.org/10.1038/nrm3510
-
FlyBase at 25: looking to the futureNucleic Acids Research 45:D663–D671.https://doi.org/10.1093/nar/gkw1016
-
Fifty years of muscle and the sliding filament hypothesisEuropean Journal of Biochemistry 271:1403–1415.https://doi.org/10.1111/j.1432-1033.2004.04044.x
-
The LIM domain: from the cytoskeleton to the nucleusNature Reviews Molecular Cell Biology 5:920–931.https://doi.org/10.1038/nrm1499
-
The ZASP-like motif in actinin-associated LIM protein is required for interaction with the alpha-actinin rod and for targeting to the muscle Z-lineJournal of Biological Chemistry 279:26402–26410.https://doi.org/10.1074/jbc.M401871200
-
New aspects of myofibrillar myopathiesCurrent Opinion in Neurology 29:628–634.https://doi.org/10.1097/WCO.0000000000000357
-
Mechanical forces during muscle developmentMechanisms of Development 144:92–101.https://doi.org/10.1016/j.mod.2016.11.003
-
SMART: recent updates, new developments and status in 2015Nucleic Acids Research 43:D257–D260.https://doi.org/10.1093/nar/gku949
-
The vertebrate muscle Z-disc: sarcomere anchor for structure and signallingJournal of Muscle Research and Cell Motility 30:171–185.https://doi.org/10.1007/s10974-009-9189-6
-
Structure, dynamics, assembly, and evolution of protein complexesAnnual Review of Biochemistry 84:551–575.https://doi.org/10.1146/annurev-biochem-060614-034142
-
Cypher and enigma homolog protein are essential for cardiac development and embryonic survivalJournal of the American Heart Association 4:e001950.https://doi.org/10.1161/JAHA.115.001950
-
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
-
Ultrastructure of developing flight muscle in Drosophila. I. assembly of myofibrilsDevelopmental Biology 160:443–465.https://doi.org/10.1006/dbio.1993.1320
-
Structure of isolated Z-disks from honeybee flight muscleJournal of Muscle Research and Cell Motility 38:241–250.https://doi.org/10.1007/s10974-017-9477-5
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Myofibrillar myopathiesCurrent Opinion in Neurology 21:585–589.https://doi.org/10.1097/WCO.0b013e32830a752b
-
Mutations in ZASP define a novel form of muscular dystrophy in humansAnnals of Neurology 57:269–276.https://doi.org/10.1002/ana.20376
-
The initial steps of myofibril assembly: integrins pave the wayNature Reviews Molecular Cell Biology 10:293–298.https://doi.org/10.1038/nrm2634
-
Alpha-actinin is a component of the Z-filament, a structural backbone of skeletal muscle Z-disksThe Journal of Biochemistry 105:529–536.https://doi.org/10.1093/oxfordjournals.jbchem.a122701
-
Echocardiographic-determined septal morphology in Z-disc hypertrophic cardiomyopathyBiochemical and Biophysical Research Communications 351:896–902.https://doi.org/10.1016/j.bbrc.2006.10.119
-
Titin and nebulin in thick and thin filament length regulationSub-Cellular Biochemistry 82:285–318.https://doi.org/10.1007/978-3-319-49674-0_10
-
Mutations in cypher/ZASP in patients with dilated cardiomyopathy and left ventricular non-compactionJournal of the American College of Cardiology 42:2014–2027.https://doi.org/10.1016/j.jacc.2003.10.021
-
Sequence-encoded material properties dictate the structure and function of nuclear bodiesCurrent Opinion in Cell Biology 46:62–71.https://doi.org/10.1016/j.ceb.2017.03.003
-
ALP/Enigma PDZ-LIM domain proteins in the heartJournal of Molecular Cell Biology 2:96–102.https://doi.org/10.1093/jmcb/mjp038
-
Ablation of cypher, a PDZ-LIM domain Z-line protein, causes a severe form of congenital myopathyThe Journal of Cell Biology 155:605–612.https://doi.org/10.1083/jcb.200107092
Article and author information
Author details
Funding
Canadian Institutes of Health Research (MOP-142475)
- Nicanor González-Morales
- Yu Shu Xiao
- Matthew Aaron Schilling
- Océane Marescal
- Kuo An Liao
Canadian Institutes of Health Research (PJT-155995)
- Nicanor González-Morales
- Yu Shu Xiao
- Matthew Aaron Schilling
- Océane Marescal
- Kuo An Liao
Natural Sciences and Engineering Research Council of Canada (RGPIN 2016-06793)
- Nicanor González-Morales
- Yu Shu Xiao
- Matthew Aaron Schilling
- Océane Marescal
- Kuo An Liao
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We appreciate the help of Tuana Correia-Mesquita and Beili Hu in making transgenic flies. We thank the Bloomington Drosophila stock center, Vienna Drosophila Resource Center, and the Drosophila Genomics Resource Center for materials and the CIAN imaging facility for providing access to confocal microscopy. This work was supported by operating grants MOP-142475 and PJT-155995 from the Canadian Institutes of Health Research and by RGPIN 2016–06793 from the Natural Sciences and Engineering Research Council of Canada.
Copyright
© 2019, González-Morales 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
-
- 2,202
- views
-
- 280
- downloads
-
- 29
- 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
Aggregation of mutant forms of Huntingtin is the underlying feature of neurodegeneration observed in Huntington's disorder. In addition to neurons, cellular processes in non-neuronal cell types are also shown to be affected. Cells expressing neurodegeneration-associated mutant proteins show altered uptake of ligands, suggestive of impaired endocytosis, in a manner as yet unknown. Using live cell imaging, we show that clathrin-mediated endocytosis (CME) is affected in Drosophila hemocytes and mammalian cells containing Huntingtin aggregates. This is also accompanied by alterations in the organization of the actin cytoskeleton resulting in increased cellular stiffness. Further, we find that Huntingtin aggregates sequester actin and actin-modifying proteins. Overexpression of Hip1 or Arp3 (actin-interacting proteins) could restore CME and cellular stiffness in cells containing Huntingtin aggregates. Neurodegeneration driven by pathogenic Huntingtin was also rescued upon overexpression of either Hip1 or Arp3 in Drosophila. Examination of other pathogenic aggregates revealed that TDP-43 also displayed defective CME, altered actin organization and increased stiffness, similar to pathogenic Huntingtin. Together, our results point to an intimate connection between dysfunctional CME, actin misorganization and increased cellular stiffness caused by alteration in the local intracellular environment by pathogenic aggregates.
-
- Cell Biology
- Developmental Biology
The bone-resorbing activity of osteoclasts plays a critical role in the life-long remodeling of our bones that is perturbed in many bone loss diseases. Multinucleated osteoclasts are formed by the fusion of precursor cells, and larger cells – generated by an increased number of cell fusion events – have higher resorptive activity. We find that osteoclast fusion and bone resorption are promoted by reactive oxygen species (ROS) signaling and by an unconventional low molecular weight species of La protein, located at the osteoclast surface. Here, we develop the hypothesis that La’s unique regulatory role in osteoclast multinucleation and function is controlled by an ROS switch in La trafficking. Using antibodies that recognize reduced or oxidized species of La, we find that differentiating osteoclasts enrich an oxidized species of La at the cell surface, which is distinct from the reduced La species conventionally localized within cell nuclei. ROS signaling triggers the shift from reduced to oxidized La species, its dephosphorylation and delivery to the surface of osteoclasts, where La promotes multinucleation and resorptive activity. Moreover, intracellular ROS signaling in differentiating osteoclasts oxidizes critical cysteine residues in the C-terminal half of La, producing this unconventional La species that promotes osteoclast fusion. Our findings suggest that redox signaling induces changes in the location and function of La and may represent a promising target for novel skeletal therapies.