Identification of genes associated with the silk gland size using multi-omics in silkworm (Bombyx mori)

Silk gland size in silkworms (Bombyx mori) affects silk output. However, the molecular mechanisms by which genes regulate silk gland size remain unclear. In this study, silk glands from three pure silkworm strains (A798, A306 and XH) with different silk gland weight phenotypes were compared using transcriptomics and proteomics to identify differentially expressed genes (DEGs) and proteins (DEPs). When comparing A798 to A306 and A798 to XH, 830 and 469 DEGs were up-regulated, respectively. These genes were related to the gene ontology terms, metabolic process, transport activity and biosynthesis process. In addition, 372 and 302 up-regulated differentially expressed proteins were detected in A798 to A306 and A798 to XH, respectively, related to the gene ontology terms, ribosome and protein export, ribosome and polypeptide biosynthesis processes. Moreover, combined transcriptomics, proteomics and weighted correlation network analyses showed that five genes (BGIBMGA002524, BGIBMGA002629, BGIBMGA005659, BGIBMGA005711 and BGIBMGA010889) were significantly associated with the silk gland weight. Reverse Transcription-quantitative real-time Polymerase Chain Reaction (RT-qPCR) and Enzyme linked immunosorbent assay (ELISA) were used to verify the mRNA and protein expression of five genes in the silk glands and tissues of 18 silkworm strains. The results showed that four genes have higher expression levels in heavier silk glands. These genes are associated with glycogen metabolism, fatty acid synthesis and branched chain amino acid metabolism, thus potentially promoting growth and silk protein synthesis. These findings provide valuable insights into the molecular mechanisms underlying the relationship between silk gland weight and silk yield in silkworms.