Isolation and Description of Catonella massiliensis sp. nov., a Novel Catonella Species, Isolated from a Stable Periodontitis Subject

The genus Catonella currently counts a unique species, C. morbi, isolated from periodontal pockets and associated with periodontitis and endodontic infections. This study contributed to the taxonomical and clinical knowledge of this genus by describing a novel species isolated from a saliva sample from a man in clinical gingival health following successful treatment of periodontitis. Morphological and chemotaxonomic characteristics were investigated using different growth conditions, pH, and temperature. Cellular fatty acid methyl ester (FAME) analysis was conducted by gas chromatography/mass spectrometry (GC/MS). Phylogenetic analysis based on 16S rRNA, orthologous average nucleotide identity (OrthoANI), and digital DNA-DNA hybridization (dDDH) relatedness were performed. Strain Marseille-Q4567T was found to be an anaerobic and non-spore-forming rod-shaped bacterium that grew at 28–41.5 °C (optimum 37 °C), pH 6.5–8.5 (optimum pH 7.5), and 5–10 g/L of NaCl (optimum 5 g/L). The predominant cellular fatty acid was C16:0 (64.2%), followed by unsaturated structures C18:1n9 (12.5%) and C18:2n6 (7.8%). Based on 16S rRNA sequence comparison, the closest phylogenetic neighbor was C. morbi ATCC 51271T (98.23% similarity). The OrthoANI and dDDH values between strain Q4567T and C. morbi ATCC 51271T were respectively 79.43% and 23.8%. Therefore, we concluded that strain Marseille-Q4567T represents a novel species of the genus Catonella, for which the name Catonella massiliensis sp. nov. is proposed (= CSUR Q4567).


Introduction
The human oral microbiome consists of a wide range of microorganisms, including bacteria, fungi, protozoa, viruses, and archaea and is thus, after the gut, the second largest human microbial community [1]. The symbiotic relationship between the oral microbiome and host ensures oral health [2]. The rupture of this balance and the establishment of a dysbiotic microbial community leads to oral disease, the most common being tooth decay and periodontal disease [3][4][5][6]. With the advent of next-generation sequencing and recent developments in bioinformatics, research on the diversity and the role of human oral microorganisms has increased dramatically [7]. The expanded Human Oral Microbiome Database (eHOMD) currently includes a total of 775 microbial species, of which 57% are officially named, 13% unnamed but cultivated, and 30% are known only as uncultivated phylotypes [8]. Several strategies have been used to cultivate "unculturable" bacteria, including media supplementation, diffusion chambers, and co-culture, as the isolation and culture of bacterial species are the only way to characterize them phenotypically and genotypically in an exhaustive way, and also to study their potential virulence [9,10]. The oral microbiome bacteria not yet cultured are all as likely to play a role in oral homeostasis or dysbiosis as cultured bacteria. It thus appears clear that the development of optimal culture means has a first order role in the investigation of oral microbial diversity and its potential implications in oral health or the development of dental and periodontal disease. Periodontitis is a chronic multifactorial inflammatory disease associated with dysbiotic plaque biofilms, which lead to the progressive destruction of the supporting structures of the teeth [11]. Periodontitis can be successfully treated, and stable periodontitis subjects have a bleeding on probing score in <10% of sites, shallow probing depths of 4 mm or less, and no 4 mm sites with bleeding on probing, optimal improvement in other clinical parameters, and lack of progressive periodontal destruction [12].
The genus Catonella (N.L. fem. dim. n. Catonella, in honor of Elizabeth P. Cato, an American microbiologist) belongs to the large family Lachnospiraceae and was first introduced into the literature by Moore and Moore (1994) [13]. C. morbi, isolated from periodontal pockets, is the only species of the genus Catonella currently described and has only been found in humans. This species was described as a putative periodontal pathogen, as it was more abundantly found in chronic periodontitis subjects [14]. Moreover, levels of C. morbi in subjects with chronic periodontitis decreased significantly after treatment in responsive sites, but not in refractory sites [15]. Recently, C. morbi has been thought to be involved in endodontic infections and was detected in 33% of the root canals associated with chronic apical periodontitis [16].
In the present study, we isolated strain Marseille-Q4567 T from a whole unstimulated saliva sample from a 51-year-old male in clinical gingival health following successful treatment of periodontitis. We used the cost-effective, accurate, and rapid routine identification by matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) for the identification of our strain. We compared strain Marseille-Q4567 T to its closely related phylogenetic neighbors and proposed that strain Marseille-Q4567 T represents a novel species of the genus Catonella, for which the name Catonella massiliensis sp. nov. is proposed (=CSUR Q4567).

