The evolutionary road from wild moth to domestic silkworm

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.

Fitness effects of synthetic and natural diet preservatives on the edible insect Bombyx mori

Silkworm pupae as widely consumed insect products are good biosources of protein and micronutrients. Silkworm rearing throughout the year can be achieved by feeding them an artificial diet instead of native plants, facilitating extensive pupa production. However, artificial diets are prone to spoilage caused by bacterial contamination. Here, we evaluated the antiseptic effect of ethylparaben (EP, chemical preservative) and medium-chain fatty acids (MCFA, natural preservative) in a silkworm artificial diet. Results showed that both preservatives effectively inhibited pathogenic bacterial growth. Furthermore, the addition of EP or MCFA did not negatively impact the production capacity of silkworms and the homeostasis of gut microbiota. However, the expression of genes involved in detoxification such as Ugt2, and immune response such as Cecropin B, were upregulated after EP consumption. Therefore, natural preservative MCFA emerges as a suitable option from a safety perspective. These findings highlight future directions for improving insect artificial diet formulation.

Multiple independent origins of the female W chromosome in moths and butterflies

Lepidoptera, the most diverse group of insects, exhibit female heterogamy (Z0 or ZW), which is different from most other insects (male heterogamy, XY). Previous studies suggest a single origin of the Z chromosome. However, the origin of the lepidopteran W chromosome remains poorly understood. Here, we assemble the genome from females down to the chromosome level of a model insect (Bombyx mori) and identify a W chromosome of approximately 10.1 megabase using a newly developed tool. In addition, we identify 3593 genes that were not previously annotated in the genomes of B. mori. Comparisons of 21 lepidopteran species (including 17 ZW and four Z0 systems) and three trichopteran species (Z0 system) reveal that the formation of Ditrysia W involves multiple mechanisms, including previously proposed canonical and noncanonical models, as well as a newly proposed mechanism called single-Z turnover. We conclude that there are multiple independent origins of the W chromosome in the Ditrysia (most moths and all butterflies) of Lepidoptera.

A long-read draft assembly of the Chinese mantis (Mantodea: Mantidae: Tenodera sinensis) genome reveals patterns of ion channel gain and loss across Arthropoda

Praying mantids (Mantodea: Mantidae) are iconic insects that have captivated biologists for decades, especially the species with cannibalistic copulatory behavior. This behavior has been cited as evidence that insects lack nociceptive capacities and cannot feel pain; however, this behaviorally driven hypothesis has never been rigorously tested at the genetic or functional level. To enable future studies of nociceptive capabilities in mantids, we sequenced and assembled a draft genome of the Chinese praying mantis (Tenodera sinensis) and identified multiple classes of nociceptive ion channels by comparison to orthologous gene families in Arthropoda. Our assembly-produced using PacBio HiFi reads-is fragmented (total size = 3.03 Gb; N50 = 1.8 Mb; 4,966 contigs), but is highly complete with respect to gene content (BUSCO complete = 98.7% [odb10_insecta]). The size of our assembly is substantially larger than that of most other insects, but is consistent with the size of other mantid genomes. We found that most families of nociceptive ion channels are present in the T. sinensis genome; that they are most closely related to those found in the damp-wood termite (Zootermopsis nevadensis); and that some families have expanded in T. sinensis while others have contracted relative to nearby lineages. Our findings suggest that mantids are likely to possess nociceptive capabilities and provide a foundation for future experimentation regarding ion channel functions and their consequences for insect behavior.

The genomic history and global migration of a windborne pest

Many insect pests, including the brown planthopper (BPH), undergo windborne migration that is challenging to observe and track. It remains controversial about their migration patterns and largely unknown regarding the underlying genetic basis. By analyzing 360 whole genomes from around the globe, we clarify the genetic sources of worldwide BPHs and illuminate a landscape of BPH migration showing that East Asian populations perform closed-circuit journeys between Indochina and the Far East, while populations of Malay Archipelago and South Asia undergo one-way migration to Indochina. We further find round-trip migration accelerates population differentiation, with highly diverged regions enriching in a gene desert chromosome that is simultaneously the speciation hotspot between BPH and related species. This study not only shows the power of applying genomic approaches to demystify the migration in windborne migrants but also enhances our understanding of how seasonal movements affect speciation and evolution in insects.

Photoperiodism of Diapause Induction in the Silkworm, Bombyx mori

The silkworm Bombyx mori exhibits a photoperiodic response (PR) for embryonic diapause induction. This article provides a comprehensive review of literature on the silkworm PR, starting from early works on population to recent studies uncovering the molecular mechanism. Makita Kogure (1933) conducted extensive research on the PR, presenting a pioneering paper on insect photoperiodism. In the 1970s and 80s, artificial diets were developed, and the influence of nutrition on PR was well documented. The photoperiodic photoreceptor has been investigated from organ to molecular level in the silkworm. Culture experiments demonstrated that the photoperiodic induction can be programmed in an isolated brain (Br)-subesophageal ganglion (SG) complex with corpora cardiaca (CC)-corpora allata (CA). The requirement of dietary vitamin A for PR suggests the involvement of opsin pigment in the photoperiodic reception, and a cDNA encoding an opsin (Boceropsin) was cloned from the brain. The effector system concerning the production and secretion of diapause hormone (DH) has also been extensively investigated in the silkworm. DH is produced in a pair of posterior cells of SG, transported to CC by nervi corporis cardiaci, and ultimately released into the hemolymph. Possible involvement of GABAergic and corazonin (Crz) signal pathways was suggested in the control of DH secretion. Knockout (KO) experiments of GABA transporter (GAT) and circadian clock genes demonstrated that GAT plays a crucial role in PR through circadian control. A model outlining the PR mechanism, from maternal photoperiodic light reception to DH secretion, has been proposed.

Cytotoxin-mediated silk gland organ dysfunction diverts resources to enhance silkworm fecundity by potentiating nutrient-sensing IIS/TOR pathways

