An update of KAIKObase, the silkworm genome database

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.

Dead ringer acts as a major regulator of juvenile hormone biosynthesis in insects

In holometabolous insects, proper control of the production of juvenile hormone (JH), which maintains larval traits, is crucial for successful metamorphosis. JH is produced specifically in the corpora allata (CA) via the functioning of a set of JH biosynthetic enzymes (JHBEs). Expression of JHBE genes in the CA is coordinated except for JH acid methyltransferase (JHAMT), which functions in the last step of JH biosynthesis. Here, we sought to determine the mechanism that enables this coordinated expression, assuming the presence of a central regulator of JHBE genes. Comparison of transcriptomes in the CA during active and inactive stages revealed the presence of 3 transcription factors, whose expression patterns matched those of JHBE genes. We propose that one of these, Dead ringer (Dri), is the central up-regulator of CA-specific JHBE genes including JHAMT, based on the following findings: (ⅰ) Knockdown of Dri in the larvae caused precocious metamorphosis, which was rescued by the exogenous application of JH analog, and (ⅱ) knockdown of Dri decreased the expression of most CA-specific JHBE genes examined. Furthermore, RNAi-based reverse genetics indicated that Dri works most upstream in the control of CA-specific JHBE genes, and that shutdown of JHAMT, which occurs independent of other JHBE genes prior to the onset of metamorphosis, can be hypothetically explained by the presence of an unidentified repressor. Our study suggests that Dri, which has been known to regulate embryonic development in a wide range of animals, is conferred a new role in holometabolous insects, i.e. central regulation of CA-specific JHBE genes.

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.

Time-course transcriptome data of silk glands in day 0-7 last-instar larvae of Bombyx mori (w1 pnd strain)

Time-course transcriptome expression data were constructed for four parts of the silk gland (anterior, middle, and posterior parts of the middle silk gland, along with the posterior silk gland) in the domestic silkworm, Bombyx mori, from days 0 to 7 of the last-instar larvae. For sample preparation, silk glands were extracted from one female and one male larva every 24 hours accurately after the fourth ecdysis. The reliability of these transcriptome data was confirmed by comparing the transcripts per million (TPM) values of the silk gene and quantitative reverse transcription PCR results. Hierarchical cluster analysis results supported the reliability of transcriptome data. These data are likely to contribute to the progress in molecular biology and genetic research using B. mori, such as elucidating the mechanism underlying the massive production of silk proteins, conducting entomological research using a meta-analysis as a model for lepidopteran insect species, and exploring medical research using B. mori as a model for disease species by utilising transcriptome data.

Effect of acute exposure of Hg on physiological parameters and transcriptome expression in silkworms (Bombyx mori)

Mercury (Hg) contamination poses a global threat to the environment, given its elevated ecotoxicity. Herein, we employed the lepidopteran model insect, silkworm (Bombyx mori), to systematically investigate the toxic effects of Hg-stress across its growth and development, histomorphology, antioxidant enzyme activities, and transcriptome responses. High doses of Hg exposure induced evident poisoning symptoms, markedly impeding the growth of silkworm larvae and escalating mortality in a dose-dependent manner. Under Hg exposure, the histomorphology of both the midgut and fat body exhibited impairments. Carboxylesterase (CarE) activity was increased in both midgut and fat body tissues responding to Hg treatment. Conversely, glutathione S-transferase (GST) levels increased in the fat body but decreased in the midgut. The transcriptomic analysis revealed that the response induced by Hg stress involved multiple metabolism processes. Significantly differently expressed genes (DEGs) exhibited strong associations with oxidative phosphorylation, nutrient metabolisms, insect hormone biosynthesis, lysosome, ribosome biogenesis in eukaryotes, and ribosome pathways in the midgut or the fat body. The findings implied that exposure to Hg might induce the oxidative stress response, attempting to compensate for impaired metabolism. Concurrently, disruptions in nutrient metabolism and insect hormone activity might hinder growth and development, leading to immune dysfunction in silkworms. These insights significantly advance our theoretical understanding of the potential mechanisms underlying Hg toxicity in invertebrate organisms.

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.