Strain Isolation and Phenotypic Characteristics
Strain Marseille-Q4567 T was isolated from a whole unstimulated saliva sample provided by a 51-year-old male in clinical gingival health following successful treatment of periodontitis, and living in Marseille, France. Strain Marseille-Q4567 T could not be identified by MALDI-TOF MS, as the score was lower than 1.8 ( Figure 1). Growth was observed on Columbia agar with 5% sheep blood (BioMérieux, Marcyl'Etoile, France) at 37 °C after 24 h of incubation in an anaerobic atmosphere. Growth was also achieved in microaerophilic but not in aerobic conditions. The temperature range of the strain was determined to be 28-41.5 °C, with an optimum growth temperature 37 °C. The bacterial cells tolerated a pH of 6.5 to 8.5 (optimum pH 7.5) and an NaCl content Growth was observed on Columbia agar with 5% sheep blood (BioMérieux, Marcyl'Etoile, France) at 37 • C after 24 h of incubation in an anaerobic atmosphere. Growth was also achieved in microaerophilic but not in aerobic conditions. The temperature range of the strain was determined to be 28-41.5 • C, with an optimum growth temperature 37 • C. The bacterial cells tolerated a pH of 6.5 to 8.5 (optimum pH 7.5) and an NaCl content ranging from 5 to 10 g/L (optimum 5 g/L). After 10 min of thermal shock at 80 • C, this bacterium did not grow at 37 • C on Columbia agar with 5% sheep blood, which indicated the absence of spore formation. Surface colonies on Columbia agar with 5% sheep blood (BioMérieux, Marcy-l'Etoile, France) incubated at 37 • C for 2 days were approximately between 0.1 and 0.5 mm in diameter, smooth, circular, and translucent. Bacterial cells were Gram-negative curved rods, round-tipped, and approximatively 0.3 to 0.5 µm wide and 3.9 to 4.1 µm long, as determined by field emission scanning electron microscopy ( Figure 2). Growth was observed on Columbia agar with 5% sheep blood (BioMérieux, Marcyl'Etoile, France) at 37 °C after 24 h of incubation in an anaerobic atmosphere. Growth was also achieved in microaerophilic but not in aerobic conditions. The temperature range of the strain was determined to be 28-41.5 °C, with an optimum growth temperature 37 °C. The bacterial cells tolerated a pH of 6.5 to 8.5 (optimum pH 7.5) and an NaCl content ranging from 5 to 10 g/L (optimum 5 g/L). After 10 min of thermal shock at 80 °C, this bacterium did not grow at 37 °C on Columbia agar with 5% sheep blood, which indicated the absence of spore formation. Surface colonies on Columbia agar with 5% sheep blood (BioMérieux, Marcy-l'Etoile, France) incubated at 37 °C for 2 days were approximately between 0.1 and 0.5 mm in diameter, smooth, circular, and translucent. Bacterial cells were Gram-negative curved rods, round-tipped, and approximatively 0.3 to 0.5 μm wide and 3.9 to 4.1 μm long, as determined by field emission scanning electron microscopy ( Figure 2).  Using an API 50 CH strip, positive reactions were obtained for glycerol, erythritol, Darabinose, L-arabinose, D-ribose, D-xylose, L-xylose, D-adonitol, methyl βD-xylopyranoside, D-galactose, D-glucose, D-fructose, D-mannose, L-sorbose, L-rhamnose, dulcitol, inositol, D-mannitol, D-sorbitol, methyl αD-mannopyranoside, Methyl αD-glucopyranoside, N-acetylglucosamine, amygdalin, arbutin, esculin ferric citrate, salicin, D-cellobiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-trehalose, Inulin, D-melezitose, D-raffinose, amidon, glycogen, xylitol, gentiobiose, D-turanose, D-lyxose, D-tagatose, D-fucose, L-fucose, D-arabitol, L-arabitol, potassium gluconate, potassium 2-ketogluconate, and potassium 5ketogluconate. Using an API ZYM strip, a positive result was shown for alkaline phosphatase, leucine arylamidase, acid phosphatase, naphthol-AS-BI-phosphohydrolase, α-galactosidase, β-galactosidase, α-glucosidase, β-glucosidase, and N-acetyl-β-glucosaminidase. Strain Marseille Q4567 T was negative for oxidase and catalase activities. Comparison of phenotypic characteristics between strain Marseille-Q4567 T and its closely related species-Cuneatibacter caecimuris, Faecalicatena orotica, Faecalicatena contorta, Faecalicatena fissicatena, Anaerocolumna cellulosilytica, Catonella morbi, and Herbinix luporum-is listed in Table 1.  [20]. L-Arabinose The most abundant fatty acid was C 16:0 (64.2%), followed by unsaturated structures C 18:1n9 (12.5%) and C 18:2n6 (7.8%). No branched fatty acids were described. Minor amounts of other unsaturated and saturated structures were also detected ( Table 2).

Genome Sequencing Information and Genome Properties
The genome size of strain Marseille-Q4567 T was 3,122,925 bp long with a 38.8% G+C content. It was assembled into 3 contigs with a mean coverage of 131.6. It was deposited into GenBank under the accession number JAEPRJ000000000. Of the 2805 predicted genes, 2747 were protein-coding genes and 58 were RNAs (four 5S rRNA, four 16S rRNA, four 23S rRNA, 42 tRNAs, and 4 ncRNA). Genes with putative function (by COGs) were 1628 for strain Marseille-Q4567 T (Table 3). Finally, 1396 genes (49%) were annotated as hypothetical proteins for strain Marseille-Q4567 T . A circular map showing a complete view of the genome of strain Marseille-Q4567 T is shown in Figure 3. The mass screening for antimicrobial and virulence genes revealed the presence of lsa(C)_1 gene (92.02% identity, 100.00% coverage) which confers resistance to lincomycin, clindamycin, dalfopristin, pristinamycin_IIA, virginiamycin_M, and tiamulin.    Figure 3. A circular map generated using the CGView Server BETA [21] showing a complete view of the genome of strain Marseille-Q45657 T .

Figure 3.
A circular map generated using the CGView Server BETA [21] showing a complete view of the genome of strain Marseille-Q45657 T .

Comparison with Closely Related Bacterial Strains
16S rDNA-based similarity analysis of strain Marseille-Q4567 T against GenBank yielded the highest nucleotide sequence similarities of 98.23% sequence identity with C. morbi ATCC 51271 T (GenBank accession no. NR_026248.1). The 16S rRNA gene sequence was deposited into GenBank under the accession number MW410925. The phylogenetic tree highlighting the position of the strain relative to other closely related species is shown in Figure 4.

Comparison with Closely Related Bacterial Strains
16S rDNA-based similarity analysis of strain Marseille-Q4567 T against GenBank yielded the highest nucleotide sequence similarities of 98.23% sequence identity with C. morbi ATCC 51271 T (GenBank accession no. NR_026248.1). The 16S rRNA gene sequence was deposited into GenBank under the accession number MW410925. The phylogenetic tree highlighting the position of the strain relative to other closely related species is shown in Figure 4. Only bootstrap values of 70% or greater are shown. Type strains are indicated with superscript T. GenBank accession numbers of 16S rRNA indicated in parentheses. Sequences were aligned using MUSCLE with default parameters, phylogenetic inference was obtained using the Maximum likelihood method and MEGA X software [22]. Bootstrap values obtained by repeating the analysis 1000 times to generate a majority consensus tree are indicated at the nodes. There was a total of 1392 positions in the final dataset.
The genome of strain Marseille-Q4567 T was compared to the available genomes of nine closely related bacterial strains. The genome size of our strain (3,122,925 bp) was larger than that of H. luporum and E. oxidoreducens but smaller than that of C. caecimuris, C. morbi, R. hominis, F. fissicatena, A. cellulosilytica, F. contorta, and F. orotica (Table 4). The G+C content of strain Marseille Q4567 T (38.8%) was smaller than that of all compared species except H. luporum (35.3%), A. cellulosilytica (36.7%), and C. morbi (37.1%).  Only bootstrap values of 70% or greater are shown. Type strains are indicated with superscript T. GenBank accession numbers of 16S rRNA indicated in parentheses. Sequences were aligned using MUSCLE with default parameters, phylogenetic inference was obtained using the Maximum likelihood method and MEGA X software [22]. Bootstrap values obtained by repeating the analysis 1000 times to generate a majority consensus tree are indicated at the nodes. There was a total of 1392 positions in the final dataset. The genome of strain Marseille-Q4567 T was compared to the available genomes of nine closely related bacterial strains. The genome size of our strain (3,122,925 bp) was larger than that of H. luporum and E. oxidoreducens but smaller than that of C. caecimuris, C. morbi, R. hominis, F. fissicatena, A. cellulosilytica, F. contorta, and F. orotica (Table 4). The G+C content of strain Marseille Q4567 T (38.8%) was smaller than that of all compared species except H. luporum (35.3%), A. cellulosilytica (36.7%), and C. morbi (37.1%).  (Table 5). These values were lower than the 70% threshold used for delineating prokaryotic species, thus confirming that strain Marseille-Q4567 T represents a new species [10]. In addition, using OrthoANI analysis, strain Marseille-Q4567 T exhibited values ranging from 64.97% with C. caecimuris to 79.43% with C. morbi ( Figure 5).

Discussion
The human oral microbiome is a complex ecosystem whose balance with the host immune system is necessary in maintaining oral health. In this study, we aimed at isolating and describing a novel bacterial species using culturomics and taxonogenomic strategies in order to contribute new knowledge on the diversity of the oral microbiome.
The phenotypic and phylogenetic analysis of our strain Marseille-Q4567 T revealed several different characteristics when compared to other members of the family Lachnospiraceae, suggesting that it could be classified as a new species. The family Lachnospiraceae is a large family currently composed of 57 genera with a validly published and correct name. Among the members of the family Lachnospiraceae, strain Marseille-Q4567 T shared the highest 16S rRNA gene sequence similarities (98.23%) with C. morbi ATCC 51271 T , an anaerobic Gram-negative bacilli isolated from diseased periodontal pockets [13]. C. morbi has been thought to be involved in periodontal and endodontic infections [14,16]. Here, our strain Marseille-Q4567 T was isolated from the saliva of a 51-year-old male in clinical gingival health following successful treatment of periodontitis. The discovery of this new species alone does not allow us to associate its presence with the patient's periodontal condition. In this aim, additional studies exploring the prevalence of our strain Catonella massiliensis Marseille-Q4567 T according to the presence or absence of periodontal diseases will be of great interest. In addition, as stable periodontitis patients remain at higher risk for recurrent disease, it would be interesting to monitor the presence of Catonella species in treated and maintained periodontitis subjects and to correlate it with potential recurrence of the disease. Further studies will also be needed to investigate the potential involvement of the genus Catonella in gingival inflammation and periodontal destruction.
The genomic content (dDDH, orthoANI, and pangenome) with biochemical characteristics clearly indicated that strain Marseille-Q4567 T could be differentiated from the closely related species C. morbi, C. caecimuris, R. hominis, F. fissicatena, A. cellulosilytica, F. contorta, H. luporum, F. orotica, and E. oxidoreducens. Based on the results from phenotypic, chemotaxonomic, genomic, and phylogenetic analyses and data, we concluded that strain Marseille-Q4567 T represents a novel species of the genus Catonella, for which the name Catonella massiliensis sp. nov. is proposed (=CSUR Q4567).

Strain Isolation and Phenotypic Tests
A sample of unstimulated whole saliva was collected from a 51-year-old male in clinical gingival health following successful treatment of periodontitis and currently referred for periodontal maintenance at the periodontal department of the Pavillon Odontologique de la Timone, Marseille, France. The patient provided signed informed consent, and the study was approved by the Comité de Protection des Personnes (C.P.P.) Sud-Ouest et Outre-Mer 1 (no. ID RCB: 2020-A01234-35-CPP 1-20-075 ID 9806). Briefly, a 2-mL volume of whole unstimulated saliva was collected into a 50-mL centrifuge tube and transported to an anaerobic workstation within 1 h. After vortexing for 30 s, a 10-fold dilution series of the sample was prepared in phosphate-buffered saline 1x. Columbia agar sheep blood plates (BioMérieux, Marcy l'Etoile, France) were inoculated with 50 µL each of a 10 −4 to 10 −8 diluted plaque suspension. After 48 h of incubation in an anaerobic atmosphere (AnaeroGen Compact; Oxoid, Thermo Scientific, Dardilly, France) at 37 • C, the culture plates were inspected using a magnifying glass, and any microcolonies or colonies showing satellitism were passaged onto a fresh Columbia agar sheep blood plate. Matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) protein analysis was performed with a Microflex LT mass spectrometer (Bruker Daltonics, Bremen, Germany; external mass spectrometer calibration accuracy ± 300 ppm), as previously reported [23]. Briefly, each isolate colony was deposited on a 96 polished steel MALDI target and then coated with 1 µL of matrix solution containing α-cyano-4-hydroxycinnamic acid diluted in 500 µL of acetonitrile, 250 µL of 10% trifluoroacetic acid, and 250 µL of HPLC-grade water. The matrix sample was then crystallized by air-drying at room temperature as previously described [24]. Two spots were systematically created for each colony, and each isolate was characterized by at least 12 spots. The obtained spectra were imported into the BioTyper-RTC TM version 3.0 software (Bruker Daltonics GmbH) and analyzed by standard pattern matching (with default parameter settings). Interpretation of the scores was carried out as previously reported [23]. One purified strain, designated Marseille-Q4567 T and deposited in the Collection de Souches de l'Unité des Rickettsies under accession number Q4567, could not be identified by MALDI-TOF MS.
Gram staining was carried out using standard Gram stain, and morphological characteristics were observed with a field emission scanning electron microscope (SU5000 FE-SEM, Hitachi High-Tech, HHT, Tokyo, Japan) using cultures grown on Columbia agar with 5% sheep blood (BioMérieux, Marcy l'Etoile, France) at 37 • C for 24 h. A colony was collected from the agar and immersed in a 2.5% glutaraldehyde fixative solution. The slide was gently washed in water, air-dried, and examined with a SU5000 (Hitachi SU5000) operated at 10.0 kV. Subculture of the strain Marseille-Q4567 T was attempted at a wide range of temperatures (25, 28, 31.5, 37, 41.5, and 56 • C) on Columbia agar with 5% sheep blood and in different conditions of pH (5.5, 6.5, 7.5, and 8.5) and salinity (5, 10, and 15g/L) on Columbia agar bases (BioMérieux, Marcy l'Etoile, France). The growth of the strain was also tested under anaerobic (AnaeroGen Compact; Oxoid, Thermo Scientific, Dardilly, France), microaerophilic (campyGEN; Oxoid, Thermo Scientific, Dardilly, France), and aerobic conditions at 37 • C for 48 h. API ZYM and API 50 CH kits (BioMérieux, Marcy l'Etoile, France) were used to perform biochemical analysis in accordance with the manufacturer's instructions. Oxidase (MASTDISCS ® ID, Mast Group Ltd., Bootle, Merseyside, UK) and catalase (BioMérieux, Marcy l'Etoile, France) assays were also performed. The ability of strain Marseille-Q4567 T to form spores was evaluated following thermal shock at 80 • C for 10 min. Cellular fatty acid methyl ester (FAME) analysis was performed by GC/MS. Two samples were prepared with approximately 3 mg of bacterial biomass per tube harvested from several culture plates. Fatty acid methyl esters were prepared as described by Sasser [25]. GC/MS analyses were carried out as previously described [26]. Briefly, fatty acid methyl esters were separated using an Elite 5-MS column and monitored by mass spectrometry (Clarus 500-SQ 8 S, Perkin Elmer, Courtaboeuf, France). A spectral database search was performed using MS Search 2.0, operated with the Standard Reference Database 1A (NIST, Gaithersburg, MD, USA) and the FAMEs mass spectral database (Wiley, Chichester, UK).

Extraction and Genome Sequencing
The genomic DNA (gDNA) of strain Marseille-Q4567 T was extracted in two steps: Mechanical treatment was first performed by acid washed glass beads (G4649-500g Sigma) using a FastPrep-24™ 5G Grinder (mpBio) at maximum speed (6.5) for 90 s. Then, after 30 min lysozyme incubation at 37 • C, DNA was extracted on the EZ1 biorobot (Qiagen) with an EZ1 DNA tissues kit. The elution volume was 50 µL. The gDNA of strain Marseille-Q4567 T was quantified by a Qubit assay with the high sensitivity kit (Life technologies, Carlsbad, CA, USA) to 0.2 ng/µL. Genomic DNA was next sequenced on MiSeq Technology (Illumina Inc, San Diego, CA, USA) with the paired-end strategy, and was barcoded in order to be mixed respectively with 20 other genomic projects prepared with the Nextera XT DNA sample prep kit (Illumina). To prepare the paired-end library, dilution was performed to require 1 ng of each genome as input to prepare the paired end library. The "tagmentation" step fragmented and tagged the DNA. Limited cycle PCR amplification (12 cycles) completed the tag adapters and introduced dual-index barcodes. After purification on AMPure XP beads (Beckman Coulter Inc, Fullerton, CA, USA), the libraries were then normalized on specific beads according to the Nextera XT protocol (Illumina). Normalized libraries were pooled into a single library for sequencing on the MiSeq. The pooled single strand library was loaded onto the reagent cartridge and then onto the instrument along with the flow cell. Automated cluster generation and paired end sequencing with dual index reads were performed in a single 39-h run in 2x250-bp. Information totaling 5.54 Gb was obtained from a 578 K/mm 2 cluster density with a cluster passing quality control filters of 95.5%. Within this run, the index representation for strain Marseille-Q4567 T was determined to index 5.97%. The 11,232,685 paired-end reads were filtered according to the read qualities.
In order to improve the genome sequence, an Oxford Nanopore approach was performed on 1D genomic DNA sequencing for the MinIon device using the SQK-LSK109 kit. A library was constructed from 1 µg genomic DNA without fragmentation or end repair. Adapters were ligated to both ends of genomic DNA. After purification on AMPure XP beads (Beckman Coulter Inc, Fullerton, CA, USA), the library was quantified by a Qubit assay with the high sensitivity kit (Life technologies, Carlsbad, CA, USA). In total, 1500 active pores were detected for the sequencing and the "What's In My Pot (WIMP)" workflow was chosen for bioinformatic analysis in real time. After 21 h as the run time and end life of the flow cell, 177,380 reads as raw data were generated.

Assembly and Annotation of the Genome Sequence
The assembly was performed with a pipeline incorporating different software (Velvet [27], Spades [28], Soap Denovo [29]), and trimmed data (MiSeq and Trimmomatic [30] software) or untrimmed data (only MiSeq software). GapCloser was used to reduce assembly gaps. Scaffolds < 800 bp and scaffolds with a depth value < 25% of the mean depth were removed. The best assembly was selected using different criteria (number of scaffolds, N50, number of N).
Prokka (Galaxy v 1.14.5) was used for prediction in the open reading frame (ORF) with the default settings [31]. Deviations in the sequencing regions predicted by ORFs were excluded. BlastP was used to predict the bacterial proteome (E value of 1e03, coverage of 70% and percent identity of 30%) according to the Orthological Group (COG) database. In the absence of a match, the search for BlastP in the database [32] was extended with an E value of 1e03, coverage of 70%, and percent identity of 30%. If the length of the sequence was less than 80 amino acids (aa), an E value of 1e05 was used. The rRNA and tRNA genes were retrieved using the Prokka (Galaxy v 1.14.5) [33,34]. CGView Server BETA [21] was used to generate a circular map showing a complete view of the genome of strain Marseille-Q4567 T . Finally, a mass screening for antimicrobial and virulence genes was performed using ABRIcate (Galaxy Version 1.0.1) [35].

Phylogenetic Analysis and Genome Comparison
The 16S rRNA gene sequence of strain Marseille-Q4567 T (1532 bp) was obtained and compared with the most closely related species retrieved using NCBI BLAST (National Center for Biotechnology Information, Basic Local Alignment Search Tool; https: //blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 10 January 2021) and then submitted to the GenBank database. Phylogenetic analyses were performed using MEGA X software [22], with genetic distances determined according to the Kimura two-parameter model [36] and phylogenies reconstructed with the maximum-likelihood method. The topology of the phylogenetic tree was conducted using the bootstrap method with 1000 repetitions. All positions containing gaps and missing data were eliminated from the dataset (complete deletion option). Digital DNA-DNA hybridization (dDDH) values between strain Marseille-Q4567 T and other closely related species were assessed using the Genome-to-Genome Distance Calculator 2.1 (GGDC) web server (http://ggdc.dsmz.de, accessed on 10 January 2021) [37]. In addition, the orthologous average nucleotide identity was calculated using OrthoANI v 0.93.1 software [38]. Finally, the pangenome distribution of strain Marseille Q4567 T and other closely related species was evaluated using Roary software (Galaxy v 3.13.0) [39].  Informed Consent Statement: Informed consent was obtained from the subject involved in the study.

Data Availability Statement:
The data presented in this study are contained within the article.