Energy reserves, primarily stored in the insect's fat body, are essential for physiological processes such as reproduction and cocoon formation. However, whether these processes are mutually constraining is unknown. Here, we showed that cocoon-free silkworms accumulate amino acid constituents of silk proteins in the hemolymph and maintain lipid and sugar reserves in the pupal fat body by repressing the expression of sericin and fibroin genes in the middle and posterior silk glands, respectively, via butterfly pierisin-1A catalytic domain expression. This, in turn, upregulates insulin/insulin-like signaling and target of rapamycin (IIS/TOR) signaling, which enhances vitellogenesis and accelerates ovarian development, thus contributing to increased fecundity. The impacts of semi-starvation on fecundity and egg hatchability were also less pronounced in cocoon-free silkworms compared with wildtype silkworms. These data uncover the resource allocation trade-off between cocoon formation and fecundity and demonstrate that nutritional signaling plays a role in regulating silkworm reproduction.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

A major endogenous glycoside hydrolase mediating quercetin uptake in Bombyx mori

Quercetin is a common plant flavonoid which is involved in herbivore-plant interactions. Mulberry silkworms (domestic silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina) take up quercetin from mulberry leaves and accumulate the metabolites in the cocoon, thereby improving its protective properties. Here we identified a glycoside hydrolase, named glycoside hydrolase family 1 group G 5 (GH1G5), which is expressed in the midgut and is involved in quercetin metabolism in the domestic silkworm. Our results suggest that this enzyme mediates quercetin uptake by deglycosylating the three primary quercetin glycosides present in mulberry leaf: rutin, quercetin-3-O-malonylglucoside, and quercetin-3-O-glucoside. Despite being located in an unstable genomic region that has undergone frequent structural changes in the evolution of Lepidoptera, GH1G5 has retained its hydrolytic activity, suggesting quercetin uptake has adaptive significance for mulberry silkworms. GH1G5 is also important in breeding: defective mutations which result in discoloration of the cocoon and increased silk yield are homozygously conserved in 27 of the 32 Japanese white-cocoon domestic silkworm strains and 12 of the 30 Chinese ones we investigated.

Single-cell transcriptomes provide insights into expansion of glial cells in Bombyx mori

The diversity of cell types in the brain and how these change during different developmental stages, remains largely unknown. The life cycle of insects is short and goes through 4 distinct stages including embryonic, larval, pupal, and adult stages. During postembryonic life, the larval brain transforms into a mature adult version after metamorphosis. The silkworm, Bombyx mori, is a lepidopteran model insect. Here, we characterized the brain cell repertoire of larval and adult B. mori by obtaining 50 708 single-cell transcriptomes. Seventeen and 12 cell clusters from larval and adult brains were assigned based on marker genes, respectively. Identified cell types include Kenyon cells, optic lobe cells, monoaminergic neurons, surface glia, and astrocyte glia. We further assessed the cell type compositions of larval and adult brains. We found that the transition from larva to adult resulted in great expansion of glial cells. The glial cell accounted for 49.8% of adult midbrain cells. Compared to flies and ants, the mushroom body kenyon cell is insufficient in B. mori, which accounts for 5.4% and 3.6% in larval and adult brains, respectively. Analysis of neuropeptide expression showed that the abundance and specificity of expression varied among individual neuropeptides. Intriguingly, we found that ion transport peptide was specifically expressed in glial cells of larval and adult brains. The cell atlas dataset provides an important resource to explore cell diversity, neural circuits and genetic profiles.

Disruption of a BTB-ZF transcription factor causes female sterility and melanization in the larval body of the silkworm, Bombyx mori

The dilute black (bd) of the silkworm Bombyx mori is a recessive mutant that produces a grayish-black color in the larval integument, instead of the characteristic white color found in wild-type larvae. In addition, eggs produced by bd females are sterile due to a deficiency in the micropylar apparatus. We identified candidate genes responsible for the bd phenotype using publicly available RNA-seq data. One of these candidate genes was homologous to the maternal gene required for meiosis (mamo) of Drosophila melanogaster, which encodes a broad-complex, tramtrack, and bric-à-brac-zinc finger (BTB-ZF) transcription factor essential for female fertility. In three independent bd strains, the expression of the B. mori mamo (Bmmamo) was downregulated in the larval integument. Using a CRISPR/Cas9-mediated knockout strategy, we found that Bmmamo knockout mutants exhibit a grayish-black color in the larval integument and female infertility. Moreover, larvae obtained from the complementation cross between bd/+ mutants and heterozygous knockouts for the Bmmamo also exhibited a grayish-black color, indicating that Bmmamo is responsible for the bd phenotype. Gene expression analysis using Bmmamo knockout mutants suggested that the BmMamo protein suppresses the expression of melanin synthesis genes. Previous comparative genome analysis revealed that the Bmmamo was selected during silkworm domestication, and we found that Bmmamo expression in the larval integument is higher in B. mori than in the wild silkworm B. mandarina, suggesting that the Bmmamo is involved in domestication-associated pigmentation changes of the silkworm.

Adaptive Changes in Detoxification Metabolism and Transmembrane Transport of Bombyx mori Malpighian Tubules to Artificial Diet

The high adaptability of insects to food sources has contributed to their ranking among the most abundant and diverse species on Earth. However, the molecular mechanisms underlying the rapid adaptation of insects to different foods remain unclear. We explored the changes in gene expression and metabolic composition of the Malpighian tubules as an important metabolic excretion and detoxification organ in silkworms (Bombyx mori) fed mulberry leaf and artificial diets. A total of 2436 differentially expressed genes (DEGs) and 245 differential metabolites were identified between groups, with the majority of DEGs associated with metabolic detoxification, transmembrane transport, and mitochondrial function. Detoxification enzymes, such as cytochrome P450 (CYP), glutathione-S-transferase (GST), and UDP-glycosyltransferase, and ABC and SLC transporters of endogenous and exogenous solutes were more abundant in the artificial diet group. Enzyme activity assays confirmed increased CYP and GST activity in the Malpighian tubules of the artificial diet-fed group. Metabolome analysis showed increased contents of secondary metabolites, terpenoids, flavonoids, alkaloids, organic acids, lipids, and food additives in the artificial diet group. Our findings highlight the important role of the Malpighian tubules in adaptation to different foods and provide guidance for further optimization of artificial diets to improve silkworm breeding.

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.

Recessive embryonic lethal mutations uncovered in heterozygous condition in silkworm semiconsomic strains

