# Illumina and PacBio DNA sequencing data, de novo assembly and annotation of the genome of Aurantiochytrium limacinum strain CCAP_4062/1

Christian Morabito, Riccardo Aiese Cigliano, Eric Maréchal, Fabrice Rébeillé, Alberto Amato

► To cite this version:

Christian Morabito, Riccardo Aiese Cigliano, Eric Maréchal, Fabrice Rébeillé, Alberto Amato. Illumina and PacBio DNA sequencing data, de novo assembly and annotation of the genome of Aurantiochytrium limacinum strain CCAP_4062/1. Data in Brief, 2020, 31, pp.105729. 10.1016/j.dib.2020.105729. hal-02969120

HAL Id: hal-02969120

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

Submitted on 6 Jan 2021

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

Illumina and PacBio DNA sequencing data, de novo assembly and annotation of the genome of Aurantiochytrium limacinum strain CCAP_4062/1

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

Christian Morabito<sup>a,1</sup>, Riccardo Aiese Cigliano<sup>b</sup>, Eric Maréchal<sup>a</sup>, Fabrice Rébeillé<sup>a,*</sup>, Alberto Amato<sup>a,*</sup>

ARTICLE INFO

**Article history:**
Received 22 April 2020
Revised 7 May 2020
Accepted 13 May 2020
Available online 21 May 2020

**Keywords:**
Genome
Third generation sequencing
Next generation sequencing
Structural annotation
Biotechnology
Thraustochytrid

ABSTRACT

The complete genome of the thraustochytrid *Aurantiochytrium limacinum* strain CCAP_4062/1 was sequenced using both Illumina Novaseq 6000 and third generation sequencing technology PacBio RSII in order to obtain trustworthy assembly and annotation. The reads from both platforms were combined at multiple levels in order to obtain a reliable assembly, then compared to the *A. limacinum* ATCC® MYA1381™ reference genome. The final assembly was annotated with the help of strain CCAP_4062/1 RNAseq data. *A. limacinum* strain CCAP_4062/1 is an industrial strain used for the production of very long chain polyunsaturated fatty acids, like the docosahexaenoic acid that is an essential fatty acid synthesised only at very low pace in humans and vertebrates. Thraustochytrids in general and *Aurantiochytrium* more specifically, are used for carotenoid and squalene production as well. Beside their biotechnological interest, thraustochytrids play a crucial role in both inshore and oceanic basins ecosystems. Genome sequences will foster biotechnological as well as ecological studies.

**Article history:**
Received 22 April 2020
Revised 7 May 2020
Accepted 13 May 2020
Available online 21 May 2020

<footnote><sup>a</sup> Laboratoire de Physiologie Cellulaire Végétale, Université Grenoble Alpes, CEA, CNRS, INRAE, IRIG-LPCV, 38054 Grenoble Cedex 9, France</footnote>

<footnote><sup>b</sup> Sequentia Biotech Carrer d'Àlaba, 61, 08005 Barcelona, Spain</footnote>

<footnote>* Corresponding authors.</footnote>

<footnote>E-mail addresses: fabrice.rebeille@cea.fr (F. Rébeillé), alberto.amato@cea.fr (A. Amato).</footnote>

