CYP303A1 has a conserved function in adult eclosion in Locusta migratoria and Drosophila melanogaster

Characterization of CYP303A1 and its potential application based on ZIF-8 nanoparticle-wrapped dsRNA in Nilaparvata lugens (Stål)

BACKGROUND

RNA interference (RNAi) technology has been put forward as a promising method for pest control and resistance management. Mining highly efficient lethal genes and constructing stable double-stranded RNA (dsRNA) delivery systems are of great significance to improve the application potential of RNAi technology.

RESULTS

In this study, we characterized a molting-related gene, NlCYP303A1, in Nilaparvata lugens that was highly expressed in the cuticle and at the end stages of each instar in nymphs. Silencing the expression of NlCYP303A1 in N. lugens resulted in a deformed phenotype and a significant increase in mortality. Furthermore, interfering with NlCYP303A1 changed the relative expression of key genes in the chitin synthesis and degradation pathway. Finally, we used the nanocarrier zeolitic imidazolate framework-8 (ZIF-8) to load dsNlCYP303A1, forming a complex denoted as dsNlCYP303A1@ZIF-8. The results of both feeding and rice-seedling dip experiments indicated that the expression of NlCYP303A1 was dramatically and persistently suppressed by the dsNlCYP303A1@ZIF-8 treatment, compared with treatment with dsNlCYP303A1, suggesting that ZIF-8 can enhance the interference efficiency as well as the stability of dsNlCYP303A1.

CONCLUSIONS

Our results demonstrate that the lethal gene NlCYP303A1 can be employed as an excellent target for RNAi technology by loading onto a nano-delivery system, and provide new insights into the creation of innovative pest control approaches. © 2024 Society of Chemical Industry.

Chromosome-level genome assembly of Scathophaga stercoraria provides new insights into the evolutionary adaptations of dung flies

The yellow dung fly Scathophaga stercoraria is a widely distributed species in high-altitude regions of the Northern Hemisphere. It plays important roles as a decomposer, predator, and pollinator in the ecosystem. As a staple model organism, S. stercoraria serves as a standard test species for assessing the toxicity of drug residues in livestock dung and has been the focus of numerous studies. The genetic mechanisms underlying the ecological adaptability of S. stercoraria remain poorly understood. To fill the gap, we first assembled a high-quality chromosome-level genome of S. stercoraria, resulting in a final assembly size of 549.64 Mb, with a contig N50 of 4.06 Mb, and 92.53 % of the sequence anchored to six chromosomes. Gene family analysis revealed an expansion of Toll (Toll1), GNBP3, Cyp303a1, Cyp4d14, Cyp6g1, OR67d, and yolk protein genes in the S. stercoraria genome. Transcriptome analysis indicated that most genes in the trypsin and carboxypeptidase gene families are predominantly expressed during the larval stage, whereas the α-Amylase gene family is mainly expressed during the adult stage. Additionally, PGRP-SC is highly expressed during the larval stage, OBPs are primarily expressed during the adult stage, and yolk protein genes exhibit female-biased expression. Our study not only provides a new resource for the dung flies genomic pool, but also identifies the expression patterns of key ecologically adaptative genes and gene families at the developmental stages, which provides new insights into the ecological adaptive evolution of dung flies.

CYP303A1 regulates molting and metamorphosis through 20E signaling in Nilaparvata lugens Stål (Hemiptera: Delphacidae)

Cytochrome P450s play a crucial role in the breakdown of external substances and perform important activities in the hormone system of insects. It has been understood that P450s were essential in the metabolism of ecdysteroids. CYP303A1 is a highly conserved CYP in most insects, but its specific physiological functions remain poorly understood in Nilaparvata lugens Stål. In this study, NlCYP303A1 was identified and highly expressed in the pre-molt stages, predominantly in the cuticle-producing tissues. Silencing of NlCYP303A1 caused a lethal phenotype with a molting defect. Moreover, the 20E titers, the expression levels of Halloween genes, and critical genes associated with the 20E signaling pathway in N. lugens nymphs were significantly decreased with the silencing NlCYP303A1. We further performed additional backfilling of 20E to rescue the RNAi effects on NlCYP303A1. The gene expression levels that were previously reduced caused by silencing NlCYP303A1 were significantly elevated. However, the molting defects of nymphs were not effectively improved. The results demonstrated NlCYP303A1 plays a crucial role in the molting and metamorphosis of N. lugens by regulating the 20E signaling pathway and cuticular formation, enhances the understanding of the functional role of CYP 2 clans, and identifies candidate gene for RNAi-based control of N. lugens.

