NADP+-dependent farnesol dehydrogenase, a corpora allata enzyme involved in juvenile hormone synthesis

Engineering of insect juvenile hormone III biosynthesis in the plant Nicotiana benthamiana

Juvenile hormones (JHs) are farnesoic acid-derived sesquiterpenoids that play a crucial role in regulating various developmental processes in insects. Based on these reported biological activities, JHs and their synthetic analogs have been utilized as insecticides with significant commercial success over the past years. Here we describe the engineering of the JH pathway of the yellow fever mosquito (Aedes aegypti) by transient gene expression in the plant Nicotiana benthamiana. This approach led to the successful production of JH III in N. benthamiana leaves at a concentration of ca. 10 μg/g fresh weight. The co-expression of a feedback-insensitive version of 3-hydroxy-3-methylglutaryl coenzyme A reductase from Arabidopsis thaliana further increased the titer eight-fold from 10 to 80 μg/g fresh weight. Our efforts also revealed that the rich endogenous metabolic background of N. benthamiana can generate farnesoic acid, a key precursor to JH III, and thus, only 3 genes need to be expressed to provide high titers of this compound. Our study demonstrates the production of high titers of JH III in N. benthamina via heterologous expression of insect JH biosynthetic genes.

Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

Evolution, structure, and drug-metabolizing activity of mammalian prenylcysteine oxidases

Prenylcysteine oxidases (PCYOXs) metabolize prenylated cysteines produced by protein degradation. They utilize oxygen as a co-substrate to produce free cysteine, an aldehyde, and hydrogen peroxide through the unusual oxidation of a thioether bond. In this study, we explore the evolution, structure, and mechanism of the two mammalian PCYOXs. A gene duplication event in jawed vertebrates originated in these two paralogs. Both enzymes are active on farnesyl- and geranylgeranylcysteine, but inactive on molecules with shorter prenyl groups. Kinetics experiments outline a mechanism where flavin reduction and re-oxidation occur rapidly without any detectable intermediates, with the overall reaction rate limited by product release. The experimentally determined three-dimensional structure of PCYOX1 reveals long and wide tunnels leading from the surface to the flavin. They allow the isoprene substrate to curl up within the protein and position its reactive cysteine group close to the flavin. A hydrophobic patch on the surface mediates membrane association, enabling direct substrate and product exchange with the lipid bilayer. Leveraging established knowledge of flavoenzyme inhibition, we designed sub-micromolar PCYOX inhibitors. Additionally, we discovered that PCYOXs bind and slowly degrade salisirab, an anti-RAS compound. This activity suggests potential and previously unknown roles of PCYOXs in drug metabolism.

Cytoskeleton-associated gelsolin responds to the midgut distention process in saline meal-fed Aedes aegypti and affects arbovirus dissemination from the midgut

The mosquito, Aedes aegypti, is the principal vector for several arboviruses. The mosquito midgut is the initial tissue that gets infected with an arbovirus acquired along with a blood meal from a vertebrate host. Blood meal ingestion leads to midgut tissue distention thereby increasing the pore size of the surrounding basal lamina. This allows newly synthesized virions to exit the midgut by traversing the distended basal lamina to infect secondary tissues of the mosquito. We conducted a quantitative label-free proteomic time course analysis with saline meal-fed Ae. aegypti females to identify host factors involved in midgut tissue distention. Around 2000 proteins were detected during each of the seven sampling time points and 164 of those were uniquely expressed. Forty-five of 97 differentially expressed proteins were upregulated during the 96-h time course and most of those were involved in cytoskeleton modulation, metabolic activity, and vesicle/vacuole formation. The F-actin-modulating Ae. aegypti (Aa)-gelsolin was selected for further functional studies. Stable knockout of Aa-gelsolin resulted in a mosquito line, which showed distorted actin filaments in midgut-associated tissues likely due to diminished F-actin processing by gelsolin. Zika virus dissemination from the midgut of these mosquitoes was diminished and delayed. The loss of Aa-gelsolin function was associated with an increased induction of apoptosis in midgut tissue indicating an involvement of Aa-gelsolin in apoptotic signaling in mosquitoes. Here, we used proteomics to discover a novel host factor, Aa-gelsolin, which affects the midgut escape barrier for arboviruses in mosquitoes and apoptotic signaling in the midgut.

Exploring the Wilderness within: An Integrative Metabolomics and Transcriptomics Study on Near-Wild and Colonized Aedes aegypti

This study examines the phenotypic differences between wild-derived F2 Central Valley mosquitoes and the insecticide-susceptible Rockefeller (Rock) lab strain of Ae. aegypti. Given the rarity of wild pyrethroid-susceptible populations, the focus of this work is to develop an understanding of the resistance physiology in this invasive mosquito population and explore the potential of metabolites as diagnostic biomarkers for metabolic resistance. This study utilizes metabolomic, gene expression, and lifespan data for a comparison between strains. The findings indicate that wild-derived mosquitoes with greater metabolic resistance have a lifespan sensitivity to restricted larval nutrition. In terms of metabolism and gene expression, Central Valley mosquitoes show increased activity in oxidoreductase, glutathione metabolism, and the pentose phosphate pathway. Conversely, Rock mosquitoes display signs of metabolic inefficiency and mitochondrial dysregulation, likely tolerated due to the consistency and nutritional abundance of a controlled lab environment. The study also examines Ae. aegypti P450 and GSTE profiles in relation to other insecticide-resistant groups. While metabolomic data can differentiate our study groups, the challenges in biomarker development arise from few detected markers meeting high fold change thresholds.

A nano-delivery system expands the insecticidal target of thiamethoxam to include a devastating pest, the fall armyworm

Nano-delivery systems have been applied to deliver various synthetic/botanical pesticides to increase the efficiency of pesticide use and reduce the volumes of pesticides applied. Previous studies have supported the hypothesis that the nanocarriers can help expand the insecticidal target of pesticides to include non-target pests. However, the potential mechanism underlying this interesting phenomenon remains unclear. Herein, a widely applied star polycation (SPc) nanocarrier was synthesized to construct a thiamethoxam (TMX) nano-delivery system. The SPc-based delivery system could promote the translocation of exogenous substances across the membrane of Sf9 cells, increase the cytotoxicity of TMX against Sf9 cells by nearly 20%, and expand the insecticidal target of TMX to include Spodoptera frugiperda (the fall armyworm), with a 27.5% mortality increase at a concentration of 0.25 mg/mL. Moreover, the RNA-seq analysis demonstrated that the SPc could upregulate various transport-related genes, such as Rab, SORT1, CYTH, and PIKfyve, for the enhanced cellular uptake of TMX. Furthermore, enhanced cell death in larvae treated with the TMX-SPc complex was observed through changes in the expression levels of death-related genes, such as Casp7, BIRC5, MSK1, and PGAM5. The SPc-based nano-delivery system improved the cellular uptake of TMX and expanded its insecticidal target by adjusting the expression levels of death-related genes. The current study mainly identified the transport and cell death genes related to nanocarrier-based insecticidal target expansion, which is beneficial for understanding the bioactivity enhancement of the nano-delivery system.

Farnesol dehydrogenase from Helicoverpa armigera (Hübner) as a promising target for pest management: molecular docking, in vitro and insect bioassay studies using geranylgeraniol as potential inhibitor

Juvenile hormone (JH) plays pivotal roles in several critical developmental processes in insects, including metamorphosis and reproduction. JH-biosynthetic pathway enzymes are considered highly promising targets for discovering novel insecticides. The oxidation of farnesol to farnesal, catalysed by farnesol dehydrogenase (FDL), represents a rate-limiting step in JH biosynthesis. Here, we report farnesol dehydrogenase (HaFDL) from H. armigera as a promising insecticidal target. The inhibitory potential of natural substrate analogue geranylgeraniol (GGol) was tested in vitro, wherein it showed a high binding affinity (kd 595 µM) for HaFDL in isothermal titration calorimetry (ITC) and subsequently exhibited dose-dependent enzyme inhibition in GC-MS coupled qualitative enzyme inhibition assay. Moreover, the experimentally determined inhibitory activity of GGol was augmented by the in silico molecular docking simulation which showed that GGol formed a stable complex with HaFDL, occupied the active site pocket and interacted with key active site residues (Ser147 and Tyr162) as well as other residues that are crucial in determining the active site architecture. Further, the diet-incorporated oral feeding of GGol caused detrimental effects on larval growth and development, exhibiting a significantly reduced rate of larval weight gain (P < 0.01), aberrant pupal and adult morphogenesis, and a cumulative mortality of ~ 63%. To the best of our knowledge, the study presents the first report on evaluating GGol as a potential inhibitor for HaFDL. Overall, the findings revealed the suitability of HaFDL as a potential insecticidal target for the management H. armigera.