Fine mapping of Green a, Ga, on chromosome 27 in Bombyx mori

Some strains of silkworms produce green cocoons of varying intensities. This results from quantitative and qualitative differences in flavonoid pigments, which are influenced by the environment and genetic background. We discovered that the appearance of a faint green cocoon is regulated by a gene (G27) located on chromosome 27. Through mating experiments, we found that G27 is identical to an essential flavonoid cocoon gene, Ga. This locus has not been previously described. Furthermore, we narrowed down the Ga region to 438 kbp using molecular markers. Within this region, several predicted genes for sugar transporters form a cluster structure, suggesting that Ga is among them.

Structure of a superoxide dismutase from a tardigrade: Ramazzottius varieornatus strain YOKOZUNA-1

Superoxide dismutase (SOD) is an essential and ubiquitous antioxidant protein that is widely present in biological systems. The anhydrobiotic tardigrades are some of the toughest micro-animals. They have an expanded set of genes for antioxidant proteins such as SODs. These proteins are thought to play an essential role in oxidative stress resistance in critical situations such as desiccation, although their functions at the molecular level have yet to be explored. Here, crystal structures of a copper/zinc-containing SOD (RvSOD15) from an anhydrobiotic tardigrade, Ramazzottius varieornatus strain YOKOZUNA-1, are reported. In RvSOD15, one of the histidine ligands of the catalytic copper center is replaced by a valine (Val87). The crystal structures of the wild type and the V87H mutant show that even though a histidine is placed at position 87, a nearby flexible loop can destabilize the coordination of His87 to the Cu atom. Model structures of other RvSODs were investigated and it was found that some of them are also unusual SODs, with features such as deletion of the electrostatic loop or β3 sheet and unusual metal-binding residues. These studies show that RvSOD15 and some other RvSODs may have evolved to lose the SOD function, suggesting that gene duplications of antioxidant proteins do not solely explain the high stress tolerance of anhydrobiotic tardigrades.

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.

Overexpression of Ultrabithorax Changes the Development of Silk Gland and the Expression of Fibroin Genes in Bombyx mori

Ultrabithorax (Ubx) is a member of the Hox gene group involved in cell fate decisions, cell proliferation and organ identity. Its function has been extensively researched in Drosophila melanogaster but little is known about it in Lepidoptera. To uncover the function of Ubx in the development of lepidopterans, we constructed the Ubx overexpression (Ubx^OE) strain based on the Nistari strain of Bombyx mori. The Ubx^OE strain showed a small body size, transparent intersegmental membrane and abnormal posterior silk gland (PSG). In the current study, we focused on the effect of Ubx overexpression on the posterior silk gland. As the major protein product of PSG, the mRNA expression of fibroin heavy chain (Fib-H) and fibroin light chain (Fib-L) was upregulated three times in Ubx^OE, but the protein expression of Fib-H and Fib-L was not significantly different. We speculated that the overexpression of Ubx downregulated the expression of Myc and further caused abnormal synthesis of the spliceosome and ribosome. Abnormalities of the spliceosome and ribosome affected the synthesis of protein in the PSG and changed its morphology.

Genome Editing of Silkworms

Silkworm is a lepidopteran insect that has been used as a model for a wide variety of biological studies. The microinjection technique is available, and it is possible to cause transgenesis as well as target gene disruption via the genome editing technique. TALEN-mediated knockout is especially effective in this species. We also succeeded in the precise and efficient integration of a donor vector using the precise integration into target chromosome (PITCh) method. Here we describe protocols for ZFN (zinc finger nuclease)-, TALEN (transcription activator-like effector nuclease)-, and CRISPR/Cas9-mediated genome editing as well as the PITCh technique in the silkworm. We consider that all of these techniques can contribute to the further promotion of various biological studies in the silkworm and other insect species.

