Silk is an important natural fiber of high economic value, and thus genetic study of the silkworm is a major area of research. Transcriptome analysis can provide guidance for genetic studies of silk yield traits. In this study, we performed a transcriptome comparison using multiple silkworms with different silk yields. A total of 22 common differentially expressed genes (DEGs) were identified in multiple strains and were mainly involved in metabolic pathways. Among these, seven significant common DEGs were verified by quantitative reverse transcription polymerase chain reaction, and the results coincided with the findings generated by RNA sequencing. Association analysis showed that BGIBMGA003330 and BGIBMGA005780 are significantly associated with cocoon shell weight and encode uridine nucleosidase and small heat shock protein, respectively. Functional annotation of these genes suggest that these play a role in silkworm silk gland development or silk protein synthesis. In addition, we performed principal component analysis (PCA) in combination with wild silkworm analysis, which indicates that modern breeding has a stronger selection effect on silk yield traits than domestication, and imply that silkworm breeding induces aggregation of genes related to silk yield.

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.

Wnt is a family of conserved glycoproteins that participate in a variety of important biological processes including embryo development, cell proliferation and differentiation, and tissue regeneration. The Wnt family is a metazoan novelty found in all animal phyla. Studies have revealed that the number of Wnt genes varies among species, presumably due to reproduction and loss of genes during evolution. However, a comprehensive inventory of Wnt genes in Lepidoptera is lacking. In this study, we identified the repertoire of Wnt genes in the silkworm and seven other species of Lepidoptera and obtained eight Wnt genes (Wnt1, Wnt5–Wnt7, Wnt9–Wnt11, and WntA) in each species. Four of these Wnt genes are clustered in two orientations (5′-Wnt9-Wnt1-Wnt6-Wnt10-3′ and 5′-Wnt10-Wnt6-Wnt1-Wnt9-3′) in both moths and butterflies. Transcript analysis of Wnt in silkworm embryonic stages showed that each BmWnt gene had a unique expression pattern during embryological development. Analysis of a larval stage revealed differential expression of Wnt family members in diverse tissues. Our study provides an overview of the Wnt family in Lepidoptera and will inspire further functional study of the Wnt genes in the silkworm.

WD40 proteins are scaffolding molecules in protein-protein interactions and play crucial roles in fundamental biological processes. Genome-wide characterization of WD40 proteins in animals has been conducted solely in humans. We retrieved 172 WD40 protein genes in silkworm (BmWD40s) and identified these genes in 7 other insects, 9 vertebrates and 5 nematodes. Comparative analysis revealed that the WD40 protein gene family underwent lineage-specific expansions during animal evolution, but did not undergo significant expansion during insect evolution. The BmWD40s were categorized into five clusters and 12 classes according to the phylogenetic classification and their domain architectures, respectively. Sequence analyses indicated that tandem and segmental duplication played minor roles in producing the current number of BmWD40s, and domain recombination events of multi-domain BmWD40s might have occurred mainly after gene duplication events. Gene Ontology (GO) analysis revealed that a higher proportion of BmWD40s was involved in processes, such as binding, transcription-regulation and cellular component biogenesis, compared to all silkworm genes annotated in GO. Microarray-based analysis demonstrated that many BmWD40s had tissue-specific expression and exhibited high and/or sex-related expression during metamorphosis. These findings contribute to a better understanding of the evolution of the animal WD40 protein family and assist the study of the functions of BmWD40s.

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.

Diapause is one of the survival strategies of insects for confronting adverse environmental conditions. Bombyx mori displays typical embryonic diapause, and offspring diapause depends on the incubation environment of the maternal embryo in the bivoltine strains of the silkworm. However, the molecular mechanisms of the diapause induction process are still poorly understood. In this study, we compared the differentially expressed miRNAs (DEmiRs) in bivoltine silkworm embryos incubated at diapause- (25 °C) and non-diapause (15 °C)-inducing temperatures during the blastokinesis (BK) and head pigmentation (HP) phases using transcriptome sequencing. There were 411 known miRNAs and 71 novel miRNAs identified during the two phases. Among those miRNAs, there were 108 and 74 DEmiRs in the BK and HP groups, respectively. By the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the predicted target genes of the DEmiRs, we found that aside from metabolism, the targets were also enriched in phototransduction-fly and insect hormone biosynthesis in the BK group and the HP group, respectively. Dual luciferase reporter assay illustrated that bmo-miR-6497-3p directly regulated Bmcycle and subsequently regulated the expression of circadian genes. These results imply that microRNAs, as vitally important regulators, respond to different temperatures and participate in the diapause induction process across species.