A comparative genomic analysis of targets of Hox protein Ultrabithorax amongst distant insect species

Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

A Micro-evolutionary Change in Target Binding Sites as a Key Determinant of Ultrabithorax Function in Drosophila

Hox genes encode Homeodomain-containing transcription factors, which specify segmental identities along the anterior-posterior axis. Functional changes in Hox genes have been directly implicated in the evolution of body plans across the metazoan lineage. The Hox protein Ultrabithorax (Ubx) is expressed and required in developing third thoracic (T3) segments in holometabolous insects studied so far, particularly, of the order Coleoptera, Lepidoptera and Diptera. Ubx function is key to specify differential development of the second (T2) and T3 thoracic segments in these insects. While Ubx is expressed in the third thoracic segment in developing larvae of Hymenopteran Apis mellifera, the morphological differences between T2 and T3 are subtle. To identify evolutionary changes that are behind the differential function of Ubx in Drosophila and Apis, which are diverged for more than 350 million years, we performed comparative analyses of genome wide Ubx-binding sites between these two insects. Our studies reveal that a motif with a TAAAT core is a preferred binding site for Ubx in Drosophila, but not in Apis. Biochemical and transgenic assays suggest that in Drosophila, the TAAAT core sequence in the Ubx binding sites is required for Ubx-mediated regulation of two of its target genes studied here; CG13222, a gene that is normally upregulated by Ubx and vestigial (vg), whose expression is repressed by Ubx in T3. Interestingly, changing the TAAT site to a TAAAT site was sufficient to bring an otherwise unresponsive enhancer of the vg gene from Apis under the control of Ubx in a Drosophila transgenic assay. Taken together, our results suggest an evolutionary mechanism by which critical wing patterning genes might have come under the regulation of Ubx in the Dipteran lineage.

A Comparative Study of Transcriptional Regulation Mechanism of Cytochrome P450 CYP6B7 between Resistant and Susceptible Strains of Helicoverpa armigera

Cytochrome P450 CYP6B7 has previously been proved to be associated with fenvalerate-resistance in Helicoverpa armigera. Here, how CYP6B7 is regulated and involved in the resistance of H. armigera is studied. Seven base differences (M1-M7) were found in CYP6B7 promoter between a fenvalerate-resistant (HDTJFR) and a susceptible (HDTJ) strain of H. armigera. M1-M7 sites in HDTJFR were mutated into the corresponding base in HDTJ, and pGL3-CYP6B7 reporter genes with different mutation sites were constructed. Fenvalerate-induced activities of reporter genes mutated at M3, M4, and M7 sites were significantly reduced. Transcription factors Ubx and Br, whose binding sites contain M3 and M7, respectively, were overexpressed in HDTJFR. Knockdown of Ubx and Br results in significant expression inhibition of CYP6B7 and other resistance-related P450 genes, and increase of sensitivity of H. armigera to fenvalerate. These results indicate that Ubx and Br regulate the expression of CYP6B7 to mediate the fenvalerate-resistance in H. armigera.

Dynamic chromatin conformation and accessibility changes mediate the spatial-specific gene regulatory network in Bombyx mori

Silk gland genes of Bombyx mori can have strict spatial expression patterns, which impact their functions and silk quality; however, our understanding of their regulation mechanisms is currently insufficient. To address this, the middle silk gland (MSG) and posterior silk gland (PSG) of the silkworm were investigated. Gene ontology annotation showed that spatially specific expressed genes were involved in the formation of H3k9me and chromatin topology. Chromatin conformation data generated by Hi-C showed that the topologically associated domain boundaries around FibL and Sericin1 genes were significantly different between MSG and PSG. Changes in chromatin conformation led to changes in chromatin activity, which significantly affected the expression of nearby genes in silkworm. Chromatin accessibility regions of MSG and PSG were analyzed using FAIRE-seq, and 1006 transcription factor motifs were identified in open chromatin regions. Furthermore, the spatial-specific expression patterns of silk gland genes were mainly associated with homeobox-contained transcription factors, such as POU-M2, which was specifically bound and relatively highly expressed in the MSG. The regulatory network mediated by POU-M2 regulated most of the spatial-specific expressed genes in MSG, such as ADH1. These results can aid in improving silk performance, optimizing silkworm breeding, and improving the gene spatial regulatory model research for insects.

Overexpression of Ultrabithorax Changes the Development of Silk Gland and the Expression of Fibroin Genes in Bombyx mori