In this study, we found two embryonic lethal mutations, t04 lethal (l-t04) and m04 lethal (l-m04), in semiconsomic strains T04 and M04, respectively. In these semiconsomic strains, the entire diploid genome, except for one chromosome 4 of the wild silkworm Bombyx mandarina, is substituted with chromosomes of the domesticated silkworm B. mori, and l-t04 and l-m04 mutations are located on B. mandarina-derived chromosome 4. To clarify the cause of the lethalities and the genes responsible for these mutations, positional cloning and CRISPR/Cas9 mediated knockout screening were performed. Finally, genetic complementation tests identified the mutations responsible for the l-t04 and l-m04 as the Bombyx homolog of imaginal discs arrested (Bmida) and TATA box binding protein-associated factor 5 (BmTaf5), respectively. Lethal stages of each knockout mutant indicated the importance of these genes in B. mori late embryogenesis. The lethal mutations responsible for l-t04 and l-m04 were not found in parental strains or wild B. mandarina collected from 39 distinct locations in Japan, indicating that both mutations were independently introduced during or after the development of the semiconsomic strains. We conclude that the recessive embryonic lethality in the T04 and M04 strains is due to deleterious mutations produced in B. mandarina-derived chromosome 4.

Mutation Rate and Spectrum of the Silkworm in Normal and Temperature Stress Conditions

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−8 (95% confidence interval, 0.33 × 10−8–0.49 × 10−8) 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.

FibH Gene Complete Sequences (FibHome) Revealed Silkworm Pedigree

The highly repetitive and variable fibroin heavy chain (FibH) gene can be used as a silkworm identification; however, only a few complete FibH sequences are known. In this study, we extracted and examined 264 FibH gene complete sequences (FibHome) from a high-resolution silkworm pan-genome. The average FibH lengths of the wild silkworm, local, and improved strains were 19,698 bp, 16,427 bp, and 15,795 bp, respectively. All FibH sequences had a conserved 5' and 3' terminal non-repetitive (5' and 3' TNR, 99.74% and 99.99% identity, respectively) sequence and a variable repetitive core (RC). The RCs differed greatly, but they all shared the same motif. During domestication or breeding, the FibH gene mutated with hexanucleotide (GGTGCT) as the core unit. Numerous variations existed that were not unique to wild and domesticated silkworms. However, the transcriptional factor binding sites, such as fibroin modulator-binding protein, were highly conserved and had 100% identity in the FibH gene's intron and upstream sequences. The local and improved strains with the same FibH gene were divided into four families using this gene as a marker. Family I contained a maximum of 62 strains with the optional FibH (Opti-FibH, 15,960 bp) gene. This study provides new insights into FibH variations and silkworm breeding.

Non-classical digestive lipase BmTGL selected by gene amplification reduces the effects of mulberry inhibitor during silkworm domestication

Efficient utilization of dietary lipids is crucial for Bombyx mori, also known as domesticated silkworms. However, the effects of domestication on the genes encoding lipases remain unknown. In this study, we investigated the expression difference of one triacylglycerol lipase (BmTGL) between B.mori and wild (ancestor) silkworm strains (Bombyx mandarina). An immunofluorescence localization analysis showed that BmTGL was present in all parts of the gut and was released into the intestinal lumen. BmTGL expression was significantly enhanced in different domesticated silkworm strains compared to that in the B. mandarina strains. The BmTGL copy numbers in the genomes of the domesticated silkworm strains were 2-to-3 folds that of the B. mandarina strains and accounted for the enhanced expression of BmTGL in the domesticated silkworm strains. The Ser144Asn substitution in the Ser-Asp-His catalytic triads of BmTGL resulted in relatively lower lipase activity and reduced sensitivity to the lipase inhibitor morachalcone A. Moreover, BmTGL overexpression significantly increased the weights of the B. mori silkworms compared to those of the non-transgenic controls. Thus, the selection of BmTGL by gene amplification may be a trade-off between maintaining high enzymatic activity and reducing the effects of mulberry inhibitors during silkworm domestication.

Effect of environmental conditions on the changes of voltinism in three lines of Bombyx mori

This work aims to estimate the voltinism differences in the studied lines and the influence of biological and environmental variables on this trait. Three Bombyx mori. lines were used. One polyvoltine (Indian, C. Nichi) and two bivoltine (Chinese6-C6 and Japanese7-J7). Silkworms were reared under controlled conditions of temperature (24±1, 26±1, 28±1°C) and photoperiod (8-16, 12-12, 16-8 light/darkness) during the life cycle. The measured biological characters were duration of larval period (DLP), larval weight on fifth day of V instar (LW(V-5)), rate of weight (RWV), duration of life cycle (DLC), fecunditiy (NE), fertility (F) and voltinism (V). Voltinism was estimated using a linear regression model; the variables that influence voltinism were found using a multiple regression analysis; and the effect of environmental conditions on change of voltinism employing an ANOVA. The results show that the variables that affect voltinism are different for each line, which demonstrates the effect of the line factor. Nichi voltinism was affected by the six measured variables, while in C6 and J7 where voltinisms depends on the duration of diapause, some variables were excluded from the model. In Indian tropical polyvoltine, the temperature was the influencing factor, while in Chinese and Japanese bivoltines lines was the photoperiod.

Effects of Microbial Transfer during Food-Gut-Feces Circulation on the Health of Bombyx mori

Change in habitual diet may negatively affect health. The domestic silkworm (Bombyx mori) is an economically important oligophagous insect that feeds on mulberry leaves. The growth, development, and immune-disease resistance of silkworms have declined under artificial dietary conditions. In this study, we used B. mori as a model insect to explore the relationship between changes in diet and balance of intestinal microbes due to its simpler guts compared with those of mammals. We found that artificial diets reduced the intestinal bacterial diversity in silkworms and resulted in a simple intestinal microbial structure. By analyzing the correlations among food, gut, and fecal microbial diversity, we found that an artificial diet was more easily fermented and enriched the lactic acid bacteria in the gut of the silkworms. This diet caused intestinal acidification and microbial imbalance (dysbiosis). When combined with the artificial diet, Enterococcus mundtii, a colonizing opportunistic pathogen, caused dysbiosis and allowed the frequent outbreak of bacterial diseases in the silkworms. This study provides further systematic indicators and technical references for future investigations of the relationship between diet-based environmental changes and intestinal microbial balance. IMPORTANCE The body often appears unwell after habitual dietary changes. The domestic silkworm (Bombyx mori) raised on artificial diets is a good model to explore the relationship between dietary changes and the balance of intestinal microbes. In this study, the food-gut-feces microbial model was established, and some potential key genera that could regulate the balance of intestinal microbiota were screened out. Our findings will provide a reference for future research to further our understanding of healthy silkworm development and may even be useful for similar research on other animals.