Isolation, purification and biochemical characterization of Conopomorpha cramerella farnesol dehydrogenase

In Southeast Asia, Conopomorpha cramerella (Snellen) which is commonly known as the cocoa pod borer (CPB) moth has been identified as the most detrimental pest of Theobroma cacao L. Apart from the various side effects on human health and non-target organisms, heavily relying on synthetic pyrethroid insecticides to control CPB infestations also increases the environmental contamination risks. Thus, developing biorational insecticides that minimally affect the non-target organism and environment by targeting the insect growth regulation process is needed to manage the pest population. In insects, juvenile hormones (JH) regulate critical biological events, especially metamorphosis, development and reproduction. Since the physiological roles of JH III vary among different organisms, the biochemical properties, especially substrate specificity and analogue inhibition, may also be different. Therefore, studies on the JH III biosynthetic pathway enzymes in both plants and insects are beneficial to discover more effective analogues. Bioinformatic analysis and biochemical characterization of a NADP⁺ -dependent farnesol dehydrogenase, an intermediate enzyme of the JH III pathway, from C. cramerella (CcFolDH), were described in this study. In addition, the farnesol analogues that may act as a potent analogue inhibitor for CcFolDH ware determined using in vitro enzymatic study. The phylogenetic analysis indicated that CcFolDH shared a close phylogenetic relationship to the honeybee's short-chain dehydrogenase/reductase. The 27 kDa CcFolDH has an NADP(H) binding domain with a typical Rossmann fold and is likely a homotetrameric protein in the solution. The enzyme had a greater preference for substrate trans, trans-farnesol and coenzyme NADP⁺ . In terms of analogue inhibitor inhibition, hexahydroxyfarnesyl acetone showed the highest inhibition (the lowest K_i ) compared to other farnesol analogues. Thus, hexahydroxyfarnesyl acetone would serve as the most potent active ingredient for future biorational pesticide management for C. cramerella infestation. Based on the bioinformatic analyses and biochemical characterizations conducted in this research, we proposed that rCcFolDH differs slightly from other reported farnesol dehydrogenases in terms of molecular weight, substrate preference, coenzymes utilization and analogue inhibitors selection.

Ensemble-based, high-throughput virtual screening of potential inhibitor targeting putative farnesol dehydrogenase of Metisa plana (Lepidoptera: Psychidae)

Metisa plana (Lepidoptera: Psychidae) bagworm is a leaf-eater caterpillar ubiquitously found as a damaging pest in oil palm plantations, specifically in Malaysia. Various strategies have been implemented, including the usage of chemical insecticides. However, the main challenges include the development of insecticide resistance and its detrimental effects on the environment and non-target organisms. Therefore, a biorational insecticide is introduced by targeting the juvenile hormone (JH) biosynthetic pathway, which is mainly present in the insect and vital for the insect's growth, diapause, metamorphosis, and adult reproduction. This study aimed to investigate the potential inhibitor for the rate-limiting enzyme involved in the JH pathway known as farnesol dehydrogenase. A 255 amino acids sequence encoded for the putative M. plana farnesol dehydrogenase (MpFolDH) open reading frame had been identified and isolated. The three-dimensional structure of MpFolDH was predicted to have seven β- sheets with α-helices at both sides, showing typical characteristics for classical short-chain dehydrogenase and associated with oxidoreductase activity. Then, the ensemble-based virtual screening was conducted based on the ZINC20 database, in which 43 768 compounds that fulfilled pesticide-likeness criteria were screened by site-specific molecular docking. After a short molecular dynamics simulation (5 ns) was conducted towards 102 compounds, only the top 10 compounds based on their most favourable binding energy were selected for a more extended simulation (100 ns). Based on the protein-ligand stability, protein compactness, residues rigidity, binding interaction, binding energy throughout the 100 ns simulation, and physicochemical analysis, ZINC000408743205 was selected as a potential inhibitor for this enzyme. Amino acids decomposition analysis indicates Ile18, Ala95, Val198 and Val202 were the critical contributor residues for MpFolDH-inhibitors(s) complex.

Transcriptome changes of Apis mellifera female embryos with fem gene knockout by CRISPR/Cas9

The sex of honey bees is decided by a regulatory cascade comprising of csd, fem and Amdsx. In order to further identify other genes involved in sex determination and differentiation of honey bees in the early stages of embryo development, the CRISPR/Cas9 method was used to knock out fem gene in the embryonic stage of diploid western honey bees, and RNA-seq was used to analyze gene expression changes in the embryo after fem knockout. Finally, we found that the bees had undergone gender changes due to fem knockout. A total of 155 differentially expressed genes (DEGs) were obtained, with 48 up-regulated and 107 down-regulated DEGs in the mutant group compared to the control group. Of them, many genes are related to sex development or differentiation. In addition, 1502 differentially expressed alternative splicing events (DEASEs) related to 1011 genes, including the main honey bee sex-determining genes csd, tra2, fem, and Amdsx, were identified between the mutant group and control group, indicating that fem regulates alternative splicing of a large number of downstream genes. Our results provide valuable clues for further investigating the molecular mechanism of sex determination and differentiation in honey bees.

Differential gene expression underlying the biosynthesis of Dufour's gland signals in Bombus impatiens

Pheromones regulating social behavior are one of the most explored phenomena in social insects. However, compound identity, biosynthesis and their genetic basis are known in only a handful of species. Here we examined the gene expression associated with pheromone biosynthesis of two main chemical classes: esters and terpenes, using the social bee Bombus impatiens. We conducted chemical and RNA-seq analyses of the Dufour's gland, an exocrine gland producing a plethora of pheromones regulating social behavior in hymenopteran species. The Dufour's gland contains mostly long-chained hydrocarbons, terpenes and esters that signal reproductive and social status in several bee species. In bumble bees, the Dufour's gland contains queen- and worker-specific esters, in addition to terpenes and terpene-esters only found in gynes and queens. These compounds are assumed to be synthesized de novo in the gland, however, their genetic basis is unknown. A whole transcriptome gene expression analysis of the gland in queens, gynes, queenless and queenright workers showed distinct transcriptomic profiles, with thousands of differentially expressed genes between the groups. Workers and queens express genes associated with key enzymes in the biosynthesis of wax esters, while queens and gynes preferentially express key genes in terpene biosynthesis. Overall, our data demonstrate gland-specific regulation of chemical signals associated with social behavior and identifies candidate genes and pathways regulating caste-specific chemical signals in social insects.

Biosynthesis of geranate via isopentenol utilization pathway in Escherichia coli

Isoprenoids are a large family of natural products with diverse structures, which allow them to play diverse and important roles in the physiology of plants and animals. They also have important commercial uses as pharmaceuticals, flavoring agents, fragrances, and nutritional supplements. Recently, metabolic engineering has been intensively investigated and emerged as the technology of choice for the production of isoprenoids through microbial fermentation. Isoprenoid biosynthesis typically originates in plants from acetyl-coA in central carbon metabolism, however, a recent study reported an alternative pathway, the isopentenol utilization pathway (IUP), that can provide the building blocks of isoprenoid biosynthesis from affordable C5 substrates. In this study, we expressed the IUP in Escherichia coli to efficiently convert isopentenols into geranate, a valuable isoprenoid compound. We first established a geraniol-producing strain in E. coli that uses the IUP. Then, we extended the geraniol synthesis pathway to produce geranate through two oxidation reactions catalyzed by two alcohol/aldehyde dehydrogenases from Castellaniella defragrans. The geranate titer was further increased by optimizing the expression of the two dehydrogenases and also parameters of the fermentation process. The best strain produced 764 mg/L geranate in 24 h from 2 g/L isopentenols (a mixture of isoprenol and prenol). We also investigated if the dehydrogenases could accept other isoprenoid alcohols as substrates.

Structural characterization and functional analysis of juvenile hormone acid methyltransferase JHAMT3 from the silkworm, Bombyx mori

Juvenile hormone acid methyltransferase (JHAMT) is a rate-limiting enzyme of juvenile hormone (JH) biosynthesis in insects. It transfers the methyl group of S-adenosyl methionine to either the carboxyl group of JH acids or farnesoic acid to produce JH. Six JHAMT paralogues have been identified in the silkworm (Bombyx mori); among them, JHAMT1 and JHAMT2 display a methyltransferase activity. Here, the three-dimensional crystal structure of inactive JHAMT3 and the binary complex of JHAMT3 with its cofactor S-adenosyl-l-homocysteine were determined through X-ray crystallization. Comparative structural analysis revealed that JHAMT3 adopted a similar structural pattern to that of functional JHAMT2, which comprised one core Rossmann fold domain and one substrate-binding domain. Similar to JHAMT2, JHAMT3 underwent a conformational change at the Rossmann fold domain because of cofactor binding, which promoted ligand accommodation. However, it exhibited a relatively rigid substrate-binding pocket compared with that of JHAMT2. JHAMT3 was also highly expressed in the silk gland of fourth- and fifth-instar B. mori larvae. The results of expression profiling combined with activity analysis suggested that JHAMT3 might function as a binding protein of JH acids for the regulation of JH acid titers. These findings provide a structural basis for enhancing the understanding of the physiological function of JHAMT3 and a rational framework for the development of potent and specific inhibitors of JHAMT family members.

