Expression of E93 provides an instructive cue to control dynamic enhancer activity and chromatin accessibility during development

Background of Insect Metamorphosis: Numerous Functions of Ecdysteroid

Insect development is mainly controlled by juvenile hormone (JH) and ecdysone, and their hemolymph titer determines the insect direction; larva, pupa, or adult. The mediators of them are Krüppel homolog 1 (Kr-h1), Broad-Complex (BR-C), and E93. They activate genes that characterize larval, pupal, and adult feature, which gives rise to the metamorphosis. Before individual ecdysis, these master factors activate target genes to produce larva, pupa, or adult. Prothoracicotropic hormone (PTTH) from the brain activates prothoracic gland (PG), resulted in the ecdysis. Other factors that activate ecdysteroid production are reported. The produced ecdysteroid was observed early stages of the last larval instar of Bombyx mori, where the ecdysteroid titer is different from previous stage. Two stages are different in JH and ecdysone titer, and the interaction of the JH and ecdysone production determines their titer. Ecdysone brings about the prominent change, which needs many gene transcriptions and the interaction of ecdysone-responsive transcription factors (ERTFs). Their target genes are successively expressed, which brings about the metamorphosis. For the activation of genes, ecdysone gives rise to chromatin remodeling and histone modification. Ecdysone and other factors bring about cell division of the wing disc of the last larval instar; for the proliferation and differentiation, which gives rise to the growth and differentiation of the wing disc for the metamorphosis. In addition, hormone-responsive miRNAs work, several ERTFs function for one gene, and the suppressive TF function along with metamorphosis. Thus, several attractive things underly around the insect metamorphosis. We will be near the understanding of the metamorphosis through these things.

Molecular basis of E93-dependent tissue morphogenesis and histolysis during insect metamorphosis

The evolution of insect metamorphosis has profoundly influenced their successful adaptation and diversification. Two key physiological processes during insect metamorphosis are notable: wing maturation and prothoracic gland (PG) histolysis. The ecdysone-induced protein 93 (E93) is a transcription factor indispensable for metamorphosis. While it has been established that both wing maturation and PG histolysis are dependent on E93, the molecular mechanisms through which E93 regulates these seemingly 'opposing' events remain poorly understood. In this study, time-course transcriptome profiles were generated for wing pads and PGs during metamorphosis in Blattella germanica, a hemimetabolous model insect. Comparative transcriptomic analyses demonstrated that E93 exerts predominant control over extensive gene transcription during wing morphogenesis and PG histolysis. During wing morphogenesis, E93 selectively enhances the expression of genes associated with cell proliferation, energy supply, signal transduction, actin cytoskeleton organization, and cell adhesion, etc. Additionally, E93 activates the transcription of the majority of genes within the wing gene network that are crucial for wing development in B. germanica. During PG histolysis, E93 preferentially promotes the expression of genes related to endocytosis, focal adhesion, the AMPK signaling pathway, adipocytokine signaling pathway, Toll and Imd signaling pathways, and autophagy, etc. The key genes involved in the aforementioned pathways were subsequently confirmed to contribute to the E93-dependent degeneration of the PG in B. germanica. In summary, our results reveal that E93 functions as a master transcriptional regulator orchestrating both tissue morphogenesis and histolysis during insect metamorphosis. These findings contribute to a deeper understanding of the genetic underpinnings of insect metamorphosis.

E93 controls adult differentiation by repressing broad in Drosophila

In Drosophila melanogaster, successful development relies on the precise coordination of both spatial and temporal regulatory axes. The temporal axis governs stage-specific identity and developmental transitions through a number of genes, collectively forming the Metamorphic Gene Network. Among these, Ecdysone inducible protein 93F (E93) serves as the critical determinant for adult specification, but its mechanism of action remains unclear. Here, we found that, rather than acting mainly as an instructive signal, E93 promotes adult differentiation through the repression of the pupal specifier broad (br). In the absence of E93, sustained high levels of Br during the pupal stage strongly represses pupal-specific enhancers that are essential for the terminal differentiation of the wing. Notably, RNA-seq analysis confirmed that the majority of E93-dependent transcriptomic changes in pupal wings are primarily driven by br repression. In addition, we also show that Br represses the pupal-enhancers during the larval and prepupal stages preventing the premature implementation of the adult genetic program, and that it also dampens the activity of larval enhancers during the latter stages of larval development. This mechanism of action seems to be a derived feature acquired in Diptera, as in the coleopteran Tribolium castaneum, repression of br by E93 is not sufficient to allow adult differentiation. In summary, our study elucidates the crucial role of the intricate interplay between E93 and Br as the governing mechanism in the process of terminal differentiation in Drosophila. This finding holds significant implications for advancing our understanding of the evolution of insect metamorphosis.

