The evolutionary road from wild moth to domestic silkworm

Structural insight into the substrate specificity of Bombyx mori β-fructofuranosidase belonging to the glycoside hydrolase family 32

Sucrose-hydrolyzing enzymes are largely divided into β-fructofuranosidase and sucrose α-glucosidase. The domestic silkworm Bombyx mori possesses both enzymes, BmSUC1 and BmSUH, belonging to the glycoside hydrolase family 32 (GH32) and GH13, respectively. BmSUC1 was presumed to be acquired by horizontal gene transfer from bacteria based on phylogenetic analysis and related to tolerance to sugar-mimic alkaloids contained in mulberry latex. Here we investigated the substrate specificity of recombinant BmSUC1 that can hydrolyze not only sucrose but also fructooligosaccharides and fructans, and revealed that the enzyme was competitively inhibited by 1,4-dideoxy-1,4-imino-D-arabinitol, one of the alkaloids. Moreover, the crystal structures of BmSUC1 in apo form and complex with sucrose were determined, and the active site pocket was shallow and suitable for shorter substrates but was related to more relaxed substrate specificity than the strict sucrose α-glucosidase BmSUH. Considering together with the distribution of BmSUC1-orthologous genes in many lepidopterans, our results suggest that BmSUC1 contributes to the digestion of fructooligosaccharides and fructans derived from feed plants.

Chromosome-level genome of the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae) provides a resource for evolutionary studies on moths

The peach fruit moth (PFM), Carposina sasakii Matsumura, is a major phytophagous orchard pest widely distributed across Northeast Asia. Here, we report the chromosome-level genome for the PFM, representing the first genome for the family Carposinidae, from the lepidopteran superfamily Copromorphoidea. The genome was assembled into 404.83 Mb sequences using PacBio long-read and Illumina short-read sequences, including 275 contigs, with a contig N50 length of 2.62 Mb. All contigs were assembled into 31 linkage groups assisted by the Hi-C technique, including 30 autosomes and a Z chromosome. BUSCO analysis showed that 98.3% of genes were complete and 0.4% of genes were fragmented, while 1.3% of genes were missing in the assembled genome. In total, 21,697 protein-coding genes were predicted, of which 84.80% were functionally annotated. Because of the importance of diapause triggered by photoperiod in PFM, five circadian genes in the PFM as well as in the other related species were annotated, and potential genes related to diapause and photoperiodic reaction were also identified from transcriptome sequencing. In addition, manual annotation of detoxification gene families was undertaken and showed a higher number of glutathione S-transferase (GST) gene in PFM than in most other lepidopterans, in contrast to a lower number of uridine diphosphate (UDP)-glycosyltransferase (UGT) gene, carboxyl/cholinesterases (CCE) gene and cytochrome P450 monooxygenase (P450) gene, suggesting different detoxication pathways in this moth. The high-quality genome provides a resource for comparative evolutionary studies of this moth and its relatives within the context of radiations across Lepidoptera.

Characterization and potential application of an α-amylase (BmAmy1) selected during silkworm domestication

Efficient resource utilization plays a central role in the high productivity of domesticated plants and animals. Whether artificial selection acts on digestive enzymes in the domesticated silkworm (Bombyx mori), which is larger than its wild ancestor, Bombyx mandarina (B. mandarina), remains unknown. In this study, we present the characteristics of a novel alpha-amylase, BmAmy1, in B. mori. The activity of recombinant BmAmy1 was maximal at 35 °C and pH 9.0, and could be suppressed by amylase inhibitors from mulberry, the exclusive food source of silkworms. Three different transposable element fragments, which were independently inserted in the 5'-upstream regulatory region, might be responsible for the enhanced expression of BmAmy1 in different domesticated silkworm strains as revealed by dual-luciferase reporter assay. The BmAmy1 overexpression increased the weight of female and male B. mori by 11.9% and 6.8%, respectively, compared with non-transgenic controls. Our results emphasize that, by exploring the genetic mechanisms of human-selected traits, the domestication process could be further accelerated through genetic engineering and targeted breeding.