Selection and validation of reference genes for RT-qPCR-based analyses of Anastatus japonicus Ashmead (Hymenoptera: Helicopteridae)

RT-qPCR remains a vital approach for molecular biology studies aimed at quantifying gene expression in a range of physiological or pathological settings. However, the use of appropriate reference genes is essential to attain meaningful RT-qPCR results. Anastatus japonicus Ashmead (Hymenoptera: Helicopteridae) is an important egg parasitoid wasp and natural enemy of fruit bugs and forest caterpillars. While recent transcriptomic studies have analyzed gene expression profiles in A. japonicus specimens, offering a robust foundation for functional research focused on this parasitoid, no validated A. japonicus reference genes have yet been established, hampering further research efforts. Accordingly, this study sought to address this issue by screening for the most stable internal reference genes in A. japonicus samples to permit reliable RT-qPCR analyses. The utility of eight candidate reference genes (ACTIN, TATA, GAPDH, TUB, RPL13, RPS6, EF1α, RPS3a) was assessed under four different conditions by comparing developmental stages (larvae, pupae, adults), tissues (abdomen, chest, head), sex (male or female adults), or diapause states (diapause induction for 25, 35, 45, or 55 days, or diapause termination). RefFinder was used to calculate gene stability based on the integration of four algorithms (BestKeeper, Normfinder, geNorm, and ΔCt method) to determine the optimal RT-qPCR reference gene. Based on this approach, RPS6 and RPL13 were found to be the most reliable reference genes when assessing different stages of development, while ACTIN and EF1α were optimal when comparing adults of different sexes, RPL13 and EF1α were optimal when analyzing different tissues, and TATA and ACTIN were optimal for different diapause states. These results provide a valuable foundation for future RT-qPCR analyses of A. japonicus gene expression and function under a range of experimental conditions.

Tachykinin-related peptides modulate immune-gene expression in the mealworm beetle Tenebrio molitor L

Tachykinins (TKs) are a group of conserved neuropeptides. In insects, tachykinin-related peptides (TRPs) are important modulators of several functions such as nociception and lipid metabolism. Recently, it has become clear that TRPs also play a role in regulating the insect immune system. Here, we report a transcriptomic analysis of changes in the expression levels of immune-related genes in the storage pest Tenebrio molitor after treatment with Tenmo-TRP-7. We tested two concentrations (10^-8 and 10^-6 M) at two time points, 6 and 24 h post-injection. We found significant changes in the transcript levels of a wide spectrum of immune-related genes. Some changes were observed 6 h after the injection of Tenmo-TRP-7, especially in relation to its putative anti-apoptotic action. Interestingly, 24 h after the injection of 10^-8 M Tenmo-TRP-7, most changes were related to the regulation of the cellular response. Applying 10^-6 M Tenmo-TRP-7 resulted in the downregulation of genes associated with humoral responses. Injecting Tenmo-TRP-7 did not affect beetle survival but led to a reduction in haemolymph lysozyme-like antibacterial activity, consistent with the transcriptomic data. The results confirmed the immunomodulatory role of TRP and shed new light on the functional homology between TRPs and TKs.

Crystal structure and molecular characterization of NADP+-farnesol dehydrogenase from cotton bollworm, Helicoverpaarmigera

Farnesol dehydrogenase (FDL) orchestrates the oxidation reaction catalyzing farnesol to farnesal, a key step in the juvenile hormone (JH) biosynthesis pathway of insects and hence, represents a lucrative target for developing insect growth regulators (IGRs). However, information on the structural and functional characterization of JH-specific farnesol dehydrogenase in insects remains elusive. Herein, we identified a transcript that encodes farnesol dehydrogenase (HaFDL) from Helicoverpa armigera, a major pest of cotton. The investigations of molecular assembly, biochemical analysis and spatio-temporal expression profiling showed that HaFDL exists as a soluble homo-tetrameric form, exhibits a broad substrate affinity and is involved in the JH-specific farnesol oxidation in H. armigera. Additionally, the study presents the first crystal structure of the HaFDL-NADP enzyme complex determined at 1.6 Å resolution. Structural analysis revealed that HaFDL belongs to the NADP-specific cP2 subfamily of the classical short-chain dehydrogenase/reductase (SDR) family and exhibits typical structural features of those enzymes including the conserved nucleotide-binding Rossman-fold. The isothermal titration calorimetry (ITC) showed a high binding affinity (dissociation constant, Kd, 3.43 μM) of NADP to the enzyme. Comparative structural analysis showed a distinct substrate-binding pocket (SBP) loop with a spacious and hydrophobic substrate-binding pocket in HaFDL, consistent with the biochemically observed promiscuous substrate specificity. Finally, based on the crystal structure, substrate modeling and structural comparison with homologs, a two-step reaction mechanism is proposed. Overall, the findings significantly impact and contribute to our understanding of farnesol dehydrogenase functional properties in JH biosynthesis in H. armigera.

Interspecies Isobaric Labeling-Based Quantitative Proteomics Reveals Protein Changes in the Ovary of Aedes aegypti Coinfected With ZIKV and Wolbachia

Zika is a vector-borne disease caused by an arbovirus (ZIKV) and overwhelmingly transmitted by Ae. aegypti. This disease is linked to adverse fetal outcomes, mostly microcephaly in newborns, and other clinical aspects such as acute febrile illness and neurologic complications, for example, Guillain-Barré syndrome. One of the most promising strategies to mitigate arbovirus transmission involves releasing Ae. aegypti mosquitoes carrying the maternally inherited endosymbiont bacteria Wolbachia pipientis. The presence of Wolbachia is associated with a reduced susceptibility to arboviruses and a fitness cost in mosquito life-history traits such as fecundity and fertility. However, the mechanisms by which Wolbachia influences metabolic pathways leading to differences in egg production remains poorly known. To investigate the impact of coinfections on the reproductive tract of the mosquito, we applied an isobaric labeling-based quantitative proteomic strategy to investigate the influence of Wolbachia wMel and ZIKV infection in Ae. aegypti ovaries. To the best of our knowledge, this is the most complete proteome of Ae. aegypti ovaries reported so far, with a total of 3913 proteins identified, were also able to quantify 1044 Wolbachia proteins in complex sample tissue of Ae. aegypti ovary. Furthermore, from a total of 480 mosquito proteins modulated in our study, we discuss proteins and pathways altered in Ae. aegypti during ZIKV infections, Wolbachia infections, coinfection Wolbachia/ZIKV, and compared with no infection, focusing on immune and reproductive aspects of Ae. aegypti. The modified aspects mainly were related to the immune priming enhancement by Wolbachia presence and the modulation of the Juvenile Hormone pathway caused by both microorganism’s infection.

Rethinking Sesquiterpenoids: A Widespread Hormone in Animals

The sesquiterpenoid hormone juvenile hormone (JH) controls development, reproduction, and metamorphosis in insects, and has long been thought to be confined to the Insecta. While it remains true that juvenile hormone is specifically synthesized in insects, other types or forms of sesquiterpenoids have also been discovered in distantly related animals, such as the jellyfish. Here, we combine the latest literature and annotate the sesquiterpenoid biosynthetic pathway genes in different animal genomes. We hypothesize that the sesquiterpenoid hormonal system is an ancestral system established in an animal ancestor and remains widespread in many animals. Different animal lineages have adapted different enzymatic routes from a common pathway, with cnidarians producing farnesoic acid (FA); non-insect protostomes and non-vertebrate deuterostomes such as cephalochordate and echinoderm synthesizing FA and methyl farnesoate (MF); and insects producing FA, MF, and JH. Our hypothesis revolutionizes the current view on the sesquiterpenoids in the metazoans, and forms a foundation for a re-investigation of the roles of this important and yet neglected type of hormone in different animals.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Genome and transcriptome analysis of the mealybug Maconellicoccus hirsutus: Correlation with its unique phenotypes

Mealybugs are aggressive pests with world-wide distribution and are suitable for the study of different phenomena like genomic imprinting and epigenetics. Genomic approaches facilitate these studies in absence of robust genetics in this system. We sequenced, de novo assembled, annotated Maconellicoccus hirsutus genome. We carried out comparative genomics it with four mealybug and eight other insect species, to identify expanded, specific and contracted gene classes that relate to pesticide and desiccation resistance. We identified horizontally transferred genes adding to the mutualism between the mealybug and its endosymbionts. Male and female transcriptome analysis indicates differential expression of metabolic pathway genes correlating with their physiology and the genes for sexual dimorphism. The significantly lower expression of endosymbiont genes in males relates to the depletion of endosymbionts in males during development.

