The evolutionary road from wild moth to domestic silkworm

Dimm targets GDAP2 to regulate larval development in the silkworm, Bombyx mori

The basic helix-loop-helix (bHLH) domain transcription factors precisely regulate various developmental processes in insects. Dimm, a specific bHLH transcription factor, integrates the insulin/insulin-like growth factor signaling (IIS) and juvenile hormone signaling (JHS) pathways to modulate larval development in silkworms. However, the molecular mechanisms underlying this regulation are not yet fully understood. This study aimed to determine the targets of Dimm through which it regulates larval development. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed ganglioside-induced differentiation-associated protein 2 (GDAP2) as a direct downstream target gene of Dimm. Further study showed that Dimm directly binds to an enhancer element located in the second intron of the GDAP2 gene to promote its transcription. GDAP2 exhibited widespread expression across different stages and tissues of silkworms, regulated by both the IIS and the JHS pathways. The systemic knockout of GDAP2 leads to delayed larval development with a significant reduction in body weight; moreover, larval development was arrested at the 4th-instar stage. Further investigation unveiled that the inhibition of the ecdysone and innate immune signaling pathways in the mutant line led to abnormal larval development. A systematic investigation of the biological functions of GDAP2 offers valuable insights into the mechanism by which Dimm integrates IIS and JHS pathways to regulate the larval development of silkworms.

Adaptive genomic signatures of globally invasive populations of the yellow fever mosquito Aedes aegypti

In the arboviral vector Aedes aegypti, adaptation to anthropogenic environments has led to a major evolutionary shift separating the domestic Aedes aegypti aegypti (Aaa) ecotype from the wild Aedes aegypti formosus (Aaf) ecotype. Aaa mosquitoes are distributed globally and have higher vectorial capacity than Aaf, which remained in Africa. Despite the evolutionary and epidemiological relevance of this separation, inconsistent morphological data and a complex population structure have hindered the identification of genomic signals distinguishing the two ecotypes. Here we assessed the correspondence between the geographic distribution, population structure and genome-wide selection of 511 Aaf and 123 Aaa specimens and report adaptive signals in 186 genes that we call Aaa molecular signatures. Our results indicate that Aaa molecular signatures arose from standing variation associated with extensive ancestral polymorphisms in Aaf populations and have been co-opted for self-domestication through genomic and functional redundancy and local adaptation. Overall, we show that the behavioural shift of Ae. aegypti mosquitoes to live in association with humans relied on the fine regulation of chemosensory, neuronal and metabolic functions, as seen in the domestication processes of rabbits and silkworms. Our results also provide a foundation for the investigation of new genic targets for the control of Ae. aegypti populations.

Insufficient wing development possibly contributes to flightlessness of the silkworm Bombyx mori during domestication

The silkworm Bombyx mori is a fully domesticated insect and has lost its flight capability during domestication. However, it does not show obvious degradation of wing, rendering the mechanism of flightlessness elusive. In this study, we demonstrated that B. mori showed insufficient wing disc developmental properties, weaker adult wings and loosely arranged flight muscles with smaller cell size and relatively weaker wing flapping frequency, compared with Bombyx mandarina. Consistently, switch of transcriptomic landscape from larval to pupal wing discs occurred earlier in B. mori, with an earlier repression of genes related to development of wing and flight muscle precursor cells, as well as of those genes functioning in chromosome remodelling during metamorphosis. Knockout of the wing-disc expressing and the well-documented flight muscle gene fln caused loosely arranged flight muscles and reduced flight capacity. Meanwhile, core mitochondrial genes CoxI and CoxIII were downregulated in B. mori, accompanied by upregulation of CoxII and CoxIV and activated adverse cellular constituents such as reactive oxygen species, suggesting a possible functional impairment in the mitochondria of B. mori wing disc cells compared with those of B. mandarina. The results indicate that the flight system of B. mori has been structurally and functionally weakened by domestication and provide new insights into understanding behavioural domestication and insect flight from the developmental view.

A chromosome-level genome assembly of wild silkmoth, Bombyx mandarina