Comparative Silk Transcriptomics Illuminates Distinctive Impact of Artificial Selection in Silkworm Modern Breeding

Early domestication and the following improvement are two important processes in the cocoon silk evolution of silkworms. In contrast to early domestication, understanding of the improvement process is still fuzzy. By systematically comparing the larval silk gland transcriptomes of the wild, early domestic, and improved silkworms, we highlighted a novel landscape of transcriptome in the silk glands of improved ones. We first clarified that silk cocoon protein genes were up-regulated in modern breeding but not in early domestication. Furthermore, we found that differentially expressed genes (DEGs) between improved and early domestic silkworms (2711), as well as between improved and wild silkworms (2264), were obviously more than those between the early domestic and wild silkworms (158), with 1671 DEGs specific in the improved silkworm (IS-DEGs). Hierarchical clustering of all the DEGs consistently indicated that improved silkworms were significantly diverged from the early domestic and wild silkworms, suggesting that modern breeding might cause prompt and drastic dynamic changes of gene expression in the silk gland. We further paid attention to these 1671 IS-DEGs and were surprised to find that down-regulated genes were enriched in basic organonitrogen compound biosynthesis, RNA biosynthesis, and ribosome biogenesis processes, which are generally universally expressed, whereas those up-regulated genes were enriched in organonitrogen compound catabolic processes and functions involving in the dynamic regulation of protein post-translation of modification. We finally highlighted one candidate improvement gene among these up-regulated IS-DEGs, i.e., GDAP2, which may play roles in silk behavior and the overall robustness of the improved silkworm. The findings strongly suggest that modern breeding may facilitate effective control of the basic consumption of nitrogen and a stronger switch of nitrogen resources from other tissues to the silk glands, for an efficient supply for silk production, and implies the importance of brain behavior and robustness in silk yield improvement of modern breeding.

Allele-specific knockouts reveal a role for apontic-like in the evolutionary loss of larval melanin pigmentation in the domesticated silkworm, Bombyx mori

The domesticated silkworm, Bombyx mori, and its wild progenitor, B. mandarina, are extensively studied as a model case of the evolutionary process of domestication. A conspicuous difference between these species is the dramatic reduction in melanin pigmentation in both larval and adult B. mori. Here we evaluate the efficiency of CRISPR/Cas9-targeted knockouts of pigment-related genes as a tool to understand their potential contributions to domestication-associated melanin pigmentation loss in B. mori. To demonstrate the efficacy of targeted knockouts in B. mandarina, we generated a homozygous CRISPR/Cas9-targeted knockout of yellow-y. In yellow-y knockout mutants, black body colour became lighter throughout the larval, pupal and adult stages, confirming a role for this gene in melanin pigment formation. Further, we performed allele-specific CRISPR/Cas9-targeted knockouts of the pigment-related transcription factor, apontic-like (apt-like) in B. mori × B. mandarina F1 hybrid individuals which exhibit B. mandarina-like larval pigmentation. Knockout of the B. mandarina allele of apt-like in F1 embryos results in white patches on the dorsal integument of larvae, whereas corresponding knockouts of the B. mori allele consistently exhibit normal F1 larval pigmentation. These results demonstrate a contribution of apt-like to the evolution of reduced melanin pigmentation in B. mori. Together, our results demonstrate the feasibility of CRISPR/Cas9-targeted knockouts as a tool for understanding the genetic basis of traits associated with B. mori domestication.

Bmmp influences wing morphology by regulating anterior-posterior and proximal-distal axes development

Insect wings are subject to strong selective pressure, resulting in the evolution of remarkably diverse wing morphologies that largely determine flight capacity. However, the genetic basis and regulatory mechanisms underlying wing size and shape development are not well understood. The silkworm Bombyx mori micropterous (mp) mutant exhibits shortened wing length and enlarged vein spacings, albeit without changes in total wing area. Thus, the mp mutant comprises a valuable genetic resource for studying wing development. In this study, we used molecular mapping to identify the gene responsible for the mp phenotype and designated it Bmmp. Phenotype-causing mutations were identified as indels and single nucleotide polymorphisms in noncoding regions. These mutations resulted in decreased Bmmp messenger RNA levels and changes in transcript isoform composition. Bmmp null mutants were generated by clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated protein 9 and exhibited changed wing shape, similar to mp mutants, and significantly smaller total wing area. By examining the expression of genes critical to wing development in wildtype and Bmmp null mutants, we found that Bmmp exerts its function by coordinately modulating anterior-posterior and proximal-distal axes development. We also studied a Drosophila mp mutant and found that Bmmp is functionally conserved in Drosophila. The Drosophila mp mutant strain exhibits curly wings of reduced size and a complete loss of flight capacity. Our results increase our understanding of the mechanisms underpinning insect wing development and reveal potential targets for pest control.

High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation

The silkworm Bombyx mori is 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.

The Bombyx mori Nucleopolyhedrovirus GP64 Retains the Transmembrane Helix of Signal Peptide to Contribute to Secretion across the Cytomembrane

