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

Abstract

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.

Introduction

Bombyx mori (B. mori) or domestic silkworm is a well-known bioreactor that produces natural silk in sericulture and is frequently used as a lepidopteran model to study insect immunology and disease resistance (Maekawa et al., 1988; Arunkumar et al., 2006). In recent decades, the genetics and genomics of B. mori in several Asian countries have been greatly elucidated (Xia et al., 2004). Its natural primal enemy, B. mori nucleopolyhedrovirus (BmNPV), can cause widespread death to the overwhelming majority of B. mori strains. Several strains bred in different laboratories display congenital BmNPV resistance (Xu et al., 2013). Similar to most higher organisms, the key prerequisite of most biological research is the assembly of whole-genome sequences. For B. mori, the first genome sequence with 3× coverage generated by shotgun sequencing technology (fosmid-end) was published in Japan in 2004, and the B. mori strain used was p50T (Koike et al., 2003; Mita et al., 2004). At almost the same time, an independent Chinese group also performed a similar whole-genome sequencing (BAC-end) with 5.9× coverage, and the B. mori strain used was Dazao (Xia et al., 2004). Both studies declared that the estimated genome coverage was more than 90%. Based on the EST database in the former project, the first version of SilkBase was released in 2006. Then, the integrated transcriptomic and genomic data were added to the project and released in the second version of SilkBase in 2015. Two years later, the p50T strain was sequenced again by making use of the combination of PacBio and Illumina sequencing technologies, and this genome assembly with increased accuracy was released as SilkBase v2.1 (Kawamoto et al., 2019). Furthermore, the Silkworm Knowledge-based database (SilkDB) was developed for data storage, retrieval, analysis, and visualization by Chinese groups in 2005 (Wang et al., 2005), and the assembly accuracy was later promoted by the International Silkworm Genome Consortium in 2010 (SilkDB 2.0) (International Silkworm Genome, 2008). Ten years later, cutting-edge PacBio and Hi-C biotechnologies were employed to assemble the latest version of the Dazao genome in 2020 (SilkDB 3.0) (Duan et al., 2010; Lu et al., 2020). To the best of our knowledge, SilkDB 3.0 and SilkBase v2.1 are the most commonly used database reference genomes. Hence, they are recognized as the basis for a number of B. mori studies and databases (Cao et al., 2017; Li et al., 2019; Zhu et al., 2019; Fujimoto et al., 2020).

It is a pity that both of the strains they used do not have BmNPV-resistance. Perhaps, as a consequence, their genome information is not perfectly accurate for the resistance studies considering the importance of the reference genome. Since countless efforts have been made for the BmNPV resistance without any acknowledged results, we think that it is time to optimize the reference genome group—for example, establish a BmNPV-resistant genome. In addition, in our opinion, the pipelines for constructing the above-mentioned two genome references could be more accurate and practical. For SilkDB 3.0, short-read sequencing data that could correct the draft genome assembly were not added to aid the assembly process (Jayakumar and Sakakibara, 2019). For SilkBase, Hi-C data should be used to help arrange and orient contigs or scaffolds (Belton et al., 2012). More importantly, neither Dazao nor p50T has resistance to BmNPV, which may increase genomic noise or even introduce errors in the resistance study of B. mori. Hence, here we introduce a high-quality de novo chromosome-level genome assembly of B. mori with BmNPV resistance by integrating data from PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing technology.

Materials and Methods

Sample Collection and DNA Isolation

Two five-instar B. mori belonging to the NB strain were collected from the School of Life Sciences, Jiangsu University, and downstream wet-laboratory experiments were performed at Novogene Co., Ltd. First, one of the imagoes was sufficiently ground with polyvinyl polypyrrolidone powder, and whole DNA was extracted for sequencing using sodium dodecyl sulfate (SDS) and cetyltrimethylammonium bromide (Chen et al., 2010; Kasajima, 2018). Then, 1% agarose gel was used to examine DNA degradation and contamination. Consequently, the DNA purity was checked using a NanoPhotometer^® spectrophotometer (IMPLEN, CA, United States), and the DNA concentration was measured using a Qubit^® DNA Assay Kit on a Qubit^® 2.0 fluorometer (Life Technologies, CA, United States).

Library Construction and Illumina Sequencing

A total of 1.5 μg extracted DNA was used as input material for the DNA sample preparations and sheared with a Covaris Focused ultrasonicator. Following the recommendations of the manufacturer, a library for Next-Gen sequencing was constructed using a Truseq^® Nano DNA Sample Preparation Kit (Illumina, San Diego, CA, United States) after end-repair and the addition of index code. Subsequently, the fragments were end-polished, A-tailed, and ligated with the full-length adapter for Illumina sequencing with further PCR amplification. After fragment selection for a size of 350 bp and PCR amplification, the products were purified with the AMPure XP system (PacBio, Menlo Park, CA, United States), and the library was analyzed for size distribution by an Agilent 2100 Bioanalyzer and quantified using real-time PCR. Finally, the library constructed as detailed above was sequenced by an Illumina HiSeq X Ten instrument, and 150 bp paired-end (PE) reads, including ∼350 bp inserts, were generated. After removing the reads with adapters, the paired reads that had more than 10% N bases or over 20% low-quality bases on either end sequence were discarded. As a result, the remaining reads passed onto the downstream analysis as clean reads.

Library Construction and PacBio Sequencing

