# De novo genome and transcriptome resources of the Adzuki bean borer Ostrinia scapulalis (Lepidoptera: Crambidae)

Bernhard Gschloessl, Franck Dorkeld, Philippe Audiot, Anthony Bretaudieu, Carole Kerdelhue, Rejane Streiff

► To cite this version:

Bernhard Gschloessl, Franck Dorkeld, Philippe Audiot, Anthony Bretaudieu, Carole Kerdelhue, et al.. De novo genome and transcriptome resources of the Adzuki bean borer Ostrinia scapulalis (Lepidoptera: Crambidae). Data in Brief, 2018, 17, pp.781-787. 10.1016/j.dib.2018.01.073. hal-02627842

HAL Id: hal-02627842

https://hal.inrae.fr/hal-02627842

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Data Article

*De novo* genome and transcriptome resources
of the Adzuki bean borer *Ostrinia scapulalis*
(Lepidoptera: Crambidae)

<figure><img src="image_3.png" /></figure>

Check for
updates

B. Gschloessl<sup>a,*</sup>, F. Dorkeld<sup>a</sup>, P. Audiot<sup>a</sup>, A. Bretaudieu<sup>b,c</sup>, C. Kerdelhué<sup>a</sup>, R. Streiff<sup>a</sup>

ARTICLE INFO

Article history:

Received 23 December 2017

Received in revised form

23 January 2018

Accepted 25 January 2018

Available online 1 February 2018

Keywords:

De novo assembly

Genome

Transcriptome

Gene prediction

Lepidoptera

Crambidae

A B S T R A C T