Molecular characterization and enzyme inhibition studies of NADP+- farnesol dehydrogenase from diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

Juvenile hormone III (JH III) plays an important role in insect reproduction, development, and behavior. The second branch of JH III production includes oxidation of farnesol to farnesal by farnesol dehydrogenase. This study reported the identification and characterization of Plutella xylostella farnesol dehydrogenase (PxFoLDH). Our results showed that PxFoLDH belongs to the short-chain dehydrogenase/reductase superfamily, consisting of a single domain with a structurally conserved Rossman fold, an NAD(P) (H)-binding region and a structurally diverse C-terminal region. The purified enzyme displayed maximum activity at 55$\ $°C with pH 9.5 and was stable in the temperature below 70$\ ^\circ $C. PxFoLDH was determined to be a monomer with a relative molecular weight of 27 kDa and highly specific for trans, trans-farnesol, and NADP+. Among analog inhibitors tested, farnesyl acetate was the most effective inhibitor with the lowest Ki value of 0.02 µm. Our findings showed this purified enzyme may represent as NADP+-farnesol dehydrogenase.

Purification, biochemical characterisation and bioinformatic analysis of recombinant farnesol dehydrogenase from Theobroma cacao

The juvenile hormones (JH) in plants are suggested to act as a form of plant defensive strategy especially against insect herbivory. The oxidation of farnesol to farnesoic acid is a key step in the juvenile hormone biosynthesis pathway. We herein present the purification and characterisation of the recombinant Theobroma cacao farnesol dehydrogenase enzyme that catalyses oxidation of farnesol to farnesal. The recombinant enzyme was purified to apparent homogeneity by affinity chromatography. The purified enzyme was characterised in terms of its deduced amino acid sequences, phylogeny, substrate specificity, kinetic parameters, structural modeling, and docking simulation. The phylogenetic analysis indicated that the T. cacao farnesol dehydrogenase (TcFolDH) showed a close relationship with A. thaliana farnesol dehydrogenase gene. The TcFolDH monomer had a large N-terminal domain which adopted a typical Rossmann-fold, harboring the GxxGxG motif (NADP(H)-binding domain) and a small C-terminal domain. The enzyme was a homotrimer comprised of subunits with molecular masses of 36 kDa. The TcFolDH was highly specific to NADP⁺ as coenzyme. The substrate specificity studies showed trans, trans-farnesol was the most preferred substrate for the TcFolDH, suggesting that the purified enzyme was a NADP⁺-dependent farnesol dehydrogenase. The docking of trans, trans-farnesol and NADP⁺ into the active site of the enzyme showed the important residues, and their interactions involved in the substrate and coenzyme binding of TcFolDH. Considering the extensive involvement of JH in both insects and plants, an in-depth knowledge on the recombinant production of intermediate enzymes of the JH biosynthesis pathway could help provide a potential method for insect control.

Identification of azadirachtin responsive genes in Spodoptera frugiperda larvae based on RNA-seq

The fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) is a polyphagous pest with 353 plant species as its hosts, including maize, sorghum, cotton, and rice. Azadirachtin is one of the most effective botanical insecticides. The effect of azadirachtin against S. frugiperda remains to be determined. Here we report strong growth inhibition of azadirachtin on S. frugiperda larvae under either 1.0 or 5.0 μg/g azadirachtin. To explore the relevant mechanisms, the larvae fed with normal artificial diet and with 1.0 μg/g azadirachtin exposure for 3 days were collected as samples for RNA-Seq. RNA-Seq on S. frugiperda larvae under different treatments identified a total of 24,153 unigenes, including 3494 novel genes, were identified. Among them, 1282 genes were affected by 1.0 μg/g azadirachtin exposure, with 672 up-regulated and 610 down-regulated. The impacted genes include 61 coding for detoxification enzymes (31 P450s, 7 GSTs, 11 CarEs, 7 UGTs and 5 ABC transporters), 31 for cuticle proteins, and several proteins involved in insect chitin and hormone biosynthesis. Our results indicated that azadirachtin could regulate the growth of S. frugiperda by affecting insect chitin and hormone biosynthesis pathway. The enhanced expression of detoxification enzymes might be related to detoxifying azadirachtin. These findings provided a foundation for further delineating the molecular mechanism of growth regulation induced by azadirachtin in S. frugiperda larvae.

Toxicity and Sublethal Effect of Farnesyl Acetate on Diamondback Moth, Plutella xylostella (L.) (Lepidoptera: Plutellidae)

The diamondback moth, Plutella xylostella (L.) (Lepidoptera: Plutellidae), is the most important pest of cruciferous vegetables worldwide. In this study, we evaluated the properties of selected farnesyl derivative compounds against P. xylostella. The toxicity and sublethal concentration (LC50) of farnesyl acetate, farnesyl acetone, farnesyl bromide, farnesyl chloride, and hexahydrofarnesyl acetone were investigated for 96 h. The leaf-dip bioassays showed that farnesyl acetate had a high level of toxicity against P. xylostella compared to other tested farnesyl derivatives. The LC50 value was 56.41 mg/L on the second-instar larvae of P. xylostella. Then, the sublethal effects of farnesyl acetate on biological parameters of P. xylostella were assessed. Compared to the control group, the sublethal concentration of farnesyl acetate decreased pupation and emergence rates, pupal weight, fecundity, egg hatching rate, female ratio, and oviposition period. Furthermore, the developmental time of P. xylostella was extended after being exposed to farnesyl acetate. Moreover, the application of farnesyl acetate on P. xylostella induced morphogenetic abnormalities in larval-pupal intermediates, adults that emerged with twisted wings, or complete adults that could not emerge from the cocoon. These results suggested that farnesyl acetate was highly effective against P. xylostella. The sublethal concentration of farnesyl acetate could reduce the population of P. xylostella by increasing abnormal pupal and adults, and by delaying its development period.

The neonicotinoid thiacloprid causes transcriptional alteration of genes associated with mitochondria at environmental concentrations in honey bees

Thiacloprid is widely used in agriculture and may affect pollinators. However, its molecular effects are poorly known. Here, we report the global gene expression profile in the brain of honey bee foragers assessed by RNA-sequencing. Bees were exposed for 72 h to nominal concentrations of 25 and 250 ng/bee via sucrose solution. Determined residue concentrations by LC-MS/MS were 0.59 and 5.49 ng/bee, respectively. Thiacloprid exposure led to 5 and 71 differentially expressed genes (DEGs), respectively. Nuclear genes encoding mitochondrial ribosomal proteins and enzymes involved in oxidative phosphorylation, as well as metabolism enzymes and transporters were altered at 5.49 ng/bee. Kyoto Encylopedia of Genes and Genomes (KEGG) analysis revealed that mitochondrial ribosome proteins, mitochondrial oxidative phosphorylation, pyrimidine, nicotinate and nicotinamide metabolism and additional metabolic pathways were altered. Among 21 genes assessed by RT-qPCR, the transcript of farnesol dehydrogenase involved in juvenile hormone III synthesis was significantly down-regulated. Transcripts of cyp6a14-like and apolipophorin-II like protein, cytochrome oxidase (cox17) and the non-coding RNA (LOC102654625) were significantly up-regulated at 5.49 ng/bee. Our findings indicate that thiacloprid causes transcriptional changes of genes prominently associated with mitochondria, particularly oxidative phosphorylation. This highlight potential effects of this neonicotinoid on energy metabolism, which may compromise bee foraging and thriving populations at environmentally relevant concentrations.

Comparative Transcriptomics Reveals Gene Families Associated with Predatory Behavior in Photuris femme fatale Fireflies

Identifying the basis of phenotypic variation is a key objective of genetics. This work has been mostly limited to model systems with a plethora of genetic manipulation and functional characterization tools. With the development of high-throughput sequencing and new computational tools, it is possible to identify candidate genes related to phenotypic variation in non-model organisms. Fireflies are excellent for studying phenotypic variation because of their diverse and well-characterized behaviors. Most adult fireflies emit a single mating flash pattern and do not eat. In contrast, adult females of many species in the genus Photuris employ multiple flash patterns and prey upon mate-seeking males of other firefly species. To investigate the genetic basis for this variation, we used comparative transcriptomics to identify positively selected genes between a predatory firefly, Photuris sp., and a non-predatory relative, Photuris frontalis, controlling for genes generally under selection in fireflies by comparing to a Photinus firefly. Nine gene families were identified under positive selection in the predatory versus non-predatory Photuris comparison, including genes involved in digestion, detoxification, vision, reproduction, and neural processes. These results generate intriguing hypotheses about the genetic basis for insect behavior and highlight the utility of comparative transcriptomic tools to investigate complex behaviors in non-model systems.

