Drosophila Kruppel homolog 1 represses lipolysis through interaction with dFOXO

Transcriptional coordination is a vital process contributing to metabolic homeostasis. As one of the key nodes in the metabolic network, the forkhead transcription factor FOXO has been shown to interact with diverse transcription co-factors and integrate signals from multiple pathways to control metabolism, oxidative stress response, and cell cycle. Recently, insulin/FOXO signaling has been implicated in the regulation of insect development via the interaction with insect hormones, such as ecdysone and juvenile hormone. In this study, we identified an interaction between Drosophila FOXO (dFOXO) and the zinc finger transcription factor Kruppel homolog 1 (Kr-h1), one of the key players in juvenile hormone signaling. We found that Kr-h1 mutants show delayed larval development and altered lipid metabolism, in particular induced lipolysis upon starvation. Notably, Kr-h1 physically and genetically interacts with dFOXO in vitro and in vivo to regulate the transcriptional activation of insulin receptor (InR) and adipose lipase brummer (bmm). The transcriptional co-regulation by Kr-h1 and dFOXO may represent a broad mechanism by which Kruppel-like factors integrate with insulin signaling to maintain metabolic homeostasis and coordinate organism growth.

FOXO interacts with a number of transcription factors within the nucleus to activate or inhibit transcription of target genes 14 . The interactions between FOXO and its binding partners contribute to the transcriptional specificity of FOXO and pleotropic functions of insulin/FOXO signaling. For instance, mouse FOXO1 interacts with PGC-1α in liver to modulate insulin-mediated gluconeogenesis 15 ; mammalian FOXO1 binds to Smad2/3 in response to TGF-beta signaling and regulates cell proliferation 16 ; mammalian FOXO transcription factors (FOXO3A and FOXO4) interacts with beta-Catenin of Wnt signaling to modulate cellular oxidative response 17 . In Drosophila, dFOXO interacts with bZIP transcription factor REPTOR of mechanistic target of rapamycin (mTOR) signaling to regulate growth and energy homeostasis 18 . Interestingly, recent studies found that FOXO interacts with Ultraspiracle (Usp), a co-factor of the ecdysone receptor, to regulate ecdysone biosynthesis and developmental timing in Drosophila 19 . To date, factors of JH signaling have not been identified to directly interact with FOXO.
Kruppel-like homolog 1 (Kr-h1) is a key regulator of insect molting and metamorphosis and a major effector in JH signaling [20][21][22] . JH strongly induces the transcription of Kr-h1 via its receptor Methoprene-tolerant (Met) [21][22][23] . During insect development, Kr-h1 functions as a transcriptional repressor on neurogenesis of mushroom body and photoreceptor maturation 24,25 . Kr-h1 belongs to Kruppel-like factors (KLFs) protein family, a group of conserved C2H2 type zinc finger transcription factors. Unlike mammalian KLFs that contain three zinc finger DNA binding domains, Drosophila Kr-h1 has eight zinc finger motifs 26 . KLFs are also closely related to transcription factor Sp1 (specificity protein 1). At least seventeen KLFs are identified in mammals 27 . Both KLFs and Sp1-like factors recognized GC-rich DNA elements or CACCC-box in the promoters of target genes 27 . While KLFs and Sp1 can function as both transcription activator and repressor, the N-terminus of KLFs contains a consensus motif PXDL(S/T) that is thought to interact with transcriptional co-repressor CtBP (C-terminal binding protein) 28,29 . Some KLFs also interact with transcriptional co-activators to enhance transcriptional activities. For instance, KLF1 is acetylated through its interaction with co-activators p300 and CREB-binding protein (CBP), which leads to elevated induction of target gene beta-globin 30 .
In this study, we identified an interaction between dFOXO and the zinc finger transcription factor Kr-h1 in Drosophila. While characterizing a Kr-h1 mutant, we found that lipolysis was elevated in fasting mutant larvae. Genetic and molecular analyses revealed that Kr-h1 physically interacts with dFOXO and represses the transcriptional activation of dFOXO target genes, such as insulin receptor (InR) and adipose triglyceride lipase (bmm or brummer). The present study suggests a mechanism by which Kruppel-like factor Kr-h1 integrates with insulin/ dFOXO signaling to control lipid metabolism and coordinate organism growth.

