A chromosome‐scale genome assembly of Antheraea pernyi (Saturniidae, Lepidoptera)

Day-night gene expression reveals circadian gene disco as a candidate for diel-niche evolution in moths

Temporal ecological niche partitioning is an underappreciated driver of speciation. While insects have long been models for circadian biology, the genes and circuits that allow adaptive changes in diel-niches remain poorly understood. We compared gene expression in closely related day- and night-active non-model wild silk moths, with otherwise similar ecologies. Using an ortholog-based pipeline to compare RNA-Seq patterns across two moth species, we find over 25 pairs of gene orthologs showing differential expression. Notably, the gene disco, involved in circadian control, optic lobe and clock neuron development in Drosophila, shows robust adult circadian mRNA cycling in moth heads. Disco is highly conserved in moths and has additional zinc-finger domains with specific nocturnal and diurnal mutations. We propose disco as a candidate gene for the diversification of temporal diel-niche in moths.

Inundative practice for screening siRNA management candidates against a notorious predatory beetle using olfactory silencing

Calosoma maximoviczi, a predatory pest beetle, poses a significant threat to wild silk farm production due to its predation on wild silkworms. Given the coexistence of this species with beneficial silkworms in the farm orchards, chemical pesticides are not an ideal solution for controlling its population. In this study, we employed a comprehensive multi-target RNA interference (RNAi) approach to disrupt the olfactory perception of C. maximoviczi through independently silencing 16 odorant receptors (ORs) in the respective genders. Specifically, gene-specific siRNAs were designed to target a panel of ORs, allowing us to investigate the specific interactions between odorant receptors and ligands within this species. Our investigation led to identifying four candidate siOR groups that effectively disrupted the beetle's olfactory tracking of various odorant ligands associated with different trophic levels. Furthermore, we observed sex-specific differences in innate RNAi responses reflected by subsequent gene expression, physiological and behavioral consequences, underscoring the complexity of olfactory signaling and emphasizing the significance of considering species/sex-specific traits when implementing pest control measures. These findings advance our understanding of olfactory coding patterns in C. maximoviczi beetles and establish a foundation for future research in the field of pest management strategies.

Exploring the proteins and metabolites associated with male antennae responses to female exposure of Antheraea pernyi (Lepidoptera: Saturniidae) moths

Detection of sex pheromones of insects relies on the antennae. The female pheromone signal transmission in the male antennae ultimately initiates the courtship and mating behaviors of males. To investigate the proteins and metabolites involved in this neural transduction, integrative proteomics and metabolomics analysis including tandem mass tag (TMT) proteomic quantification and liquid chromatography tandem mass spectrometry (LC/MS)-based metabolomics was adopted for comparing proteomic and metabolic changes between the antennae of male moths following stimulation by females and the non-stimulated males of Antheraea pernyi (Guérin-Méneville, Lepidoptera: Saturniidae) in this study. A total of 92 differentially expressed proteins (DEPs) containing 52 upregulated and 40 downregulated proteins and 545 differentially expressed metabolites (DEMs) including 218 upregulated and 327 downregulated metabolites were identified from the antennae of female-stimulated male moths based on the proteome and metabolome data, respectively. Bioinformatics analysis was performed for the 45 DEPs and 160 DEMs, including Gene Ontology (GO), Clusters of Orthologous Groups (COG), and Kyoto Encylopaedia of Genes and Genomes (KEGG) enrichment analysis and Human Metabolome Database (HMDB) annotation. A number of DEPs and DEMs related to neural transmission of female pheromone signals in the male antennae of A. pernyi were screened, including tyrosine hydroxylase, cryptochrome-1, tachykinin, arylalkylamine N-acetyltransferase, cadherin-23, glutathione S-transferase delta 3, tyramine, tryptamine, n-oleoyl dopamine, n-stearoyl dopamine, and n-stearoyl tyrosine. The altered expression levels of those proteins or metabolites were speculated involved in regulating the neuron activity for enhanced transmission of neural impulses and continuous perception, reception, and transduction of female pheromone signals. Our findings yielded novel insights into the potential mechanisms in the antennae of male A. pernyi responding to female attraction.