We present a draft genome assembly with a *de novo* prediction and automated functional annotation of coding genes, and a reference transcriptome of the Adzuki bean borer, *Ostrinia scapulalis*, based on RNA sequencing of various tissues and developmental stages. The genome assembly spans 419 Mb, has a GC content of 37.4% and includes 26,120 predicted coding genes. The reference transcriptome holds 33,080 unigenes and contains a high proportion of a set of genes conserved in eukaryotes and arthropods, used as quality assessment of the reconstructed transcripts. The new genomic and transcriptomic data presented here significantly enrich the public sequence databases for the Crambidae and Lepidoptera, and represent useful resources for future researches related to the evolution and the adaptation of phytophagous moths. The genome and transcriptome assemblies have been deposited and made accessible *via* a NCBI BioProject (id PRJNA390510) and the LepidoDB database ([http://bipaa.genouest.org/sp/ostrinia_scapulalis/](http://bipaa.genouest.org/sp/ostrinia_scapulalis/)).

© 2018 The Authors. Published by Elsevier Inc. This is an open

access article under the CC BY license

([http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)).

<footnote><sup>a</sup> CBGP, INRA, CIRAD, IRD, Montpellier SupAgro, Univ Montpellier, Montpellier, France</footnote>

<footnote><sup>b</sup> INRA, UMR Institut de Génétique, Environnement et Protection des Plantes (IGEPP), BioInformatics Platform for Agroecosystems Arthropods (BIPAA), Campus Beaulieu, Rennes, France</footnote>

<footnote><sup>c</sup> INRIA, IRISA, GenOuest Core Facility, Campus de Beaulieu, Rennes, France</footnote>

<footnote>* Correspondence to: CBGP, INRA, CIRAD, IRD, Montpellier SupAgro, Univ., 755 avenue du Campus Agropolis, CS30016, F-34988 Montferrier-sur-Lez cedex, France. Fax: +33 4 99 62 33 45.</footnote>

<footnote>E-mail address: Bernhard.Gschloessl@inra.fr (B. Gschloessl).</footnote>
**Specifications Table** Genome and Transcriptome

<table><tr><td>Subject area</td><td>Biology</td></tr><tr><td>More specific subject area</td><td>Lepidoptera, Genomics</td></tr><tr><td>Type of data</td><td>DNA and cDNA sequence reads, genome assembly and transcript assembly</td></tr><tr><td>How data was acquired</td><td>Shotgun whole genome and cDNA sequencing using Illumina HiSeq 2000</td></tr><tr><td>Data format</td><td>Analyzed: i.e. raw data and assembled sequences</td></tr><tr><td>Experimental factors</td><td>Genome: total DNA extraction from male larvae of wild samples<br/>Transcriptome: total RNA extraction from various tissues, developmental stages and of males and females</td></tr><tr><td>Experimental features</td><td>Genome: DNA sequencing<br/>Transcriptome: RNA sequencing of various tissues from eggs to adults, controlled conditions</td></tr><tr><td>Data source location</td><td>Genome: Amiens, Picardie/France (49°54'0.01"N, 2°18'0"E)<br/>Transcriptome: Nadarzin, Poland (52°4'2.05"N, 20°47'33.00"E)</td></tr><tr><td>Data accessibility</td><td>All raw sequence reads are accessible as NCBI BioProject (id PRJNA390510).<br/>The OSCA v1.2 draft genome assembly, the reference transcriptome assembly and automatic functional annotations can be found in the LepidoDB database (http://bipaa.genouest.org/sp/ostrinia_scapulalis/).</td></tr></table>

**Value of the data**

* The draft genome represents the first available genome assembly for *O. scapulalis*.

* The reference transcriptome of *O. scapulalis* will allow comparative expression studies.

* The new genomic and transcriptomic data enrich the public sequence databases for the Crambidae and Lepidoptera.

* The data represent pangenomic resources for future researches related to the evolution and the adaptation of phytophagous moths.

**1. Data**

The Adzuki bean borer, *Ostrinia scapulalis* (hereafter OSCA), is a palaearctic phytophagous moth feeding on various dicotyledons, including hop (*Humulus lupulus*), mugwort (*Artemisia vulgaris*) and hemp (*Cannabis sativa*) [1]. In Europe, it partly co-occurs with its sibling species, the European corn borer, *Ostrinia nubilalis*, which is a major pest of maize (*Zea mays*). Previous studies demonstrated that *O. scapulalis* and *O. nubilalis* are specialized to their respective host plants [2–7] and that their genetic divergence is rather low so that they can be considered as sibling species [8]. Yet, a few genomic sequences and rearrangements are much more divergent than the rest of the genomic background [1,9–11]. These genomic regions are of particular interest to understand the divergence process between *O. scapulalis* and *O. nubilalis*. To further investigate the host adaptation and divergence between these two sibling species at a pangenomic scale, we have elaborated new genomic and transcriptomic resources consisting of an OSCA draft genome and a related reference transcriptome. The latter extends a published transcriptomic set generated with Roche 454 sequencing technology [12].
<table><caption>**Table 1**

Statistic features of sequence reads issued from different genomic libraries and used for the OSCA v1.2 genome assembly.</caption><thead><tr><th>Library</th><th>PE300</th><th>Mate3kb</th><th>Mate8kb</th></tr></thead><tbody><tr><td>Raw read count</td><td>350,636,628</td><td>234,843,432</td><td>295,988,124</td></tr><tr><td>Read count after clipping<sup>a</sup></td><td>349,359,410</td><td>233,657,586</td><td>287,067,720</td></tr><tr><td>Total Size [Gb]</td><td>34.5</td><td>11.9</td><td>23.9</td></tr><tr><td>Minimum Length [bp]</td><td>20</td><td>20</td><td>20</td></tr><tr><td>Maximum Length [bp]</td><td>100</td><td>51</td><td>100</td></tr><tr><td>Mean Length [bp]</td><td>99</td><td>51</td><td>83</td></tr><tr><td>Expected insert size [bp]<sup>b</sup></td><td>300</td><td>2500</td><td>7500</td></tr><tr><td>Mapped on OSCA v1.2 genome</td><td>207,000,258</td><td>101,565,815</td><td>146,139,301</td></tr><tr><td>Mapped as paired-end on OSCA v1.2 genome</td><td>170,567,848</td><td>16,877,194</td><td>31,116,832</td></tr></tbody></table>

<table><caption>**Table 2**

Genome features of the OSCA v1.2 assembly. The coverage is defined as the average read count per assembled bp.</caption><thead><tr><th></th><th>OSCA v1.2 assembly</th></tr></thead><tbody><tr><td>Contig count</td><td>163,703</td></tr><tr><td>Scaffold (scf) count</td><td>50,738</td></tr><tr><td>N50 scf [bp]</td><td>29,308</td></tr><tr><td>N50 scf sequence count</td><td>3395</td></tr><tr><td>N90 scf [bp]</td><td>3051</td></tr><tr><td>N90 scf sequence count</td><td>21,872</td></tr><tr><td>Minimum scf length [bp]</td><td>883</td></tr><tr><td>Maximum scf length [kb]</td><td>619.8</td></tr><tr><td>Illumina PE300 coverage [reads/bp]</td><td>49.7</td></tr><tr><td>GC content [%]</td><td>37.4</td></tr><tr><td>N base content in assembly [%]</td><td>30.9</td></tr><tr><td>Total Length [Mb]</td><td>419.2</td></tr><tr><td>Total Length without Ns [Mb]</td><td>289.8</td></tr><tr><td>GenomeScope PE genome size estimation [Mb]</td><td>302.9</td></tr><tr><td>CEGMA identified [%] (count of 248)</td><td>152 (61.3)</td></tr><tr><td>CEGMA full-length [%] (count of 248)</td><td>83 (33.5)</td></tr><tr><td>BUSCO2 euk identified [%] (count of 303)</td><td>211 (69.6)</td></tr><tr><td>BUSCO2 euk full-length [%] (count of 303)</td><td>156 (51.5)</td></tr><tr><td>BUSCO2 arthropod identified [%] (count of 2675)</td><td>1363 (50.9)</td></tr><tr><td>BUSCO2 arthropod full-length [%] (count of 2675)</td><td>842 (31.5)</td></tr></tbody></table>

## 2. Experimental design, materials and methods

### 2.1. De novo draft genome

Diapausing larvae were collected in mugwort stems in 2008 near Amiens (Picardie, France) and stored in 95% ethanol at –20 °C. Whole genomic DNA extracts were obtained from a CTAB-based method [13]. DNA quality and integrity was evaluated through migration on an agarose gel and nanodrop technology. The sex of each sampled larvae was determined with a molecular coamplification of markers specific to each heterochromosome (Z and W in Lepidoptera) as described in Orsucci et al. [6]. Only samples of the ZZ homogametic sex (males in Lepidoptera), were retained for the libraries construction. A 2 × 100 bp shot-gun paired-end library and a 3 (2 × 50 bp) and an 8 kb (2 × 100 bp) mate-pair library were generated using the DNA extract of one larva for each library and the Illumina TruSeq TM and Nextera Mate Pair Library Preparation kits, respectively. All libraries were sequenced by LGC Genomics GmbH (Berlin, Germany) on an Illumina HiSeq 2000 platform using the paired-end protocol. Between 234 and 351 million DNA raw reads were generated per library (Table 1). Assembly and scaffolding of the cleaned reads were done with the software Allpaths-LG

<footnote><sup>a</sup> Only paired-end reads taken into account</footnote>

<footnote><sup>b</sup> Based on *in silico* measurements</footnote>
<table><caption>**Table 3**
Number of repeated elements found in the OSCA v1.2 draft genome assembly and corresponding genome ratio.</caption><thead><tr><th>Family</th><th>Fragments</th><th>Total length [Mb]</th><th>% of genome</th></tr></thead><tbody><tr><td>LTR</td><td>413,393</td><td>45.8</td><td>10.9</td></tr><tr><td>LINE</td><td>95,506</td><td>11.8</td><td>2.8</td></tr><tr><td>SINE</td><td>34,673</td><td>3.8</td><td>0.9</td></tr><tr><td>DNA</td><td>64,646</td><td>8</td><td>1.9</td></tr><tr><td>Total</td><td>608,218</td><td>69.4</td><td>16.6</td></tr></tbody></table>

<table><caption>**Table 4**
Characteristics of the different transcriptome assemblies and the genes predicted from the genome.</caption><thead><tr><th></th><th>HiSeq transcriptome</th><th>454 transcriptome</th><th>MAKER genes</th></tr></thead><tbody><tr><td>Raw read count</td><td>325,008,948</td><td>322,504</td><td>N/A</td></tr><tr><td>Cleaned read count</td><td>267,359,188</td><td>287,429</td><td>N/A</td></tr><tr><td>Mapped Reads</td><td>198,962,467</td><td>145,588</td><td>N/A</td></tr><tr><td>Transcriptome size [Mb]</td><td>49.2</td><td>10.4</td><td>21.6</td></tr><tr><td>Coverage (mean read count per bp)</td><td>339.2</td><td>9.8</td><td>N/A</td></tr><tr><td>Transcript count</td><td>44,564</td><td>11,231</td><td>26,120</td></tr><tr><td>Unigene count</td><td>33,080</td><td>8892</td><td>26,120</td></tr><tr><td>Mean CDS/transcript length [bp]</td><td>1103</td><td>922</td><td>829</td></tr><tr><td>Median transcript length [bp]</td><td>591</td><td>693</td><td>498</td></tr><tr><td>N50 transcript length [bp]</td><td>2006</td><td>1036</td><td>1296</td></tr><tr><td>N50 sequence count</td><td>7061</td><td>2858</td><td>4475</td></tr><tr><td>N90 transcript length [bp]</td><td>429</td><td>485</td><td>353</td></tr><tr><td>N90 sequence count</td><td>26,722</td><td>8736</td><td>17,539</td></tr><tr><td>Minimum length</td><td>201</td><td>111</td><td>66</td></tr><tr><td>Maximum length</td><td>27,559</td><td>10,991</td><td>53,685</td></tr><tr><td>CEGMA identified [%] (count of 248)</td><td>209 (84.3)</td><td>56 (22.6)</td><td>N/A</td></tr><tr><td>CEGMA full-length [%] (count of 248)</td><td>189 (76.2)</td><td>22 (8.9)</td><td>N/A</td></tr><tr><td>BUSCO2 euk identified [%] (count of 303)</td><td>268 (88.4)</td><td>70 (23.1)</td><td>226 (74.6)</td></tr><tr><td>BUSCO2 euk full-length [%] (count of 303)</td><td>256 (84.5)</td><td>17 (5.6)</td><td>162 (53.5)</td></tr><tr><td>BUSCO2 arthropod identified [%] (count of 2675)</td><td>1892 (70.7)</td><td>265 (9.9)</td><td>1453 (54.3)</td></tr><tr><td>BUSCO2 arthropod full-length [%] (count of 2675)</td><td>2109 (78.8)</td><td>119 (4.4)</td><td>923 (34.5)</td></tr><tr><td>Transcripts with predicted CDS (%)</td><td>18,494 (41.5)</td><td>4016 (35.8)</td><td>26,120 (100)</td></tr><tr><td>Transcripts with full-length CDS (%)</td><td>12,515 (28.1)</td><td>1826 (16.3)</td><td>26,120 (100)</td></tr><tr><td>Located on OSCA v1.2 genome [count] (%)</td><td>11,010 (24.7)</td><td>3331 (29.7)</td><td>26,120 (100)</td></tr><tr><td>Split on two OSCA v1.2 scaffolds [count] (%)</td><td>11,721 (26.3)</td><td>2350 (20.9)</td><td>0 (0)</td></tr></tbody></table>

[14], followed by GapCloser for closing gaps in the assembled scaffolds. The resulting scaffolding was then further improved by integrating independent reconstructed transcript data after RNA sequencing (see Supplementary material for details). The OSCA v1.2 draft genome assembly consisted of 50,738 scaffolds, representing 419 Mb with a mean read coverage of 50 reads per base (Table 2). N50 and N90 were 29,308 and 3051 bp, respectively. Scaffold lengths ranged between 883 bp and 619.8 kb. The genome assembly had a GC content of 37.4% and a proportion of repeated elements of 16.6% (Table 3). A total of 8372 short duplicated regions (average length: 1813 bp, min.: 808 bp, max.: 11,745 bp) were identified on 7009 scaffolds. The quality of the OSCA v1.2 genome assembly was evaluated by the recovery rate of three sets of genes highly conserved in eukaryotes and arthropods [15,16]. Among both eukaryotic gene sets, 61% of the CEGMA genes and 70% of the BUSCO genes were recovered in the OSCA v1.2 genome, while 51% of the BUSCO arthropod conserved genes could be identified. Using the MAKER pipeline [see Supplementary material, [17]] on the OSCA v1.2 nuclear genome 26,120 coding genes were predicted (Table 4). Of these coding genes, 80.3% could be functionally annotated. Furthermore, 19,023 OSCA genes were assigned to 9785 ortholog groups (Fig. 1) of which 93% were shared with at least one of the three Lepidoptera species *Bombyx mori*, *Danaus plexippus* or *Spodoptera frugiperda*.
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Venn diagram showing all OrthoMCL ortholog groups among the MAKER-predicted *O. scapulalis* proteins (OSCA) and the proteomes of *Spodoptera frugiperda* (SFRU), *Bombyx mori* (BMOR), *Danaus plexippus* (DPLE) and *Drosophila melanogaster* (DMEL).</figcaption></figure>

<caption>**Table 5**</caption>

Developmental stages and tissues of 7 RNA extracts issued from F1 individuals obtained in the laboratory after rearing diapausing larvae collected in the field.

<table><thead><tr><th>Extract/library ID</th><th>Developmental stage</th><th>Tissue</th><th>n</th><th>Raw read count</th></tr></thead><tbody><tr><td>Lib1</td><td>egg</td><td>whole egg</td><td>3 egg masses (ca. 60 eggs)</td><td>58,941,438</td></tr><tr><td>Lib2</td><td>L5</td><td>whole body</td><td>9</td><td>54,020,124</td></tr><tr><td>Lib3</td><td>L5</td><td>hemolymph</td><td>31</td><td>43,406,048</td></tr><tr><td>Lib4</td><td>Female adult</td><td>Head/thorax</td><td>4</td><td>43,465,884</td></tr><tr><td>Lib5</td><td>Female adult</td><td>Abdomen</td><td>4</td><td>54,635,026</td></tr><tr><td>Lib6</td><td>Male adult</td><td>Head/thorax</td><td>4</td><td>38,226,332</td></tr><tr><td>Lib7</td><td>Male adult</td><td>Abdomen</td><td>4</td><td>32,314,096</td></tr><tr><td>Total</td><td></td><td></td><td></td><td>325,008,948</td></tr></tbody></table>

## 2.2. *De novo* transcriptome

In March 2011, diapausing larvae were collected in mugwort stems from Nadarzin (Poland) and then reared in the laboratory to obtain fresh tissues from the following developmental stages: eggs and larval whole body and hemolymph from the fifth instar. At the adult stage we sampled and separated heads/thorax from abdomens and males from females. In total, we prepared 7 RNA extracts
corresponding to these various stages and tissues (Table 5). RNA quality and concentration were evaluated using the RNA 6000 Nano kit with an Agilent 2100 Bioanalyser (Agilent Technologies, Palo Alto, CA, USA). Indexed cDNA libraries with an insert size of 150–200 bp were constructed for each developmental stage and tissue extract using the Illumina TruSeq RNA sample preparation kit. Subsequently, the libraries were sequenced by an Illumina HiSeq 2000 System at GATC Biotech (Konstanz, Germany) using the paired-end protocol. Between 32 and 59 million 100 bp raw reads were generated per library (Table 5). After a series of read cleaning and normalization steps, transcripts were reconstructed with Trinity, CD-HIT-EST and CAP3 [see Supplementary material, [18–20]]. The de novo transcriptome assembly had an overall size of 49.2 Mb and resulted in 44,564 transcripts, grouped into 33,080 unigenes (Table 4). Transcript lengths ranged from 201 to 27,559 bp. The N50 and N90 lengths were 2006 and 429 bp, respectively. Regarding the conserved eukaryotic gene sets, 84% of the CEGMA and 88% of the BUSCO genes were identified. Furthermore, 71% of the conserved arthropod BUSCO genes were present within the reference transcriptome. Coding sequences (CDS) were predicted with FrameDP [21] for 18,494 transcripts of which 12,515 were complete. An additional analysis of the 26,070 transcripts without predicted CDS with FEELnc [22] identified 10,883 potential long non-coding RNAs and sequences for which the CDS was either too fragmented (n = 15,132) or not present at all (n = 55). Further comparative analysis of the reference transcripts with the Lepbase [23] reference protein set detected probable homologs for 20,835 transcripts. OHR analyzes [12] on the best matches indicated that 61% of the CDS were reconstructed at least at 60% of the corresponding reference lepidopteran protein homolog, whereas 43% of the transcript CDS were assembled at full length.

**Acknowledgements**

This work was supported by a grant from INRA (AIP Bioressources 2011 (grant id SCOOP, INRA)). Most of the bioinformatics analyses were performed on the CBGP HPC computational platform (Montferrier-sur-Lez, France, http://www6.montpellier.inra.fr/cbgp_eng/Platforms/Computing-platform). We are grateful to the bioinformatics platforms ABiMS (Roscoff, France, http://abims.sb-roscoff.fr/), GenOuest (Rennes, France, https://www.genouest.org/) and Genotoul (Toulouse, France, http://bioinfo.genotoul.fr/) for providing access to their computational resources. We also thank Denis Bourguet for his comments on an earlier version of the manuscript.

**Transparency document. Supporting information**

Supplementary data associated with this article can be found in the online version at http://dx.doi.org/10.1016/j.dib.2018.01.073.

**Appendix A. Supporting information**

Supplementary data associated with this article can be found in the online version at http://dx.doi.org/10.1016/j.dib.2018.01.073.

**References**

[1] A.N. Frolov, et al., Reconsidering the taxonomy of several Ostrinia species in the light of reproductive isolation: a tale for Ernst Mayr, Biol. J. Linn. Soc. 91 (2007) 49–72.

[2] M.T. Bethenod, et al., Genetic isolation between two sympatric host plant races of the European corn borer, *Ostrinia nubilalis* Hubner. II: assortative mating and host-plant preferences for oviposition, Heredity 94 (2005) 264–270.

[3] V. Calcagno, et al., Sympatric host races of the European corn borer: adaptation to host plants and hybrid performance, J Evol. Biol. 20 (2007) 1720–1729.
[4] V. Calcagno, et al., Parallel evolution of behaviour during independent host-shifts following maize introduction into Asia and Europe, Evol. Appl. 10 (2017) 881–889.

[5] T. Malausa, et al., Differences in oviposition behaviour of two sympatric sibling species of the genus Ostrinia, Bull. Entomol. Res. 98 (2008) 193–201.

[6] M. Orsucci, et al., Host specialization involving attraction, avoidance and performance, in two phytophagous moth species, J. Evol. Biol. 29 (2016) 114–125.

[7] Y. Thomas, et al., Genetic isolation between two sympatric host-plant races of the European corn borer, *Ostrinia nubilalis* Hubner. *I. Sex pheromone, moth emergence timing, and parasitism, Evolution* 57 (2003) 261–273.

[8] T. Malausa, et al., Genetic structure and gene flow in French populations of two Ostrinia taxa: host races or sibling species? Mol. Ecol. 16 (2007) 4210–4222.

[9] A. Midamegbe, et al., Scanning the European corn borer (Ostrinia spp.) genome for adaptive divergence between host-affiliated sibling species, Mol. Ecol. 20 (2011) 1414–1430.

[10] H. Alexandre, et al., When history repeats itself: exploring the genetic architecture of host-plant adaptation in two closely related lepidopteran species, PLoS One 8 (2013) e69211.

[11] C.B. Wadsworth, et al., A recombination suppressor contributes to ecological speciation in OSTRINIA moths, Heredity 114 (2015) 593–600.

[12] B. Gschloessl, et al., De novo transcriptomic resources for two sibling species of moths: *Ostrinia nubilalis* and *O. scapulalis*, BMC Res. Notes 6 (2013).

[13] J.J. Doyle, J.L. Doyle, A rapid DNA isolation procedure for small quantities of fresh leaf tissue, Phytochem. Bull. 19 (1987) 11–15.

[14] S. Gnerre, et al., High-quality draft assemblies of mammalian genomes from massively parallel sequence data, Proc. Natl. Acad. Sci. USA 108 (2011) 1513–1518.

[15] G. Parra, et al., CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes, Bioinformatics 23 (2007) 1061–1067.

[16] F.A. Simao, et al., BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs, Bioinformatics 31 (2015) 3210–3212.

[17] B.L. Cantarel, et al., MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes, Genome Res. 18 (2008) 188–196.

[18] L. Fu, et al., CD-HIT: accelerated for clustering the next-generation sequencing data, Bioinformatics 28 (2012) 3150–3152.

[19] M.G. Grabherr, et al., Full-length transcriptome assembly from RNA-Seq data without a reference genome, Nat. Biotechnol. 29 (2011) 644–652.

[20] X. Huang, A. Madan, CAP3: a DNA sequence assembly program, Genome Res. 9 (1999) 868–877.

[21] J. Gouzy, et al., FrameDP: sensitive peptide detection on noisy matured sequences, Bioinformatics 25 (2009) 670–671.

[22] V. Wucher, et al., FEELnc: a tool for long non-coding RNA annotation and its application to the dog transcriptome, Nucleic Acids Res. 45 (2017) e57.

[23] R.J. Challis, et al., Lepbase: The Lepidopteran genome database, bioRxiv, 2016.