Feeding preference is critical for insect adaptation and survival. However, little is known regarding the determination of insect feeding preference, and the genetic basis is poorly understood. As a model lepidopteran insect with economic importance, the domesticated silkworm, Bombyx mori, is a well-known monophagous insect that predominantly feeds on fresh mulberry leaves. This species-specific feeding preference provides an excellent model for investigation of host-plant selection of insects, although the molecular mechanism underlying this phenomenon remains unknown. Here, we describe the gene GR66, which encodes a putative bitter gustatory receptor (GR) that is responsible for the mulberry-specific feeding preference of B. mori. With the aid of a transposon-based, clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 nuclease (Cas9) system, the GR66 locus was genetically mutated, and homozygous mutant silkworm strains with truncated gustatory receptor 66 (GR66) proteins were established. GR66 mutant larvae acquired new feeding activity, exhibiting the ability to feed on a number of plant species in addition to mulberry leaves, including fresh fruits and grain seeds that are not normally consumed by wild-type (WT) silkworms. Furthermore, a feeding choice assay revealed that the mutant larvae lost their specificity for mulberry. Overall, our findings provide the first genetic and phenotypic evidences that a single bitter GR is a major factor affecting the insect feeding preference. The molecular mechanism underlying species-specific feeding preference in insects is poorly understood. The silkworm, Bombyx mori, is a typical monophagous plant-eating insect, but the genetic basis for its famous mulberry-specific feeding preference is unknown. Here, we identify gustatory receptor 66 (GR66) as a determinant of the silkworm’s mulberry-specific monophagy. GR66-mutant larvae generated by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 nuclease (Cas9) acquired new feeding activity and showed the ability to feed on various plant species that are not normally consumed by the wild-type (WT) animals; a two-choice assay demonstrated that the mutant larvae had lost their feeding preference for mulberry. Our genetic and phenotypic evidence therefore demonstrates that GR66 is a major factor affecting the feeding preference of the silkworm.

Diapause is a common seasonal adaptive strategy that regulates annual timing in insects. Very few causal loci underlying diapause variation have yet been identified. By leveraging cross-mapping and genome-wide association analysis, we identified the N terminus of the clock protein CYCLE as a major causal effector underlying embryonic diapause differences in the silk moth. We found that the nondiapause phenotype in polyvoltine strains results from a specific deletion that disrupts an alternative isoform of CYCLE. We further demonstrated that different CYCLE isoforms contribute to a functional diversity in modulating circadian rhythms and diapause, which has been preserved in Lepidoptera for at least 110 million years. Our study proposes a model that explains how adaptive phenotypes can evolve rapidly without affecting related essential functions.

Cocoon and silk yields are the most important characteristics of sericulture. However, few studies have examined the genes that modulate these features. Further studies of these genes will be useful for improving the products of sericulture. JingSong (JS) and Lan10 (L10) are two strains having significantly different cocoon and silk yields. In the current study, RNA-Seq and quantitative polymerase chain reaction (qPCR) were performed on both strains in order to determine divergence of the silk gland, which controls silk biosynthesis in silkworms. Compared with L10, JS had 1375 differentially expressed genes (DEGs; 738 up-regulated genes and 673 down-regulated genes). Nine enriched gene ontology (GO) terms were identified by GO enrichment analysis based on these DEGs. KEGG enrichment analysis results showed that the DEGs were enriched in three pathways, which were mainly associated with the processing and biosynthesis of proteins. The representative genes in the enrichment pathways and ten significant DEGs were further verified by qPCR, the results of which were consistent with the RNA-Seq data. Our study has revealed differences in silk glands between the two silkworm strains and provides a perspective for understanding the molecular mechanisms determining silk yield.

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the primary pathogens of the silkworm. Cytochrome c (cytc) showed a significant response to BmNPV infection in our previous transcriptome study. However, little is known about the role of Bombyx mori cytc (Bmcytc) in resistance to BmNPV infection. In this study, the expression levels analysis of Bmcytc showed stable expression levels in selected tissues of the resistant strain AN following BmNPV infection, while there was downregulation in the susceptible strain p50, except in the malpighian tubule. To further study the role of Bmcytc in viral infection, Bmcytc was knocked down with siRNA in vitro, resulting in significant downregulation of selected downstream genes of the mitochondrial pathway, including Bmapaf, Bmcaspase-Nc, and Bmcaspase-1; this was also confirmed by overexpression of Bmcytc using the pIZT/V5-His-mCherry insect vector, except Bmcaspase-1. Moreover, knockdown of Bmcytc significantly promoted the infection process of BmNPV in vitro, while the infection was inhibited by overexpression of Bmcytc at the early stage and subsequently increased rapidly. Based on these results, we concluded that Bmcytc plays a vital role in BmNPV infection by regulating the mitochondrial apoptosis pathway. Our work provides valuable data for the clarification of the mechanism of silkworm resistance to BmNPV infection.