More than 5 μg of concentrated genomic DNA (gDNA) that passed the inspections was used for the size-select DNA fragment step. g-TUBE was utilized to shear the gDNA into ∼20-kb fragments to construct a 20-kb SMRTbell library. After shearing, AMPure PB beads were used to concentrate the sheared gDNA. Then, the fragments were treated to remove the single-strand overhangs and repair DNA damage with ExoVII enzyme. Subsequently, the ends of the double-strand fragments were polished with T4 DNA Polymerase and then ligated to single-molecule real-time (SMRT) hairpin adapters. Through EXOIII and VII enzyme treatment, the SMRTbell templates were washed out. The library was again concentrated and purified with AMPure PB beads. In this case, we employed the BluePippin system (Sage Science, Inc.) and set a 20-kb cutoff threshold for size selection. The sublibrary was also concentrated and purified with AMPure PB beads. After annealing the sequencing primer to the SMRTbell template, polymerase was bound to both ends of the SMRTbell templates using a binding kit for efficient loading into ZWMs. Finally, a total of 14 SMRT cells were run on the PacBio Sequel system with the P6-C4 sequencing reagent (Kingan et al., 2019). Briefly, a SMRTbell library was constructed with the SMRTbell Express Template Prep Kit 2.0 according to the latest protocol published by Pacific Biosciences of California, Inc.

Hi-C Experiment and Sequencing

The Hi-C experiment was performed as described previously (Belton et al., 2012; Yardimci et al., 2019). Cell suspensions were made with the other B. mori imago and treated with paraformaldehyde to fix the three-dimensional structures of the nuclei to retain the relationship between genomic and physical distance. Therefore, after cell lysis with RIPA lysis buffer (1 M Tris–HCl, pH 8, 1 M NaCl, 10% CA-630, and 13 units of protease inhibitor), the exposed chromatin was cross-linked in situ to trap sequence interactions across the entire genome and between different chromosomes. Then, the supernatant was centrifuged at 5,000 rpm at room temperature for 10 min. The obtained pellet was washed twice with 100 μl ice cold 1× NEB buffer, followed by centrifugation at 5,000 rpm for 6 min.

The chromatin was resuspended in 100 μl NEB buffer and solubilized with dilute SDS, followed by incubation at 65°C for 10 min. After quenching SDS with Triton X-100, overnight digestion was performed with the four-cutter endonuclease MboI at 37°C on a rocking platform. Subsequently, the fragmented chromatin was digested with the restriction enzymes HindIII and DpnII at 37°C for 16 h. Next, for proximity ligation, the fragmented loci were ligated and marked with biotin to create chimeric junctions between adjacent sequences after incubation at 37°C for 45 min. The enzymatic system was inactivated with 20% SDS solution. Finally, proteinase K was added for reverse cross-linking overnight at 65°C, and the resulting sample was purified and dissolved in 90 μl of double-distilled water. The purified fragments were sheared to a size of ∼350 bp, and then their ends were repaired. The biotin-containing fragments were isolated with Dynabeads^TM M-280 Streptavidin purchased from Thermo Fisher Scientific Inc. After adding A-tails to the fragment ends and following ligation by Illumina PE sequencing adapters, the Hi-C library was amplified by 12–14 PCR cycles and sequenced on the Illumina NovaSeq 6000 platform with the PE set to 150 bp.

RNAseq Experiment and Transcriptome Sequencing

Transcriptome sequencing was performed as described by Conesa et al. (2016). The RNA pool was made from a whole one-instar B. mori using the RNeasy Maxi Kit purchased from QIAGEN, and raw sequencing data were obtained on the Illumina HiSeq X sequence platform. The subsequent quality control steps were performed with FastQC, FASTX-Toolkit, and Trimmomatic (Bolger et al., 2014).

Results

De novo Assembly of the B. mori Genome and Assessment

The raw Illumina sequencing data first underwent quality control after image recognition, contamination, and adapter elimination. Then, the initial characterization of the B. mori genome was estimated through k-mer (k = 17 in this study) analysis of the clean data by jellyfish and GenomeScope (Vurture et al., 2017). As a result, the estimated genome size (475.39 and 464.90 Mbp after later correction), heterozygosity (0.23%), and repeat content (43.78.09%) were basically consistent with previous findings (Mita et al., 2004; Xia et al., 2004; International Silkworm Genome, 2008; Kawamoto et al., 2019). All the sequencing data produced in this study are listed in Table 1. As the table shows, long-read sequencing on the PacBio Sequel system yielded 50.08 G of data with a high coverage of 107.72× short-read sequencing on the Illumina NovaSeq 6000 system yielded 64.81 G of data with an average coverage of 139.41× and Hi-C sequencing on the HiSeq X system yielded 83.78 G of data with an average coverage of 180.21×.

During the assembly process, the long reads downloaded from the PacBio Sequel platform were subjected to self-correction by the NextCorrect module of NextDenovo.¹ Taking the corrected reads as input materials, FALCON was employed to assemble the draft framework of the B. mori genome (Chin et al., 2016). To correct the high error rate of the long-read sequencing technology (Rang et al., 2018), two genome sequence polishing steps were performed: the Quiver algorithm was first used to polish the genome using PacBio long reads (Chin et al., 2013), and another round of genome-wide base-level correction of the Illumina clean reads was performed using Pilon (Walker et al., 2014). For chromosome-level scaffolding (Oluwadare et al., 2019), the clean Hi-C reads were mapped to the assembled genome using BWA (Li and Durbin, 2010), and the repeats and unmapped reads were removed by SAMtools (Li et al., 2009) to obtain high-quality reads and information about the restriction sites for enzyme cutting. As a result, only uniquely mapped read pairs were considered for subsequent analysis, and the Hi-C heatmap did not show misassembly during scaffolding (Supplementary Figure 1). LACHESIS was finally used to cluster, order, and orient the assembled contigs (Burton et al., 2013). The assembly statistics of the completed genome are shown in Table 2 and Supplementary Tables 1–3. Only scaffolds larger than 100 bp were selected to perform the assembly; the contig N50 was 3.75 Mbp, and the scaffold N50 was 17.26 Mbp.

