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.

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.

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

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

Analysis of the Toll and Spaetzle Genes Involved in Toll Pathway-Dependent Antimicrobial Gene Induction in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae).

Insects rely only on their innate immune system to protect themselves from pathogens. Antimicrobial peptide (AMP) production is the main immune reaction in insects. In Drosophila melanogaster, the reaction is regulated mainly by the Toll and immune deficiency (IMD) pathways. Spaetzle proteins, activated by immune signals from upstream components, bind to Toll proteins, thus, activating the Toll pathway, which in turn, induces AMP genes. Previous studies have shown the difference in immune systems related to Toll and IMD pathways between D. melanogaster and Tribolium castaneum. In T. castaneum, nine Toll and seven spaetzle (spz) genes were identified. To extend our understanding of AMP production by T. castaneum, we conducted functional assays of Toll and spaetzle genes related to Toll-pathway-dependent AMP gene expression in T. castaneum under challenge with bacteria or budding yeast. The results revealed that Toll3 and Toll4 double-knockdown and spz7 knockdown strongly and moderately reduced the Toll-pathway-dependent expression of AMP genes, respectively. Moreover, Toll3 and Toll4 double-knockdown pupae more rapidly succumbed to entomopathogenic bacteria than the control pupae, but spz7 knockdown pupae did not. The results suggest that Toll3 and Toll4 play a large role in Toll-pathway-dependent immune reactions, whereas spz7 plays a small part.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae).

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

Meta-Analysis of Transcriptomes in Insects Showing Density-Dependent Polyphenism.

With increasing public data, a statistical analysis approach called meta-analysis, which combines transcriptome results obtained from multiple studies, has succeeded in providing novel insights into targeted biological processes. Locusts and aphids are representative of insect groups that exhibit density-dependent plasticity. Although the physiological mechanisms underlying density-dependent polyphenism have been identified in aphids and locusts, the underlying molecular mechanisms remain largely unknown. In this study, we performed a meta-analysis of public transcriptomes to gain additional insights into the molecular underpinning of density-dependent plasticity. We collected RNA sequencing data of aphids and locusts from public databases and detected differentially expressed genes (DEGs) between crowded and isolated conditions. Gene set enrichment analysis was performed to reveal the characteristics of the DEGs. DNA replication (GO:0006260), DNA metabolic processes (GO:0006259), and mitotic cell cycle (GO:0000278) were enriched in response to crowded conditions. To date, these processes have scarcely been the focus of research. The importance of the oxidative stress response and neurological system modifications under isolated conditions has been highlighted. These biological processes, clarified by meta-analysis, are thought to play key roles in the regulation of density-dependent plasticity.

Meta-Analysis of the Public RNA-Seq Data of the Western Honeybee Apis mellifera to Construct Reference Transcriptome Data.

The Western honeybee (Apis mellifera) is valuable in biological research and agriculture. Its genome sequence was published before those for other insect species. RNA-Seq data for A. mellifera have been applied in several recently published studies. Nevertheless, these data have not been prepared for use in subsequent meta-analyses. To promote A. mellifera transcriptome analysis, we constructed reference transcriptome data using the reference genome sequence and RNA-Seq data curated from about 1,000 runs of public databases. The new reference transcriptome data construct comprised 149,685 transcripts, and 194,174 protein sequences were predicted. Approximately 50-60% of the predicted protein sequences were functionally annotated using the protein sequence data for several model and insect species. Novel candidate immune-related transcripts were searched by meta-analysis using immune-response-related RNA-Seq and reference transcriptome data. Three to twenty candidate transcripts including autophagy-related protein 3 were upregulated or downregulated in response to both viral and bacterial infections. The constructed reference transcriptome data may facilitate future transcriptome analyses of A. mellifera.

Revealing Landscapes of Transposable Elements in Apis Species by Meta-Analysis.