Comparative Genomics Uncovers the Evolutionary Dynamics of Detoxification and Insecticide Target Genes Across 11 Phlebotomine Sand Flies

Sand flies infect more than 1 million people annually with Leishmania parasites and other bacterial and viral pathogens. Progress in understanding sand fly adaptations to xenobiotics has been hampered by the limited availability of genomic resources. To address this gap, we sequenced, assembled, and annotated the transcriptomes of 11 phlebotomine sand fly species. Subsequently, we leveraged these genomic resources to generate novel evolutionary insights pertaining to their adaptations to xenobiotics, including those contributing to insecticide resistance. Specifically, we annotated over 2,700 sand fly detoxification genes and conducted large-scale phylogenetic comparisons to uncover the evolutionary dynamics of the five major detoxification gene families: cytochrome P450s (CYPs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs), carboxyl/cholinesterases (CCEs), and ATP-binding cassette (ABC) transporters. Using this comparative approach, we show that sand flies have evolved diverse CYP and GST gene repertoires, with notable lineage-specific expansions in gene groups evolutionarily related to known xenobiotic metabolizers. Furthermore, we show that sand flies have conserved orthologs of (i) CYP4G genes involved in cuticular hydrocarbon biosynthesis, (ii) ABCB genes involved in xenobiotic toxicity, and (iii) two primary insecticide targets, acetylcholinesterase-1 (Ace1) and voltage gated sodium channel (VGSC). The biological insights and genomic resources produced in this study provide a foundation for generating and testing hypotheses regarding the molecular mechanisms underlying sand fly adaptations to xenobiotics.

Expression and Functional Analysis of Two Cytochrome P450 Monooxygenase Genes and a UDP-Glycosyltransferase Gene Linked with Thiamethoxam Resistance in the Colorado Potato Beetle

Cytochrome P450 monooxygenases (P450s) and UDP-glycosyltransferases (UGTs) are involved in the evolution of insecticide resistance. Leptinotarsa decemlineata (Say), the Colorado potato beetle (CPB), is a notorious insect that has developed resistance to various insecticides including neonicotinoids. This study investigated whether the differentially expressed P450 genes CYP9Z140 and CYP9AY1 and UGT gene UGT321AP1, found in our transcriptome results, conferred resistance to thiamethoxam in L. decemlineata. Resistance monitoring showed that the sampled field populations of L. decemlineata adults collected from Urumqi City and Qapqal, Jimsar, and Mulei Counties of Xinjiang in 2021-2023 developed low levels of resistance to thiamethoxam with resistance ratios ranging from 6.66- to 9.52-fold. Expression analyses indicated that CYP9Z140, CYP9AY1, and UGT321AP1 were significantly upregulated in thiamethoxam-resistant populations compared with susceptible populations. The expression of all three genes also increased significantly after thiamethoxam treatment compared with the control. Spatiotemporal expression patterns showed that the highest expression of CYP9Z140 and CYP9AY1 occurred in pupae and the midgut, whereas UGT321AP1 was highly expressed in adults and Malpighian tubules. Knocking down all three genes individually or simultaneously using RNA interference increased the sensitivity of adult L. decemlineata to thiamethoxam. These results suggest that overexpression of CYP9Z140, CYP9AY1, and UGT321AP1 contributes to the development of thiamethoxam resistance in L. decemlineata and provides a scientific basis for improving new resistance management of CPB.

Cyp6g2 is the major P450 epoxidase responsible for juvenile hormone biosynthesis in Drosophila melanogaster

BACKGROUND

Juvenile hormones (JH) play crucial role in regulating development and reproduction in insects. The most common form of JH is JH III, derived from MF through epoxidation by CYP15 enzymes. However, in the higher dipterans, such as the fruitfly, Drosophila melanogaster, a bis-epoxide form of JHB3, accounted most of the JH detected. Moreover, these higher dipterans have lost the CYP15 gene from their genomes. As a result, the identity of the P450 epoxidase in the JH biosynthesis pathway in higher dipterans remains unknown.

RESULTS