Ultrabithorax (Ubx) is a member of the Hox gene group involved in cell fate decisions, cell proliferation and organ identity. Its function has been extensively researched in Drosophila melanogaster but little is known about it in Lepidoptera. To uncover the function of Ubx in the development of lepidopterans, we constructed the Ubx overexpression (Ubx^OE) strain based on the Nistari strain of Bombyx mori. The Ubx^OE strain showed a small body size, transparent intersegmental membrane and abnormal posterior silk gland (PSG). In the current study, we focused on the effect of Ubx overexpression on the posterior silk gland. As the major protein product of PSG, the mRNA expression of fibroin heavy chain (Fib-H) and fibroin light chain (Fib-L) was upregulated three times in Ubx^OE, but the protein expression of Fib-H and Fib-L was not significantly different. We speculated that the overexpression of Ubx downregulated the expression of Myc and further caused abnormal synthesis of the spliceosome and ribosome. Abnormalities of the spliceosome and ribosome affected the synthesis of protein in the PSG and changed its morphology.

Hox dosage and morphological diversification during development and evolution

Hox genes encode for evolutionary conserved transcription factors that have long fascinated biologists since the observation of the first homeotic transformations in flies. Hox genes are developmental architects that instruct the formation of various and precise morphologies along the body axes in cnidarian and bilaterian species. In contrast to these highly specific developmental functions, Hox genes encode for proteins that display poorly selective DNA-binding properties in vitro. This "Hox paradox" has been partially solved with the discovery of the TALE-class cofactors, which interact with all Hox members and form versatile Hox/TALE protein complexes on DNA. Here, we describe the role of the Hox dosage as an additional molecular strategy contributing to further resolve the Hox paradox. We present several cases where the Hox dosage is involved in the formation of different morphologies in invertebrates and vertebrates, with a particular emphasis on flight appendages in insects. We also discuss how the Hox dosage could be interpreted in different types of target enhancers within the nuclear environment in vivo. Altogether our survey underlines the Hox dosage as a key mechanism for shaping Hox molecular function during development and evolution.

CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis

Simple Summary

Homeotic genes encode transcription factors that coordinated the anatomical structure formation during the early embryonic development of organisms. In this study, we functionally characterized two homeotic genes, Abdominal-A (Abd-A) and Ultrabithorax (Ubx), in the Asian corn borer, Ostrinia furnacalis (a maize pest that has devastated the Asia-Pacific region) by using a CRISPR/Cas9 genome editing system. Our results show that the mutagenesis of OfAbd-A and OfUbx led to severe morphological defects in O. furnacalis, which included fused segments and segmental twist during the larval stage, and hollowed and incision-like segments during the pupal stage in OfAbd-A mutants, as well as defects in the wing-pad development in pupal and adult OfUbx mutants. Overall, knocking out Abd-A and Ubx in O. furnacalis resulted in the embryonic lethality to, and pleiotropic impact on, other homeotic genes. This study not only confirms the conserved body planning functions in OfAbd-A and OfUbx, but it also strengthens the control implications of these homeotic genes for lepidopteran pests.

Abstract

(1) Background: Abdominal-A (Abd-A) and Ultrabithorax (Ubx) are homeotic genes that determine the identity and morphology of the thorax and abdomen in insects. The Asian corn borer, Ostrinia furnacalis (Guenée) (Lepidoptera: Pyralidae), is a devastating maize pest throughout Asia, the Western Pacific, and Australia. Building on previous knowledge, we hypothesized that the knockout of Abd-A and Ubx would disrupt the abdominal body planning in O. furnacalis. (2) Methods: CRISPR/Cas9-targeted mutagenesis was employed to decipher the functions of these homeotic genes. (3) Results: Knockout insects demonstrated classical homeotic transformations. Specifically, the mutagenesis of OfAbd-A resulted in: (1) Fused segments and segmental twist during the larval stage; (2) Embryonic lethality; and (3) The pleiotropic upregulation of other homeotic genes, including Lab, Pd, Dfd, Antp, and Abd-B. The mutagenesis of OfUbx led to: (1) Severe defects in the wing pads, which limited the ability of the adults to fly and mate; (2) Female sterility; and (3) The pleiotropic upregulation of other homeotic genes, including Dfd, Abd-B, and Wnt1. (4) Conclusions: These combined results not only support our hypothesis, but they also strengthen the potential of using homeotic genes as molecular targets for the genetic control of this global insect pest.