Chronologically inappropriate morphogenesis (Chinmo) is required for maintenance of larval stages of fall armyworm

Broad complex (Br-C) and eip93F (E93) transcription factors promote insect metamorphosis from larva to pupa and from pupa to adult, respectively. Recently, chronologically inappropriate morphogenesis (Chinmo) has been proposed as a larval specifier in Drosophila melanogaster. However, whether Chinmo is required for larval maintenance in lepidopteran insects, the underlying mechanisms involved in maintaining the larval stage, and its interactions with the JH signaling pathway are not well understood. Here, we used a binary transgenic CRISPR/Cas9 system to knockout Chinmo and Kr-h1 (primary response gene in the JH signaling pathway) in the fall armyworm (FAW). Kr-h1 knockout induced premature metamorphosis only after L5 (penultimate), whereas Chinmo and Kr-h1 double knockout induced premature metamorphosis in L3. Sequencing and differential gene expression (DEG) analysis of RNA isolated from mutants and single-cell multiome ATAC analysis of Chinmo, Kr-h1, and Chinmo and Kr-h1 double knockout Sf9 cells revealed that Chinmo participates in chromatin modifications that prevent the promoter accessibility and expression of metamorphosis promoting genes. These results suggest that Chinmo is a larval specifier that plays a major role in preventing metamorphosis in early larval stages by controlling chromatin accessibility near the promoters of genes such as Br-C and E93 required for pupal and adult development.

Redesigning the Drosophila histone gene cluster: an improved genetic platform for spatiotemporal manipulation of histone function

Mutating replication-dependent (RD) histone genes is an important tool for understanding chromatin-based epigenetic regulation. Deploying this tool in metazoans is particularly challenging because RD histones in these organisms are typically encoded by many genes, often located at multiple loci. Such gene arrangements make the ability to generate homogenous histone mutant genotypes by site-specific gene editing quite difficult. Drosophila melanogaster provides a solution to this problem because the RD histone genes are organized into a single large tandem array that can be deleted and replaced with transgenes containing mutant histone genes. In the last ∼15 years several different RD histone gene replacement platforms were developed using this simple strategy. However, each platform contains weaknesses that preclude full use of the powerful developmental genetic capabilities available to Drosophila researchers. Here we describe the development of a newly engineered platform that rectifies many of these weaknesses. We used CRISPR to precisely delete the RD histone gene array (HisC), replacing it with a multifunctional cassette that permits site-specific insertion of either one or two synthetic gene arrays using selectable markers. We designed this cassette with the ability to selectively delete each of the integrated gene arrays in specific tissues using site-specific recombinases. We also present a method for rapidly synthesizing histone gene arrays of any genotype using Golden Gate cloning technologies. These improvements facilitate the generation of histone mutant cells in various tissues at different stages of Drosophila development and provide an opportunity to apply forward genetic strategies to interrogate chromatin structure and gene regulation.

Redesigning the Drosophila histone gene cluster: An improved genetic platform for spatiotemporal manipulation of histone function