Transcriptome sequencing of the central nervous system to identify the neuropeptides and neuropeptide receptors of Antheraea pernyi

Neuropeptides and neuropeptide receptors are crucial regulators for the behavior, lifecycle, and physiology of insects and are mainly produced and released from the neurosecretory cells of the central nervous system (CNS). In this study, RNA-seq was employed to investigate the transcriptome profile of the CNS which is composed of the brain and ventral nerve cord (VNC) of Antheraea pernyi. From the data sets, a total of 18 and 42 genes were identified, which respectively encode the neuropeptides and neuropeptide receptors involved in regulating multiple behaviors including feeding, reproductive behavior, circadian locomotor, sleep, and stress response and physiological processes such as nutrient absorption, immunity, ecdysis, diapause, and excretion. Comparison of the patterns of expression of those genes between the brain and VNC showed that most had higher levels of expression in the brain than VNC. Besides, 2760 differently expressed genes (DEGs) (1362 up-regulated and 1398 down-regulated ones between the B and VNC group) were also screened and further analyzed via gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses. The results of this study could provide comprehensive profiles of the neuropeptides and neuropeptide receptors of A. pernyi CNS and lay the foundation for further research into their functions.

Genome-Wide Identification of Detoxification Genes in Wild Silkworm Antheraea pernyi and Transcriptional Response to Coumaphos

For a half-century, the commercial wild silkworm, Antheraea pernyi, has been protected by coumaphos, which is an internal organophosphorus insecticide used to kill the potential parasitic fly larvae inside. Knowledge about the detoxification genes of A. pernyi as well as the detoxification mechanism for this species remains severely limited. In this study, we identified 281 detoxification genes (32 GSTs, 48 ABCs, 104 CYPs, and 97 COEs) in the genome of this insect, which are unevenly distributed over 46 chromosomes. When compared to the domesticated silkworm, Bombyx mori, a lepidopteran model species, A. pernyi has a similar number of ABCs, but a greater number of GSTs, CYPs, and COEs. By transcriptome-based expression analysis, we found that coumaphos at a safe concentration level significantly changed the pathways related to ATPase complex function and the transporter complex in A. pernyi. KEGG functional enrichment analysis indicated that protein processing in the endoplasmic reticulum was the most affected pathway after coumaphos treatment. Finally, we identified four significantly up-regulated detoxification genes (ABCB1, ABCB3, ABCG11, and ae43) and one significantly down-regulated detoxification gene (CYP6AE9) in response to coumaphos treatment, suggesting that these five genes may contribute to detoxification of coumaphos in A. pernyi. Our study provides the first set of detoxification genes for wild silkworms from Saturniidae and highlights the importance of detoxification gene repertoire in insect pesticide tolerance.

Genome-Wide Identification and Transcriptome-Based Expression Profile of Cuticular Protein Genes in Antheraea pernyi

Antheraea pernyi is one of the most famous edible and silk-producing wild silkworms of Saturniidae. Structural cuticular proteins (CPs) are the primary component of insect cuticle. In this paper, the CPs in the genome of A. pernyi were identified and compared with those of the lepidopteran model species Bombyx mori, and expression patterns were analyzed based on the transcriptomic data from the larval epidermis/integument (epidermis in the following) and some non-epidermis tissues/organs of two silkworm species. A total of 217 CPs was identified in the A. pernyi genome, a comparable number to B. mori (236 CPs), with CPLCP and CPG families being the main contribution to the number difference between two silkworm species. We found more RR-2 genes expressed in the larval epidermis of fifth instar of A. pernyi than B. mori, but less RR-2 genes expressed in the prothoracic gland of A. pernyi than B. mori, which suggests that the hardness difference in the larval epidermis and prothoracic gland between the two species may be caused by the number of RR-2 genes expressed. We also revealed that, in B. mori, the number of CP genes expressed in the corpus allatum and prothoracic gland of fifth instar was higher than that in the larval epidermis. Our work provided an overall framework for functional research into the CP genes of Saturniidae.