The juvenile hormone described in Rhodnius prolixus by Wigglesworth is juvenile hormone III skipped bisepoxide

Juvenile hormones (JHs) are sesquiterpenoids synthesized by the corpora allata (CA). They play critical roles during insect development and reproduction. The first JH was described in 1934 as a “metamorphosis inhibitory hormone” in Rhodnius prolixus by Sir Vincent B. Wigglesworth. Remarkably, in spite of the importance of R. prolixus as vectors of Chagas disease and model organisms in insect physiology, the original JH that Wigglesworth described for the kissing-bug R. prolixus remained unidentified. We employed liquid chromatography mass spectrometry to search for the JH homologs present in the hemolymph of fourth instar nymphs of R. prolixus. Wigglesworth’s original JH is the JH III skipped bisepoxide (JHSB3), a homolog identified in other heteropteran species. Changes in the titer of JHSB3 were studied during the 10-day long molting cycle of 4^th instar nymph, between a blood meal and the ecdysis to 5^th instar. In addition we measured the changes of mRNA levels in the CA for the 13 enzymes of the JH biosynthetic pathway during the molting cycle of 4^th instar. Almost 90 years after the first descriptions of the role of JH in insects, this study finally reveals that the specific JH homolog responsible for Wigglesworth’s original observations is JHSB3.

A differential transcriptional profile by Culex quinquefasciatus larvae resistant to Lysinibacillus sphaericus IAB59 highlights genes and pathways associated with the resistance phenotype

Background

The study of the mechanisms by which larvae of the Culex quinquefasciatus mosquito survive exposure to the entomopathogen Lysinibacillus sphaericus has benefited substantially from the generation of laboratory-selected colonies resistant to this bacterium. One such colony, RIAB59, was selected after regular long-term exposure of larvae to the L. sphaericus IAB59 strain. This strain is characterized by its ability to produce the well known Binary (Bin) toxin, and the recently characterized Cry48Aa/Cry49Aa toxin, able to kill Bin-resistant larvae. Resistance to Bin is associated with the depletion of its receptor, Cqm1 α-glucosidase, from the larvae midgut. This study aimed to identify novel molecules and pathways associated with survival of the RIAB59 larvae and the resistance phenotype.

Methods

A transcriptomic approach and bioinformatic tools were used to compare the profiles derived from the midguts of larvae resistant and susceptible to L. sphaericus IAB59.

Results

The RNA-seq profiles identified 1355 differentially expressed genes (DEGs), with 673 down- and 682 upregulated transcripts. One of the most downregulated DEGs was cqm1, which validates the approach. Other strongly downregulated mRNAs encode the enzyme pantetheinase, apolipoprotein D, lipases, heat-shock proteins and a number of lesser known and hypothetical polypeptides. Among the upregulated DEGs, the top most encodes a peroxisomal enzyme involved in lipid metabolism, while others encode enzymes associated with juvenile hormone synthesis, ion channels, DNA binding proteins and defense polypeptides. Further analyses confirmed a strong downregulation of several enzymes involved in lipid catabolism while the assignment of DEGs into metabolic pathways highlighted the upregulation of those related to DNA synthesis and maintenance, confirmed by their clustering into related protein networks. Several other pathways were also identified with mixed profiles of down- and upregulated transcripts. Quantitative RT-PCR confirmed the changes in levels seen for selected mRNAs.

Conclusions

Our transcriptome-wide dataset revealed that the RIAB59 colony, found to be substantially more resistant to Bin than to the Cry48Aa/Cry49Aa toxin, developed a differential expression profile as well as metabolic features co-selected during the long-term adaptation to IAB59 and that are most likely linked to Bin resistance.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3661-y) contains supplementary material, which is available to authorized users.

Reconstruction of insect hormone pathways in an aquatic firefly, Sclerotia aquatilis (Coleoptera: Lampyridae), using RNA-seq

Insect hormones: ecdysteroids and juvenile hormones have crucial functions during the regulation of different developmental pathways in insects. Insect metamorphosis is one of the primary pathways regulated by these hormones. The insect hormone biosynthetic pathway is conserved among arthropods, including insects, with some variations in the form of hormones used among each group of insects. In this study, the candidate genes involved in the insect hormone pathways and their functional roles were assessed in an aquatic firefly, Sclerotia aquatilis using a high-throughput RNA sequencing technique. Illumina next-generation sequencing (NGS) was used to generate transcriptome data for the different developmental stages (i.e., larva, pupa, and adult) of S. aquatilis. A total of 82,022 unigenes were generated across all different developmental stages. Functional annotation was performed for each gene, based on multiple biological databases, generating 46,230 unigenes. These unigenes were subsequently mapped using KEGG pathways. Accordingly, 221 protein-encoding genes involved in the insect hormone pathways were identified, including, JHAMT, CYP15A1, JHE, and Halloween family genes. Twenty potential gene candidates associated with the biosynthetic and degradation pathways for insect hormones were subjected to real-time PCR, reverse transcriptase PCR (RT-PCR) and sequencing analyses. The real-time PCR results showed similar expression patterns as those observed for transcriptome expression profiles for most of the examined genes. RT-PCR and Sanger sequencing confirmed the expressed coding sequences of these gene candidates. This study is the first to examine firefly insect hormone pathways, facilitating a better understanding of firefly growth and development.

Varroa destructor parasitism has a greater effect on proteome changes than the deformed wing virus and activates TGF-β signaling pathways

Honeybee workers undergo metamorphosis in capped cells for approximately 13 days before adult emergence. During the same period, Varroa mites prick the defenseless host many times. We sought to identify proteome differences between emerging Varroa-parasitized and parasite-free honeybees showing the presence or absence of clinical signs of deformed wing virus (DWV) in the capped cells. A label-free proteomic analysis utilizing nanoLC coupled with an Orbitrap Fusion Tribrid mass spectrometer provided a quantitative comparison of 2316 protein hits. Redundancy analysis (RDA) showed that the combination of Varroa parasitism and DWV clinical signs caused proteome changes that occurred in the same direction as those of Varroa alone and were approximately two-fold higher. Furthermore, proteome changes associated with DWV signs alone were positioned above Varroa in the RDA. Multiple markers indicate that Varroa activates TGF-β-induced pathways to suppress wound healing and the immune response and that the collective action of stressors intensifies these effects. Furthermore, we indicate JAK/STAT hyperactivation, p53-BCL-6 feedback loop disruption, Wnt pathway activation, Wnt/Hippo crosstalk disruption, and NF-κB and JAK/STAT signaling conflict in the Varroa–honeybee–DWV interaction. These results illustrate the higher effect of Varroa than of DWV at the time of emergence. Markers for future research are provided.

Epigenetic modifications acetylation and deacetylation play important roles in juvenile hormone action

BACKGROUND

Epigenetic modifications including DNA methylation and post-translational modifications of histones are known to regulate gene expression. Antagonistic activities of histone acetyltransferases (HATs) and histone deacetylases (HDACs) mediate transcriptional reprogramming during insect development as shown in Drosophila melanogaster and other insects. Juvenile hormones (JH) play vital roles in the regulation of growth, development, metamorphosis, reproduction and other physiological processes. However, our current understanding of epigenetic regulation of JH action is still limited. Hence, we studied the role of CREB binding protein (CBP, contains HAT domain) and Trichostatin A (TSA, HDAC inhibitor) on JH action.

RESULTS

Exposure of Tribolium castaneum cells (TcA cells) to JH or TSA caused an increase in expression of Kr-h1 (a known JH-response gene) and 31 or 698 other genes respectively. Knockdown of the gene coding for CBP caused a decrease in the expression of 456 genes including Kr-h1. Interestingly, the expression of several genes coding for transcription factors, nuclear receptors, P450 and fatty acid synthase family members that are known to mediate JH action were affected by CBP knockdown or TSA treatment.

CONCLUSIONS

These data suggest that acetylation and deacetylation mediated by HATs and HDACs play an important role in JH action.

RNA interference against the putative insulin receptor substrate gene chico affects metamorphosis in Leptinotarsa decemlineata