Transposable elements (TEs) are grouped into several families with diverse sequences. Owing to their diversity, studies involving the detection, classification, and annotation of TEs are difficult tasks. Moreover, simple comparisons of TEs among different species with different methods can lead to misinterpretations. The genome data of several honey bee (Apis) species are available in public databases. Therefore, we conducted a meta-analysis of TEs, using 11 sets of genome data for Apis species, in order to establish data of "landscape of TEs". Consensus TE sequences were constructed and their distributions in the Apis genomes were determined. Our results showed that TEs belonged to four to seven TE families among 13 and 15 families of TEs detected in classes I and II respectively mainly consisted of Apis TEs and that more DNA/TcMar-Mariner consensus sequences and copies were present in all Apis genomes tested. In addition, more consensus sequences and copy numbers of DNA/TcMar-Mariner were detected in Apis mellifera than in other Apis species. These results suggest that TcMar-Mariner might exert A. mellifera-specific effects on the host A. mellifera species. In conclusion, our unified approach enabled comparison of Apis genome sequences to determine the TE landscape, which provide novel evolutionary insights into Apis species.

Co-Expression Network and Time-Course Expression Analyses to Identify Silk Protein Regulatory Factors in Bombyx mori.

Bombyx mori is an important economic insect and an animal model in pharmacomedical research. Although its physiology has been studied for many years, the mechanism via which silk protein genes are regulated remains unclear. In this study, we performed two-step expression screening, namely co-expression network and time-course expression analyses to screen silk protein regulation factors. A co-expression network analysis using RNA-seq data that were obtained from various tissues, including the silk glands of B. mori, was performed to identify novel silk protein regulatory factors. Overall, 91 transcription factors, including some known ones, were found to be co-expressed with silk protein genes. Furthermore, time-course expression analysis during the fifth instar larvae stage revealed that the expression pattern of 13 novel transcription factors was highly relevant to that of silk protein genes and their known regulatory factor genes. In particular, the expression peak of several transcription factors (TFs) was detected before the expression of silk protein genes peak. These results indicated that a larger number of genes than expected may be involved in silk protein regulation in B. mori. Functional analyses of function-unknown transcription factors should enhance our understanding of this system.

A simple method for ex vivo honey bee cell culture capable of in vitro gene expression analysis.

Cultured cells are a very powerful tool for investigating biological events in vitro; therefore, cell lines have been established not only in model insect species, but also in non-model species. However, there are few reports on the establishment of stable cell lines and development of systems to introduce genes into the cultured cells of the honey bee (Apis mellifera). We describe a simple ex vivo cell culture system for the honey bee. Hemocyte cells obtained from third and fourth instar larvae were cultured in commercial Grace's insect medium or MGM-450 insect medium for more than two weeks maintaining a normal morphology without deterioration. After an expression plasmid vector bearing the enhanced green fluorescent protein (egfp) gene driven by the immediate early 2 (IE2) viral promoter was transfected into cells, EGFP fluorescence was detected in cells for more than one week from one day after transfection. Furthermore, double-stranded RNA corresponding to a part of the egfp gene was successfully introduced into cells and interfered with egfp gene expression. A convenient and reproducible method for an ex vivo cell culture that is fully practicable for gene expression assays was established for the honey bee.

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.

Bumble bee queens activate dopamine production and gene expression in nutritional signaling pathways in the brain.

To explore the neuroendocrine mechanisms underlying caste-specific behavior and its evolution from primitive to advanced eusocial bees, the monoamine levels and expression of genes involved in monoamine production and signaling in the brain were compared between the castes of Bombus ignitus. Higher levels of dopamine and its related substances were found in the brains of newly emerged queens than in the brains of emerged workers. The degree of caste differences in B. ignitus was smaller than that reported in Apis mellifera, indicating a link to different social stages in the two species. There was no differential expression in genes involved in dopamine biosynthesis between castes, suggesting that the high dopamine production in queens was not largely influenced by the expression of these genes at emergence, rather it might be influenced by tyrosine supply. Genome-wide analyses of gene expression by RNA-sequencing indicated that a greater number of genes involved in nutrition were actively expressed in the brains of newly emerged queens in comparison to the emerged workers. Some of the expression was confirmed by real-time quantitative PCR. The signaling pathways driven by the expression of these genes may be associated with dopamine signaling or the parallel activation of dopamine production.

An update of KAIKObase, the silkworm genome database.