Identification and effect of Zf-AD-containing C2H2 zinc finger genes on BmNPV replication in the silkworm (Bombyx mori)

Zf-AD-containing C2H2 zinc -finger genes (ZAD) are uniquely present and have lineage-specific expansion in arthropods. Arthropods are also the hosts of Baculoviruses. We studied the possible relationship between the lineage-specific expansion of ZAD genes and arthropod-Baculovirus co-evolution. We used the silkworm (Bombyx mori) as a model. We identified 73 ZAD genes (BmZAD) in the silkworm. Sequence-based similarity analysis showed that nine clusters involving 28 BmZADs may have undergone species-specific expansion in the silkworm. Expression pattern analysis showed that the BmZADs were divided into five groups. Group I comprised 10 genes with high expression in multiple tissues, suggesting that BmZADs may play roles in the development of various tissues. We identified six BmZADs that could be induced by the Nucleopolyhedrovirus (BmNPV). Among them, BmZAD69 expression is capable of responding to BmNPV infection, and the ZAD domain is indispensable for the function of BmZAD69 in BmNPV replication. We also detected a 3 bp deletion at 1.7 kb upstream of BmZAD69, which may make it more sensitive to BmNPV infection, and thus elevate the BmNPV resistance in Qiufeng_N, a strain with strong virus resistance. These data suggest that BmZADs may be involved in BmNPV infection and that ZAD genes may play a role in arthropod-Baculovirus co-evolution.

Identification, expression, and artificial selection of silkworm epigenetic modification enzymes

Background

Understanding the genetic basis of phenotype variations during domestication and breeding is of great interest. Epigenetics and epigenetic modification enzymes (EMEs) may play a role in phenotypic variations; however, no comprehensive study has been performed to date. Domesticated silkworm (Bombyx mori) may be utilized as a model in determining how EMEs influence domestication traits.

Results

We identified 44 EMEs in the genome of silkworm (Bombyx mori) using homology searching. Phylogenetic analysis showed that genes in a subfamily among different animals were well clustered, and the expression pattern of EMEs is constant among Bombyx mori, Drosophila melanogaster, and Mus musculus. These are most highly expressed in brain, early embryo, and internal genitalia. By gene-related selective sweeping, we identified five BmEMEs under artificial selection during the domestication and breeding of silkworm. Among these selected genes, BmSuv4–20 and BmDNMT2 harbor selective mutations in their upstream regions that alter transcription factor-binding sites. Furthermore, these two genes are expressed higher in the testis and ovary of domesticated silkworm compared to wild silkworms, and correlations between their expression pattern and meiosis of the sperm and ova were observed.

Conclusions

The domestication of silkworm has induced artificial selection on epigenetic modification markers that may have led to phenotypic changes during domestication. We present a novel perspective to understand the genetic basis underlying animal domestication and breeding.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07155-z.

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

Holometabolous insects have distinct larval, pupal, and adult stages. The pupal stage is typically immobile and can be subject to predation, but cocoon offers pupal protection for many insect species. The cocoon provides a space in which the pupa to adult metamorphosis occurs. It also protects the pupa from weather, predators and parasitoids. Silk protein is a precursor of the silk used in cocoon construction. We used the silkworm as a model species to identify genes affecting silk protein synthesis and cocoon construction. We used quantitative genetic analysis to demonstrate that β-1,4-N-acetylglucosaminidase 1 (BmGlcNase1) is associated with synthesis of sericin, the main composite of cocoon. BmGlcNase1 has an expression pattern coupled with silk gland development and cocoon shell weight (CSW) variation, and CSW is an index of the ability to synthesize silk protein. Up-regulated expression of BmGlcNase1 increased sericin content by 13.9% and 22.5% while down-regulation reduced sericin content by 41.2% and 27.3% in the cocoons of females and males, respectively. Genomic sequencing revealed that sequence variation upstream of the BmGlcNase1 transcriptional start site (TSS) is associated with the expression of BmGlcNase1 and CSW. Selective pressure analysis showed that GlcNase1 was differentially selected in insects with and without cocoons (ω₁ = 0.044 vs. ω₂ = 0.154). This indicates that this gene has a conserved function in the cocooning process of insects. BmGlcNase1 appears to be involved in sericin synthesis and silkworm cocooning.

The cocoon provides a protected space for the metamorphosis of many insect species. Silk protein is a precursor of the fiber used for cocoon construction. Deciphering the genetic basis underlying silk protein synthesis will improve our understanding of cocoon construction and the adaptations of species that construct cocoons. We used the silkworm (Bombyx mori) as a model to identify genes affecting silk protein synthesis and cocoon construction. Quantitative genetic analysis was used to show that β-1,4-N-acetylglucosaminidase 1 (BmGlcNase1), a gene selected during silkworm domestication and breeding, is associated with sericin synthesis. Transgenic-based functional validation confirmed that BmGlcNase1 positively regulates sericin content in the silkworm cocoon. The selective pressure of GlcNase1 in the evolution of insects with cocoons is higher than those without cocoons. This indicates that it has a conserved function in the cocooning process. These results reveal aspects of the genetic basis of silk protein synthesis and the cocoon construction of insects.

SilkDB 3.0: visualizing and exploring multiple levels of data for silkworm

SilkDB is an open-accessibility database and powerful platform that provides comprehensive information on the silkworm (Bombyx mori) genome. Since SilkDB 2.0 was released 10 years ago, vast quantities of data about multiple aspects of the silkworm have been generated, including genome, transcriptome, Hi-C and pangenome. To visualize data at these different biological levels, we present SilkDB 3.0 (https://silkdb.bioinfotoolkits.net), a visual analytic tool for exploring silkworm data through an interactive user interface. The database contains a high-quality chromosome-level assembly of the silkworm genome, and its coding sequences and gene sets are more accurate than those in the previous version. SilkDB 3.0 provides a view of the information for each gene at the levels of sequence, protein structure, gene family, orthology, synteny, genome organization and gives access to gene expression information, genetic variation and genome interaction map. A set of visualization tools are available to display the abundant information in the above datasets. With an improved interactive user interface for the integration of large data sets, the updated SilkDB 3.0 database will be a valuable resource for the silkworm and insect research community.

SGID: a comprehensive and interactive database of the silkworm

Although the domestic silkworm (Bombyx mori) is an important model and economic animal, there is a lack of comprehensive database for this organism. Here, we developed the silkworm genome informatics database (SGID). It aims to bring together all silkworm-related biological data and provide an interactive platform for gene inquiry and analysis. The function annotation in SGID is thorough and covers 98% of the silkworm genes. The annotation details include function description, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes pathway, subcellular location, transmembrane topology, protein secondary/tertiary structure, homologous group and transcription factor. SGID provides genome-scale visualization of population genetics test results based on high-depth resequencing data of 158 silkworm samples. It also provides interactive analysis tools of transcriptomic and epigenomic data from 79 NCBI BioProjects. SGID will be extremely useful to silkworm research in the future.

Gut microbiota metabolic potential correlates with body size between mulberry-feeding lepidopteran pest species

BACKGROUND

Many insect pests rely on microbial symbionts to obtain nutrients or for defence, thereby allowing them to exploit novel food sources and degrade environmental xenobiotics, including pesticides. Although Lepidoptera is one of the most diverse insect taxa and includes important agricultural pests, lepidopteran microbiotas, particularly functional traits, have not been studied widely. Here, we provide a comprehensive characterization of the gut microbiota across multiple mulberry-feeding lepidopteran species, resolving both community structure and metabolic potential.

RESULTS

Our results indicate abundant bacteria inside the gut of larval Lepidoptera. However, even though they were fed the same diet, the structures of the bacterial communities differed in four major mulberry pest species, suggesting host-specific effects on microbial associations. Community-level metabolic reconstructions further showed that although taxonomic composition varied greatly, carbohydrate and amino acid metabolism and membrane transporter were key functional capabilities of the gut bacteria in all samples, which may play basic roles in the larval gut. In addition, principal coordinate analysis (PCoA) of gut bacterial-predicted gene ontologies revealed specialized features of the microbiota associated with these mulberry pests, which were divided into two distinct clusters (macrolepidopterans and microlepidopterans). This pattern became even more prominent when further Lepidoptera species were involved.

CONCLUSIONS

A suite of gut microbiota metabolic functions significantly correlated with larval size; the metabolism of terpenoids and polyketides, xenobiotics biodegradation and metabolism were specifically enriched in large species, while small larvae had enhanced nucleotide metabolism. Our report paves the way for uncovering the correlation between host phenotype and microbial symbiosis in this notorious insect pest group. © 2019 Society of Chemical Industry.

Genetic and genomic analysis for cocoon yield traits in silkworm

Domestic species provides a powerful model for examining genetic mechanisms in the evolution of yield traits. The domestic silkworm (Bombyx mori) is an important livestock species in sericulture. While the mechanisms controlling cocoon yield are largely unknown. Here, using B. mori and its wild relative B. mandarina as intercross parents, 100 BC₁ individuals were sequenced by restriction site-associated DNA sequencing (RAD-Seq). The linkage map contained 9,632 markers was constructed. We performed high-resolution quantitative trait locus (QTL) mapping for four cocoon yield traits. A total of 11 QTLs were identified, including one yield-enhancing QTL from wild silkworm. By integrating population genomics and transcriptomic analysis with QTLs, some favourable genes were revealed, including 14 domestication-related genes and 71 differentially expressed genes (DEGs) in the fifth-instar larval silk gland transcriptome between B. mori and B. mandarina. The relationships between the expression of two important candidate genes (KWMTBOMO04917 and KWMTBOMO12906) and cocoon yield were supported by quantitative real-time PCR (qPCR). Our results provide some new insights into the molecular mechanisms of complex yield traits in silkworm. The combined method might be an efficient approach for identifying putative causal genes in domestic livestock and wild relatives.

A Comparison of Co-expression Networks in Silk Gland Reveals the Causes of Silk Yield Increase During Silkworm Domestication

Long-term domestication and selective breeding have increased the silk yield of the domestic silkworm (Bombyx mori) by several times the amount of the silk yield of its wild ancestor (Bombyx mandarina). However, little is known about the molecular mechanisms behind the increase in silk yield during domestication. Based on dynamic patterns of functional divergence in the silk gland between domestic and wild silkworms, we found that at early and intermediate stages of silk gland development, the up-regulated genes of the domestic silkworm were mainly involved in DNA integration, nucleic acid binding, and transporter activity, which are related to the division and growth of cells. This has led to the posterior silk gland (PSG) of the domestic silkworm having significantly more cells (“factories” of fibroin protein synthesis) than that of the wild silkworm. At the late stage of silk gland development, the up-regulated genes in the domestic silkworm was enriched in protein processing and ribosome pathways, suggesting protein synthesis efficiency is greatly improved during silkworm domestication. While there was an increase in fibroin protein synthesis, the production of sericin protein was simultaneously reduced in the silk gland of the domestic silkworm. This reflects that domestic and wild silkworms have been under different selection pressures. Importantly, we found that the network co-expressed with the silk-coding genes of the domestic silkworm was larger than that of the wild silkworm. Furthermore, many more genes co-expressed with silk-coding genes in the domestic silkworm were subjected to artificial selection than those in the wild silkworm. Our results revealed that the increase of silk yield during silkworm domestication is involved in improvement of a biological system which includes not only expansion of “factories” (cells of PSG) of protein synthesis, but also a high expression of silk-coding genes and silk production-related genes such as biological energy, transport, and ribosome pathway genes.

Most lepidopteran neuroparsin genes seem functional, but in some domesticated silkworm strains it has a fatal mutation

The primary sequence of the Arthropod neurohormone neuroparsin is so variable that so far no orthologs from moths and butterflies have been characterized, even though classical neurosecretory stains identify cells that are homologous to those producing this hormone in other insect species. Here Lepidopteran cDNAs showing limited sequence similarity to other insect neuroparsins are described. That these cDNAs do indeed code for authentic neuroparsins was confirmed by in situ hybridization in the wax moth, Galleria mellonella, which labeled the neuroparsin neuroendocrine cells. Although in virtually all genome assemblies from Lepidoptera a neuroparsin gene could be identified, the genome assembly from the silkworm, Bombyx mori, has a neuroparsin gene containing a 16 nucleotide deletion that renders this gene nonfunctional. Although only a small number of all silkworm strains carry this deletion, it suggests that the domestication of the silkworm has rendered the function of this neurohormone dispensable.

Chromosomal-level reference genome of Chinese peacock butterfly (Papilio bianor) based on third-generation DNA sequencing and Hi-C analysis

BACKGROUND

Papilio bianor Cramer, 1777 (commonly known as the Chinese peacock butterfly) (Insecta, Lepidoptera, Papilionidae) is a widely distributed swallowtail butterfly with a wide number of geographic populations ranging from the southeast of Russia to China, Japan, India, Vietnam, Myanmar, and Thailand. Its wing color consists of both pigmentary colored scales (black, reddish) and structural colored scales (iridescent blue or green dust). A high-quality reference genome of P. bianor is an important foundation for investigating iridescent color evolution, phylogeography, and the evolution of swallowtail butterflies.

FINDINGS

We obtained a chromosome-level de novo genome assembly of the highly heterozygous P. bianor using long Pacific Biosciences sequencing reads and high-throughput chromosome conformation capture technology. The final assembly is 421.52 Mb on 30 chromosomes (29 autosomes and 1 Z sex chromosome) with 13.12 Mb scaffold N50. In total, 15,375 protein-coding genes and 233.09 Mb of repetitive sequences were identified. Phylogenetic analyses indicated that P. bianor separated from a common ancestor of swallowtails ∼23.69-36.04 million years ago. Demographic history suggested that the population expansion of this species from the last interglacial period to the last glacial maximum possibly resulted from its decreased natural enemies and its adaptation to climate change during the glacial period.

CONCLUSIONS

We present a high-quality chromosome-level reference genome of P. bianor using long-read single-molecule sequencing and Hi-C-based chromatin interaction maps. Our results lay the foundation for exploring the genetic basis of special biological features of P. bianor and also provide a useful data source for comparative genomics and phylogenomics among butterflies and moths.

Synergism of open chromatin regions involved in regulating genes in Bombyx mori

The dynamic variability of transcription factors (TFs) and their binding sites makes it challenging to conduct genome-wide transcription regulation research. The silkworm Bombyx mori, which produces silk, is one of the most valuable model insects in the order Lepidoptera. The "opening" and "closing" of chromatin in different silk yield strains is associated with changes in silk production, making this insect a good model for studying the transcriptional regulation of genes. However, few studies have examined the open chromatin regions (OCRs) of silkworms, and studying OCR synergism and their function in silk production remains challenging. Here, we performed formaldehyde-assisted isolation of regulatory elements (FAIRE) to isolate OCRs from the silk glands of fifth-instar larvae of the DaZao and D872 strains. In total, 128,908 high confidence OCRs were identified and approximately 80% of OCRs were located in non-coding regions. OCRs upregulated adjacent genes and showed signal-dependent vulnerability to single-nucleotide polymorphisms. Mid- and low-signal OCRs were more likely to have single-nucleotide polymorphisms (SNP). Further, OCRs interacted with each other within a distance of 5 kb. We named the OCR interaction complex as the "cluster of related regions" (COREs). The functions of the CORE and its harbored OCRs showed some differences. Additionally, COREs enriched many silk protein synthesis-associated genes, some of which were upregulated. This study identified numerous high confidence regulation sites and synergistic regulatory modes of OCRs that affect adjacent genes. These results provide insight into silkworm transcriptional regulation and improve our understanding of cis-element cooperation.

Artificial selection on storage protein 1 possibly contributes to increase of hatchability during silkworm domestication

Like other domesticates, the efficient utilization of nitrogen resources is also important for the only fully domesticated insect, the silkworm. Deciphering the way in which artificial selection acts on the silkworm genome to improve the utilization of nitrogen resources and to advance human-favored domestication traits, will provide clues from a unique insect model for understanding the general rules of Darwin's evolutionary theory on domestication. Storage proteins (SPs), which belong to a hemocyanin superfamily, basically serve as a source of amino acids and nitrogen during metamorphosis and reproduction in insects. In this study, through blast searching on the silkworm genome and further screening of the artificial selection signature on silkworm SPs, we discovered a candidate domestication gene, i.e., the methionine-rich storage protein 1 (SP1), which is clearly divergent from other storage proteins and exhibits increased expression in the ova of domestic silkworms. Knockout of SP1 via the CRISPR/Cas9 technique resulted in a dramatic decrease in egg hatchability, without obvious impact on egg production, which was similar to the effect in the wild silkworm compared with the domestic type. Larval development and metamorphosis were not affected by SP1 knockout. Comprehensive ova comparative transcriptomes indicated significant higher expression of genes encoding vitellogenin, chorions, and structural components in the extracellular matrix (ECM)-interaction pathway, enzymes in folate biosynthesis, and notably hormone synthesis in the domestic silkworm, compared to both the SP1 mutant and the wild silkworm. Moreover, compared with the wild silkworms, the domestic one also showed generally up-regulated expression of genes enriched in the structural constituent of ribosome and amide, as well as peptide biosynthesis. This study exemplified a novel case in which artificial selection could act directly on nitrogen resource proteins, further affecting egg nutrients and eggshell formation possibly through a hormone signaling mediated regulatory network and the activation of ribosomes, resulting in improved biosynthesis and increased hatchability during domestication. These findings shed new light on both the understanding of artificial selection and silkworm breeding from the perspective of nitrogen and amino acid resources.

Like other domesticates, nitrogen resources are also important for the only fully domesticated insect, the silkworm. Deciphering the way in which artificial selection acts on the silkworm genome to improve the utilization of nitrogen resources, thereby advancing human-favored domestication traits, will provide clues from a unique insect model for understanding the general rules of Darwin's theory on artificial selection. However, the mechanisms of domestication in the silkworm remain largely unknown. In this study, we focused on one important nitrogen resource, the storage protein (SP). We discovered that the methionine-rich storage protein 1 (SP1), which is divergent from other SPs, is the only target of artificial selection. Based on functional evidence, together with key findings from the comprehensive comparative transcriptome, we propose that artificial selection favored higher expression of SP1 in the domestic silkworm, which would influence the genes or pathways vital for egg development and eggshell formation. Artificial selection also consistently favored activated ribosome activities and improved amide and peptide biosynthesis in the ova, like what they may act in the silk gland to increase silk-cocoon yield. We highlighted a novel case in which artificial selection could directly act on a nitrogen resource protein associated with a human-desired domestication trait.

Decrease of gene expression diversity during domestication of animals and plants

BACKGROUND

The genetic mechanisms underlying the domestication of animals and plants have been of great interest to biologists since Darwin. To date, little is known about the global pattern of gene expression changes during domestication.

RESULTS

We generated and collected transcriptome data for seven pairs of domestic animals and plants including dog, silkworm, chicken, rice, cotton, soybean and maize and their wild progenitors and compared the expression profiles between the domestic and wild species. Intriguingly, although the number of expressed genes varied little, the domestic species generally exhibited lower gene expression diversity than did the wild species, and this lower diversity was observed for both domestic plants and different kinds of domestic animals including insect, bird and mammal in the whole-genome gene set (WGGS), candidate selected gene set (CSGS) and non-CSGS, with CSGS exhibiting a higher degree of decreased expression diversity. Moreover, different from previous reports which found 2 to 4% of genes were selected by human, we identified 6892 candidate selected genes accounting for 7.57% of the whole-genome genes in rice and revealed that fewer than 8% of the whole-genome genes had been affected by domestication.

CONCLUSIONS

Our results showed that domestication affected the pattern of variation in gene expression throughout the genome and generally decreased the expression diversity across species, and this decrease may have been associated with decreased genetic diversity. This pattern might have profound effects on the phenotypic and physiological changes of domestic animals and plants and provide insights into the genetic mechanisms at the transcriptome level other than decreased genetic diversity and increased linkage disequilibrium underpinning artificial selection.