AbstractInsect development is primarily controlled by juvenile hormone (JH) and 20-hydroxyecdysone (20E), which regulate gene cascades leading to changes in phenotype, physiology, and behavior. Besides these hormones, microRNAs play a crucial role in insect development by regulating gene expression at the post-transcriptional level. To advance the molecular understanding of holometabolous developmental events, we investigate the pupal phase in the honeybee, Apis mellifera. In this study, we assessed the expression profiles of genes components of JH and 20E cascades – Usp, ftz-f1, EcR, Met, Chd64, InR-2, Kr-h1 and Tai – as well as the microRNAs miRNA-34 and miRNA-281 during pupal development of A. mellifera. We then analyzed the impact of JH and 20E treatments on the expression of these developmental genes and their putative regulators, the microRNAs. Overall, the selected genes and miRNAs remained stable or were downregulated following 20E treatment, while treatments with JH, upregulated most of our candidate developmental genes and microRNAs. Notably, the expression profile of Met, an intracellular receptor of JH, showed a strong correlation with fluctuations in 20E titers during pupal development. Furthermore, a computational analysis, followed by experimental assays, points to both miR-34 and miR-281 as potential regulators of pupal development in A. mellifera. This study paves the way for a better understanding of how JH and 20E hormones interact with developmental genes and microRNAs (miR-34 and miR-281) to regulate pupal development in honeybees, elucidating a piece of this complex network of interactions.
IntroductionIn insects, development is characterized by growth periods intercalated by molts, transitioning from immature stages to adults1,2,3. Holometabolous insects, to which our study model belongs, develop from embryo to adult, passing through larval, pre-pupal, and pupal phases. These developmental stages are achieved through a series of molts, primarily regulated by the levels of juvenile hormone (JH) and 20-hydroxyecdysone (20E)4,5.The hormones JH and 20E regulate insect development by exerting opposing effects at specific developmental stages; while JH prevents premature molting, 20E stimulates it. Together, they regulate the precise timing of metamorphosis6 by activating gene cascades that result in phenotypic, physiological, and behavioral changes1,4,5.Among the genes whose expression is triggered by JH and 20E, ultraspiracle (Usp), ecdysone receptor (EcR), Methoprene-tolerant (Met), kruppel homolog 1 (Kr-h1), calponin-like (Chd64), insulin-like receptor-like (InR-2), nuclear hormone receptor FTZ-F1 (ftz-f1) and Taiman (Tai) stand out as key regulators of insect development1,7,8,9,10,11,12. USP and ECR are nuclear receptors for 20E, promoting the transcription of 20E target genes13. Additionally, USP also plays as JH receptor14,15; in Drosophila melanogaster and Apis mellifera it was characterized as component of the JH transcription factor complex16. MET acts as an intracellular receptor for JH, and upon binding to the hormone, it translocates to the nucleus where it associates with transcription factors, such TAI and USP. This transcriptional complex then binds to JH response elements in the promoter regions of JH target genes, regulating their expression17,18,19,20,21. Downstream of MET and modulated by JH titers, the zinc finger transcription factor Kr-h1 inhibits precocious metamorphosis and adult morphogenesis1,22. Playing in the crosstalk between JH and 20E with role in both hormone pathways, is the calponin-like protein, identified as a chaperone, CHD648. When phosphorylated by PKC, triggered by 20E, CHD64 translocates to the nucleus and participates in the 20E signaling pathway. In the absence of phosphorylation, induced by JH, CHD64 binds to USP in the nucleus, integrating into the JH transcription factor complex, mediating JH signal transduction8. The insulin receptor 2 activates the expression of 3-hydroxy-3-methyl-glutaryl-CoA reductase, an enzyme that participates in the JH biosynthesis process7,23. Finally, the transcription factor FTZ-F1 has been shown to be essential for corpora allata activation, participating in JH biosynthesis in Blattella germanica24. Moreover, it acts as a competence factor at the end of each molt, as a consequence of the 20E peak25, allowing expression of downstream genes involved in developmental events (e.g., cuticular genes in pharate-adult stages).In addition to protein-coding genes, microRNAs (miRNAs) are regulatory non-coding RNAs known to modulate biological processes, including metamorphosis26,27,28,29. In A. mellifera, 197 miRNAs are expressed in the wing discs during metamorphosis (L5-PP2 stage), including miR-281 and miR-3430. The miR-281 regulates EcR levels during the development of Malpighian tubules in Bombyx mori and in a cross-species interaction between the endoparasitoid Cotesia vestalis and the moth Plutella xylostella. In this latter case, miR-281 is synthesized by C. vestalis and modulates EcR expression, thereby inhibiting the growth of its host, P. xylostella. This highlights the conserved role of miR-281 in targeting EcR transcripts across species31,32. In honey bees, the expression levels of miR-34 are perturbed by the knockdown of EcR33.In this study, we examined the expression profile of the genes Usp, EcR, Met, Kr-h1, Chd64, InR-2, ftz-f1, and Tai, and how their expression is affected by alterations in JH and 20E levels. We also investigated whether miR-34 and miR-281 can bind to the predicted miRNA binding sites on the 3′ UTR of these developmental genes, and assessed their responsiveness to JH and 20E. Our results show that during the pupal development Met responds to 20E treatment, while expression of Usp is affected by both JH and 20E treatments. miR-34 binds to the 3′ UTR of InR-2, Chd64, Kr-h1, and ftz-f1, while miR-281 targets the 3′ UTR of EcR and ftz-f1. Both miR-34 and miR-281 are upregulated in response to JH treatment and downregulated in response to 20E treatment. Altogether, our study helps to unveil the complex network of interactions between genes, hormones, and miRNAs that orchestrate the honeybee pupal development.Material and methodsSampling and age controlHoneybees (A. mellifera) larvae and pupae were obtained from the colonies kept in the experimental apiary of the Department of Genetics at the University of São Paulo (Ribeirão Preto, Brazil). To control the age of the samples, a fertilized queen was caged in an empty area of the frame for oviposition, for a period of 6 h. To avoid the effects of circadian rhythms on the gene expression of target genes, we implemented a schedule to cage the queen at different times of the day, ensuring that each specific stage could be collected at 10 AM (Supplementary Figure S1).Gene expression profile of the candidate genes during preimaginal developmentWe collected worker bees at each of the following stages: L5S1, L5S2, L5S3, PP1, PP2, PP3, Pw, Pp, Pdp, Pb, Pbl, Pbm, Pbd and Ne (Supplementary Figure S1). The bees were staged according to the classification by Michelette and Soares34. To ensure genetic variability, each stage consisted of three worker bees from distinct colonies. The heads were removed from the samples to exclude the primary sites of morphogenetic hormone production, specifically the corpora allata and corpora cardiaca glands. The remaining body was used to extract total RNA with TRIzol reagent, according to the manufacturer’s instructions. For each sample, 3 µg of total RNA were incubated with 0.5 unit of DNase I Amplification Grade (Invitrogen) and then proceeded to reverse transcription reactions using 100 units SuperScript II® (Invitrogen), and 0.5 µg/µL oligo(dT)12–18 (Invitrogen). The resulting cDNA was used as template for qPCR reactions using 10 µL of PowerUp™ Sybr™ Green Master Mix (Applied Biosystems), 0.8 µL of each primer (forward and reverse—10 µM), 6.4 µL of water and 2 µL of cDNA (1:10 dilution). Each biological replicate was processed in three technical replicates. The gene sequences were retrieved from the National Center for Biotechnology Information (NCBI) database (version Amel_HAv3.135,36). Gene annotation was performed using the Artemis platform37 to precisely define exon boundaries for primer design (Supplementary Figure S2). The sequences and fragment lengths are listed in Table 1. Amplification was carried out in a StepOnePlus™ Real-Time PCR System (Applied Biosystems) following the conditions: 50 °C for 2 min (UDG activation) and 95 °C for 2 min (enzyme activation), followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. The dissociation cycle for each primer pair was: from 95 to 60 °C (15 s each degree), 60 °C (1 min), from 60 to 95 °C (15 s each degree).Table 1 Primers used for gene expression analysis of larvae and pupae of Apis mellifera.Full size tableEffects of hormonal treatments on candidate genes and microRNAs in pupal stagesWe selected white-eyed pupae (Pw) for Juvenile Hormone-III (Sigma-Aldrich) treatment and brown-eyed pupae (Pb) for 20-hydroxyecdysone (Sigma-Aldrich) treatment (for details about the selected stages, see Supplementary Figure S3). Three points were sampled: 1 h, 1.5 h, and 24 h after treatment. The Pw group received a topical application of 1 µL of JH III (3 µg/µL) dissolved in acetone, while the control group was treated with 1 µL of acetone. The Pb group was injected with 1 µL of 20E (3 µg/µL) diluted in saline solution, with the control group receiving only saline (Supplementary Material and Methods). The injections were performed between the 3rd and 4th tergites. The pupae were kept in an incubator at 34 °C with 70% humidity. The development was monitored daily and documented using Zeiss Stereo Discovery V12 (Supplementary Figure S4). Five biological replicates were used in this step of the study.For the miRNAs, total RNA extraction, cDNA synthesis and qPCR, were conducted following steps mentioned in the previous section. For cDNA synthesis of the miRNAs we used the NCode™ miRNA First-Strand cDNA synthesis Kit (Invitrogen). For this, 2.5 µg of total RNA were used.Gene expression analysisFor the expression analysis of candidate genes, we tested reference genes as suggested by Lourenço et al.39, the most suitable for normalization was Rpl32 (see Supplementary Figure S5). In the case of miRNA expression analysis, miR-252a was used as the reference gene, as it exhibited the most stable expression profile across the tested conditions. The relative expression of target genes and miRNAs were assessed following the model (2−∆∆Cq method) proposed by Livak and Schmittigen40.Prediction and test of miRNAs regulation sites in the candidate genesFor prediction of miRNA binding sites on the developmental genes, the 3′ untranslated region (3′UTR) of the genes: Usp, EcR, Met, Kr-h1, Chd64, InR-2 and ftz-f1 were retrieved from the NCBI (Amel_HAv3.1, release 104) and the sequences of miR-34 and miR-281 from miRBase41. The regulatory elements were searched using RNAhybrid software42. Interactions with binding free energy below − 20 kcal/mol and a p-value < 0.05 were considered.To test if the predicted microRNA regulatory elements (MREs) were functional, we used the Dual-Luciferase Reporter Assay System (Promega) as described by Cristino et al.43. In summary, the MREs were inserted in the 3′UTR of a luciferase coding region yielding a chimeric PSYCHECK 2 (Promega) vector. HEK293T cells were co-transfected with chimeric vector and mimic miRNA for 24 h and left to rest until next day. Cells lysates were used to measure luciferase activity in a FlexStation 3 Multi-Mode Microplate Reader. Luciferase activity was compared between treatment and control groups using t-test (p-value < 0.05): 1) chimeric vector and mimic miRNA; 2) chimeric vector and scrambled miRNA; 3) chimeric vector only. Six biological replicates were used for each condition.Statistical analysesStatistical analyses were conducted using GraphPad Prism software (version 7, GraphPad Software, Inc., San Diego, CA, USA). Differences between the control and treated groups in the hormone treatment experiments were evaluated using a t-test, with significance set at p < 0.05. Correlations between gene expression profiles and hormone levels during pupal development were assessed by calculating Pearson’s correlation coefficient.ResultsExpression profile of Usp, EcR, Met, Kr-h1, Chd64, InR-2, ftz-f1 and Tai, and correlation with 20E and JH titersTo test if the expression profile of the genes is correlated with the hormones titers (JH titers obtained by Rembold44, and 20E titers from Rachinsky et al.45 and Pinto et al.46) we calculated the Pearson’s correlation coefficient (r), between the expression profile and circulating hormone titers in the hemolymph (Supplementary Table S1). We established a correlation threshold of ± 0.3, considering moderate correlations or stronger47. The genes were then categorized based on their responsiveness: JH-responsive or 20E-responsive.In the set of genes responsive to JH and 20E (Supplementary Table S1), two subclasses are identified: those with direct and inverse correlations. Among the JH-responsive genes, Kr-h1 (r = 0.9) and Usp (r = 0.6) show a direct correlation, where increased hormone titers lead to elevated gene expression (Fig. 1A). In contrast, InR-2 (r = − 0.31) exhibits an inverse correlation, where high hormone titers inhibit receptor expression (Fig. 1A). Regarding 20E-responsive genes, EcR (r = 0.74) and Met (r = 0.49) show direct correlations (Fig. 1B), while Chd64 (r = − 0.48), ftz-f1 (r = − 0.42), Tai (− 0.55) and Usp (r = − 0.43) display inverse correlations. Notably, Usp presents positive correlation with JH and a negative correlation with 20E during pre-imaginal development.Fig. 1Expression profile of the genes Usp, EcR, Met, Kr-h1, Chd64, InR-2, ftz-f1 and Tai during development of Apis mellifera workers, from the fifth larval instar (L5S1) to newly emerged (Ne), and their correlation with Juvenile Hormone (JH) and 20-hydroxyecdysone (20E). (A) JH responsive genes: Usp and Kr-h1 exhibit a positive correlation, while InR-2 shows a negative correlation with hormone fluctuations. (B) 20E responsive genes: EcR and Met shows a positive correlation, while Usp, Tai, Chd64, ftz-f1 and Tai exhibit a negative correlation with hormone fluctuations. The JH hormone titers were redrawn from Rembold44 and 20E titers redrawn from Rachinsky et al.45 and Pinto et al.46. Bars show means + -SEM (n = 3, samples processed individually). * EcR showed strong direct correlation with 20E when removing Pbm point (see discussion).Full size imageEffects of hormonal treatment on developmental genes and microRNAsThe effects of altered hormonal regimes on developmental genes and miRNAs expression were investigated by qPCR following the application of JH III and 20E on Pw and Pb pupae, respectively (Supplementary Figure S3). The treatment with 20E led to a decrease in the transcript levels of Usp (p = 0.01), EcR (p = 0.006), Met (p = 0.000001), Chd64 (p = 0.02), ftz-f1 (p = 0.003), miR-34 (p = 0.007), miR-281 (p = 0.03) and an increase of transcripts levels of Kr-h1 (p = 0.007) and InR-2 (p = 0.04) (Fig. 2A). The treatment with JH III increased the transcripts levels of Usp (p = 0.004), Kr-h1 (p = 0.0003), InR-2 (p = 0.01) and miR-34 (p = 0.009) (Fig. 2B). Transcripts levels of both EcR (p = 0.03) and miR-281 (p = 0.04) dropped in JH III-treated pupae 1 h after hormone application, rapidly increased 1h30 after the treatment (p = 0.006 and p = 0.003) and showed no difference to non-treated pupae 24 h after hormone application.Fig. 2Effects of hormonal treatment on the expression of the genes Usp, EcR, Kr-h1, Chd64, InR-2, ftz-f1, and the miRNAs miR-34 and miR-281 in Apis mellifera workers. (A) Gene expression in brown-eyed pupae (Pb) treated with 20-hydroxyecdysone (20E). (B) Gene expression in white-eyed pupae (Pw) treated with Juvenile Hormone III (JH III). We used 3 µg of each hormone, JH III and 20E, for the treatment of Pw and Pb pupae respectively. The pupae were sampled at 1, 1.5 and 24 h after treatment. Bars show means + -SEM (n = 5, samples processed individually). T-test: * < 0.05; ** < 0.01.Full size imageInteraction between developmental genes and their regulators miR-34 and miR-281Both JH III and 20E treatments affected the transcript levels of miR-34 and miR-281 (Fig. 2). To further explore the role of these regulatory RNAs, we searched for MREs in the 3′ UTR of developmental genes. MREs for miR-34 were predicted in the 3′ UTRs of the genes Kr-h1, Chd64, ftz-f1, EcR, and USP (Supplementary Figures S6 and S7), while MREs for miR-281 were identified in the 3′ UTRs of ftz-f1, EcR, USP, and InR-2 (Supplementary Figure S8). To test if the predicted MREs were functional, we constructed chimeric vectors by inserting the MREs into the 3′ end of the coding sequence of a reporter-gene (Fig. 3A–J). The activity of the reporter-gene was measured following the transfection of HEK293T cells with both the chimeric vector and mimic miRNAs. Mimic miR-34 decreased the gene-reporter activity when co-transfected with chimeric vectors bearing MREs predicted in the 3′ UTR of genes Chd64 (12.5%) (Fig. 3C), InR-2 (33.3%) (Fig. 3D) and Kr-h1 (25.6%) (Fig. 3E). Mimic miR-281 decreased the gene-reporter activity when co-transfected with chimeric vectors containing MREs predicted in the 3′ UTR of genes ftz-f1 (13.3%) (Fig. 3H) and EcR (39%) (Fig. 3I). Based on MREs results, a miRNA-target interaction network was proposed (Fig. 3J).Fig. 3Validation of the predicted binding sites for miR-34 and miR-281 in the 3′UTR of the genes EcR, Kr-h1, Chd64, InR-2 and ftz-f1, by luciferase assay. Changes in the ratio of chimeric Renilla luciferase (Ren) and firefly luciferase (FF) activity is shown for each gene. The chimeric Ren luciferase is linked to the 3′UTR of target genes containing the predicted binding sites for miR-34 and miR-281. Schematic representation of the positive control, perfect pairing with the miR-34 and miR-281 sequence (A and F). The positive control (B and G) and the 3′UTRs of the Chd64 (B), InR2 (C), Kr-h1 (D) ftz-f1 (H), and EcR (I) genes were directly regulated by miR-281. (J) Network of interactions based on the validated binding sites between the miR-34 and miR-281, in the 3′UTR of the target genes: EcR, Kr-h1, Chd64, InR-2 and ftz-f1. Bars show means + -SEM (n = 6) T-test: * < 0.05.Full size imageDiscussionIn A. mellifera, the genes Usp38,48,49,50, EcR33,51,52, InR-253,54 and ftz-f155,56,57, have been studied in the context of pre-imaginal development. In contrast, research on Kr-h1, has focused on the nervous system58,59,60,61, and for Chd64 the information available is behavior-related62,63,64.Our results provide a comprehensive analysis of key genes involved in the control of late pre-imaginal development in honeybees (Fig. 4). We show that Usp responds to both morphogenetic hormones, and that Kr-h1 responds very precisely to JH. Meanwhile, Met, which is known to participate only in the JH pathway, followed the 20E titers in the assessed stages. The expression of the genes ftz-f1 and Chd64 showed to be reliable biomarkers for the metamorphic molt and imaginal molt. Finally, we expand the regulatory network by demonstrating that miR-34 and miR-281 also interact with these genes, pointing to them as regulators of pupal development in honeybees.Fig. 4Conceptual model of interactions between the morphogenetic hormones (JH/20E) and the genes Usp, EcR, Met, Kr-h1, Chd64, InR-2, ftz-f1, Tai, as well as the miRNAs miR-34 and miR-281 during the pupal development of Apis mellifera workers. The modulation of JH titers was obtained from Rembold44, while the modulation of 20E titers was based on the studies of Rachinsky et al.45 and Pinto et al.46. The classification of developmental stages followed Michelette and Soares34.Full size imageMet expression is affected by 20E and presents stage-specific response to JH-III stimuliOur results point to a novel role for Met in the pupal development of honeybees, independent of the JH pathway, as its expression responds to 20E rather than JH stimuli. Transcript levels of Met increase during the pupal stages (Pdp and Pb) in response to the 20E peak at the Pp stage (Fig. 1B), and when 20E levels are maintained high (by injections of 20E in Pb), Met expression is also affected. On the other hand, treatment with JH had no effect on Met transcript levels during pupal development, when JH levels are basal in honeybees44 (Fig. 1B).Met is widely recognized for its role as JH receptor65,66,67 being described in distinct JH-dependent biological processes such reproduction68, embryonic development69 and development of follicular ephitelial cells70. However, its significant expression during the final developmental stages (Fig. 1B), when JH levels are low, raises intriguing questions. In Tribolium castaneum, in a study focused on characterizing Kr-h1 function, high Met expression levels can be similarly observed during the pupal phase, when JH is minimal71, suggesting Met may adopt a JH-independent function at this stage.We also provide evidence that Met responds differentially to JH stimuli depending on the developmental stage. Supplementary Figure S9 offers additional findings, showing that JH III influences Met expression in L2 larvae but does not affect Met expression in pupal stages, prompting the question: Do tissues lose their competence to respond to JH during pupal development due to the hormone’s sharp decline? At this stage, the organism is preparing for the imaginal molt, making it crucial that no effects from JH interfere with this process. Thus, Met is JH-responsive in early stages and 20E-responsive in late development. Bitra and Palli72 demonstrated Met interacting with ECR and USP in Drosophila cell line (two-hybrid interaction assay – L57 cells) while Guo et al.73 showed Met responding to 20E during pupation in Bombyx mori, again by interacting with ECR and USP. Although there is no evidence of interaction between ECR/USP and MET in honeybees, it still supports that MET plays a role in the 20E pathway during pupal development in this species. Additionally, the highest expression levels of Met (Pdp and Pb) coincide precisely with a decrease in Usp and EcR expression levels (Fig. 1A,B), demonstrating that Met and Usp/EcR have opposing expression profiles in these moments. This indicates a potential interaction—an assertion that requires further experimental validation.Usp, a rapid response gene, plays a role in the JH and 20E pathwaysUSP was identified as component of the ecdysone response complex in 199374. Years later, Jones and Sharp14 used fluorescence-based binding assays to demonstrate the USP binding affinity for JH III, indicating a dual role for USP in both hormones pathways. In honeybees, Barchuk et al.38 demonstrated that Usp expression quickly responds to JH, confirming USP as primary JH binding protein in this species. Our findings show that Usp has a positive correlation with JH while displaying a negative correlation with 20E, supporting that during late development of A. mellifera, Usp plays a role in both hormones signalling pathways. In 2008, Barchuck et al.49 first proposed a negative regulation of Usp by the morphogenetic hormone 20E in A. mellifera, based on expression observations. Our hormonal treatment confirms this, as Usp was downregulated by exogenous 20E at all three tested points. This regulatory pattern is consistent across other bee species, the stingless bees Melipona scutellaris and Scaptotrigona depilis75.The hormonal treatments reveal that while some genes are quickly responsive, others respond more slowly. Usp for example responds to JH treatment within 1 h, but its original levels are restored 30 min after that. In contrast, the expression of genes responsive to 20E are not restored even 24 h post-treatment. Barchuk et al.38 also observed a rapid response of Usp to JH treatment, however, the original expression levels did not restore even 4 h and 30 min post-treatment. The authors used a high dosage of JH III (10 µg/µL), approximately 3.3 times higher than we used here, which probably explains the different results obtained.Additionally, the expression profile of Usp closely resembles that of Tai, suggesting a functional partnership. TAI is a highly versatile component involved in regulating both 20E and JH gene cascades. It functions as a steroid coactivator for the 20E receptor complex EcR/USP76 and also forms heterodimers with MET in response to JH to repress premature molting19. Our results reveal a strong negative correlation between Tai expression (r = − 0.55) and 20E titer fluctuations, similar to Usp (r = − 0.43), during late honeybee development, suggesting an interaction between TAI and USP in this species. However, our approach focuses on the pupal phase, with only the last larval stages (spinning stages) represented in our data. Tai may also correlate with JH, as it mediates JH’s antimetamorphic action through dimerization with MET. However, this effect is expected to occur during the larval molts, and our results do not include this stages.Hormonal treatments validate the expression profile correlationInteractions between hormones and the genes regulating development create a complex and intricate network, where the perturbation of a single component can significantly influence the entire system. Zhang et al.77 demonstrated that Kr-h1, a well-established target gene for JH, affects the synthesis of 20E by repressing the expression of genes involved in ecdysteroidogenesis in Drosophila and Bombyx. Given the cross-regulation that exists between hormonal pathways, where one hormone can modulate the effects of another, we correlated the earliest time point of hormonal treatment (1 h post-treatment) with gene expression profile data obtained from specific developmental stages that align with these treatments. For JH III, we analyzed data from the Pw onwards, while for 20E, we focused on the expression data starting from the Pb onwards. This integrated approach enabled us to validate 10 out of the 14 hypothesized relationships presented in this study (see Supplementary Figure S10).Kr-h1 and InR-2 have opposite patterns under JH stimuliThe transcriptional regulation of JH target genes, such as Kr-h1, is mediated by interaction between USP, MET/GCE and Taiman for insects in general6,19,78. The strong positive correlation found between Kr-h1 expression pattern and JH levels show Kr-h1 as a reliable biomarker for JH titers in honeybees. Furthermore, this relationship is highly conserved across different insect species, spanning from hemimetabolous79,80 to holometabolous insects71,78.The inverse pattern observed between InR-2 and JH is consistent with the findings of de Azevedo and Hartfelder54. This relationship was also identified in D. melanogaster, where the insulin-like peptide 6 (DILP6), another component of the insuling pathway, negatively regulates JH titers81. However, this poses a contradiction, as InR-2 levels are low during the feeding stages of A. mellifera larvae and peak during the pupal stages, when the organism does not feed. Perez-Hedo et al.82 demonstrated that during starvation in Aedes aegypti, the expression of the insulin receptor increases alongside expression of Forkhead-box O (FoxO), a terminal transcription factor of the insulin signaling pathway. FoxO coordinates developmental timing and tissue differentiation83,84 being active during development, especially during pupal development when feeding ceases. The peaks in InR-2 expression observed at the Pp and Pbm stages (Fig. 1A) can likely be attributed to the activity of FoxO during these periods. In another aspect, FoxO is also known to regulate 20E biosynthesis85 and to play a role in JH degradation84,86. Regarding the modulation of these two key morphogenetic hormones, during honeybee development, FoxO is probably involved in the metamorphic transition, in which JH levels are decreasing44 and the 20E levels are increasing45,46, underscoring the critical importance of FoxO during these developmental transitions. Despite the characterization of the FoxO gene family in the genome of A. mellifera87, its function in the context of honeybee development remains to be elucidated.Expression profile of ftz-f1 and Chd64 indicates metamorphic and imaginal moltThe genes ftz-f1 and Chd64 exhibited similar expression patterns in response to 20E stimuli. As shown in their expression profiles (Fig. 1B), an inverse relationship with hormone titers is observed, which is further confirmed by hormonal treatments, where both genes display the same response to both hormonal stimuli (Fig. 2). This inverse correlation between gene expression and hormone levels has also been documented by Mello et al.57 for the ftz-f1 transcription factor, which requires a decay in 20E levels to increase its expression in honey bees, as it does in D. melanogaster88. This delayed response coincides with the molting process, occurring after the decay of 20E. The stages examined in this study encompass both metamorphic and imaginal molts, specifically the transitions from PP3 to Pw and from Pbd to Ne. Since the expression of ftz-f1 and Chd64 peaks during these stages, their expression profiles are biological markers that distinctly indicate the occurrence of metamorphic and imaginal molts in A. mellifera.EcR expression peaks after 20E titers dropEcR expression peaks immediately after the plummet of 20E titers in the Pbm stage (Fig. 1B). Truman et al.89 demonstrated in Manduca sexta that a rise and subsequent fall in ecdysteroid titers is a prerequisite to turn tissues responsive to the hormone. Similarly, in honey bees, Mello et al.33 showed that exogenous 20E treatments repress EcR expression when 20E levels are high, but EcR expression increases once 20E levels decline. This dynamic explains the peak in EcR expression profile (Pbm stage) we observed after the significant decline in 20E (Fig. 1B). Due to this specific regulatory mechanism between 20E and EcR, we did not include the Pbm value (highest 20E peak) when calculating the Pearson’s correlation between the EcR expression profile and 20E titers (Supplementary Table S1). The expression profile obtained for EcR and its response to 20E treatments aligns with the findings of Mello et al.33, with exception to the results observed under JH treatment. However, the dosage and the chosen stage used for JH treatment is different, Mello and colleagues used Pbl, close to the final pupal development, while we aimed at the beginning of pupal development using Pw. In this study, we assessed isoform A of the gene, which predominates during post-embryonic development, while isoform B is more highly expressed during embryonic stages33. Additionally, Apis mellifera ligustica also exhibits a similar positive correlation between 20E and EcR52.miR-34 and miR-281 are involved in the regulation of developmental genesThe treatment with JH III notably increased the transcript levels of miR-34, along with its downstream targets, Kr-h1 and InR-2. This observation raises the question of whether miR-34 is directly involved in regulating the decrease of Kr-h1 and InR-2 or if the increase in transcript levels observed is merely a response to JH III. In a study by Liu and colleagues29, the authors demonstrated that miR-34-KO silkworms (with miRNA ablation via the CRISPR/Cas9 system) exhibited a severe developmental delay during the larval stages, highlighting miR-34 as a regulator of larval growth. This regulatory axis seems to be conserved, as in honeybees, in addition to the increase in miR-34 expression induced by JH III treatment (Fig. 2B), we also confirmed interaction between miR-34 and Kr-h1 (Fig. 3E). Although the relationship between JH III and miR-34 has not been extensively documented in the literature, our findings suggest a regulatory connection. Furthermore, several other miRNAs have been reported with binding sites in the 3′ UTR of Kr-h122, such as miR-2 (B. germanica)90, miR-927 (D. melanogaster)91, and let-7 and miR-278 (Locusta migratoria)92, reinforcing the cooperative role of hormones and miRNAs in regulating gene expression.While establishing evolutionarily and functionally conserved patterns remains a complex task due to the vast diversity of species and their life histories, our data aligns with prior research. MicroRNAs are known to integrate molecular regulatory networks modulating insect developmental timing and stage transitions93. The scientific literature provides conclusive evidence of miRNA-hormone interactions associated with body transformations in Hexapoda. In particular, aspects of organismal growth, including cell proliferation and differentiation, are controled by some JH- or 20E-responsive miRNAs as well as by miRNA-mediated activation of JH- or 20E-signaling94. An increasing number of microRNAs crucial for the fine-tuned regulation of insect ontogeny have been well documented, comprising (but not limited to) bantam, let-7, iab-8, mir-13a, miR-13b, miR-100, miR-125, miR-14, miR-2, miR-252a, miR-252b, miR-263, miR-278, miR-281, miR-34, miR-71, miR-8, miR-9, miR-927, miR-96528,93,94,95,96,97,98.Both miR-34 and miR-281 were upregulated by JH and downregulated by 20E treatment in the pupal development of honeybees. The expression data for miR-34 and EcR following 20E treatment (Fig. 2A) align with the findings of Mello et al.33, which reported a significant downregulation of miR-34 expression in EcR-knockdown bees. The highly conserved miR-34 exhibited a similar response to JH and 20E treatments in D. melanogaster embryonic cells96, reinforcing that the interplay between these hormones and miR-34 is critical for regulating development and supporting the idea of conserved hormonal regulatory mechanism across insect species.MiR-281 and EcR displayed similar response to JH and 20E treatments in the pupal development of honeybees (Fig. 2A,B) and miR-281 directly targets EcR transcripts (Fig. 3I,J). In Bombyx mori, the expression of miR-281 was downregulated by 20E but not affected by JH treatments. In this species, miR-281 and EcR are co-expressed in the Malpighian tubules from larvae to pupation, with miR-281 overexpression downregulating EcR protein levels, while miR-281 suppression leads to EcR protein upregulation31. Our findings reveal a conserved response of miR-34 and miR-281 to variations in JH and 20E levels, and by validating the interactions of these microRNAs with the developmental genes, we add a new layer of complexity to the understanding of hormonal regulation of gene expression, in this context, fine-tuning the pre-imaginal development process. Combining these findings with the presented discussion and evidence from the literature, we propose a model (Fig. 4) that illustrates a potential regulatory mechanism underlying pupal development in honeybees.Here, we characterized gene expression dynamics during the late development of A. mellifera. Notably, an increase in the transcriptional levels of key genes (Usp, Kr-h1, InR-2, Chd64, ftz-f1 and Tai) as development progresses to the later pupal stages, takes place immediately after the decay in 20E titers. Furthermore, a well-established JH receptor, Met, appears to assume a distinct role in pupal development, independent of this hormone. The findings regarding the physical interaction between microRNAs and critical developmental genes further advance our understanding of the relationship between these complexes and the two primary hormones, 20E and JH, which govern pupal progression.
Data availability
All data generated or analysed during this study are included in this published article and its supplementary information files.
ReferencesBellés, X. The metamorphosis of insects and their regulation. C.R. Biol. 342(7–8), 254–256. https://doi.org/10.1016/j.crvi.2019.09.009 (2019).Article
MATH
Google Scholar
Gilbert, S. F. & Barresi, M. J. F. Developmental Biology (Sinauer Associates, 2016).MATH
Google Scholar
Dubrovsky, E. B. Hormonal cross talk in insect development. Trends Endocrinol. Metab. 16(1), 6–11. https://doi.org/10.1016/j.tem.2004.11.003 (2005).Article
MathSciNet
PubMed
MATH
Google Scholar
Truman, J. W. & Riddiford, L. M. The evolution of insect metamorphosis: A developmental and endocrine view. Philos. Trans. R. Soc. B 374(1783), 20190070. https://doi.org/10.1098/rstb.2019.0070 (2019).Article
MATH
Google Scholar
Truman, J. W. The evolution of insect metamorphosis. Curr. Biol. 29(23), R1252–R1268. https://doi.org/10.1016/j.cub.2019.10.009 (2019).Article
PubMed
MATH
Google Scholar
Jindra, M., Bellés, X. & Shinoda, T. Molecular basis of juvenile hormone signaling. Curr. Opin. Insect Sci. 11, 39–46. https://doi.org/10.1016/j.cois.2015.08.004 (2015).Article
PubMed
MATH
Google Scholar
Belgacem, Y. H. & Martin, J. R. Hmgcr in the corpus allatum controls sexual dimorphism of locomotor activity and body size via the insulin pathway in Drosophila. PLoS ONE 2(1), e187. https://doi.org/10.1371/journal.pone.0000187 (2007).Article
ADS
PubMed
PubMed Central
MATH
Google Scholar
Liu, P. C., Wang, J. X., Song, Q. S. & Zhao, X. F. The participation of calponin in the cross talk between 20-hydroxyecdysone and juvenile hormone signaling pathways by phosphorylation variation. PLoS ONE 6(5), e19776. https://doi.org/10.1371/journal.pone.0019776 (2011).Article
ADS
PubMed
PubMed Central
Google Scholar
Xu, Q. Y. et al. Functional characterization of ultraspiracle in Leptinotarsa decemlineata using RNA interference assay. Insect Mol. Biol. 28(5), 676–688. https://doi.org/10.1111/imb.12580 (2019).Article
PubMed
MATH
Google Scholar
Ali, M. S. & Takaki, K. Transcriptional regulation of cuticular genes during insect metamorphosis. Front. Biosci. (Landmark Edition) 25, 106–117. https://doi.org/10.2741/4796 (2020).Article
MATH
Google Scholar
Bellés, X. Krüppel homolog 1 and E93: The doorkeeper and the key to insect metamorphosis. Arch. Insect Biochem. Physiol. 103(3), e21609. https://doi.org/10.1002/arch.21609 (2020).Article
PubMed
Google Scholar
Okude, G. et al. Molecular mechanisms underlying metamorphosis in the most-ancestral winged insect. Proc. Natl. Acad. Sci. https://doi.org/10.1073/pnas.2114773119 (2022).Article
PubMed
PubMed Central
Google Scholar
Truman, J. W. & Riddiford, L. M. Drosophila postembryonic nervous system development: A model for the endocrine control of development. Genetics 223(3), iyac184. https://doi.org/10.1093/genetics/iyac184 (2023).Article
PubMed
PubMed Central
Google Scholar
Jones, G. & Sharp, P. A. Ultraspiracle: An invertebrate nuclear receptor for juvenile hormones. Proc. Natl. Acad. Sci. 94(25), 13499–13503. https://doi.org/10.1073/pnas.94.25.13499 (1997).Article
ADS
PubMed
PubMed Central
MATH
Google Scholar
Jones, G. et al. Ligand binding pocket function of Drosophila USP is necessary for metamorphosis. Gen. Comp. Endocrinol. 182, 73–82. https://doi.org/10.1016/j.ygcen.2012.11.009 (2013).Article
PubMed
MATH
Google Scholar
Li, Y., Zhang, Z., Robinson, G. E. & Palli, S. R. Identification and characterization of a juvenile hormone response element and its binding proteins. J. Biol. Chem. 282(52), 37605–37617. https://doi.org/10.1074/jbc.M704595200 (2007).Article
PubMed
Google Scholar
Palli, S. R. Juvenile hormone receptor Methoprene tolerant: Functions and applications. Vitam. Horm. 123, 619–644. https://doi.org/10.1016/bs.vh.2023.03.002 (2023).Article
PubMed
MATH
Google Scholar
Li, K., Jia, Q. Q. & Li, S. Juvenile hormone signaling–a mini review. Insect Sci. 26(4), 600–606. https://doi.org/10.1111/1744-7917.12614 (2019).Article
PubMed
MATH
Google Scholar
Lozano, J., Kayukawa, T., Shinoda, T. & Bellés, X. A role for Taiman in insect metamorphosis. PLoS Genet. 10(10), e1004769. https://doi.org/10.1371/journal.pgen.1004769 (2014).Article
PubMed
PubMed Central
Google Scholar
Li, M., Mead, E. A. & Zhu, J. Heterodimer of two bHLH-PAS proteins mediates juvenile hormone-induced gene expression. Proc. Natl. Acad. Sci. 108(2), 638–643. https://doi.org/10.1073/pnas.1013914108 (2011).Article
ADS
PubMed
Google Scholar
Zhang, Z., Xu, J., Sheng, Z., Sui, Y. & Palli, S. R. Steroid receptor co-activator is required for juvenile hormone signal transduction through a bHLH-PAS transcription factor, methoprene tolerant. J. Biol. Chem. 286(10), 8437–8447. https://doi.org/10.1074/jbc.M110.191684 (2011).Article
PubMed
Google Scholar
He, Q. & Zhang, Y. Kr-h1, a cornerstone gene in insect life history. Front. Physiol. 13, 905441. https://doi.org/10.3389/fphys.2022.905441 (2022).Article
PubMed
PubMed Central
Google Scholar
Bellés, X., Martín, D. & Piulachs, M. D. The mevalonate pathway and the synthesis of juvenile hormone in insects. Annu. Rev. Entomol. 50(1), 181–199. https://doi.org/10.1146/annurev.ento.50.071803.130356 (2005).Article
PubMed
Google Scholar
Borras-Castells, F. et al. Juvenile hormone biosynthesis in adult Blattella germanica requires nuclear receptors Seven-up and FTZ-F1. Sci. Rep. 7, 40234. https://doi.org/10.1038/srep40234 (2017).Article
ADS
PubMed
PubMed Central
MATH
Google Scholar
Broadus, J., McCabe, J. R., Endrizzi, B., Thummel, C. S. & Woodard, C. T. The Drosophila βFTZ-F1 orphan nuclear receptor provides competence for stage-specific responses to the steroid hormone ecdysone. Mol. Cell 3(2), 143–149. https://doi.org/10.1016/S1097-2765(00)80305-6 (1999).Article
PubMed
Google Scholar
Li, C. et al. Insect metamorphosis is regulated differently between sexes by members of a microRNA cluster. eLife 13, RP100499. https://doi.org/10.1101/2024.06.25.600558 (2024).Article
Google Scholar
Ma, X. et al. Large-scale annotation and evolution analysis of MiRNA in insects. Genome Biol. Evol. 13(5), e083. https://doi.org/10.1093/gbe/evab083 (2021).Article
MATH
Google Scholar
Song, J. & Zhou, S. Post-transcriptional regulation of insect metamorphosis and oogenesis. Cell. Mol. Life Sci. 77(10), 1893–1909. https://doi.org/10.1007/s00018-019-03361-5 (2020).Article
PubMed
MATH
Google Scholar
Liu, Z. et al. miR-34 regulates larval growth and wing morphogenesis by directly modulating ecdysone signalling and cuticle protein in Bombyx mori. RNA Biol. 17(9), 1342–1351. https://doi.org/10.1080/15476286.2020.1767953 (2020).Article
PubMed
PubMed Central
MATH
Google Scholar
Soares, M. P. M., Pinheiro, D. G., de Paula Freitas, F. C., Simões, Z. L. P. & Bitondi, M. M. G. Transcriptome dynamics during metamorphosis of imaginal discs into wings and thoracic dorsum in Apis mellifera castes. BMC Genom. 22, 1–24. https://doi.org/10.1186/s12864-021-08040-z (2021).Article
Google Scholar
Jiang, J. et al. MicroRNA-281 regulates the expression of ecdysone receptor (EcR) isoform B in the silkworm, Bombyx mori. Insect Biochem. Mol. Biol. 43(8), 692–700. https://doi.org/10.1016/j.ibmb.2013.05.002 (2013).Article
PubMed
MATH
Google Scholar
Wang, Z. Z. et al. Parasitic insect-derived miRNAs modulate host development. Nat. Commun. 9(1), 1–9. https://doi.org/10.1038/s41467-018-04504-1 (2018).Article
ADS
Google Scholar
Mello, T. R. et al. Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera). Front. Genet. 5, 445. https://doi.org/10.3389/fgene.2014.00445 (2014).Article
PubMed
PubMed Central
Google Scholar
Michelette, E. R. & Soares, A. E. E. Characterization of preimaginal developmental stages in Africanized honey bee workers (Apis mellifera L.). Apidologie 24(4), 431–440. https://doi.org/10.1051/apido:19930410 (1993).Article
MATH
Google Scholar
Honeybee Genome Sequencing Consortium. Insights into social insects from the genome of the honeybee Apis mellifera. Nature 443(7114), 931. https://doi.org/10.1038/nature05260 (2006).Article
ADS
Google Scholar
Elsik, C. G. et al. Hymenoptera genome database: Integrating genome annotations in HymenopteraMine. Nucleic Acids Res. 44(D1), D793–D800. https://doi.org/10.1093/nar/gkv1208 (2015).Article
PubMed
PubMed Central
Google Scholar
Carver, T., Harris, S. R., Berriman, M., Parkhill, J. & McQuillan, J. A. Artemis: An integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics 28(4), 464–469. https://doi.org/10.1093/bioinformatics/btr703 (2012).Article
PubMed
Google Scholar
Barchuk, A. R., Maleszka, R. & Simões, Z. L. P. Apis mellifera ultraspiracle: cDNA sequence and rapid up-regulation by juvenile hormone. Insect Mol. Biol. 13(5), 459–467. https://doi.org/10.1111/j.0962-1075.2004.00506.x (2004).Article
PubMed
Google Scholar
Lourenço, A. P., Mackert, A., dos Santos Cristino, A. & Simões, Z. L. P. Validation of reference genes for gene expression studies in the honey bee, Apis mellifera, by quantitative real-time RT-PCR. Apidologie 39(3), 372–385. https://doi.org/10.1051/apido:2008015 (2008).Article
Google Scholar
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25(4), 402–408. https://doi.org/10.1006/meth.2001.1262 (2001).Article
PubMed
MATH
Google Scholar
Kozomara, A., Birgaoanu, M. & Griffiths-Jones, S. miRBase: From microRNA sequences to function. Nucleic Acids Res. 47(D1), D155–D162. https://doi.org/10.1093/nar/gky1141 (2019).Article
PubMed
Google Scholar
Krüger, J. & Rehmsmeier, M. RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 34(suppl_2), W451–W454. https://doi.org/10.1093/nar/gkl243 (2006).Article
PubMed
PubMed Central
MATH
Google Scholar
Cristino, A. S. et al. Neuroligin-associated microRNA-932 targets actin and regulates memory in the honeybee. Nat. Commun. 5(1), 5529. https://doi.org/10.1038/ncomms6529 (2014).Article
ADS
PubMed
MATH
Google Scholar
Rembold, H. Caste specific modulation of juvenile hormone titers in Apis mellifera. Insect Biochem. 17(7), 1003–1006. https://doi.org/10.1016/0020-1790(87)90110-7 (1987).Article
MATH
Google Scholar
Rachinsky, A., Strambi, C., Strambi, A. & Hartfelder, K. Caste and metamorphosis: Hemolymph titers of juvenile hormone and ecdysteroids in last instar honeybee larvae. Gen. Comp. Endocrinol. 79(1), 31–38. https://doi.org/10.1016/0016-6480(90)90085-Z (1990).Article
PubMed
MATH
Google Scholar
Pinto, L. Z., Hartfelder, K., Bitondi, M. G. & Simões, Z. L. P. Ecdysteroid titers in pupae of highly social bees relate to distinct modes of caste development. J. Insect. Physiol. 48(8), 783–790. https://doi.org/10.1016/S0022-1910(02)00103-8 (2002).Article
PubMed
Google Scholar
Cohen, J. Statistical Power Analysis for the Behavioral Sciences 2nd edn. (Lawrence Erlbaum Associates, 1988).MATH
Google Scholar
Barchuk, A. R. et al. Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera. BMC Dev. Biol. 7, 1–19. https://doi.org/10.1186/1471-213X-7-70 (2007).Article
Google Scholar
Barchuk, A. R., Figueiredo, V. L. C. & Simões, Z. L. Downregulation of ultraspiracle gene expression delays pupal development in honeybees. J. Insect. Physiol. 54(6), 1035–1040. https://doi.org/10.1016/j.jinsphys.2008.04.006 (2008).Article
PubMed
MATH
Google Scholar
Corby-Harris, V., Snyder, L. & Meador, C. Fat body lipolysis connects poor nutrition to hypopharyngeal gland degradation in Apis mellifera. J. Insect. Physiol. 116, 1–9. https://doi.org/10.1016/j.jinsphys.2019.04.001 (2019).Article
PubMed
Google Scholar
Paul, R. K., Takeuchi, H. & Kubo, T. Expression of two ecdysteroid-regulated genes, Broad-Complex and E75, in the brain and ovary of the honeybee (Apis mellifera L.). Zool. Sci. 23(12), 1085–92. https://doi.org/10.2108/zsj.23.1085 (2006).Article
Google Scholar
Yu, J. et al. 20-hydroxyecdysone upregulates ecdysone receptor (ECR) gene to promote pupation in the honeybee, Apis mellifera Ligustica. Integr. Comp. Biol. 63(2), 288–303. https://doi.org/10.1093/icb/icad077 (2023).Article
PubMed
Google Scholar
Wheeler, D. E., Buck, N. O. R. M. A. N. & Evans, J. D. Expression of insulin pathway genes during the period of caste determination in the honey bee Apis mellifera. Insect Mol. Biol. 15(5), 597–602. https://doi.org/10.1111/j.1365-2583.2006.00681.x (2006).Article
PubMed
PubMed Central
Google Scholar
de Azevedo, S. V. & Hartfelder, K. The insulin signaling pathway in honey bee (Apis mellifera) caste development—differential expression of insulin-like peptides and insulin receptors in queen and worker larvae. J. Insect. Physiol. 54(6), 1064–1071. https://doi.org/10.1016/j.jinsphys.2008.04.009 (2008).Article
PubMed
Google Scholar
Hepperle, C. & Hartfelder, K. Differentially expressed regulatory genes in honey bee caste development. Naturwissenschaften 88, 113–116. https://doi.org/10.1007/s001140000196 (2001).Article
ADS
PubMed
Google Scholar
Soares, M. P. M. et al. Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera. BMC Genom. 14, 1–17. https://doi.org/10.1186/1471-2164-14-576 (2013).Article
Google Scholar
Mello, T. R. P. et al. Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: A positive loop linking juvenile hormone, ftz-f1, and vitellogenin. Insect Mol. Biol. 28(1), 145–159. https://doi.org/10.1111/imb.12541 (2019).Article
PubMed
MATH
Google Scholar
Grozinger, C. M. & Robinson, G. E. Endocrine modulation of a pheromone-responsive gene in the honeybee brain. J. Comp. Physiol. A 193, 461–470. https://doi.org/10.1007/s00359-006-0202-x (2007).Article
MATH
Google Scholar
Fussnecker, B. & Grozinger, C. Dissecting the role of Kr-h1 brain gene expression in foraging behavior in honeybees (Apis mellifera). Insect Mol. Biol. 17(5), 515–522. https://doi.org/10.1111/j.1365-2583.2008.00819.x (2008).Article
PubMed
Google Scholar
Moda, L. M. et al. Nutritionally driven differential gene expression leads to heterochronic brain development in honeybee castes. PLoS ONE 8(5), e64815. https://doi.org/10.1371/journal.pone.0064815 (2013).Article
ADS
PubMed
PubMed Central
Google Scholar
Gospocic, J. et al. Kr-h1 maintains distinct caste-specific neurotranscriptomes in response to socially regulated hormones. Cell 184(23), 5807–5823. https://doi.org/10.1016/j.cell.2021.10.006 (2021).Article
PubMed
PubMed Central
Google Scholar
Ament, S. A. et al. The transcription factor ultraspiracle influences honeybee social behavior and behavior-related gene expression. PLoS Genet. 8(3), e1002596. https://doi.org/10.1371/journal.pgen.1002596 (2012).Article
PubMed
PubMed Central
MATH
Google Scholar
Pandey, A. & Bloch, G. Juvenile hormone and ecdysteroids as major regulators of brain and behavior in bees. Curr. Opin. Insect Sci. 12, 26–37. https://doi.org/10.1016/j.cois.2015.09.006 (2015).Article
MATH
Google Scholar
Lino, S., Oya, S., Kakutani, T., Kohno, H. & Kubo, T. Identification of ecdysone receptor target genes in the worker honey bee brains during foraging behavior. Sci. Rep. 13(1), 10491. https://doi.org/10.1038/s41598-023-37001-7 (2023).Article
Google Scholar
Baumann, A., Barry, J., Wang, S., Fujiwara, Y. & Wilson, T. G. Paralogous genes involved in juvenile hormone action in Drosophila melanogaster. Genetics 185(4), 1327–1336. https://doi.org/10.1534/genetics.110.116962 (2010).Article
PubMed
PubMed Central
Google Scholar
Lozano, J. & Bellés, X. Role of Methoprene-tolerant (Met) in adult morphogenesis and in adult ecdysis of Blattella germanica. PLoS ONE 9(7), e103614. https://doi.org/10.1371/journal.pone.0103614 (2014).Article
ADS
PubMed
PubMed Central
Google Scholar
Jindra, M., Uhlirova, M., Charles, J. P., Smykal, V. & Hill, R. J. Genetic evidence for function of the bHLH-PAS protein Gce/Met as a juvenile hormone receptor. PLoS Genet. 11(7), e1005394. https://doi.org/10.1371/journal.pgen.1005394 (2015).Article
PubMed
PubMed Central
Google Scholar
Zou, Z. et al. Juvenile hormone and its receptor, methoprene-tolerant, control the dynamics of mosquito gene expression. Proc. Natl. Acad. Sci. 110(24), E2173–E2181. https://doi.org/10.1073/pnas.1305293110 (2013).Article
PubMed
PubMed Central
MATH
Google Scholar
Naruse, S., Washidu, Y., Miura, K., Shinoda, T. & Minakuchi, C. Methoprene-tolerant is essential for embryonic development of the red flour beetle Tribolium castaneum. J. Insect. Physiol. 121, 104017. https://doi.org/10.1016/j.jinsphys.2020.104017 (2020).Article
PubMed
Google Scholar
HuangFu, N. et al. Dynamic transcriptome analysis and Methoprene-tolerant gene knockdown reveal that juvenile hormone regulates oogenesis and vitellogenin synthesis in Propylea Japonica. Genomics 113(4), 2877–2889. https://doi.org/10.1016/j.ygeno.2021.06.002 (2021).Article
PubMed
MATH
Google Scholar
Minakuchi, C., Namiki, T. & Shinoda, T. Krüppel homolog 1, an early juvenile hormone-response gene downstream of Methoprene-tolerant, mediates its anti-metamorphic action in the red flour beetle Tribolium castaneum. Dev. Biol. 325(2), 341–350. https://doi.org/10.1016/j.ydbio.2008.10.016 (2009).Article
PubMed
Google Scholar
Bitra, K. & Palli, S. R. Interaction of proteins involved in ecdysone and juvenile hormone signal transduction. Arch. Insect Biochem. Physiol. 70(2), 90–105. https://doi.org/10.1002/arch.20281 (2009).Article
PubMed
Google Scholar
Guo, E. et al. MET is required for the maximal action of 20-hydroxyecdysone during Bombyx metamorphosis. PLoS ONE 7(12), e53256. https://doi.org/10.1371/journal.pone.0053256 (2012).Article
ADS
PubMed
PubMed Central
Google Scholar
Yao, T. P. et al. Functional ecdysone receptor is the product of EcR and Ultraspiracle genes. Nature 366(6454), 476–479. https://doi.org/10.1038/366476a0 (1993).Article
ADS
PubMed
MATH
Google Scholar
Teles, A. C. A. S., Mello, T. R. P., Barchuk, A. R. & Simões, Z. L. P. Ultraspiracle of the stingless bees Melipona scutellaris and Scaptotrigona depilis: cDNA sequence and expression profiles during pupal development. Apidologie 38(5), 462–471. https://doi.org/10.1051/apido:2007035 (2007).Article
Google Scholar
Xu, Q. Y., Du, J. L., Mu, L. L., Guo, W. C. & Li, G. Q. Importance of Taiman in larval-pupal transition in Leptinotarsa decemlineata. Front. Physiol. 10, 724. https://doi.org/10.3389/fphys.2019.00724 (2019).Article
PubMed
PubMed Central
Google Scholar
Zhang, T. et al. Krüppel homolog 1 represses insect ecdysone biosynthesis by directly inhibiting the transcription of steroidogenic enzymes. Proc. Natl. Acad. Sci. 115(15), 3960–3965. https://doi.org/10.1073/pnas.1800435115 (2018).Article
ADS
PubMed
PubMed Central
Google Scholar
Kayukawa, T. et al. Transcriptional regulation of juvenile hormone-mediated induction of Krüppel homolog 1, a repressor of insect metamorphosis. Proc. Natl. Acad. Sci. 109(29), 11729–11734. https://doi.org/10.1073/pnas.1204951109 (2012).Article
ADS
PubMed
PubMed Central
Google Scholar
Lozano, J. & Bellés, X. Conserved repressive function of Krüppel homolog 1 on insect metamorphosis in hemimetabolous and holometabolous species. Sci. Rep. 1, 163. https://doi.org/10.1038/srep00163 (2011).Article
PubMed
PubMed Central
MATH
Google Scholar
Gujar, H. & Palli, S. R. Krüppel homolog 1 and E93 mediate Juvenile hormone regulation of metamorphosis in the common bed bug Cimex lectularius. Sci. Rep. 6, 26092. https://doi.org/10.1038/srep26092 (2016).Article
ADS
PubMed
PubMed Central
Google Scholar
Rauschenbach, I. Y. et al. Insulin-like peptide DILP6 regulates juvenile hormone and dopamine metabolism in Drosophila females. Gen. Compar. Endocrinol. 243, 1–9. https://doi.org/10.1016/j.ygcen.2016.11.004 (2017).Article
Google Scholar
Perez-Hedo, M., Rivera-Perez, C. & Noriega, F. G. Starvation increases insulin sensitivity and reduces juvenile hormone synthesis in mosquitoes. PLoS ONE 9(1), e86183. https://doi.org/10.1371/journal.pone.0086183 (2014).Article
ADS
PubMed
PubMed Central
MATH
Google Scholar
Barthel, A., Schmoll, D. & Unterman, T. G. FoxO proteins in insulin action and metabolism. Trends Endocrinol. Metabol. 16(4), 183–189. https://doi.org/10.1016/j.tem.2005.03.010 (2005).Article
MATH
Google Scholar
Zeng, B. et al. The FOXO transcription factor controls insect growth and development by regulating juvenile hormone degradation in the silkworm, Bombyx mori. J. Biol. Chem. 292(28), 11659–11669. https://doi.org/10.1074/jbc.M117.777797 (2017).Article
PubMed
PubMed Central
Google Scholar
Lin, X., Yu, N. & Smagghe, G. FoxO mediates the timing of pupation through regulating ecdysteroid biosynthesis in the red flour beetle, Tribolium castaneum. Gen. Comp. Endocrinol. 258, 149–156. https://doi.org/10.1016/j.ygcen.2017.05.012 (2018).Article
PubMed
Google Scholar
Süren-Castillo, S., Abrisqueta, M. & Maestro, J. L. FoxO inhibits juvenile hormone biosynthesis and vitellogenin production in the German cockroach. Insect Biochem. Mol. Biol. 42(7), 491–498. https://doi.org/10.1016/j.ibmb.2012.03.006 (2012).Article
PubMed
Google Scholar
Nie, H. et al. Genome-wide identification and characterization of Fox genes in the honeybee, Apis cerana, and comparative analysis with other bee Fox genes. Int. J. Genom. 2018(1), 5702061. https://doi.org/10.1155/2018/5702061 (2018).Article
Google Scholar
Bond, N. D. et al. ssFTZ-F1 and Matrix metalloproteinase 2 are required for fat-body remodeling in Drosophila. Dev. Biol. 360(2), 286–296. https://doi.org/10.1016/j.ydbio.2011.09.015 (2011).Article
PubMed
Google Scholar
Truman, J. W., Rountree, D. B., Reiss, S. E. & Schwartz, L. M. Ecdysteroids regulate the release and action of eclosion hormone in the tobacco hornworm, Manduca sexta (L.). J. Insect. Physiol. 29(12), 895–900. https://doi.org/10.1016/0022-1910(83)90052-5 (1983).Article
Google Scholar
Lozano, J., Montañez, R. & Belles, X. MiR-2 family regulates insect metamorphosis by controlling the juvenile hormone signaling pathway. Proc. Natl. Acad. Sci. 112(12), 3740–3745. https://doi.org/10.1073/pnas.1418522112 (2015).Article
ADS
PubMed
PubMed Central
Google Scholar
He, Q., Zhang, Y. & Dong, W. MicroRNA miR-927 targets the juvenile hormone primary response gene Krüppel homolog1 to control Drosophila developmental growth. Insect Mol. Biol. 29(6), 545–554. https://doi.org/10.1111/imb.12662 (2020).Article
PubMed
Google Scholar
Song, J. et al. The microRNAs let-7 and miR-278 regulate insect metamorphosis and oogenesis by targeting the juvenile hormone early-response gene Krüppel-homolog 1. Development 145(24), dev170670. https://doi.org/10.1242/dev.170670 (2018).Article
PubMed
Google Scholar
Reynolds, J. A. MicroRNAs in the developmental toolbox-a comparative approach to understanding their role in regulating insect development. Curr. Opin. Insect Sci. 101256. https://doi.org/10.1016/j.cois.2024.101256 (2024).Belles, X. MicroRNAs and the evolution of insect metamorphosis. Ann. Rev. Entomol. 62(1), 111–125. https://doi.org/10.1146/annurev-ento-031616-034925 (2017).Article
MATH
Google Scholar
Bashirullah, A. et al. Coordinate regulation of small temporal RNAs at the onset of Drosophila metamorphosis. Dev. Biol. 259(1), 1–8. https://doi.org/10.1016/S0012-1606(03)00063-0 (2003).Article
PubMed
Google Scholar
Sempere, L. F., Sokol, N. S., Dubrovsky, E. B., Berger, E. M. & Ambros, V. Temporal regulation of microRNA expression in Drosophila melanogaster mediated by hormonal signals and broad-Complex gene activity. Dev. Biol. 259(1), 9–18. https://doi.org/10.1016/S0012-1606(03)00208-2 (2003).Article
PubMed
Google Scholar
Asgari, S. MicroRNA functions in insects. Insect Biochem. Mol. Biol. 43(4), 388–397. https://doi.org/10.1016/j.ibmb.2012.10.005 (2013).Article
PubMed
MATH
Google Scholar
Lucas, K. & Raikhel, A. S. Insect microRNAs: Biogenesis, expression profiling and biological functions. Insect Biochem. Mol. Biol. 43(1), 24–38. https://doi.org/10.1016/j.ibmb.2012.10.009 (2013).Article
PubMed
Google Scholar
Download referencesAcknowledgementsThis research was supported by the São Paulo Research Foundation – FAPESP, which provided a fellowship to T.S.D. (Process No. 2016/13854-6) and funded projects No. 2011/03171-5 and No. 2016/06657-0 awarded to Z.L.P.S. We also acknowledge support from CAPES (finance code 001) for institutional resources.Author informationAuthors and AffiliationsDepartment of Genetics, Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, SP, BrazilT. S. Depintor, F. C. P. Freitas, N. Hernandes, F. M. F. Nunes & Z. L. P. SimõesDepartment of Genetics and Evolution, Center for Biological and Health Sciences, Federal University of São Carlos, São Carlos, SP, BrazilF. M. F. NunesDepartment of Biochemistry and Tissue Biology, Institute of Biology, University of Campinas, Campinas, SP, BrazilF. C. P. FreitasDepartment of Biology, Faculty of Philosophy, Sciences and Letters at Ribeirão Preto, University of São Paulo, Ribeirão Preto, SP, BrazilZ. L. P. SimõesAuthorsT. S. DepintorView author publicationsYou can also search for this author inPubMed Google ScholarF. C. P. FreitasView author publicationsYou can also search for this author inPubMed Google ScholarN. HernandesView author publicationsYou can also search for this author inPubMed Google ScholarF. M. F. NunesView author publicationsYou can also search for this author inPubMed Google ScholarZ. L. P. SimõesView author publicationsYou can also search for this author inPubMed Google ScholarContributionsT.S.D., N.H., F.C.P.F., designed the experiments and analyses. T.S.D., F.C.P.F., Z.L.P.S. wrote the paper. Z.L.P.S. provided reagents, materials, and analysis tools. F.C.P.F., F.M.F.N., Z.L.P.S. directed the study. All authors reviewed and approved the manuscript.Corresponding authorCorrespondence to
T. S. Depintor.Ethics declarations
Competing interests
The authors declare no competing interests.
Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationSupplementary Information.Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissionsAbout this articleCite this articleDepintor, T.S., Freitas, F.C.P., Hernandes, N. et al. Interactions of juvenile hormone, 20-hydroxyecdysone, developmental genes, and miRNAs during pupal development in Apis mellifera.
Sci Rep 15, 10354 (2025). https://doi.org/10.1038/s41598-025-93580-7Download citationReceived: 31 October 2024Accepted: 07 March 2025Published: 26 March 2025DOI: https://doi.org/10.1038/s41598-025-93580-7Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy to clipboard
Provided by the Springer Nature SharedIt content-sharing initiative
Keywords
Usp
EcR
Kr-h1
Chd64
InR-2
Met
Ftz-f1
Tai
Metamorphosis regulationmiR-34 and miR-281