ABSTRACT
Myofibrils are long intracellular cables specific to muscles, composed mainly of actin and myosin filaments. The actin and myosin filaments are organized into repeated units called sarcomeres, which form the myofibrils. Muscle contraction is achieved by the simultaneous shortening of sarcomeres, which requires all sarcomeres to be the same size. Muscles have a variety of ways to ensure sarcomere homogeneity. We have previously shown that the controlled oligomerization of Zasp proteins sets the diameter of the myofibril. Here, we looked for Zasp-binding proteins at the Z-disc to identify additional proteins coordinating myofibril growth and assembly. We found that the E1 subunit of the oxoglutarate dehydrogenase complex localizes to both the Z-disc and the mitochondria, and is recruited to the Z-disc by Zasp52. The three subunits of the oxoglutarate dehydrogenase complex are required for myofibril formation. Using super-resolution microscopy, we revealed the overall organization of the complex at the Z-disc. Metabolomics identified an amino acid imbalance affecting protein synthesis as a possible cause of myofibril defects, which is supported by OGDH-dependent localization of ribosomes at the Z-disc.
INTRODUCTION
Striated muscles are long contractile cables that bridge two rigid structures, such as bones or regions of the exoskeleton. Muscle contraction is responsible for providing the energy required to displace those rigid units and thus animal movement. Muscles are formed by long intracellular cables, called myofibrils, that bridge the ends of muscles; the coordinated shortening of these myofibrils causes muscle contraction (Ahmed et al., 2022). Myofibrils are themselves composed of tandemly repeated units called sarcomeres. Sarcomeres are built up from a complex array of antiparallel actin and myosin filaments, where actin filaments are anchored to the flanks of the sarcomere at a protein complex called the Z-disc, whereas the myosin filaments are anchored at the center of the sarcomere at another protein complex called the M-line. Coordinated shortening of the sarcomeres produces myofibril contraction (Ahmed et al., 2022; Dos Remedios and Gilmour, 2017; Gunage et al., 2017; Lemke and Schnorrer, 2017; Luis and Schnorrer, 2021; Nikonova et al., 2020).
Because all sarcomeres contract in synchrony, their sizes are identical, and muscles use a variety of ways to ensure sarcomere size homogeneity. This is particularly true for the indirect flight muscle (IFM), a special muscle that evolved in insect lineages to sustain high-frequency contractions for prolonged periods. One of the adaptations of the IFM is very regular sarcomeres with a very small contractile range. In Drosophila, the IFM develops during the early pupal stages and then rapidly grows to fill most of the thoracic space during the late pupal stages (Reedy and Beall, 1993). Muscle growth is a very coordinated process. Myofibrils first form very thin longitudinal cables that stably grow by recruiting cytoplasmic proteins (Katzemich et al., 2013; Loison et al., 2018; Reedy and Beall, 1993; Spletter et al., 2018). The M-line and Z-disc grow together with the myofibril, actively mediating myofibril growth (González-Morales et al., 2019b; Katzemich et al., 2012, 2013; Orfanos et al., 2015).
Local and global mechanisms control muscle growth and sarcomere homogeneity. Local mechanisms acting on individual myofibrils or sarcomeres control sarcomere length, myofibril width and length of the I-band (the region that contains Z-discs and is devoid of myosin filaments). Sarcomere length is controlled by the length of the connecting protein titin (Tskhovrebova and Trinick, 2017) and the fine regulation of actin and myosin dynamics (Molnar et al., 2014; Shwartz et al., 2016). The length of the I-bands is controlled by the function of the Lasp proteins (Fernandes and Schöck, 2014). Myofibril width is set in place by controlled oligomerization of Zasp proteins at the Z-disc (González-Morales et al., 2019b). In contrast to these, global mechanisms act on the whole muscle, containing thousands of sarcomere units. One example of a global mechanism is the muscle growth coordination imposed by the continuous increase in tissue tension. As muscles grow, tension builds because of premature muscle contractions, and the increase in tendon and cuticle stiffness (Chu and Hayashi, 2021; Lemke and Schnorrer, 2017; Weitkunat et al., 2017, 2014). This continuous increase in tension coordinates the growth and the shape of myofibrils and mitochondria (Avellaneda et al., 2021). Other global mechanisms include growth regulation by the Hippo pathway (Kaya-Copur et al., 2021), the role of the E2F–DP heterodimeric transcription factor (Zappia and Frolov, 2016; Zappia et al., 2019) and the role of insulin (Demontis and Perrimon, 2009). The Hippo pathway and E2F–DP separately coordinate myofibril and mitochondrial growth rates by promoting the expression of myofibril and mitochondria proteins (Kaya-Copur et al., 2021; Zappia et al., 2019). Despite all these well-described mechanisms, the interconnections between the global growth cues and the local growth mechanisms have yet to be deciphered.
To gain insights into the mechanisms of Zasp proteins promoting myofibril growth, we screened for Z-disc proteins recruited by Zasp involved in myofibril diameter size regulation. Zasp proteins are members of the ALP/Enigma family; they consist of a PDZ, a ZM and up to 4 C-terminal LIM domains. They localize to the Z-disc, with their LIM domains at the very center of the disc and their PDZ domain at the periphery (Katzemich et al., 2013; Szikora et al., 2020). Through the PDZ domain, Zasp binds actinin and establishes the structural core of the Z-disc (Liao et al., 2016). Actinin anchors actin filaments from opposing sarcomeres. Zasp exists in two forms, a blocking and a growing form, with opposite roles during myofibril diameter growth (González-Morales et al., 2019a,b; Katzemich et al., 2013). The blocking forms prevent the recruitment of the growing isoforms to the Z-disc, whereas the growing isoforms recruit Zasp proteins to the Z-disc. The self-association of Zasp is mediated by a physical interaction between its LIM domains and ZM domain (González-Morales et al., 2019b).
We report that the E1 subunit of 2-Oxoglutarate Dehydrogenase (OGDH/E1), a crucial enzyme in the tricarboxylic acid (TCA) cycle, is recognized by the LIM domains of Zasp. OGDH/E1 is recruited to the growing Z-disc in addition to its mitochondrial localization, and its function is required for myofibril growth and assembly. The TCA cycle is a loop of chemical reactions and constitutes a metabolic buffering system. Anaplerotic metabolic pathways replenish the TCA cycle, whereas cataplerotic reactions use the TCA cycle metabolites (Martinez-Reyes and Chandel, 2020; Owen et al., 2002). An important step in the TCA cycle is the conversion of 2-oxoglutarate into succinyl-CoA, which is catalyzed by the OGDH complex – a giant enzyme cluster composed of multiples of three subunits (Tretter and Adam-Vizi, 2005). The dihydrolipoyllysine-residue succinyltransferase subunit (DLST, also known as E2 or CG5214) serves as a structural core unit (Skalidis et al., 2020). The OGDH (also known as E1 or Nc73EF) and dihydrolipoyl dehydrogenase (DLD, also known as E3 or CG7430) subunits sit in the periphery of the E2 core (Larkin et al., 2021; Skalidis et al., 2020; Tretter and Adam-Vizi, 2005) (Chen et al., 2015; Gruntenko et al., 1998; Yoon et al., 2017). The OGDH (also known as E1) subunit recognizes the 2-oxoglutarate substrate and provides specificity to the complex (Tretter and Adam-Vizi, 2005). This enzyme complex, therefore, fulfills crucial functions in mitochondria, providing metabolites and reduced electron carriers for oxidative phosphorylation. Here, we propose an additional role for the OGDH complex in providing amino acids for protein synthesis to sustain myofibril growth and Z-disc assembly. Together, these data suggest a novel link between a local Zasp-dependent myofibril growth control mechanism and a global myofibril growth control mechanism: amino acid availability.
RESULTS
Because Zasp proteins are often insoluble, we used a bioinformatic approach to look for proteins with similar evolutionary rates to the Zasp proteins in order to find proteins recruited by Zasp to the growing Z-disc. Evolutionary rate covariation (ERC) is a measurement of shared evolutionary history between proteins (Clark and Aquadro, 2010; Raza et al., 2019). We retrieved the ERC values for Zasp52, Zasp66 and Actn (α-actinin), and matched them to all Drosophila proteins with available ERC values (roughly 11,100 proteins) from a previously characterized project (Findlay et al., 2014). We found 16 proteins with ERC values greater than 0.5, with at least two of the three bait proteins. We then used an automatic clustering method to group the candidate proteins. Two clusters were obtained, the Zasp52/Zasp66 and the Zasp52/actinin groups (Fig. 1A), possibly reflecting the two roles of Zasp52: stabilization of actinin and recruitment of other Zasp proteins. As we primarily aimed for Z-disc proteins, we investigated which of the candidate proteins colocalizes with Zasp and actinin at the Z-disc. To do so, we obtained tagged versions of the candidate proteins and analyzed their intracellular localization. To avoid the signal coming from the mitochondria, we used a glycerol washing step to remove the mitochondria (Xiao et al., 2017). From the candidates, OGDH (also known as Nc73EF) had the most obvious Z-disc localization compared with the negative control (Fig. 1B). We then confirmed the presence of OGDH at the Z-disc by co-immunoprecipitation (IP). Because Zasp52 is highly insoluble in traditional immunoprecipitation (IP) buffers, we used the well-characterized Zasp-binding protein α-actinin. OGDH-GFP purified from adult thoraces precipitates with α-actinin, confirming the presence of OGDH at the Z-disc (Fig. S1). Oxoglutarate dehydrogenase (OGDH) is a crucial enzyme in the TCA cycle (Martinez-Reyes and Chandel, 2020).
Because OGDH localization at the Z-disc was unexpected, we decided to further study its role in myofibril assembly. First, we tested whether the GFP tag faithfully reflects OGDH localization by using a smaller FLAG tag that is less likely to cause localization artifacts. Using the UAS/Gal4 system, we expressed OGDH-FLAG and human OGDH-FLAG in the IFM, and analyzed their localization in glycerinated muscles using immunohistochemistry. Both forms localize mainly to the Z-disc, but also slightly to the M-line, suggesting that the Z-disc localization we observed in OGDH-GFP flies is not an artifact caused by the GFP (Fig. S1).
OGDH is the E1 subunit of the OGDH complex, the other two subunits are DLD and DLST. We therefore tested the localization at the Z-disc of all three subunits. We used GFP-tagged genomic versions of OGDH, DLD and DLST, and to determine the precise localization of the subunits within the Z-disc, we used a previously validated dSTORM (direct stochastic optical reconstruction microscopy) super-resolution microscopy method (Szikora et al., 2020). Overall, we found that all the subunits are present at the Z-disc. The small differences in distribution patterns provide clues about the structure of the complex at the Z-disc. (Fig. 1C-E). The distribution of DLST, the core subunit, and the OGDH molecules are largely overlapping (Fig. 1C,E,G). The third subunit, DLD, also accumulates at the Z-disc, but it localizes into two discrete bands alongside the Z-disc (Fig. 1D). The double-band pattern of DLD still exhibits a significant overlap with that of OGDH (Fig. 1F); therefore, these super-resolution data suggest that two OGDH complexes assemble at either side of the Z-disc center, with the DLD subunits concentrated on the outer periphery (Fig. S2).
We then analyzed the consequences of inactivating OGDH. We removed OGDH from the developing indirect flight muscles using Act88f-Gal4 in combination with an RNAi directed against all OGDH isoforms (HMS00554; we refer to these flies as OGDH-HM). We noted that all the OGDH-HM flies were flightless, and all their myofibrils were abnormal when compared with controls (Fig. 2A,B). OGDH-HM muscles had very small Z-discs and occasionally protein aggregates were visible in the cytoplasm (Fig. 2B, asterisk and arrow). These protein aggregates are Z-disc aggregates, as shown by the presence of Zasp66 and actinin, and the absence of the M-line marker obscurin (Fig. S3). Importantly, the OGDH defects are specific to the Z-disc, because the M-line marker obscurin crosses the entire myofibril (Fig. 2D). In addition, we used transmission electron microscopy to characterize the phenotype in more detail, and we noted that the Z-discs in OGDH-HM flies, which are located at the center of the fibrils, were severely reduced in size (Fig. 2F,G), suggesting a Z-disc growth or assembly defect. This results in unattached thin filaments that invade the H zone; however, the overall width of myofibrils is largely maintained, presumably because obscurin still tethers myosin filaments at the M-line (Fig. 2D).
To test the specificity of the RNAi knockdown, we used other approaches to reduce the function of OGDH in the muscles. We used two other RNAi lines that target different sequences of the OGDH gene (GD12778 and GD50393), and an indirect flight muscle-specific CRISPR-Cas9-based method targeting the OGDH gene (OGDH-CRISPR, TKO.GS00550). All these conditions had muscles with smaller Z-discs and protein aggregates (Fig. S3), confirming the OGDH requirement for proper myofibril and Z-disc growth. The Z-discs were smaller in the case of OGDH-HM and OGDH-CRISPR than in the two other RNAi conditions, whereas the aggregates were most common in OGDH-CRISP, followed by OGDH-HM, GD50393 and GD12778 (Fig. S3). Thus, overall, all four methods resulted in similar muscle defects, and the slight differences observed likely reflect the effectiveness of reducing OGDH levels in different conditions.
After the initial phenotypic characterization of OGDH, we explored the physical association between Zasp52 and OGDH proteins. First, we used a yeast two-hybrid assay to test for protein-protein binding. Yeast expressing Zasp52 and OGDH grow on selective media, suggesting a physical interaction, whereas yeast expressing Zasp66 or actinin together with OGDH did not grow (Fig. 3 and data not shown). We then used this assay to find the binding site mediating the interaction between OGDH and Zasp52. First, we paired OGDH with all the possible individual domains of Zasp52. Only Zasp52-LIM2a and -LIM2b together with OGDH restored yeast growth (Fig. 3A). We have previously shown that the LIM domains of Zasp52 bind to the ZM domains of Zasp proteins (González-Morales et al., 2019b). We noted a region in OGDH that weakly aligns with the ZM domain of Zasp66 and asked whether this sequence would also mediate OGDH-Zasp52 binding (Fig. 3B,C). We made a mutant OGDH version (referred to as OGDH-BM, short for binding mutant) that lacks this sequence, and we found that Zasp52 was unable to bind OGDH-BM (Fig. 3D). The LIM2a or LIM2b domains are also unable to bind OGDH-BM in Y2H assays (Fig. 3E).
Because both forms localize to the Z-disc when overexpressed, we used bimolecular fluorescence complementation assays to test the direct interaction between Zasp and OGDH at the Z-disc (Marescal et al., 2020). Muscles expressing OGDH fused to the N-terminal region of YFP, and Zasp52 fused to the complementary C-terminal region of YFP showed a clear fluorescent signal at the Z-disc (Fig. 4A and Fig. S4), suggesting OGDH and Zasp bind at the Z-disc. Importantly, the signal is lost in OGDH-BM conditions (Fig. 4B,C and Fig. S4). We then tested the role of Zasp52 in recruiting OGDH to the Z-disc. We imaged OGDH-GFP in two Zasp52 mutants. Zasp52MI02988 affects the PDZ domain but not the LIM domains, and Zasp52MI00979 forces a stop before the last three LIM domains but leaves the PDZ domain intact. OGDH localizes to the Z-disc in control and Zasp52MI02988 mutant muscles but not in Zasp52MI00979 mutants (Fig. 4D-F), further supporting the involvement of the LIM domains in OGDH binding.
Because overexpression studies are prone to subcellular localization artifacts, we also analyzed the endogenous OGDH protein. To this end, we created OGDH mutant alleles by incorporating a mCherry tag either into the wild-type OGDH or the OGDH-BM form (OGDH-WT-mCh and OGDH-BM-mCh, respectively). First, we tested the correct mitochondrial localization of both OGDH forms in the IFM by avoiding the glycerinated step during sample preparation. As expected, both OGDH-WT-mCh and OGDH-BM-mCh strongly localize to the mitochondria (Fig. S5). Homozygous OGDH-WT-mCh flies are viable and their myofibrils develop properly. In contrast, homozygous OGDH-BM-mCh is lethal. As heterozygotes, both alleles develop normal myofibrils. Because the mitochondrial signal is strong, it masks the Z-disc signal, so we used glycerinated muscles, to test the Z-disc localization of OGDH-BM-mCh. The Z-disc localization of OGDH-BM-mCh is diminished roughly by half compared with the OGDH-WT-mCh control (Fig. 4G), indicating that Zasp binding is at least partially required for Z-disc localization.
Because OGDH normally functions as one of the three subunits of the OGDH protein complex, we used tissue-specific CRISPR-Cas9 targeted mutagenesis to examine the function of the three OGDH complex subunits. We expressed Cas9 together with gRNAs targeting OGDH, DLST or DLD in the indirect flight muscles and observed a very similar sarcomere phenotype in all three cases, characterized by a strong impairment of myofibrillar organization and sarcomere structure (Fig. 5). Collectively, these data strongly argue that the OGDH complex is crucial for proper sarcomere arrangement, including Z-disc organization.
Given the crucial role of OGDH in the TCA cycle, we next asked whether the depletion of other TCA cycle components would affect sarcomere assembly. We expressed RNAi transgenes targeting the major components of the TCA cycle (Fig. 6A) and, strikingly, we observed sarcomere phenotypes in 83% of them (Fig. 6B-H, Table S1). Most notably, Aconitase, Isocitrate Dehydrogenase, OGDH complex E2 subunit and Succinyl CoA-Synthetase, which catalyze sequential steps in the cycle, had the most dramatic effects (Fig. 6B), often exhibiting a complete disintegration of the myofibril structure (Fig. 6E,F,H, red arrows) or myofibrils with more subtle defects (Fig. 6D,F-H, orange arrows). In the case of the few TCA components that do not affect sarcomere structure, enzyme redundancy is a likely explanation for the lack of phenotype. As the TCA cycle fuels the electron transport chain, which generates ATP, we also tested the role of ATP by analyzing muscles with a compromised electron transport chain. Removing Cox5a, an essential subunit of the Cytochrome c Oxidase (Mandal et al., 2005), with two different RNAi constructs, does not result in myofibril defects, suggesting that the electron transport chain does not account for the defects observed in the TCA cycle enzymes (Fig. 6B). We asked whether other TCA cycle enzymes localize to the Z-disc. We focused on Aconitase and Isocitrate Dehydrogenase because they are the two enzymes that directly precede OGDH and tagged versions are available. We used a GFP version of Aconitase (Acon1CC00758) and a Flag-tagged version of Isocitrate Dehydrogenase (UAS-Idh-FLAG) to test their localization. Acon1 has a faint Z-disc fluorescence signal, while Idh has a strong Z-disc signal, suggesting that both are present at the Z-disc (Fig. S6).
The TCA cycle is a metabolic hub that connects many aspects of cellular homeostasis including the biosynthesis of many amino acids (Martinez-Reyes and Chandel, 2020). To investigate how the loss of OGDH affects muscle metabolism, we performed a metabolomic analysis of OGDH-HM and control muscles (Fig. 7A). Consistent with a role in the TCA cycle, we observed a 2.7-fold increase in oxoglutarate accumulation, indicating the sensitivity of the approach. Interestingly, this change was paralleled with abnormal levels of numerous amino acids. The most extreme cases are Histidine and β-alanine, which are almost absent compared to control muscles (Fig. 7A and B). Homoserine and aspartic acid are also less abundant in OGDH-HM muscles (Fig. 7B), whereas valine, tyrosine, leucine, lysine, isoleucine, phenylalanine and sarcosine, a glycine intermediate, accumulate in OGDH-HM muscles (Fig. 7B). Overall, these data suggest a widespread effect on amino acid metabolism.
We then hypothesized that problems in protein synthesis could cause the myofibril and Z-disc defects in OGDH-depleted muscles. It is well known that when amino acids are missing, cells actively block protein synthesis by reducing ribosome biogenesis or by actively degrading them (Destefanis et al., 2020; Iadevaia et al., 2014). To test this possibility, we artificially blocked global protein synthesis by expressing the A subunit of the ricin toxin in muscles during different developmental periods (Moffat et al., 1992). Ricin expression during early development results in a complete absence of sarcomere structure (Fig. 8A, ∼28 h). In contrast, ricin expression slightly after the first appearance of sarcomeres and before sarcomere growth blocks sarcomere growth (Fig. 8A, ∼32 h and ∼56 h). Importantly, this results in normal myofibrils that are only reduced in width, not in the Z-disc-specific defects seen in OGDH mutants. Finally, ricin expression after the growth period of ∼80 h has little effect on sarcomere structure or size (Fig. 8A, ∼80 h). We then asked whether the number of ribosomes was affected by OGDH depletion. We used a GFP-tagged version of the ribosomal protein RpS5a (Kong et al., 2019). In control muscles, RpS5a localized to the perinuclear region and the Z-discs (Fig. 8B,D); however, in OGDH-depleted muscles, RpS5a fluorescence was strongly reduced, especially at the Z-disc (Fig. 8C,E,F). This suggests partial ribosome degradation or a termination in ribosome biogenesis.
DISCUSSION
Through a bioinformatic and Z-disc localization screen for myofibril growth regulators, we found the E1 subunit of the OGDH complex. We found that OGDH-E1 is a Z-disc and mitochondrial component required for the growth and assembly of myofibrils. Through a metabolomic analysis, we found that OGDH is required for the synthesis of some amino acids that we propose are required for myofibril growth.
Unexpected localization of the OGDH complex at the Z-disc
Although metabolic enzymes are not typically recognized as myofibril components, their presence at the Z-disc is not entirely unprecedented. Six glycolytic enzymes that catalyze consecutive reactions along the glycolytic pathway localize to the Z-disc (Sullivan et al., 2003; Wojtas et al., 1997). Our data, together with previous data on glycolytic enzymes, suggest that the Z-disc may be a common space for metabolic reactions that take place in the cytoplasm of muscles. Similarly, an interactome study using titin-BioID knock-in mice showed the presence of several enzymes at the Z-disc (Rudolph et al., 2020). Overall, our work demonstrates that the OGDH complex and at least two other TCA cycle enzymes localize to the Z-disc.
Using dSTORM super-resolution imaging, we show that the three subunits of the OGDH complex localize at the Z-disc and assemble into an asymmetric configuration in which the E3 subunit lies at the periphery of the complex. This is consistent with asymmetries observed by cryo-EM models of native α-keto acid dehydrogenase complexes from the Chaetomium thermophilum fungus (Kyrilis et al., 2021). Our work shows for the first time the asymmetry of the OGDH complex in animals.
OGDH is a novel Zasp-binding protein
Zasp proteins exist in two isoforms: the ones that contain LIM domains were designated as growing isoforms, while the ones that do not are called blocking isoforms (Katzemich et al., 2013; Liao et al., 2016). The growing isoforms recruit both growing and blocking forms through their LIM domains. The LIM domains bind the ZM domain present in both forms (González-Morales et al., 2019b). Here, we show that the LIM2 domain of Zasp52 recognizes the OGDH protein through a sequence very similar to the ZM domain. Interestingly, although all LIM domains can bind Zasp66, only the LIM2 domain binds OGDH, indicating that the other LIM domains might be used exclusively for Zasp proteins or for yet unidentified Z-disc proteins. Screening for proteins with ZM-like regions might be a strategy to find other Zasp-binding proteins. OGDH-BM-mCherry, which lacks the Zasp52 binding site, shows reduced localization at the Z-disc.
Control of myofibril growth
As sarcomere homogeneity is crucial for muscle function, several pathways coordinate Drosophila IFM growth. At the global muscle level, the Hippo pathway controls the rapid growth burst of the flight muscles by controlling the expression of many of the major sarcomere proteins (Kaya-Copur et al., 2021). The activating function of RBf with the E2F–DP heterodimeric transcription factor promotes the postmitotic expression of sarcomere proteins required for myofibril growth (Zappia and Frolov, 2016; Zappia et al., 2019), and the Insulin and mTOR signaling pathways regulate endoreplication, which is required for postmitotic growth (Demontis and Perrimon, 2009). At the local myofibril level, the Zasp oligomerization pathway then controls myofibril diameter (González-Morales et al., 2019b), and a combination of titin filaments and actin regulators set the sarcomere length (Molnar et al., 2014; Shwartz et al., 2016; Tskhovrebova and Trinick, 2017). Finally, some myofibril proteins are translated locally at the Z-disc (Denes et al., 2021). Here, we propose a link between the Zasp oligomerization system and the OGDH complex, which lies at the core of the TCA cycle. We propose a feedback system between Zasp oligomerization and the local amino acid pool that provides robustness to the myofibril diameter size control.
Mitochondrial versus Z-disc functions
The majority of OGDH and other TCA cycle enzymes localize to mitochondria where they are involved in two main functions: generating reduced electron carriers for eventual ATP production; and providing metabolic intermediates, mostly amino acids for protein synthesis. A smaller fraction surprisingly localizes to the Z-disc, where OGDH may locally provide amino acids for protein synthesis or may ‘moonlight’ in a different function, e.g. a structural role at the Z-disc.
Based on the currently available data, we propose that OGDH has dual roles in myofibril growth and Z-disc assembly: mitochondrial OGDH contributes to global amino acid production required for myofibril growth, whereas Z-disc OGDH is necessary to provide amino acids at the Z-disc for local protein synthesis that results in proper Z-disc growth and assembly. We can exclude a role for mitochondrial ATP production, because two different RNAi lines depleting cytochrome c oxidase do not affect myofibril growth or assembly (Fig. 6B). Metabolomics confirms, as expected, that amino acid levels are severely misregulated in OGDH mutants, which could therefore cause the observed phenotypes (Fig. 7). When OGDH function is impaired, many amino acids appear to be affected, some accumulate and others have highly reduced levels. The most affected amino acid is histidine. Histidine availability is linked to ribosome biogenesis. In fast-growing epithelial cells, high levels of histidine are required for ribosome biogenesis (Froldi et al., 2019). In addition, the loss of SLC15, a histidine amino acid transporter, correlates with low cytoplasmic histidine levels and with disruption of the mTOR pathway, which in turn controls ribosome biogenesis (Iadevaia et al., 2014; Kobayashi et al., 2014). Although we have not explored the requirements of specific amino acids, we speculate that low levels of histidine, like the ones in OGDH-HM muscles, are sufficient to reduce ribosome biogenesis. Many myofibrils have a smaller diameter in OGDH and other TCA cycle mutants, which is consistent with a global downregulation of histidine that results in reduced protein synthesis. However, the main OGDH phenotype cannot be explained by global downregulation of amino acids, because OGDH specifically affects Z-disc growth and assembly, but not M-line growth (Fig. 2). In contrast, global shutdown of protein synthesis disrupts only myofibril growth: myofibrils stop growing whenever protein synthesis is shut down, but sarcomeres look normal except for their smaller diameter (Fig. 8). We believe that a local production of amino acids at the Z-disc is more likely than a ‘moonlighting’ function of OGDH for three reasons: first, ribosomes are also enriched at the Z-disc, and their Z-disc localization depends on OGDH (Fig. 8); second, all three components of the OGDH complex and other TCA cycle components are found at the Z-disc, which would not be necessary if OGDH-E1 ‘moonlighted’ as a structural protein; and third, it appears unlikely that recruitment and precise localization to the Z-disc of the OGDH complex evolves without any associated function in amino acid production. To better dissect these contributions in the future will require OGDH mutants that cleanly disrupt enzymatic functions versus localization. Such mutants are currently unavailable, e.g. the OGDH-BM mutant identifies the area required for Z-disc localization (Figs 3 and 4), but may also affect enzymatic functions because it lies within the transketolase domain (Fig. 3B). An additional complication is that genetic mutants in these essential metabolic components are typically embryonic lethal and can therefore not be analyzed in adult muscles.
In conclusion, we observe a global alteration of amino acid homeostasis in OGDH mutants and specific localization of the OGDH complex to Z-discs, which aligns well with perturbed myofibril growth and specific Z-disc defects. This suggests important functions of the OGDH complex that go beyond its classical role in mitochondria.
MATERIALS AND METHODS
Experimental model
As a model organism, we used Drosophila melanogaster. Flies were raised at 25°C on standard cornmeal glucose media. A comprehensive list of all strains used and generated can be found in Table S2. Briefly, we used Saccharomyces cerevisiae for the yeast two-hybrid assays. The UAS/Gal4 system was used for transgene expression. The Act88F-Gal4 transgene was used to direct expression in the indirect flight muscles (Bryantsev et al., 2012). OGDH-TRiP.HMS00554 expresses RNAi against OGDH under UAS control (BDSC:33686). OGDHMI06026-GFSTF.1 is a GFP trap allele of OGDH (BDSC:59416). UAS-LacZ (BDSC:3356) was used to express LacZ. OGDH-GD.50393 expresses RNAi against OGDH under UAS control (VDRC:50393). Zasp52MI02988-mCherry, is a mCherry trap allele of Zasp52 (Xiao et al., 2017). Zasp52-GFP Zasp52ZCL423 is a GFP trap allele of Zasp52 (BDSC:58790). Zasp66-GFP Zasp66ZCL0663 is a GFP trap allele of Zasp66 (BDSC:6824). Obscurin-GFP is a fosmid duplication from the fTRG library of the obscurin locus with a GFP at the C terminus. Sls-GFP is a GFP trap allele of Sls (Orfanos et al., 2015). Zasp52MI02988 (BDSC_41034) and Zasp52MI00979 (BDSC:33099) are MIMIC-based alleles of Zasp52 that introduce early stop codons. They have been described previously (González-Morales et al., 2019b). DLD-GFP is a fosmid duplication from the fTRG library of the DLD locus with a GFP at the C terminus (VDRC:318906). DLST-GFP is a fosmid duplication from the fTRG library of the DLST locus. UAS-RA.CS2 expresses a cold-sensitive ricin toxin under UAS control (BDSC:538624). UAS-RpS5a-Venus was used as Ribosome-GFP and is a gift from Paul Lasko (McGill University, Quebec, Canada) (Kong et al., 2019). The following RNAi strains targeting different TCA cycle enzymes or cytochrome c oxidase were used: mAcon1-KK103809 (VDRC:103809), Idh3b-KK102960 (VDRC:102960), Idh3b-GD6219 (VDRC:14443), Idh3a-PGD5222 (VDRC:41191), Idh3a-GD5222 (VDRC:41192), Idh3a-GD16641 (VDRC:50828), Idh3a-KK107912 (VDRC:107912), CG5028-GD6271 (VDRC:52043), CG5028-KK102781 (VDRC:102781), Fum1-KK108008 (VDRC:108008), Fum1-HMC03334 (BDSC:51779), Cox5A-HMJ22367 (BDSC:58282), Cox5A-JF02700 (BDSC:27548), DLST-HMC03051 (BDSC:50650), DLST-KK109081 (VDRC:109081), DLD-KK102614 (VDRC:102614), kdn-KK107737 (VDRC:107737), mdh2-KK109040 (VDRC:109040), mdh1-KK108844 (VDRC:108844), skap-KK101171 (VDRC:101171), Scsα1-KK102542 (VDRC:102542), ScsβG-KK109063 (VDRC:109063), UAS-Idh-Flag (BDSC:56202) and Acon1CC00758 (BDSC:51542).
Evolutionary rate covariation
We obtained the ERC values from a previously characterized project (Findlay et al., 2014). We used Zasp66, Zasp52 and actinin as baits. This set contains pair-wise ERC values from 11,100 proteins calculated form multiple alignments of 12 Drosophila species. We then selected all the proteins with ERC values above 0.5 when compared with the bait proteins. We retrieved 94 proteins for Zasp66, 32 for Actinin and 25 for Zasp52. We used R to plot the values from the subset as a heat map and selected the proteins common to at least two of the bait proteins. Sixteen proteins were selected. From these, only OGDH and TER94 localized to the Z-disc.
Confocal microscopy imaging of flight muscles
The muscles were prepared for confocal imaging as described previously (Xiao et al., 2017). Briefly, the thoraces were dissected in half and incubated overnight at −20°C in relaxing-glycerol solution [20 mM sodium phosphate (pH 7.2), 2 mM MgCl2, 2 mM EGTA, 5 mM DTT, 0.5% Triton X-100 and 50% glycerol]. We then fixed the muscles in 4% paraformaldehyde and dissected them. For visualizing actin filaments, we used 488-phalloidin or 555-phalloidin (1:1000; Cytoskeleton) in PBS. Finally, we mounted the samples in Mowiol 4-88 mounting media (Sigma, 9002-89-5). All images were acquired using a 63×1.4 NA HC Plan Apochromat oil objective on a Leica SP8 confocal microscope. We used more than 10 flies for each experiment and randomly picked the muscle area to image. Control and experimental samples were prepared and imaged simultaneously, and imaged with comparable parameters. We used muscles from very young flies – 1-2 days old.
Bimolecular fluorescence complementation assay
BiFC assays were carried out as previously described (Marescal et al., 2020). The UAS-OGDH-NYFP construct was made using Gateway cloning using the OGDH-GEO09867 donor vector that contains the PA isoform as a donor and pBIDUAS-GV, pUAST-RfB-myc-NYFP as destination vector (Gohl et al., 2010). To make the UAS-OGDH-BM-NYFP construct, we first deleted the coding sequence for amino acids 741-769 in the OGDH-GEO09867 donor vector. The resulting vector was then transferred to pUAST-RfB-myc-NYFP using Gateway cloning. The resulting vectors were sequence verified and then inserted into the ZH-58A attp landing site. The UAS-Zasp52-PK-CYFP and the control lines have been described previously (Gohl et al., 2010; González-Morales, 2019b). At least 10 samples were used for each condition. We normalized the data to the basal noise levels and made plots in R software.
Tissue-specific CRISPR mutants
Tissue-specific CRISPR disruption works by expressing the Cas9 endonuclease in a specific tissue, using the UAS/Gal4 system together with a gene targeting gRNA expressed ubiquitously. The tissue containing the Cas9 protein generates small insertion or deletion mutations in the gene targeted by the gRNA (Port et al., 2014). To express the Cas9 protein in the IFM, we used Act88F-Gal4 with UAS-Cas9.P2 (BDSC:58986). As gRNA constructs, we used TKO.GS03432 targeting DLST (CG5214), TKO.GS00548 targeting DLD (CG7430), and TKO.GS00550 targeting OGDH (Nc73EF). The muscle defects were observed in 1- to 2-day-old flies.
Construction of DLST-GFP line
The fosmid carrying the GFP-tagged version of DLST (SourceBioscience: CBGtg9060A03104D) is part of the Flyfos library, a collection of fosmids that contain C-terminally GFP-tagged versions of genes at their genomic locations (Sarov et al., 2016). These constructs are then introduced into the fly genome by site-directed integration using the PhiC31 integrase (Bischof et al., 2007). We used P[CaryP]attP40 as the landing site for the DLST-GFP fosmid. Genome ProLab did the microinjections and Px3-RFP was used to screen for successful transformants.
Precise genome engineering at the OGDH locus
To precisely modify the OGDH locus, we used the recombination-mediated cassette exchange method using OGDHMI06026, which carries a MiMIC transposon between exons 5 and 6 (Venken et al., 2011). The rationale was to replace the MiMIC transposon with wild-type or mutant versions of the OGDH gene, starting with the sequence where OGDHMI06026 is inserted. In all cases, we added a C-terminal tag consisting of 6XHis and mCherry. First, we gene-synthesized the wild-type replacement construct and then mutagenized that construct using site-directed mutagenesis in bacteria. Gene synthesis and mutagenesis were carried out by Genscript. We created OGDHΔ741-769 (OGDH-BM) by deleting residues 741 to 769. Residue numbering is in accordance with the OGDH-PA isoform. All constructs were then introduced into OGDHMI06026 as a landing site. GenetiVision carried out the microinjections and initial confirmation of the mutants.
Yeast two-hybrid assays
Yeast two-hybrid assays were carried out as described previously (González-Morales et al., 2019b).
GC-MS sample preparation and metabolite measurements
The thoraces were dissected then flash frozen and crushed by mortar and pestle on liquid nitrogen (40 thoraces per sample). Frozen tissue powder was placed in pre-chilled Eppendorf brand tubes to which 1 ml of 80% methanol in water was added along with four 2.8 mm ceramic beads. Samples were subjected to 45 s of bead beating at 50 Hz (SpeedMill Plus homogenizer) four times. Samples were kept on ice between bead-beating sessions. Samples were then centrifuged at 1°C for 10 min at 21,130 g. Supernatants were transferred to fresh pre-chilled tubes containing 1 µl of 800 ng/µl 2H27-Myristic in pyridine. The protein concentration of the pellets was estimated and used for normalization. Samples were then dried by vacuum centrifugation operating at a sample temperature of −4°C (LabConco).
After drying, samples were subjected to a two-step derivatization: First, the samples were resuspended in 30 µl of 10 mg/ml methoxyamine:HCl in anhydrous pyridine (MOX). They were sonicated and vortexed for 15 s three times then centrifuged for 3 min at room temperature at 21,130 g. Incubation for methoximation was 30 min at room temperature. The samples were then centrifuged for 2 min at 21,130 g and the supernatants were transferred to GC-MS sample vials containing 250 µl glass inserts pre-filled with 70 µl of N-tert-butyldimethylsilyl-N-methyltrifluoroacetamide (MTBSTFA) and incubated at 70°C for 60 min.
An Agilent 5975C GC-MS equipped with a DB-5MS+DG (30 m×250 µm×0.25 µm) capillary column (Agilent J&W) was used for all GC-MS measurements, and data were collected by electron impact set at 70 eV both in scan (50-1000 m/z) and single ion monitoring modes. A volume of 1 ml of derivatized sample was injected in splitless mode with an inlet temperature set to 280°C, using helium as a carrier gas, and the flow rate was adjusted to 18 min for 2H27-myristic acid. The quadrupole was set at 150°C and the GC-MS interface at 285°C. The oven program for all metabolite analyses started at 60°C held for 1 min, then increased at a rate of 10°C/min until 320°C. Bake-out was at 320°C for 10 min. Sample data were acquired in scan mode (50-1000 m/z) or in single ion monitoring (SIM) with a 5 ms dwell time where the M-57 [M+•-C4H9•]+ fragment was used for quantitation (area under the curve) in both modes of data acquisition. Citrate and isocitrate used the m/z 459 ion for quantification as described previously (Mamer et al., 2013). The spectra and retention times of all metabolites reported were confirmed by methoxylamine–tert-butyldimethylsilylated authentic standards. For saturating metabolites, samples were diluted 1:25 with the same ratio of derivatization reagents and run in scan mode. Metabolite area under the curve was normalized to tissue weight.
Transmission electron microscopy
Muscles samples from 1- to 2-day-old flies were prepared for transmission electron microscopy imaging as described previously with slight modifications (González-Morales et al., 2017). Briefly, the thoraces were dissected in half and were treated with 5 mM MOPS (pH 6.8), 150 mM KCl, 5 mM EGTA, 5 mM ATP and 1% Triton X-100 for 2 h at 4°C. Samples were then washed in rigor solution [5 mM MOPS (pH 6.8), 40 mM KCl, 5 mM EGTA, 5 mM MgCl2 and 5 mM NaN3] and fixed in 3% glutaraldehyde, 0.2% tannic acid in 20 mM MOPS (pH 6.8), 5 mM EGTA, 5 mM MgCl2 and 5 mM NaN3 for 2 h at 4°C. Images were acquired on a Tecnai 12 BioTwin 120 kV transmission electron microscope with an AMT XR80C CCD camera (FEI).
Super-resolution dSTORM microscopy
Super-resolution imaging was carried out essentially as described previously (Szikora et al., 2020). Briefly, all the dSTORM images were captured under EPI illumination (Nikon CFI Apo 100×, NA=1.49) on a custom-made inverted microscope based on a Nikon Eclipse Ti-E frame. The laser (MPB Communication; 647 nm, Pmax=300 mW) intensity was controlled via an acousto-optic tunable filter (AOTF) set to 2-4 kW/cm2 on the sample plane. An additional laser (Nichia: 405 nm, Pmax=60 mW) was used for reactivation. Images were captured by an Andor iXon3 897 BV EMCCD digital camera (512×512 pixels with 16 μm pixel size). Frame stacks for dSTORM super-resolution imaging were captured at a reduced image size. A fluorescence filter set (Semrock, LF405/488/561/635-A-000) with an additional emission filter (AHF, 690/70 H Bandpass) was used to select and separate the excitation and emission lights in the microscope. During the measurements, the perfect focus system of the microscope was used to keep the sample in focus with a precision of <30 nm. Immediately before the measurement, the storage buffer of the sample was replaced with a GLOX switching buffer (van de Linde et al., 2011), and the sample was mounted onto a microscope slide. Typically, 20,000-50,000 frames were captured with an exposure time of 20 or 30 ms. The captured and stored image stacks were evaluated and analyzed with the rainSTORM localization software (Rees et al., 2013). Individual images of single molecules were fitted with a Gaussian point spread function and their center positions were associated with the position of the fluorescent molecule. Localizations were filtered via their intensity, precision and standard deviation values. Only localizations with precisions of <20 nm and standard deviation (σ) of 0.8≤σ≤1.0 were used to form the final image and for further analysis. Mechanical drift introduced by either the mechanical movement of the sample or thermal effects was analyzed and reduced using a correlation-based blind drift correction algorithm. Spatial coordinates of the localized events were stored, and the final super-resolved image was visualized with a pixel size of 10 nm. We used GFP-tagged lines to determine the nanoscopic localization of OGDH, DLD and DLST. Individual myofibrils were isolated from the IFM of anesthetized adult (∼24 h after eclosion) Drosophila as described previously (Burkart et al., 2007), with minor modifications. In brief, bisected hemithoraces were incubated in relaxing solution [100 mM NaCl, 20 mM NaPi (pH, 7.0), 5 mM MgCl2, 5 mM EGTA and 5 mM ATP] supplemented with 50% glycerol for 2 h at 4°C. Afterwards, the dorsal longitudinal muscles were isolated from the hemithoraces and dissociated by gently pipetting them in an Eppendorf tube in the presence of 0.5% Triton X-100. Dissociated myofibrils were centrifuged at 12,300 g for 2 min. Myofibrils were washed and centrifuged two more times in relaxing solution. Myofibrils were resuspended in a relaxing solution, and 20 µl of the sample was dropped on a glass coverslip and fixed with 4% paraformaldehyde (Alfa Aesar) in relaxing solution for 15 min. After washing three times in relaxing solution, the samples were blocked in blocking solution [5% goat serum (Sigma) and 0.1% Triton X-100 in relaxing solution] for 30 min in a humidity chamber. To detect OGDH, DLD and DLST-GFP, an anti-GFP antibody (1:1000; Abcam; ab13970) was applied overnight at 4°C in a blocking solution. After washing, goat anti-chicken secondary antibody coupled to AlexaFluor 647 (1:600; Invitrogen; A21449) was applied for 2 h at room temperature. F-actin was labeled with AlexaFluor 488-phalloidin (1:200; Thermo Fisher Scientific; A12379). The samples were thoroughly washed and stored in PBS before imaging. Experimental spatial resolution and localization precision were determined by the Fourier Ring Correlation and the Nearest Neighbor approaches (Endesfelder et al., 2014; Nieuwenhuizen et al., 2013). Drift-corrected measurement data generated by rainSTORM were first preconditioned and then evaluated by the FIRE and Coordinate Based Localization Precision Estimator codes. Spatial resolution of 52.8±6.3 nm and localization precision of 10.5±1.5 nm were achieved by the evaluation of ten randomly selected datasets. The values did not show any correlation with the samples, they were instead specified by the dSTORM microscope system and the data acquisition process. Based on these experimental results, localizations with theoretical (Thompson) localization precisions of less than 20 nm and a standard deviation (σ) of 0.8≤σ≤1.0 were used to form the final image and for further analysis.
Ribosome block
To block protein synthesis, we used UAS-RA.cs2, a cold-sensitive version of the ricin-A toxin subunit. Ricin-A inactivates ribosomes by the specific depurination of the 28S rRNA (Moffat et al., 1992). Ricin-A-TS is a temperature-sensitive allele that is active at 30°C but not at 20°C. We raised Act88F-Gal4 UAS-Ricin-A-TS flies at 20°C then transferred them to a 30°C incubator for 48 h, and then back into 20°C. The muscles were analyzed 2 days after emergence.
Acknowledgements
We appreciate the help of Beili Hu in making transgenic flies. We appreciate the support from the community resources such as the Bloomington Drosophila Stock Center, FlyBase and the Vienna Drosophila Resource Center. Confocal imaging was carried out at the Advanced Bioimaging Facility. Transmission electron microscopy was carried out with the assistance of Jeannie Mui at the Facility for Electron Microscopy Research. All GC-MS data were collected at the Metabolomics Innovation Resource.
Footnotes
Author contributions
Conceptualization: N.G.M., J.M., F.S.; Methodology: N.G.M., S.S., A.K., T.C.-M.; Software: N.G.M., S.S.; Validation: N.G.M., O.M., A.K., F.S.; Formal analysis: N.G.M., O.M., S.S., A.K., P.B., M.E., T.C.-M., J.M., F.S.; Investigation: N.G.M., A.K., P.B., F.S.; Resources: N.G.M., M.E., T.C.-M., F.S.; Data curation: N.G.M., F.S.; Writing - original draft: N.G.M., F.S.; Writing - review & editing: N.G.M., F.S.; Visualization: N.G.M., O.M., S.S., A.K., P.B., M.E.; Supervision: N.G.M., M.E., J.M., F.S.; Project administration: N.G.M., J.M., F.S.; Funding acquisition: N.G.M., J.M., F.S.
Funding
This work was supported by operating grants from the Canadian Institutes of Health Research (MOP-142475 and PJT-155995), by the Hungarian Science Foundation (OTKA) (K132782 to J.M.; FK138894 and PD128623 to S.S.), by the Hungarian National Research, Development and Innovation Office (NKFIH-871-3/2020 to J.M.), and by the János Bolyai Research Scholarship of the Magyar Tudományos Akadémia to S.S. Open Access funding provided by McGill University. Deposited in PMC for immediate release.
Data availability
All relevant data can be found within the article and its supplementary information.
Peer review history
The peer review history is available online at https://journals.biologists.com/jcs/lookup/doi/10.1242/jcs.260717.reviewer-comments.pdf
References
Competing interests
The authors declare no competing or financial interests.