- Research article
- Open Access
Drosophila insulin and target of rapamycin (TOR) pathways regulate GSK3 beta activity to control Myc stability and determine Myc expression in vivo
BMC Biology volume 9, Article number: 65 (2011)
Genetic studies in Drosophila melanogaster reveal an important role for Myc in controlling growth. Similar studies have also shown how components of the insulin and target of rapamycin (TOR) pathways are key regulators of growth. Despite a few suggestions that Myc transcriptional activity lies downstream of these pathways, a molecular mechanism linking these signaling pathways to Myc has not been clearly described. Using biochemical and genetic approaches we tried to identify novel mechanisms that control Myc activity upon activation of insulin and TOR signaling pathways.
Our biochemical studies show that insulin induces Myc protein accumulation in Drosophila S2 cells, which correlates with a decrease in the activity of glycogen synthase kinase 3-beta (GSK3β ) a kinase that is responsible for Myc protein degradation. Induction of Myc by insulin is inhibited by the presence of the TOR inhibitor rapamycin, suggesting that insulin-induced Myc protein accumulation depends on the activation of TOR complex 1. Treatment with amino acids that directly activate the TOR pathway results in Myc protein accumulation, which also depends on the ability of S6K kinase to inhibit GSK3β activity. Myc upregulation by insulin and TOR pathways is a mechanism conserved in cells from the wing imaginal disc, where expression of Dp110 and Rheb also induces Myc protein accumulation, while inhibition of insulin and TOR pathways result in the opposite effect. Our functional analysis, aimed at quantifying the relative contribution of Myc to ommatidial growth downstream of insulin and TOR pathways, revealed that Myc activity is necessary to sustain the proliferation of cells from the ommatidia upon Dp110 expression, while its contribution downstream of TOR is significant to control the size of the ommatidia.
Our study presents novel evidence that Myc activity acts downstream of insulin and TOR pathways to control growth in Drosophila. At the biochemical level we found that both these pathways converge at GSK3β to control Myc protein stability, while our genetic analysis shows that insulin and TOR pathways have different requirements for Myc activity during development of the eye, suggesting that Myc might be differentially induced by these pathways during growth or proliferation of cells that make up the ommatidia.
Genetic studies in Drosophila have identified Myc as well as components of insulin and target of rapamycin (TOR) signaling pathways as key regulators of growth [1, 2]. Insulin and TOR signaling pathways are highly conserved. In Drosophila, binding of insulin-like peptides to the insulin receptor (InR) results in the activation and phosphorylation of chico, the ortholog of the insulin receptor substrates 1-4 (IRS1-4) . This event leads to the production of phosphatidylinositol-3,4,5-triphosphate (PIP3) by phosphoinositide 3-kinase (PI3K), a reaction that is counteracted by the lipid phosphatase and tensin homolog (PTEN) [4, 5]. PIP3 recruits several Ser/Thr kinases to the plasma membrane including Akt/PKB (Protein Kinase B) and PDK1 (3'-phosphoinositide-dependent protein kinase-1) . Activation of Akt results in the inhibition of glycogen synthase kinase 3-beta (GSK3β), a conserved kinase that not only controls energy metabolism by inactivation of glycogen synthase, but also regulates Wnt signaling by controlling βcatenin/armadillo  and Myc stability [9, 10]. Activation of Akt also inhibits tuberous sclerosis complex 1 and 2 (TSC1/2), a binary-complex that negatively regulates Rheb, a GTPase upstream of TOR kinase responsible for activation of TOR complex 1 . TOR is found in two complexes: TOR complex 1, which includes Raptor and LSt8 adaptor molecules, is sensitive to amino acids and is inhibited by rapamycin; and TOR complex 2, which is composed of LSt8 and Rictor adaptor molecules, and does not respond to amino acids or rapamycin [12, 13]. Activation of TOR complex 1 results in phosphorylation of ribosomal protein kinase p-70-S6 (S6K) on threonine 398, and of eukaryotic translation initiation factor 4E-binding protein 1(4E-BP1), thereby triggering protein synthesis and initiation of translation [14, 15]. Insulin and TOR activities are also balanced by a negative feedback mechanism that is activated when S6K is hyper-activated to counteract insulin activity. Under this condition, S6K phosphorylates IRS1-4/chico triggering its internalization and subsequent proteasomal degradation [16, 17]. This feedback mechanism is reduced in pathological conditions, such as TSC syndromes where cells carrying mutations in tsc1 or tsc2 display an abnormal increase in size and exhibit constitutive phosphorylation of S6K . In these cells, hyper-activation of S6K correlates with inactivation of GSK3β by phosphorylation of Serine 9, which results in c-Myc protein accumulation .
The Drosophila dmyc gene, called diminutive (dm), was identified as the sole ortholog of the human c-myc gene [19, 20]. Analysis of dm target genes revealed a prominent role for Myc in the regulation of genes controlling ribosomal biogenesis and protein synthesis [21–24], thus animals carrying dm hypomorphic mutations are smaller due to a reduction in their protein synthesis and cell size [25–28]. This phenotype is reminiscent of hypomorphic mutants for tor or s6k, components of the TOR signaling pathway [29–31]. Signaling through the InR controls growth; mutants for the Drosophila-insulin-like peptides or of the adaptor molecule chico/IRS are also smaller because of fewer and smaller cells . Recent genomic analysis in whole larvae showed a strong correlation between the targets of Myc and those of the TOR pathway; however, less overlap was found between the targets of Myc and those of PI3K signaling . Whether Myc acts downstream or in parallel to PI3K is not totally clear and previous observations in vivo, in cells of the imaginal discs, indicate that activation of PI3K, even though it influenced growth and decreased the G1 phase of the cell cycle, did not significantly alter Myc protein levels . However, in salivary glands, Myc expression was able to partially rescue the growth defect caused by lowering PI3K activity  suggesting that Myc activity might be acting downstream of PI3K signaling. Together, these data suggest that Myc could participate in the regulation of growth in response to insulin and TOR signaling, however a molecular link between Myc and these signaling pathways has not been clearly identified yet.
In this report, we provide evidence of a molecular mechanism for the Myc protein to be stabilized by insulin and amino acid signaling in Drosophila S2 cells, which converge to decrease the activity of GSK3β, a kinase responsible for Myc protein degradation. We found Myc protein regulated in vivo both by insulin and TOR pathways in epithelial cells of the imaginal discs. Using genetic analysis we demonstrate that Myc functions downstream of insulin and TOR to sustain proper growth of the eye during development.
Insulin induces Myc protein accumulation in DrosophilaS2 cells, with a mechanism that results in inactivation of GSK3β and is dependent on TOR complex 1 activity
Our previous observation that Myc protein stability in Drosophila S2 cells is reduced by GSK3β activity  led us to investigate if stimulation of insulin signaling, which inhibits GSK3β via Akt phosphorylation, could result in Myc protein accumulation. Treatment of Drosophila S2 cells with insulin induced an increase in Myc protein levels visible after 30 minutes of stimulation that was still detectable after 180 minutes of treatment (Figure 1A). This event was accompanied by a small increase in dmyc-RNA that peaked after 30 minutes and rapidly returned to baseline levels (Additional file 1). Myc protein accumulation by insulin was accompanied by phosphorylation of Akt on Ser 505, an event that correlated with phosphorylation of GSK3β on Ser 9 (Figure 1B), and was inhibited in the presence of the PI3K inhibitor wortmannin (Figure 1C). In order to analyze if GSK3β signaling contributes to insulin-induced Myc protein upregulation, S2 cells were treated with insulin in the presence of the GSK3β inhibitor lithium chloride (LiCl), or insulin was added to S2 cells expressing the GSK3β- kinase dead (KD) mutant, which was shown previously to reduce GSK3β activity .
In these experiments, addition of insulin increased Myc protein levels and this was accompanied by increased phosphorylation of GSK3β on Ser 9 (Figure 1D and 1E, lane 2). LiCl or expression of GSK3β-KD also increased endogenous Myc protein levels (Figure 1D and 1E, lane 3) as a result of inhibition of the endogenous GSK3β activity, which is known to control Myc protein stability . Blocking GSK3β activity with LiCl or using GSK3β-KD, along with the addition of insulin, led to an increase in Myc, which was comparable to that of insulin alone (Figure 1D and 1E, compare lane 2 and 4).
Activation of insulin signaling enhances TOR activity by releasing the negative feedback of Akt on TSC1/2 [34–36]. We therefore treated S2 cells with rapamycin, an inhibitor of TOR complex 1 (TORC1), to analyze the contribution of TOR to insulin-mediated Myc upregulation. These experiments showed that rapamycin suppresses Myc protein accumulation by insulin (Figure 1F), and similar results were obtained for the regulation of dmyc-RNA (Additional file 2). These data also showed that rapamycin reduced endogenous Myc protein levels without affecting dmyc-RNA, suggesting that TOR signaling might regulate Myc protein stability. To better understand this mechanism we blocked the proteasome using MG132 and analyzed Myc protein level upon rapamycin treatment. These data showed that Myc protein degradation in the presence of rapamycin was completely suppressed by MG132 (Figure 1G) suggesting that TOR activity regulates Myc protein stability by mediating its degradation through the ubiquitin-proteosomal pathway.
Amino acids and/or TOR signaling increases Myc protein stability by GSK3β inhibition in DrosophilaS2 cells
TOR signaling is activated by insulin and also by amino acids (AAs) [35, 36]. To determine whether AAs directly controlled Myc protein accumulation, we serum-starved S2 cells and then performed a complete amino acid starvation by subsequently bathing them in an AA-free medium without serum. After 30 minutes, AAs were added back to the cells for the indicated times and Myc protein levels were analyzed by western blotting. As shown in Figure 2A, treatment with AAs increased Myc protein levels, which peaked between 60 and 90 minutes after treatment. Quantitative RT-PCR analysis showed that dmyc-mRNA was not significantly affected (Figure 1B in Additional file 1), suggesting that TOR signaling regulates Myc protein mainly at its post-translational level. AA starvation resulted in a reduction of Myc protein levels (Figure 2B, compare lane 1 and lane 2), which was increased by adding AAs back to the medium (Figure 2B, lane 3). This correlated with an increase in GSK3β phosphorylation of Ser 9. Myc upregulation by AAs was significantly reduced in the presence of rapamycin (Figure 2B, lane 4). In order to analyze if GSK3β activity contributes to TOR-induced Myc protein upregulation, we stimulated S2 cells with AAs in the presence of LiCl and analyzed whether AAs could still induce Myc protein accumulation in this condition.
These experiments showed that addition of LiCl to an AA-free medium slightly increased Myc protein level (Figure 2C, lane 2) but to a lesser extent than in medium with amino acids (compare with Figure 1D and 1E, lane 3). Addition of AAs increased Myc protein levels and was accompanied by an increase in phosphorylation of S6K on Thr 398 that correlated with phosphorylation of GSK3β on Ser 9 (Figure 2C, lane 3). Addition of LiCl together with AAs did not further increase Myc protein levels, an outcome similar to that observed with AA treatment alone (Figure 2C, compare lane 3 and 4). Similar results were obtained for the phosphorylation of S6K on Thr 398 or GSK3β on Ser 9.
We then analyzed if ectopic expression of Rheb and S6K resulted in accumulation of Myc protein. In these experiments the epitope-tagged form of Rheb (MYC-Rheb) and of S6K (HA-S6K)  were expressed in S2 cells and endogenous Myc protein was analyzed by western blot. These experiments showed that while S6K alone had a very small effect on Myc protein levels (Figure 2D), co-expression of S6K with Rheb was able to substantially induce Myc protein accumulation, an event that correlated with phosphorylation of S6K on Thr 398 and GSK3β on Ser 9. In order to understand if accumulation of Myc protein induced by Rheb and S6K was regulated at the post-transcriptional level, we co-expressed Rheb and S6K together with an epitope HA-tagged form of Myc . In these experiments the UAS/Gal4 system  was used to express UAS-HA-Myc in a stable S2-line that constitutively expressed Gal4 under the tubulin promoter (S2-tub > Gal4). Rheb and S6K were co-expressed with HA-Myc and their protein levels were analyzed using anti-tag antibodies. These experiments showed that expression of S6K alone was not sufficient to induce significant changes in Myc protein levels, however expression of Rheb alone resulted in the accumulation of HA-Myc protein and this effect was accompanied by phosphorylation of GSK3β on Ser 9 (Figure 2E).
Taken together, these data show that TOR signaling controls Myc protein stability at the post-transcriptional level, which correlates with inhibition of GSK3β by its phosphorylation on Ser 9.
Modulation of insulin and TOR signaling regulates Myc protein levels in epithelial cells of the wing imaginal discs
Our next step was to analyze if components of the insulin and TOR signaling pathways regulate Myc protein levels in vivo in cells from the wing imaginal discs. Because it was previously shown that reduction of components of the insulin and of TOR pathways in vivo resulted in clones of reduced size [3, 31, 39], we decided to take advantage of a conditional gene expression technique where the Act > CD2 > PR-Gal4 construct was used to temporally induce the expression of UAS-transgenes under the actin promoter through the addition of the synthetic steroid mifepristone . Under these conditions, we generated clones expressing UAS-Dp110 or UAS-PTEN to modulate insulin signaling, or we activated the TOR pathway by expressing rheb using the UAS-RhebAV4 line, derived from insertion of a P element in the rheb locus . The Myc protein level was analyzed in clones of cells from wing imaginal discs of third instar larvae after five hours of induction. Myc protein at this stage of development (Figure 3A, in red) is predominantly expressed within the cells of the presumptive notum and in the wing pouch, with the exception of the area of the hinge and a stripe of cells in the zone of non-proliferative cells (ZNC) located along the dorsal-ventral boundary, where Myc is transcriptionally repressed by Wingless activity [26, 42, 43]. These experiments showed that clones expressing Dp110, marked by co-expression of GFP, accumulate Myc protein (Figure 3B, red), which was also visible in cells from a clone that was generated within the ZNC, where normally Myc expression is repressed (arrow). In contrast, Myc protein was visibly reduced in clones overexpressing PTEN, a negative regulator of the insulin pathway (Figure 3C). Upregulation of TOR signaling, using UAS-RhebAV , also induced the accumulation of Myc protein (Figure 3D); on the contrary Myc protein was reduced in clones expressing TOR TED , a mutant of TOR that functions as a dominant negative for TORC1 activity  (Figure 3E). Interestingly, these experiments revealed a strong non-autonomous accumulation of Myc protein in cells neighboring the clones, particularly visible in Dp110 and PTEN clones and to a lesser extent in Rheb AV4 and TOR TED clones (arrowhead). This was more evident when clones where positioned along the dorsal-ventral axis of the wing disc, suggesting that differences in Myc expression in cells at the border of the dorsal-ventral axis may induce changes that resulted in the regulation of Myc levels non-autonomously.
Genetic interaction of Myc with members of the insulin and TOR signaling pathways in the adult eye
Given our biochemical and cellular evidence indicating that both the insulin and TOR pathways regulate Myc protein levels, we used a genetic approach to analyze the contribution of Myc to the growth exerted by these pathways. Using the UAS/Gal4 system in combination with flip-out techniques [45, 46], we expressed components of the insulin or TOR pathways in the eye and analyzed their relative effect on growth by comparing the size and number of ommatidia from animals with different dm genetic backgrounds. In these experiments we recombined the tubulin-FRT-dmyc-FRT-Gal4 line, which we previously used to rescue viability of dmP0 and lethality of dm4 mutants, with the eyeless-Flp transgene . This line contains the dmyc-cDNA located in a removable cassette between the tubulin promoter and Gal4. Upon FLP-mediated excision of dmyc-cDNA, the tubulin promoter drives Gal4 expression allowing for the expression of the UAS transgenes in the eyeless compartment (eye and antenna). This chromosome (hereafter called ey) was recombined in wild-type animals (ey-dm+ ) and in flies carrying the hypomorphic dmP0 or the null dm4 alleles (hereafter called ey-dmP0 or ey-dm4 ). The size and number of the ommatidia in the eyes from control wild-type (ey-dm+ ) was compared to that from ey-dmP0 or ey-dm4 mutant animals expressing the different transgenes. The statistical significance of these data, across the different dm genetic backgrounds was calculated using the two-tailed z test (Appendix 1 in Additional file 3). Our first approach was to analyze if Myc activity modulates ommatidial size changes induced by components of the insulin pathway. This analysis revealed that even though expression of UAS-Dp110 increased the size of the ommatidia by 38% in a wild-type dm+ background (Figure 4 and Table 1; P < 0.001, using a Student's t test) this effect was not dependent on Myc expression (Table S1 in Additional file 4; P = 0.4292, using a two-tailed z test) and only in dm4/Y null animals did Dp110 exhibit a weak but significant inhibition of its effect in increasing the size of the ommatidia (Table S1 in Additional file 4; P = 0.0367, using a two-tailed z test). By contrast, the increase in the total number of the ommatidia induced by Dp110 in ey-dm+/Y animals (Table 1; P < 0.001, t test) was significantly reduced in dmP0 /Y and dm4 /Y flies (Table S1 in Additional file 4; P < 0.001, two-tailed z test). Reduction of insulin signaling by PTEN showed a significant decrease in the size of the ommatidia in ey-dm+ /Y animals (92%) (Figure 4, Table 1, P < 0.001, t test), and this effect was more pronounced in ey-dmP0/Y and ey-dm4/Y animals, where the size of the ommatidia was reduced to 64% and 67%, respectively (Table 1; P < 0.001, t test). PTEN also significantly reduced the total number of the ommatidia in the eyes of ey-dm+ /Y and dmP0 /Y flies (Table 1; P < 0.001, t test). Statistical analysis of relative growth across the different dm genetic backgrounds showed that this reduction was highly significant for both dm backgrounds (Table S1 in Additional file 4; P < 0.001, two-tailed z test). These data suggest that Myc activity significantly contributes to insulin signaling-induced growth in the eye.
Genetic analysis using components of TOR signaling showed that, while expression of p70-S6K moderately affected the size of the ommatidia (Table 1, P < 0.05, t test), expression of the UAS-RhebAV4 allele  showed an 89% increase of ommatidia size in ey-dm+ /Y flies (Table 1, P < 0.001, t test) and this effect was drastically reduced in ey-dmP0/Y and ey-dm4/Y animals (Table 1, P < 0.001, t test). This event was statistically significant as demonstrated by our analysis of the relative size increases across the different dm genetic backgrounds (Table S1 in Additional file 4; P < 0.001, z test).
In summary these data suggests that Myc functions downstream of TOR signaling to control the growth of the size of the ommatidia.
While performing these experiments we noticed that expression of components of TOR signaling, and in particular of the strong RhebAV4 allele, had a significant negative effect on the total number of ommatidia (Table 1; P < 0.001, t test). Moreover, this effect was rescued by reducing dmyc levels (Figure 4, Table 1). To understand the molecular mechanisms that caused Rheb to reduce the ommatidia number, imaginal discs from third instar larvae expressing UAS- RhebAV4 transgenes, were examined for defects in cell proliferation or for increased cell death. Imaginal eye discs from ey-dmP0/Y or wild-type ey-dm+/Y animals carrying the UAS-RhebAV4 transgene were subjected to bromodeoxyuridine (BrdU) labeling to detect DNA replication (S phase), or immunostained with anti-active caspase-3 to detect apoptotic cells. This analysis revealed that, while no significant changes were observed in the pattern of BrdU labeling between the different genotypes (Additional file 6), a significant increase in the number of caspase-3 positive cells in the antennal and eye imaginal discs of ey-dm+/Y; UAS-RhebAV4 /+ larvae was seen, which was significantly reduced in ey-dmP0/Y; UAS-RhebAV4 /+ animals (Additional file 7). This highlights a potential mechanism for TOR signaling to induce cell death when growth is in excess.
Previous studies in vertebrates have indicated a critical function for Myc downstream of growth factor signaling including insulin-like growth factor, insulin and TOR pathways [18, 48–50]. In Drosophila, despite a few notes that Myc transcriptional activity acts downstream of insulin and TOR pathways [23, 24], no clear molecular mechanisms linking these pathways to Myc have been elucidated yet.
We previously demonstrated that inhibition of GSK3β prevents Myc degradation by the proteasome pathway . In this report, we further unravel the pathways that control Myc protein stability and show that signaling by insulin and TOR induce Myc protein accumulation by regulating GSK3β activity in S2 cells. GSK3β is a constitutively active kinase that is regulated by multiple signals and controls numerous cellular processes . With our biochemical data we propose that GSK3β acts as a common point where insulin and TOR signaling converge to regulate Myc protein stability (Figure 5). In particular, we showed that activation of insulin signaling induces activation of Akt, an event that is accompanied by GSK3β phosphorylation on Ser 9 that causes its inactivation and Myc protein to stabilize (Figure 1B). Interestingly, insulin-induced Myc protein accumulation, when GSK3β activity was blocked by the presence of LiCl or by expression of GSK3β-KD, was similar to that obtained with insulin alone. Since we showed that activation of insulin signaling leads to GSK3β inhibition and to an increase in Myc protein, if insulin and GSK3β signaling were acting independently, we would expect that activation of insulin signaling concomitantly with the inhibition of GSK3β activity would result in a higher level of Myc than that obtained with insulin or LiCl alone. Our results instead showed a similar level of Myc protein accumulation with insulin in the presence of GSK3β inhibitors as compared to insulin alone (Figure 1D and 1E, compare lane 2 and 4), supporting the hypothesis that GSK3β and insulin signaling, at least in our experimental condition, depend on each other in the mechanism that regulates Myc protein stability.
In a similar biochemical approach, we analyzed the effect of AAs on Myc protein stability and how TOR signaling is linked to mechanisms that inactivate GSK3β to stabilize Myc protein in S2 cells. In these experiments we were able to demonstrate that AAs increased Myc protein stability, and we also showed that treatment with rapamycin, an inhibitor of TORC1, reduced insulin-induced Myc upregulation. The reduction of Myc protein accumulation by rapamycin was blocked by inhibition of the proteasome pathway, linking TOR signaling to the pathway that controls Myc protein stability (Figure 1F). TORC1 is a central node for the regulation of anabolic and catabolic processes and contains the central enzyme Rheb-GTPase, which responds to amino acids by activating TOR kinase to induce phosphorylation of p70-S6K and 4E-BP1 [14, 15]. Our analysis of the molecular mechanisms that act downstream of TOR to regulate Myc stability shows that AA treatment induces p70-S6K to phosphorylate GSK3β on Ser 9, an event that results in its inactivation and accumulation of Myc protein (Figure 2).
Reducing GSK3β activity with LiCl, in medium lacking AAs, resulted in a slight increase in Myc protein levels (Figure 2C, lane 1 and 2). Adding back AAs lead to a substantial increase in Myc protein levels, which did not further increase when AAs where added to cells in the presence of the GSK3β inhibitor LiCl (Figure 2C, lane 3 and 4). These events were accompanied by phosphorylation of S6K on Thr 398, which correlated with phosphorylation of GSK3β on Ser 9. From these experiments we conclude that TOR signaling also converges to inhibit GSK3β activity to regulate Myc protein stability (Figure 5). However, we need to point out that since AAs alone increased Myc protein levels to a higher extent than that observed with LiCl alone (Figure 2C, compare lane 2 and 3), our experiments also suggest that Myc protein stability by TOR signaling is not solely directed through the inhibition of GSK3β activity, but other events and/or pathways contribute to Myc regulation. In conclusion, our biochemical experiments demonstrate that GSK3β acts downstream of insulin and TOR pathways to control Myc stability, however we do not exclude that other pathways may control Myc protein stability upon insulin and amino acids stimulation in S2 cells.
Reduction of insulin and TOR signaling in vivo reduces cell size and proliferation, and clones mutant for chico, the Drosophila orthologue of IRS1-4, or for components of TOR signaling, are smaller due a reduction in size and the number of cells [29, 31, 39, 51]. Our experiments showed that reducing insulin signaling by expression of PTEN or using TOR TED , a dominant negative form of TOR, decreased Myc protein levels in clones of epithelial cells of the wing imaginal discs, while the opposite was true when these signals were activated using Dp110 or RhebAV4 (Figure 3). Those experiments suggested that the mechanism of regulation of Myc protein by insulin and TOR pathways was conserved also in vivo in epithelial cells of the larval imaginal discs.
During these experiments we also noted that Myc protein was induced in the cells surrounding and bordering the clones (non-autonomously), particularly when clones where positioned along the dorsal-ventral axis of the wing disc. This upregulation of Myc protein was not restricted to components of the insulin signaling pathway since we also observed it in cells surrounding the clones mutant for components of the Hippo pathway  or for the tumor suppressor lethal giant larvae (lgl), which upregulates Myc protein cell-autonomously . We suspect that this non-autonomous regulation of Myc may be induced by a novel mechanism that controls proliferation of cells when 'growth' is unbalanced. We can speculate that clones with different growth rates, caused by different Myc levels, might secrete factors to induce Myc expression in neighboring cells. As a consequence, these Myc-expressing cells will speed up their growth rate in an attempt to maintain proliferation and tissue homeostasis. Further analysis is required to identify the mechanisms responsible for this effect.
In order to distinguish if Myc activity was required downstream of insulin and TOR signaling to induce growth, we performed a genetic analysis. The ability to induce growth and proliferation was measured in the eye by measuring the size and number of the ommatidia from animals expressing members of the insulin and TOR pathways in different dm genetic background (dm+, dmP0 and dm4 ). Our data showed that Dp110 increased the size and number of the ommatidia, however only the alteration in the total number was dependent on dm levels. These data suggest that Myc is required downstream of insulin pathway to achieve the proper number of ommatidia. However, when insulin signaling was reduced by PTEN, a significant decrease in the size of ommatidia was seen and it was dependent on dm expression levels, suggesting that Myc activity is limiting for ommatidial size and number. Activation of TOR signaling induces growth , and our genetic analysis showed that Myc significantly contributes to the size of the ommatidial cells thus suggesting that Myc acts downstream of TOR pathway to control growth.
Recent genomic analysis showed a strong correlation between the targets of Myc and those of the TOR pathway , implying that they may share common targets. In support of this observation our mosaic analysis with a repressible cell marker (MARCM) experiments in the developing wing disc showed that overexpression of Myc partially rescues the growth disadvantage of clones mutant for the hypomorphic Rheb7A1 allele (Additional file 8), further supporting the idea that Myc acts downstream of TOR to activate targets that control growth in these clones.
Our genetic interaction revealed a stronger dependence on Myc expression when Rheb was used as opposed to S6K (Table 1). A possible explanation for this difference could lie in the fact that S6K is not capable of auto-activation of its kinase domain unless stimulated by TOR kinase. TOR activity is dependent on its upstream activator Rheb; consequently the enzymatic activity of the Rheb/GTPase is the limiting factor that influences S6K phosphorylation and therefore capable of maximizing its activity .
Interestingly, these experiments also showed that activation of TOR signaling has a negative effect on the number of ommatidia, and this correlates with the ability of Rheb AV4 to induce cell death during the development of the eye imaginal disc. Rheb-induced cell death was rescued in a dmP0 mutant background, which led us to speculate that 'excessive' protein synthesis, triggered by overexpression of TOR signaling, could elicit a Myc-dependent stress response, which induces apoptosis. Alternatively, high protein synthesis could result in an enrichment of misfolded proteins  that may result in a stress response and induces cell death. Further analysis is required to delineate the mechanisms underlying this process.
Our analyses provide novel genetic and biochemical evidences supporting a role for Myc in the integration of the insulin and TOR pathway during the control of growth, and highlights the role of GSK3β in this signaling. We found that insulin signaling inactivates GSK3β to control Myc protein stability, and a similar biochemical regulation is also shared by activation of the TOR pathways. In support of this data, a recent genomic analysis in whole larvae showed a strong correlation between the targets of Myc and those of the TOR pathway; however, less overlap was found between the targets of Myc and those of PI3K signaling .
Statistical analysis applied to our genetic interaction experiments revealed that, in the Drosophila eye, proliferation induced by activation of the insulin pathway is sensitive to variations in Myc levels, while a significant interaction was seen mostly when TOR increased cell size. Our data therefore suggests that there is a correlation between Myc and the InR signaling and it is expected that the InR pathway also shares some transcriptional targets with Myc. Indeed, we found an overlap between the targets induced by insulin and Myc in Drosophila S2 cells (PB, unpublished data) and these targets have also been reported in transcriptome analyses in the fat body upon nutritional stress , suggesting that Myc acts downstream of InR/PI3K and TOR signaling and that this interaction might be specific to some tissues or in a particular metabolic state of the cell.
Fly stocks were obtained from the Bloomington stock center including w; UAS-RhebAV4/TM6b which is derived from P-element inserted in the rheb locus , UAS-Dp110 , UAS-PTEN , with the exception of w; UAS-TORTED a mutant that expresses the toxic extended domain of TOR protein (Thomas Neufeld, UMN, MN, USA), UAS-HA-Myc , yw; UAS-Rheb and w; FRT82, Rheb7A1/TM6b (Hugo Stocker, ETH-ZH, Switzerland), w; hs-Flp; tub-Gal4, UAS-GFP; FRT82 (hs-CD2 y+)-tub-Gal80 (w+) (Myriam Zecca, Columbia University, NY, USA); UAS-GFP;Ay-Gal4:PR/TM6b (Kenneth Irvine, Rutgers, NJ, USA ), ey-dm4 = yw dm4 tubulin-FRT-dmyc-cDNA-FRT-Gal4 ey-Flp/Y , ey-dmP0 = yw dmP0tubulin-FRT-dmyc-cDNA-FRT-Gal4, ey-Flp/Y, and ey-dm+ = yw dm+ tubulin-FRT-dmyc-cDNA-FRT-Gal4, ey-Flp/Y.
Cell culture and western blot
Drosophila Schneider S2 cells or the line expressing the plasmid tubulin-Gal4 (S2-tub-Gal4)  were grown at 25°C using Schneider medium (Invitrogen, GIBCO, Carlsbad, CA, USA) supplemented with 10% heat-inactivated FCS and 100 I.U. of penicillin/streptomycin (Invitrogen, GIBCO, Carlsbad, CA, USA). Insulin and AA treatments: S2 cells were serum-starved in 0.5% serum for 12 h and then various chemical inhibitors were added to the medium with or without insulin (from porcine, Sigma, St Louis, MO, USA 1 μM final concentration). For AA treatment, S2 cells were serum-starved as described above and further incubated for 30 min in medium without amino acids, made from the recipe of Drosophila Schneider medium (Invitrogen, GIBCO, Carlsbad, CA, USA) containing 0.5% dialyzed FCS. A 2× solution containing the same concentration of AAs as the Drosophila Schneider medium was added to induce an AA response for the indicated times. Cells were washed in PBS and lysed in a buffer containing 50 mM HEPES, pH 7.4; 250 mM NaCl; 1 mM EDTA, 1% Triton and protease and phosphatases inhibitors (Roche, Mannheim, Germany). Protein concentration was measured using the Bio-Rad protein assay. Western blot analysis was performed using the following primary antibodies: anti-Drosophila Myc mAb ; anti-HA (Roche, Mannheim, Germany), anti-actin mAb and anti vinculin mAb (Sigma, St. Louis, MO, USA); anti-Drosophila phospho-Akt, anti-phospho-GSK3β anti-phospho-S6K and anti GSK3β (Cell Signaling Technology inc. Danvers, MA, USA). Chemical inhibitors were used at 1 μM for rapamycin, 100 nM for wortmannin, and 50 mM for LiCl. Plasmids encoding HA-S6K (pAct-HA-s6k) and Myc-Rheb (pAct-Myc-Rheb), were a gift from Duojia J. Pan . Transfection of S2 cells or of S2-tub-Gal4 was performed using Cellfectin reagent (Invitrogen, Carlsbad, CA, USA). The efficiency of transfection in S2 cells was analyzed by co-transfecting tub-Gal4 with UAS-GFP plasmids, and was found in the range of 30% to 60% depending on the experiment. Experiments were repeated at least three times with similar results.
Generation of inducible Flp-out clones
Flp-out clones were generated using UAS-GFP;AyGal4:PR/TM6b flies  by heat shock for 20 min at 37°C. Heat shock was performed 48 h after egg laying (AEL). Gal4:PR was activated by transferring larvae to instant food (Instant Drosophila Medium, Connecticut Valley Biological, MA, USA) containing RU486 (mifepristone, Sigma, St. Louis, MO, USA). Two grams of instant food were mixed with 7 mL of water previously supplemented with RU486, resulting in a final concentration of 20 μg/ml of mifepristone. Immunostaining protocol is reported in the Additional file 9: Supplementary Material and Methods.
Analysis of ommatidial size and number
Flies were reared at 25°C under reproducible growth conditions and age matched (three-day-old males) before determining the ommatidial size and number. Total ommatidial number was counted on scanning electron micrographs from animals of the indicate genotypes. Size of the ommatidia was calculated by measuring the area of 20 ommatidia located at the center of the eye using Adobe Photoshop 7.0, as described previously . At least eight animals of each genotype were measured.
To perform this statistical analysis we considered two components: an estimation of the relative size and number of the ommatidia and their variability, and a standard two-sided z test to establish the significance of the differences of our data within the different dm genetic backgrounds. Means comparison between groups of animals was made using a Student's t test. Data in Table 1 are presented as means. Standard deviations and the number of animals are represented in parenthesis. A two-tailed z test analysis was performed to calculate the p values for the hypothesis that the relative growth of the ommatidial size and number for the various phenotypes are equal (see Appendix 1 in Additional file 3 and Table S1 in Additional file 4).
after egg laying
fetal calf serum
green fluorescent protein
glycogen synthase kinase 3-beta
GSK3β kinase dead
mosaic analysis with a repressible cell marker
phosphate buffered saline
phosphatase and tensin homolog
Ras homolog enriched in brain
reverse transcription polymerase chain reaction
p70-S6 ribosomal protein kinase
target of rapamycin
TOR complex 1
tuberous sclerosis complex
zone of non-proliferative cells
eukaryotic translation initiation factor 4E-binding protein 1.
Bellosta P, Gallant P: Myc Function in Drosophila. Genes Cancer. 2010, 1 (6): 542-546. 10.1177/1947601910377490.
Grewal SS: Insulin/TOR signaling in growth and homeostasis: a view from the fly world. Int J Biochem Cell Biol. 2009, 41 (5): 1006-1010. 10.1016/j.biocel.2008.10.010.
Bohni R, Riesgo-Escovar J, Oldham S, Brogiolo W, Stocker H, Andruss BF, Beckingham K, Hafen E: Autonomous control of cell and organ size by CHICO, a Drosophila homolog of vertebrate IRS1-4. Cell. 1999, 97 (7): 865-875. 10.1016/S0092-8674(00)80799-0.
Huang H, Potter CJ, Tao W, Li DM, Brogiolo W, Hafen E, Sun H, Xu T: PTEN affects cell size, cell proliferation and apoptosis during Drosophila eye development. Development. 1999, 126 (23): 5365-5372.
Gao X, Neufeld TP, Pan D: Drosophila PTEN regulates cell growth and proliferation through PI3K-dependent and -independent pathways. Dev Biol. 2000, 221 (2): 404-418. 10.1006/dbio.2000.9680.
Stocker H, Andjelkovic M, Oldham S, Laffargue M, Wymann MP, Hemmings BA, Hafen E: Living with lethal PIP3 levels: viability of flies lacking PTEN restored by a PH domain mutation in Akt/PKB. Science. 2002, 295 (5562): 2088-2091. 10.1126/science.1068094.
Alessi DR, James SR, Downes CP, Holmes AB, Gaffney PR, Reese CB, Cohen P: Characterization of a 3-phosphoinositide-dependent protein kinase which phosphorylates and activates protein kinase Balpha. Curr Biol. 1997, 7 (4): 261-269. 10.1016/S0960-9822(06)00122-9.
Xu C, Kim NG, Gumbiner BM: Regulation of protein stability by GSK3 mediated phosphorylation. Cell Cycle. 2009, 8 (24): 4032-4039. 10.4161/cc.8.24.10111.
Sears R, Nuckolls F, Haura E, Taya Y, Tamai K, Nevins JR: Multiple Ras-dependent phosphorylation pathways regulate Myc protein stability. Genes Dev. 2000, 14 (19): 2501-2514. 10.1101/gad.836800.
Galletti M, Riccardo S, Parisi F, Lora C, Saqcena MK, Rivas L, Wong B, Serra A, Serras F, Grifoni D, Pelicci P, Jiang J, Bellosta P: Identification of domains responsible for ubiquitin-dependent degradation of dMyc by glycogen synthase kinase 3beta and casein kinase 1 kinases. Mol Cell Biol. 2009, 29 (12): 3424-3434. 10.1128/MCB.01535-08.
Sengupta S, Peterson TR, Sabatini DM: Regulation of the mTOR complex 1 pathway by nutrients, growth factors, and stress. Mol Cell. 2010, 40 (2): 310-322. 10.1016/j.molcel.2010.09.026.
Loewith R, Jacinto E, Wullschleger S, Lorberg A, Crespo JL, Bonenfant D, Oppliger W, Jenoe P, Hall MN: Two TOR complexes, only one of which is rapamycin sensitive, have distinct roles in cell growth control. Mol Cell. 2002, 10 (3): 457-468. 10.1016/S1097-2765(02)00636-6.
Kim E, Goraksha-Hicks P, Li L, Neufeld TP, Guan KL: Regulation of TORC1 by Rag GTPases in nutrient response. Nat Cell Biol. 2008, 10 (8): 935-945. 10.1038/ncb1753.
Inoki K, Ouyang H, Li Y, Guan KL: Signaling by target of rapamycin proteins in cell growth control. Microbiol Mol Biol Rev. 2005, 69 (1): 79-100. 10.1128/MMBR.69.1.79-100.2005.
Mamane Y, Petroulakis E, LeBacquer O, Sonenberg N: mTOR, translation initiation and cancer. Oncogene. 2006, 25 (48): 6416-6422. 10.1038/sj.onc.1209888.
Haruta T, Uno T, Kawahara J, Takano A, Egawa K, Sharma PM, Olefsky JM, Kobayashi M: A rapamycin-sensitive pathway down-regulates insulin signaling via phosphorylation and proteasomal degradation of insulin receptor substrate-1. Mol Endocrinol. 2000, 14 (6): 783-794. 10.1210/me.14.6.783.
Radimerski T, Montagne J, Rintelen F, Stocker H, van der Kaay J, Downes CP, Hafen E, Thomas G: dS6K-regulated cell growth is dPKB/dPI(3)K-independent, but requires dPDK1. Nat Cell Biol. 2002, 4 (3): 251-255. 10.1038/ncb763.
Zhang HH, Lipovsky AI, Dibble CC, Sahin M, Manning BD: S6K1 regulates GSK3 under conditions of mTOR-dependent feedback inhibition of Akt. Mol Cell. 2006, 24 (2): 185-197. 10.1016/j.molcel.2006.09.019.
Gallant P, Shiio Y, Cheng PF, Parkhurst SM, Eisenman RN: Myc and Max homologs in Drosophila. Science. 1996, 274 (5292): 1523-1527. 10.1126/science.274.5292.1523.
Schreiber-Agus N, Stein D, Chen K, Goltz JS, Stevens L, DePinho RA: Drosophila Myc is oncogenic in mammalian cells and plays a role in the diminutive phenotype. Proc Natl Acad Sci USA. 1997, 94 (4): 1235-1240. 10.1073/pnas.94.4.1235.
Orian A, van Steensel B, Delrow J, Bussemaker HJ, Li L, Sawado T, Williams E, Loo LW, Cowley SM, Yost C, Pierce S, Edgar BA, Parkhurst SM, Eisenman RN: Genomic binding by the Drosophila Myc, Max, Mad/Mnt transcription factor network. Genes Dev. 2003, 17 (9): 1101-1114. 10.1101/gad.1066903.
Hulf T, Bellosta P, Furrer M, Steiger D, Svensson D, Barbour A, Gallant P: Whole-genome analysis reveals a strong positional bias of conserved dMyc-dependent E-boxes. Mol Cell Biol. 2005, 25 (9): 3401-3410. 10.1128/MCB.25.9.3401-3410.2005.
Li L, Edgar BA, Grewal SS: Nutritional control of gene expression in Drosophila larvae via TOR, Myc and a novel cis-regulatory element. BMC Cell Biol. 2010, 11: 7-10.1186/1471-2121-11-7.
Teleman AA, Hietakangas V, Sayadian AC, Cohen SM: Nutritional control of protein biosynthetic capacity by insulin via Myc in Drosophila. Cell Metab. 2008, 7 (1): 21-32. 10.1016/j.cmet.2007.11.010.
Benassayag C, Montero L, Colombie N, Gallant P, Cribbs D, Morello D: Human c-Myc isoforms differentially regulate cell growth and apoptosis in Drosophila melanogaster. Mol Cell Biol. 2005, 25 (22): 9897-9909. 10.1128/MCB.25.22.9897-9909.2005.
Johnston LA, Prober DA, Edgar BA, Eisenman RN, Gallant P: Drosophila myc regulates cellular growth during development. Cell. 1999, 98 (6): 779-790. 10.1016/S0092-8674(00)81512-3.
Pierce SB, Yost C, Britton JS, Loo LW, Flynn EM, Edgar BA, Eisenman RN: dMyc is required for larval growth and endoreplication in Drosophila. Development. 2004, 131 (10): 2317-2327. 10.1242/dev.01108.
Grewal SS, Li L, Orian A, Eisenman RN, Edgar BA: Myc-dependent regulation of ribosomal RNA synthesis during Drosophila development. Nat Cell Biol. 2005, 7 (3): 295-302. 10.1038/ncb1223.
Oldham S, Montagne J, Radimerski T, Thomas G, Hafen E: Genetic and biochemical characterization of dTOR, the Drosophila homolog of the target of rapamycin. Genes Dev. 2000, 14 (21): 2689-2694. 10.1101/gad.845700.
Montagne J, Stewart MJ, Stocker H, Hafen E, Kozma SC, Thomas G: Drosophila S6 kinase: a regulator of cell size. Science. 1999, 285 (5436): 2126-2129. 10.1126/science.285.5436.2126.
Zhang H, Stallock JP, Ng JC, Reinhard C, Neufeld TP: Regulation of cellular growth by the Drosophila target of rapamycin dTOR. Genes Dev. 2000, 14 (21): 2712-2724. 10.1101/gad.835000.
Prober DA, Edgar BA: Interactions between Ras1, dMyc, and dPI3K signaling in the developing Drosophila wing. Genes Dev. 2002, 16 (17): 2286-2299. 10.1101/gad.991102.
Yanagawa S, Lee JS, Haruna T, Oda H, Uemura T, Takeichi M, Ishimoto A: Accumulation of Armadillo induced by Wingless, Dishevelled, and dominant-negative Zeste-White 3 leads to elevated DE-cadherin in Drosophila clone 8 wing disc cells. J Biol Chem. 1997, 272 (40): 25243-25251. 10.1074/jbc.272.40.25243.
Avruch J, Hara K, Lin Y, Liu M, Long X, Ortiz-Vega S, Yonezawa K: Insulin and amino-acid regulation of mTOR signaling and kinase activity through the Rheb GTPase. Oncogene. 2006, 25 (48): 6361-6372. 10.1038/sj.onc.1209882.
Tee AR, Fingar DC, Manning BD, Kwiatkowski DJ, Cantley LC, Blenis J: Tuberous sclerosis complex-1 and -2 gene products function together to inhibit mammalian target of rapamycin (mTOR)-mediated downstream signaling. Proc Natl Acad Sci USA. 2002, 99 (21): 13571-13576. 10.1073/pnas.202476899.
Lizcano JM, Alrubaie S, Kieloch A, Deak M, Leevers SJ, Alessi DR: Insulin-induced Drosophila S6 kinase activation requires phosphoinositide 3-kinase and protein kinase B. Biochem J. 2003, 374 (Pt 2): 297-306.
Zhang Y, Gao X, Saucedo LJ, Ru B, Edgar BA, Pan D: Rheb is a direct target of the tuberous sclerosis tumour suppressor proteins. Nat Cell Biol. 2003, 5 (6): 578-581. 10.1038/ncb999.
Brand AH, Manoukian AS, Perrimon N: Ectopic expression in Drosophila. Methods Cell Biol. 1994, 44: 635-654.
Stocker H, Radimerski T, Schindelholz B, Wittwer F, Belawat P, Daram P, Breuer S, Thomas G, Hafen E: Rheb is an essential regulator of S6K in controlling cell growth in Drosophila. Nat Cell Biol. 2003, 5 (6): 559-565. 10.1038/ncb995.
Rogulja D, Irvine KD: Regulation of cell proliferation by a morphogen gradient. Cell. 2005, 123 (3): 449-461. 10.1016/j.cell.2005.08.030.
Patel PH, Thapar N, Guo L, Martinez M, Maris J, Gau CL, Lengyel JA, Tamanoi F: Drosophila Rheb GTPase is required for cell cycle progression and cell growth. J Cell Sci. 2003, 116 (Pt 17): 3601-3610.
Herranz H, Perez L, Martin FA, Milan M: A Wingless and Notch double-repression mechanism regulates G1-S transition in the Drosophila wing. EMBO J. 2008, 27 (11): 1633-1645. 10.1038/emboj.2008.84.
Duman-Scheel M, Johnston LA, Du W: Repression of dMyc expression by Wingless promotes Rbf-induced G1 arrest in the presumptive Drosophila wing margin. Proc Natl Acad Sci USA. 2004, 101 (11): 3857-3862. 10.1073/pnas.0400526101.
Hennig KM, Neufeld TP: Inhibition of cellular growth and proliferation by dTOR overexpression in Drosophila. Genesis. 2002, 34 (1-2): 107-110. 10.1002/gene.10139.
Brand AH, Perrimon N: Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Development. 1993, 118 (2): 401-415.
Golic KG: Site-specific recombination between homologous chromosomes in Drosophila. Science. 1991, 252 (5008): 958-961. 10.1126/science.2035025.
Bellosta P, Hulf T, Balla Diop S, Usseglio F, Pradel J, Aragnol D, Gallant P: Myc interacts genetically with Tip48/Reptin and Tip49/Pontin to control growth and proliferation during Drosophila development. Proc Natl Acad Sci USA. 2005, 102 (33): 11799-11804. 10.1073/pnas.0408945102.
Conover CA, Bale LK: Insulin-like growth factor I induction of c-myc expression in bovine fibroblasts can be blocked by antecedent insulin receptor activation. Exp Cell Res. 1998, 238 (1): 122-127. 10.1006/excr.1997.3815.
Bravo R, Burckhardt J, Curran T, Muller R: Stimulation and inhibition of growth by EGF in different A431 cell clones is accompanied by the rapid induction of c-fos and c-myc proto-oncogenes. EMBO J. 1985, 4 (5): 1193-1197.
Gulati P, Gaspers LD, Dann SG, Joaquin M, Nobukuni T, Natt F, Kozma SC, Thomas AP, Thomas G: Amino acids activate mTOR complex 1 via Ca2+/CaM signaling to hVps34. Cell Metab. 2008, 7 (5): 456-465. 10.1016/j.cmet.2008.03.002.
Gao X, Pan D: TSC1 and TSC2 tumor suppressors antagonize insulin signaling in cell growth. Genes Dev. 2001, 15 (11): 1383-1392. 10.1101/gad.901101.
Ziosi M, Baena-Lopez LA, Grifoni D, Froldi F, Pession A, Garoia F, Trotta V, Bellosta P, Cavicchi S: dMyc functions downstream of Yorkie to promote the supercompetitive behavior of hippo pathway mutant cells. PLoS Genet. 2010, 6 (9):
Froldi F, Ziosi M, Garoia F, Pession A, Grzeschik NA, Bellosta P, Strand D, Richardson HE, Grifoni D: The lethal giant larvae tumour suppressor mutation requires dMyc oncoprotein to promote clonal malignancy. BMC Biol. 2010, 8 (1): 33-10.1186/1741-7007-8-33.
Tee AR, Manning BD, Roux PP, Cantley LC, Blenis J: Tuberous sclerosis complex gene products, Tuberin and Hamartin, control mTOR signaling by acting as a GTPase-activating protein complex toward Rheb. Curr Biol. 2003, 13 (15): 1259-1268. 10.1016/S0960-9822(03)00506-2.
Zhou X, Ikenoue T, Chen X, Li L, Inoki K, Guan KL: Rheb controls misfolded protein metabolism by inhibiting aggresome formation and autophagy. Proc Natl Acad Sci USA. 2009, 106 (22): 8923-8928. 10.1073/pnas.0903621106.
Leevers SJ, Weinkove D, MacDougall LK, Hafen E, Waterfield MD: The Drosophila phosphoinositide 3-kinase Dp110 promotes cell growth. EMBO J. 1996, 15 (23): 6584-6594.
Schwinkendorf D, Gallant P: The conserved Myc box 2 and Myc box 3 regions are important, but not essential, for Myc function in vivo. Gene. 2009, 436 (1-2): 90-100. 10.1016/j.gene.2009.02.009.
We thank Sheri Zola, Peter Gallant and Laura A Johnston for helpful discussion and critical reading of the manuscript. This work was supported by a Public Health Service grant from the NIH-DK085047 (PB), NIH-GM078464 (LA Johnston), and "Fondazione Cassa di Risparmio in Bologna" (FB, SR). The authors declare no conflict of interest.
FP, SR, MD, MS and NK performed the experiments. AP, DG and HS participated in the discussion and conceiving part of the experiments. ET formulated the mathematics and statistics analysis. PB participated in the discussion and conceiving the experiments and wrote the manuscript. All authors read and approved the final manuscript.