Bombyx mori nucleopolyhedrovirus (BmNPV) is the primary pathogen of silkworms that causes severe economic losses in sericulture. GP64 is the key membrane fusion protein that mediates budded virus (BV) fusion with the host cell membrane. Previously, we found that the n-region of the GP64 signal peptide (SP) is required for protein secretion and viral pathogenicity; however, our understanding of BmNPV GP64 remains limited. Here, we first reported that BmNPV GP64 retained its SP in the mature protein and virion in only host cells but did not retain in nonhost cells. Uncleaved SP mediates protein targeting to the cytomembrane or secretion in Bombyx mori cells. The exitance of the n-region extended the transmembrane helix length, which resulted in the cleavage site to be located in the helix structure and thus blocked cleavage from signal peptidase (SPase). Without the n-region, the protein fails to be transported to the cytomembrane, but this failure can be rescued by the cleavage site mutation of SP. Helix-breaking mutations in SP abolished protein targeting to the cytomembrane and secretion. Our results revealed a previously unrecognized mechanism by which SP of membrane fusion not only determines protein localization but also determines viral pathogenicity, which highlights the escape mechanism of SP from the cleavage by SPase. IMPORTANCE BmNPV is the primary pathogen of silkworms, which causes severe economic losses in sericulture. BmNPV and Autographa californica multiple nucleopolyhedrovirus (AcMNPV) are closely related group I alphabaculoviruses, but they exhibit nonoverlapping host specificity. Recent studies suppose that GP64 is a determinant of host range, while knowledge remains limited. In this study, we revealed that BmNPV GP64 retained its SP in host cells but not in nonhost cells, and the SP retention is required for GP64 secretion across the cytomembrane. This is the first report that a type I membrane fusion protein retained its SP in mature proteins and virions. Our results unveil the mechanism by which SP GP64 escapes cleavage and the role of SP in protein targeting. This study will help elucidate an important mechanistic understanding of BmNPV infection and host range specificity.

Use of Bombyx mori silk fibroin in tissue engineering: From cocoons to medical devices, challenges, and future perspectives

Silk fibroin has become a prominent material in tissue engineering (TE) over the last 20 years with almost 10,000 published works spanning in all the TE applications, from skeleton to neuronal regeneration. Fibroin is an extremely versatile biopolymer that, due to its ease of processing, has enabled the development of an entire plethora of materials whose properties and architectures can be tailored to suit target applications. Although the research and development of fibroin TE materials and devices is mature, apart from sutures, only a few medical products made of fibroin are used in the clinical routines. <40 clinical trials of Bombyx mori silk-related products have been reported by the FDA and few of them resulted in a commercialized device. In this review, after explaining the structure and properties of silk fibroin, we provide an overview of both fibroin constructs existing in the literature and fibroin devices used in clinic. Through the comparison of these two categories, we identified the burning issues faced by fibroin products during their translation to the market. Two main aspects will be considered. The first is the standardization of production processes, which leads both to the standardization of the characteristics of the issued device and the correct assessment of its failure. The second is the FDA regulations, which allow new devices to be marketed through the 510(k) clearance by demonstrating their equivalence to a commercialized medical product. The history of some fibroin medical devices will be taken as a case study. Finally, we will outline a roadmap outlining what actions we believe are needed to bring fibroin products to the market.

Bombyx mori Vps13d is a key gene affecting silk yield

Bombyx mori is an important economic insect, its economic value mainly reflected in the silk yield. The major functional genes affecting the silk yield of B. mori have not been determined yet. Bombyx mori vacuolar protein sorting-associated protein 13d (BmVps13d) has been identified, but its function is not reported. In this study, BmVps13d protein shared 30.84% and 34.35% identity with that of in Drosophila melanogaster and Homo. sapiens, respectively. The expressions of BmVps13d were significantly higher in the midgut and silk gland of JS (high silk yield) than in that of L10 (low silk yield). An insertion of 9 bp nucleotides and two deficiencies of adenine ribonucleotides in the putative promoter region of BmVps13d gene in L10 resulted in the decline of promoter activity was confirmed using dual luciferase assay. Finally, the functions of BmVps13d in B. mori were studied using the CRISPR/Cas9 system, and the mutation of BmVps13d resulted in a 24.7% decline in weight of larvae, as well as a 27.1% (female) decline and a 11.8% (male) decline in the silk yield. This study provides a foundation for studying the molecular mechanism of silk yield and breeding the silkworm with high silk yield.

Molecular mechanism and potential application of bacterial infection in the silkworm, Bombyx mori

As a representative species of Lepidoptera, Bombyx mori has been widely studied and applied. However, bacterial infection has always been an important pathogen threatening the growth of silkworms. Bombyx mori can resist various pathogenic bacteria through their own physical barrier and innate immune system. However, compared with other insects, such as Drosophila melanogaster, research on the antibacterial mechanism of silkworms is still in its infancy. This review systematically summarized the routes of bacterial infection in silkworms, the antibacterial mechanism of silkworms after ingestion or wounding infection, and the intestinal bacteria and infection of silkworms. Finally, we will discuss silkworms as a model animal for studying bacterial infectious diseases and screening antibacterial drugs.

Identification and Characterization of Genes Related to Resistance of Autographa californica Nucleopolyhedrovirus Infection in Bombyx mori

Simple Summary

Autographa californica nucleopolyhedrovirus (AcMNPV) is a kind of baculovirus that was initially found and named for its host, but the previous study reveals several silkworm strains are preferentially susceptible to AcMNPV through intrahemocelical injection method. In the following study, genetics analysis showed that a set of potential genes which controlled resistance of AcMNPV was located on chromosome 3. In the present research, we performed Genome-Wide Association Studies to identify the gene that controls the resistance of AcMNPV, results show that the Niemann-Pick C1 (NPC-1) gene is strongly associated with this resistance. Then we found that there are several amino acid mutations in the protein sequence of BmNPC1 between two different resistance strains of Bombyx mori. RNAi results showed that BmNPC1 successfully suppressed virus infection ability and changed the expression pattern of viral genes.

Abstract

In Bombyx mori, as an important economic insect, it was first found that some strains were completely refractory to infection with Autographa californica nucleopolyhedrovirus (AcMNPV) through intrahemocelical injection; whereas almost all natural strains had difficulty resisting Bombyx mori nucleopolyhedrovirus (BmNPV), which is also a member of the family Baculoviridae. Previous genetics analysis research found that this trait was controlled by a potentially corresponding locus on chromosome 3, but the specific gene and mechanism was still unknown. With the help of the massive silkworm strain re-sequencing dataset, we performed the Genome-Wide Association Studies (GWAS) to identify the gene related to the resistance of AcMNPV in this study. The GWAS results showed that the Niemann-Pick type C1 (NPC-1) gene was the most associated with the trait. The knockdown experiments in BmN cells showed that BmNPC1 has a successful virus suppression infection ability. We found a small number of amino acid mutations among different resistant silkworms, which indicates that these mutations contributed to the resistance of AcMNPV. Furthermore, inhibition of the BmNPC1 gene also changed the viral gene expression of the AcMNPV, which is similar to the expression profile in the transcriptome data of p50 and C108 strains.

