High-quality genome assembly of the silkworm, Bombyx mori. (April 2019)

• Record Type: Journal Article
• Authors: Kawamoto, Munetaka; Jouraku, Akiya; Toyoda, Atsushi; Yokoi, Kakeru; Minakuchi, Yohei; Katsuma, Susumu; Fujiyama, Asao; Kiuchi, Takashi; Yamamoto, Kimiko; Shimada, Toru
• Is Part Of: Insect biochemistry and molecular biology. Volume 107(2019)
• Page Start: 53
• Page End: 62
• Publication Date: 2019-04
• Subjects: Silkworm (Bombyx mori) -- Genome assembly -- Long-read sequencing -- Gene prediction; ABC ATP-binding cassette -- COE carboxylesterase -- GST glutathione S-transferase -- piRNA PIWI-interacting RNA -- TE transposable element; Insect biochemistry -- Periodicals; Insects -- Physiology -- Periodicals; Insects -- Molecular aspects -- Periodicals; Biochemistry -- Periodicals; Insectes -- Biochimie -- Périodiques; Insectes -- Composition -- Périodiques; Insectes -- Physiologie -- Périodiques; Insectes -- Aspect moléculaire -- Périodiques; Biochimie -- Périodiques; Biochemistry; Insect biochemistry; Insects -- Molecular aspects; Insects -- Physiology; Periodicals; 572.8157
• DOI: 10.1016/j.ibmb.2019.02.002
• Languages: English
• ISSNs: 0965-1748

Abstract

Abstract: 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 ). Graphical abstract: Image 1 Highlights: Whole silkworm genome resequencing was performed using PacBio long-read and Illumina short-read sequencers. High quality genome assembly (460.3 Mb) containing 696 scaffolds with an N50 size of 16.8 Mb was obtained. Gene prediction using the newly assembled genome identified 16, 880 genes with an N50 of 2154 bp. The new genome assembly and the gene models were evaluated. The new genome assembly and the gene models are available in SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp ).