Proteotranscriptomic analyses of the midgut and Malpighian tubules after a sublethal concentration of Cry1Ab exposure on Spodoptera litura.

BACKGROUND: Cry1Ab has emerged as a bio-insecticide to control Spodoptera litura (Lepidoptera: Noctuidae). However, the sublethal effects of Cry1Ab on the physiological changes and molecular level of S. litura have not been well documented. Our aims in this study were to assess the sublethal effect of Cry1Ab on S. litura, including midgut and Malpighian tubules as targets. RESULTS: After sublethal Cry1Ab exposure, distinct histological alterations were mainly observed in the midgut. Furthermore, the results of comparative RNA sequencing and tandem mass tag-based proteomics showed that, in the midgut, most differential expression genes (DEGs) were up-regulated and significantly enriched in the serine protease activity pathway, and up-regulated differential expression proteins (DEPs) were mainly associated with the oxidative phosphorylation pathway, whereas the down-regulated involved in the ribosome pathways. In the Malpighian tubules, DEGs and DEPs were significantly enriched in the ribosome pathway. We proposed that ribosome may act as a universal target in energy metabolism with other pathways via the results of protein-protein interaction analysis. Further, by verification of the mRNA expression of some Cry protein receptor and detoxification genes after Cry1Ab treatment, it was suggested that the ribosomal proteins (RPs) possibly participate in influencing the Bt-resistance of S. litura larvae under sublethal Cry1Ab exposure. CONCLUSION: Under sublethal Cry1Ab exposure, the midgut of S. litura was damaged, and the proteotranscriptomic analysis elucidated that Cry1Ab disrupted the energy homeostasis of larvae. Furthermore, we emphasized the potential role of ribosomes in sublethal Cry1Ab exposure. © 2024 Society of Chemical Industry.

Lepidopteran wing scales contain abundant cross-linked film-forming histidine-rich cuticular proteins.

Scales are symbolic characteristic of Lepidoptera; however, nothing is known about the contribution of cuticular proteins (CPs) to the complex patterning of lepidopteran scales. This is because scales are resistant to solubilization, thus hindering molecular studies. Here we succeeded in dissolving developing wing scales from Bombyx mori, allowing analysis of their protein composition. We identified a distinctive class of histidine rich (His-rich) CPs (6%-45%) from developing lepidopteran scales by LC-MS/MS. Functional studies using RNAi revealed CPs with different histidine content play distinct and critical roles in constructing the microstructure of the scale surface. Moreover, we successfully synthesized films in vitro by crosslinking a 45% His-rich CP (BmorCPR152) with laccase2 using N-acetyl- dopamine or N-ß-alanyl-dopamine as the substrate. This molecular study of scales provides fundamental information about how such a fine microstructure is constructed and insights into the potential application of CPs as new biomaterials.

Circadian regulation of night feeding and daytime detoxification in a formidable Asian pest Spodoptera litura.

Voracious feeding, trans-continental migration and insecticide resistance make Spodoptera litura among the most difficult Asian agricultural pests to control. Larvae exhibit strong circadian behavior, feeding actively at night and hiding in soil during daytime. The daily pattern of larval metabolism was reversed, with higher transcription levels of genes for digestion (amylase, protease, lipase) and detoxification (CYP450s, GSTs, COEs) in daytime than at night. To investigate the control of these processes, we annotated nine essential clock genes and analyzed their transcription patterns, followed by functional analysis of their coupling using siRNA knockdown of interlocked negative feedback system core and repressor genes (SlituClk, SlituBmal1 and SlituCwo). Based on phase relationships and overexpression in cultured cells the controlling mechanism seems to involve direct coupling of the circadian processes to E-boxes in responding promoters. Additional manipulations involving exposure to the neonicotinoid imidacloprid suggested that insecticide application must be based on chronotoxicological considerations for optimal effectiveness.

A defective prostaglandin E synthase could affect egg formation in the silkworm Bombyx mori.

We had previously reported a prostaglandin E synthase (bmPGES) in the silkworm Bombyx mori that catalyzes the isomerization of PGH2 to PGE2. The present study aimed to provide a genome-editing characterization of bmPGES in B. mori. Results showed bmPGES gene disruption to result in a reduced content of PGE2. The change affected the expression of chorion genes and egg formation in silkworms. Collectively, the results indicated that bmPGES could be involved in reproduction of B. mori. Therefore, this study provides insights into the physiological role of bmPGES and PGE2 in silkworms.

Genome-wide annotation and comparative analysis of cuticular protein genes in the noctuid pest Spodoptera litura.

Insect cuticle is considered an adaptable and versatile building material with roles in the construction and function of exoskeleton. Its physical properties are varied, as the biological requirements differ among diverse structures and change during the life cycle of the insect. Although the bulk of cuticle consists basically of cuticular proteins (CPs) associated with chitin, the degree of cuticular sclerotization is an important factor in determining its physical properties. Spodoptera litura, the tobacco cutworm, is an important agricultural pest in Asia. Compared to the domestic silkworm, Bombyx mori, another lepidopteran whose CP genes have been well annotated, S. litura has a shorter life cycle, hides in soil during daytime beginning in the 5th instar and is exposed to soil in the pupal stage without the protection of a cocoon. In order to understand how the CP genes may have been adapted to support the characteristic life style of S. litura, we searched its genome and found 287 putative cuticular proteins that can be classified into 9 CP families (CPR with three groups (RR-1, RR-2, RR-3), CPAP1, CPAP3, CPF, CPFL, CPT, CPG, CPCFC and CPLCA), and a collection of unclassified CPs named CPH. There were also 112 cuticular proteins enriched in Histidine residues with content varying from 6% to 30%, comprising many more His-rich cuticular proteins than B. mori. A phylogenetic analysis between S. litura, M. sexta and B. mori uncovered large expansions of RR-1 and RR-2 CPs, forming large gene clusters in different regions of S. litura chromosome 9. We used RNA-seq analysis to document the expression profiles of CPs in different developmental stages and tissues of S. litura. The comparative genomic analysis of CPs between S. litura and B. mori integrated with the unique behavior and life cycle of the two species offers new insights into their contrasting ecological adaptations.

Genome-wide patterns of copy number variations in Spodoptera litura.

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.

A single amino acid substitution in the Bombyx-specific mucin-like membrane protein causes resistance to Bombyx mori densovirus.

Bombyx mori densovirus type 1 (BmDV) is a pathogen that causes flacherie disease in the silkworm. The absolute nonsusceptibility to BmDV among certain silkworm strains is determined independently by two genes, nsd-1 and Nid-1. However, neither of these genes has been molecularly identified to date. Here, we isolated the nsd-1 gene by positional cloning and characterized the properties of its product, NSD-1. Sequence and biochemical analyses revealed that this gene encodes a Bombyx-specific mucin-like glycoprotein with a single transmembrane domain. The NSD-1 protein was specifically expressed in the larval midgut epithelium, the known infection site of BmDV. Sequence analysis of the nsd-1 gene from 13 resistant and 12 susceptible strains suggested that a specific arginine residue in the extracellular tail of the NSD-1 protein was common among susceptible strains. Germline transformation of the susceptible-type nsd-1 (with a single nucleotide substitution) conferred partial susceptibility to resistant larvae, indicating that the + nsd-1 gene is required for the susceptibility of B. mori larvae to BmDV and the susceptibility is solely a result of the substitution of a single amino acid with arginine. Taken together, our results provide striking evidence that a novel membrane-bound mucin-like protein functions as a cell-surface receptor for a densovirus.

A new approach for comprehensively describing heterogametic sex chromosomes.

Notwithstanding the rapid developments in sequencing techniques, Y and W sex chromosomes have still been mostly excluded from whole genome sequencing projects due to their high repetitive DNA content. Therefore, Y and W chromosomes are poorly described in most species despite their biological importance. Several methods were developed for identifying Y or W-linked sequences among unmapped scaffolds. However, it is not enough to discover functional regions from short unmapped scaffolds. Here, we provide a new and simple strategy based on k-mer comparison for comprehensive analysis of the W chromosome in Bombyx mori. Using this novel method, we effectively assembled de novo 1281 W-derived genome contigs (totaling 1.9 Mbp), and identified 156 W-linked transcript RNAs and 345 W-linked small RNAs. This method will help in the elucidation of mechanisms of sexual development and exploration of W chromosome biological functions, and provide insights into the evolution of sex chromosomes. Moreover, we showed this method can be employed in identifying heterogametic sex chromosomes (W and Y chromosomes) in many other species where genomic information is still scarce.

Enzymatic characterization of two epsilon-class glutathione S-transferases of Spodoptera litura.

Two cDNAs encoding glutathione S-transferase (GST) of the tobacco cutworm, Spodoptera litura, were cloned by reverse transcriptase-polymerase chain reaction. The deduced amino acid sequences of the resulting clones revealed 32-51% identities to the epsilon-class GSTs from other organisms. The recombinant proteins were functionally overexpressed in Escherichia coli cells in soluble form and were purified to homogeneity. The enzymes were capable of catalyzing the bioconjugation of glutathione with 1-chloro-2,4-dinitrobenzene, 1,2-epoxy-3-(4-nitrophenoxy)-propane, and ethacrynic acid. A competition assay revealed that the GST activity was inhibited by insecticides, suggesting that it could be conducive to insecticide tolerance in the tobacco cutworm.

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest.

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect's natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India-South China-Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

Osiris9a is a major component of silk fiber in lepidopteran insects.

In a previous high-throughput proteomics study, it was found that the silkworm cocoon contains hundreds of complex proteins, many of which have unknown functions, in addition to fibroins, sericins, and some protease inhibitors. Osiris was one of the proteins with no known function. In this study, we identified the Osiris gene family members and constructed a phylogenetic tree based on the sequences from different species. Our results indicate that the Osiris9 gene subfamily contains six members; it is specifically expressed in lepidopteran insects and has evolved by gene duplication. An Osiris gene family member from Bombyx mori was designated as BmOsiris9a (BmOsi9a) on the basis of its homology to Drosophila melanogaster Osiris9. The expression pattern of BmOsi9a showed that it was highly expressed only in the middle silk gland of silkworm larvae, similar to Sericin1 (Ser1). BmOsi9a was visualized as two bands in western blot analysis, implying that it probably undergoes post-translational modifications. Immunohistochemistry analysis revealed that BmOsi9a was synthesized and secreted into the lumen of the middle silk gland, and was localized in the sericin layer in the silk fiber. BmOsi9a was found in the silk fibers of not only three Bombycidae species, viz. B. mori, B. mandarina, and B. huttoni, but also in the fibers collected from Saturniidae species, including Antheraea assama, Antheraea mylitta, and Samia cynthia. Although the exact biological function of Osi9a in the silk fibers is unknown, our results are important because they demonstrate that Osi9a is a common structural component of silk fiber and is expressed widely among the silk-producing Bombycidae and Saturniidae insects. Our results should help in understanding the role of Osi9a in silk fibers.

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori.

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO2 receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

The angiotensin-converting enzyme (ACE) gene family of Bombyx mori.

We previously reported regarding an ecdysone-inducible angiotensin-converting enzyme (ACE) gene. We found another four ACE genes in the Bombyx genome. The present study was undertaken to clarify the evolutionally changed function of the ACE of Bombyx mori. Core regions of deduced amino acid sequences of ACE genes were compared with those of other insect ACE genes. Five Bombyx genes have the conserved Zn2+-binding-site motif (HEXXH); however, BmAcer4 has only one and BmAcer3 has no catalytic ligand. BmAcer1 and BmAcer2 were expressed in several organs. BmAcer3 was expressed in testes, and BmAcer4 and BmAcer5 were expressed in compound eyes; however, the transcription levels of these three genes were very low. Quantitative RT-PCR and Western analysis were conducted to determine the tissue distribution and developmental expression of BmAcer1and BmAcer2. Transcripts of BmAcer1 and BmAcer2 were found in the reproductive organs during the larval and pupal stages. BmAcer1 was dominant in fat bodies during the feeding stage and showed high expression in the epidermis, wing discs, and pupal wing tissues after the wandering stage. Its expression patterns in epidermis, wing discs, and wing tissues resembled the hemolymph ecdysteroid titer in the larval and pupal stages. Acer1 was observed in the hemolymph at all stages, appearing to be the source of it are fat bodies, wings, and epidermis, and functioning after being secreted into the hemolymph. BmAcer2 was abundant in the midgut during the feeding stage and after the wandering stage and in silk glands after the pupal stage. We conclude that the evolution of BmAcer occurred through duplication, and, thereafter, functional diversification developed.