Comparative Proteomic Analysis Provides New Insights into the Molecular Basis of Thermal-Induced Parthenogenesis in Silkworm (Bombyx mori)

Artificial parthenogenetic induction via thermal stimuli in silkworm is an important technique that has been used in sericultural production. However, the molecular mechanism underlying it remains largely unknown. We have created a fully parthenogenetic line (PL) with more than 85% occurrence and 80% hatching rate via hot water treatment and genetic selection, while the parent amphigenetic line (AL) has less than 30% pigmentation rate and less than 1% hatching rate when undergoing the same treatment. Here, isobaric tags for relative and absolute quantitation (iTRAQ)-based analysis were used to investigate the key proteins and pathways associated with silkworm parthenogenesis. We uncovered the unique proteomic features of unfertilized eggs in PL. In total, 274 increased abundance proteins and 211 decreased abundance proteins were identified relative to AL before thermal induction. Function analysis displayed an increased level of translation and metabolism in PL. After thermal induction, 97 increased abundance proteins and 187 decreased abundance proteins were identified. An increase in stress response-related proteins and decrease in energy metabolism suggested that PL has a more effective response to buffer the thermal stress than AL. Cell cycle-related proteins, including histones, and spindle-related proteins were decreased in PL, indicating an important role of this decrease in the process of ameiotic parthenogenesis.

SilkBase: an integrated transcriptomic and genomic database for Bombyx mori and related species

We introduce SilkBase as an integrated database for transcriptomic and genomic resources of the domesticated silkworm Bombyx mori and related species. SilkBase is the oldest B. mori database that was originally established as the expressed sequence tag database since 1999. Here, we upgraded the database by including the datasets of the newly assembled B. mori complete genome sequence, predicted gene models, bacterial artificial chromosome (BAC)-end and fosmid-end sequences, complementary DNA (cDNA) reads from 69 libraries, RNA-seq data from 10 libraries, PIWI-interacting RNAs (piRNAs) from 13 libraries, ChIP-seq data of 9 histone modifications and HP1 proteins and transcriptome and/or genome data of four B. mori-related species, i.e. Bombyx mandarina, Trilocha varians, Ernolatia moorei and Samia ricini. Our new integrated genome browser easily provides a snapshot of tissue- and stage-specific gene expression, alternative splicing, production of piRNAs and histone modifications at the gene locus of interest. Moreover, SilkBase is useful for performing comparative studies among five closely related lepidopteran insects.

Database URL: https://silkbase.ab.a.u-tokyo.ac.jp

Transcriptome analysis in the silkworm Bombyx mori overexpressing piRNA-resistant Masculinizer gene

Dosage compensation is a process that produces a similar expression of sex-linked and autosomal genes. In the silkworm Bombyx mori with a WZ sex-determination system, the expression from the single Z in WZ females matches that of ZZ males due to the suppression of Z-linked genes in males. A primary maleness determinant gene, Masculinizer (Masc), is also required for dosage compensation. In females, silkworm Piwi is complexed with the W chromosome-derived female-specific Feminizer (Fem) PIWI-interacting RNA (piRNA) and cleaves Masc mRNA. When Fem piRNA-resistant Masc cDNA (Masc-R) is overexpressed in both sexes, only female larvae are dead during the larval stage. In this study, transcriptome analysis was performed in neonate larvae to examine the effects of Masc-R overexpression on a global gene expression profile. Z-linked genes were globally repressed in Masc-R-overexpressing females due to force-driven dosage compensation. In contrast, Masc-R overexpression had little effect on the expression of Z-linked genes and the male-specific isoform of B. mori insulin-like growth factor II mRNA-binding protein in males, indicating that excessive Masc expression strengthens neither dosage compensation nor maleness in males. Fourteen genes were differentially expressed between Masc-R-overexpressing and control neonate larvae in both sexes, suggesting Masc functions other than dosage compensation and masculinization.

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.

InsectBase 2.0: a comprehensive gene resource for insects