Mutating replication-dependent (RD) histone genes is an important tool for understanding chromatin-based epigenetic regulation. Deploying this tool in metazoan models is particularly challenging because RD histones in these organisms are typically encoded by many genes, often located at multiple loci. Such RD histone gene arrangements make the ability to generate homogenous histone mutant genotypes by site-specific gene editing quite difficult. Drosophila melanogaster provides a solution to this problem because the RD histone genes are organized into a single large tandem array that can be deleted and replaced with transgenes containing mutant histone genes. In the last ∼15 years several different RD histone gene replacement platforms have been developed using this simple strategy. However, each platform contains weaknesses that preclude full use of the powerful developmental genetic capabilities available to Drosophila researchers. Here we describe the development of a newly engineered platform that rectifies many of these weaknesses. We used CRISPR to precisely delete the RD histone gene array ( HisC ), replacing it with a multifunctional cassette that permits site-specific insertion of either one or two synthetic gene arrays using selectable markers. We designed this cassette with the ability to selectively delete each of the integrated gene arrays in specific tissues using site-specific recombinases. We also present a method for rapidly synthesizing histone gene arrays of any genotype using Golden Gate cloning technologies. These improvements facilitate generation of histone mutant cells in various tissues at different stages of Drosophila development and provide an opportunity to apply forward genetic strategies to interrogate chromatin structure and gene regulation.

The SWI/SNF nucleosome remodeler constrains enhancer activity during Drosophila wing development

Chromatin remodeling is central to the dynamic changes in gene expression that drive cell fate determination. During development, the sets of enhancers that are accessible for use change globally as cells transition between stages. While transcription factors and nucleosome remodelers are known to work together to control enhancer accessibility, it is unclear how the short stretches of DNA that they individually unmask yield the kilobase-sized accessible regions characteristic of active enhancers. Here, we performed a genetic screen to investigate the role of nucleosome remodelers in control of dynamic enhancer activity. We find that the Drosophila Switch/Sucrose Non-Fermenting complex, BAP, is required for repression of a temporally dynamic enhancer, brdisc. Contrary to expectations, we find that the BAP-specific subunit Osa is dispensable for mediating changes in chromatin accessibility between the early and late stages of wing development. Instead, we find that Osa is required to constrain the levels of brdisc activity when the enhancer is normally active. Genome-wide profiling reveals that Osa directly binds brdisc as well as thousands of other developmentally dynamic regulatory sites, including multiple genes encoding components and targets of the Notch signaling pathway. Transgenic reporter analyses demonstrate that Osa is required for activation and for constraint of different sets of target enhancers in the same cells. Moreover, Osa loss results in hyperactivation of the Notch ligand Delta and development of ectopic sensory structures patterned by Notch signaling early in development. Together, these findings indicate that proper constraint of enhancer activity is necessary for regulation of dose-dependent developmental events.

Transcriptomic exploration of the Coleopteran wings reveals insight into the evolution of novel structures associated with the beetle elytron

The acquisition of novel traits is central to organismal evolution, yet the molecular mechanisms underlying this process are elusive. The beetle forewings (elytra) are evolutionarily modified to serve as a protective shield, providing a unique opportunity to study these mechanisms. In the past, the orthologs of genes within the wing gene network from Drosophila studies served as the starting point when studying the evolution of elytra (candidate genes). Although effective, candidate gene lists are finite and only explore genes conserved across species. To go beyond candidate genes, we used RNA sequencing and explored the wing transcriptomes of two Coleopteran species, the red flour beetle (Tribolium castaneum) and the Japanese stag beetle (Dorcus hopei). Our analysis revealed sets of genes enriched in Tribolium elytra (57 genes) and genes unique to the hindwings, which possess more "typical" insect wing morphologies (29 genes). Over a third of the hindwing-enriched genes were "candidate genes" whose functions were previously analyzed in Tribolium, demonstrating the robustness of our sequencing. Although the overlap was limited, transcriptomic comparison between the beetle species found a common set of genes, including key wing genes, enriched in either elytra or hindwings. Our RNA interference analysis for elytron-enriched genes in Tribolium uncovered novel genes with roles in forming various aspects of morphology that are unique to elytra, such as pigmentation, hardening, sensory development, and vein formation. Our analyses deepen our understanding of how gene network evolution facilitated the emergence of the elytron, a unique structure critical to the evolutionary success of beetles.

Increased chromatin accessibility promotes the evolution of a transcriptional silencer in Drosophila