KAIKObase was established in 2009 as the genome database of the domesticated silkworm Bombyx mori. It provides several gene sets and genetic maps as well as genome annotation obtained from the sequencing project of the International Silkworm Genome Consortium in 2008. KAIKObase has been used widely for silkworm and insect studies even though there are some erroneous predicted genes due to misassembly and gaps in the genome. In 2019, we released a new silkworm genome assembly, showing improvements in gap closure and covering more and longer gene models. Therefore, there is a need to include new genome and new gene models to KAIKObase. In this article, we present the updated contents of KAIKObase and the methods to generate, integrate and analyze the data sets. Database URL: https://kaikobase.dna.affrc.go.jp.

Effectiveness of orally-delivered double-stranded RNA on gene silencing in the stinkbug Plautia stali.

Development of a reliable method for RNA interference (RNAi) by orally-delivered double-stranded RNA (dsRNA) is potentially promising for crop protection. Considering that RNAi efficiency considerably varies among different insect species, it is important to seek for the practical conditions under which dsRNA-mediated RNAi effectively works against each pest insect. Here we investigated RNAi efficiency in the brown-winged green stinkbug Plautia stali, which is notorious for infesting various fruits and crop plants. Microinjection of dsRNA into P. stali revealed high RNAi efficiency-injection of only 30 ng dsRNA into last-instar nymphs was sufficient to knockdown target genes as manifested by their phenotypes, and injection of 300 ng dsRNA suppressed the gene expression levels by 80% to 99.9%. Knockdown experiments by dsRNA injection showed that multicopper oxidase 2 (MCO2), vacuolar ATPase (vATPase), inhibitor of apoptosis (IAP), and vacuolar-sorting protein Snf7 are essential for survival of P. stali, as has been demonstrated in other insects. By contrast, P. stali exhibited very low RNAi efficiency when dsRNA was orally administered. When 1000 ng/µL of dsRNA solution was orally provided to first-instar nymphs, no obvious phenotypes were observed. Consistent with this, RT-qPCR showed that the gene expression levels were not affected. A higher concentration of dsRNA (5000 ng/µL) induced mortality in some cohorts, and the gene expression levels were reduced to nearly 50%. Simultaneous oral administration of dsRNA against potential RNAi blocker genes did not improve the RNAi efficiency of the target genes. In conclusion, P. stali shows high sensitivity to RNAi with injected dsRNA but, unlike the allied pest stinkbugs Halyomorpha halys and Nezara viridula, very low sensitivity to RNAi with orally-delivered dsRNA, which highlights the varied sensitivity to RNAi across different species and limits the applicability of the molecular tool for controlling this specific insect pest.

Genome-wide assessment and development of molecular diagnostic methods for imidacloprid-resistance in the brown planthopper, Nilaparvata lugens (Hemiptera; Delphacidae).

BACKGROUND: The brown planthopper, Nilaparvata lugens (Stål), is one of the most notorious pests of rice throughout Asia. The brown planthopper has developed high resistance to imidacloprid, a member of neonicotinoid insecticides. Several genes and mutations conferring imidacloprid resistance in N. lugens, especially in eastern and southeastern Asia populations, have been reported. Thus, the key mechanisms of imidacloprid resistance need to be examined. RESULTS: RNA-seq analyses revealed that only one cytochrome P450 monooxygenase gene, CYP6ER1, was commonly upregulated in the five resistant strains tested. Sequences of CYP6ER1, which were highly expressed in the imidacloprid-resistant strains, contained a three-nucleotide deletion in the coding region, and amino acid substitutions and deletion, compared to that in an imidacloprid-susceptible strain. RNAi-mediated gene knockdown of CYP6ER1 increased imidacloprid susceptibility in a resistant strain. Further, we established two simple and convenient PCR-based molecular diagnostic methods to detect the CYP6ER1 locus with the three-nucleotide deletion. Using these methods, the resistance of F2 progenies derived from the crosses of F1 siblings from susceptible and resistant parents was analyzed, showing that the imidacloprid resistance had a relationship to the CYP6ER1 locus with the three-nucleotide deletion. CONCLUSION: The overexpression of a variant CYP6ER1 with amino acid substitutions and deletion was involved in imidacloprid resistance in N. lugens. Based on these findings, molecular diagnostic methods have been developed and are promising tools for monitoring imidacloprid resistance in paddy fields. © 2020 Society of Chemical Industry.