In this study, we show that Cyp6g2 serves as the major JH epoxidase responsible for the biosynthesis of JHB3 and JH III in D. melanogaster. The Cyp6g2 is predominantly expressed in the corpus allatum (CA), concurring with the expression pattern of jhamt, another well-studied gene that is crucial in the last steps of JH biosynthesis. Mutation in Cyp6g2 leads to severe disruptions in larval-pupal metamorphosis and exhibits reproductive deficiencies, exceeding those seen in jhamt mutants. Notably, Cyp6g2-/-::jhamt2 double mutants all died at the pupal stage but could be rescued through the topical application of JH analogs. JH titer analyses revealed that both Cyp6g2-/- mutant and jhamt2 mutant lacking JHB3 and JH III, while overexpression of Cyp6g2 or jhamt caused a significant increase in JHB3 and JH III titer.

CONCLUSIONS

These findings collectively established that Cyp6g2 as the major JH epoxidase in the higher dipterans and laid the groundwork for the further understanding of JH biosynthesis. Moreover, these findings pave the way for developing specific Cyp6g2 inhibitors as insect growth regulators or insecticides.

Nuclear Receptor FTZ-F1 Controls Locust Molt by Regulating the Molting Process of Locusta migratoria

Fushi-tarazu factor 1 (FTZ-F1) is a class of transcription factors belonging to the nuclear receptor superfamily and an important molting regulator in insects; however, its detailed function in the molting process of Locusta migratoria is still unclear. This study identified two FTZ-F1 transcripts (LmFTZ-F1-X1 and LmFTZ-F1-X2) in L. migratoria. The classical domains of FTZ-F1 were present in their protein sequences and distinguished based on their variable N-terminal domains. Reverse-transcription quantitative polymerase chain reaction analysis revealed that LmFTZ-F1-X1 and LmFTZ-F1-X2 were highly expressed in the integument. RNA interference (RNAi) was used to explore the function of LmFTZ-F1s in the molting of the third-instar nymph. Separate LmFTZ-F1-X1 or LmFTZ-F1-X2 silencing did not affect the normal development of third-instar nymphs; however, the simultaneous RNAi of LmFTZ-F1-X1 and LmFTZ-F1-X2 caused the nymphs to be trapped in the third instar stage and finally die. Furthermore, the hematoxylin-eosin and chitin staining of the cuticle showed that the new cuticles were thickened after silencing the LmFTZ-F1s compared to the controls. RNA-seq analysis showed that genes encoding four cuticle proteins, two chitin synthesis enzymes, and cytochrome P450 303a1 were differentially expressed between dsGFP- and dsLmFTZ-F1s-injected groups. Taken together, LmFTZ-F1-X1 and LmFTZ-F1-X2 are involved in the ecdysis of locusts, possibly by regulating the expression of genes involved in cuticle formation, chitin synthesis, and other key molting processes.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Sterol Regulation of Development and 20-Hydroxyecdysone Biosynthetic and Signaling Genes in Drosophila melanogaster

Ecdysteroids are crucial in regulating the growth and development of insects. In the fruit fly Drosophila melanogaster, both C27 and C28 ecdysteroids have been identified. While the biosynthetic pathway of the C27 ecdysteroid 20-hydroxyecdysone (20E) from cholesterol is relatively well understood, the biosynthetic pathway of C28 ecdysteroids from C28 or C29 dietary sterols remains unknown. In this study, we found that different dietary sterols (including the C27 sterols cholesterol and 7-dehydrocholesterol, the C28 sterols brassicasterol, campesterol, and ergosterol, and the C29 sterols β-sitosterol, α-spinasterol, and stigmasterol) differentially affected the expression of 20E biosynthetic genes to varying degrees, but similarly activated 20E primary response gene expression in D. melanogaster Kc cells. We also found that a single dietary sterol was sufficient to support D. melanogaster growth and development. Furthermore, the expression levels of some 20E biosynthetic genes were significantly altered, whereas the expression of 20E signaling primary response genes remained unaffected when flies were reared on lipid-depleted diets supplemented with single sterol types. Overall, our study provided preliminary clues to suggest that the same enzymatic system responsible for the classical C27 ecdysteroid 20E biosynthetic pathway also participated in the conversion of C28 and C29 dietary sterols into C28 ecdysteroids.

Identification and characterization of CYP307A1 as a molecular target for controlling the small hive beetle, Aethina tumida

BACKGROUND