The loss of discrete morphological traits, the most common evolutionary transition, is typically driven by changes in developmental gene expression. Mutations accumulating in regulatory elements of these genes can disrupt DNA binding sites for transcription factors patterning their spatial expression, or delete entire enhancers. Regulatory elements, however, may be silenced through changes in chromatin accessibility or the emergence of repressive elements. Here, we show that increased chromatin accessibility at the gene yellow, combined with the gain of a repressor site, underlies the loss of a wing spot pigmentation pattern in a Drosophila species. The gain of accessibility of this repressive element is regulated by E93, a transcription factor governing the progress of metamorphosis. This convoluted evolutionary scenario contrasts with the parsimonious mutational paths generally envisioned and often documented for morphological losses. It illustrates how evolutionary changes in chromatin accessibility may directly contribute to morphological diversification.

Reduction of embryonic E93 expression as a hypothetical driver of the evolution of insect metamorphosis

The early embryo of the cockroach Blattella germanica exhibits high E93 expression. In general, E93 triggers adult morphogenesis during postembryonic development. Here we show that E93 is also crucial in early embryogenesis in the cockroach, as a significant number of E93-depleted embryos are unable to develop the germ band under maternal RNAi treatment targeting E93. Moreover, transcriptomic analysis indicates that E93 depletion results in important gene expression changes in the early embryo, and many of the differentially expressed genes are involved in development. Then, using public databases, we gathered E93 expression data in embryo and preadult stages, finding that embryonic expression of E93 is high in hemimetabolan species (whose juveniles, or nymphs, are similar to the adult) and low in holometabolans (whose juveniles, or larvae, are different from the adult). E93 expression is also low in Thysanoptera and in Hemiptera Sternorrhyncha, hemimetabolans with postembryonic quiescent stages, as well as in Odonata, the nymph of which is very different from the adult. In ametabolans, such as the Zygentoma Thermobia domestica, E93 transcript levels are very high in the early embryo, whereas during postembryonic development they are medium and relatively constant. We propose the hypothesis that during evolution, a reduction of E93 expression in the embryo of hemimetabolans facilitated the larval development and the emergence of holometaboly. Independent decreases of E93 transcripts in the embryo of Odonata, Thysanoptera, and different groups of Hemiptera Sternorrhyncha would have allowed the development of modified juvenile stages adapted to specific ecophysiological conditions.

CRISPR-Cas9 Genome Editing Uncovers the Mode of Action of Methoprene in the Yellow Fever Mosquito, Aedes aegypti

Methoprene, a juvenile hormone (JH) analog, is widely used for insect control, but its mode of action is not known. To study methoprene action in the yellow fever mosquito, Aedes aegypti, the E93 (ecdysone-induced transcription factor) was knocked out using the CRISPR-Cas9 system. The E93 mutant pupae retained larval tissues similar to methoprene-treated insects. These insects completed pupal ecdysis and died as pupa. In addition, the expression of transcription factors, broad complex and Krüppel homolog 1 (Kr-h1), increased and that of programmed cell death (PCD) and autophagy genes decreased in E93 mutants. These data suggest that methoprene functions through JH receptor, methoprene-tolerant, and induces the expression of Kr-h1, which suppresses the expression of E93, resulting in a block in PCD and autophagy of larval tissues. Failure in the elimination of larval tissues and the formation of adult structures results in their death. These results answered long-standing questions on the mode of action of methoprene.

Knockdown of Ecdysone-Induced Protein 93F Causes Abnormal Pupae and Adults in the Eggplant Lady Beetle

Ecdysone-induced protein 93F (E93) plays triple roles during post-embryonic development in insects whose juvenile instars are more than four. However, it only acts as a specifier of adult structures in Drosophila flies whose larval instars are fixed at three. In this study, we determined the functions of E93 in the eggplant lady beetle (Henosepilachna vigintioctopunctata), which has four larval instars. We uncovered that E93 was abundantly expressed at the prepupal and pupal stages. A precocious inhibition of the juvenile hormone signal by RNA interference (RNAi) of HvKr-h1 or HvHairy, two vital downstream developmental effectors, at the penultimate instar larval stage increased the expression of E93, Conversely, ingestion of JH by the third-instar larvae stimulated the expression of HvKr-h1 but repressed the transcription of either HvE93X1 or HvE93X2. However, disturbance of the JH signal neither drove premature metamorphosis nor caused supernumerary instars. In contrast, depletion of E93 at the third- and fourth-instar larval and prepupal stages severely impaired pupation and caused a larval-pupal mixed phenotype: pupal spines and larval scoli were simultaneously presented on the cuticle. RNAi of E93 at the pupal stage affected adult eclosion. When the beetles had suffered from a dsE93 injection at the fourth-instar larval and pupal stages, a few resultant adults emerged, with separated elytra, abnormally folded hindwings, a small body size and short appendages. Taken together, our results suggest the larval instars are fixed in H. vigintioctopunctata; E93 serves as a repressor of larval characters and a specifier of adult structures during the larval-pupal-adult transition.

