The beta-1, 4-N-acetylglucosaminidase 1 gene, selected by domestication and breeding, is involved in cocoon construction of Bombyx mori

QTL analysis to identify genes involved in the trade-off between silk protein synthesis and larva-pupa transition in silkworms

BACKGROUND

Insect-based food and feed are increasingly attracting attention. As a domesticated insect, the silkworm (Bombyx mori) has a highly nutritious pupa that can be easily raised in large quantities through large-scale farming, making it a highly promising source of food. The ratio of pupa to cocoon (RPC) refers to the proportion of the weight of the cocoon that is attributed to pupae, and is of significant value for edible utilization, as a higher RPC means a higher ratio of conversion of mulberry leaves to pupa. In silkworm production, there is a trade-off between RPC and cocoon shell ratiao(CSR), which refers the ratio of silk protein to the entire cocoon, during metamorphosis process. Understanding the genetic basis of this balance is crucial for breeding edible strains with a high RPC and further advancing its use as feed.

RESULTS

Using QTL-seq, we identified a quantitative trait locus (QTL) for the balance between RPC and CSR that is located on chromosome 11 and covers a 9,773,115-bp region. This locus is an artificial selection hot spot that contains ten non-overlapping genomic regions under selection that were involved in the domestication and genetic breeding processes. These regions include 17 genes, nine of which are highly expressed in the silk gland, which is a vital component in the trade-off between RPC and CSR. These genes are annotate with function related with epigenetic modifications and the regulation of DNA replication et al. We identified one and two single nucleotide polymorphisms (SNPs) in the exons of teh KWMTBOMO06541 and KWMTBOMO06485 genes that result in amino acid changes in the protein domains. These SNPs have been strongly selected for during the domestication process. The KWMTBOMO06485 gene encodes the Bombyx mori (Bm) tRNA methyltransferase (BmDnmt2) and its knockout results in a significant change in the trade-off between CSR and RPC in both sexes.

CONCLUSIONS

Taken together, our results contribute to a better understanding of the genetic basis of RPC and CSR. The identified QTL and genes that affect RPC can be used for marker-assisted and genomic selection of silkworm strains with a high RPC. This will further enhance the production efficiency of silkworms and of closely-related insects for edible and feed purposes.

Mutation in the Bombyx mori BmGMC2 gene impacts silk production and silk protein synthesis

Silk is a natural protein fiber that is predominantly comprised of fibroin and sericin. In addition, it contains seroins, protease inhibitors, enzymes, and other proteins. We found an ecdysone oxidase BmGMC2, notably, which is specifically and highly expressed only in the silk glands of silkworms (Bombyx mori L.). It is also one of the main components of non-cocoon silk, however, its precise function remains unclear. In this study, we examined the spatiotemporal expression pattern of this protein and obtained a homozygous mutant strain (K-GMC2) using the CRISPR-Cas9 system. Compared to the wild-type strain (WT), the silk production and main silk proteins significantly decreased in the larval stage, and the adhesive strength of native silk proteins decreased in the final instar. Proteomic data indicated the abundance of ribosomal proteins decreased significantly in K-GMC2, differentially expressed proteins (DEPs) were enriched in pathways related to neurodegenerative diseases and genetic information processing, indicating that knockout may lead to a certain degree of cell stress, affecting the synthesis of silk proteins. This study investigated the expression pattern and gene function of ecdysone oxidase BmGMC2 in silk and silk glands, laying the groundwork for understanding the role of enzymes in the production of silk fibers.

Fatty Acid Desaturase Bmdesat5, Suppressed in the Salivary Glands by Domestication, is Involved in Regulation of Food Intake in Silkworms

Understanding the evolutionary genetics of food intake regulation in domesticated animals has relevance to evolutionary biology, animal improvement, and obesity treatment. Here, we observed that the fatty acid desaturase gene (Bmdesat5), which regulates food intake, is suppressed in domesticated silkworms, but expressed in the salivary glands of the wild silkworm Bombyx mandarina. The content of its catalytic product, cis-vaccenic acid, was related to the expression levels of Bmdesat5 in the salivary glands of domesticated and wild silkworm strains. These two strains also showed significant differences in food intake. Using orally administering cis-vaccenic acid and transgenic-mediated overexpression, we verified that cis-vaccenic acid functions as a satiation signal, regulating food intake and growth in silkworms. Selection analysis showed that Bmdesat5 experienced selection, especially in the potential promoter, 5'-untranslated, and intron regions. This study highlights the importance of the decrement of satiety in silkworm domestication and provides new insights into the potential involvement of salivary glands in the regulation of satiety in animals, by acting as a supplement to gut-brain nutrient signaling.

Cold-stored mulberry leaves affect antioxidant system and silk proteins of silkworm (Bombyx mori) larva

BACKGROUND

Cold storage has been widely used to maintain the quality of vegetables, but whether eating cold-stored vegetables affects health remains unknown.

RESULTS