Briefly, we assembled a high-quality chromosome-scale genome for B. mori de novo by making use of the 198.67 G of sequencing data (Figure 1). The overall length and N50 value of the contigs were 455.46 and 3.75 Mbp, respectively, while the overall length and N50 value of the scaffolds were 455.46 and 17.26 Mbp. To evaluate the quality of the newly assembled genome, its completeness was first assessed by the core eukaryotic genes mapping approach (CEGMA), whose core gene database includes 248 genes derived from six model eukaryotes (Parra et al., 2007). Supplementary Table 4 shows that 223 core eukaryotic genes were retrieved, which was up to 89.92% of the total number of genes. In addition to CEGMA, Benchmarking Universal Single-Copy Ortholog (BUSCO) software, invoking TBLASTN (Gertz et al., 2006), AUGUSTUS (Stanke and Morgenstern, 2005), and HMMER (Wheeler and Eddy, 2013), was employed (Seppey et al., 2019). Of the 978 orthologous genes, 98.1% complete BUSCOs (C + D) were detected in its report (see Supplementary Tables 4, 5). CEGMA and BUSCO revealed that we obtained a high-quality and complete B. mori genome.

FIGURE 1

Moreover, the reference genomes from SilkDB and SilkBase were used to check the collinearity with our genome assembly with JCVI (Tang et al., 2008). As expected, great collinearity existed among the three chromosome-level genomes (Figure 2). The distortions on several chromosomes, such as Chr1, Chr10, and Chr23, were attributable to scaffold perversions, which is quite reasonable, as large chromosomal structure changes may occur among closely related strains. In addition, optical mapping methods and approaches for the ordering and orientation of scaffolds can also cause a systematic deviation (Ekblom and Wolf, 2014; Belser et al., 2018; Deschamps et al., 2018). Nevertheless, crevices with different sizes indicate small INDELs that may be the origins of BmNPV resistance; this will be analyzed in detail in future studies in our laboratory.

FIGURE 2

Genome Annotation

Supplementary Figure 2 shows three routes in the road map of the genome annotation phase. From the left, an integration of homology-based and de novo approaches was used to identify the repeat sequences. First, LTR_FINDER² (Xu and Wang, 2007), RepeatScout,³ and RepeatModeler⁴ were employed to create a de novo repeat sequence dataset. Then, based on the RepBase database, RepeatMasker (see text footnote 3) and RepeatProteinMask⁵ were used to recognize the sequences that had a high similarity in RepBase. Moreover, Tandem Repeats Finder (Benson, 1999) and RepeatModeler2 (Flynn et al., 2020) were also utilized to detect tandem repeats and transposable element (TE) families in the assembled genome with default settings. Consequently, all the above-mentioned repetitive elements were combined to annotate the genome assembly using RepeatMasker (Table 3 and Supplementary Figure 3; Supplementary Table 6).

For the prediction of gene structures, de novo, homology-based, and transcriptome-based strategies were integrated to predict genes in the B. mori genome (Table 4 and Supplementary Figure 2). De novo prediction was performed with AUGUSTUS (Stanke and Morgenstern, 2005), GlimmerHMM (Majoros et al., 2004), SNAP (Johnson et al., 2008), Geneid,⁶ and Genscan (Burge and Karlin, 1997), while homology-based prediction was conducted with BLAT, which aligns the protein sequences of Ame, Pra, Bmo, Pxy, Tca, and Dme to the above-mentioned genome assembly (Kent, 2002). Then, the transcript reads were aligned to the genome assembly using the TopHat package (Kim et al., 2013), followed by gene structure prediction with Cufflinks (Ghosh and Chan, 2016). Finally, all gene models were merged by EVidenceModeler (EVM)⁷ (Haas et al., 2008) after removing redundancy. To enrich the contents, the RNAseq sequence data were assembled de novo with Trinity (Grabherr et al., 2011; Haas et al., 2013) and applied to insert untranslated regions and alternative splicing information to the EVM dataset by PASA (Haas et al., 2003). In conclusion, 9,113 gene structures were synchronously supported by de novo, homology-based, and transcriptome-based prediction strategies (Supplementary Figure 5). The same prediction pipeline was also applied to several proximal species (Bmj, Ame, Bmo, Dme, Pra, Pxy, and Tca; see Supplementary Table 10). Notably, the number of predicted genes for Bmo was 13,850, which was very close to the 13,103 predicted genes for Bmj.

For gene function annotation, the final protein sequences obtained above were used as inputs for the SwissProt,⁸ Nr,⁹ Kyoto Encyclopedia of Genes and Genomes (KEGG),¹⁰ InterPro,¹¹ Gene Ontology,¹² and Pfam¹³ databases. After alignment to these published protein databases, 99.6% of the genes were annotated with functional entities (Supplementary Table 8). In addition, non-coding RNAs were also annotated (Supplementary Table 9). tRNAs were detected with tRNAscan-SE¹⁴ (Lowe and Chan, 2016; Chan and Lowe, 2019), while miRNAs and snRNAs were detected with the INFERNAL module of Rfam¹⁵ (Kalvari et al., 2018, 2021).

Genomic Comparison

Genomic comparison in bioinformatics analysis, including gene family clustering and phylogenetic and divergence analysis, was performed by comparing our in-house B. mori genome assembly with the assemblies of a set of 17 other representative species. In addition to three closely related species (Bmo, Aya, and Bma) and five other lepidopteran species (Sfr, Pxy, Pra, Har, and Dpl), three dipteran species (Aya, Aae, and Dme), three hymenopteran species (Ame, Lhu, and Nvi), and two coleopteran species (Lde and Tca) were also included in the comparison list. Moreover, Caenorhabditis elegans, belonging to Rhabditia, was also included. For convenience, their abbreviations are summarized in Supplementary Table 10.