Insects are the largest group of animals on the planet and have a huge impact on human life by providing resources, transmitting diseases, and damaging agricultural crop production. Recently, a large amount of insect genome and gene data has been generated. A comprehensive database is highly desirable for managing, sharing, and mining these resources. Here, we present an updated database, InsectBase 2.0 (http://v2.insect-genome.com/), covering 815 insect genomes, 25 805 transcriptomes and >16 million genes, including 15 045 111 coding sequences, 3 436 022 3'UTRs, 4 345 664 5'UTRs, 112 162 miRNAs and 1 293 430 lncRNAs. In addition, we used an in-house standard pipeline to annotate 1 434 653 genes belonging to 164 gene families; 215 986 potential horizontally transferred genes; and 419 KEGG pathways. Web services such as BLAST, JBrowse2 and Synteny Viewer are provided for searching and visualization. InsectBase 2.0 serves as a valuable platform for entomologists and researchers in the related communities of animal evolution and invertebrate comparative genomics.

Reference Transcriptome Data in Silkworm Bombyx mori

Herein, we performed RNA-seq analysis of ten major tissues/subparts of silkworm larvae. The sequences were mapped onto the reference genome assembly and the reference transcriptome data were successfully constructed. The reference data provided a nearly complete sequence for sericin-1, a major silk gene with a complex structure. We also markedly improved the gene model for other genes. The transcriptomic expression was investigated in each tissue and a number of transcripts were identified that were exclusively expressed in tissues such as the testis. Transcripts strongly expressed in the midgut formed tight genomic clusters, suggesting that they originated from tandem gene duplication. Transcriptional factor genes expressed in specific tissues or the silk gland subparts were also identified. We successfully constructed reference transcriptome data in the silkworm and found that a number of transcripts showed unique expression profiles. These results will facilitate basic studies on the silkworm and accelerate its applications, which will contribute to further advances in lepidopteran and entomological research as well as the practical use of these insects.

The Essential and Enigmatic Role of ABC Transporters in Bt Resistance of Noctuids and Other Insect Pests of Agriculture

Simple Summary

The insect family, Noctuidae, contains some of the most damaging pests of agriculture, including bollworms, budworms, and armyworms. Transgenic cotton and maize expressing Cry-type insecticidal proteins from Bacillus thuringiensis (Bt) are protected from such pests and greatly reduce the need for chemical insecticides. However, evolution of Bt resistance in the insects threatens the sustainability of this environmentally beneficial pest control strategy. Understanding the interaction between Bt toxins and their targets in the insect midgut is necessary to evaluate the risk of resistance evolution. ABC transporters, which in eukaryotes typically expel small molecules from cells, have recently been proposed as a target for the pore-forming Cry toxins. Here we review the literature surrounding this hypothesis in noctuids and other insects. Appreciation of the critical role of ABC transporters will be useful in discovering counterstrategies to resistance, which is already evolving in some field populations of noctuids and other insects.

Abstract

In the last ten years, ABC transporters have emerged as unexpected yet significant contributors to pest resistance to insecticidal pore-forming proteins from Bacillus thuringiensis (Bt). Evidence includes the presence of mutations in resistant insects, heterologous expression to probe interactions with the three-domain Cry toxins, and CRISPR/Cas9 knockouts. Yet the mechanisms by which ABC transporters facilitate pore formation remain obscure. The three major classes of Cry toxins used in agriculture have been found to target the three major classes of ABC transporters, which requires a mechanistic explanation. Many other families of bacterial pore-forming toxins exhibit conformational changes in their mode of action, which are not yet described for the Cry toxins. Three-dimensional structures of the relevant ABC transporters, the multimeric pore in the membrane, and other proteins that assist in the process are required to test the hypothesis that the ATP-switch mechanism provides a motive force that drives Cry toxins into the membrane. Knowledge of the mechanism of pore insertion will be required to combat the resistance that is now evolving in field populations of insects, including noctuids.