This study used silkworms as an animal model to evaluate the effects of nutrient changes in cold-stored mulberry leaves (CSML) on health. Compared with fresh mulberry leaves (FML), CSML contained lower vitamin C, soluble sugars and proteins, and higher H2 O2 , suggesting decreased antioxidant ability and nutrition. The CSML did not obviously affect larval survival rate, body weight or dry matter rate, cocoon shape, weight and size, or final rates of cluster and cocooning relative to the FML, suggesting CSML did not alter overall growth and development. However, the CSML increased the initial rates of cluster and cocooning and upregulated BmRpd3, suggesting CSML shortened larval lifespan and enhanced senescence. CSML upregulated BmNOX4, downregulated BmCAT, BmSOD and BmGSH-Px and increased H2 O2 in silkworms, suggesting CSML caused oxidative stress. CSML upregulated ecdysone biosynthesis and inactivation genes and elevated ecdysone concentration in silkworms, suggesting that CSML affected hormone homeostasis. CSML upregulated apoptosis-related genes, downregulated sericin and silk fibroin genes and decreased sericin content rate in silkworms, suggesting oxidative stress and protein deficiency.

CONCLUSION

Cold storage reduced nutrition and antioxidant capability of mulberry leaves. CSML did not influence growth and development of silkworm larva, but affected health by causing oxidative stress and reducing protein synthesis. The findings show that the ingredient changes in CSML had negative effects on health of silkworms. © 2023 Society of Chemical Industry.

Highly accurate genome assembly of an improved high-yielding silkworm strain, Nichi01

The silkworm (Bombyx mori) is an important lepidopteran model insect and an industrial domestic animal traditionally used for silk production. Here, we report the genome assembly of an improved Japanese strain Nichi01, in which the cocoon yield is comparable to that of commercial silkworm strains. The integration of PacBio Sequel II long-read and ddRAD-seq-based high-density genetic linkage map achieved the highest quality genome assembly of silkworms to date; 22 of the 28 pseudomolecules contained telomeric repeats at both ends, and only four gaps were present in the assembly. A total of 452 Mbp of the assembly with an N50 of 16.614 Mbp covered 99.3% of the complete orthologs of the lepidopteran core genes. Although the genome sequence of Nichi01 and that of the previously reported low-yielding tropical strain p50T assured their accuracy in most regions, we corrected several regions, misassembled in p50T, in our assembly. A total of 18,397 proteins were predicted using over 95 Gb of mRNA-seq derived from 10 different organs, covering 96.9% of the complete orthologs of the lepidopteran core genes. The final assembly and annotation files are available in KAIKObase (https://kaikobase.dna.affrc.go.jp/index.html) along with a genome browser and BLAST searching service, which would facilitate further studies and the breeding of silkworms and other insects.

Analysis of histomorphometric and proteome dynamics inside the silk gland lumen of Bombyx mori revealed the dynamic change of silk protein during the molt stage

Silkworms spin different silks at different growth stages for specific purposes. The silk spun before the end of each instar is stronger than that at the beginning of each instar and cocoon silk. However, the compositional changes in silk proteins during this process are unknown. Consequently, we performed histomorphological and proteomic analyses of the silk gland to characterize changes from the instar end to the next instar beginning. The silk glands were collected on day 3 of third- and fourth-instar larvae (III-3 and IV-3) and the beginning of fourth-instar larvae (IV-0). Proteomic analysis identified 2961 proteins from all silk glands. Silk proteins P25 and Ser5 were significantly more abundant in III-3 and IV-3 than in IV-0, and many cuticular proteins and protease inhibitors increased significantly in IV-0 compared with III-3 and IV-3. This shift may cause mechanical property differences between the instar end and beginning silk. Using section staining, qPCR, and western blotting, we found for the first time that silk proteins were degraded first and then resynthesized during the molting stage. Furthermore, we revealed that fibroinase mediated the changes of silk proteins during molting. Our results provide insights into the molecular mechanisms of silk proteins dynamic regulation during molting.

The Role of Chitooligosaccharidolytic β-N-Acetylglucosamindase in the Molting and Wing Development of the Silkworm Bombyx mori

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.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

Gene mapping reveals the association between tyrosine protein kinase Abl1 and the silk yield of Bombyx mori

The Z chromosome of the silkworm contains a major gene that influences silk yield. This major locus on chromosome Z accounts for 35.10% of the phenotypic variance. The location and identification of the gene have been a focus of silkworm genetics research. Unfortunately, identification of this gene has been difficult. We used extreme phenotype subpopulations and selected from a backcross population, BC₁ M, which was obtained using the high-yield strain 872B and the low-yield strain IS-Dazao as parents, for mapping the gene on the chromosome Z. The candidate region was narrowed down to 134 kb at the tip of the chromosome. BmAbl1 in this region correlated with silk gland development by spatiotemporal expression analysis. This gene was differentially expressed in the posterior silk glands of the high- and low-yield strains. In BmAbl1, an insertion-deletion (indel) within the 10th exonic region and an SNP within the 6th intronic region were detected and shown to be associated with cocoon shell weight in 84 Bombyx mori strains with different yields. Nucleotide diversity analysis of BmAbl1 and its 50 kb flanking regions indicated that BmAbl1 has experienced strong artificial selection during silkworm domestication. This study is the first to identify the genes controlling silk yield in the major QTL of the Z chromosome using forward genetics.