Kr-h1 mutants delay larval development and have reduced triglyceride.
Here we study the role of Drosophila Kruppel-like factor Kr-h1 in larval development and metabolic control using a P-element-induced hypomorphic allele Kr-h1 [7] (also known as Kr-h1[k04411]) 20,31 . The P-element insertion is located within the first intron of the Kr-h1 locus and is reported to interfere with the transcription of Kr-h1 isoforms. Kr-h1 [7] homozygous mutants are partially viable during embryonic and larval development 31 . We backcrossed this Kr-h1 [7] allele into a yw R background for seven generations. The heterozygotes of cleaned Kr-h1 [7] mutants showed prolonged developmental time to pupariation (Fig. 1A), while the homozygotes arrest at either second or third instar larval stage. Kr-h1 mRNA is largely reduced in homozygous animals based on primers for the common region of all three isoforms (Fig. 1B). Using a newly generated anti-Kr-h1 antibody, three major bands were detected in larval samples from wild-type and Kr-h1 [7] heterozygotes (Fig. 1C). Each of these bands was significantly reduced in homozygous animals (Fig. 1C). In the Drosophila genome, there are three Kr-h1 isoforms. The predict molecular weight of Kr-h1 is 85.6 kDa for the α and γ isoforms, and 91.5 kDa for the β isoform 20,31 . In Fig. 1C, top band is probably the α isoform (the most abundant isoform at larval stages) because this band corresponds to the recombinant Kr-h1α protein purified using a HaloTag Protein Expression System (Fig. 1D). The nature of two lower bands in Fig. 1C is unclear, probably reflecting cross-reactivity of our Kr-h1 antibodies with other Drosophila proteins, or proteolysis events of the Kr-h1 protein. The latter might be more likely because we observed a reduction of the intensity of all three bands in Kr-h1 [7] mutants (Fig. 1C).
Defects in metabolic regulation also occur in the developmentally delayed Kr-h1 mutants. Using Nile red staining, we found that the fat body lipid droplet was slightly enlarged in Kr-h1 [7] homozygous mutants, despite their smaller nuclear size (Fig. 1D). Upon fasting (16 hours in 1x PBS), the size of lipid droplet was reduced in both wild-type and Kr-h1 mutants. To quantify fat reserves, we performed a colorimetric assay of triglycerides (TAG), the major form of lipid storage in the fat body, using larvae at 90 hours after egg laying (AEL). Among fed animals, triglycerides were reduced 2-fold by Kr-h1 mutation, while glycogen was similar among genotypes ( Fig. 1F and G). TAG is a major stored nutrient mobilized during fasting. Fasting reduced TAG stores in both genotypes, but to a significantly greater extent in Kr-h1 [7] homozygous larvae (2.9-fold vs 1.4-fold in wild-type) (Two-way ANOVA, interaction p < 0.047). Fasting reduced stored glycogen to the same extent in both genotypes ( Fig. 1F and G).
In flies, adipose triglyceride lipase brummer (bmm) is a key lipase involved in TAG mobilization 11,32 . While transcripts of bmm were slightly up-regulated by fasting in wild-type larvae (Fig. 1H), its expression was dramatically increased in fasted Kr-h1 [7] homozygous larvae (4.3-fold vs. 1.8-fold in wild-type) (Two-way ANOVA, interaction p < 0.0196) (Fig. 1H). This result is consistent with the greater TAG mobilization in fasted Kr-h1 mutants as shown in Fig. 1F, suggesting that lipase activities might be enhanced in Kr-h1 mutants, especially upon fasting. In parallel, transcripts of fly perilipin Lsd-1 were upregulated in Kr-h1 mutants and down-regulated in both genotypes upon fasting (Fig. 1I). Perilipin proteins (PLINs) are a group of lipid droplet-associated proteins that act as protective coating factors to prevent lipid breakdown by triglyceride lipases 33,34 . The elevated Lsd-1 transcripts in Kr-h1 mutants is consistent with the enlarged lipid droplet observed in Fig. 1E. Notably, repression of Lsd-1 by fasting was significantly enhanced in Kr-h1 mutants (10.2-fold vs. 2.4-fold in wild-type) (Two-way ANOVA, interaction p < 0.0001). In addition, we observed a reduced expression of FASN1 (Fatty acid synthase 1) in Kr-h1 mutants (Fig. 1J). FASN1 is the major enzyme in condensation of acetyl-CoA and malonyl-CoA to palmitic acid and palmitoyl-CoA during fatty acid biosynthesis. Since the expression of TAG lipase bmm was  [7] heterozygous mutants delayed pupation and homozygotes arrested at early larval stages. Percentage of pupariation at different developmental time points is shown. Data are represented as mean ± SE of three trials. Student t-test (**p < 0.01, *p < 0.05). (B) Kr-h1 transcripts were significantly down-regulated in Kr-h1 [7] mutants. Primers targeting common regions among three isoforms were used in qRT-PCR. Each bar represents mean ± SE of three biological replicates. Statistical significance between wild-type and mutants is assessed by student t-test (***p < 0.001). (C) Reduced Kr-h1 protein expression in Kr-h1 [7] homozygous mutants. Larvae at 90 hr AEL (after egg laid) were used in western blots. (D) Western blot for recombinant Kr-h1α isoform (second band). The top band is uncleaved Kr-h1α-Halo fusion protein. (E) Nile red staining of fat body lipid droplet in wildtype and Kr-h1 [7] homozygous mutants. Nuclear staining is in blue. Scale bar: 20 µm. (F) Kr-h1 mutant larvae have reduced TAG level. Upon starvation, TAG mobilization was faster in Kr-h1 mutants than in wild-type larvae. Larvae at 90 hr AEL were fasted for 16 hr in culture vials with wet kimwipe soaked with PBS. Each bar represents mean ± SE of three biological replicates. Statistical significance is assessed by two-way ANOVA followed by Tukey multiple comparisons test (***p < 0.001, **p < 0.01, *p < 0.05). (G) Glycogen contents and the utilization rate were not affected by Kr-h1 mutation. (H) Transcripts of TAG lipase brummer (bmm) were up-regulated by fasting and Kr-h1 mutation. The fasting-induced bmm expression was further enhanced by Kr-h1 mutation. (I) Fasting-triggered fly perilipin Lsd-1 repression was significantly enhanced in Kr-h1 mutants. Statistical significance is assessed by two-way ANOVA followed by Tukey multiple comparisons test (***p < 0.001, **p < 0.01, *p < 0.05). (J) The expression of FASN1 (Fatty acid synthase 1) was reduced in Kr-h1 mutants. Student t-test (*p < 0.05). not altered in fed Kr-h1 mutants (Fig. 1H), it is likely that the reduction of TAG in Kr-h1 mutants is through the regulation of FASN1 expression and lipogenesis, rather than lipolysis.
Collectively, these results suggest that Kr-h1 plays an important role in lipid metabolism. Upon fasting, Kr-h1 regulations lipolysis through the transcriptional regulation of triglyceride lipase bmm, and lipid droplet-associated protein Lsd-1. In contrast, Kr-h1 acts on lipogenesis under normal fed condition. The present study focuses on the role of Kr-h1 in lipolysis under nutrition deprivation.
Kr-h1 mutants have reduced insulin signaling. One way Kr-h1 might modulate TAG is through interactions with insulin/IGF signaling. Insulin/IGF signaling is a metabolic master regulator that controls lipase gene expression through its downstream transcription factor dFOXO 35 . Here we see that phosphorylation of IIS-regulated kinase AKT was reduced in Kr-h1 [7] homozygotes ( Fig. 2A). Furthermore, Kr-h1 [7] homozygotes had reduced expression of two insulin-like peptides (dilp2 and dilp5), which are the major DILPs produced from brain neurosecretory cells, known as insulin producing cells (IPCs) ( Fig. 2B and C).
Reduced insulin signaling is expected to activate forkhead transcription factor dFOXO 6 . Accordingly, mRNA expression of two key dFOXO target genes, 4ebp (eukaryotic translation initiation factor 4E binding protein) and InR were significantly induced in Kr-h1 mutants ( Fig. 2D and E), and InR expression was further increased in fasted Kr-h1 [7] homozygotes (5.2-fold vs. 3.4-fold in wild-type) (Two-way ANOVA, interaction p = 0.1023) (Fig. 2F). It is known that the transcription of InR is controlled by dFOXO through a negative feedback mechanism 36 . Thus, in Kr-h1 mutant larvae, insulin signaling is inhibited and dFOXO is activated. Kr-h1 may produce these effects in two ways: Kr-h1 regulates the expression of dFOXO target genes via transcriptional co-regulation and direct interaction with dFOXO, or through indirect modulation of insulin/AKT signaling and the expression of dilp2 and dilp5.