Roles of Drosophila Hox Genes in the Assembly of Neuromuscular Networks and Behavior

Hox genes have been known for specifying the anterior-posterior axis (AP) in bilaterian body plans. Studies in vertebrates have shown their importance in developing region-specific neural circuitry and diversifying motor neuron pools. In Drosophila, they are instrumental for segment-specific neurogenesis and myogenesis early in development. Their robust expression in differentiated neurons implied their role in assembling region-specific neuromuscular networks. In the last decade, studies in Drosophila have unequivocally established that Hox genes go beyond their conventional functions of generating cellular diversity along the AP axis of the developing central nervous system. These roles range from establishing and maintaining the neuromuscular networks to controlling their function by regulating the motor neuron morphology and neurophysiology, thereby directly impacting the behavior. Here we summarize the limited knowledge on the role of Drosophila Hox genes in the assembly of region-specific neuromuscular networks and their effect on associated behavior.

Transcriptome dynamics during metamorphosis of imaginal discs into wings and thoracic dorsum in Apis mellifera castes

BACKGROUND

Much of the complex anatomy of a holometabolous insect is built from disc-shaped epithelial structures found inside the larva, i.e., the imaginal discs, which undergo a rapid differentiation during metamorphosis. Imaginal discs-derived structures, like wings, are built through the action of genes under precise regulation.

RESULTS

We analyzed 30 honeybee transcriptomes in the search for the gene expression needed for wings and thoracic dorsum construction from the larval wing discs primordia. Analyses were carried out before, during, and after the metamorphic molt and using worker and queen castes. Our RNA-seq libraries revealed 13,202 genes, representing 86.2% of the honeybee annotated genes. Gene Ontology analysis revealed functional terms that were caste-specific or shared by workers and queens. Genes expressed in wing discs and descendant structures showed differential expression profiles dynamics in premetamorphic, metamorphic and postmetamorphic developmental phases, and also between castes. At the metamorphic molt, when ecdysteroids peak, the wing buds of workers showed maximal gene upregulation comparatively to queens, thus underscoring differences in gene expression between castes at the height of the larval-pupal transition. Analysis of small RNA libraries of wing buds allowed us to build miRNA-mRNA interaction networks to predict the regulation of genes expressed during wing discs development.

CONCLUSION

Together, these data reveal gene expression dynamics leading to wings and thoracic dorsum formation from the wing discs, besides highlighting caste-specific differences during wing discs metamorphosis.

Hox genes are essential for the development of eyespots in Bicyclus anynana butterflies

The eyespot patterns found on the wings of nymphalid butterflies are novel traits that originated first in hindwings and subsequently in forewings, suggesting that eyespot development might be dependent on Hox genes. Hindwings differ from forewings in the expression of Ultrabithorax (Ubx), but the function of this Hox gene in eyespot development as well as that of another Hox gene Antennapedia (Antp), expressed specifically in eyespots centers on both wings, are still unclear. We used CRISPR-Cas9 to target both genes in Bicyclus anynana butterflies. We show that Antp is essential for eyespot development on the forewings and for the differentiation of white centers and larger eyespots on hindwings, whereas Ubx is essential not only for the development of at least some hindwing eyespots but also for repressing the size of other eyespots. Additionally, Antp is essential for the development of silver scales in male wings. In summary, Antp and Ubx, in addition to their conserved roles in modifying serially homologous segments along the anterior-posterior axis of insects, have acquired a novel role in promoting the development of a new set of serial homologs, the eyespot patterns, in both forewings (Antp) and hindwings (Antp and Ubx) of B. anynana butterflies. We propose that the peculiar pattern of eyespot origins on hindwings first, followed by forewings, could be due to an initial co-option of Ubx into eyespot development followed by a later, partially redundant, co-option of Antp into the same network.

Hox dosage contributes to flight appendage morphology in Drosophila

Flying insects have invaded all the aerial space on Earth and this astonishing radiation could not have been possible without a remarkable morphological diversification of their flight appendages. Here, we show that characteristic spatial expression profiles and levels of the Hox genes Antennapedia (Antp) and Ultrabithorax (Ubx) underlie the formation of two different flight organs in the fruit fly Drosophila melanogaster. We further demonstrate that flight appendage morphology is dependent on specific Hox doses. Interestingly, we find that wing morphology from evolutionary distant four-winged insect species is also associated with a differential expression of Antp and Ubx. We propose that variation in the spatial expression profile and dosage of Hox proteins is a major determinant of flight appendage diversification in Drosophila and possibly in other insect species during evolution.