The molting hormone 20-hydroxyecdysone (20E) plays a key role in insect development, metamorphosis, and reproduction. Previous studies have shown that ecdysteroid metabolism is regulated by a series of CYP genes in most of the insect species. However, the roles of these CYP genes in a Coleopteran beetle, Aethina tumida (small hive beetle, SHB) have not yet been explored.

RESULTS

In the current study, we identified seven CYP genes (six Halloween genes and one AtCYP18A1 gene) related to 20E metabolism. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) showed that AtCYP307A1 and AtCYP307B1 were primarily expressed in the embryonic stage and in the cephalothorax of larvae. RNA interference (RNAi) screening revealed that suppression of AtCYP307A1 expression caused a lethal phenotype during the larval-pupal metamorphosis. Furthermore, Hematoxylin and Eosin staining of the integument showed that the RNAi of AtCYP307A1 inhibited the apolysis and degradation of the old cuticle. In addition, silencing of AtCYP307A1 resulted in significant down-regulation of 20E titers and the expression levels of 20E signaling pathway genes. Finally, the AtCYP307A1 RNAi phenotype was rescued by topical application of 20E.

CONCLUSION

Our studies suggest that AtCYP307A1 involved in 20E synthesis is indispensable during the larval-pupal metamorphosis of beetles, which could serve as a putative insecticide target for pest control. © 2022 Society of Chemical Industry.

Differential expression profiling of Oxycarenus laetus Kirby (Hemiptera: Lygaeidae) upon exposure to gossypol

BACKGROUND

Gossypium hirsutum seeds are rich in gossypol. In addition to its diverse beneficial properties, it is a known anti-fertility inducing agent in humans. Oxycarenus laetus feeds on the cottonseeds and yet its courtship, mating and reproduction is unaffected.

METHODS AND RESULTS

In this study, we performed a transcriptome profiling of O. laetus fed on Abutilon indicum (AB-no gossypol), G. hirsutum (GH-natural gossypol) and 1400 ppm commercial gossypol-soaked GH seeds (GHGO). Illumina NextSeq-500 paired-end 75 bp reads were generated and de novo assembled (48,214 genes) to identify the differentially expressed transcripts (DET) between the samples. Gene enrichment, KEGG pathway and cluster profiling of the DETs resulted in the identification of vital genes involved in the detoxification, pheromone biosynthesis, cuticle protein in the GHGO sample. Cyp4C1, Cyp6a13, Cyp6a14, Cyp4g15, Cyp4em8, Cyp303a1 were the detoxification related genes identified. Similarly, SDR dehydrogenase family 11 and fatty acid synthase in pheromone biosynthesis and cuticle proteins (RR1 and RR2) coding transcripts were found to be differentially expressed.

CONCLUSION

This is the first study to report the expression of genes induced by gossypol in O. laetus. Based on the findings from the DET analysis, we conclude that the detoxification related genes of gossypol treated samples were affected.

Involvement of CYP2 and mitochondrial clan P450s of Helicoverpa armigera in xenobiotic metabolism

Insect CYP2 and mitochondrial clan P450s are relatively conserved genes encoding enzymes generally thought to be involved in biosynthesis or metabolism of endobiotics. However, emerging evidence argues they have potential roles in chemical defense as well, but their actual detoxification functions remain largely unknown. Here, we focused on the full complement of 8 CYP2 and 10 mitochondrial P450s in the generalist herbivore, Helicoverpa armigera. Their varied spatiotemporal expression profiles were analyzed and reflected their specific functions. For functional study of the mitochondrial clan P450s, the redox partners, adrenodoxin reductase (AdR) and adrenodoxin (Adx), were identified from genomes of eight insects and an efficient in vitro electron transfer system of mitochondrial P450 was established by co-expression with Adx and AdR of H. armigera. All CYP2 clan P450s and 8 mitochondrial P450s were successfully expressed in Sf9 cells and compared functionally. In vitro metabolism assays showed that two CYP2 clan P450s (CYP305B1 and CYP18A1) and CYP333B3 (mito clan) could epoxidize aldrin to dieldrin, while CYP305B1 and CYP339A1 (mito clan) have limited but significant hydroxylation capacities to esfenvalerate. CYP303A1 of the CYP2 clan exhibits high metabolic efficiency to 2-tridecanone. Screening the xenobiotic metabolism competence of CYP2 and mitochondrial clan P450s not only provides new insights on insect chemical defense but also can give indications on their physiological functions in H. armigera and other insects.