Genome-Wide Identification of M14 Family Metal Carboxypeptidases in Antheraea pernyi (Lepidoptera: Saturniidae)

The M14 family metal carboxypeptidase genes play an important role in digestion and pathogenic infections in the gut of insects. However, the roles of these genes in Antheraea pernyi (Guérin-Méneville, 1855) remain to be analyzed. In the present study, we cloned a highly expressed M14 metal carboxypeptidase gene (ApMCP1) found in the gut and discovered that it contained a 1,194 bp open reading frame encoding a 397-amino acid protein with a predicted molecular weight of 45 kDa. Furthermore, 14 members of the M14 family metal carboxypeptidases (ApMCP1-ApMCP14) were identified in the A. pernyi genome, with typical Zn_pept domains and two Zn-anchoring motifs, and were further classified into M14A, M14B, and M14D subfamilies. Expression analysis indicated that ApMCP1 and ApMCP9 were mainly expressed in the gut. Additionally, we observed that ApMCP1 and ApMCP9 displayed opposite expression patterns after starvation, highlighting their functional divergence during digestion. Following natural infection with baculovirus NPV, their expression was significantly upregulated in the gut of A. pernyi. Our results suggest that the M14 family metal carboxypeptidase genes are conservatively digestive enzymes and evolutionarily involved in exogenous pathogenic infections.

Pupal Diapause Termination and Transcriptional Response of Antheraea pernyi (Lepidoptera: Saturniidae) Triggered by 20-Hydroxyecdysone

The pupal diapause of univoltine Antheraea pernyi hampers sericultural and biotechnological applications, which requires a high eclosion incidence after artificial diapause termination to ensure production of enough eggs. The effect of pupal diapause termination using 20-hydroxyecdysone (20E) on the eclosion incidence has not been well-documented in A. pernyi. Here, the dosage of injected 20E was optimized to efficiently terminate pupal diapause of A. pernyi, showing that inappropriate dosage of 20E can cause pupal lethality and a low eclosion incidence. The optimal ratio of 20E to 1-month-old pupae was determined as 6 μg/g. Morphological changes showed visible tissue dissociation at 3 days post-injection (dpi) and eye pigmentation at 5 dpi. Comprehensive transcriptome analysis identified 1,355/1,592, 494/203, 584/297, and 1,238/1,404 upregulated and downregulated genes at 1, 3, 6, and 9 dpi, respectively. The 117 genes enriched in the information processing pathways of “signal transduction” and “signaling molecules and interaction” were upregulated at 1 and 3 dpi, including the genes involved in FOXO signaling pathway. One chitinase, three trehalase, and five cathepsin genes related to energy metabolism and tissue dissociation showed high expression levels at the early stage, which were different from the upregulated expression of four other chitinase genes at the later stage. Simultaneously, the expression of several genes involved in molting hormone biosynthesis was also activated between 1 and 3 dpi. qRT-PCR further verified the expression patterns of two ecdysone receptor genes (EcRB1 and USP) and four downstream response genes (E93, Br-C, βFTZ-F1, and cathepsin L) at the pupal and pharate stages, respectively. Taken together, these genes serve as a resource for unraveling the mechanism underlying pupal-adult transition; these findings facilitate rearing of larvae more than once a year and biotechnological development through efficient termination of pupal diapause in A. pernyi in approximately half a month.

A chromosome-level genome assembly of Paracymoriza distinctalis (Lepidoptera: Crambidae: Acentropinae)

