The silkwormBombyx moriis an important economic insect for producing silk, the “queen of fabrics”. The currently available genomes limit the understanding of its genetic diversity and the discovery of valuable alleles for breeding. Here, we deeply re-sequence 1,078 silkworms and assemble long-read genomes for 545 representatives. We construct a high-resolution pan-genome dataset representing almost the entire genomic content in the silkworm. We find that the silkworm population harbors a high density of genomic variants and identify 7308 new genes, 4260 (22%) core genes, and 3,432,266 non-redundant structure variations (SVs). We reveal hundreds of genes and SVs that may contribute to the artificial selection (domestication and breeding) of silkworm. Further, we focus on four genes responsible, respectively, for two economic (silk yield and silk fineness) and two ecologically adaptive traits (egg diapause and aposematic coloration). Taken together, our population-scale genomic resources will promote functional genomics studies and breeding improvement for silkworm.

Insect cuticle is critical for the environmental adaptability and insecticide resistance of insects. However, there is no clear understanding of the structure and protein components of the cuticle during each developmental stage of holometabolous insects, and knowledge about the protein components within each layer is vague. We conducted serial sectioning, cuticular structure analysis, and transcriptome sequencing of the larval, pupal, and adult cuticles of Bombyx mori. The deposition processes of epicuticle, exocuticle, and endocuticle during larval, pupal, and adult cuticle formation were similar. Transcriptome analysis showed that these cuticle formations share 74% of the expressed cuticular protein (CP) genes and 20 other structural protein genes, such as larval serum protein and prisilkin. There are seven, six, and eleven stage-specific expressed CP genes in larval, pupal, and adult cuticles, respectively. The types and levels of CP genes may be the key determinants of the properties of each cuticular layer. For example, the CPs of the RR-2 protein family with high contents of histidine (His) are more essential for the exocuticle. Functional analysis suggested that BmorCPAP1-H is involved in cuticle formation. This study not only offers an in-depth understanding of cuticle morphology and protein components but also facilitates the elucidation of molecular mechanisms underlying cuticle formation in future studies.

Flight loss has occurred in many winged insect taxa. The flightless silkmoth Bombyx mori, is domesticated from the wild silkmoth, Bombyx mandarina, which can fly. In this paper, we studied morphological characteristics attributed to flightlessness in silkmoths. Three domestic flightless B. mori strains and one B. mandarina population were used to compare morphological components of the flight apparatus, including wing characteristics (shape, forewing area, loading, and stiffness), flight muscle (weight, ratio, and microscopic detail) and body mass. Compared with B. mandarina, B. mori strains have a larger body, greater wing loading, more flexible wings and a lower flight muscle ratio. The arrangement in microscopy of dorsal longitudinal flight muscles (DLFMs) of B. mori was irregular. Comparative analysis of the sexes suggests that degeneration of flight muscles and reduction of wing mechanical properties (stiffness) are associated with silkmoth flightlessness. The findings provide important clues for further research of the molecular mechanisms of B. mori flight loss.

The insect glycoside hydrolase family 20 β-N-acetylhexosaminidases (HEXs) are key enzymes involved in chitin degradation. In this study, nine HEX genes in Bombyx mori were identified by genome-wide analysis. Bioinformatic analysis based on the transcriptome database indicated that each gene had a distinct expression pattern. qRT-PCR was performed to detect the expression pattern of the chitooligosaccharidolytic β-N-acetylglucosaminidase (BmChiNAG). BmChiNAG was highly expressed in chitin-rich tissues, such as the epidermis. In the wing disc and epidermis, BmChiNAG has the highest expression level during the wandering stage. CRISPR/Cas9-mediated BmChiNAG deletion was used to study the function. In the BmChiNAG-knockout line, 39.2% of female heterozygotes had small and curly wings. The ultrastructure of a cross-section showed that the lack of BmChiNAG affected the stratification of the wing membrane and the formation of the correct wing vein structure. The molting process of the homozygotes was severely hindered during the larva to pupa transition. Epidermal sections showed that the endocuticle of the pupa was not degraded in the mutant. These results indicate that BmChiNAG is involved in chitin catabolism and plays an important role in the molting and wing development of the silkworm, which highlights the potential of BmChiNAG as a pest control target.

Mutation rate is a crucial parameter in evolutionary genetics. However, the mutation rate of most species as well as the extent to which the environment can alter the genome of multicellular organisms remain poorly understood. Here, we used parents–progeny sequencing to investigate the mutation rate and spectrum of the domestic silkworm (Bombyx mori) among normal and two temperature stress conditions (32 °C and 0 °C). The rate of single-nucleotide mutations in the normal temperature rearing condition was 0.41 × 10⁻⁸ (95% confidence interval, 0.33 × 10⁻⁸–0.49 × 10⁻⁸) per site per generation, which was up to 1.5-fold higher than in four previously studied insects. Moreover, the mutation rates of the silkworm under the stresses are significantly higher than in normal conditions. Furthermore, the mutation rate varies less in gene regions under normal and temperature stresses. Together, these findings expand the known diversity of the mutation rate among eukaryotes but also have implications for evolutionary analysis that assumes a constant mutation rate among species and environments.

The spectrin family genes play critical roles in cytoskeletal organization and cellular integrity, yet their evolutionary and functional significance in non-classical model organisms remains poorly explored. Here, we systematically identified and characterized spectrin family genes in the silkworm Bombyx mori. Genome-wide analysis identified 17 predicted spectrin genes, which were integrated into eight optimized members through transcriptome-guided structural refinement. Multi-species genomic analysis revealed 8, 23, and 24 spectrin family genes in Drosophila melanogaster, Mus musculus, and Homo sapiens, respectively. Phylogenetic analysis revealed conserved clades across insects and mammals, with gene family expansions in vertebrates. Spatiotemporal expression profiling demonstrated ubiquitous expression of these genes during silkworm development. Population genomic analyses detected strong selection signatures in BmTrio during domestication and implicated BmBeta_spc as a candidate gene for silk yield enhancement in Chinese-improved strains (CHN-I). Expression profiles of parental strains and F₁ offspring from a commercial hybrid cross (Jingsong × Haoyue) revealed BmBeta_spc expression correlating with heterosis in silk yield traits. This study elucidates the characterization and functional relevance of silkworm spectrin genes, providing insights into their roles in domestication and breeding.