First, OrthoMCL¹⁶ was used for gene family clustering analysis (Li et al., 2003; Chen et al., 2006) using the following parameters: mode, 3; inflation, 1.5 blast_file, gg_file. As a result, a total of 21,904 gene families were identified, and 650 strict single-copy orthologs were recovered in the 18 genomes (Figure 3A). In addition, a core set of 6,897 gene families was shared by Bmj (for shipshape, use Bmj instead of B. mori), Bmo, Bma, and Aya (Figure 3B), and only 220 gene families uniquely belonged to Bmj. A consequent large-scale analysis found 84 gene families specific to Bmj but not the 17 other selected species. Most of these specific genes were involved in receptor activity, transmembrane signaling receptor activity, and glycosphingolipid biosynthesis (Supplementary Figure 7).

FIGURE 3

Second, CAFÉ v3.0¹⁷ and CAFE_fig¹⁸ were used to analyze changes in the gene family size using the following parameters (-i -t 8 -l -p 0.05). As Figure 3C shows, 132 gene families were expanded, while 1,021 gene families were contracted from the Bmj genome. In comparison with the closely interlinked Bmo, which gained 439 gene families and missed 461 gene families, Bmj lost many more gene families. These missed genes were rich in diverse functions, such as catalytic activity, fatty acid biosynthesis, fatty acid metabolism, and AMPK signaling pathway (Supplementary Figure 8). Corresponding to the phylogenetic analysis performed with the mcmctree program of the PAML package¹⁹ (parameters: rootage, 1,000; clock, 3; alpha, 0.603820) (Yang, 1997), Figure 3A shows that Bmj has a very close relationship to Bmo and phylogenetically diverged from the common ancestor ∼3 million years ago. However, the history of human breeding is up to 5,000 years long (Tabunoki et al., 2016; Meng et al., 2017), and the best explanation is manual intervention in sericulture, which increases their divergence. By using the phylogenetic analysis as a supplement, the above-mentioned 650 single-copy gene families were subjected to internal alignment with MUSCLE,²⁰ and then the aligned results were merged to form a super alignment matrix (Edgar, 2004a, b). After that, the ML TREE algorithm in RAxML²¹ was used with default parameters to create a phylogenetic tree (Figure 3C; Stamatakis, 2014).

To detect positively selected genes related to different phenotypes of Bmj, a positive selection analysis was executed in three diverse combinations with the other 17 species. MUSLCE was again employed to align the protein sequences of the single-copy gene families between foreground branches and background branches. Then, the aligned results were cleaned by removing the low-quality regions using Gblocks and then reversed to CDS²² (Castresana, 2000). For each gene family, the branch-site model in the codeml program of the PAML package was applied to check whether it was positively selected in the Bmj branch (Yang, 1997). Instead of simply seeking genes with ka/ks > 1, PAML makes use of two hypothesized likelihood ratios to double-check the positive selection (Yang and Nielsen, 2002; Zhang et al., 2005) (Figures 4A,B). For the eating pattern of mulberry (Ge et al., 2018a, b), the foreground branches were Bmj, Bmo, and Bma, and the background branches were Aya, Har, Sfr, Pra, Dpl, and Pxy in the first group. Sixty-two targeted genes that were involved in metabolic processes, catalytic activity, and Huntington’s disease were obtained (Supplementary Figure 9). For the capability of silk production, the foreground branches were Bmj, Bmo, Bma, and Aya, and the background branches were Lhu, Ame, Nvi, Lde, Tca, Aae, Aga, and Dme in the second group. A total of 119 targeted genes were enriched in metal ion binding, zinc ion binding, and transmembrane transport (Supplementary Figure 10). In the third group, for resistance to BmNPV, the foreground branch was Bmj, while the background branches were Bmo, Har, Sfr, Pra, Dpl, Pxy, Bma, and Aya. The downstream analysis indicated that the BmNPV resistance of our NB species may be related to pathways of transmembrane signaling receptor activity and lysosome and signal transducer activity (Figures 4C,D and Supplementary Figure 11).

FIGURE 4

Conclusion and Discussion

Since approximately 5,000 years ago, B. mori has been one of the most important agricultural economic insects for silk production in Asia (Xiang et al., 2018). As it is completely domesticated in sericulture, its survival and reproduction totally depend on human beings (Jiang and Xia, 2014). With advances of biotechnology, B. mori has been treated as an important bioreactor for producing recombinant protein (Tomita et al., 2000, 2003). Moreover, because of their intermediate genome size, short life cycle, affordability, and ability to be easily used for drug screening, B. mori is a perfect model organism for scientific discovery, especially for lepidopterans, and is useful for economic research (Kaito et al., 2002; International Silkworm Genome, 2008; Tabunoki et al., 2016; Meng et al., 2017). To maintain the consistency of biological background among different laboratories, the B. mori strains of Dazao and p50T were first separately sequenced and annotated in China and Japan by whole-genome shotgun sequencing in 2004, and the genome sequence has been updated every few years, which is considered the gold standard (Mita et al., 2004; Xia et al., 2004; Meng et al., 2017; Kawamoto et al., 2019).