The wild silk moth, Bombyx mandarina, is the closest relative of the domesticated silk moth, Bombyx mori. National BioResource Project of Japan (NBRP) maintains a B. mandarina strain derived from individuals captured at Sakado (Saitama, Japan) in 1982. Now, NBRP has developed chromosome replacement strains of B. mori. In each strain, one autosome of B. mori was replaced with the corresponding chromosome of B. mandarina. To facilitate the use of chromosome replacement strains and B. mandarina itself, we constructed a chromosome-level genome assembly of B. mandarina. Furthermore, we performed functional annotations of the genome assembly, i.e., transcriptome-based gene prediction, Assay for Transposase-Accessible Chromatin (ATAC)-seq, and PIWI-interacting RNA (piRNA)-targeted small RNA-seq. The assembly harbors 14,859 protein-coding genes and 628 piRNA clusters across three tissues: ovaries, testis, and embryos. ATAC-seq data comprehensively detected open chromosome regions, which will benefit when CRISPR/Cas9-mediated genome editing is conducted.

The T2T Genome of the Domesticated Silkworm Bombyx mori

Genome sequences contain the fundamental genetic information that largely determines the biology of a species. Over the past 20 years, advancements in high-throughput sequencing technologies and bioinformatics tools have matured, facilitating genome assembly and ushering in the telomere-to-telomere (T2T) era. Bombyx mori is renowned as a silk-producing insect and serves as an important model organism extensively studied across various fields of biology. In this study, we present the first assembled T2T genome by integrating HiFi, ultra-long ONT, NGS, and Hi-C data. This assembly comprises 450,267,439 base pairs from 28 chromosomes and includes annotations for a total of 18,253 protein-coding genes. A completeness evaluation revealed that 99.1% of conserved single-copy genes were included, as determined by a BUSCO analysis. Furthermore, the consensus quality (QV) assessed through Merqury was recorded at 59.88. The proportion of repeat sequence achieved 60.77%, marking it as the highest reported value for B. mori to date. In comparison to previously published genomes, our assembly offers a more complete and higher quality representation, particularly concerning highly homologous tandem regions such as telomeres, rDNA clusters, and Gr family regions. Furthermore, our extensive experience in genome assembly, including sample preparation experience and assembly strategies to reduce complexity, will provide valuable references for other species aiming to achieve their own T2T genome assemblies.

Sharp decline in male fertility in F2 hybrids of the female-heterogametic silk moth Bombyx

Sexual selection drives rapid evolution of morphological, physiological, and behavioral traits, especially in males, and it may also drive the rapid evolution of hybrid male sterility. Indeed, the faster male theory of speciation was once viewed as a major cause of Haldane's rule in male-heterogametic XY taxa, but is increasingly being replaced by the genetic conflict hypothesis partly because it cannot explain the faster evolution of hybrid female sterility in female-heterogametic ZW taxa. The theory nonetheless predicts that there should be more genes for hybrid male sterility than for hybrid female sterility even in such taxa, but this remains untested. Thus, finding evidence for the faster male theory of reproductive isolation beyond the F1 generation in ZW systems still represents a challenge to studying the impact of sexual selection. In this study, we examined F2 hybrids between the domesticated silkworm Bombyx mori and the wild silk moth Bombyx mandarina, which have ZW sex determination. We found that although only females showed reduced fertility in the F1 generation, the F2 hybrid males had a significant reduction in fertility compared with the parental and F1 males. Importantly, 27% of the F2 males and 15% of the F2 females were completely sterile, suggesting the presence of recessive incompatibilities causing male sterility in female-heterogametic taxa.

PI3K-Akt-SGF1-Dimm pathway mediates the nutritional regulation of silk protein synthesis in Bombyx mori

The efficient synthesis of silk protein is heavily reliant on the ingestion of massive nutrients during the peak growth phase in the silkworm. However, the molecular mechanism of nutritional regulation of silk protein synthesis remains unknown. In this study, we investigated the impact of nutrient deficiency on the synthesis of silk protein. Nutritional deficiency led to a reduction in silk yield, accompanied by decreased levels of silk proteins and fibroin heavy chain (FibH)-activating transcription factors SGF1 and Dimm. Furthermore, insulin enhanced the protein levels of SGF1 and Dimm, which can be attenuated by specific inhibitors of PI3K. Co-immunoprecipitation analysis showed that the nutrient pathway factor protein kinase B (Akt) could interact with SGF1 protein. Knockdown of Akt reduced the phosphorylation level of SGF1 and impedes its nuclear translocation. Further studies revealed that SGF1 was directly bound to Fkh site in the 22-43 region upstream of ATG of Dimm gene to activate its transcription. In conclusion, during the peak growth phase, nutrition promotes the massive synthesis of silk protein through the PI3K-Akt-SGF1-Dimm pathway. This study offers valuable insights into the efficient synthesis of silk proteins and establishes a theoretical foundation for improving silk yield.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

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.

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.

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.

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.

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.