Gut transcriptome of two bark beetle species stimulated with the same kairomones reveals molecular differences in detoxification pathways

Dendroctonus bark beetles are the most destructive agents in coniferous forests. These beetles come into contact with the toxic compounds of their host's chemical defenses throughout their life cycle, some of which are also used by the insects as kairomones to select their host trees during the colonization process. However, little is known about the molecular mechanisms by which the insects counteract the toxicity of these compounds. Here, two sibling species of bark beetles, D. valens and D. rhizophagus, were stimulated with vapors of a blend of their main kairomones (α-pinene, β-pinene and 3-carene), in order to compare the transcriptional response of their gut. A total of 48 180 unigenes were identified in D. valens and 43 704 in D. rhizophagus, in response to kairomones blend. The analysis of differential gene expression showed a transcriptional response in D. valens (739 unigenes, 0.58–10.36 Log2FC) related to digestive process and in D. rhizophagus (322 unigenes 0.87–13.08 Log2FC) related to xenobiotics metabolism. The expression profiles of detoxification genes mainly evidenced the up-regulation of COEs and GSTs in D. valens, and the up-regulation of P450s in D. rhizophagus. Results suggest that terpenes metabolism comes accompanied by an integral hormetic response, result of compensatory mechanisms, including the activation of other metabolic pathways, to ensure the supply of energy and the survival of organisms which is specific for each species, according to its life history and ecological strategy.

Detoxification gene families in Phylloxera: Endogenous functions and roles in response to the environment

Phylloxera, Daktulosphaira vitifoliae, is an agronomic pest that feeds monophagously on grapevine, Vitis spp. host plants. Phylloxera manipulates primary and secondary plant metabolism to establish either leaf or root galls. We manually annotated 198 detoxification genes potentially involved in plant host manipulation, including cytochrome P450 (66 CYPs), carboxylesterase (20 CCEs), glutathione-S-transferase (10 GSTs), uridine diphosphate-glycosyltransferase (35 UGTs) and ABC transporter (67 ABCs) families. Transcriptomic expression patterns of these detoxification genes were analyzed for root and leaf galls. In addition to these transcriptomic analyses, we reanalyzed recent data from L1 and L2-3 stages feeding on tolerant and resistant rootstock. Data from two agricultural pest aphids, the generalist Myzus persicae and the Fabaceae specialist Acyrthosiphon pisum, and from the true bug vector of Chagas disease, Rhodnius prolixus, were used to perform phylogenetic analyses for each detoxification gene family. We found expansions of several gene sub-families in the genome of D. vitifoliae. Phylogenetically close genes were found to be organized in clusters in the same genomic position and orientation suggesting recent successive duplications. These results highlight the roles of the phylloxera detoxification gene repertoire in insect physiology and in adaptation to plant secondary metabolites, and provide gene candidates for further functional analyses.

The microRNA miR-184 regulates the CYP303A1 transcript level to control molting of Locusta migratoria

Cytochrome P450 monooxygenases (CYPs) play essential physiological functions in insects. CYP303A1 is highly conserved in insect species studied to date, and shows an indispensable role for adult eclosion in both Locusta migratoria and Drosophila melanogaster. However, how CYP303A1 is regulated to control insect developmental processes remains uninvestigated. In this study, we discovered functional binding sites for miR-184 in the coding sequence of LmCYP303A1. The luciferase reporter assay showed that miR-184 could target LmCYP303A1 and regulate its expression in vitro. Furthermore, overexpression of miR-184 through microinjection of agomir to locusts reduced the transcripts of LmCYP303A1 and led to abnormal molting, which is similar to the phenotype of silencing LmCYP303A1 by direct injection of dsLmCYP303A1 to locusts. Meanwhile, down-regulation of miR-184 by injection of antagomir increased the LmCYP303A1 transcript and caused molting defects. These findings suggested that miR-184 could target LmCYP303A1 to regulate the molting process in L. migratoria, which might be considered as a novel target for pest control.

Matrix metalloproteinases are involved in eclosion and wing expansion in the American cockroach, Periplaneta americana