According to a search of NCBI PubMed with the keywords “(Bombyx mori) OR (silkworm),” 7,024 relevant works have been published since 2004. Among these, transcriptomics studies have inevitably made use of the above-mentioned reference genomes to identify transcripts and quantify gene expression. To obtain the brain transcriptome profiles of BmNPV-infected and non-infected silkworm larvae, Wang et al. (2015) mapped their clean reads to the SilkDB reference genome. Moreover, Li et al. (2016) comprehensively investigated the transcriptomic changes between susceptible and resistant B. mori strains after BmNPV infection. Three years later, they (Li et al., 2019) updated the analysis of the different alternative splicing events based on the same reference genome. Recently, Sun Q. et al. (2020) also downloaded genome sequences and annotation files from the SilkDB website for the transcriptome analysis of the immune response of B. mori after bidensovirus infection at the early stage. Based on the SilkBase genome sequence, Shoji et al. (2013) characterized a novel chromodomain-containing gene and later (Shoji et al., 2014) identified a number of genes whose expression can be enhanced by heterochromatin protein 1. Therefore, it follows that the two reference genomes with non-resistance to BmNPV are being widely used as the B. mori “Bible.” However, with the increasing demand for precise reference genomes, a personalized edition is needed for B. mori research. With our more refined reference genome, one can find precise bioinformatics information such as SNPs and INDELs related to BmNPV resistance.

Thanks to the rapid advances in sequencing technology and bioinformatic tools in the past decade, transcriptome sequencing (∼$130/sample) has become a cost-effective technology for obtaining biological information. However, the emerging third-generation sequencing technologies, such as Nanopore/PacBio long-read sequencing and Hi-C sequencing, require a larger budget than that of standard laboratories (Makałowski and Shabardina, 2019; Wang et al., 2019). Even the use of whole-genome sequencing on the popular Illumina platform to generate short reads for draft assembly correction is expensive, let alone optical genome mapping (i.e., Bionano maps), which can correct chimeric contigs and scaffolding errors caused by Hi-C sequencing (Chan et al., 2018; Bocklandt et al., 2019; Choi et al., 2020; Ning et al., 2020; Sun L. et al., 2020). In addition, the entire pipeline of genome assembly is a tedious and time-consuming process, as ideal computing modules should be constructed for each individual species. As such, methods and algorithms have been developed or improved by scientists with different professional backgrounds (Shin et al., 2019; Nakabayashi and Morishita, 2020). Thus, there is a need for many users to find a single universal pipeline for datasets from different species for more accurate results (Minervini et al., 2020). In short, for the further analysis of species, an updated reference genome with a high resolution is needed to enhance the robustness of the final results.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  ArunkumarK. P.MettaM.NagarajuJ. (2006). Molecular phylogeny of silkmoths reveals the origin of domesticated silkmoth, Bombyx mori from Chinese Bombyx mandarina and paternal inheritance of Antheraea proylei mitochondrial DNA.Mol. Phylogenet. Evol.40419–427. 10.1016/j.ympev.2006.02.023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  BelserC.IstaceB.DenisE.DubarryM.BaurensF. C.FalentinC.et al (2018). Chromosome-scale assemblies of plant genomes using nanopore long reads and optical maps.Nat. Plants4879–887. 10.1038/s41477-018-0289-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BeltonJ. M.McCordR. P.GibcusJ. H.NaumovaN.ZhanY.DekkerJ. (2012). Hi-C: a comprehensive technique to capture the conformation of genomes.Methods58268–276. 10.1016/j.ymeth.2012.05.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BensonG. (1999). Tandem repeats finder- a program to analyze DNA sequences.Nucleic Acids Res.27573–580. 10.1093/nar/27.2.573

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BocklandtS.HastieA.CaoH. (2019). Bionano genome mapping: high-throughput, ultra-long molecule genome analysis system for precision genome assembly and haploid-resolved structural variation discovery.Adv. Exp. Med. Biol.112997–118. 10.1007/978-981-13-6037-4_7

  • CrossRef
  • Google Scholar

• 6

  BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data.Bioinformatics302114–2120. 10.1093/bioinformatics/btu170

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  BurgeC.KarlinS. (1997). Prediction of complete gene structures in human genomic DNA.J. Mol. Biol.26878–94. 10.1006/jmbi.1997.0951

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  BurtonJ. N.AdeyA.PatwardhanR. P.QiuR.KitzmanJ. O.ShendureJ. (2013). Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions.Nat. Biotechnol.311119–1125. 10.1038/nbt.2727

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  CaoX.HuangY.XiaD.QiuZ.ShenX.GuoX.et al (2017). BmNPV-miR-415 up-regulates the expression of TOR2 via Bmo-miR-5738.Saudi J. Biol. Sci.241614–1619. 10.1016/j.sjbs.2015.09.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  CastresanaJ. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis.Mol. Biol. Evol.17540–552. 10.1093/oxfordjournals.molbev.a026334

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  ChanP. P.LoweT. M. (2019). tRNAscan-SE: searching for tRNA genes in genomic sequences.Methods Mol. Biol.19621–14. 10.1007/978-1-4939-9173-0_1

  • CrossRef
  • Google Scholar

• 12

  ChanS.LamE.SaghbiniM.BocklandtS.HastieA.CaoH.et al (2018). Structural variation detection and analysis using bionano optical mapping.Methods Mol. Biol.1833193–203. 10.1007/978-1-4939-8666-8_16

  • CrossRef
  • Google Scholar

