The total DNA of the genomes was recovered using the EZ1 biorobot (Qiagen, Courtaboeuf, France) and the EZ1 DNA tissue kit. Sequencing was performed using MiSeq technology (Illumina, San Diego, CA, USA) with the Nextera Mate Pair sample prep kit and Nextera XT paired end (Illumina), as previously described (Morel et al. 2015 (link)). Several bioinformatic tools, including Velvet (Zerbino and Birney 2008 (link)), Spades (Bankevich et al. 2012 (link)), and SOAPdenovo (Luo et al. 2012 (link)) were used to assemble the genomes. GapCloser software (Xu et al. 2019 (link)) was used to reduce assembly gaps. Scaffolds which had fewer than 800 base pairs (bp) or had a depth value lower than 25% of the mean depth were removed. The best assembly was selected using different criteria (number of scaffolds, N50, number of N). The degree of genomic similarity of each strain was evaluated by processing sequences using the Orthologous ANI Tool (OAT) software (Lee et al. 2016 (link)). Furthermore, the Genome-to-Genome Distance Calculator (GGDC) web server, which is available online (http://ggdc.dsmz.de), was used to calculate digital DNA–DNA hybridisation (dDDH) values between the genomes of closest species, as previously described (Meier-Kolthoff et al. 2013 (link)).
Free full text: Click here