Role of Bmznf-2, a Bombyx mori CCCH zinc finger gene, in masculinisation and differential splicing of Bmtra-2.

Deciphering the regulatory factors involved in Bombyx mori sex determination has been a puzzle, challenging researchers for nearly a century now. The pre-mRNA of B. mori doublesex (Bmdsx), a master regulator gene of sexual differentiation, is differentially spliced, producing Bmdsxm and Bmdsxf transcripts in males and females respectively. The putative proteins encoded by these differential transcripts orchestrate antagonistic functions, which lead to sexual differentiation. A recent study in B. mori illustrated the role of a W-derived fem piRNA in conferring femaleness. In females, the fem piRNA was shown to suppress the activity of a Z-linked CCCH type zinc finger (znf) gene, Masculiniser (masc), which indirectly promotes the Bmdsxm type of splicing. In this study, we report a novel autosomal (Chr 25) CCCH type znf motif encoding gene Bmznf-2 as one of the potential factors in the Bmdsx sex specific differential splicing, and we also provide insights into its role in the alternative splicing of Bmtra2 by using ovary derived BmN cells. Over-expression of Bmznf-2 induced Bmdsxm type of splicing (masculinisation) with a correspondingly reduced expression of Bmdsxf type isoform in BmN cells. Further, the site-directed mutational studies targeting the tandem CCCH znf motifs revealed their indispensability in the observed phenotype of masculinisation. Additionally, the dual luciferase assays in BmN cells using 5' UTR region of the Bmznf-2 strongly implied the existence of a translational repression over this gene. From these findings, we propose Bmznf-2 to be one of the potential factors of masculinisation similar to Masc. From the growing number of Bmdsx splicing regulators, we assume that the sex determination cascade of B. mori is quite intricate in nature; hence, it has to be further investigated for its comprehensive understanding.

Gene expression analysis in the larval silk gland of the eri silkworm Samia ricini.

Insects produce silk for a range of purposes. In the Lepidoptera, silk is utilized as a material for cocoon production and serves to protect larvae from adverse environmental conditions or predators. Species in the Saturniidae family produce an especially wide variety of cocoons, for example, large, golden colored cocoons and those with many small holes. Although gene expression in the silk gland of the domestic silkworm (Bombyx mori L.) has been extensively studied, considerably fewer investigations have focused on members of the saturniid family. Here, we established expression sequence tags from the silk gland of the eri silkworm (Samia ricini), a saturniid species, and used these to analyze gene expression. Although we identified the fibroin heavy chain gene in the established library, genes for other major silk proteins, such as fibroin light chain and fibrohexamerin, were absent. This finding is consistent with previous reports that these latter proteins are lacking in saturniid silk. Recently, a series of fibrohexamerin-like genes were identified in the Bombyx genome. We used this information to conduct a detailed analysis of the library established here. This analysis identified putative homologues of these genes. We also found several genes encoding small silk protein molecules that are also present in the silk of other Lepidoptera. Gene expression patterns were compared between eri and domestic silkworm, and both conserved and nonconserved expression patterns were identified for the tested genes. Such differential gene expression might be one of the major causes of the differences in silk properties between these species. We believe that our study can be of value as a basic catalogue for silk gland gene expression, which will yield to the further understanding of silk evolution.

Construction, complete sequence, and annotation of a BAC contig covering the silkworm chorion locus.

The silkmoth chorion was studied extensively by F.C. Kafatos' group for almost 40 years. However, the complete structure of the chorion locus was not obtained in the genome sequence of Bombyx mori published in 2008 due to repetitive sequences, resulting in gaps and an incomplete view of the locus. To obtain the complete sequence of the chorion locus, expressed sequence tags (ESTs) derived from follicular epithelium cells were used as probes to screen a bacterial artificial chromosome (BAC) library. Seven BACs were selected to construct a contig which covered the whole chorion locus. By Sanger sequencing, we successfully obtained complete sequences of the chorion locus spanning 871,711 base pairs on chromosome 2, where we annotated 127 chorion genes. The dataset reported here will recruit more researchers to revisit one of the oldest model systems which has been used to study developmentally regulated gene expression. It also provides insights into egg development and fertilization mechanisms and is relevant to applications related to improvements in breeding procedures and transgenesis.