• 13

  ChenF.MackeyA. J.StoeckertC. J.Jr.RoosD. S. (2006). OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups.Nucleic Acids Res.34D363–D368. 10.1093/nar/gkj123

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  ChenH.RangasamyM.TanS. Y.WangH.SiegfriedB. D. (2010). Evaluation of five methods for total DNA extraction from western corn rootworm beetles.PLoS One5:e11963. 10.1371/journal.pone.0011963

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  ChinC. S.AlexanderD. H.MarksP.KlammerA. A.DrakeJ.HeinerC.et al (2013). Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data.Nat. Methods10563–569. 10.1038/nmeth.2474

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  ChinC. S.PelusoP.SedlazeckF. J.NattestadM.ConcepcionG. T.ClumA.et al (2016). Phased diploid genome assembly with single-molecule real-time sequencing.Nat. Methods131050–1054. 10.1038/nmeth.4035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  ChoiJ. H.LeeJ. H.ChoiJ. W. (2020). Applications of bionano sensor for extracellular vesicles analysis.Materials13:3677. 10.3390/ma13173677

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  ConesaA.MadrigalP.TarazonaS.Gomez-CabreroD.CerveraA.McPhersonA.et al (2016). A survey of best practices for RNA-seq data analysis.Genome Biol.17:13. 10.1186/s13059-016-0881-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  DeschampsS.ZhangY.LlacaV.YeL.SanyalA.KingM.et al (2018). A chromosome-scale assembly of the sorghum genome using nanopore sequencing and optical mapping.Nat. Commun.9:4844. 10.1038/s41467-018-07271-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  DuanJ.LiR.ChengD.FanW.ZhaX.ChengT.et al (2010). SilkDB v2.0: a platform for silkworm (Bombyx mori) genome biology.Nucleic Acids Res.38D453–D456. 10.1093/nar/gkp801

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  EdgarR. C. (2004a). MUSCLE: a multiple sequence alignment method with reduced time and space complexity.BMC Bioinformatics5:113. 10.1186/1471-2105-5-113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  EdgarR. C. (2004b). MUSCLE: multiple sequence alignment with high accuracy and high throughput.Nucleic Acids Res.321792–1797. 10.1093/nar/gkh340

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  EkblomR.WolfJ. B. (2014). A field guide to whole-genome sequencing, assembly and annotation.Evol. Appl.71026–1042. 10.1111/eva.12178

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  FlynnJ. M.HubleyR.GoubertC.RosenJ.ClarkA. G.FeschotteC.et al (2020). RepeatModeler2 for automated genomic discovery of transposable element families.Proc. Natl. Acad. Sci. U.S.A.1179451–9457. 10.1073/pnas.1921046117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  FujimotoS.KawamotoM.ShojiK.SuzukiY.KatsumaS.IwanagaM. (2020). Whole-genome sequencing and comparative transcriptome analysis of Bombyx mori nucleopolyhedrovirus La strain.Virus Genes56249–259. 10.1007/s11262-019-01727-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  GeQ.ChenL.TangM.ZhangS.LiuL.GaoL.et al (2018a). Analysis of mulberry leaf components in the treatment of diabetes using network pharmacology.Eur. J. Pharmacol.83350–62. 10.1016/j.ejphar.2018.05.021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  GeQ.ZhangS.ChenL.TangM.LiuL.KangM.et al (2018b). Mulberry leaf regulates differentially expressed genes in diabetic mice liver based on RNA-Seq analysis.Front. Physiol.9:1051. 10.3389/fphys.2018.01051

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  GertzE. M.YuY. K.AgarwalaR.SchafferA. A.AltschulS. F. (2006). Composition-based statistics and translated nucleotide searches: improving the TBLASTN module of BLAST.BMC Biol.4:41. 10.1186/1741-7007-4-41

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  GhoshS.ChanC. K. (2016). Analysis of RNA-seq data using tophat and cufflinks.Methods Mol. Biol.1374339–361. 10.1007/978-1-4939-3167-5_18

  • CrossRef
  • Google Scholar