Whole-genome sequences of 37 breeding line Bombyx mori strains and their phenotypes established since 1960s

Bombyx mori is a key insect in the sericulture industry and one of the very important economic animals that are responsible for not only the livelihood of many farmers internationally but also expended biomedical use. The National Institute of Agricultural Sciences of the Rural Development Administration of Korea (NIAS, RDA, Korea) has been collecting silkworm resources with various phenotypic traits from the 1960s and established breeding lines for using them as genetic resources. And these breeding line strains have been used to develop suitable F1 hybrid strains for specific use. In this study, we report the whole-genome sequences of 37 breeding line B. mori strains established over the past 60 years, along with the description of their phenotypic characteristics with photos of developmental stages. In addition, we report the example phenotypic characteristics of the F1-hybrid strain using these breeding line strains. We hope this data will be used as valuable resources to the related research community for studying B. mori and similar other insects.

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.

Mitochondrial Genome Architecture and Evolutionary Origin of the Yao Silkworm, a Living Fossil of the Domestic Silkworm Bombyx mori (Lepidoptera: Bombycidae)

The Yao silkworm is a unique silkworm resource producing yellow flat plate silk that has only been reared by the Baiku Yao ethnic group in Nandan County, Guangxi Province, China for a thousand years. Here, we report the mitochondrial genomes (mitogenomes) of five Yao silkworm strains and 10 local Guangxi strains of the domestic silkworm (Bombyx mori) L. (Lepidoptera: Bombycidae), and use the resulting mitogenomes and the available Bombyx mitogenomes to characterize their genome architecture and trace the evolutionary origin of the Yao silkworm. The five Yao silkworm mitogenomes exhibited genome architectures identical to typical set of 37 mitochondrial genes (13 protein-coding genes, 22 transfer RNAs, and two ribosomal RNAs) and a high level of genome sequence similarity with the domestic silkworm. Mitogenome-based phylogenetic reconstruction provided solid evidence that the Yao silkworm shares a common ancestor with the domestic silkworm. Sliding window analysis uncovered a distinct variation pattern in the mitogenome between the Yao silkworm and the other domestic silkworm strains. The phylogenetic analyses revealed a basal placement of the Yao silkworm among all available domestic silkworm strains, indicating that the Yao silkworm is an ancient population of the domestic silkworm. Our data indicated that the Yao silkworm (B. mori) is a lineage of the domestic silkworm, which for the first time provides insights into the origin of the Yao silkworm.

Phylogeographic Relationships among Bombyx mandarina (Lepidoptera: Bombycidae) Populations and Their Relationships to B. mori Inferred from Mitochondrial Genomes

We report 37 mitochondrial genome (mitogenome) sequences of Bombyx mori strains (Lepidoptera: Bombycidae) and four of B. mandarina individuals, each preserved and collected, respectively, in South Korea. These mitogenome sequences combined with 45 public data showed a substantial genetic reduction in B. mori strains compared to the presumed ancestor B. mandarina, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, concordant to the Qinling-Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses, as has previously been detected. However, one individual that was in close proximity to the south Qinling-Huaihe line was exceptional, belonging to the northern group. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. Although a substantial genetic divergence is present between B. mandarina in southern China and Japan, a highly supported sister relationship between the two regional populations may suggest the potential origin of Japanese B. mandarina from southern China instead of the Korean peninsula.

Involvement of the Clock Gene Period in the Photoperiodism of the Silkmoth Bombyx mori

We established a knockout strain of a clock gene, period (per), by using TALEN in a bivoltine strain (Kosetsu) of Bombyx mori (Insecta, Lepidoptera), and examined the effect of per knockout on the circadian rhythm and photoperiodism. The generated per knockout allele was considered to be null, because a new stop codon was present in the insertion allele. The wild type (Kosetsu) showed clear circadian rhythms in eclosion and hatching, whereas the per knockout strain showed arrhythmic eclosion and hatching under constant darkness. In this strain, moreover, temporal expression changes of clock genes per and timeless were disrupted. The wild type showed a clear long-day response for induction of embryonic diapause: when larvae were reared under long-day and short-day conditions at 25°C, adults produced nondiapause and diapause eggs, respectively. However, the per knockout strain lost the sensitivity to photoperiod and laid nondiapause eggs under both conditions. We conclude that per plays an important role both in circadian rhythms and in photoperiodism of B. mori, indicating the involvement of the circadian clock consisting of per in the photoperiodism.

Transgenic genome editing-derived antiviral therapy to nucleopolyhedrovirus infection in the industrial strain of the silkworm

The silkworm (Bombyx mori) is a domesticated and economically important insect. During the whole growth period, silkworm suffers various pathogen infection. Nearly 80% of silk cocoon crop losses are attributed to viral diseases. The circular double-stranded DNA virus Bombyx mori nuclepolyhedrovirus (BmNPV) is the major viral pathogen responsible for massive silkworm death and huge economic losses in the sericulture industry. Up to now, almost all the new strategies for developing immunity against BmNPV are in laboratory strains because of the lack of transgenic technology in industrial silkworm strains. We previously demonstrated that modification of BmNPV genome DNA with the antivirus transgenic CRISPR/Cas9 system effectively improved the resistance of laboratory silkworm strains to viral pathogens. The industrial strains are monovoltine or bivoltine. It is very difficult to break the diapause before microinjection for genetic transformation. Here, we show that the anti-BmNPV transgenic CRISPR/Cas9 system also works in the industrial silkworm strain Jingsong. In this study, we successfully broke diapause and obtained generation-skipping embryos and constructed two TG Jingsong lines. Both lines exhibited significantly enhanced immunity to BmNPV without significant changes in most of the commercially important traits. These results demonstrate that the construction of TG silkworm lines carrying a heritable immune defense system against BmNPV could be an effective strategy to enhance the resistance of industrial silkworm strains to the most devastating DNA virus. The research opened a window for genetic modification of the important strains from laboratory strains to industrial strains.

Drivers of insect consumption across human populations