Ultrabithorax Is a Micromanager of Hindwing Identity in Butterflies and Moths

The forewings and hindwings of butterflies and moths (Lepidoptera) are differentiated from each other, with segment-specific morphologies and color patterns that mediate a wide range of functions in flight, signaling, and protection. The Hox geneUltrabithorax(Ubx) is a master selector gene that differentiates metathoracic from mesothoracic identities across winged insects, and previous work has shown this role extends to at least some of the color patterns from the butterfly hindwing. Here we used CRISPR targeted mutagenesis to generateUbxloss-of-function somatic mutations in two nymphalid butterflies (Junonia coenia,Vanessa cardui) and a pyralid moth (Plodia interpunctella). The resulting mosaic clones yielded hindwing-to-forewing transformations, showingUbxis necessary for specifying many aspects of hindwing-specific identities, including scale morphologies, color patterns, and wing venation and structure. These homeotic phenotypes showed cell-autonomous, sharp transitions between mutant and non-mutant scales, except for clones that encroached into the border ocelli (eyespots) and resulted in composite and non-autonomous effects on eyespot ring determination. In the pyralid moth, homeotic clones converted the folding and depigmented hindwing into rigid and pigmented composites, affected the wing-coupling frenulum, and induced ectopic scent-scales in male androconia. These data confirmUbxis a master selector of lepidopteran hindwing identity and suggest it acts on many gene regulatory networks involved in wing development and patterning.

Ultrabithorax is a key regulator for the dimorphism of wings, a main cause for the outbreak of planthoppers in rice

Rice planthoppers, the most devastating rice pests, occur in two wing forms: the short-wing form for rapid population growth and long-wing form for long-distance migration, which together create the mechanism for outbreak. Here we show that Ultrabithorax (Ubx) is a key regulator for switching between the long- and short-wing forms of rice planthoppers. Ubx is expressed in both forewing and hindwing pads, which is different from the canonical model of Ubx expression. In brown planthoppers, expression of Ubx (NlUbx) is regulated by nutritional status of the rice host. High-quality young plants induce NlUbx expression leading to the short-wing form; low-quality ripe plants reduce NlUbx expression resulting in long-wing form. We also showed that NlUbx is regulated by the insulin receptors NlInR1 and NlInR2. The default expression of NlInR1 inhibits NlUbx resulting in long-wings, while high-quality hosts induce NlInR2 expression, which represses NlInR1 thus promoting NlUbx expression to produce short-wings.

How to study enhancers in non-traditional insect models

Transcriptional enhancers are central to the function and evolution of genes and gene regulation. At the organismal level, enhancers play a crucial role in coordinating tissue- and context-dependent gene expression. At the population level, changes in enhancers are thought to be a major driving force that facilitates evolution of diverse traits. An amazing array of diverse traits seen in insect morphology, physiology and behavior has been the subject of research for centuries. Although enhancer studies in insects outside of Drosophila have been limited, recent advances in functional genomic approaches have begun to make such studies possible in an increasing selection of insect species. Here, instead of comprehensively reviewing currently available technologies for enhancer studies in established model organisms such as Drosophila, we focus on a subset of computational and experimental approaches that are likely applicable to non-Drosophila insects, and discuss the pros and cons of each approach. We discuss the importance of validating enhancer function and evaluate several possible validation methods, such as reporter assays and genome editing. Key points and potential pitfalls when establishing a reporter assay system in non-traditional insect models are also discussed. We close with a discussion of how to advance enhancer studies in insects, both by improving computational approaches and by expanding the genetic toolbox in various insects. Through these discussions, this Review provides a conceptual framework for studying the function and evolution of enhancers in non-traditional insect models.

Synergism of open chromatin regions involved in regulating genes in Bombyx mori