Opportunistic binding of EcR to open chromatin drives tissue-specific developmental responses

Steroid hormones perform diverse biological functions in developing and adult animals. However, the mechanistic basis for their tissue specificity remains unclear. In Drosophila, the ecdysone steroid hormone is essential for coordinating developmental timing across physically separated tissues. Ecdysone directly impacts genome function through its nuclear receptor, a heterodimer of the EcR and ultraspiracle proteins. Ligand binding to EcR triggers a transcriptional cascade, including activation of a set of primary response transcription factors. The hierarchical organization of this pathway has left the direct role of EcR in mediating ecdysone responses obscured. Here, we investigate the role of EcR in controlling tissue-specific ecdysone responses, focusing on two tissues that diverge in their response to rising ecdysone titers: the larval salivary gland, which undergoes programmed destruction, and the wing imaginal disc, which initiates morphogenesis. We find that EcR functions bimodally, with both gene repressive and activating functions, even at the same developmental stage. EcR DNA binding profiles are highly tissue-specific, and transgenic reporter analyses demonstrate that EcR plays a direct role in controlling enhancer activity. Finally, despite a strong correlation between tissue-specific EcR binding and tissue-specific open chromatin, we find that EcR does not control chromatin accessibility at genomic targets. We conclude that EcR contributes extensively to tissue-specific ecdysone responses. However, control over access to its binding sites is subordinated to other transcription factors.

Chinmo is the larval member of the molecular trinity that directs Drosophila metamorphosis

The genome of insects with complete metamorphosis contains the instructions for making three distinct body forms, that of the larva, of the pupa, and of the adult. However, the molecular mechanisms by which each gene set is called forth and stably expressed are poorly understood. A half century ago, it was proposed that there was a set of three master genes that inhibited each other’s expression and enabled the expression of genes for each respective stage. We show that the transcription factor chinmo is essential for maintaining the larval stage in Drosophila, and with two other regulatory genes, broad and E93, makes up the trinity of mutually repressive master genes that underlie insect metamorphosis.

The molecular control of insect metamorphosis from larva to pupa to adult has long been a mystery. The Broad and E93 transcription factors, which can modify chromatin domains, are known to direct the production of the pupa and the adult, respectively. We now show that chinmo, a gene related to broad, is essential for the repression of these metamorphic genes. Chinmo is strongly expressed during the formation and growth of the larva and its removal results in the precocious expression of broad and E93 in the first stage larva, causing a shift from larval to premetamorphic functions. This trinity of Chinmo, Broad, and E93 regulatory factors is mutually inhibitory. The interaction of this network with regulatory hormones likely ensures the orderly progression through insect metamorphosis.

Drosophila E93 promotes adult development and suppresses larval responses to ecdysone during metamorphosis

Pulses of the steroid hormone ecdysone act through transcriptional cascades to direct the major developmental transitions during the Drosophila life cycle. These include the prepupal ecdysone pulse, which occurs 10 ​hours after pupariation and triggers the onset of adult morphogenesis and larval tissue destruction. E93 encodes a transcription factor that is specifically induced by the prepupal pulse of ecdysone, supporting a model proposed by earlier work that it specifies the onset of adult development. Although a number of studies have addressed these functions for E93, little is known about its roles in the salivary gland where the E93 locus was originally identified. Here we show that E93 is required for development through late pupal stages, with mutants displaying defects in adult differentiation and no detectable effect on the destruction of larval salivary glands. RNA-seq analysis demonstrates that E93 regulates genes involved in development and morphogenesis in the salivary glands, but has little effect on cell death gene expression. We also show that E93 is required to direct the proper timing of ecdysone-regulated gene expression in salivary glands, and that it suppresses earlier transcriptional programs that occur during larval and prepupal stages. These studies support the model that the stage-specific induction of E93 in late prepupae provides a critical signal that defines the end of larval development and the onset of adult differentiation.