A comprehensive analysis of the chorion locus in silkmoth.

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as "middle", and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Molecular architecture of silk fibroin of Indian golden silkmoth, Antheraea assama.

The golden silk spun by Indian golden silkmoth Antheraea assama, is regarded for its shimmering golden luster, tenacity and value as biomaterial. This report describes the gene coding for golden silk H-fibroin (AaFhc), its expression, full-length sequence and structurally important motifs discerning the underlying genetic and biochemical factors responsible for its much sought-after properties. The coding region, with biased isocodons, encodes highly repetitious crystalline core, flanked by a pair of 5' and 3' non-repetitious ends. AaFhc mRNA expression is strictly territorial, confined to the posterior silk gland, encoding a protein of size 230 kDa, which makes homodimers making the elementary structural units of the fibrous core of the golden silk. Characteristic polyalanine repeats that make tight ß-sheet crystals alternate with non-polyalanine repeats that make less orderly antiparallel ß-sheets, ß-turns and partial α-helices. Phylogenetic analysis of the conserved N-terminal amorphous motif and the comparative analysis of the crystalline region with other saturniid H-fibroins reveal that AaFhc has longer, numerous and relatively uniform repeat motifs with lower serine content that assume tighter ß-crystals and denser packing, which are speculated to be responsible for its acclaimed properties of higher tensile strength and higher refractive index responsible for golden luster.

Functional Loss of Bmsei Causes Thermosensitive Epilepsy in Contractile Mutant Silkworm, Bombyx mori.

The thermoprotective mechanisms of insects remain largely unknown. We reported the Bombyx mori contractile (cot) behavioral mutant with thermo-sensitive seizures phenotype. At elevated temperatures, the cot mutant exhibit seizures associated with strong contractions, rolling, vomiting, and a temporary lack of movement. We narrowed a region containing cot to ~268 kb by positional cloning and identified the mutant gene as Bmsei which encoded a potassium channel protein. Bmsei was present in both the cell membrane and cytoplasm in wild-type ganglia but faint in cot. Furthermore, Bmsei was markedly decreased upon high temperature treatment in cot mutant. With the RNAi method and injecting potassium channel blockers, the wild type silkworm was induced the cot phenotype. These results demonstrated that Bmsei was responsible for the cot mutant phenotype and played an important role in thermoprotection in silkworm. Meanwhile, comparative proteomic approach was used to investigate the proteomic differences. The results showed that the protein of Hsp-1 and Tn1 were significantly decreased and increased on protein level in cot mutant after thermo-stimulus, respectively. Our data provide insights into the mechanism of thermoprotection in insect. As cot phenotype closely resembles human epilepsy, cot might be a potential model for the mechanism of epilepsy in future.

Identification of two juvenile hormone inducible transcription factors from the silkworm, Bombyx mori.

Juvenile hormone (JH) regulates many physiological processes in insects. However, the signal cascades in which JH is active have not yet been fully elucidated, particularly in comparison to another major hormone ecdysteroid. Here we identified two JH inducible transcription factors as candidate components of JH signaling pathways in the silkworm, Bombyx mori. DNA microarray analysis showed that expression of two transcription factor genes, E75 and Enhancer of split mß (E(spl)mß), was induced by juvenile hormone I (JH I) in NIAS-Bm-aff3 cells. Real time RT-PCR analysis confirmed that expression of four E75 isoforms (E75A, E75B, E75C and E75D) and E(spl)mß was 3-8 times greater after JH I addition. Addition of the protein synthesis inhibitor cycloheximide did not suppress JH-induced expression of the genes, indicating that they were directly induced by JH. JH-induced expression of E75 and E(spl)mß was also observed in four other B. mori cell lines and in larval hemocytes of final instar larvae. Notably, E75A expression was induced very strongly in larval hemocytes by topical application of the JH analog fenoxycarb; the level of induced expression was comparable to that produced by feeding larvae with 20-hydroxyecdysone. These results suggest that E75 and E(spl)mß are general and direct target genes of JH and that the transcription factors encoded by these genes play important roles in JH signaling.