Matrix metalloproteinases (MMPs) are the major proteinases that process or degrade numerous extracellular matrix (ECM) components and are evolutionarily conserved from nematodes to humans. During molting in insects, the old cuticle is removed and replaced by a new counterpart. Although the regulatory mechanisms of hormones and nutrients in molting have been well studied, very little is known about the roles of ECM-modifying enzymes in this process. Here, we found that MMPs are necessary for imaginal molting of the American cockroach, Periplaneta americana. Inhibition of Mmp activity via inhibitor treatment led to the failure of eclosion and wing expansion. Five Mmps genes were identified from the P. americana genome, and PaMmp2 played the dominant roles during molting. Further microscopic investigations showed that newly formed adult cuticles were attenuated and that then chitin content was reduced upon Mmp inhibition. Transcriptomic analysis of the integument demonstrated that multiple signaling and metabolic pathways were changed. Microscopic investigation of the wings showed that epithelial cells were restrained together because they were incapable of degrading the ECM upon Mmp inhibition. Transcriptomic analysis of the wing identified dozens of possible genes functioned in wing expansion. This is the first study to show the essential roles of Mmps in the nymph-adult transition of hemimetabolous insects.

Gene expression during larval caste determination and differentiation in intermediately eusocial bumblebees, and a comparative analysis with advanced eusocial honeybees