It is noted that insect insulin/insulin-like growth factor/target of rapamycin signaling is critical for the regulation of metamorphosis in holometabolous insects. However, the molecular mechanism remains undetermined. Our previous findings reveal that RNA interference (RNAi)-mediated knockdown of an insulin gene (LdILP2) in Leptinotarsa decemlineata disturbs both 20-hydroxyecdysone (20E) and juvenile hormone (JH) signaling, and impairs pupation. In the present paper, we further observed that the expression of the insulin receptor substrate gene chico (Ldchico) and the phosphoinositide-3-kinase gene pi3k (Ldpi3k92E) was repressed in LdILP2 depleted larvae. Moreover, RNAi of Ldchico or Ldpi3k92E decreased food consumption, affected absorption and metabolism of amino acids and sugars, and reduced expression of several 20E (LdEcR, LdHR3 and LdE75) and JH (LdJHAMT, LdKr-h1 and LdHairy) signaling genes. As a result, larval development was postponed and larval growth was inhibited. Intriguingly, knockdown of Ldchico, rather than Ldpi3k92E, impaired larval-pupal and pupal-adult ecdysis, and specifically repressed transcription of another 20E signaling gene LdUSP. Ingestion of 20E rescued the expression of LdEcR, LdHR3 and LdE75, whereas 20E feeding restored neither the decreased LdUSP mRNA level, nor the reduced pupation and adult emergence rates in Ldchico RNAi larvae. Therefore, Chico is critical for the regulation of larval-pupal-adult transition by a PI3K-independent pathway, perhaps through activation of USP in L. decemlineata.

Ecdysteroid and juvenile hormone biosynthesis, receptors and their signaling in the freshwater microcrustacean Daphnia

The two essential insect hormones, ecdysteroids and juvenile hormones, are possessed not only by insects, but also widely by arthropods, and regulate various developmental and physiological processes. In contrast to the abundant information about molecular endocrine mechanisms in insects, the knowledge of non-insect arthropod endocrinology is still limited. In this review, we summarize recent reports about the molecular basis of these two major insect hormones in the freshwater microcrustacean Daphnia, a keystone taxon in limnetic ecology and a bioindicator in environmental studies. Comprehensive comparisons of endocrine signaling pathways between insects and daphnids may shed light on the regulatory mechanisms of various biological phenomena and, moreover, evolutionary processes of arthropod species.

An Alcohol Dehydrogenase from the Short-Chain Dehydrogenase/Reductase Family of Enzymes for the Lactonization of Hexane-1,6-diol

Biocatalytic production of lactones, and in particular ϵ-caprolactone (CL), have gained increasing interest as a greener route to polymer building blocks, especially through the use of Baeyer-Villiger monooxygenases (BVMOs). Despite several advances in the field, BVMOs, however, still suffer several practical limitations. Alcohol dehydrogenase (ADH)-mediated lactonization of diols in turn has received far less attention and very few enzymes have been identified for the conversion of diols to lactones, with horse-liver ADH (HLADH) remaining the catalyst of choice. Screening of a diverse panel of ADHs, AaSDR-1, a member of the short-chain dehydrogenase/reductase family, was found to produce ϵ-caprolactone from hexane-1,6-diol. Moreover, cofactor regeneration by an NADH oxidase eliminated the requirement of co-substrates, yielding water as the sole by-product. Despite lower turnover frequencies as compared to HLADH, higher selectivity was found for the production of CL, with HLADH forming significant amounts of 6-hydroxyhexanoic acid and adipic acid through aldehyde dehydrogenation/oxidation of the gem-diol intermediates. Also, CL yield were shown to be dependent on buffer choice, as structural elucidation of a Tris adduct confirmed the buffer amine to react with aliphatic aldehydes forming a Schiff-base intermediate which through further ADH oxidation, forms a tricyclic acetal product.

Omics approaches to study juvenile hormone synthesis

The juvenile hormones (JHs) are a family of insect acyclic sesquiterpenoids produced by the corpora allata (CA), a pair of endocrine glands connected to the brain. They are involved in the regulation of development, reproduction, behavior, caste determination, diapause, stress response, and numerous polyphenisms. In the post-genomics era, comprehensive analyses using functional 'omics' technologies such as transcriptomics, proteomics and metabolomics have increased our understanding of the activity of the minute CA. This review attempts to summarize some of the 'omics' studies that have contributed to further understand JH synthesis in insects, with an emphasis on our own research on the mosquito Aedes aegypti.

RNA-seq analysis of gene expression changes during pupariation in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

Background

The oriental fruit fly, Bactrocera dorsalis (Hendel) has been considered to be one of the most important agricultural pest around the world. As a holometabolous insect, larvae must go through a metamorphosis process with dramatic morphological and structural changes to complete their development. To better understand the molecular mechanisms of these changes, RNA-seq of B. dorsalis from wandering stage (WS), late wandering stage (LWS) and white puparium stage (WPS) were performed.

Results

In total, 11,721 transcripts were obtained, out of which 1914 genes (578 up-regulated and 1336 down-regulated) and 2047 genes (655 up-regulated and 1392 down-regulated) were found to be differentially expressed between WS and LWS, as well as between WS and WPS, respectively. Of these DEGs, 1862 and 1996 genes were successfully annotated in various databases. The analysis of RNA-seq data together with qRT-PCR validation indicated that during this transition, the genes in the oxidative phosphorylation pathway, and genes encoding P450s, serine protease inhibitor, and cuticular proteins were down-regulated, while the serine protease genes were up-regulated. Moreover, we found some 20-hydroxyecdysone (20E) biosynthesis and signaling pathway genes had a higher expression in the WS, while the genes responsible for juvenile hormone (JH) synthesis, degradation, signaling and transporter pathways were down-regulated, suggesting these genes might be involved in the process of larval pupariation in B. dorsalis. For the chitinolytic enzymes, the genes encoding chitinases (chitinase 2, chitinase 5, chitinase 8, and chitinase 10) and chitin deacetylase might play the crucial role in the degradation of insect chitin with their expressions significantly increased during the transition. Here, we also found that chitin synthase 1A might be involved in the chitin synthesis of cuticles during the metamorphosis in B. dorsalis.

Conclusions

Significant changes at transcriptional level were identified during the larval pupariation of B. dorsalis. Importantly, we also obtained a vast quantity of RNA-seq data and identified metamorphosis associated genes, which would all help us to better understand the molecular mechanism of metamorphosis process in B. dorsalis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5077-z) contains supplementary material, which is available to authorized users.

Assembling large genomes: analysis of the stick insect (Clitarchus hookeri) genome reveals a high repeat content and sex-biased genes associated with reproduction

BACKGROUND

Stick insects (Phasmatodea) have a high incidence of parthenogenesis and other alternative reproductive strategies, yet the genetic basis of reproduction is poorly understood. Phasmatodea includes nearly 3000 species, yet only the genome of Timema cristinae has been published to date. Clitarchus hookeri is a geographical parthenogenetic stick insect distributed across New Zealand. Sexual reproduction dominates in northern habitats but is replaced by parthenogenesis in the south. Here, we present a de novo genome assembly of a female C. hookeri and use it to detect candidate genes associated with gamete production and development in females and males. We also explore the factors underlying large genome size in stick insects.

RESULTS

The C. hookeri genome assembly was 4.2 Gb, similar to the flow cytometry estimate, making it the second largest insect genome sequenced and assembled to date. Like the large genome of Locusta migratoria, the genome of C. hookeri is also highly repetitive and the predicted gene models are much longer than those from most other sequenced insect genomes, largely due to longer introns. Miniature inverted repeat transposable elements (MITEs), absent in the much smaller T. cristinae genome, is the most abundant repeat type in the C. hookeri genome assembly. Mapping RNA-Seq reads from female and male gonadal transcriptomes onto the genome assembly resulted in the identification of 39,940 gene loci, 15.8% and 37.6% of which showed female-biased and male-biased expression, respectively. The genes that were over-expressed in females were mostly associated with molecular transportation, developmental process, oocyte growth and reproductive process; whereas, the male-biased genes were enriched in rhythmic process, molecular transducer activity and synapse. Several genes involved in the juvenile hormone synthesis pathway were also identified.

CONCLUSIONS

The evolution of large insect genomes such as L. migratoria and C. hookeri genomes is most likely due to the accumulation of repetitive regions and intron elongation. MITEs contributed significantly to the growth of C. hookeri genome size yet are surprisingly absent from the T. cristinae genome. Sex-biased genes identified from gonadal tissues, including genes involved in juvenile hormone synthesis, provide interesting candidates for the further study of flexible reproduction in stick insects.

Modeling the flux of metabolites in the juvenile hormone biosynthesis pathway using generalized additive models and ordinary differential equations

Juvenile hormone (JH) regulates development and reproductive maturation in insects. The corpora allata (CA) from female adult mosquitoes synthesize fluctuating levels of JH, which have been linked to the ovarian development and are influenced by nutritional signals. The rate of JH biosynthesis is controlled by the rate of flux of isoprenoids in the pathway, which is the outcome of a complex interplay of changes in precursor pools and enzyme levels. A comprehensive study of the changes in enzymatic activities and precursor pool sizes have been previously reported for the mosquito Aedes aegypti JH biosynthesis pathway. In the present studies, we used two different quantitative approaches to describe and predict how changes in the individual metabolic reactions in the pathway affect JH synthesis. First, we constructed generalized additive models (GAMs) that described the association between changes in specific metabolite concentrations with changes in enzymatic activities and substrate concentrations. Changes in substrate concentrations explained 50% or more of the model deviances in 7 of the 13 metabolic steps analyzed. Addition of information on enzymatic activities almost always improved the fitness of GAMs built solely based on substrate concentrations. GAMs were validated using experimental data that were not included when the model was built. In addition, a system of ordinary differential equations (ODE) was developed to describe the instantaneous changes in metabolites as a function of the levels of enzymatic catalytic activities. The results demonstrated the ability of the models to predict changes in the flux of metabolites in the JH pathway, and can be used in the future to design and validate experimental manipulations of JH synthesis.