<footnote n="1">¹ Present address: Université Paris-Saclay, INRAE, MGP, 78350 Jouy-en-Josas, France.</footnote>
<table><caption>**Specifications Table**</caption><tr><td>Subject</td><td>Applied Microbiology and Biotechnology</td></tr><tr><td>Specific subject area</td><td>Marine eukaryotic microbiology</td></tr><tr><td>Type of data</td><td>DNA Sequencing Data</td></tr><tr><td>How data were acquired</td><td>The data were acquired by Next-Generation Sequencing technology using Illumina Novaseq 6000 and third generation sequencing technology using PacBio RSII platforms</td></tr><tr><td>Data format</td><td>Raw reads were deposited in GenBank. Analysed: The assembly (fasta file), the gene descriptions (txt file), and the gene models (gtf file) were deposited in Mendeley</td></tr><tr><td>Parameters for data collection</td><td>DNA was extracted from six day-old cultures.</td></tr><tr><td>Description of data collection</td><td>Whole-genome sequencing, genome assembly, and annotation</td></tr><tr><td>Data source location</td><td>Institution: LPCV-IRIG<br>City: Grenoble<br>Country: France<br>Aurantiochytrium limacinum strain was collected in Mayotte (Indian Ocean, 12°48'51.8"S, 45°14'21.7"E)</td></tr><tr><td>Data accessibility</td><td>Repository name: NCBI BioProjects<br>Data identification number: PRJNA612804<br>Direct URL to data: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA612804<br>Repository name: Mendeley<br>Data identification number: 10.17632/v3w485jnz5.2<br>Direct URL to data: http://dx.doi.org/10.17632/v3w485jnz5.2</td></tr></table>

**Value of the Data**

* The biotechnology based on thraustochytrids has been gaining importance in the last decades. Physiology and life cycle traits can be investigated as well under a molecular point of view.

* The dataset presented here can provide information to both academic and private laboratories for reverse genetics studies.

* The genome can be advantageous for biotechnological as well as physiological studies aiming at improving growth or lipid production in thraustochytrids.

**1. Data Description**

Thraustochytrids are marine non-photosynthetic protists whose ability to produce high amounts of lipids, like long-chain polyunsaturated fatty acids [1] used in nutraceutical, and some terpene derivatives [2] like astaxanthin, a potent antioxidant agent, and squalene [3], has attracted industrial interest [4]. Although ecologically relevant [5], *Aurantiochytrium limacinum* biotechnological attractiveness is the main trigger for transcriptomic and genomic studies. To date an increasing number of thraustochytrid genomes [6,7] and transcriptomes [8,9] have been sequenced, and this has fostered reverse genetics studies [10, 11]. Here we present the genome sequenced by Illumina and PacBio of an *A. limacinum* strain that groups in the same 18S rDNA clade as the type species of the genus *Aurantiochytrium* [1], *A. limacinum* strain SR21 stored at the American Type Culture Collection under the entry ATCC® MYA1381TM. The taxonomic and systematic environment of thraustochytrids is very complex and requires profound rearrangements [12], but the phylogenetic position of strain CCAP_4062/1 was recently confirmed [1].
<table><caption>**Table 1**

Description of the genomics datasets used in this study.</caption><thead><tr><th></th><th>PacBio Sequel RSII</th><th>Illumina NovaSeq 6000</th></tr></thead><tbody><tr><td>Sequenced Bases</td><td>12,160,726,429 bp</td><td>9,045,844,656 bp</td></tr><tr><td>Number of Reads</td><td>770,207</td><td>59,906,256</td></tr><tr><td>Sequencing Layout</td><td>Single End Long Reads</td><td>Paired End 2 × 150 bp</td></tr><tr><td>Max Read Length</td><td>120,000 bp</td><td>150 bp</td></tr><tr><td>Read N50</td><td>37,932 bp</td><td>150 bp</td></tr><tr><td>Estimate Genome Coverage</td><td>202 ×</td><td>150 ×</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Schematic representation of the genome assembly pipeline. Raw PacBio reads were corrected using the Illumina data and using three iterations of the program LoRDEC. The corrected PacBio reads were used to create a draft assembly with the tool wtdbg2, thus producing the first raw assembly. The latter was then polished using the Illumina reads and performing five iterations of Pilon corrections and one run of REAPR to remove misassemblies. The polished assembly was used together with the Illumina reads to perform an assembly with Spades. The obtained assembly was polished with 10 iterations of Pilon, then gap closing was performed with LR_GapCloser.</figcaption></figure>

## 2. Genome assembly

The starting dataset included PacBio data obtained with the RSII long read technology and Illumina Novaseq 6000 reads (Table 1). The PacBio data included 770,000 polymerase reads with an N50 of 37,932 bp and a total of 12.16 Gbp. The Illumina data consisted of 59.9 million paired-end 2 × 150 bp reads corresponding to 9 Gbp [13]. Considering that the estimated genome size of *A. limacinum* is about 60 Mbp, the PacBio and Illumina data corresponded to a 202 × and 150 × coverage, respectively.

In order to obtain a highly accurate genome assembly, the Illumina and PacBio datasets were combined at multiple levels following the bioinformatics pipeline showed in Fig. 1. The Illumina dataset was quality-checked and trimmed in order to remove adapters and low-quality sequences thus obtaining 47 million high quality reads. PacBio reads were corrected using the high-quality Illumina reads and used to perform a draft assembly with the software wtdbg2. The obtained genome sequence was polished using the Illumina reads and then used as a guide
<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Dotplot obtained by aligning the Aurli1 reference genome [13] assembly from *A. limacinum* ATCC® MYA1381™ (X-axis) against the *Aurantiochyrium limacinum* strain CCAP_4062/1 scaffolds.</figcaption></figure>

for a new assembly with Spades. The newly obtained genome was subjected to error correction, scaffolding and gap-closing. The procedure led to the final assembly.

The obtained assembly was aligned against the *A. limacinum* ATCC® MYA1381™ reference genome [13] assembly as shown in Fig. 2. The analysis highlighted a high degree of collinearity between the two genomes with more than 93% of the sequences having an identity higher than 75%. Indeed, the Average Nucleotide Identity (ANI) between the two genomes was 98.89%.

Given the high collinearity and similarity of the two genomes, a scaffolding step was performed with the software DGenies using the *A. limacinum* ATCC® MYA1381™ assembly as a reference. The final assembly was then generated by creating a super scaffold with the smallest unplaced contig. The obtained final assembly included 478 sequences with a total assembled genome of 62 Mbp (Table 2).

The quality of the assembly was evaluated by mapping back the Illumina and PacBio reads to measure the percentage of alignment and its quality. About 97.7% of the Illumina reads were mapped back to the assembly with a mean mapping quality of 52 and 95% of the reads were properly paired. About 95% of the PacBio reads were mapped to the assembly with a mean mapping quality of 48. In addition, the tool BUSCO was used to predict the presence of single
<table><caption>**Table 2**
Genome assembly statistics.</caption><thead><tr><th></th><th>Aurantiochytrium limacinum strain CCAP_4062/1</th></tr></thead><tbody><tr><td>Number of Contigs</td><td>478</td></tr><tr><td>Genome Size</td><td>62,086,374 bp</td></tr><tr><td>Number of Contigs larger than 50 Kbp</td><td>210</td></tr><tr><td>N50</td><td>358,008 bp</td></tr><tr><td>L50</td><td>51</td></tr><tr><td>Largest Contig</td><td>2029,424 bp</td></tr><tr><td>GC Content</td><td>45.66%</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Fig. 3.** Results of the BUSCO analysis highlighting the presence of complete and single copy eukaryotic genes in the assembly. Letters indicate the BUSCO categories presented in the figure, numbers indicate the number of genes composing a category. 'n' indicate the total number of genes in all BUSCO categories.</figcaption></figure>

copy conserved Eukaryotic genes in the assembly (Fig. 3), highlighting that 86% of the genes were complete.

### 3. RNA-seq guided annotation of the genome

Once the final genome assembly was produced, genome annotation was performed. For this purpose, RNA-seq dataset [8] was used to assist in the gene prediction. The bioinformatics pipeline consisted in mapping the RNA-seq reads on the genome assembly followed by a com-
<figure><img src="image_7.png" /><figcaption>**Fig. 4.** Histograms showing the distribution of alignment and identity percentage between *Aurantiochytrium limacinum* ATCC® MYA1381TM (reference) and *Aurantiochytrium limacinum* strain CCAP_4062/1 predicted proteins (present study).</figcaption></figure>

bination of GeneMark and Augustus in order to predict the gene models. The Braker2 pipeline was used for this purpose. As a first step, the aligned RNA-seq reads were used to detect introns and transcribed loci. More than 95% of the reads could be correctly mapped on the obtained assembly. Then, the information used in the previous step was used to feed GeneMark to create an HMM model of splicing sites which, in turn, was used to train Augustus. In addition, the proteins from the *A. limacinum* ATCC® MYA1381TM genome were used as input to identify coding genes by sequence similarity. The generated Augustus model was used to create the final annotation. About 14,470 protein coding genes were predicted, of which 12,610 (87%) had a match with an *A. limacinum* ATCC® MYA1381TM protein. Gene descriptions and Gene Ontology annotations were associated to 11,208 proteins using the Pannzer 2 pipeline.

To assess the quality of the annotation, the RNA-seq reads were used to detect the expression of the annotated genes and about 75% of the reads could be properly associated to the annotated genes. In addition, the length of the predicted proteins was compared with the length of *A. limacinum* ATCC® MYA1381TM proteins. About 90% of the predicted proteins covered at least 90% of the *A. limacinum* ATCC® MYA1381TM proteins, showing that the genes were complete (Fig. 4). Finally, a BUSCO analysis was performed on the predicted proteins using the Eukaryotic
gene-set which highlighted the presence of 90% of complete genes, of which 87% were in single copy.

# 4. Experimental design, materials, and methods

## 4.1. The strain

By means of pollen grain bait method [14] carried out on samples from coastal seawater gathered in Mayotte (Indian Ocean, 12°48′51.8″S, 45°14′21.7″E), a cell was isolated and an axenic culture established. Once the culture was proved to be monospecific and axenic, it was deposited at the Culture Collection of Algae and Protozoa (CCAP) under the accession number CCAP_4062/1.

## 4.2. Culture conditions and DNA extraction

In order to accumulate biomass for DNA extraction, a six-day old culture was inoculated in 50 mL fresh R medium [1,5] at $5 \times 10^5$ cells mL$^{-1}$ in sterile 250 mL-Erlenmeyer glass flasks and incubated for six days at 20°C and 100 rpm. Genomic DNA was extracted by phenol:chloroform:isoamyl alcohol (25:24:1) on lysed cells and precipitated with Na-Acetate 3 M pH 5 and absolute ethanol. Quality and concentration were estimated using a Thermo Scientific™ NanoDrop 2000 Spectrophotometer and Qubit Flex Fluorometer. Genomic DNA was sent to Macrogen (Korea) for both Illumina and PacBio sequencing. The sample was prepared according to a guide for preparing SMRTbell template for sequencing on the PacBio Sequel System. The templates were sequenced using SMRT Sequencing. Illumina TruSeq Nano DNA Kit was used to generate the Illumina library according to manufacturer's specifications. Illumina sequencing was performed on a Novaseq-6000 producing paired-end 2 × 150 bp reads.

## 4.3. Bioinformatics methods

Raw Illumina reads were analysed with FASTQC [15] to obtain quality statistics, then BB-Duk v38.75 [16] was used to perform trimming and clipping (minimum base quality 35 and minimum read length 35 bp). PacBio reads were corrected using three iterations of the software LoRDEC v0.3 [17] together with the trimmed Illumina reads, the three iterations were performed with three different K-mer lengths: 19, 31 and 41 bp. The software wtdbg2 v2.5 [18] was used to obtain the first draft assembly with the PacBio corrected reads, an estimated genome size of 60 Mbp was indicated. The Illumina reads were then mapped against the obtained assembly with minimap2 version 2.17-r954-dirty [19] and the results treated with Pilon v1.23 [20] and REAPR v1.0.18 [21] in order to fix mismatches and assembly rearrangements. This process was iterated five times obtaining a polished first assembly. The obtained assembly was used as input, together with the Illumina reads, for Spades v3.14.0 [22] to perform a second genome assembly, which was then polished by 10 iterations of Pilon corrections. Gap closing was performed using the LR_GapCloser algorithm [23] using the polished first assembly as input. The obtained scaffolds were aligned with DGenies [24] against the *Aurantiochytrium limacinum* ATCC® MYA1381™ genome [13] using minimap2 as an aligner. The FASTA file containing the scaffolds ordered according to the *A. limacinum* ATCC® MYA1381™ genome alignment was finally downloaded. The CCAP_4062/1 scaffolds with no match against the *A. limacinum* ATCC® MYA1381™ genome were concatenated introducing 40 Ns between each scaffold to generate a FASTA of unplaced sequences. Genome assembly statistics were produced with the software
QUAST [25] whereas BUSCO analyses were performed with the software version 3 [26] and the eukaryote_odb9 dataset. Average Nucleotide Identity between A. limacinum ATCC® MYA1381™ and A. limacinum strain CCAP_4062/1 genomes was calculated with FastANI [27].

For genome structural annotations, RNA-seq reads [8] were trimmed and clipped using BB-Duk using the same parameters mentioned above. Reads were mapped against the genome assembly using STAR v2.7.3a [28] in double pass mode. The obtained mapping was used as input for Braker2 v2.1.0 [29] also providing a FASTA file containing the protein sequences of A. limacinum ATCC® MYA1381™. The genome assembly was repeat masked with RepeatMasker version open-4.0.9 [30] before performing the annotation selecting the option -species stramenopiles. Gene expression levels were obtained using Kallisto v0.46.0 [31] against the predicted transcript sequences. Finally, gene functional annotations were obtained using the software PANNZER2 [32] with the following options: Minimum query coverage 0.4 or minimum sbjct coverage 0.4, and Minimum alignment length 50. The ARGOT [33] scoring function implemented in PANNZER2 as default advanced option for Gene Ontology Annotation that proved to be the best [32], was chosen. The functional annotation table (available at doi:http://dx.doi.org/10.17632/v3w485jnz5.2) contains the description, Gene Ontology and KEGG Enzymes identified through PANNZER2. The four columns contain 1. **Locus**: The transcript ID created during the gene prediction step; 2. **Annotation Type**, which can include: DE (general description), MFARGOT (Gene Ontology Molecular Function), BPARGOT (Gene Ontology Biological Process), CCARGOT (Gene Ontology Cellular Component), ECARGOT (EC number); 3. **Annotation ID**, which contains the accession number in the case of MFARGOT, CCARGOT, BPARGOT and ECARGOT and the score of the prediction in the case of DE; 4. **Description** of the annotation ID, which includes the gene description in the case of DE, the Gene Ontology description for MFARGOT, CCARGOT, BPARGOT and the associated GO ID for the ECARGOT. The sequences of CCAP_4062/1 predicted proteins were aligned with the A. limacinum ATCC® MYA1381™ proteins using the BLASTp algorithm setting a minimum evalue of 0.001.

**Declaration of Competing Interest**

The authors declare that they have no known competing financial interests or personal relationships which have, or could be perceived to have, influenced the work reported in this article.

**Acknowledgments**

CM, EM, FR, and AA were supported by the French National Research Agency (ANR-10-LABEX-04, GRAL Labex; ANR-11-BTBR-0008, Océanomics; ANR-17-EURE-0003, EUR CBS) and by the Trans'Alg Bpifrance PSPC partnership.

**Author contribution**

AA, EM, FR conceptualization; AA, CM, RAC formal analysis and investigation; AA, RAC writing and visualization; EM, FR funding acquisition; AA, EM, FR project administration.

**Supplementary materials**

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.dib.2020.105729.
## References

[1] Y. Dellero, O. Cagnac, S. Rose, K. Seddiki, M. Cussac, C. Morabito, J. Lupette, R. Aiese Cigliano, W. Sanseverino, M. Kuntz, J. Jouhet, E. Maréchal, F. Rébeillé, A. Amato, Proposal of a new thraustochytrid genus *Hondaea* gen. nov. and comparison of its lipid dynamics with the closely related pseudo-cryptic genus *Aurantiochytrium*. Algal Res. 35 (2018) 125-141. doi:10.1016/j.algal.2018.08.018.

[2] Y. Xie, B. Sen, G. Wang, Mining terpenoids production and biosynthetic pathway in thraustochytrids. Bioresour. Technol. 244 (2017) 1269-1280. doi:10.1016/j.biortech.2017.05.002.

[3] Patel, A., Liefeldt, S., Rova, U., Christakopoulos, P., Matsakas, L., 2020. Co-production of DHA and squalene by thraustochytrid from forest biomass. Sci. Rep. 10, e1992. doi:10.1038/s41598-020-58728-7.

[4] Morabito, C., Bournaud, C., Maës, C., Schuler, M., Aiese Cigliano, R., Dellero, Y., Maréchal, E., Amato, A., Rébeillé, F., 2019. The lipid metabolism in thraustochytrids. Prog. Lipid Res. 76, e101007. doi:10.1010/j.plipres.2019.101007.

[5] Y. Dellero, S. Rose, C. Metton, C. Morabito, J. Lupette, J. Jouhet, E. Maréchal, F. Rébeillé, A. Amato, Ecophysiology and lipid dynamics of a eukaryotic mangrove decomposer, Environ. Microbiol. 20 (2018) 3057-3068. doi:10.1111/1462-2920.14346.

[6] Seddiki, K., Godart, F., Aiese Cigliano, R., Sanseverino, W., Barakat, M., Ortet, P., Rébeillé, F., Maréchal, E., Cagnac, O., Amato, A., 2018. Sequencing, de novo assembly, and annotation of the complete genome of a new thraustochytrid species, strain CCAP_4062/3. Genome Announc. 6, e01335-17. doi:10.1128/genomeA.01335-17.

[7] B. Liu, H. Ertesvåg, I.M. Aasen, O. Vadstein, T. Brautaset, T.M. Heggeset, Draft genome sequence of the docosahexaenoic acid producing thraustochytrid *Aurantiochytrium* sp. T66. Genom. Data 8 (2016) 115-116. doi:10.1016/j.gdata.2016.04.013.

[8] Dellero, Y., Maës, C., Morabito, C., Schuler, M., Bournaud, C., Aiese Cigliano, R., Maréchal, E., Amato, A., Rébeillé, F., 2020. The zoospores of the thraustochytrid *Aurantiochytrium limacinum*: transcriptional reprogramming and lipid metabolism associated to their specific functions. Environ. Microbiol. doi:10.1111/1462-2920.14978.

[9] Heggeset, T.M.B., Ertesvåg, H., Liu, B., Ellingsen, T.E., Vadstein, O., Aasen, I.M., 2019. Lipid and DHA-production in *Aurantiochytrium* sp. - responses to nitrogen starvation and oxygen limitation revealed by analyses of production kinetics and global transcriptomes. Sci. Rep. 9, e19470. doi:10.1038/s41598-019-55902-4.

[10] T. Watanabe, R. Sakiyama, Y. Limi, S. Sekine, E. Abe, K.H. Nomura, K. Nomura, Y. Ishibashi, N. Okino, M. Hayashi, M. Ito, Regulation of TG accumulation and lipid droplet morphology by the novel TLDP1 in *Aurantiochytrium limacinum* F26-b. J. Lipid Res. 58 (2017) 2334-2347. https://www.jlr.org/content/58/12/2334.

[11] Nutahara, E., Abe, E., Uno, S., Ishibashi, Y., Watanabe, T., Hayashi, M., Okino, N., Ito, M., 2019. The glycerol-3-phosphate acyltransferase PLAT2 functions in the generation of DHA-rich glycerolipids in *Aurantiochytrium limacinum* F26-b. PLoS ONE 14, e0211164. doi:10.1371/journal.pone.0211164.

[12] R.A. Andersen, E. Ganuza, Nomenclatural errors in the Thraustochytridiales (Heterokonta/Staminipilina), especially with regard to the type species of *Schizochytrium*. Not. Algar. 64 (2018) 1-8.

[13] Info • Aurantiochytrium limacinum ATCC MYA-1381 https://phycocosm.ijgi.doe.gov/Aurli1/Aurli1.info.html.

[14] A. Gupta, S. Wilkens, J.L. Adcock, M. Puri, C.J. Barrow, Pollen baiting facilitates the isolation of marine thraustochytrids with potential in omega-3 and biodiesel production. J. Ind. Microbiol. Biotechnol. 40 (2013) 1231-1240. doi:10.1007/s10295-013-1324-0.

[15] Babraham Bioinformatics https://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

[16] Bushnell, B., Rood, J., Singer, E., 2017. BBMerge – accurate paired shotgun read merging via overlap. PLoS ONE 12, e0185056. doi:10.1371/journal.pone.0185056.

[17] L. Salmela, E. Rivals, LoRDEC: accurate and efficient long read error correction. Bioinformatics 30 (2014) 3506-3514. doi:10.1093/bioinformatics/btu538.

[18] J. Ruan, H. Li, Fast and accurate long-read assembly with wtdbg2. Nat. Methods 17 (2020) 155-158. doi:10.1038/s41592-019-0669-3.

[19] H. Li, Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34 (2018) 3094-3100. doi:10.1093/bioinformatics/bty191.

[20] Walker, B.J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., Cuomo, C.A., Zeng, Q., Wortman, J., Young, S.K., Earl, A.M., 2014. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9, e112963. doi:10.1371/journal.pone.0112963.

[21] Hunt, M., Kikuchi, T., Sanders, M., Newbold, C., Berriman, M., Otto, T.D., 2013. REAPR: a universal tool for genome assembly evaluation. Genome Biol. 14, R47. doi:10.1186/gb-2013-14-5-r47.

[22] CAB Center for Algorithmic Biotechnology St. Petersburg State University http://cab.spbu.ru/software/spades/.

[23] Xu, G.-C., Xu, T.-J., Zhu, R., Zhang, Y., S.-Q., Wang, H.-W., Li, J.-T., 2019. LR_Gapcloser: a tiling path-based gap closer that uses long reads to complete genome assembly. Gigascience 8, giy157. doi:10.1093/gigascience/giy157.

[24] Cabanettes, F., Klopp, C., 2018. D-GENIES: dot plot large genomes in an interactive, efficient and simple way. PeerJ. 6, e4958. doi:10.7717/peerj.4958.

[25] A. Gurevich, V. Saveliev, N. Vyahgi, G. Tesler, QUAST: quality assessment tool for genome assemblies, Bioinformatics 29 (2013) 1072-1075. doi:10.1093/bioinformatics/btt086.

[26] M. Seppey, M. Manni, E.M. Zdobnov, BUSCO: assessing genome assembly and annotation completeness, in: M. Kollmar (Ed.) Gene Prediction. Methods in Molecular Biology, Vol 1962. Humana, New York, 2019, pp. 227-245. doi:10.1007/978-1-4939-9173-0_14.

[27] Jain, C., Rodriguez-R, L.M., Phillippy, A.M., Konstantinidis, K.T., Aluru, S. 2018. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9, e5114. doi:10.1038/s41467-018-07641-9.

[28] STAR https://github.com/alexdobin/STAR.

[29] Braker2 https://github.com/Gaius-Augustus/BRAKER.

[30] RepeatMasker http://www.repeatmasker.org/.

[31] N.L. Bray, H. Pimentel, P. Melsted, L. Pachter, Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 34 (2016) 525-527. doi:10.1038/nbt.3519.
[32] P. Törönen, A. Medlar, L. Holm, PANNZER2: a rapid functional annotation web server, Nucl. Acids Res. 46 (2018) W84–W88. doi:10.1093/nar/gky350.

[33] Falda, M., Toppo, S., Pescarolo, A., Lavezzo, E., Di Camillo, B., Facchinetti, A., Cilia, E., Velasco, R., Fontana, P. 2012. Argot2: a large scale function prediction tool relying on semantic similarity of weighted Gene Ontology terms. BMC Bioinform. 13, S14. doi:10.1186/1471-2105-13-S4-S14.