The queen‐worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally‐induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste‐associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen‐dependent period during which their caste fate as adults is determined followed by a nutrition‐sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA‐seq and qRT‐PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen‐dependent period. Relatively few novel (i.e., taxonomically‐restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late‐instar larvae in the worker pathway. We compared sets of caste‐associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste‐associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

Both LmCYP4G genes function in decreasing cuticular penetration of insecticides in Locusta migratoria

BACKGROUND

Cuticular hydrocarbons (CHCs) have a critical role in preventing desiccation and penetration of xenobiotics in insects. Previous studies have shown that cytochrome P450 subfamily 4G (CYP4G) enzymes are oxidative decarbonylases, essential for CHC biosynthesis. However, it is unclear whether there are functional differences between the two CYP4G genes in most insects. In Locusta migratoria, we identified two CYP4G genes (LmCYP4G62 and LmCYP4G102). LmCYP4G102 plays a critical role in the synthesis of CHCs, but the function of LmCYP4G62 is unknown.

RESULTS

We identified, characterized, and compared two LmCYP4G genes, based on L. migratoria transcriptomic and genomic databases. RT-qPCR showed that both were highly expressed in tissues with which oenocytes are associated, the integument and fat body. Immunostaining indicated that LmCYP4G62 and LmCYP4G102 were highly abundant in oenocytes in these tissues. However, the two enzymes had a different subcellular distribution, with LmCYP4G62 localized on the plasma membrane and LmCYP4G102 dispersed throughout the oenocyte cytoplasm, presumably on the endoplasmic reticulum. RNA interference-mediated gene silencing against each of the two genes resulted in reduced CHC contents, in all classes for LmCYP4G102, but mostly shorter chain CHCs for LmCYP4G62. Silencing of both genes resulted in increased insecticide penetration through the cuticle, and increased locust susceptibility to desiccation and insecticides.

CONCLUSION

Our studies suggest that both LmCYP4G62 and LmCYP4G102 contribute to hydrocarbon biosynthesis and play key roles in protecting locusts from water loss and insecticide penetration, but they are not fully redundant. Further, the two LmCYP4G genes might be used as new targets for insect pest management. © 2020 Society of Chemical Industry.

Knockdown of LmCYP303A1 alters cuticular hydrocarbon profiles and increases the susceptibility to desiccation and insecticides in Locusta migratoria

Cytochrome P450 monooxygenases (CYPs) serve many functions in insects, from the regulation of development to xenobiotic detoxification. Several conserved CYPs have been shown to play a role in insect growth and development. CYP303A1 is a highly conserved CYP with a single ortholog in most insects, but its underlying molecular characteristics and specific physiological functions remain poorly understood. In Drosophila melanogaster and Locusta migratoria, CYP303A1 is indispensable for eclosion to adult. Here, we report additional functions of the locust gene LmCYP303A1 in nymphal molts, cuticular lipid deposition and insecticide penetration. RT-qPCR revealed that LmCYP303A1 had a high expression level before ecdysis and was highly expressed in integument, wing pads, foregut and hindgut. Suppression of LmCYP303A1 expression by RNA interference (RNAi) caused a lethal phenotype with molting defect from nymph to nymph. In addition, LmCYP303A1 RNAi resulted in locusts being more susceptible to desiccation and to insecticide toxicity. Furthermore, knockdown of LmCYP303A1 efficiently suppressed the transcript level of key genes (ELO7, FAR15 and CYP4G102) responsible for cuticular hydrocarbon (CHC) synthesis, which led to a decrease in some CHC levels. Taken together, our results suggest that one of the functions of LmCYP303A1 is to regulate the biosynthesis of CHC, which plays critical roles in protecting locusts from water loss and insecticide penetration.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum

In the postgenomics era, comparative genomics have advanced the understanding of evolutionary processes of neuropeptidergic signaling systems. The evolutionary origin of many neuropeptidergic signaling systems can be traced date back to early metazoan evolution based on the conserved sequences. Insect parathyroid hormone receptor (iPTHR) was previously described as an ortholog of vertebrate PTHR that has a well-known function in controlling bone remodeling. However, there was no sequence homologous to PTH sequence in insect genomes, leaving the iPTHR as an orphan receptor. Here, we identified the authentic ligand insect PTH (iPTH) for the iPTHR. The taxonomic distribution of iPTHR, which is lacking in Diptera and Lepidoptera, provided a lead for identifying the authentic ligand. We found that a previously described orphan ligand known as PXXXamide (where X is any amino acid) described in the cuttlefish Sepia officinalis has a similar taxonomic distribution pattern as iPTHR. Tests of this peptide, iPTH, in functional reporter assays confirmed the interaction of the ligand-receptor pair. Study of a model beetle, Tribolium castaneum, was used to investigate the function of the iPTH signaling system by RNA interference followed by RNA sequencing and phenotyping. The results suggested that the iPTH system is likely involved in the regulation of cuticle formation that culminates with a phenotype of defects in wing exoskeleton maturation at the time of adult eclosion. Moreover, RNAi of iPTHRs also led to significant reductions in egg numbers and hatching rates after parental RNAi.

Vertebrate parathyroid hormone (PTH) and its receptors have been extensively studied with respect to their function in bone remodeling and calcium metabolism. Insect parathyroid hormone receptors (iPTHRs) have been previously described as counterparts of vertebrate PTHRs, however, they are still orphan receptors for which the authentic ligands and biological functions remain unknown. We describe an insect form of parathyroid hormone (iPTH) by analyzing its interactions with iPTHRs. Identification of this new insect peptidergic system proved that the PTH system is an ancestral signaling system dating back to the evolutionary time before the divergence of protostomes and deuterostomes. We also investigated the functions of the iPTH system in a model beetle Tribolium castaneum by using RNA interference. RNA interference of iPTHR resulted in defects in wing exoskeleton maturation and fecundity. Based on the differential gene expression patterns and the phenotype induced by RNAi, we propose that the iPTH system is likely involved in the regulation of exoskeletal cuticle formation and fecundity in insects.

Aryl hydrocarbon receptor regulates the expression of LmGSTd7 and is associated with chlorpyrifos susceptibility in Locusta migratoria

BACKGROUND

The aryl hydrocarbon receptor (AhR) belongs to the bHLH-PAS (basic Helix-Loop-Helix - Period/ARNT/Single minded) family of transcription factors. AhR is a ligand-activated transcription factor, which participates in the sensing and transmitting stimuli of endogenous and exogenous chemicals, and subsequently activates the transcription of genes related to various physiological and detoxification functions.

RESULT

In this study, a single full-length LmAhR sequence was cloned and characterized. RNA interference (RNAi) and insecticide bioassays showed that LmAhR plays a vital role in chlorpyrifos susceptibility. To better identify aryl hydrocarbon receptor from locusta migratoria (LmAhR)-regulated genes involved in chlorpyrifos susceptibility, a comparative transcriptome analysis was performed using double-stranded (ds)GFP- and dsLmAhR-injected Locusta migratoria. Differential gene expression analysis identified 145 down-regulated and 67 up-regulated genes (P ≤ 0.05 and fold change ≥2) in dsLmAhR-knockdown insects. We selected 27 down-regulated genes and verified their expression levels using reverse transcription quantitative PCR. Finally, one glutathione S-transferase (GST) gene (LmGSTd7) was selected as a candidate detoxification gene and was further validated via RNAi and chlorpyrifos bioassays.

CONCLUSION

Our data suggest that AhR is associated with chlorpyrifos susceptibility via the regulation of LmGSTd7 expression in L. migratoria. © 2019 Society of Chemical Industry.