Paracymoriza distinctalis is a semiaquatic lepidopteran insect, which is of great value for studying the differentiation of the Pyraloidea superfamily. However, the understanding of heredity, evolution, and functional genomics of P. distinctalis are limited by few genome-wide resources. Here, we applied PacBio sequencing and the chromosome capture technique to assemble the first P. distinctalis genome from a single female individual. The genome size is 1.2 Gb with 32 chromosomes and the N50 is 38.91 Mb. Approximately 576.37 Mb, accounting for 48.93% of the genome, was identified as repeats. The genome comprises 39,003 protein-coding genes, 66.56% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of P. distinctalis and Chilo suppressalis lived ~83.5 million years ago. This chromosome-level genome assembly work is not only conducive to the understanding of P. distinctalis, but also may promote the study of the genomes of other lepidopteran species.

Behavioural regulator and molecular reception of a double-edge-sword hunter beetle

BACKGROUND

The black carabid beetle Calosoma maximoviczi is a successful predator that serves as both a beneficial insect and a severe threat to economic herbivores. Its hunting technique relies heavily on olfaction, but the underlying mechanism has not been studied. Here, we report the electrophysiological, ecological and molecular traits of bioactive components identified from a comprehensive panel of natural odorants in the beetle-prey-plant system. The aim of this work was to investigate olfactory perceptions and their influence on the behaviours of C. maximoviczi.

RESULTS

Among the 200 identified volatiles, 18 were concentrated in beetle and prey samples, and 14 were concentrated in plants. Insect feeding damage to plants led to a shift in the emission fingerprint. Twelve volatiles were selected using successive electrophysiological tests. Field trials showed that significant sex differences existed when trapping with a single chemical or chemical mixture. Expression profiles indicated that sex-biased catches were related to the expression of 15 annotated CmaxOBPs and 40 CmaxORs across 12 chemosensory organs. In silico evaluations were conducted with 16 CmaxORs using modelling and docking. Better recognition was predicted for the pairs CmaxOR5-(Z)-3-hexenyl acetate, CmaxOR6-β-caryophyllene, CmaxOR18-(E)-β-ocimene and CmaxOR18-tetradecane, with higher binding affinity and a suitable binding pocket. Lastly, 168Y in CmaxOR6 and 142Y in CmaxOR18 were predicted as key amino acid residues for binding β-caryophyllene and tetradecane, respectively.

CONCLUSION

This work provides an example pipeline for de novo investigation in C. maximoviczi baits and the underlying olfactory perceptions. The results will benefit the future development of trapping-based integrated pest management strategies and the deorphanization of odorant receptors in ground beetles. © 2022 Society of Chemical Industry.

A Gut-Specific LITAF-Like Gene in Antheraea pernyi (Lepidoptera: Saturniidae) Involved in the Immune Response to Three Pathogens

Antheraea pernyi (Guérin-Méneville 1855) is an important resource for silk, food, and biohealth products; however, exogenous pathogens largely affect the commercial application potential of this species. Since the gut is a key organ for the digestion and absorption of nutrients as well as for immune defense, we used comparative transcriptome analysis to screen for a gut-specific molecular tool for further functional research in A. pernyi. In total, 3,331 differentially expressed genes (DEGs) were identified in the gut compared with all other pooled tissues of A. pernyi, including 1,463 upregulated genes in the gut. Among these, we further focused on a lipopolysaccharide-induced tumor necrosis factor-α factor (LITAF) gene because of its high gut-specific expression and the presence of a highly conserved SIMPLE-like domain, which is related to the immune response to pathogenic infections in many species. The cDNA sequence of ApLITAF was 447-bp long and contained a 243-bp open reading frame encoding an 80-amino acid protein. Immune challenge assays indicated that ApLITAF expression was significantly upregulated in the gut of A. pernyi naturally infected with nucleopolyhedrovirus (NPV) or fed leaves infected with the gram-negative bacterium Escherichia coli (Migula 1895) and the gram-positive bacterium Bacillus subtilis (Ehrenberg 1835). Cell transfection showed that ApLITAF localized to the lysosome. Collectively, these results suggested that ApLITAF played a role in the immune response of A. pernyi and could facilitate the future research and breeding application in this species.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.