Kr-h1 genetically interacts with dfoxo to regulate the transcription of InR and bmm, and lipid metabolism.
To determine if dFOXO is required for Kr-h1 to mediated lipid metabolism, we generated a double mutant by combining Kr-h1 [7] and dfoxo [21] 37 . Interestingly, we found that dfoxo [21] heterozygotes rescued the delayed pupariation of Kr-h1 [7] heterozygous mutants (Fig. 2G). As well, dfoxo [21] mutants suppressed the elevated InR and bmm expression found in Kr-h1 [7] mutants ( Fig. 2H and I), confirming that these transcription factors co-regulate key metabolic genes. Furthermore, the reduction of TAG in Kr-h1 [7] homozygous mutants was rescued by dfoxo [21]−/− (Fig. 2J). Together, these results reveal a genetic interaction between Kr-h1 and dFOXO in the control of the transcription of metabolic genes and lipid metabolism.

Kr-h1 physically interacts with dFOXO.
Kr-h1 and dFOXO may interact directly or indirectly to regulate the expression of InR and bmm. To test the possibility of direct interaction, we co-immunoprecipitated (Co-IP) Kr-h1 and dFOXO in cultured Drosophila cells. We were able to pull down endogenous dFOXO from nuclear and cytoplasmic extracts using an anti-dFOXO antibody. Interestingly, Kr-h1 was detected in the protein complex from the nuclear extracts, but not from the cytoplasmic extracts ( Fig. 3A), suggesting that Kr-h1 can form a protein complex with dFOXO in the nuclei.
To identify the protein interaction site between these transcriptional factors, we cloned a series of deletion fragments that contained different protein domains into the Gateway expression vectors. Both the DNA binding domain and transactivation domain of dFOXO bound to full-length Kr-h1 proteins (Fig. 3B). On the other hand, the Kr-h1 fragments that contain transactivation/repression domain (a Q-rich domain) bound to full-length dFOXO proteins, while Kr-h1 fragments with no Q-rich domain showed no binding (Fig. 3C). Therefore, the transactivation/repression domain of Kr-h1 is responsible for the interaction between Kr-h1 and dFOXO.
The direct interaction between dFOXO and Kr-h1 may provide a mechanism for the transcriptional repression of dFOXO target genes by Kr-h1. To test this idea, we co-expressed dFOXO with Kr-h1 (or Kr-h1 Q-rich domain) in Kc167 cells. The mRNA expression of bmm was significantly induced by dFOXO alone, and this induction was blocked by co-expressing either full-length of Kr-h1 or Q-rich domain (Fig. 3D). Thus, Kr-h1 appears to repress dFOXO transcriptional activity through direct protein-protein interactions.