Novel NAD+-Farnesal Dehydrogenase from Polygonum minus Leaves. Purification and Characterization of Enzyme in Juvenile Hormone III Biosynthetic Pathway in Plant

Juvenile Hormone III is of great concern due to negative effects on major developmental and reproductive maturation in insect pests. Thus, the elucidation of enzymes involved JH III biosynthetic pathway has become increasing important in recent years. One of the enzymes in the JH III biosynthetic pathway that remains to be isolated and characterized is farnesal dehydrogenase, an enzyme responsible to catalyze the oxidation of farnesal into farnesoic acid. A novel NAD⁺-farnesal dehydrogenase of Polygonum minus was purified (315-fold) to apparent homogeneity in five chromatographic steps. The purification procedures included Gigacap S-Toyopearl 650M, Gigacap Q-Toyopearl 650M, and AF-Blue Toyopearl 650ML, followed by TSK Gel G3000SW chromatographies. The enzyme, with isoelectric point of 6.6 is a monomeric enzyme with a molecular mass of 70 kDa. The enzyme was relatively active at 40°C, but was rapidly inactivated above 45°C. The optimal temperature and pH of the enzyme were found to be 35°C and 9.5, respectively. The enzyme activity was inhibited by sulfhydryl agent, chelating agent, and metal ion. The enzyme was highly specific for farnesal and NAD⁺. Other terpene aldehydes such as trans- cinnamaldehyde, citral and α- methyl cinnamaldehyde were also oxidized but in lower activity. The K_m values for farnesal, citral, trans- cinnamaldehyde, α- methyl cinnamaldehyde and NAD⁺ were 0.13, 0.69, 0.86, 1.28 and 0.31 mM, respectively. The putative P. minus farnesal dehydrogenase that’s highly specific towards farnesal but not to aliphatic aldehydes substrates suggested that the enzyme is significantly different from other aldehyde dehydrogenases that have been reported. The MALDI-TOF/TOF-MS/MS spectrometry further identified two peptides that share similarity to those of previously reported aldehyde dehydrogenases. In conclusion, the P. minus farnesal dehydrogenase may represent a novel plant farnesal dehydrogenase that exhibits distinctive substrate specificity towards farnesal. Thus, it was suggested that this novel enzyme may be functioning specifically to oxidize farnesal in the later steps of JH III pathway. This report provides a basic understanding for recombinant production of this particular enzyme. Other strategies such as adding His-tag to the protein makes easy the purification of the protein which is completely different to the native protein. Complete sequence, structure and functional analysis of the enzyme will be important for developing insect-resistant crop plants by deployment of transgenic plant.

Identification of Multiple Loci Associated with Social Parasitism in Honeybees

In colonies of the honeybee Apis mellifera, the queen is usually the only reproductive female, which produces new females (queens and workers) by laying fertilized eggs. However, in one subspecies of A. mellifera, known as the Cape bee (A. m. capensis), worker bees reproduce asexually by thelytoky, an abnormal form of meiosis where two daughter nucleii fuse to form single diploid eggs, which develop into females without being fertilized. The Cape bee also exhibits a suite of phenotypes that facilitate social parasitism whereby workers lay such eggs in foreign colonies so their offspring can exploit their resources. The genetic basis of this switch to social parasitism in the Cape bee is unknown. To address this, we compared genome variation in a sample of Cape bees with other African populations. We find genetic divergence between these populations to be very low on average but identify several regions of the genome with extreme differentiation. The regions are strongly enriched for signals of selection in Cape bees, indicating that increased levels of positive selection have produced the unique set of derived phenotypic traits in this subspecies. Genetic variation within these regions allows unambiguous genetic identification of Cape bees and likely underlies the genetic basis of social parasitism. The candidate loci include genes involved in ecdysteroid signaling and juvenile hormone and dopamine biosynthesis, which may regulate worker ovary activation and others whose products localize at the centrosome and are implicated in chromosomal segregation during meiosis. Functional analysis of these loci will yield insights into the processes of reproduction and chemical signaling in both parasitic and non-parasitic populations and advance understanding of the process of normal and atypical meiosis.

Knockdown of juvenile hormone acid methyl transferase severely affects the performance of Leptinotarsa decemlineata (Say) larvae and adults

BACKGROUND

Juvenile hormone (JH) plays a critical role in the regulation of metamorphosis in Leptinotarsa decemlineata, a notorious defoliator of potato. JH acid methyltransferase (JHAMT) is involved in one of the final steps of JH biosynthesis.

RESULTS

A putative JHAMT cDNA (LdJHAMT) was cloned. Two double-stranded RNAs (dsRNAs) (dsJHAMT1 and dsJHAMT2) against LdJHAMT were constructed and bacterially expressed. Experiments were conducted to investigate the effectiveness of RNAi in both second- and fourth-instar larvae. Dietary introduction of dsJHAMT1 and dsJHAMT2 successfully knocked down the target gene, lowered JH titre in the haemolymph and reduced the transcript of Krüppel homologue 1 gene. Ingestion of dsJHAMT caused larval death and weight loss, shortened larval developmental period and impaired pupation. Moreover, the dsJHAMT-fed pupae exhibited lower adult emergence rates. The resulting adults weighed an average of 50 mg less than the control group, and the females did not deposit eggs. Application of pyriproxyfen to the dsJHAMT-fed insects rescued all the negative effects.

CONCLUSIONS

LdJHAMT expresses functional JHAMT enzyme. The RNAi targeting LdJHAMT could be used for control of L. decemlineata. © 2015 Society of Chemical Industry.

Exploring complex pheromone biosynthetic processes in the bumblebee male labial gland by RNA sequencing

Male marking pheromones (MPs) are used by the majority of bumblebee species (Hymenoptera: Apidae), including a commercially important greenhouse pollinator, the buff-tailed bumblebee (Bombus terrestris), to attract conspecific females. MP biosynthetic processes in the cephalic part of the bumblebee male labial gland (LG) are of extraordinary complexity, involving enzymes of fatty acid and isoprenoid biosynthesis, which jointly produce more than 50 compounds. We employed a differential transcriptomic approach to identify candidate genes involved in MP biosynthesis by sequencing Bombus terrestris LG and fat body (FB) transcriptomes. We identified 12 454 abundantly expressed gene products (reads per kilobase of exon model per million mapped reads value > 1) that had significant hits in the GenBank nonredundant database. Of these, 876 were upregulated in the LG (> 4-fold difference). We identified more than 140 candidate genes potentially involved in MP biosynthesis, including esterases, fatty acid reductases, lipases, enzymes involved in limited fatty acid chain shortening, neuropeptide receptors and enzymes involved in biosynthesis of triacylglycerols, isoprenoids and fatty acids. For selected candidates, we confirmed their abundant expression in LG using quantitative real-time reverse transcription-PCR (qRT-PCR). Our study shows that the Bombus terrestris LG transcriptome reflects both fatty acid and isoprenoid MP biosynthetic processes and identifies rational gene targets for future studies to disentangle the molecular basis of MP biosynthesis. Additionally, LG and FB transcriptomes enrich the available transcriptomic resources for Bombus terrestris.

Comparative Transcriptomics Reveals Key Gene Expression Differences between Diapausing and Non-Diapausing Adults of Culex pipiens

Diapause is a critical eco-physiological adaptation for winter survival in the West Nile Virus vector, Culex pipiens, but little is known about the molecular mechanisms that distinguish diapause from non-diapause in this important mosquito species. We used Illumina RNA-seq to simultaneously identify and quantify relative transcript levels in diapausing and non-diapausing adult females. Among 65,623,095 read pairs, we identified 41 genes with significantly different transcript abundances between these two groups. Transcriptome divergences between these two phenotypes include genes related to juvenile hormone synthesis, anaerobic metabolism, innate immunity and cold tolerance.

Mevalonate-Farnesal Biosynthesis in Ticks: Comparative Synganglion Transcriptomics and a New Perspective