• 30

  GrabherrM. G.HaasB. J.YassourM.LevinJ. Z.ThompsonD. A.AmitI.et al (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome.Nat. Biotechnol.29644–652. 10.1038/nbt.1883

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  HaasB. J.DelcherA. L.MountS. M.WortmanJ. R.SmithR. K.Jr.HannickL. I.et al (2003). Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies.Nucleic Acids Res.315654–5666. 10.1093/nar/gkg770

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HaasB. J.PapanicolaouA.YassourM.GrabherrM.BloodP. D.BowdenJ.et al (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis.Nat. Protoc.81494–1512. 10.1038/nprot.2013.084

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  HaasB. J.SalzbergS. L.ZhuW.PerteaM.AllenJ. E.OrvisJ.et al (2008). Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments.Genome Biol.9:R7. 10.1186/gb-2008-9-1-r7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  International Silkworm Genome. (2008). The genome of a lepidopteran model insect, the silkworm Bombyx mori.Insect Biochem. Mol. Biol.381036–1045. 10.1016/j.ibmb.2008.11.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  JayakumarV.SakakibaraY. (2019). Comprehensive evaluation of non-hybrid genome assembly tools for third-generation PacBio long-read sequence data.Brief Bioinform.20866–876. 10.1093/bib/bbx147

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  JiangL.XiaQ. (2014). The progress and future of enhancing antiviral capacity by transgenic technology in the silkworm Bombyx mori.Insect Biochem. Mol. Biol.481–7. 10.1016/j.ibmb.2014.02.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  JohnsonA. D.HandsakerR. E.PulitS. L.NizzariM. M.O’DonnellC. J.de BakkerP. I. (2008). SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap.Bioinformatics242938–2939. 10.1093/bioinformatics/btn564

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  KaitoC.AkimitsuN.WatanabeH.SekimizuK. (2002). Silkworm larvae as an animal model of bacterial infection pathogenic to humans.Microb. Pathog.32183–190. 10.1006/mpat.2002.0494

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  KalvariI.NawrockiE. P.ArgasinskaJ.Quinones-OlveraN.FinnR. D.BatemanA.et al (2018). Non-coding RNA analysis using the Rfam database.Curr. Protoc. Bioinformatics62:e51. 10.1002/cpbi.51

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  KalvariI.NawrockiE. P.Ontiveros-PalaciosN.ArgasinskaJ.LamkiewiczK.MarzM.et al (2021). Rfam 14: expanded coverage of metagenomic, viral and microRNA families.Nucleic Acids Res.49D192–D200. 10.1093/nar/gkaa1047

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  KasajimaI. (2018). Successful tips of DNA extraction and PCR of plants for beginners.Trends Res.11–2. 10.15761/tr.1000115

  • CrossRef
  • Google Scholar

• 42

  KawamotoM.JourakuA.ToyodaA.YokoiK.MinakuchiY.KatsumaS.et al (2019). High-quality genome assembly of the silkworm, Bombyx mori.Insect Biochem. Mol. Biol.10753–62. 10.1016/j.ibmb.2019.02.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  KentW. J. (2002). BLAT–the BLAST-like alignment tool.Genome Res.12656–664. 10.1101/gr.229202

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  KimD.PerteaG.TrapnellC.PimentelH.KelleyR.SalzbergS. L. (2013). TopHat2- accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.Genome Biol.14:R36.

  • Google Scholar

• 45

  KinganS. B.HeatonH.CudiniJ.LambertC. C.BaybayanP.GalvinB. D.et al (2019). A high-quality de novo genome assembly from a single mosquito using PacBio sequencing.Genes10:62. 10.3390/genes10010062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  KoikeY.MitaK.SuzukiM. G.MaedaS.AbeH.OsoegawaK.et al (2003). Genomic sequence of a 320-kb segment of the Z chromosome of Bombyx mori containing a kettin ortholog.Mol. Genet. Genomics269137–149. 10.1007/s00438-003-0822-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  LiG.QianH.LuoX.XuP.YangJ.LiuM.et al (2016). Transcriptomic analysis of resistant and susceptible Bombyx mori strains following BmNPV infection provides insights into the antiviral mechanisms.Int. J. Genomics2016:2086346. 10.1155/2016/2086346

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  LiG.ZhouK.ZhaoG.QianH.XuA. (2019). Transcriptome-wide analysis of the difference of alternative splicing in susceptible and resistant silkworm strains after BmNPV infection.3 Biotech9:152. 10.1007/s13205-019-1669-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  LiH.DurbinR. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform.Bioinformatics26589–595. 10.1093/bioinformatics/btp698

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al (2009). The sequence alignment/map format and SAMtools.Bioinformatics252078–2079. 10.1093/bioinformatics/btp352

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  LiL.StoeckertC. J.Jr.RoosD. S. (2003). OrthoMCL: identification of ortholog groups for eukaryotic genomes.Genome Res.132178–2189. 10.1101/gr.1224503

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  LoweT. M.ChanP. P. (2016). tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes.Nucleic Acids Res.44W54–W57. 10.1093/nar/gkw413

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  LuF.WeiZ.LuoY.GuoH.ZhangG.XiaQ.et al (2020). SilkDB 3.0: visualizing and exploring multiple levels of data for silkworm.Nucleic Acids Res.48D749–D755. 10.1093/nar/gkz919

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  MaekawaH.TakadaN.MikitaniK.OguraT.MiyajimaN.FujiwaraH.et al (1988). Nucleolus organizers in the wild silkworm Bombyx mandarina and the domesticated silkworm B. mori.Chromosoma96263–269. 10.1007/bf00286912

  • CrossRef
  • Google Scholar

• 55

  MajorosW. H.PerteaM.SalzbergS. L. (2004). TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders.Bioinformatics202878–2879. 10.1093/bioinformatics/bth315

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  MakałowskiW.ShabardinaV. (2019). Bioinformatics of nanopore sequencing.J. Hum. Genet.6561–67. 10.1038/s10038-019-0659-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  MengX.ZhuF.ChenK. (2017). Silkworm: a promising model organism in life science.J. Insect Sci.17:97. 10.1093/jisesa/iex064

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  MinerviniC. F.CumboC.OrsiniP.AnelliL.ZagariaA.SpecchiaG.et al (2020). Nanopore sequencing in blood diseases: a wide range of opportunities.Front. Genet.11:76. 10.3389/fgene.2020.00076

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  MitaK.KasaharaM.SasakiS.NagayasuY.YamadaT.KanamoriH.et al (2004). The genome sequence of silkworm, Bombyx mori.DNA Res.1127–35. 10.1093/dnares/11.1.27

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  NakabayashiR.MorishitaS. (2020). HiC-Hiker: a probabilistic model to determine contig orientation in chromosome-length scaffolds with Hi-C.Bioinformatics363966–3974. 10.1093/bioinformatics/btaa288

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  NingD. L.WuT.XiaoL. J.MaT.FangW. L.DongR. Q.et al (2020). Chromosomal-level assembly of Juglans sigillata genome using Nanopore, BioNano, and Hi-C analysis.Gigascience9:giaa006. 10.1093/gigascience/giaa006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  OluwadareO.HighsmithM.ChengJ. (2019). An overview of methods for reconstructing 3-D chromosome and genome structures from Hi-C data.Biol. Proced. Online21:7. 10.1186/s12575-019-0094-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  ParraG.BradnamK.KorfI. (2007). CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes.Bioinformatics231061–1067. 10.1093/bioinformatics/btm071

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  RangF. J.KloostermanW. P.de RidderJ. (2018). From squiggle to basepair: computational approaches for improving nanopore sequencing read accuracy.Genome Biol.19:90. 10.1186/s13059-018-1462-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  SeppeyM.ManniM.ZdobnovE. M. (2019). BUSCO: assessing genome assembly and annotation completeness.Methods Mol. Biol.1962227–245. 10.1007/978-1-4939-9173-0_14

  • CrossRef
  • Google Scholar

• 66

  ShinS. C.KimH.LeeJ. H.KimH. W.ParkJ.ChoiB. S.et al (2019). Nanopore sequencing reads improve assembly and gene annotation of the Parochlus steinenii genome.Sci. Rep.9:5095. 10.1038/s41598-019-41549-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  ShojiK.HaraK.KawamotoM.KiuchiT.KawaokaS.SuganoS.et al (2014). Silkworm HP1a transcriptionally enhances highly expressed euchromatic genes via association with their transcription start sites.Nucleic Acids Res.4211462–11471. 10.1093/nar/gku862

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  ShojiK.KiuchiT.HaraK.KawamotoM.KawaokaS.ArimuraS.et al (2013). Characterization of a novel chromodomain-containing gene from the silkworm, Bombyx mori.Gene527649–654. 10.1016/j.gene.2013.06.071

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  StamatakisA. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.Bioinformatics301312–1313. 10.1093/bioinformatics/btu033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  StankeM.MorgensternB. (2005). AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints.Nucleic Acids Res.33W465–W467. 10.1093/nar/gki458

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  SunL.GaoT.WangF.QinZ.YanL.TaoW.et al (2020). Chromosome-level genome assembly of a cyprinid fish Onychostoma macrolepis by integration of nanopore sequencing, Bionano and Hi-C technology.Mol. Ecol. Resour.201361–1371. 10.1111/1755-0998.13190

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  SunQ.GuoH.XiaQ.JiangL.ZhaoP. (2020). Transcriptome analysis of the immune response of silkworm at the early stage of Bombyx mori bidensovirus infection.Dev. Comp. Immunol.106:103601. 10.1016/j.dci.2019.103601

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  TabunokiH.BonoH.ItoK.YokoyamaT. (2016). Can the silkworm (Bombyx mori) be used as a human disease model?Drug Discov. Ther.103–8. 10.5582/ddt.2016.01011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  TangH.BowersJ. E.WangX.MingR.AlamM.PatersonA. H. (2008). Synteny and collinearity in plant genomes.Science320486–488. 10.1126/science.1153917

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  TomitaM.MunetsunaH.SatoT.AdachiT.HinoR.HayashiM.et al (2000). Germline transformation of the silkworm Bombyx mori L. using a piggyBac transposon-derived vector.Nat. Biotechnol.1881–84. 10.1038/71978

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  TomitaM.MunetsunaH.SatoT.AdachiT.HinoR.HayashiM.et al (2003). Transgenic silkworms produce recombinant human type III procollagen in cocoons.Nat. Biotechnol.2152–56. 10.1038/nbt771

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  VurtureG. W.SedlazeckF. J.NattestadM.UnderwoodC. J.FangH.GurtowskiJ.et al (2017). GenomeScope: fast reference-free genome profiling from short reads.Bioinformatics332202–2204. 10.1093/bioinformatics/btx153

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  WalkerB. J.AbeelT.SheaT.PriestM.AbouellielA.SakthikumarS.et al (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement.PLoS One9:e112963. 10.1371/journal.pone.0112963

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  WangG.ZhangJ.ShenY.ZhengQ.FengM.XiangX.et al (2015). Transcriptome analysis of the brain of the silkworm Bombyx mori infected with Bombyx mori nucleopolyhedrovirus: a new insight into the molecular mechanism of enhanced locomotor activity induced by viral infection.J. Invertebr. Pathol.12837–43. 10.1016/j.jip.2015.04.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  WangJ.XiaQ.HeX.DaiM.RuanJ.ChenJ.et al (2005). SilkDB: a knowledgebase for silkworm biology and genomics.Nucleic Acids Res.33D399–D402. 10.1093/nar/gki116

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  WangJ.YangJ.YingY. L.LongY. T. (2019). Nanopore-based confined spaces for single-molecular analysis.Chem. Asian J.14389–397. 10.1002/asia.201801648

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  WheelerT. J.EddyS. R. (2013). nhmmer: DNA homology search with profile HMMs.Bioinformatics292487–2489. 10.1093/bioinformatics/btt403

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  XiaQ.ZhouZ.LuC.ChengD.DaiF.LiB.et al (2004). A draft sequence for the genome of the domesticated silkworm (Bombyx mori).Science3061937–1940. 10.1126/science.1102210

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  XiangH.LiuX.LiM.ZhuY.WangL.CuiY.et al (2018). The evolutionary road from wild moth to domestic silkworm.Nat. Ecol. Evol.21268–1279. 10.1038/s41559-018-0593-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  XuY. P.ChengR. L.XiY.ZhangC. X. (2013). Genomic diversity of Bombyx mori nucleopolyhedrovirus strains.Genomics10263–71. 10.1016/j.ygeno.2013.04.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  XuZ.WangH. (2007). LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons.Nucleic Acids Res.35W265–W268. 10.1093/nar/gkm286

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  YangZ. (1997). PAML: a program package for phylogenetic analysis by maximum likelihood.Comput. Appl. Biosci.13555–556. 10.1093/bioinformatics/13.5.555

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  YangZ.NielsenR. (2002). Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages.Mol. Biol. Evol.19908–917. 10.1093/oxfordjournals.molbev.a004148

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  YardimciG. G.OzadamH.SauriaM. E. G.UrsuO.YanK. K.YangT.et al (2019). Measuring the reproducibility and quality of Hi-C data.Genome Biol.20:57. 10.1186/s13059-019-1658-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  ZhangJ.NielsenR.YangZ. (2005). Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level.Mol. Biol. Evol.222472–2479. 10.1093/molbev/msi237

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91

  ZhuZ.GuanZ.LiuG.WangY.ZhangZ. (2019). SGID: a comprehensive and interactive database of the silkworm.Database2019:baz134. 10.1093/database/baz134

  • Pubmed Abstract
  • CrossRef
  • Google Scholar