Discussions regarding entomophagy in humans have been typically led by entomologists. While anthropologists devote much time to understanding diverse human subsistence practices, historical and cultural variation in insect consumption remains largely unexplained. This review explores the relation between variable ecologies, subsistence strategies, and social norms on insect consumption patterns across past and contemporary human populations. Ecological factors, such as the nutritional contribution of edible insects relative to those of other foraged or farmed resources available, may help explain variation in their consumption. Additionally, our evolved social learning strategies may help propagate social norms that prohibit or prioritize the consumption of some or all edible insects, independent of their profitability. By adopting a behavioral ecological and cultural evolutionary approach, this review aims to resolve current debates on insect consumption and provide directions for future research.

High-Quality de novo Chromosome-Level Genome Assembly of a Single Bombyx mori With BmNPV Resistance by a Combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing

The reference genomes of Bombyx mori (B. mori), Silkworm Knowledge-based database (SilkDB) and SilkBase, have served as the gold standard for nearly two decades. Their use has fundamentally shaped model organisms and accelerated relevant studies on lepidoptera. However, the current reference genomes of B. mori do not accurately represent the full set of genes for any single strain. As new genome-wide sequencing technologies have emerged and the cost of high-throughput sequencing technology has fallen, it is now possible for standard laboratories to perform full-genome assembly for specific strains. Here we present a high-quality de novo chromosome-level genome assembly of a single B. mori with nuclear polyhedrosis virus (BmNPV) resistance through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. In addition, regular bioinformatics analyses, such as gene family, phylogenetic, and divergence analyses, were performed. The sample was from our unique B. mori species (NB), which has strong inborn resistance to BmNPV. Our genome assembly showed good collinearity with SilkDB and SilkBase and particular regions. To the best of our knowledge, this is the first genome assembly with BmNPV resistance, which should be a more accurate insect model for resistance studies.

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.

Role of juvenile hormone receptor Methoprene-tolerant 1 in silkworm larval brain development and domestication

The insect brain is the central part of the neurosecretory system, which controls morphology, physiology, and behavior during the insect's lifecycle. Lepidoptera are holometabolous insects, and their brains develop during the larval period and metamorphosis into the adult form. As the only fully domesticated insect, the Lepidoptera silkworm Bombyx mori experienced changes in larval brain morphology and certain behaviors during the domestication process. Hormonal regulation in insects is a key factor in multiple processes. However, how juvenile hormone (JH) signals regulate brain development in Lepidoptera species, especially in the larval stage, remains elusive. We recently identified the JH receptor Methoprene tolerant 1 ( Met1) as a putative domestication gene. How artificial selection on Met1 impacts brain and behavioral domestication is another important issue addressing Darwin's theory on domestication. Here, CRISPR/Cas9-mediated knockout of Bombyx Met1 caused developmental retardation in the brain, unlike precocious pupation of the cuticle. At the whole transcriptome level, the ecdysteroid (20-hydroxyecdysone, 20E) signaling and downstream pathways were overactivated in the mutant cuticle but not in the brain. Pathways related to cell proliferation and specialization processes, such as extracellular matrix (ECM)-receptor interaction and tyrosine metabolism pathways, were suppressed in the brain. Molecular evolutionary analysis and in vitro assay identified an amino acid replacement located in a novel motif under positive selection in B. mori, which decreased transcriptional binding activity. The B. mori MET1 protein showed a changed structure and dynamic features, as well as a weakened co-expression gene network, compared with B. mandarina. Based on comparative transcriptomic analyses, we proposed a pathway downstream of JH signaling (i.e., tyrosine metabolism pathway) that likely contributed to silkworm larval brain development and domestication and highlighted the importance of the biogenic amine system in larval evolution during silkworm domestication.

Large-scale whole-genome resequencing unravels the domestication history of Cannabis sativa

Cannabis sativa has long been an important source of fiber extracted from hemp and both medicinal and recreational drugs based on cannabinoid compounds. Here, we investigated its poorly known domestication history using whole-genome resequencing of 110 accessions from worldwide origins. We show that C. sativa was first domesticated in early Neolithic times in East Asia and that all current hemp and drug cultivars diverged from an ancestral gene pool currently represented by feral plants and landraces in China. We identified candidate genes associated with traits differentiating hemp and drug cultivars, including branching pattern and cellulose/lignin biosynthesis. We also found evidence for loss of function of genes involved in the synthesis of the two major biochemically competing cannabinoids during selection for increased fiber production or psychoactive properties. Our results provide a unique global view of the domestication of C. sativa and offer valuable genomic resources for ongoing functional and molecular breeding research.

The global significance of biodiversity science in China: an overview

Biodiversity science in China has seen rapid growth over recent decades, ranging from baseline biodiversity studies to understanding the processes behind evolution across dynamic regions such as the Qinghai-Tibetan Plateau. We review research, including species catalogues; biodiversity monitoring; the origins, distributions, maintenance and threats to biodiversity; biodiversity-related ecosystem function and services; and species and ecosystems' responses to global change. Next, we identify priority topics and offer suggestions and priorities for future biodiversity research in China. These priorities include (i) the ecology and biogeography of the Qinghai-Tibetan Plateau and surrounding mountains, and that of subtropical and tropical forests across China; (ii) marine and inland aquatic biodiversity; and (iii) effective conservation and management to identify and maintain synergies between biodiversity and socio-economic development to fulfil China's vision for becoming an ecological civilization. In addition, we propose three future strategies: (i) translate advanced biodiversity science into practice for biodiversity conservation; (ii) strengthen capacity building and application of advanced technologies, including high-throughput sequencing, genomics and remote sensing; and (iii) strengthen and expand international collaborations. Based on the recent rapid progress of biodiversity research, China is well positioned to become a global leader in biodiversity research in the near future.

Low levels of genetic differentiation with isolation by geography and environment in populations of Drosophila melanogaster from across China

The fruit fly, Drosophila melanogaster, is a model species in evolutionary studies. However, population processes of this species in East Asia are poorly studied. Here we examined the population genetic structure of D. melanogaster across China. There were 14 mitochondrial haplotypes with 10 unique ones out of 23 known from around the globe. Pairwise F_ST values estimated from 15 novel microsatellites ranged from 0 to 0.11, with geographically isolated populations showing the highest level of genetic uniqueness. STRUCTURE analysis identified high levels of admixture at both the individual and population levels. Mantel tests indicated a strong association between genetic distance and geographical distance as well as environmental distance. Full redundancy analysis (RDA) showed that independent effects of environmental conditions and geography accounted for 62.10% and 31.58% of the total explained genetic variance, respectively. When geographic variables were constrained in a partial RDA analysis, the environmental variables bio2 (mean diurnal air temperature range), bio13 (precipitation of the wettest month), and bio15 (precipitation seasonality) were correlated with genetic distance. Our study suggests that demographic history, geographical isolation, and environmental factors have together shaped the population genetic structure of D. melanogaster after its introduction into China.

Proteomic Identification of Bombyx mori Organelles Using the Engineered Ascorbate Peroxidase APEX and Development of Silkworm Organelle Proteome Database (SilkOrganPDB)

Silkworm Bombyx mori is an economically important insect and a lepidopteran model. Organelle proteome is vital to understanding gene functions; however, it remains to be identified in silkworm. Here, using the engineered ascorbate peroxidase APEX, we constructed transgenic B. mori embryo cells (BmE) expressing APEX-NLS, COX4-APEX, APEX-Rev, and APEX-KDEL in nucleus, mitochondrial matrix (MM), cytosol, and endoplasmic reticulum (ER), and isolated the biotin-labeled proteins using streptavidin-affinity purification, respectively. The isolated proteins were determined using LC-MS/MS and annotated by searching B. mori genomes downloaded from GenBank, SilkBase, SilkDB 2.0, and SilkDB 3.0, resulting in 842, 495, 311, and 445 organelle proteins identified, respectively. We mapped the 296 MM proteins annotated in the GenBank data to mitochondrial protein databases of the fly, human, and mouse, and found that 140 (47%) proteins are homologous to 80 fly proteins, and 65 (22%) proteins match to 31 and 29 human and mouse proteins, respectively. Protein orthology was predicted in multiple insects using OrthoMCL, producing 460 families containing 839 proteins we identified. Out of 460 families, 363 were highly conserved and found in all insects, leaving only three proteins without orthology in other insects, indicating that the identified proteins are highly conserved and probably play important roles in insects. A gene ontology enrichment analysis by clusterProfiler revealed that the nucleus proteins significantly enriched in cellular component terms of nucleus and nucleolus, the MM proteins markedly enriched in molecular function terms of nucleotide binding, and the cytosol proteins mainly enriched in biological process terms of small molecule metabolism. To facilitate the usage and analysis of our data, we developed an open-access database, Silkworm Organelle Proteome Database (SilkOrganPDB), which provides multiple modules for searching, browsing, downloading, and analyzing these proteins, including BLAST, HMMER, Organelle Proteins, Protein Locations, Sequences, Gene Ontology, Homologs, and Phylogeny. In summary, our work revealed the protein composition of silkworm BmE organelles and provided a database resource helpful for understanding the functions and evolution of these proteins.

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.

The Number of Larval Molts Is Controlled by Hox in Caterpillars

Animals with exoskeletons molt for further growth. In insects, the number of larval (or nymphal) molts varies inter- and intra-specifically, and it is widely accepted that the variation in the number of larval molts is an adaptive response to diverse environmental conditions.^1-5 However, the molecular mechanism that underlies the variety and plasticity in the number of larval molts is largely unknown. In the silkworm, Bombyx mori, there are strains that molt three, four, or five times, and these numbers are determined by allelic variation at a single autosomal locus, Moltinism (M).^6-9 Here, we demonstrate that the Hox gene Sex combs reduced (Scr) is responsible for the phenotypes of the M locus. Scr is selectively expressed in the larval prothoracic gland (PG), an endocrine organ that produces molting hormones.²Scr represses the biosynthesis of molting hormones in the PG, thereby regulating the incremental increase in body size during each larval instar. Our experiments consistently suggest that the differential expression levels of Scr among the three M alleles result in different growth ratios that ultimately lead to the different number of larval molts. Although the role of Hox genes in conferring segmental identity along the body axis and in molding segment-specific structure later in development has been well established,^10-13 the present study identifies an unexpected role of Hox gene in hormone biosynthesis. This new role means that, in addition to shaping segment-specific morphology, Hox genes also drive the evolution of life history traits by regulating animal physiology.

Cytochrome P450 monooxygenase genes in the wild silkworm, Bombyx mandarina

Wild (Bombyx mandarina) and domestic silkworms (B. mori) are good models for investigating insect domestication, as 5000 years of artificial breeding and selection have resulted in significant differences between B. mandarina and B. mori. In this study, we improved the genome assemblies to the chromosome level and updated the protein-coding gene annotations for B. mandarina. Based on this updated genome, we identified 68 cytochrome P450 genes in B. mandarina. The cytochrome P450 repository in B. mandarina is smaller than in B. mori. Certain currently unknown key genes, rather than gene number, are critical for insecticide resistance in B. mandarina, which shows greater resistance to insecticides than B. mori. Based on the physical maps of B. mandarina, we located 66 cytochrome P450s on 18 different chromosomes, and 27 of the cytochrome P450 genes were concentrated into seven clusters. KEGG enrichment analysis of the P450 genes revealed the involvement of cytochrome P450 genes in hormone biosynthesis. Analyses of the silk gland transcriptome identified candidate cytochrome P450 genes (CYP306A) involved in ecdysteroidogenesis and insecticide metabolism in B. mandarina.

Genetic Mapping of Climbing and Mimicry: Two Behavioral Traits Degraded During Silkworm Domestication

Behavioral changes caused by domestication in animals are an important issue in evolutionary biology. The silkworm, Bombyx mori, is an ideal fully domesticated insect model for studying both convergent domestication and behavior evolution. We explored the genetic basis of climbing for foraging and mimicry, two degraded behaviors during silkworm domestication, in combination of bulked segregant analysis (BSA) and selection sweep screening. One candidate gene, ASNA1, located in the 3-5 Mb on chromosome 19, harboring a specific non-synonymous mutation in domestic silkworm, might be involved in climbing ability. This mutation was under positive selection in Lepidoptera, strongly suggesting its potential function in silkworm domestication. Nine candidate domesticated genes related to mimicry were identified on chromosomes 13, 21, and 27. Most of the candidate domesticated genes were generally expressed at higher levels in the brain of the wild silkworm. This study provides valuable information for deciphering the molecular basis of behavioral changes associated with silkworm domestication.

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.