Kr-h1 binds to the promoters of insulin receptor and brummer lipase adjacent to dFOXO binding sites.
Kr-h1 and dFOXO physically interact and may thus transcriptionally co-regulate metabolic genes. It has been previously shown that dFOXO binds to the promoter regions near transcriptional start sites of InR and bmm 35,36 , although our recent ChIP-Seq analysis (unpublished) suggests that dFOXO also strongly bound the promoter region near the 5′-UTR of InR (P1 region as shown in Fig. 4A) that contains a canonical FOXO binding motif (GTAAATAA). To identify potential Kr-h1 response elements of InR and bmm, we searched their promoters using mammalian KLF motifs in the Jaspar database (http://jaspar.genereg.net). Three putative KLF binding sites denoted P1~P3 in each promoter were identified including sites in 5′-UTR and intronic regions ( Fig. 4A and B) (Supplementary Table S1). We did not find any sites corresponding to the Bombyx Kr-h1 response element (GACCTACGCTAACGCTAAATAGAGTTCCGA) reported by Kayukawa et al. 23 .
Binding of dFOXO and Kr-h1 to these putative sites was determined by ChIP-PCR analysis in fasted animals. At InR, dFOXO binding was strongest in the P1 region located at the 5′-UTR region (Fig. 4C), while Kr-h1 bound most strongly to the P3 regions, about 15k bp away from P1 region (Fig. 4D). At bmm lipase, both dFOXO and Kr-h1 bound with highest affinity in the P1 region ( Fig. 4E and F). The co-localization of Kr-h1 and dFOXO binding suggests these factors could interact at promoters to control the transcriptional activation of the key metabolic genes (e.g., bmm lipase).
Kr-h1 represses dFOXO binding to the promoter of InR and bmm. Kr-h1 may repress dFOXO activity by inhibiting its binding at response elements in bmm and InR. We performed a ChIP-PCR to test this possibility using anti-dFOXO antibody and Kr-h1 [7] mutants. dFOXO binding to the InR P1 region was increased from SCIEnTIfIC RepoRts | 7: 16369 | DOI:10.1038/s41598-017-16638-1 2.9-fold relative to negative control (Act5C) in fasted wild-type to 8.95-fold in fasted Kr-h1 mutants (Two-way ANOVA, interaction p < 0.0001) (Fig. 4G). In contrast, dFOXO binding to the bmm P1 region was slightly but
non-significantly increased from 2.1-fold in fasted wild-type to 2.8-fold in fasted Kr-h1 mutants (Two-way ANOVA, interaction p = 0.5862) (Fig. 4H). At the InR promoter in particular, inhibition of dFOXO-DNA interaction may be one mechanism by which Kr-h1 modulates dFOXO transcriptional activity. Notably, in a reciprocal experiment with anti-Kr-h1 antibody, the binding of Kr-h1 to InR and bmm promoters was abolished in dfoxo [21] mutants ( Fig. 4I and J). These data suggest that Kr-h1 may be recruited after dFOXO binds to the promoters of target genes, and Kr-h1 subsequently modulates the transcriptional activities of dFOXO through interfering with dFOXO-DNA interactions.
Kr-h1 expresses in adipose tissue to control larval development and lipid metabolism. To determine where Kr-h1 and dFOXO interact in vivo, we examined the nuclear co-localization of Kr-h1 and dFOXO in larval fat body using our anti-Kr-h1 antibodies. Previous studies showed that Kr-h1 expressed broadly in many tissues during late larval and prepupal stages, including imaginal discs, trachea, central nervous system, ring gland, salivary gland, muscle, gut, and fat body 38 , and the expression was primarily detected in nuclei. We find that in early L3 larval development, Kr-h1 is expressed at low level and primarily in cytosol (Fig. 5A). We also observed cytoplasmic expression of Kr-h1 in Kc167 cells under normal culture condition (Supplementary figure S6). Interestingly, nuclear co-localization of dFOXO and Kr-h1 was increased in fasted larval fat body (90 hr AEL) (Fig. 5A). We verified the specificity of our Kr-h1 antibody and found no nuclear expression of Kr-h1 in Kr-h1 [7] homozygous mutants (Fig. 5B). Thus, dFOXO and Kr-h1 may interact in fat body nuclei to co-regulate the transcriptional activation of target genes.
dFOXO expressed in fat body regulates local lipid metabolism and the production of Dilp2 from brain IPCs 39 . To determine where Kr-h1 acts to regulate lipid metabolism and larval development, we knocked down Kr-h1 message through RNA interference (RNAi) with tissue-specific Gal4 drivers. Knockdown of Kr-h1 in fat body (r4-gal4) and muscle (Mhc-gal4) delayed pupariation, while knockdown in gut, IPCs and Corpus allatum  [21] mutants. (J) Kr-h1 binding to bmm promoter is abolished in fasted dfoxo [21] mutants. Each bar represents mean ± SE of three biological replicates. Statistical significance is assessed by one-way ANOVA, followed by Dunnett's multiple comparisons (*p < 0.05, ns: not significant).
SCIEnTIfIC RepoRts | 7: 16369 | DOI:10.1038/s41598-017-16638-1 (CA) showed no effects on larval development (Fig. 6A). Since fat body is the major site for triglyceride storage in Drosophila, we further examined the role of Kr-h1 in the regulation of lipid metabolism in this tissue. Fat body-specific knockdown of Kr-h1 induced bmm transcription, while overexpression of Kr-h1 repressed it (Fig. 6B). Kr-h1 expressed in fat body also increased whole larval TAG levels (Fig. 6C). Thus, adipose-expressed Kr-h1 is essential for larval development and metabolic regulation.

Juvenile hormone signaling regulates lipase bmm through dFOXO. The interaction between Kr-h1
and dFOXO has the potential to integrate development and nutrient signaling. Nutrient signaling through FOXO involves insulin, AMPK, SIRT and JNK in both insect and mammals alike 6,14 . On the other hand, the upstream regulators of Kruppel-like factors are poorly characterized in vertebrates, but among insects Kr-h1 is decisively regulated by JH, a key hormonal signal involved in molting and metamorphosis 21 . In particular, JH induces the transcription of Kr-h1 via the JH receptor Methoprene-tolerant (Met) 22,40 . In this capacity, recent studies suggest that JH and Met are involved in not only development programming, but also in metabolic control 9,41-43 , although how JH affects metabolism is largely unknown.
Given that Kr-h1 and dFOXO functionally interact to control lipid metabolism, we examined if this feature provides a way for JH to affect metabolic regulation through bmm transcription. Consistent with previous studies 9 , triglyceride levels were reduced in flies where the corpora allata were genetically ablated (CAX) (Fig. 7A). Conversely, wild-type flies exposed to the JH analog (JHA) methoprene had elevated TAG contents compared to controls (Fig. 7B). Additionally, Met mutations down-regulated TAG levels (Fig. 7C) and up-regulated bmm mRNA (Fig. 7D). Met also genetically interacts with dFOXO to regulate the mRNA expression of bmm (Fig. 7D). JH may therefore regulate bmm through a Met-mediated interaction with dFOXO. Supporting this prediction, methoprene treatment inhibited the expression of bmm in wild-type flies, but not in dfoxo [21] mutants (Fig. 7E). Furthermore, fasting reduced JH titers about 2-fold in both female and male flies (Fig. 7F). But while it is known that JH positively regulates Kr-h1 transcription 21,23 , Kr-h1 mRNA was not reduced upon fasting (Fig. 7G). On the other hand, methoprene treatment was sufficient to induce Kr-h1 transcription (Fig. 7H). Overall these results indicate that JH signaling interacts with dFOXO to regulate lipid metabolism and lipase gene expression, but the mechanistic role of Kr-h1 in this process remains to be elucidated.

Discussion
Transcriptional coordination is a key process contributing to metabolic homeostasis 44 . Multiple transcription factors interact at their genomic binding sites to enhance transcriptional specificity and pleiotropic functions of metabolic pathways. As a key node in the metabolic network, forkhead transcription factor FOXO has been shown to interact with diverse transcription co-factors and thereby integrate signals to control metabolism and oxidative stress 6,14 . Intriguingly, in recent genomic studies [45][46][47] , enriched FOXO binding at specific genes does not always correlate to elevated transcriptional output, suggesting there exists inhibitory or inertial mechanisms to repress the transcriptional activity of DNA-bound FOXO.
Here we find that Drosophila Kruppel-like factor Kr-h1 acts as a repressor of dFOXO to modulate induction of two dFOXO target genes, InR and bmm. Like other FOXO interacting partners, Kr-h1 physically binds to dFOXO and inhibits the expression of dFOXO targets by influencing the binding affinity of dFOXO to DNA. The transcriptional activity of FOXO is typically regulated in two layers. The first and probably most important regulation is through PTM, including phosphorylation, acetylation and ubiquitination 6 . PTM of FOXO proteins can affect its subcellular localization (by phosphorylation), DNA binding affinity (by acetylation) and protein degradation (by ubiquitination). Interestingly, the effects of acetylation on FOXO factors seem to be quite different from those by acetylation on KLFs. Acetylation of FOXO by co-factor CBP/p300 weakens the FOXO binding to its DNA targets 48 , while CBP/p300 acetylated KLF1 shows increased transcriptional activation of target gene beta-globin 30 . The second mechanism for the regulation of FOXO activity is through the interaction between FOXO and other transcription factors or co-factors. FOXO factors have been shown to interact with diverse transcription factors (e.g. Smad3/4, PGC-1α, STAT3) that often potentiate the expression of FOXO target genes 14 . Kr-h1 identified in our study presents another example for this type of modulatory regulation, although the interaction between Kr-h1 and dFOXO results in transcription repression, instead of activation. While we do not fully resolve how Kr-h1 blocks dFOXO activity, it seems that Kr-h1 can inhibit dFOXO binding to its DNA targets. This result is similar to previous studies showing reduced FOXO-DNA binding upon interaction with androgen receptor (AR) 49 and with peroxisome proliferator-activated receptor-γ (PPARγ) 50 . Alternatively, Kr-h1 may act by inhibiting recruitment of dFOXO-coactivators (e.g. SIRT or CBP/p300) or by sequestering these coactivators away from dFOXO. Kr-h1 may also recruit additional co-repressors (e.g. CtBP or Sin3-HDAC) to the dFOXO transactivation sites to block the transcriptional activation of target genes. The N-terminal Q-rich domain of KLFs is crucial for the recruitment of co-repressors CtBP and Sin3A 51 . In our co-immunoprecipitation assays, the Q-rich domain of Kr-h1 strongly binds to dFOXO, suggesting Kr-h1 might inhibit dFOXO activity through recruiting co-repressors. One possible candidate is an Sds3-like gene (CG14220), which was previously found to co-immunoprecipitate with Kr-h1 52 . Sds3-like gene family proteins can form Each bar represents mean ± SE of three biological replicates. Student t-test (*p < 0.05). (B) Flies exposed to JH analog (JHA) methoprene show increased TAG levels. Each bar represents mean ± SE of three biological replicates. One-way ANOVA (*p < 0.05). (C) Met mutants have reduced TAG levels. Student t-test (*p < 0.05). (D) Genetic interaction between Met and dfoxo in the regulation of bmm transcripts. bmm transcription is up-regulated in Met mutants, which was rescued by dfoxo 21 mutants. Oneway ANOVA (*p < 0.05, ns: not significant). (E) JH analog (JHA) methoprene treatment led to reduced bmm expression in wild-type female flies, but not in dfoxo 21 mutant flies. Each bar represents mean ± SE of three biological replicates. Statistical significance is assessed by two-way ANOVA (*p < 0.05, ns: not significant) (F) JH titer is decreased upon fasting. 10-day-old adult flies were fasted (in culture vial with wet kimwipe soaked with PBS) for 16 hours before collected for JH quantification. Each bar represents mean ± SE of 5~7 biological replicates. Statistical significance is assessed by student t-test (*p < 0.05). (G) The mRNA expression of Kr-h1 did not change upon fasting. (H) Methoprene treatment induced Kr-h1 transcription. Each bar represents mean ± SE of three biological replicates. One-way ANOVA (*p < 0.05, ns: not significant).
co-repressor complex with Sin3A and HDAC to inhibit gene transcription via interactions with sequence-specific transcription factors 53 .
KLFs have well documented roles in cell proliferation, differentiation and apoptosis 51 . The function of Kr-h1 in lipid metabolism and insulin signaling identified in the present study complements recent studies where KLFs function in cellular metabolic regulation, such as gluconeogenesis 54,55 . Likewise, KLF15 deletion in mice produces hypoglycemia and impaired amino acid catabolism upon fasting 55 . Additionally, KLF5 heterozygous mice are resistant to high-fat diet-induced obesity. SUMOylation modulates the transcriptional activities of KLF5 and its association with peroxisome proliferator-activated receptor-delta (PPAR-delta) to control the expression of carnitine-palmitoyl transferase-1b (Cpt1b), uncoupling proteins 2 and 3 56 . Interestingly, KLF4 has recently been identified as a direct target gene of FOXO-mediated transcription during B cell development 57 , suggesting a potential interaction between KLF and FOXO transcriptional regulatory network.
Our ChIP-PCR studies suggest that Drosophila KLF Kr-h1 transcriptionally controls many metabolic genes, including key dFOXO targets (e.g. InR and triglyceride lipase bmm). While it is not known whether Drosophila Kr-h1 could broadly interact with dFOXO across the genome, such a genome-wide interaction between mammalian FOXO factors and KLFs has been suggested by a recent meta-analysis 58 . Because Kr-h1 plays an important role in morphogenesis during Drosophila development 21 , the interplay between Kr-h1 and dFOXO raises the possibility that Kr-h1 coordinates growth and development through insulin/dFOXO-mediated metabolic regulation.
We also noticed that the direct interaction between Kr-h1 and dFOXO occurs strongly upon fasting. In the normal fed state, transcriptional regulation of lipase bmm is weakly regulated by Kr-h1, although the TAG is significantly reduced in fed Kr-h1 [7] mutants. We speculate that decreased TAG in Kr-h1 [7] mutants is probably due to the reduction in lipogenesis and fatty acid synthesis, rather than to increased lipolysis. Our preliminary data reveal that the expression of fatty acid synthase is down-regulated in Kr-h1 [7] mutants. The potential connection between Kr-h1 and lipogenesis needs to be further examined.
In insects, Kr-h1 is one of the key effectors of JH signaling, an important hormonal pathway governing insect molting, metamorphosis and reproduction. Recent studies reveal that JH also participates in the regulation of carbohydrate and lipid metabolism 9,41-43,59,60 . In Tsetse flies, JH and insulin co-regulate the expression of TAG lipase and inhibit lipolysis 9 , which is consistent with our observation that TAG metabolism and lipase brummer expression are regulated by JH and its receptor Met. Similarly, JH III treatment reduced FOXO levels in diapause-destined female mosquitoes, Culex pipiens 13 . Although we did not test the paralog of Met, germ cell-expressed (GCE) [61][62][63] , previous studies have demonstrated GCE is not involved in lipid metabolism during lactation state of Tsetse flies 9 . Finally, JH signaling was recently found to genetically interact with insulin/dFOXO to control larval growth rate and define final body size in Drosophila 10 . Thus, the transcriptional co-regulation of lipid metabolism by Drosophila Kr-h1 and dFOXO may contribute to a novel mechanism through which JH interacts with insulin signaling to integrate metabolism and growth during larval development.
Since both Kr-h1 and dFOXO express highly in metabolic tissues (fat body and muscle) of Drosophila, these transcription factors are likely to co-regulate many key metabolic genes. As well, these metabolic tissues contribute significantly to other insect physiology and organismal functions, such as stress resistance and aging that are tightly regulated by insulin/dFOXO signaling 64 and by JH signaling 41 . We predict that the interplay between Kr-h1 and dFOXO will also contribute to the regulation of these adult physiological processes. In fact, a recent study reported that loss of kruppel-like factors (KLFs) abolished the lifespan extension in multiple longevity paradigms in C. elegans, such as dietary restriction and insulin receptor mutant daf-2 65 . Thus, identifying the key co-factors, mechanisms, and downstream events of the Kr-h1/dFOXO transcriptional network may advance our understanding of how JH and insulin signaling coordinately regulate metabolism, development, and aging.

Materials and Methods
Fly Husbandry and Stocks. Flies were maintained at 25 °C, 40% relative humidity and 12-hour light/dark. Adults were reared on agar-based diet with 0.8% cornmeal, 10% sugar, and 2.5% yeast (unless otherwise noted). Fly stocks used in the present study are: Kr-h1 [7] 72 , S106-GS-gal4 73 . Double mutants were made by crossing Kr-h1 [7] or Met [27] to foxo [21] respectively. Corpus allatum (CA) ablation flies (named CAX flies) are generated in our laboratory as previously described 41 . yw R flies were used as wild-type flies in most of the experiments. For methoprene treatment, adult flies were exposed for 24~48 hours to various concentrations of methoprene applied to the side of culture vials.
Kr-h1 Antibody, Kr-h1 recombinant protein synthesis and Western Blot. Kr-h1 polyclonal antibody was generated in rabbits against the short peptide sequence 'LIEHFKRGDLARHG' (Covance, Dedham, MA, USA) and affinity purified (Thermo Fisher Scientific, Waltham, MA, USA). The specificity of Kr-h1 antibody was verified by western blots (Fig. 1) and immunostaining (Fig. 5) using Kr-h1 [7] mutants. Recombinant Kr-h1α protein was produced and purified using the HaloTag Protein Expression System (Promega, Madison, WI, USA). Kr-h1α cDNA clone (LD32311) was obtained from DGRC (Drosophila Genomics Resource Center) and cloned into pFN29K His 6  Quantitative RT-PCR. Total RNA was extracted using Trizol reagent (Thermo Fisher Scientific, Waltham, MA, USA) from 10~15 synchronously staged larvae or whole adult flies. DNase-treated total RNA was quantified and about 500 ng of total RNA was reverse transcribed to cDNA using iScript cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA). QPCR was performed with an ABI prism 7300 Sequence Detection System (Thermo Fisher Scientific, Waltham, MA, USA). Three to five biological replicates were used for each experimental treatment. mRNA abundance of each gene was normalized to the expression of ribosomal protein L32 (RpL32 or rp49) by the method of comparative C T . Primer sequences are listed in Supplementary Table S2.
Pupariation timing analysis. Synchronized eggs were placed on 35 × 10 mm petri dishes containing standard medium (see above) at 20~30 eggs per dish. The numbers of pupae were recorded 2~3 times every day around 120 hr AEL till all larvae molt into pupae.
Nile red staining. Larval fat body were dissected in 1x PBS and stained with 0.00005% Nile Red (Sigma-Aldrich, St Louis, MO, USA) and 1 µg/ml Hoechst 33342 (ImmunoChemistry Technologies, Bloomington, MN, USA) in 75% glycerol for 10 min at room temperature. Fat body was then mounted in 75% glycerol and imaged using an Olympus BX51WI upright epifluorescence microscope.
Metabolic assays. All metabolic analyses were performed as previously described 73,74 . For TAG assay, 25 staged larvae or six adult flies were collected and homogenized in 1xPBS containing 0.1% Tween 20 and TAG was quantified using Thermo Scientific ™ Triglycerides Reagent (Thermo Fisher Scientific, Waltham, MA, USA).
For glycogen measurement, samples were digested with amyloglucosidase (Sigma-Aldrich, St Louis, MO, USA) and glucose contents were quantified using Thermo Scientific ™ Glucose Hexokinase Reagents (Thermo Fisher Scientific, Waltham, MA, USA). The relative level of each metabolite was obtained by normalizing the metabolites to total protein.
Immunoprecipitation and pull-down. All the immunoprecipitation and pull-down experiments were conducted in Drosophila Kc167 cells adapted to serum-free culture medium (Drosophila Schneider Medium). Either full-length (Kr-h1 α-isoform) or partial gene products were cloned into Drosophila Gateway Vectors with N-terminal tags (FLAG and HA) following Drosophila Gateway Vectors protocols (https://emb.carnegiescience. edu/Drosophila-gateway-vector-collection). About 1 µg of constructs were transfected to 2 × 10 6 Kc167 cells using Effectene reagent (Qiagen, Hilden, Germany). Two days after transfection, cells were harvested and lysed in NP-40 lysis buffer (Thermo Fisher Scientific, Waltham, MA, USA) with proteinase inhibitors (Sigma-Aldrich, St Louis, MO, USA). To pull-down target proteins, total protein extracts were incubated with proper antibodies and Dynabeads Protein A (Thermo Fisher Scientific, Waltham, MA, USA). Following pull-down, western blotting was performed to examine protein complex. Antibodies used in pull-down and western blots include rabbit anti-Kr-h1 and anti-dFOXO produced in our laboratory, rabbit anti-HA (Covance, Dedham, MA, USA), and mouse anti-FLAG (Sigma-Aldrich, St Louis, MO, USA). Nuclear extracts for immunoprecipitation were conducted with a nuclear extraction kit (Active motif, Carlsbad, CA, USA).

Immunostaining and imaging.
To examine the co-localization of Kr-h1 and dFOXO, larval fat body was dissected from fed or fasted 3 rd instar larvae (90 hr AEL) (For fasting, larvae were placed onto wet kimwipe soaked with 1 × PBS for 16 hours). Tissue immunostaining were performed as previously described 47  Chromatin immunoprecipitation (ChIP). ChIP was conducted as previously described 47 . About 50 staged larvae were used in each sample. Flies were homogenized and cross-linked in 1xPBS containing 1% formaldehyde. The fly nuclear extractions were sonicated using a Branson 450 sonicator to break down the chromatins. Immunoprecipitation was performed using Dynabeads Protein A and anti-Kr-h1 and anti-dFOXO antibodies. Following the wash with LiCl and TE buffer, the DNA-protein complex was eluted, reverse cross-linked, digested with Proteinase K and RNase. Kr-h1-bound or dFOXO-bound DNA fragments were purified and used as templates in qPCR analysis. Binding enrichment was calculated as the fold change between ChIP DNA vs. input DNA (Chromatin extracts before immunoprecipitation). The binding to the coding region of Actin (Act5C) was used as negative controls.
Juvenile hormone quantification. For each sample, 197-200 individual flies (7~10-day-old) were placed in 500 µl hexane in a glass vial with a Teflon cap insert and stored at −80 °C prior to analysis. To extract the hormone, the flies were crushed with a Teflon tissue grinder. The resultant homogenate was centrifuged at 3500 rpm for 5 min, and the supernatant was removed to clean vial. Extraction was conducted three times, combining the resultant supernatant from each sample. The gas chromatography/mass spectrometry (GC-MS) method 76 , as modified 77,78 , was used to quantify juvenile hormone (JH). Samples were eluted through aluminum oxide columns successively with hexane, 10% ethyl ether-hexane and 30% ethyl ether-hexane. Samples were subjected to a second series of aluminum oxide elutions (30% ethyl ether-hexane then 50% ethyl-acetate-hexane) after derivatization with methyl-d alcohol (Sigma-Aldrich, St Louis, MO, USA) and trifluoroacetic acid (Sigma-Aldrich, St Louis, MO, USA). Purified samples were analyzed on an HP 7890A Series GC (Agilent Technologies, Santa Clara, CA, USA) equipped with a 30 m × 0.25 mm Zebron ZB-WAX column (Phenomenex, Torrence, CA, USA) and coupled to an HP 5975C inert mass selective detector with helium as the carrier gas. MS analysis occurred in the SIM mode, monitoring at m/z 76 and 225 to ensure specificity for the d3-methoxyhydrin derivative of JH III. Total abundance was quantified against a standard curve of derivatized JH III and using farnesol (Sigma-Aldrich, St Louis, MO, USA) as an internal standard. The detection limit is approximately 1 pg. Statistical analysis. GraphPad Prism 6 (GraphPad Software, La Jolla, CA) was used for statistical analysis.
To compare the mean value of treatment groups versus that of control, either student t-test or one-way ANOVA was performed using Dunnett's test for multiple comparison. The effects of mutants on starvation responses was analyzed by two-way ANOVA, including Tukey multiple comparisons test.