Diversity and function of multicopper oxidase genes in the stinkbug Plautia stali.

Multicopper oxidase (MCO) genes comprise multigene families in bacteria, fungi, plants and animals. Two families of MCO genes, MCO1 (laccase1) and MCO2 (laccase2), are conserved among diverse insects and relatively well-characterized, whereas additional MCO genes, whose biological functions have been poorly understood, are also found in some insects. Previous studies reported that MCO1 participates in gut immunity and MCO2 plays important roles in cuticle sclerotization and pigmentation of insects. In mosquitoes, MCO2 was reported to be involved in eggshell sclerotization and pigmentation, on the ground that knockdown of MCO2 caused deformity and fragility of the eggshell. Here we identified a total of 7 MCO genes, including PsMCO1 and PsMCO2, and investigated their expression and function in the brown-winged green stinkbug Plautia stali. RNA interference (RNAi) knockdown of MCO genes by injecting double-stranded RNA (dsRNA) into nymphs revealed that MCO2, but not the other 6 MCOs, is required for cuticle sclerotization and pigmentation, and also for survival of P. stali. Trans-generational knockdown of MCO2 by injecting dsRNA into adult females (maternal RNAi) resulted in the production of unhatched eggs despite the absence of deformity or fragility of the eggshell. These results suggested that MCO2 plays an important role in sclerotization and pigmentation of the cuticle but not in eggshell integrity in P. stali. Maternal RNAi of any of the other 6 MCO genes and 3 tyrosinase genes affected neither survival nor eggshell integrity of P. stali. Contrary to the observations in the red flour beetle and the brown rice planthopper, RNAi knockdown of MCO6 (MCORP; Multicopper oxidase related protein) exhibited no lethal effects on P. stali. Taken together, our findings provide insight into the functional diversity and commonality of MCOs across hemipteran and other insect groups.

Genome-wide SNP marker discovery and phylogenetic analysis of mulberry varieties using double-digest restriction site-associated DNA sequencing.

There has been long taxonomic debate on mulberry species (genus Morus) because the classification of mulberry species has relied on morphological characteristics. Although attempts for classifying mulberry species using molecular markers have been performed, phylogenetic relationships among diploid mulberry species remain unclear. In this study, we aim to investigate the genetic relationship between 54 diploid mulberry varieties belonging to seven different Morus species (M. alba, M. indica, M. bombycis, M. acidosa, M. latifolia, M. kagayamae, and M. rotundiloba) and one unspecified Morus species ('Enbu') using genome-wide SNP discovery and phylogenetic analysis via double-digest restriction site-associated DNA sequencing (ddRAD-seq). Genome-wide 2229 homozygous SNPs of 54 mulberry varieties in the eight species were identified by ddRAD-seq. Results of the phylogenetic analysis identified only three clear monophyletic clades in two Japanese native species, M. acidosa and M. kagayamae, which are found on different geographically isolated islands and a Thai species, M. rotundiloba, whereas the other species were non-monophyletic. Varieties of M. bombycis, another Japanese native species, were roughly classified into three groups. Of these, two M. bombycis groups were monophyletic with M. acidosa and M. kagayamae, respectively, while another M. bombycis group was not monophyletic. Varieties of M. indica, an Indian native species, were classified into two different monophyletic clades. Of these, one clade was clearly monophyletic with an indigenous variety in Kenya, 'Enbu', while another clade was monophyletic with M. rotundiloba and one M. latifolia variety. There were no clear monophyletic clades within M. alba and M. latifolia varieties, which could be a result of several hybridization events after their introductions from China to Japan. Our results suggested that it was difficult to clearly classify the hybridized mulberry varieties even with genome-wide DNA markers. In addition to phylogenetic analysis, we also evaluated morphological characteristics of mulberry leaves for each variety. The results of morphological evaluation indicated that leaf tip ratio may correlate to genetic difference among the two M. bombycis groups in monophyletic clades and another M. bombycis group in non-monophyletic clades. These results suggested that leaf tip ratio might be used for evaluating hybridization of M. bombycis varieties. Over all, our results may provide new insights into taxonomic debate of mulberry species.