Juvenile hormone (JH) controls the growth, development, metamorphosis, and reproduction of insects. For many years, the general assumption has been that JH regulates tick and other acarine development and reproduction the same as in insects. Although researchers have not been able to find the common insect JHs in hard and soft tick species and JH applications appear to have no effect on tick development, it is difficult to prove the negative or to determine whether precursors to JH are made in ticks. The tick synganglion contains regions which are homologous to the corpora allata, the biosynthetic source for JH in insects. Next-gen sequencing of the tick synganglion transcriptome was conducted separately in adults of the American dog tick, Dermacentor variabilis, the deer tick, Ixodes scapularis, and the relapsing fever tick, Ornithodoros turicata as a new approach to determine whether ticks can make JH or a JH precursor. All of the enzymes that make up the mevalonate pathway from acetyl-CoA to farnesyl diphosphate (acetoacetyl-CoA thiolase, HMG-S, HMG-R, mevalonate kinase, phosphomevalonate kinase, diphosphomevalonate decarboxylase, and farnesyl diphosphate synthase) were found in at least one of the ticks studied but most were found in all three species. Sequence analysis of the last enzyme in the mevalonate pathway, farnesyl diphosphate synthase, demonstrated conservation of the seven prenyltransferase regions and the aspartate rich motifs within those regions typical of this enzyme. In the JH branch from farnesyl diphosphate to JH III, we found a putative farnesol oxidase used for the conversion of farnesol to farnesal in the synganglion transcriptome of I. scapularis and D. variabilis. Methyltransferases (MTs) that add a methyl group to farnesoic acid to make methyl farnesoate were present in all of the ticks studied with similarities as high as 36% at the amino acid level to insect JH acid methyltransferase (JHAMT). However, when the tick MTs were compared to the known insect JHAMTs from several insect species at the amino acid level, the former lacked the farnesoic acid binding motif typical in insects. The P450s shown in insects to add the C10,11 epoxide to methyl farnesoate, are in the CYP15 family; this family was absent in our tick transcriptomes and in the I. scapularis genome, the only tick genome available. These data suggest that ticks do not synthesize JH III but have the mevalonate pathway and may produce a JH III precursor.

High substrate specificity of ipsdienol dehydrogenase (IDOLDH), a short-chain dehydrogenase from Ips pini bark beetles

Ips spp. bark beetles use ipsdienol, ipsenol, ipsdienone and ipsenone as aggregation pheromone components and pheromone precursors. For Ips pini, the short-chain oxidoreductase ipsdienol dehydrogenase (IDOLDH) converts (-)-ipsdienol to ipsdienone, and thus likely plays a role in determining pheromone composition. In order to further understand the role of IDOLDH in pheromone biosynthesis, we compared IDOLDH to its nearest functionally characterized ortholog with a solved structure: human L-3-hydroxyacyl-CoA dehydrogenase type II/ amyloid-β binding alcohol dehydrogenase (hHADH II/ABAD), and conducted functional assays of recombinant IDOLDH to determine substrate and product ranges and structural characteristics. Although IDOLDH and hHADH II/ABAD had only 35% sequence identity, their predicted tertiary structures had high identity. We found IDOLDH is a functional homo-tetramer. In addition to oxidizing (-)-ipsdienol, IDOLDH readily converted racemic ipsenol to ipsenone, and stereo-specifically reduced both ketones to their corresponding (-)-alcohols. The (+)-enantiomers were never observed as products. Assays with various substrate analogs showed IDOLDH had high substrate specificity for (-)-ipsdienol, ipsenol, ipsenone and ipsdienone, supporting that IDOLDH functions as a pheromone-biosynthetic enzyme. These results suggest that different IDOLDH orthologs and or activity levels contribute to differences in Ips spp. pheromone composition.

Nuclear receptor ecdysone-induced protein 75 is required for larval-pupal metamorphosis in the Colorado potato beetle Leptinotarsa decemlineata (Say)

20-hydroxyecdysone (20E) and juvenile hormone (JH) are key regulators of insect development. In this study, three Leptinotarsa decemlineata Ecdysone-induced protein 75 (LdE75) cDNAs (LdE75A, B and C) were cloned from L. decemlineata. The three LdE75 isoforms were highly expressed just before or right after each moult. Within the fourth larval instar, they showed a small rise and a big peak 40 and 80 h after ecdysis. The expression peaks of the three LdE75s coincided with the peaks of circulating 20E levels. In vitro midgut culture and in vivo bioassay revealed that 20E and an ecdysteroid agonist halofenozide (Hal) enhanced LdE75 expression in the day 1 final larval instars. Conversely, a decrease in 20E by feeding a double-stranded RNA (dsRNA) against an ecdysteroidogenesis gene, Shade (LdSHD), repressed the expression of LdE75. Moreover, Hal upregulated the expression of the three LdE75s in LdSHD-silenced larvae. Thus, 20E pulses activate the transcription of LdE75s. Furthermore, ingesting dsE75-1 and dsE75-2 from a common fragment of the three isoforms successfully knocked down these LdE75s, and caused developmental arrest. Finally, knocking down LdE75s significantly repressed the transcription of three ecdysteroidogenesis genes, lowered the 20E titre and affected the expression of two 20E-response genes. Silencing LdE75s also induced the expression of a JH biosynthesis gene, increased JH titre and activated the transcription of a JH early-inducible gene. Thus, Ld E75s are required for larval-pupal metamorphosis and act mainly by modulating 20E and JH titres and mediating their signalling pathways.

RNA interference suppression of the receptor tyrosine kinase Torso gene impaired pupation and adult emergence in Leptinotarsa decemlineata

In Drosophila melanogaster prothoracic gland (PG) cells, Torso mediates prothoracicotropic hormone (PTTH)-triggered mitogen activated protein kinase (MAPK) pathway (consisting of four core components Ras, Raf, MEK and ERK) to stimulate ecdysteroidogenesis. In this study, LdTorso, LdRas, LdRaf and LdERK were cloned in Leptinotarsa decemlineata. The four genes were highly or moderately expressed in the larval prothoracic glands. At the first- to third-instar stages, their expression levels were higher just before and right after the molt, and were lower in the mid instars. At the fourth-instar stage, their transcript levels were higher before prepupal stage. RNA interference-mediated knockdown of LdTorso delayed larval development, increased pupal weight, and impaired pupation and adult emergence. Moreover, knockdown of LdTorso decreased the mRNA levels of LdRas, LdRaf and LdERK, repressed the transcription of two ecdysteroidogenesis genes (LdPHM and LdDIB), lowered 20E titer, and downregulated the expression of several 20E-response genes (LdEcR, LdUSP, LdHR3 and LdFTZ-F1). Furthermore, silencing of LdTorso induced the expression of a JH biosynthesis gene LdJHAMT, increased JH titer, and activated the transcription of a JH early-inducible gene LdKr-h1. Thus, our results suggest that Torso transduces PTTH-triggered MAPK signal to regulate ecdysteroidogenesis in the PGs in a non-drosophiline insect.

Purification and Characterization of a Novel NAD(P)+-Farnesol Dehydrogenase from Polygonum minus Leaves

Juvenile hormones have attracted attention as safe and selective targets for the design and development of environmentally friendly and biorational insecticides. In the juvenile hormone III biosynthetic pathway, the enzyme farnesol dehydrogenase catalyzes the oxidation of farnesol to farnesal. In this study, farnesol dehydrogenase was extracted from Polygonum minus leaves and purified 204-fold to apparent homogeneity by ion-exchange chromatography using DEAE-Toyopearl, SP-Toyopearl, and Super-Q Toyopearl, followed by three successive purifications by gel filtration chromatography on a TSK-gel GS3000SW. The enzyme is a heterodimer comprised of subunits with molecular masses of 65 kDa and 70 kDa. The optimum temperature and pH were 35°C and pH 9.5, respectively. Activity was inhibited by sulfhydryl reagents, metal-chelating agents and heavy metal ions. The enzyme utilized both NAD⁺ and NADP⁺ as coenzymes with K_m values of 0.74 mM and 40 mM, respectively. Trans, trans-farnesol was the preferred substrate for the P. minus farnesol dehydrogenase. Geometrical isomers of trans, trans-farnesol, cis, trans-farnesol and cis, cis-farnesol were also oxidized by the enzyme with lower activity. The K_m values for trans, trans-farnesol, cis, trans-farnesol and cis, cis-farnesol appeared to be 0.17 mM, 0.33 mM and 0.42 mM, respectively. The amino acid sequences of 4 tryptic peptides of the enzyme were analyzed by MALDI-TOF/TOF-MS spectrometry, and showed no significant similarity to those of previously reported farnesol dehydrogenases. These results suggest that the purified enzyme is a novel NAD(P)⁺-dependent farnesol dehydrogenase. The purification and characterization established in the current study will serve as a basis to provide new information for recombinant production of the enzyme. Therefore, recombinant farnesol dehydrogenase may provide a useful molecular tool in manipulating juvenile hormone biosynthesis to generate transgenic plants for pest control.