Drosophila TRIM32 cooperates with glycolytic enzymes to promote cell growth
Abstract
Cell growth and/or proliferation may require the reprogramming of metabolic pathways, whereby a switch from oxidative to glycolytic metabolism diverts glycolytic intermediates towards anabolic pathways. Herein, we identify a novel role for TRIM32 in the maintenance of glycolytic flux mediated by biochemical interactions with the glycolytic enzymes Aldolase and Phosphoglycerate mutase. Loss of Drosophila TRIM32, encoded by thin (tn), shows reduced levels of glycolytic intermediates and amino acids. This altered metabolic profile correlates with a reduction in the size of glycolytic larval muscle and brain tissue. Consistent with a role for metabolic intermediates in glycolysis-driven biomass production, dietary amino acid supplementation in tn mutants improves muscle mass. Remarkably, TRIM32 is also required for ectopic growth - loss of TRIM32 in a wing disc-associated tumor model reduces glycolytic metabolism and restricts growth. Overall, our results reveal a novel role for TRIM32 for controlling glycolysis in the context of both normal development and tumor growth.
Introduction
The metabolism of all cells must adapt to meet the energetic and biosynthetic needs of growth and homeostasis (Lloyd, 2013; Zhu and Thompson, 2019). For example, tissues composed of non-dividing, differentiated cells must strike a balance between catabolic pathways that provide energy for cellular homeostasis and anabolic pathways that repair the cell and generate cell-type-specific molecules (Lloyd, 2013). In contrast, the metabolic requirements of cell growth and proliferation often require a shift toward anabolic pathways that favors the synthesis of macromolecules, such as proteins, lipids, nucleic acids, and complex carbohydrates (Zhu and Thompson, 2019). Striking this delicate balance between degradative and biosynthetic processes requires the integration of extracellular and intracellular information by complex signaling networks.
The mechanisms by which cell proliferation and tissue growth rewire metabolism to enhance biosynthesis are diverse and complex (Lunt and Vander Heiden, 2011). These changes in metabolic flux involve pathways such as the pyrimidine and purine biosynthesis, one carbon metabolism, and the interplay between the citric acid cycle and amino acids pools. However, the pathway most commonly associated with enhanced biosynthesis is glycolysis, where in biological systems ranging from yeast to human T-cells, glycolytic flux is often elevated in the context of cell growth and proliferation (Zhu and Thompson, 2019). This observation is particularly apparent in the fruit fly Drosophila melanogaster, where the onset of larval development is preceded by a metabolic switch that induces the coordinate upregulation of genes involved in glycolysis, the pentose phosphate pathway, and lactate dehydrogenase (LDH) (Tennessen et al., 2011). The resulting metabolic program allows larvae to use dietary carbohydrates for both energy production and biomass accumulation. Moreover, studies of Drosophila larval muscles reveal that this metabolic transition is essential for muscle growth and development, suggesting that glycolysis serves a key role in controlling growth (Tennessen et al., 2014b). The mechanisms that control glycolysis specifically in larval muscle, however, remain relatively unexplored. As a result, Drosophila larval development provides an excellent model for understanding how glycolysis and biomass production are regulated in a rapidly growing tissue. Moreover, since larval muscle increases in size without cell divisions, larval muscle provides an unusual opportunity to understand how glycolytic metabolism promotes growth independent of cell division.
Of the known factors that promote muscle development, TRIM32 is an intriguing candidate for coordinating metabolism with cell growth. This protein is a member of the Tripartite motif (TRIM)-containing family of proteins defined by an N-terminal RING domain, one or two B-boxes, a coiled-coil domain, and a variable C-terminal region (Tocchini and Ciosk, 2015; Watanabe and Hatakeyama, 2017). In TRIM32, six Ncl-1, HT2A, Lin-41 (NHL) repeats comprise the C-terminus and are proposed to mediate the diverse functions of TRIM32, including cell proliferation, neuronal differentiation, muscle physiology and regeneration, and tumorigenesis (Lazzari and Meroni, 2016; Tocchini and Ciosk, 2015; Watanabe and Hatakeyama, 2017). A single mutation in the B-box region of TRIM32 causes the multisystemic disorder Bardet-Biedl syndrome (BBS) (Chiang et al., 2006), while multiple mutations that cluster in the NHL domains result in the muscle disorders Limb-girdle muscular dystrophy type 2H (LGMD2H) and Sarcotubular Myopathy (STM) (Borg et al., 2009; Frosk et al., 2005; Lazzari et al., 2019; Nectoux et al., 2015; Neri et al., 2013; Schoser et al., 2005; Servián-Morilla et al., 2019).
A complete understanding of TRIM32 function is confounded by its ubiquitous expression and multitude of potential substrates for E3 ligase activity via the RING domain. Many known TRIM32 target substrates include proteins implicated in muscle physiology (Albor et al., 2006; Cohen et al., 2014; Cohen et al., 2012; Kudryashova et al., 2005; Locke et al., 2009; Volodin et al., 2017) or the prevention of satellite cell senescence (Kudryashova et al., 2012; Mokhonova et al., 2015; Servián-Morilla et al., 2019), consistent with a role for TRIM32 in LGMD2H. However, additional polyubiquitinated substrates, including p53, Abi2, Piasy, XIAP, and MYCN, are implicated in tumorigenesis (Albor et al., 2006; Izumi and Kaneko, 2014; Kano et al., 2008; Liu et al., 2014; Ryu et al., 2011). Importantly, TRIM32 protein levels are upregulated in multiple tumor types, suggesting that TRIM32 is a key player in growth regulation (Horn et al., 2004; Ito et al., 2017; Zhao et al., 2018). There is precedence for NHL function in controlling cell proliferation as two other Drosophila NHL-containing proteins, Brat and Mei-P26, act as tumor suppressors in the larval brain and female germline, respectively (Arama et al., 2000; Edwards et al., 2003).
Here, we provide a novel mechanism for TRIM32 in cell growth. Our data show that TRIM32 promotes glucose metabolism through the stabilization of glycolytic enzyme levels. This increased rate of TRIM32-mediated glycolytic flux generates precursors that are utilized for biomass production. Surprisingly, this mechanism operates in both non-dividing muscle cells as well as in proliferating larval brain cells, demonstrating a universal metabolic function for TRIM32 in growth control.
Results
TRIM32 binds to glycolytic enzymes
While NHL domain-containing proteins can interact with both RNAs and proteins (Tocchini and Ciosk, 2015; Watanabe and Hatakeyama, 2017), few bona fide TRIM32 binding partners have been identified. Causative mutations in human LGMD2H cluster in the NHL repeats (Figure 1A; Borg et al., 2009; Cossée et al., 2009; Frosk et al., 2005; Lazzari et al., 2019; Nectoux et al., 2015; Neri et al., 2013; Schoser et al., 2005), suggesting that this region may mediate protein-protein interactions important in disease prevention. We previously showed that mutations in Drosophila TRIM32 also show progressive larval muscle degeneration, despite modest sequence identity (~42%) across the NHL domains (LaBeau-DiMenna et al., 2012). To gain molecular-level insight into the functional conservation between the NHL regions in fly and mammalian TRIM32, we obtained a crystal structure of TRIM32_NHL at 2.6 Å resolution (Rwork/Rfree = 18.2/22.8%) (Figure 1—figure supplement 1A–C; Figure 1—source data 1). The structure was solved by molecular replacement, using the NHL-repeat region of the Drosophila Brat protein as a search model (PDB Code 1Q7F) (Edwards et al., 2003). This structure features a six-bladed β-propeller, where each NHL repeat is comprised of four antiparallel β-strands (Figure 1B; Figure 1—figure supplement 1C). We then used the TRIM32_NHL structure to construct a model for mouse TRIM32_NHL using SWISS-MODEL. Indeed, the superimposition of these two structures demonstrate that the NHL region of Drosophila TRIM32 is a faithful model to understand NHL function (Figure 1B).
To identify TRIM32_NHL-interacting proteins, we performed in vivo pulldowns followed by mass spectrometry (MS) to uncover associated peptides. Briefly, the C-terminal region of Drosophila TRIM32 containing all six NHL domains (AA 1061–1353) was fused in-frame with an N-terminal Glutathione S-transferase (GST) tag. This GST_TRIM32_NHL protein was expressed, purified, and incubated with third instar larval (L3) lysates. These resulting protein complexes were subjected to MS to detect peptide fragments that co-purify with TRIM32_NHL. In addition to the recovery of our bait protein, we detected peptides corresponding to Tropomyosin (Tm) and Troponin T (TnT) (Figure 1—source data 2), which are both known polyubiquitinated substrates of mammalian TRIM32 (Cohen et al., 2012). Surprisingly, glycolytic enzymes were also enriched in the list of possible NHL-binding proteins (Figure 1C; Figure 1—source data 2) and were further evaluated for physical interactions with TRIM32.
We utilized two independent methods to validate candidate protein interactions with TRIM32. First, an in vitro binding assay using purified proteins was developed with Tropomyosin 2 (Tm2) as a positive control. Indeed, incubation of TRIM32_NHL together with His-tagged Tm2 confirmed a physical interaction between these two proteins (Figure 2—figure supplement 1A). This assay was then used to test for physical interactions with candidate glycolytic enzymes. To this end, we expressed and purified His-tagged Aldolase (Ald) or His-tagged Phosphoglycerate mutase 78 (Pglym), both of which were low or absent in control pulldowns. Both Ald and Pglym were found to directly interact with the NHL domains of TRIM32 (Figure 2A,B). Note that no binding was detected between TRIM32_NHL and the His-tagged control protein SCIN, ruling out non-specific binding of the His tag. To extend these findings, we utilized immunoprecipitation experiments to test for in vivo interactions with full-length TRIM32. As expected, immunoprecipitation of TRIM32 successfully pulled-down Tm as an interacting protein at the predicted molecular weight of ~37 kD (Figure 2—figure supplement 1B). We also confirmed that both Ald and Pglym co-immunoprecipitate with anti-TRIM32 antibody. Interestingly, the bands corresponding to these proteins migrated at a higher molecular weight than the expected size present in input lysates (Figure 2C,D; asterisk), suggestive of a post-translational modification. Two approaches were taken to confirm the specificity of the slower migrating bands that correspond to Ald or Pglym. First, we concentrated L3 lysates before SDS-PAGE analysis. After overexposure of the western blots, faint bands corresponding to the higher molecular weight forms of Ald or Pglym were present (Figure 2C,D; asterisk). We also immunoprecipitated TRIM32 from larval lysates and probed for Ald or Pglym using antibodies raised against the orthologous human proteins (ALD or PGAM1). While these antibodies did not cross react well in input lysates, there was an obvious band corresponding to the higher molecular weight forms of both Ald and Pglym after pulling down TRIM32 (Figure 2—figure supplement 1D,E). Collectively, these data demonstrate that Ald and Pglym interact with the NHL domain of TRIM32 and in vivo, a subpopulation of these post-translationally modified glycolytic enzymes can be found in a complex with TRIM32.
The majority of TRIM32 studies have focused on the biological impact of polyUb chain addition onto substrates (Cohen et al., 2012; Kudryashova et al., 2005; Liu et al., 2014; Locke et al., 2009; Mokhonova et al., 2015). This modification targets a protein for proteasomal degradation and usually increases protein substrate levels in the absence of its E3 counterpart. Indeed, the known mammalian substrate Tm binds to TRIM32_NHL (Cohen et al., 2012) and accordingly, Drosophila Tm was modestly increased in larvae mutant for tn (Figure 2—figure supplement 1C). Strikingly, this same loss of TRIM32 resulted in a ~ 50% decrease in the protein levels of Ald or Pglym78 (Figure 2E,F). This downregulation was not due to a decrease in TRIM32-mediated transcription as mRNA levels of Ald or Pglym were not altered in tn-/- (Figure 2—figure supplement 1F,G). Taken together, these data highlight a unique role for TRIM32 in the stabilization of glycolytic enzyme levels.
Loss of TRIM32 disrupts glycolytic metabolism in larval muscles
To understand the biological implications for the physical interaction between TRIM32 and the glycolytic enzymes Ald and Pglym, we profiled 85 metabolites present in WT or tn-/- L3 larvae. Principal component analysis (PCA) revealed a distinct separation between WT and tn-/- experimental groups (Figure 3—figure supplement 1A). Examination of individual metabolites upon loss of TRIM32 showed a decrease in the terminal glycolytic products pyruvate and lactate as well as significant depletion of the glucose-derived metabolites glycerol-3-phosphate and 2-hydroxyglutarate (Figure 3A,B). These metabolic changes suggest that TRIM32 is required for sustaining the conversion of glucose to pyruvate and lactate in L3 larvae.
The production of pyruvate at the end of glycolysis has two fates, either a reduction to lactate or oxidation to CO2 in the mitochondrion (Mookerjee et al., 2015; TeSlaa and Teitell, 2014). To determine if reduced lactate levels in tn-/- results in increased respiration and therefore elevated ATP synthesis, we assessed these two outputs of overall energy metabolism in control or mutant whole larvae. Respirometry analysis revealed a modest increase in CO2 production upon loss of TRIM32 (Figure 3—figure supplement 1B). However, ATP levels were significantly decreased in tn-/- larvae (Figure 3—figure supplement 1C). Given that different larval tissues may have altered metabolic profiles that influence outputs of whole body metabolism, we sought to assess glycolytic activity in individual tissues.
A unique aspect of Drosophila larval development is the 200-fold increase in body size from first instar (L1) larvae through the L3 stage. In order to support high levels of biomass accumulation, glycolytic metabolism is increased prior to larval growth (Tennessen et al., 2011; Tennessen et al., 2014b). One larval tissue that requires high glycolytic activity is the somatic musculature, where WT body wall muscles undergo dramatic growth as development proceeds from second instar larvae (L2) (Figure 3C,G) to the L3 stage (Figure 3D,G). We, and others, previously reported a loss of structural integrity in TRIM32-deficient muscles (Domsch et al., 2013; LaBeau-DiMenna et al., 2012) and assumed that the associated reduction in muscle size was a secondary consequence of progressive tissue degeneration. While multiple independent alleles of tn have been shown to produce smaller muscles (Domsch et al., 2013; LaBeau-DiMenna et al., 2012), quantification of tn-/- confirmed a decrease in larval muscle diameter in both L2 (Figure 3E,G) and L3 individuals (Figure 3F,G). Thus, loss of TRIM32 compromises the glycolytic-driven growth of muscles during larval development.
As an initial assessment of glycolytic activity, we measured the extracellular acidification rate (ECAR) in WT or tn-/- muscle carcasses using the Agilent Seahorse XFe96 Analyzer. Since a major contributor to extracellular acidification is lactate excretion into medium during glycolysis (Mookerjee and Brand, 2015; TeSlaa and Teitell, 2014), we anticipated, and indeed, observed a reduction in ECAR upon loss of TRIM32 (Figure 3H). While quantitative, end-point assays such as ECAR, only measure a population of muscles at a single time point. For real-time measurements of glycolytic flux, the Xe96 Analyzer was used to calculate the proton efflux rate (PER). This assay specifically provides a measurement of protons derived from glycolysis by correcting for both the production of protons derived from mitochondrial reactions and for the buffering capacity of the assay medium. Briefly, the assay measures ECAR of the system under three conditions: (1) without manipulation, (2) following addition of the electron transport chain inhibitors rotenone and antimycin A, and (3) after addition of the glycolytic inhibitor 2-deoxy-D-glucose (2-DG). The end result of the assay is a measurement of acidification rate largely due to glycolytic sources. Our analysis revealed that glycolytic metabolism produces fewer protons in tn-/- muscle carcasses as compared with WT controls (Figure 3I), thus supporting our hypothesis that TRIM32 is required to support glycolytic flux.
To independently validate that loss of TRIM32 results in decreased glycolytic flux, we directly measured NADH levels within muscle tissues using fluorescence lifetime imaging (FLIM), which exploits the fluorescent characteristic of this cofactor to visualize NADH levels in cellular microenvironments (Martin et al., 2018; Provenzano et al., 2009; Szaszák et al., 2011; Yaseen et al., 2017). Since the fluorescent lifetime of NADH (the time required for NADH to decay when exposed to 740 nm wavelength light) is longer when NADH is bound to mitochondrial enzymes than in the free state, this assay can distinguish intracellular NADH pools (Bird et al., 2005; Skala et al., 2007). This analysis revealed that WT muscles analyzed by FLIM showed a markedly shorter lifetime than tn-/- muscles (Figure 3J). Since free, unbound NADH is predominate in highly glycolytic cells, here we confirm a reduction in the glycolytic profile upon loss of TRIM32. These data, taken together, strongly support the hypothesis that TRIM32 maintains glycolytic flux in larval muscle tissue.
One consequence of analyzing tn mutant alleles is the possibility that the observed muscle defects result from loss of TRIM32 in other tissues. To rule out TRIM32-mediated systemic defects, we utilized tissue-specific RNAi approaches (Brand and Perrimon, 1993). First, we reconfirmed that induction of three independent RNAi constructs targeting tn mRNA transcripts with the mef2-Gal4 muscle driver produced smaller muscles (Figure 4A,B,D; Brooks et al., 2016; Domsch et al., 2013; LaBeau-DiMenna et al., 2012). Knockdown of TRIM32 in a single muscle (5053> >tn RNAi, asterisk) within each hemisegment was sufficient to reduce muscle cell size (Figure 4C). As a proxy to monitor glycolytic activity, we assayed L-lactate levels in muscle carcasses. Consistent with our metabolomics data, the relative concentration of lactate levels was decreased in tn-/- muscles (Figure 4E). Induction of tn RNAi in muscle tissue (mef2 >tn RNAi) decreased lactate levels, while this same reduction did not occur upon TRIM32 knockdown in neurons (elav >tn RNAi) (Figure 4E). Here, we conclude that TRIM32-mediated regulation of glycolysis in muscle tissue is cell autonomous.
Estrogen-related receptor (ERR) is a nuclear hormone receptor that acts as a transcriptional switch in embryogenesis to induce genes required for aerobic glycolysis during larval growth (Tennessen et al., 2011). Therefore, we posited that genetic upregulation of carbohydrate metabolism genes via ERR may improve muscle growth and function. As a positive control, we expressed a cDNA encoding for TRIM32 in tn-/- muscle tissue (tn-/-; mef >TRIM32) and found that the muscle diameter was restored to WT (Figure 4F). This result also confirmed the cell autonomy of TRIM32. Overexpression of ERR in tn-/- muscles not only improved muscle size (Figure 4F), but also corrected the functional deficit associated with the inability of tn-/- to contract body wall muscles during pupal morphogenesis (Figure 4G,H; Domsch et al., 2013; LaBeau-DiMenna et al., 2012). Importantly, protein levels of Ald and Pglym were stabilized upon expression of ERR in tn-/- muscles (Figure 4I,J), indicating that restoration of glycolytic protein levels is sufficient to recover TRIM32-mediated growth defects.
TRIM32 maintains amino acid pools
Our metabolomics analysis revealed that loss of TRIM32 not only disrupts the production of glycolytic intermediates, but also induces a significant decrease in eleven of the twenty amino acids (Figure 5A). These changes in amino acid abundance are likely due to both decreased synthesis and increased catabolism (Figure 3—figure supplement 1D), as we observed that tn-/- exhibited a>1.5 decrease in not only serine and glycine levels which are normally synthesized from glucose, but also a reduction in the anaplerotic amino acids proline and aspartic acid. Moreover, loss of TRIM32 also induced a significant depletion of alanine, which is both synthesized from and catabolized into pyruvate. Overall, the metabolomic profile of tn-/- indicates that disruption of glucose catabolism results in depletion of larval amino acid pools, raising the possibility that decreased amino acid availability contributes to the tn-/- muscle defects. We tested this possibility by supplementing the diets of both mutant and control larvae with amino acid sources. As an initial approach, we first determined if the tn-/- larval phenotype exhibited enhanced sensitivity to nutrient deprivation. Indeed, when reared on starvation media, the muscle diameter of tn mutants was smaller than muscles of control larvae raised on the same media (Figure 5B,C,F). In contrast, mutant larvae fed a diet consisting of only yeast extract or supplemented with all 20 amino acids exhibited no decrease in muscle diameter (Figure 5D,E,F). We observed similar results in the context of L3 body mass and muscle contraction during the larval to pupal transition, whereby tn-/- raised on an agar-only diet exhibited dramatically more severe phenotypes compared with WT control larvae (Figure 5G,H,I). Supplementation of the larval diet with yeast extract or amino acids, however, suppressed both phenotypes in tn-/- (Figure 5G,H,I). These results suggest that tn-/- mutants are uniquely sensitive to dietary amino acids, consistent with the smaller amino acid pool size observed upon loss of TRIM32.
TRIM32 regulates glycolysis in a diversity of tissues
LDH expression generally correlates with LDH activity and increased glycolytic flux (Tanner et al., 2018; Ždralević et al., 2017). Only a few larval tissues show high LDH activity and exhibit elevated glycolytic rates (Eichenlaub et al., 2018; Li et al., 2017), including muscle and the larval brain (Figure 6—figure supplement 1A). Since we already established that tn-/- muscles were smaller, we wondered if growth of the larval brain also requires TRIM32. The overall size of larval brains dissected from WT (Figure 6A) or tn-/- (Figure 6B) L3 individuals showed a dramatic size difference, whereby loss of TRIM32 reduced the average area of the larval brain by ~40% (Figure 6G). To determine if a brain-specific reduction in TRIM32 is cell autonomous, we induced two independent tn RNAi lines in both neuronal (elav >tn RNAi) and muscle tissue (mef2 >tn RNAi). As expected, reduced larval brain size was observed upon RNAi knockdown of TRIM32 with elav-Gal4 (Figure 6C,D,G), but not with the mef2 driver (Figure 6E–G).
Increased biomass accumulation is a primary mechanism for the hypertrophic growth of post-mitotic larval muscles (Demontis and Perrimon, 2009), while larval brain development requires both cell growth and cell proliferation (Hartenstein et al., 2008). Thus, we wondered whether the TRIM32-mediated reduction in brain size was associated with altered glycolytic activity that reduced cell growth. Two experiments were performed to assess the glycolytic state in larval brain tissue. PER analysis (Neville et al., 2018) of individual larval brains isolated from tn-/- showed a reduced glycolytic rate compared to their WT counterparts (Figure 6H). Consistent with reduced substrate flux through the glycolytic pathway, L-lactate levels were also lower in dissected tn-/- larval brains (Figure 6I). Tissue-specific knockdown experiments verified that decreased lactate levels in larval brains resulted from inducing tn RNAi in neuronal, but not muscle tissue (Figure 6I). Next, we assessed how lower glycolytic activity affects overall larval brain size. Defective cell proliferation (assayed by EdU incorporation) or increased cell death (assessed using TUNEL labeling and cleaved-Caspase3 immunoreactivity) did not account for the overall brain size reduction (Figure 6—figure supplement 1B–I). Quantitation of individual brain cell size revealed a marked reduction upon loss of TRIM32 (Figure 6—figure supplement 1J–M) consistent with an increase in cell size as the primary mechanism for TRIM32-mediated tissue growth.
The elevated glycolytic rate that operates in larval muscle and brain tissue is analogous to the Warburg effect in rapidly proliferating cancer cells (Li et al., 2017; Tennessen et al., 2014b), whereby glucose metabolism is used to synthesize amino acids and other metabolites required for cell growth (Liberti and Locasale, 2016; Lunt and Vander Heiden, 2011). Thus, we hypothesized that TRIM32 could be a general regulator of highly glycolytic tumor cells. Unlike muscle or brain tissue, neither LDH activity (Wang et al., 2016) nor LDH-GFP expression (Figure 7A) were detectable in control wing discs, suggesting that this tissue does not exhibit elevated glycolytic activity. Accordingly, loss of TRIM32 did not alter the volume (Figure 7D) or area (Figure 7—figure supplement 1A–C) of mutant discs compared to WT control discs. Moreover, measurements of glycolytic flux using either PER assays (Figure 7—figure supplement 1D) or FILM analysis confirmed that TRIM32-deficient wing discs (Figure 7F,G) did not show altered glycolytic activity.
To examine the possibility that TRIM32 regulates the growth of highly glycolytic tumor cells, we confirmed that overexpression of the activated platelet PDGF/VEGF receptor (Pvract) in dpp-expressing wing disc cells increased LDH-GFP expression and promoted tissue overgrowth (Figure 7B,D; Figure 7—figure supplement 1C; Wang et al., 2016). Remarkably, removal of TRIM32 caused smaller tumors, effectively reducing the overall size of the wing disc (Figure 7C,D; Figure 7—figure supplement 1C). Analysis of LDH-GFP expression revealed biological variability as ~50% of TRIM32-deficient wing discs lost LDH-GFP expression (Figure 7E). FILM analysis validated a decrease in the glycolytic activity of tumors grown in tn-/- larvae compared to Pvr-induced tumors in WT individuals (Figure 7F,G). Collectively, our data show that TRIM32 is required for cell growth in both normal and tumorous glycolytic tissues, thus providing a novel molecular explanation for the upregulation of TRIM32 in multiple types of cancer cells.
Discussion
A unique feature of Drosophila larval development is the inherent glycolytic nature of muscle and brain tissue (Li et al., 2017; Tennessen et al., 2014b; Tixier et al., 2013), which promotes biomass synthesis during this stage of rapid organismal growth. Maintenance of such a high metabolic rate predicts that enzymes are present at sufficient concentrations in the cell to mediate the rapid shunting of intermediates through the pathway (Menard et al., 2014). We show here that TRIM32 directly interacts with and maintains the levels of two glycolytic enzymes. Decreased protein levels of both Ald and Pglym (and possibly other glycolytic enzymes) cripple this rapid flux, effectively blunting the generation of metabolic intermediates that contribute to anabolic synthesis necessary to sustain cell growth (Figure 8).
An alternative mechanism that limits cell and tissue growth is nutrient deprivation (Ahmad et al., 2018). In many organisms, the insulin/target of rapamycin (TOR) pathways integrate nutritional signals to physiologically control body size (Hyun, 2013). Implicit in this mechanism is the scaling of individual organs. Four pieces of evidence refute nutritional status as a mechanism for TRIM32-mediated tissue growth. First, the overall body size of tn mutant pupae is not smaller than their WT counterparts under starvation conditions, but instead elongated due to defective muscle contraction (Domsch et al., 2013; LaBeau-DiMenna et al., 2012; Figure 4G,H; Figure 5H,I). Second, the size of the wing disc and midgut, both tissues that show reduced growth in poorly fed larvae, is not altered upon loss of TRIM32 (Figure 7D; Figure 7—figure supplement 1A–C; Figure 7—figure supplement 2F–H). Surprisingly, cellular glucose uptake assayed by the non-metabolizable fluorescent glucose analog 2-[N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino]−2-deoxyglucose (2-NBDG), was normal in tn-/- larval brains and wing disc-derived tumors (Figure 7—figure supplement 2A–D), demonstrating that glucose is not a limiting substrate for glycolysis in isolated tn-/- tissues. Finally, even though systemic effects compromise feeding in tn-/- whole larvae (Figure 7—figure supplement 2E), food intake remained constant upon tissue-specific knockdown of TRIM32 in muscle or brains that show reduced tissue growth (Figure 7—figure supplement 2I,J). Collectively, these data demonstrate that TRIM32 functions in a cell autonomous manner to regulate tissue growth, independent of systemic nutritional status.
Li, et al., recently reported that pathways controlling lactate and glycerol-3-phosphate metabolism function redundantly in larval growth (Li et al., 2019). Removal of LDH and hence lactate production caused an increase in glycerol-3-phosphate, which was sufficient to maintain larval redox balance. Since both LDH and GPDH1 regulate redox balance necessary for maintaining high glycolytic flux to promote biomass accumulation, Ldh, Gpdh1 double mutants exhibit severe growth defects with reduced brain size. Our results show that loss of TRIM32 decreases both lactate and glycerol-3-phosphate levels (Figure 3A,B), thus mimicking the reduced carbohydrate metabolism in Ldh, Gpdh1 double mutants.
It is not clear how mutations in a ubiquitously expressed protein such as TRIM32 result in tissue-specific diseases. One prediction is that TRIM32 differentially interacts with proteins in diverse cell types to elicit distinct biological outputs. There is strong evidence to support this hypothesis in the context of LGMD2H. TRIM32 is upregulated in proliferating satellite cells and loss of this protein prevents myotube regeneration, partially through the misregulation of NDRG and c-Myc (Kudryashova et al., 2012; Mokhonova et al., 2015; Nicklas et al., 2012; Servián-Morilla et al., 2019). Muscle-specific targets that contribute to disease progression are less clear. TRIM32-mediated deregulation of key muscle substrates, including actin, α-actinin, tropomyosin, and desmin, contribute to muscle atrophy (Cohen et al., 2014; Cohen et al., 2012; Kudryashova et al., 2005; Locke et al., 2009), but studies have not been performed to directly test this model in the context of LGMD2H. Interestingly, mammalian glycolytic type II fibers are preferentially affected over oxidative type I fibers in muscle atrophy induced by aging/starvation, as well as in Duchenne’s and Becker muscular dystrophies (DMD) (Ciciliot et al., 2013; Pant et al., 2015). TRIM32 KO muscles also show a decrease in the glycolytic proteins GAPDH and PyK (Mokhonova et al., 2015), just as we observe a reduction in Ald and Pglym levels in tn-/- muscles, suggesting that TRIM32-mediated regulation of glycolysis may be a general mechanism that underlies some muscular dystrophies.
The multi-faceted roles exhibited by TRIM32 in muscle physiology and cancer seem quite different on the surface. However, control of glycolytic flux may be a common mode of regulation that has been overlooked. As in muscle tissue, the majority of studies on TRIM32 and cancer have focused on identifying substrates that are subject to poly-ubiquitination and subsequent proteasomal degradation. Piasy, p53 and Abi2 are known targets of TRIM32 E3 activity that regulate the proliferative balance in cancer cells (Albor et al., 2006; Kano et al., 2008; Liu et al., 2014). The proteolytic turnover of these proteins may affect signaling pathways independent of glycolytic TRIM32 regulation or may be a compensatory mechanism in response to metabolic shifts in which normal cells can transiently adopt cancer-like metabolism during periods of rapid proliferation.
Limitations of this study
How does this loss of TRIM32 lead to a reduction in glycolytic enzymes? Glycolytic proteins may be substrates for TRIM32 E3 ligase activity. It seems unlikely that TRIM32 polyubiquitinates Ald or Pglym for proteasomal degradation as protein levels are not elevated upon loss of this putative E3 activity. Furthermore, co-immunoprecipitaion of higher molecular weight forms of Ald or Pglym with TRIM32 suggests a yet unidentified post-translational modification. Another possibility, which we favor, is that the NHL domain of TRIM32 serves as a scaffold for the subcellular localization of glycolytic proteins to limit diffusion of substrates during glycolysis. This does not negate, but rather expands the repertoire of TRIM32 functions.
Materials and methods
Fly genetics
Request a detailed protocolFly stocks were reared on standard cornmeal media at 25°C unless otherwise indicated. The WT control strain was w1118 [Bloomington (BL) Drosophila Stock Center; BL3605). The tnΔA mutation (LaBeau-DiMenna et al., 2012) was maintained over the CyO, Tb/Sco (BL36335) balancer chromosome. Non-Tb individuals from the tnΔA/CyO, Tb1 stock were used for tn-/- analysis. The following Gal4 and/or RNAi lines were used: mef2-Gal4 (BL27390), elav-Gal4 (BL458), 5053-Gal4 (BL2702), UAS-tn RNAi-A [Vienna Drosophila Resource Center (VDRC); v19290), UAS-tn RNAi-B (BL31588), and UAS-tn RNAi-C (v19291). UAS-TRIM32 was previously described (LaBeau-DiMenna et al., 2012). The UAS-2XFLAG-dERR transgene was generated by amplifying the dERR cDNA using the oligos listed in Table 1. The resulting PCR product was sequenced, inserted into the NotI site of pUAST-attP, and injected to a strain containing the attP40 site by BestGene Inc (Chino Hills, CA). The transgenic dpp-Gal4, UAS-mCherry, LDH-GFP and UAS-Pvract flies were a kind gift from U. Banerjee (Wang et al., 2016). The LDH-GFP line (Figure 6—figure supplement 1A) was generated using a previously described pLdh genomic rescue construct (Li et al., 2017). Briefly, GFP was inserted at the 3’ end of the Ldh coding region using a PCR based method. The plasmid was injected into the strain attP40w by Rainbow Transgenics (Camarillo, CA) and the F1 generation was screened for transgene integration at the attP40 docking site. This transgene is capable of rescuing the Ldh mutant larval lethal phenotype (Chawla and Tennessen, in preparation).
Immunostaining and microscopy
Request a detailed protocolMuscle. Synchronized L2 or wandering L3 larvae were rinsed with 0.7% (w/v) NaCl/0.1% Triton, filleted in ice cold 1X PBS on Sylgard plates, and fixed in 4% (v/v) formaldehyde (Fisher) followed by three washes with 0.5% PBT. Phalloidin 488 or 594 was used to label F-actin (1:400, Molecular Probes). Other tissues. Wing discs, larval brains, or midguts were isolated from wandering L3 larvae and fixed for 25 min in 4% formaldehyde in PBS. Both discs and brains were stained overnight with DAPI (1:400) at 4°C. Midguts were stained with Phalloidin 488 (1:400, Molecular Probes). Anti-Caspase-3 (1:100, Cell Signaling Technology, Danvers, MA) staining on brains was performed overnight at 4°C followed by labeling with Alexa Fluor anti-rabbit 568 secondary antibody (1:400). Either muscle fillets, isolated brains or wing discs were mounted in anti-fade mounting medium (190% glycerol, 0.5% n-propyl gallate in 20 mM Tris buffer, pH = 8.0) and imaged using a Zeiss 700 confocal microscope.
Molecular biology
Request a detailed protocolThe NHL region of Drosophila TRIM32 (nucleotides 3231–4062) was PCR amplified using Phusion polymerase (ThermoFisher Scientific), digested with EcoRI and NotI, and ligated into the pGEX-5X-2 expression vector containing an N-terminal GST tag. Tm2, Ald, and Pglym were amplified by PCR and subcloned into the pT7HMT protein expression vector with either SalI/NotI or BamHI/NotI. Primers used for PCR amplification are included in Table 1.
Protein expression and purification
Request a detailed protocolProtein expression was performed in E. coli BL21 cells. A single colony was inoculated into 100 mL of LB media supplemented with the appropriate antibiotic at 37°C, incubated overnight and then diluted into 1L of Terrific Broth. Protein expression was induced at OD600 = 0.6 after the addition of Isopropyl β-D-1-thiogalactopyranoside (IPTG) to a final concentration of 1 mM (Geisbrecht et al., 2006). After 18 hr of incubation at 18°C, cells were centrifuged and were lysed using a microfluidizer for cell lysis. After centrifugation, the GST_TRIM32_NHL tagged protein was purified using GST affinity chromatography (GE Healthcare), further purified using size exclusion and ion-exchange Chromatography, and concentrated using 30K Amicon centrifuge columns (Millipore). The large scale expression and purification of Tm2, Ald, and Pglym was similar but used Ni2+ column for the initial purification step for His-tag binding. The purified proteins were stored at −80°C until use.
Crystallization, X-ray diffraction, structure solution, and refinement
Request a detailed protocolCrystals of a recombinant form of the NHL-repeat region of Drosophila TRIM32 were obtained by vapor diffusion of hanging drops. Briefly, a sample of purified protein was buffer exchanged into double-deionized water and concentrated to 5 mg/ml. Crystals were obtained from 2 µl droplets that had been established by mixing 1 µl of protein with 1 µl of precipitant solution [0.1M HEPES (pH 7.8), 0.2 M NaCl, and 25% (v/v) PEG-3350] and incubating over 500 µl of precipitant solution at 20°C. Rod-shaped crystals appeared within 2–3 days and grew to their full size within 1–2 weeks. Single crystals were harvested and briefly soaked in a cryopreservation buffer consisting of precipitation solution supplemented with an additional 20% (v/v) glycerol. Monochromatic X-ray diffraction data were collected at beamline 22-BM of the Advanced Photon Source at Argonne National Laboratory using incident radiation of λ = 1.000 Å. Reflections were indexed, integrated, and scaled using the HKL2000 software suite (Otwinowski and Minor, 1997). Initial phases were obtained by molecular replacement using program PHASER (McCoy et al., 2007) and a poly-alanine model derived from chain A of PDB entry 1Q7F. The model was constructed and refined through an iterative process consisting of automated building and refinement in PHENIX (Adams et al., 2002; Zwart et al., 2008), coupled with manual inspection and modifications. The final model consists of 292 protein residues, 85 ordered solvent molecules, and six ligand molecules. 94% of the protein residues occupy favored regions in the Ramachandran plot, with an additional 4% in allowed regions. Additional information on data collection and model statistics may be found in Figure 1—source data 1. The final model and structure factors have been deposited in the PDB under accession code 6D69.
Mass spectrometry (MS) analysis
Request a detailed protocolL3 larvae were homogenized in lysis buffer [50 mM Tris-HCl (pH 7.5), 100 mM NaCl, 10% (v/v) glycerol and 1 mM EDTA] plus inhibitors [1 mM Na3VO4, 5 mM NPPS, 2 mM PMSF, 2 ug/ml Leupeptin, 10 µM MG132, 1x Halt Pro inhibitor cocktail (Roche)]. 10 µg of purified TRIM32_NHL protein was coupled to Cyanogen bromide-activated-Sepharose 4B beads (Sigma-Aldrich), incubated for 2 hr at 4°C with 100 mg of lysates prepared from L3 larvae. Bead-protein complexes were washed 3X with wash buffer (50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 1 mM EDTA and 1% (v/v) Triton]. Control pulldowns were performed using beads alone. Beads containing protein complexes were sent to Oklahoma State University for MS analysis. Statistical analysis was performed using Perseus MaxQuant (Cox and Mann, 2012; Tyanova et al., 2016).
In vitro binding assay
Request a detailed protocolTo assess the interaction with TRIM32_NHL, 10 µg of purified candidate proteins (Tm2, Ald or Pglym) or the negative control (SCIN) were immobilized on Ni2+-NTA magnetic beads (ThermoFisher Scientific) for 1.5 hr followed by incubation with 10 µg of dTRIM32 NHL for 30 min at 4°C on a rotating platform. The complexes were washed 6x [300 mM NaCl, PBS (pH = 7.0) + 1% (v/v) Triton]. The bound proteins were eluted by boiling at 100°C in 6X Laemmli buffer for 10 min. The binding was analyzed on 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) under reducing conditions followed by Coomassie Blue Staining.
Co-Immunoprecipitation
Request a detailed protocolFor lysate preparation, WT third-instar larvae were homogenized in lysis buffer (50 mM Tris-HCl pH 7.5, 100 mM NaCl, 10% glycerol, 1% Triton, 1 mM EDTA) containing 1X Halt Protease Inhibitor (Thermo fisher Scientific) and 10 μM MG132. The lysate was precleared via centrifugation at 13,000 x for 20 min at 4◦C and the resulting protein concentration was measured using the Bradford Assay (Biorad, Hercules, CA). 500 μg of the lysate was set aside for input analysis. For immunoprecipitations, 20 μL of Protein A Sepharose 4b beads (Thermo Fisher Scientific) were washed twice with the lysis buffer. The beads were then incubated with 15 μg of Drosophila anti-TRIM32 antibody overnight at 4°C. The antibody conjugated beads were then mixed with 2 mg of lysate and incubated for 2 hr at 4◦C on a rotating wheel. Post incubation, the resin was centrifuged at 1000 x g for 3 min washed five times with wash buffer (1M NaCl, 1%Triton, PBS). 30 μL of 6X Laemmli buffer was used to elute bound proteins off the resin and these samples were denatured by heating at 100°C for 10 min. Eluted protein samples were subjected to 10% SDS PAGE gel electrophoresis and transferred to polyvinyl difluoride (PVDF) membranes (Pierce Biotechnology, Inc, Waltham, MA) for Western blotting with the following primary antibodies: rabbit anti-Drosophila Ald [1:1000, (Sullivan, 2003), rabbit anti- Drosophila Pglym [1:1000, (Sullivan, 2003), guinea pig anti-Drosophila TRIM32 [1:500, (LaBeau-DiMenna et al., 2012), rabbit anti-human ALD [1:1000, Biorad, Hercules, CA], rabbit anti-human PGAM1 [1:1000, Cell Signaling, Danvers, MA]. Post-treatment with the primary antibodies, blots were washed thoroughly with the wash buffer (TBS+ 0.1%Tween) and incubated with HRP conjugated and fluorescent secondary antibodies (dilution 1:5000) for 2 hr at room temperature. Protein detection was performed using the Prometheus ProSignal Pico western blotting detection kit (Genesee Scientific).
Metabolomics analysis
Request a detailed protocolPooled WT or tn-/- L3 larvae (25 each) were selected and washed with NaCl/0.1% Triton to remove food or debris. Each batch of pooled larvae were flash frozen in liquid nitrogen and sent to the University of Utah Metabolomics Core. Metabolites were extracted and derivatized before gas chromatography-mass spectrometry (GC-MS) analysis with an Agilent 5977B GC-MS. Data was collected using MassHunter software and metabolite identity was determined using MassHunter Quant. The metabolite data was normalized to standards, parsed, and Metaboanalyst 3.0 was used for statistical analysis and data processing. Independently, groups of 10 larvae were weighed on an analytical scale to determine the difference in body mass between WT and mutant cohorts before analysis. Six biological replicates were performed for each genotype.
In vivo feeding assays
Request a detailed protocolTo assess the impact of amino acid supplements on muscle mass in WT and tn-/-, embryos were reared from hatching until analysis at the L3 stage at 25°C on three different food conditions in 60 mm petri dishes: (1) Agar only control: 2.25% (w/v) agar in distilled water; (2) Yeast extract powder (mixture of amino acids, carbohydrates, peptides and soluble vitamins): 15% (w/v) yeast extract powder added to 2.25% (w/v) agar; or (3) Amino acids: 20 individual amino acids were individually weighed and added to 2.25% (w/v) agar (Lee and Micchelli, 2013). Both WT and tn-/- L3 larvae were dissected, stained with phalloidin, and analyzed using confocal microscopy for each condition. Note that agar contains two polysaccharides, agarose and agaropectin. Thus, the larvae are starved, but not completely devoid of nutrients.
ATP assay
Request a detailed protocol25 L3 larvae were pooled and homogenized in 100 µl of extraction buffer [6 M guanidine HCL, 100 mM Tris (pH 7.8), 4 mM EDTA], boiled for 5 min, and centrifuged at 13,000 x g for 5 min at 4°C. Protein concentrations were measured using the Bradford Assay kit (Biorad, Hercules, CA). ATP levels were determined using an ATP Determination Kit (Molecular Probes) as described in Tennessen et al. (2014a). 100 µl assays were performed in a 96 well plate and the luminescence was measured using Perkin Elmer EnSpire Multimode plate reader. Each sample was processed in triplicate and read in triplicate. The amount of ATP was normalized to total protein concentration.
Quantitation and statistical analysis
Request a detailed protocolMuscle diameter measurements. Z-stack images for each fillet were converted to a maximum intensity projection. The Polylines plugin from ImageJ was used to measure the muscle width of all VL3 muscles from dissected L2 and L3 individuals. Pupal case axial ratio determination. Pupa of the appropriate genotype were removed from vials, oriented dorsal side up, and attached to slides using a small drop of nail polish. Images were taken with a Leica M165 FC Stereomicroscope. Length and width measurements for each pupae were performed in ImageJ using the line and measure functions. Values were put into an Excel spreadsheet and the axial ratio (length/width) was calculated for each individual. The raw data was imported into Graphpad Prism 6.0 and graphed as a box and whiskers plot. N ≥ 10 for each genotype. Brain and midgut analysis. Z-stack images of isolated brain or midgut tissue areas were measured in ImageJ using the outline and measure functions N ≥ 8 for each genotype. Brain lobe cell size was determined using the analyze particle function after thresholding in ImageJ. N > 300 cells in 10 individual brain lobes for each genotype. Larval mass measurements. For each genotype and condition tested, pools of 10 larvae weighed on a digital scale. The average of each pool (N = 3) was plotted. Wing disc-associated tumor analysis. Z-stack images of each wing disc were used to measure the area and volume. For area measurements, the single plane that contained the maximum area for each disc was fully outlined using the free draw tool followed by the measure command in ImageJ. For volume quantitation, each disc in each Z-section was outlined, subjected to thresholding, and the stack command was used to measure the stacked volume. Statistics. Raw data was imported to GraphPad Prism 6.0 for statistical analysis and graph generation. All error bars represent mean ± standard deviation (SD). Statistical significances were determined using either student t-tests, Mann-Whitney tests or one-way ANOVA. Differences were considered significant if p<0.05 and are indicated in each figure legend.
Western blots
Request a detailed protocolFive whole L3 larvae were homogenized in 3X SDS sample buffer [188 mM Tris-HCl (pH 6.8), 3% (w/v) SDS, 30% (v/v) glycerol, 0.01% (w/v) bromophenol-blue, and 15% (v/v) β-mercaptoethanol], boiled at 95°C for 10 min, and centrifuged at 15,000 x g to remove cellular debris. To analyze overall protein levels, lysates were subjected to SDS–PAGE, transferred to polyvinyl difluoride (PVDF) membranes (Pierce Biotechnology, Inc, Waltham, MA), and probed with the appropriate primary antibodies: rabbit-Pglym [1:1000, (Sullivan, 2003)], rabbit-Ald, [1:1000, (Sullivan, 2003)], rat-Tm (1:500, Babraham Institute, Cambridge, UK) and mouse anti-ATPase 5α (1:10000, Abcam, Cambridge, MA). Horseradish Peroxidase (HRP) conjugated secondary antibodies (1:3000–1:5000, GE Healthcare, Chicago, IL) were used to detect primary antibodies. Protein detection was carried out using the ECL Plus western blotting detection kit (ThermoFisher Scientific, Waltham, MA). Densitometry analysis was performed by calculating the relative band intensities of candidate proteins to ATPase 5α loading control using ImageJ software.
Quantitative PCR (qPCR)
Request a detailed protocolRNA was isolated from a pool of five whole L3 larvae for each genotype using the RNAeasy Mini Kit (QIAGEN, Valencia, CA). After elution, RNA concentrations were determined and single strand complementary DNA (cDNA) was synthesized from 100 ng of RNA using the SuperScript VILO cDNA Synthesis Kit (Invitrogen, Carlsbad, CA). For qPCR, each cDNA sample was diluted to 1:50 and mixed with Power UP SYBR Green Master mix also mixed with the appropriate primers (Applied Biosystems, Foster City, CA). rp49 was used as the reference gene. Primers were synthesized by Integrated DNA Technologies (IDT):
rp49: F5’;-GCCCAAGGGTATCGACAACA-3’, R5’-GCGCTTGTTCGATCCGTAAC-3’
Ald: F5’- GGCCGCCGTCTACAAGGC-3’, R5’-GTTCTCCTTCTTGCCAGC-3’
Pglym78: F5-AGTCCGAGTGGAACCAGAAGA-3’, R5’-GGCTTGAAGTTCTCGTCCAG-3’
Three independent biological replicates were processed for each genotype and reactions were run in triplicate using the Quant Studio 3 Applied Biosystem with Quant studio design and analysis software. The average of the triplicates was used to calculate the 2-ΔΔCt values (normalized fold expression). Quantification of mRNA levels between different genotypes at the same time point was performed using the student t-test.
Agilent seahorse XFe96 analyzer experiments
Request a detailed protocolLarval brain, muscle, and wing disc preparation. The Agilent Seahorse XFe96 metabolic analyzer was placed in an incubator set to 12°C, and analyzer was set to 25°C with the heat on (running Wave software version 2.4). An Agilent Seahorse XFe96 cartridge (Agilent, Santa Clara, CA) was hydrated with 200 μl of calibrant solution (Agilent, Santa Clara, CA) overnight at 25°C. The next day, brains, muscles, or discs were dissected in phosphate buffered solution (PBS) and added to an Agilent 96-well cell plate (Agilent, Santa Clara, CA) containing 50 μl of Agilent Seahorse assay media with supplements required for specific assay (see glycolytic rate assay method). Tissue was sunk to the bottom of the well and centered in the middle between the three raised spheres. Forceps were used to lower the tissue restraint such that the plastic ring is facing toward the bottom of the well and the nylon screen is facing the top of the well. A probe was used to gently push the edge of the tissue restraint down toward the bottom of the well until the restraint did not move or float in the well. 130 μl of assay media was added to each well; resulting in a total of 180 μl final in each well. The cell plate was placed on the tray of the XFe96 analyzer. The instrument was for basal and glycolytic rate assays with all cycle procedures consisting of one-minute mixing, zero-minutes waiting, and 3-min measuring. Basal ECAR measurements using the XFe96. Basal levels of extracellular acidification (ECAR) were measured for a minimum of six cycles. Tissue restraints were measured without tissue as a control. Agilent Seahorse XF assay medium (Agilent, Santa Clara, CA) supplemented with 10 mM glucose, and 10 mM sodium pyruvate was used for all basal measurement assays. A minimum of 4 biological replicates were used to analyze the basal rate of WT and tn-/- larval muscle. Standard error of the mean (SEM) was used in analyzing metabolic measurement levels. Statistical significance was determined using the Holm-Sidak method with alpha = 0.05. Glycolytic Rate Assay. Analysis of glycolytic rate with mitochondrial-produced acidification subtracted was conducted using base medium without phenol red (Agilent, Santa Clara, CA), supplemented with 5 mM Hepes (Agilent, Santa Clara, CA), 2 mM glutamine, 10 mM glucose and 1 mM sodium pyruvate. 20 μl of 50 μM rotenone and antimycin-A was added to port A and injected at the 7th cycle, resulting in a final concentration of 5 μM rotenone and antimycin-A. 22 μl of 1M 2-deoxyglutarate (2-DG) was added to port B and injected at the 12th cycle, resulting in a final concentration of 100 mM 2-DG. The software package included with this kit analyzes the oxygen consumption and extracellular acidification rates, while factoring in the buffer capacity of the media. It also calculates the acidification caused by the mitochondria and subtracts this from the data. This method produces the proton efflux rate (PER). A minimum of four biological replicates were used to analyze the glycolytic rate of WT and tn-/- larval brains. Standard error of the mean (SEM) was used in analyzing metabolic measurement levels. Statistical significance was determined using the Holm-Sidak method with alpha = 0.05. Normalization of XFe96 measurements by protein concentration. XFe96 analysis data was normalized by protein concentration. Protein normalization was conducted using the Pierce 660 nm protein assay reagent (Thermo Fisher Scientific, Waltham, MA). Each sample was homogenized in default lysis buffer (50 mM Tris (pH 7.5), 125 mM NaCl, 5% glycerol, 0.2% IGEPAL, 1.5 mM MgCl2, 1 mM DTT, 25 mM NaF, 1 mM Na3VO4, 1 mM EDTA and 2 × Complete protease inhibitor (Roche, Indianapolis, IN) on ice, incubated on ice for 15 min, and centrifuged at full speed for 15 min at 4°C. Supernatant was collected and measured as indicated in the Pierce 660 nm protein assay reagent manual. Samples were assayed using a Biotek Cytation three plate reader. PER and muscle ECAR values were divided by μg amount of protein measured to determine the normalized pmol/min/μg rate of proton efflux and mpH/min/μg rate of extracellular acidification. Standard error of the mean (SEM) was used in analyzing metabolic measurement levels. Statistical significance was determined using the Holm-Sidak method with alpha = 0.05.
Respirometry
Request a detailed protocolLarval metabolic rates were assessed by indirect calorimetry, measuring CO2 production with flow-through respirometry using the Multiple Animal Versatile Energetics platform for metabolic phenotyping (MAVEn, Sable Systems International, Las Vegas, NV). Baseline ultra zero air was provided from a compressed air cylinder and regulated at approximately 20 mL/min through each of 16 respirometry chambers simultaneously. Larvae were measured in groups of five in each chamber. The MAVEn automated flow switching between chambers and baseline (interleave ratio 16:1; dwell time per chamber 2 min). Differential carbon dioxide concentration was measured with a LiCor 7000 infrared gas analyzer. The wet mass of insects prior to and following respirometry was measured to the nearest 0.01 mg (Mettler Toledo XS22SDU Analytical Balance). N = 6 chambers with five larvae in each chamber for each genotype.
FLIM
Request a detailed protocolThe multiphoton based (Denk et al., 1990) lifetime and intensity imaging was performed on a custom multiphoton laser scanning system built around an inverted Nikon Eclipse TE2000U at the Laboratory for Optical and Computational Instrumentation (Yan et al., 2006). A 20x air immersion objective (Nikon Plan Apo VC, 0.75 NA) (Melville, NY, USA) was used for all imaging. For NADH imaging, data was collected using an excitation wavelength of 740 nm, and the emission was filtered at 457 ± 50 nm (Semrock, Rochester, NY) for the spectral peak for NADH/NADPH. For intensity imaging, the excitation was set at 980 nm, and an emission 520 ± 35 filter was used (Semrock, Rochester, NY). The FLIM fitting process was done according to the methods sections describing FLIM analysis performed for the same scope (Ghanbari et al., 2019). For each sample, around eight neighboring fields were randomly selected, and the average value of lifetime and free NADH ratio were calculated.
L-lactate assays
Request a detailed protocolL-lactate levels in the brain and muscle were measured in the indicated genotypes. L3 larval brains and muscle carcasses were dissected and placed in ice-cold 1X PBS. Brains (n = 20 each, total 60 brains per genotype) or muscle tissue (N = 8, total 24 carcasses per genotype) were pooled together and homogenized in 50 μL 1X PBS. Bradford Assay was used to quantitate protein concentrations. Each lysate was transferred to a 96-well plate and incubated with lactate dehydrogenase and NAD/MTT for 2 hr at room temperature (EnzyChrom Glycolysis Assay Kit, BioAssay Systems). The intensity of the reduced dye was measured at 565 nm, which is directly proportional to the concentration of the L-Lactate in the sample.
EdU incorporation assay
Request a detailed protocolDissected larval brains were incubated in Drosophila Schneider's medium containing 10 uM EdU (ThermoFisher Scientific, Waltham, MA) for 2 hr at room temperature. Tissues were fixed and Click-iT EdU staining was performed according to the manufacturer’s protocol. For quantification, images of the whole larval brain were captured using a Zeiss 700 confocal microscope. ImageJ was used to manually count the number of EdU positive cells.
Glucose uptake assay
Request a detailed protocolWandering L3 larvae of the following genotypes (dpp-Gal4 UAS-mCherry, UAS-LDHGFP >UAS Pvract, tn-/-; dpp-Gal4 UAS-mCherry, UAS-LDHGFP >UAS Pvract) were placed in Drosophila HL3 buffer (NaCl 70 mM, KCl 5 mM, CaCl21.5 mM, MgCl220 mM, NaHCO310 mM, trehalose 5 mM, sucrose 115 mM, and Hepes 5 mM (pH 7.2) on ice. Wing discs were dissected in ice cold PBS followed by incubation with 2-NBDG (2.5 mg/mL, Cayman Chemical, Ann Arbor, MI) for 15 min at room temperature. The tissues were washed with 1X PBS twice, fixed and imaged. Similarly, for analyzing glucose uptake in larval brain, WT and tn-/- mutant brains were dissected and incubated with 2-NBDG as indicated above.
Food intake assay
Request a detailed protocolApproximately 15–35 early L3 larvae of the indicated genotypes were fed on yeast paste with 0.16% Erioglaucine dye for 24 hr (Aditi et al., 2016). Post feeding larvae were pooled together and washed thrice with distilled water. The pooled larvae were homogenized in 250 μL of ddH2O, centrifuged at 15,000 rpm for 20 min. 225 μL of supernatant was carefully transferred to a 1.5 mL tube containing 50 μL 100% ethanol and centrifuged for 10 min. 250 μL of the supernatant was placed in a new tube and centrifuged again for 5 min. Following centrifugation, 200 μL of the supernatant was placed in a 96-well crystal plate and the OD was measured at 633 nm. Three groups of 15 larvae of each genotype were analyzed.
TUNEL assay
Request a detailed protocolFor TUNEL labeling, fixed brains were treated with 20 µg/mL Proteinase K in PBS for 20 min. For positive control, WT brains were incubated with 100 µl of DNase (1 U/ µl) for 10 min at room temperature post Proteinase K treatment. The tissues were rinsed three times with PBT, followed by a wash in 100 µl terminal deoxynucleotidyltransferase (TdT) equilibration buffer (from the Kit, DeadEnd Fluorometric TUNEL System Promega). For the enzymatic reaction, a working solution of TdT enzyme (30%) with the reaction buffer (70%) was made and added to the tissues. The samples were then incubated for 2 hr in a 37°C incubator and flickered every 30 min to mix the contents. The reaction was terminated by immersing the tissues in stop buffer (0.3 M NaCl/0.03 M sodium citrate) for 15 min at room temperature. DAPI (1:400) was used as a counterstain.
Data availability
Diffraction data have been deposited in PDB under the accession code 6D69. All other data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for all raw data.
-
RCSB Protein Data BankID 6D69. Crystal Structure of the NHL Repeat Region of D. melanogaster Thin.
References
-
PHENIX: building new software for automated crystallographic structure determinationActa Crystallographica Section D Biological Crystallography 58:1948–1954.https://doi.org/10.1107/S0907444902016657
-
Altered lipid metabolism in Drosophila model of Huntington's diseaseScientific Reports 6:31411.https://doi.org/10.1038/srep31411
-
Evolutionary responses ofFrontiers in Ecology and Evolution 6:00095.https://doi.org/10.3389/fevo.2020.00095
-
Targeted gene expression as a means of altering cell fates and generating dominant phenotypesDevelopment 118:401–415.
-
Optimization of wrMTrck to monitorDrosophilalarval locomotor activityJournal of Insect Physiology 94:11–17.https://doi.org/10.1016/j.jinsphys.2016.07.007
-
Muscle type and fiber type specificity in muscle wastingThe International Journal of Biochemistry & Cell Biology 45:2191–2199.https://doi.org/10.1016/j.biocel.2013.05.016
-
Ubiquitylation by Trim32 causes coupled loss of desmin, Z-bands, and thin filaments in muscle atrophyThe Journal of Cell Biology 198:575–589.https://doi.org/10.1083/jcb.201110067
-
Trim32 reduces PI3K-Akt-FoxO signaling in muscle atrophy by promoting plakoglobin-PI3K dissociationThe Journal of Cell Biology 204:747–758.https://doi.org/10.1083/jcb.201304167
-
Abba is an essential TRIM/RBCC protein to maintain the integrity of sarcomeric cytoarchitectureJournal of Cell Science 126:3314–3323.https://doi.org/10.1242/jcs.122366
-
Model of the brain tumor-Pumilio translation repressor complexGenes & Development 17:2508–2513.https://doi.org/10.1101/gad.1119403
-
Hutterite brothers both affected with two forms of limb girdle muscular dystrophy: lgmd2h and LGMD2IEuropean Journal of Human Genetics 13:978–982.https://doi.org/10.1038/sj.ejhg.5201436
-
An optimized system for expression and purification of secreted bacterial proteinsProtein Expression and Purification 46:23–32.https://doi.org/10.1016/j.pep.2005.09.003
-
Cortex-wide neural interfacing via transparent polymer skullsNature Communications 10:1500.https://doi.org/10.1038/s41467-019-09488-0
-
The development of the Drosophila larval brainAdvances in Experimental Medicine and Biology 628:1–31.https://doi.org/10.1007/978-0-387-78261-4_1
-
Body size regulation and insulin-like growth factor signalingCellular and Molecular Life Sciences 70:2351–2365.https://doi.org/10.1007/s00018-013-1313-5
-
Satellite cell senescence underlies myopathy in a mouse model of limb-girdle muscular dystrophy 2HJournal of Clinical Investigation 122:1764–1776.https://doi.org/10.1172/JCI59581
-
TRIM32 ubiquitin E3 ligase, one enzyme for several pathologies: from muscular dystrophy to tumoursThe International Journal of Biochemistry & Cell Biology 79:469–477.https://doi.org/10.1016/j.biocel.2016.07.023
-
The warburg effect: how does it benefit Cancer cells?Trends in Biochemical Sciences 41:211–218.https://doi.org/10.1016/j.tibs.2015.12.001
-
E3 ubiquitin ligase TRIM32 negatively regulates tumor suppressor p53 to promote tumorigenesisCell Death & Differentiation 21:1792–1804.https://doi.org/10.1038/cdd.2014.121
-
TRIM32 is an E3 ubiquitin ligase for dysbindinHuman Molecular Genetics 18:2344–2358.https://doi.org/10.1093/hmg/ddp167
-
Aerobic glycolysis: meeting the metabolic requirements of cell proliferationAnnual Review of Cell and Developmental Biology 27:441–464.https://doi.org/10.1146/annurev-cellbio-092910-154237
-
GSK3β regulates brain energy metabolismCell Reports 23:1922–1931.https://doi.org/10.1016/j.celrep.2018.04.045
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
The E3 ubiquitin ligase TRIM32 regulates myoblast proliferation by controlling turnover of NDRG2Human Molecular Genetics 24:2873–2883.https://doi.org/10.1093/hmg/ddv049
-
The contributions of respiration and glycolysis to extracellular acid productionBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1847:171–181.https://doi.org/10.1016/j.bbabio.2014.10.005
-
Measurement and analysis of extracellular acid production to determine glycolytic rateJournal of Visualized Experiments e53464.https://doi.org/10.3791/53464
-
Detection of TRIM32 deletions in LGMD patients analyzed by a combined strategy of CGH array and massively parallel sequencingEuropean Journal of Human Genetics 23:929–934.https://doi.org/10.1038/ejhg.2014.223
-
A novel ex vivo method for measuring whole brain metabolism in model systemsJournal of Neuroscience Methods 296:32–43.https://doi.org/10.1016/j.jneumeth.2017.12.020
-
Processing of X-ray diffraction data collected in oscillation modeMethods in Enzymology 276:307–326.https://doi.org/10.1016/S0076-6879(97)76066-X
-
Multiphoton microscopy and fluorescence lifetime imaging microscopy (FLIM) to monitor metastasis and the tumor microenvironmentClinical & Experimental Metastasis 26:357–370.https://doi.org/10.1007/s10585-008-9204-0
-
A molecular insight into complement evasion by the staphylococcal complement inhibitor protein familyThe Journal of Immunology 183:2565–2574.https://doi.org/10.4049/jimmunol.0901443
-
Commonality of TRIM32 mutation in causing sarcotubular myopathy and LGMD2HAnnals of Neurology 57:591–595.https://doi.org/10.1002/ana.20441
-
Altered myogenesis and premature senescence underlie human TRIM32-related myopathyActa Neuropathologica Communications 7:30.https://doi.org/10.1186/s40478-019-0683-9
-
Analysis of glycolytic enzyme co-localization in Drosophila flight muscleJournal of Experimental Biology 206:2031–2038.https://doi.org/10.1242/jeb.00367
-
Coordinated metabolic transitions during Drosophila embryogenesis and the onset of aerobic glycolysisG3: Genes, Genomes, Genetics 4:839–850.https://doi.org/10.1534/g3.114.010652
-
Techniques to monitor glycolysisMethods in Enzymology 542:91–114.https://doi.org/10.1016/B978-0-12-416618-9.00005-4
-
TRIM-NHL proteins in development and diseaseSeminars in Cell & Developmental Biology 47-48:52–59.https://doi.org/10.1016/j.semcdb.2015.10.017
-
Fluorescence lifetime microscopy of NADH distinguishes alterations in cerebral metabolism in vivoBiomedical Optics Express 8:2368–2385.https://doi.org/10.1364/BOE.8.002368
-
Metabolic regulation of cell growth and proliferationNature Reviews Molecular Cell Biology 20:436–450.https://doi.org/10.1038/s41580-019-0123-5
-
Automated structure solution with the PHENIX suiteMethods in Molecular Biology 426:419–435.https://doi.org/10.1007/978-1-60327-058-8_28
Article and author information
Author details
Funding
National Institute of Arthritis and Musculoskeletal and Skin Diseases (R01AR060788)
- Erika R Geisbrecht
National Institute of General Medical Sciences (R35GM119557)
- Jason M Tennessen
National Institute of Arthritis and Musculoskeletal and Skin Diseases (R21AR073373)
- Erika R Geisbrecht
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jim Vigoreaux for sharing antibodies and Kasra X Ramyar, Joe McWhorter, and Samantha Gameros for technical assistance. We also thank James Cox for services at the University of Utah Metabolomics Core as well as Steve Hartson and Janet Rogers for services and guidance at the Oklahoma State University Proteomics and Mass Spectrometry Core. Special appreciation to the VDRC and BDSC for fly lines used in this study. This work was supported by grant through the National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS) to ERG (R21AR073373, R01AR060788, and R56AR060788). JMT is supported by a MIRA award from NIGMS (R35GM119557). X-ray diffraction data were collected at Southeast Regional Collaborative Access Team 22-BM beamline at the Advanced Photon Source, Argonne National Laboratory. Supporting institutions may be found at www.ser-cat.org/members.html. Use of the Advanced Photon Source was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract W-31–109-Eng-38*.
Copyright
© 2020, Bawa 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
-
- 3,564
- views
-
- 378
- downloads
-
- 25
- 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
-
- Cancer Biology
- Genetics and Genomics
Solitary fibrous tumors have gene expression signatures similar to those of neuroendocrine tumors.
-
- Cancer Biology
Gastric cancer (GC) is a major cause of cancer-related mortality worldwide. Despite the widespread recognition of tumor immunotherapy in treating unresectable GC, challenges, including ineffective immunotherapy and drug resistance, persist. Therefore, understanding the regulatory mechanisms of PD-L1, particularly in the context of super-enhancers (SEs) and zinc finger protein 36 ring finger protein-like 1 (ZFP36L1) RNA-binding protein, is crucial. In this study, we performed H3K27ac Cleavage Under Targets and Tagmentation (CUT&Tag) sequencing, investigated the heterogeneity of SEs between two GC subtypes with differential growth patterns, and revealed the immune escape signatures driven by ZFP36L1-SE in infiltrative GC through SEs inhibitors treatment. The regulation of ZFP36L1 to PD-L1 was evaluated by quantitative PCR, western blot, flow cytometry, and immunohistochemistry. Furthermore, we explored its regulatory mechanisms using a combination of molecular biology techniques, including luciferase reporter assay, GST/RNA pull-down, chromatin immunoprecipitation (ChIP)/RIP experiments, and in vivo functional assays. We demonstrated that ZFP36L1, driven by an SE, enhances IFN-γ-induced PD-L1 expression, with SPI1 identified as the specific transcription factor binding to ZFP36L1-SE. Mechanistically, ZFP36L1 binds to the adenylate uridylate-rich element in the 3ʹ untranslated region (3ʹUTR) of HDAC3 mRNA, exacerbating its mRNA decay, and thereby facilitating PD-L1 abnormal transcriptional activation. Collectively, our findings provide mechanistic insights into the role of the SPI1-ZFP36L1-HDAC3-PD-L1 signaling axis in orchestrating immune escape mechanisms in GC, thereby offering valuable insights into the potential targets for immune checkpoint therapy in GC management.