Memes: A motif analysis environment in R using tools from the MEME Suite

Identification of biopolymer motifs represents a key step in the analysis of biological sequences. The MEME Suite is a widely used toolkit for comprehensive analysis of biopolymer motifs; however, these tools are poorly integrated within popular analysis frameworks like the R/Bioconductor project, creating barriers to their use. Here we present memes, an R package that provides a seamless R interface to a selection of popular MEME Suite tools. memes provides a novel “data aware” interface to these tools, enabling rapid and complex discriminative motif analysis workflows. In addition to interfacing with popular MEME Suite tools, memes leverages existing R/Bioconductor data structures to store the multidimensional data returned by MEME Suite tools for rapid data access and manipulation. Finally, memes provides data visualization capabilities to facilitate communication of results. memes is available as a Bioconductor package at https://bioconductor.org/packages/memes, and the source code can be found at github.com/snystrom/memes.

Biologically active molecules such as DNA, RNA, and proteins are polymers composed of repeated subunits. For example, nucleotides are the subunits of DNA and RNA, and amino acids are the subunits of proteins. Functional properties of biopolymers are determined by short, recurring stretches of subunits known as motifs. Motifs can serve as binding sites between molecules, they can influence the structure of molecules, and they can contribute to enzymatic activities. For these reasons, motif analysis has become a key step in determining the function of biopolymers and in elucidating their roles in biological networks. The MEME Suite is a widely used compilation of tools used for identifying and analyzing motifs found in biological sequences. Here, we describe a new piece of software named “memes” that connects MEME Suite tools to R, the statistical analysis environment. By providing an interface between the MEME Suite and R, memes allows for improved motif analysis workflows and easy access to a wide range of existing data visualization tools, further expanding the utility of MEME Suite tools.

The ecdysone-induced protein 93 is a key factor regulating gonadotrophic cycles in the adult female mosquito Aedes aegypti

Repeated blood feedings are required for adult female mosquitoes to maintain their gonadotrophic cycles, enabling them to be important pathogen carriers of human diseases. Elucidating the molecular mechanism underlying developmental switches between these mosquito gonadotrophic cycles will provide valuable insight into mosquito reproduction and could aid in the identification of targets to disrupt these cycles, thereby reducing disease transmission. We report here that the transcription factor ecdysone-induced protein 93 (E93), previously implicated in insect metamorphic transitions, plays a key role in determining the gonadotrophic cyclicity in adult females of the major arboviral vector Aedes aegypti Expression of the E93 gene in mosquitoes is down-regulated by juvenile hormone (JH) and up-regulated by 20-hydroxyecdysone (20E). We find that E93 controls Hormone Receptor 3 (HR3), the transcription factor linked to the termination of reproductive cycles. Moreover, knockdown of E93 expression via RNAi impaired fat body autophagy, suggesting that E93 governs autophagy-induced termination of vitellogenesis. E93 RNAi silencing prior to the first gonadotrophic cycle affected normal progression of the second cycle. Finally, transcriptomic analysis showed a considerable E93-dependent decline in the expression of genes involved in translation and metabolism at the end of a reproductive cycle. In conclusion, our data demonstrate that E93 acts as a crucial factor in regulating reproductive cycle switches in adult female mosquitoes.

Mechanisms underlying the control of dynamic regulatory element activity and chromatin accessibility during metamorphosis

Cis-regulatory modules of metazoan genomes determine the when and where of gene expression during development. Here we discuss insights into the genetic and molecular mechanisms behind cis-regulatory module usage that have come from recent application of genomics assays to insect metamorphosis. Assays including FAIRE-seq, ATAC-seq, and CUT&RUN indicate that sequential changes in chromatin accessibility play a key role in mediating stage-specific cis-regulatory module activity and gene expression. We review the current understanding of what controls precisely coordinated changes in chromatin accessibility during metamorphosis and describe evidence that points to systemic hormone signaling as a primary signal to trigger genome-wide shifts in accessibility patterns and cis-regulatory module usage.