The dynamic variability of transcription factors (TFs) and their binding sites makes it challenging to conduct genome-wide transcription regulation research. The silkworm Bombyx mori, which produces silk, is one of the most valuable model insects in the order Lepidoptera. The "opening" and "closing" of chromatin in different silk yield strains is associated with changes in silk production, making this insect a good model for studying the transcriptional regulation of genes. However, few studies have examined the open chromatin regions (OCRs) of silkworms, and studying OCR synergism and their function in silk production remains challenging. Here, we performed formaldehyde-assisted isolation of regulatory elements (FAIRE) to isolate OCRs from the silk glands of fifth-instar larvae of the DaZao and D872 strains. In total, 128,908 high confidence OCRs were identified and approximately 80% of OCRs were located in non-coding regions. OCRs upregulated adjacent genes and showed signal-dependent vulnerability to single-nucleotide polymorphisms. Mid- and low-signal OCRs were more likely to have single-nucleotide polymorphisms (SNP). Further, OCRs interacted with each other within a distance of 5 kb. We named the OCR interaction complex as the "cluster of related regions" (COREs). The functions of the CORE and its harbored OCRs showed some differences. Additionally, COREs enriched many silk protein synthesis-associated genes, some of which were upregulated. This study identified numerous high confidence regulation sites and synergistic regulatory modes of OCRs that affect adjacent genes. These results provide insight into silkworm transcriptional regulation and improve our understanding of cis-element cooperation.

The cis-regulatory logic underlying abdominal Hox-mediated repression versus activation of regulatory elements in Drosophila

During development diverse transcription factor inputs are integrated by cis-regulatory modules (CRMs) to yield cell-specific gene expression. Defining how CRMs recruit the appropriate combinations of factors to either activate or repress gene expression remains a challenge. In this study, we compare and contrast the ability of two CRMs within the Drosophila embryo to recruit functional Hox transcription factor complexes. The DCRE CRM recruits Ultrabithorax (Ubx) and Abdominal-A (Abd-A) Hox complexes that include the Extradenticle (Exd) and Homothorax (Hth) transcription factors to repress the Distal-less leg selector gene, whereas the RhoA CRM selectively recruits Abd-A/Exd/Hth complexes to activate rhomboid and stimulate Epidermal Growth Factor secretion in sensory cell precursors. By swapping binding sites between these elements, we found that the RhoA Exd/Hth/Hox site configuration that mediates Abd-A specific activation can convey transcriptional repression by both Ubx and Abd-A when placed into the DCRE. We further show that the orientation and spacing of Hox sites relative to additional binding sites within the RhoA and DCRE is critical to mediate cell- and segment-specific output. These results indicate that the configuration of Exd, Hth, and Hox site within RhoA is neither Abd-A specific nor activation specific. Instead Hox specific output is largely dependent upon the presence of appropriately spaced and oriented binding sites for additional TF inputs. Taken together, these studies provide insight into the cis-regulatory logic used to generate cell-specific outputs via recruiting Hox transcription factor complexes.

Genome-wide open chromatin regions and their effects on the regulation of silk protein genes in Bombyx mori

Nucleosome-depleted open chromatin regions (OCRs) often harbor transcription factor (TF) binding sites that are associated with active DNA regulatory elements. To investigate the regulation of silk-protein genes, DNA molecules isolated from the silk glands of third-day fifth-instar silkworm larvae and embryo-derived (BmE) cells were subjected to formal dehyde-assisted isolation of regulatory elements (FAIRE) and high-throughput sequencing. In total, 68,000 OCRs were identified, and a number of TF-binding motifs were predicted. In particular, OCRs located near silk-protein genes contained potential binding sites for functional TFs. Moreover, many TFs were found to bind to clusters of OCRs upstream of silk-protein genes, and to regulate the expression of these genes. The expression of silk protein genes may be related not only to regulating TFs (such as fkh, Bmdimm, and Bmsage), but also to developmental and hormone-induced TFs (such as zen, eve, Br, and eip74ef). Elucidation of genome-wide OCRs and their regulatory motifs in silk protein genes will provide valuable data and clues for characterizing the mechanisms of transcriptional control of silk protein genes.

Ultrabithorax and the evolution of insect forewing/hindwing differentiation

Decades have passed since the stunning four-winged phenotype of the Drosophila Ultrabithorax (Ubx) mutant was reported, and accumulating knowledge obtained from studies on Ubx in Drosophila has provided a framework to investigate the role of Ubx during insect wing evolution. With several new insights emerging from recent studies in non-Drosophila insects, along with the outcomes of genomic studies focused on identifying Ubx targets, it appears to be an appropriate time to revisit the Drosophila paradigm regarding insect wing development and evolution. Here, I review the recent findings related to Ubx during wing development, and discuss the impact of these findings on the current view of how Ubx came to regulate wing differentiation in the evolution of insect flight structures.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.