Functional crosstalk across IMD and Toll pathways: insight into the evolution of incomplete immune cascades.

In insects, antimicrobial humoral immunity is governed by two distinct gene cascades, IMD pathway mainly targeting Gram-negative bacteria and Toll pathway preferentially targeting Gram-positive bacteria, which are widely conserved among diverse metazoans. However, recent genomic studies uncovered that IMD pathway is exceptionally absent in some hemipteran lineages like aphids and assassin bugs. How the apparently incomplete immune pathways have evolved with functionality is of interest. Here we report the discovery that, in the hemipteran stinkbug Plautia stali, both IMD and Toll pathways are present but their functional differentiation is blurred. Injection of Gram-negative bacteria and Gram-positive bacteria upregulated effector genes of both pathways. Notably, RNAi experiments unveiled significant functional permeation and crosstalk between IMD and Toll pathways: RNAi of IMD pathway genes suppressed upregulation of effector molecules of both pathways, where the suppression was more remarkable for IMD effectors; and RNAi of Toll pathway genes reduced upregulation of effector molecules of both pathways, where the suppression was more conspicuous for Toll effectors. These results suggest the possibility that, in hemipterans and other arthropods, IMD and Toll pathways are intertwined to target wider and overlapping arrays of microbes, which might have predisposed and facilitated the evolution of incomplete immune pathways.

Genome-wide Identification of Tebufenozide Resistant Genes in the smaller tea tortrix, Adoxophyes honmai (Lepidoptera: Tortricidae).

The smaller tea tortrix, Adoxophyes honmai, has developed strong resistance to tebufenozide, a diacylhydrazine-type (DAH) insecticide. Here, we investigated its mechanism by identifying genes responsible for the tebufenozide resistance using various next generation sequencing techniques. First, double-digest restriction site-associated DNA sequencing (ddRAD-seq) identified two candidate loci. Then, synteny analyses using A. honmai draft genome sequences revealed that one locus contained the ecdysone receptor gene (EcR) and the other multiple CYP9A subfamily P450 genes. RNA-seq and direct sequencing of EcR cDNAs found a single nucleotide polymorphism (SNP), which was tightly linked to tebufenozide resistance and generated an amino acid substitution in the ligand-binding domain. The binding affinity to tebufenozide was about 4 times lower in in vitro translated EcR of the resistant strain than in the susceptible strain. RNA-seq analyses identified commonly up-regulated genes in resistant strains, including CYP9A and choline/carboxylesterase (CCE) genes. RT-qPCR analysis and bioassays showed that the expression levels of several CYP9A and CCE genes were moderately correlated with tebufenozide resistance. Collectively, these results suggest that the reduced binding affinity of EcR is the main factor and the enhanced detoxification activity by some CYP9As and CCEs plays a supplementary role in tebufenozide resistance in A. honmai.

High-quality genome assembly of the silkworm, Bombyx mori.

In 2008, the genome assembly and gene models for the domestic silkworm, Bombyx mori, were published by a Japanese and Chinese collaboration group. However, the genome assembly contains a non-negligible number of misassembled and gap regions due to the presence of many repetitive sequences within the silkworm genome. The erroneous genome assembly occasionally causes incorrect gene prediction. Here we performed hybrid assembly based on 140 × deep sequencing of long (PacBio) and short (Illumina) reads. The remaining gaps in the initial genome assembly were closed using BAC and Fosmid sequences, giving a new total length of 460.3 Mb, with 30 gap regions and an N50 comprising 16.8 Mb in scaffolds and 12.2 Mb in contigs. More RNA-seq and piRNA-seq reads were mapped on the new genome assembly compared with the previous version, indicating that the new genome assembly covers more transcribed regions, including repetitive elements. We performed gene prediction based on the new genome assembly using available mRNA and protein sequence data. The number of gene models was 16,880 with an N50 of 2154 bp. The new gene models reflected more accurate coding sequences and gene sets than old ones. The proportion of repetitive elements was also reestimated using the new genome assembly, and was calculated to be 46.8% in the silkworm genome. The new genome assembly and gene models are provided in SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp).