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

Angéline Antezack, Manon Boxberger, Bernard La Scola, Virginie Monnet-Corti

► To cite this version:

Angéline Antezack, Manon Boxberger, Bernard La Scola, Virginie Monnet-Corti. Isolation and Description of Catonella massiliensis sp. nov., a Novel Catonella Species, Isolated from a Stable Periodontitis Subject. Pathogens, 2021, 10 (3), pp.367. 10.3390/pathogens10030367. hal-03243237

HAL Id: hal-03243237

https://amu.hal.science/hal-03243237

Submitted on 31 May 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
Article

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

Angéline Antezack<sup>1,2,3,4</sup>, Manon Boxberger<sup>3,4 ID</sup>, Bernard La Scola<sup>3,4 ID</sup> and Virginie Monnet-Corti<sup>1,2,3,4,* ID</sup>

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

check for
updates

**Citation:** Antezack, A.; Boxberger, M.; La Scola, B.; Monnet-Corti, V. *Isolation and Description of *Catonella massiliensis** sp. nov., a Novel *Catonella* Species, Isolated from a Stable Periodontitis Subject. *Pathogens* **2021**, *10*, 367. https://doi.org/10.3390/pathogens10030367

Academic Editor: Lawrence S. Young

Received: 23 February 2021

Accepted: 16 March 2021

Published: 19 March 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** 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).

**Keywords:** *Catonella*; gingival health; saliva; culturomics; sp. nov.

## 1. 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–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

<footnote n="1">¹ Ecole de Médecine Dentaire, Faculté des Sciences Médicales et Paramédicales, Aix-Marseille Université, 27 boulevard Jean Moulin, 13385 Marseille, France; angeline.antezack@univ-amu.fr</footnote>

<footnote n="2">² Assistance Publique-Hôpitaux de Marseille (AP-HM), Hôpital Timone, Service de Parodontologie, 264, rue Saint Pierre, 13385 Marseille, France</footnote>

<footnote n="3">³ Institut de Recherche pour le Développement (IRD), Aix-Marseille Université, Assistance Publique-Hôpitaux de Marseille (AP-HM), MEPHI, 27 boulevard Jean Moulin, 13005 Marseille, France; manon.boxberger@hotmail.fr (M.B.); bernard.la-scola@univ-amu.fr (B.L.S.)</footnote>

<footnote n="4">⁴ IHU Méditerranée Infection, 19–21 boulevard Jean Moulin, 13005 Marseille, France</footnote>

<footnote>* Correspondence: virginie.corti@univ-amu.fr</footnote>
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-Q4567T 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-Q4567T to its closely related phylogenetic neighbors and proposed that strain Marseille-Q4567T represents a novel species of the genus *Catonella*, for which the name *Catonella massiliensis* sp. nov. is proposed (=CSUR Q4567).

**2. Results**

**2.1. Strain Isolation and Phenotypic Characteristics**

Strain Marseille-Q4567T 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-Q4567T could not be identified by MALDI-TOF MS, as the score was lower than 1.8 (Figure 1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** MALDI-TOF MS reference mass spectrum for strain Marseille-Q4567T. Spectra from 12 individual colonies were compared and a reference spectrum was generated.</figcaption></figure>
Growth was observed on Columbia agar with 5% sheep blood (BioMérieux, Marcy-l'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).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Micrograph electron microscopy of strain Marseille-Q4567T. Scales and acquisition settings are shown in the figure.</figcaption></figure>

Using an API 50 CH strip, positive reactions were obtained for glycerol, erythritol, D-arabinose, 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-acetyl-glucosamine, 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 5-ketogluconate. 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 Q4567T was negative for oxidase and catalase activities. Comparison of phenotypic characteristics between strain Marseille-Q4567T 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.
<table><caption>**Table 1.** Phenotypic characterization of *Catonella massiliensis* sp. nov., compared to other closely related species. +, Positive; w, weak; -, negative; ND, no data available. Strains: 1, Marseille-Q4567T [17]; 2, *Cuneatibacter caecimuris* DSM 29486T [17]; 3, *Faecalicatena orotica* JCM 1429T [18]; 4, *Faecalicatena contorta* JCM 6483T [18]; 5, *Faecalicatena fissicatena* JCM 31501T [18]; 6, *Anaerocolumna cellulosilytica* strain SN021T [19]; 7, *Catonella morbi* ATCC 51271T [13]; 8, *Herbinix luporum* strain SD1D* [20].</caption><thead><tr><th>Characteristic</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th></tr></thead><tbody><tr><td>Gram stain</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>Spore formation</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td></tr><tr><td>Temperature (°C)</td><td>28-41.5</td><td>37</td><td>20-40</td><td>20-40</td><td>20-40</td><td>15-40</td><td>37</td><td>40-65</td></tr><tr><td>pH</td><td>6.5-8.5</td><td>7.0</td><td>6.0-8.0</td><td>6.0-8.0</td><td>6.0-8.0</td><td>6.2-8.5</td><td>5.0-5.6</td><td>6.5-8.5</td></tr><tr><td>Catalase activity</td><td>-</td><td>ND</td><td>ND</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Oxidase activity</td><td>-</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>-</td><td>ND</td><td>ND</td></tr><tr><td>Acid produced from:</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>L-Arabinose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>Cellobiose</td><td>+</td><td>-</td><td>+</td><td>+</td><td>-</td><td>+</td><td>+</td><td>+</td></tr><tr><td>Lactose</td><td>+</td><td>-</td><td>+</td><td>+</td><td>-</td><td>+</td><td>+</td><td>-</td></tr><tr><td>Maltose</td><td>+</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>ND</td></tr><tr><td>D-Mannitol</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>ND</td></tr><tr><td>D-Mannose</td><td>+</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>Raffinose</td><td>+</td><td>-</td><td>+</td><td>+</td><td>-</td><td>-</td><td>+</td><td>ND</td></tr><tr><td>L-Rhamnose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td><td>ND</td></tr><tr><td>D-Sorbitol</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Enzyme activity:</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Alkaline phosphatase</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>4</sub> esterase</td><td>-</td><td>ND</td><td>w</td><td>w</td><td>w</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>8</sub> esterase lipase</td><td>-</td><td>ND</td><td>w</td><td>w</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>14</sub> lipase</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td></tr><tr><td>Leucine arylamidase</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>Valine arylamidase</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>Cystine arylamidase</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>Trypsin</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>α-Chymotrypsin</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>Acid phosphatase</td><td>+</td><td>ND</td><td>+</td><td>+</td><td>w</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>Naphthol-AS-BI-phosphohydrolase</td><td>+</td><td>ND</td><td>w</td><td>w</td><td>w</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>α-Galactosidase</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>β-Galactosidase</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>β-Glucuronidase</td><td>-</td><td>+</td><td>+</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>α-Glucosidase</td><td>+</td><td>+</td><td>+</td><td>w</td><td>+</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>β-Glucosidase</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>N-Acetyl-β-glucosaminidase</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>α-Mannosidase</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>α-Fucosidase</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr></tbody></table>

The most abundant fatty acid was C<sub>16:0</sub> (64.2%), followed by unsaturated structures C<sub>18:1n9</sub> (12.5%) and C<sub>18:2n6</sub> (7.8%). No branched fatty acids were described. Minor amounts of other unsaturated and saturated structures were also detected (Table 2).

## 2.2. Genome Sequencing Information and Genome Properties

The genome size of strain Marseille-Q4567T 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-Q4567T (Table 3). Finally, 1396 genes (49%) were annotated as hypothetical proteins for strain Marseille-Q4567T. A circular map showing a complete view of the genome of strain Marseille-Q4567T 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.

<table><caption>**Table 2.** Cellular fatty acid compositions of strain Marseille-Q4567T and its closely related species. DMA, dimethyl acetal; ALDE, aldehyde. -, not detected; TR, trace amounts (<1%); ND, no data available. Strains: 1, Marseille-Q4567T; 2, *Cuneatibacter caecimuris* DSM 29486T [17]; 3, *Faecalicatena orotica* JCM 1429T [18]; 4, *Faecalicatena contorta* JCM 6483T [18]; 5, *Faecalicatena fissicatena* JCM 31501T [18]; 6, *Anaerocolumna cellulosilytica* strain SN021T [19]; 7, *Catonella morbi* ATCC 51271T [13]; 8, *Herbinix luporum* strain SD1D T [20].</caption><thead><tr><th>Fatty Acid</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th></tr></thead><tbody><tr><td>C<sub>12:0</sub></td><td>TR</td><td>ND</td><td>1.6</td><td>3.3</td><td>2.6</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>13:0</sub></td><td>TR</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>14:0</sub></td><td>6.4</td><td>6.2</td><td>6.0</td><td>9.3</td><td>9.2</td><td>0.6</td><td>42</td><td>14.0</td></tr><tr><td>C<sub>14:0</sub> DMA</td><td>ND</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>1.5</td><td>14</td><td>9.1</td></tr><tr><td>C<sub>15:0</sub></td><td>TR</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>16:0</sub></td><td>64.2</td><td>43.4</td><td>12.0</td><td>15.6</td><td>15.6</td><td>5.5</td><td>12</td><td>19.9</td></tr><tr><td>C<sub>16:0</sub> ALDE</td><td>-</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>20.6</td><td>ND</td><td>2.9</td></tr><tr><td>C<sub>16:0</sub> DMA</td><td>ND</td><td>ND</td><td>2.0</td><td>3.4</td><td>3.4</td><td>17.4</td><td>ND</td><td>5.3</td></tr><tr><td>C<sub>16:0</sub> 3-OH</td><td>-</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>1.2</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>16:1</sub>ω<sub>9c</sub> DMA</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>3.8</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>16:1n9</sub></td><td>-</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>2</td><td>1.2</td></tr><tr><td>C<sub>17:0</sub></td><td>TR</td><td>ND</td><td>-</td><td>-</td><td>-</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>17:1</sub>ω<sub>8c</sub></td><td>ND</td><td>8.3</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>18:0</sub></td><td>6.6</td><td>ND</td><td>1.4</td><td>-</td><td>3.9</td><td>ND</td><td>4</td><td>1.2</td></tr><tr><td>C<sub>18:1n7</sub></td><td>1.3</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>18:1n9</sub></td><td>12.5</td><td>ND</td><td>25.9</td><td>21.9</td><td>22.9</td><td>ND</td><td>5</td><td>ND</td></tr><tr><td>C<sub>18:1</sub>ω<sub>9c</sub> DMA</td><td>ND</td><td>ND</td><td>37.0</td><td>36.3</td><td>28.6</td><td>3.1</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>18:1</sub>ω<sub>7c</sub> DMA</td><td>ND</td><td>ND</td><td>2.4</td><td>-</td><td>4.0</td><td>9.0</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>18:2n6</sub></td><td>7.8</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td></tr><tr><td>C<sub>19:0</sub> cyclo 9,10 DMA</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>38.3</td></tr><tr><td>C<sub>19:0</sub> cyclo 11,12 DMA</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>ND</td><td>1.4</td><td>ND</td><td>ND</td></tr></tbody></table>

<table><caption>**Table 3.** Number of genes associated with the general clusters of orthologous group (COG) functional categories of strain Marseille-Q4567T.</caption><thead><tr><th>Code</th><th>Strain Marseille-Q4567 T</th><th>Description</th></tr></thead><tbody><tr><td>[J]</td><td>145</td><td>Translation, ribosomal structure, and biogenesis</td></tr><tr><td>[A]</td><td>0</td><td>RNA processing and modification</td></tr><tr><td>[K]</td><td>150</td><td>Transcription</td></tr><tr><td>[L]</td><td>137</td><td>Replication, recombination, and repair</td></tr><tr><td>[B]</td><td>0</td><td>Chromatin structure and dynamics</td></tr><tr><td>[D]</td><td>29</td><td>Cell cycle control, cell division, chromosome partitioning</td></tr><tr><td>[Y]</td><td>0</td><td>Nuclear structure</td></tr><tr><td>[V]</td><td>76</td><td>Defense mechanisms</td></tr><tr><td>[T]</td><td>82</td><td>Signal transduction mechanisms</td></tr><tr><td>[M]</td><td>64</td><td>Cell wall/membrane/envelope biogenesis</td></tr><tr><td>[N]</td><td>49</td><td>Cell motility</td></tr><tr><td>[Z]</td><td>0</td><td>Cytoskeleton</td></tr><tr><td>[W]</td><td>0</td><td>Extracellular structures</td></tr><tr><td>[U]</td><td>38</td><td>Intracellular trafficking, secretion, and vesicular transport</td></tr><tr><td>[O]</td><td>60</td><td>Posttranslational modification, protein turnover, chaperones</td></tr><tr><td>[X]</td><td>0</td><td>Mobilome: Prophages, transposons</td></tr><tr><td>[C]</td><td>74</td><td>Energy production and conversion</td></tr><tr><td>[G]</td><td>183</td><td>Carbohydrate transport and metabolism</td></tr><tr><td>[E]</td><td>141</td><td>Amino acid transport and metabolism</td></tr><tr><td>[F]</td><td>51</td><td>Nucleotide transport and metabolism</td></tr><tr><td>[H]</td><td>31</td><td>Coenzyme transport and metabolism</td></tr><tr><td>[I]</td><td>47</td><td>Lipid transport and metabolism</td></tr><tr><td>[P]</td><td>86</td><td>Inorganic ion transport and metabolism</td></tr><tr><td>[Q]</td><td>21</td><td>Secondary metabolites biosynthesis, transport, and catabolism</td></tr><tr><td>[R]</td><td>214</td><td>General function prediction only</td></tr><tr><td>[S]</td><td>134</td><td>Function unknown</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** A circular map generated using the CGView Server<sup>BETA</sup> [21] showing a complete view of the genome of strain Marseille-Q45657<sup>T</sup>.</figcaption></figure>

## 2.3. *Comparison with Closely Related Bacterial Strains*

16S rDNA-based similarity analysis of strain Marseille-Q4567T against GenBank yielded the highest nucleotide sequence similarities of 98.23% sequence identity with *C. morbi* ATCC 51271T (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.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Maximum likelihood tree based on the comparison of 16S rRNA gene sequences showing the phylogenetic relationships of *Catonella massiliensis* Marseille-Q4567T and other closely related species. Bootstrap values (expressed as percentages of 1000 replications) are displayed at the nodes. 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.</figcaption></figure>
The genome of strain Marseille-Q4567T 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 Q4567T (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><caption>**Table 4.** Genomic comparison (sequence size, number of contigs, G+C contents, and gene content) between strain Marseille-Q4567T and nine closely related bacterial strains.</caption><thead><tr><th>Species</th><th>Genome Size (bp)</th><th>Number of Contigs</th><th>G+C Content (%)</th><th>Gene Content</th></tr></thead><tbody><tr><td><i>Catonella massiliensis</i><br>Marseille-Q4567T</td><td>3,122,925</td><td>3</td><td>38.8</td><td>2849</td></tr><tr><td><i>Cuneatibacter caecimuris</i> DSM<br>29486T (NZ_SGF01000001.1)</td><td>3,462,725</td><td>16</td><td>49.1</td><td>3268</td></tr><tr><td><i>Eubacterium oxidoreducens</i> DSM<br>3217T (NZ_FMXR01000031.1)</td><td>2,912,287</td><td>33</td><td>39.8</td><td>2700</td></tr><tr><td><i>Faecalicatena orotica</i> strain<br>NLAE-zl-C242<br>(NZ_QGDL01000001.1)</td><td>5,717,637</td><td>33</td><td>44.7</td><td>5174</td></tr><tr><td><i>Roseburia hominis</i> A2-183T<br>(NC_015977.1)</td><td>3,592,125</td><td>1</td><td>48.5</td><td>3349</td></tr><tr><td><i>Faecalicatena contorta</i><br>2789STDY5834876<br>(NZ_CYZU01000001.1)</td><td>5,545,490</td><td>139</td><td>45.9</td><td>5148</td></tr><tr><td><i>Faecalicatena fissicatena</i> KCTC<br>15010 (NZ_LDAQ01000001.1)</td><td>5,014,239</td><td>184</td><td>45.6</td><td>4411</td></tr><tr><td><i>Anaerocolumna cellulosilytica</i><br>strain SN021T<br>(NZ_AP023367.1)</td><td>5,430,627</td><td>1</td><td>36.7</td><td>4545</td></tr><tr><td><i>Catonella morbi</i> ATCC 51271T<br>(NZ_KI535366.1)</td><td>3,479,204</td><td>8</td><td>37.1</td><td>3155</td></tr><tr><td><i>Herbinix luporum</i> strain SD1D T<br>(NZ_LN879430.1)</td><td>2,609,352</td><td>1</td><td>35.3</td><td>2430</td></tr></tbody></table>

Using dDDH analysis, strain Marseille-Q4567T exhibited values ranging from 22.9% [20.6–25.3%] with *H. luporum* to 28.4% [26.1–30.9%] with *C. caecimuris* (Table 5). These values were lower than the 70% threshold used for delineating prokaryotic species, thus confirming that strain Marseille-Q4567T represents a new species [10]. In addition, using OrthoANI analysis, strain Marseille-Q4567T exhibited values ranging from 64.97% with *C. caecimuris* to 79.43% with *C. morbi* (Figure 5).

<table><caption>**Table 5.** Numerical DNA-DNA hybridization values (%) obtained by comparison between strain Marseille-Q4567T and other closely related species using GGDC formula 2 software (DDH estimates based on HSP identities/length) [https://ggdc.dsmz.de/ggdc.php#, accessed on 10 January 2021]. Type strains are indicated with superscript T.</caption><thead><tr><th>Species</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>1 <i>Catonella massiliensis</i> Marseille-Q4567T</td><td>100.00</td><td>28.4<br>[26.1-30.9]</td><td>27.0<br>[24.6-29.5]</td><td>26.0<br>[23.7-28.5]</td><td>24.8<br>[22.5-27.3]</td><td>24.7<br>[22.3-27.1]</td><td>24.60<br>[22.3-27.1]</td><td>24.5<br>[22.2-27]</td><td>23.80<br>[21.5-26.2]</td><td>22.9<br>[20.6-25.3]</td></tr><tr><td>2 <i>Cuneatibacter caecimuris</i> DSM 29486T<br>(NZ_SGF01000001.1)</td><td></td><td>100.00</td><td>25.7<br>[23.3-28.1]</td><td>23.0<br>[20.7-25.4]</td><td>25.8<br>[23.5-28.3]</td><td>19.2<br>[17-21.6]</td><td>19.0<br>[16.8-21.4]</td><td>25.4<br>[23.1-27.9]</td><td>23.1<br>[20.8-25.5]</td><td>30.1<br>[27.7-32.6]</td></tr><tr><td>3 <i>Eubacterium oxidoreducens</i> DSM 3217T<br>(NZ_FMXR01000031.1)</td><td></td><td></td><td>100.00</td><td>22.6<br>[20.3-25]</td><td>23.2<br>[20.9-25.7]</td><td>26.1<br>[23.8-28.6]</td><td>22.9<br>[20.6-25.3]</td><td>28.8<br>[26.4-31.3]</td><td>29.3<br>[26.9-31.8]</td><td>20.7<br>[18.4-23.1]</td></tr><tr><td>4 <i>Faecalicatena orotica</i> strain NLAE-zl-C242<br>(NZ_QGDL01000001.1)</td><td></td><td></td><td></td><td>100.00</td><td>22.0<br>[19.7-24.4]</td><td>21.6<br>[19.4-24.1]</td><td>21.7<br>[19.4-24.1]</td><td>22.2<br>[20-24.7]</td><td>22.2<br>[19.9-24.6]</td><td>25.5<br>[23.2-28]</td></tr><tr><td>5 <i>Roseburia hominis</i> A2-183T (NC_015977.1)</td><td></td><td></td><td></td><td></td><td>100.00</td><td>24.0<br>[21.7-26.5]</td><td>24.6<br>[22.3-27.1]</td><td>32.9<br>[30.5-35.4]</td><td>23.70<br>[21.4-26.1]</td><td>29.2<br>[26.9-31.7]</td></tr><tr><td>6 <i>Faecalicatena contorta</i> 2789STDY5834876<br>(NZ_CYZU01000001.1)</td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>28.5<br>[26.2-31]</td><td>28.5<br>[26.1-31]</td><td>26.4<br>[24-28.8]</td><td>25.8<br>[23.5-28.3]</td></tr><tr><td>7 <i>Faecalicatena fissicatena</i> KCTC 15010<br>(NZ_LDAQ01000001.1)</td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>26.4<br>[24.1-28.9]</td><td>25.3<br>[22.9-27.8]</td><td>23.1<br>[20.9-25.6]</td></tr><tr><td>8 <i>Anaerocolumna cellulosilytica</i> strain SN021T<br>(NZ_AP023367.1)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>24.1<br>[21.8-26.5]</td><td>26.5<br>[24.1-29]</td></tr><tr><td>9 <i>Catonella morbi</i> ATCC 51271T (NZ_KI535366.1)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>21.0<br>[18.8-23.4]</td></tr><tr><td>10 <i>Herbinix luporum</i> strain SD1D T<br>(NZ_LN879430.1)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td></tr></tbody></table>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Heatmap generated with orthologous average nucleotide identity (OrthoANI) values calculated using the OAT software for strain Marseille-Q4567T with other closely related species validly described.</figcaption></figure>

Finally, pangenome analysis of strain Marseille-Q4567T showed a total of 34,912 gene clusters distributed as follows: (Core genes = 0), (soft core genes = 0), (shell genes = 1171), and (cloud genes = 33,741), respectively (Figure 6).

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Pangenome analysis strain of Marseille-Q4567T whole-genome sequences. A maximum likelihood tree was constructed from the accessory genome elements (left). The presence (blue) and the absence (white) of accessory genome elements are presented on the right.</figcaption></figure>
## 2.4. Description of *Catonella Massiliensis* sp. nov.

*Catonella massiliensis* (mas.si.li.en'sis L. masc./fem. adj. massiliensis, pertaining to Massilia, the ancient Roman name of Marseille, France, where the organism was isolated).

Cells were Gram-negative, anaerobic, curved rods, round-tipped, and approximatively 0.3 to 0.5 µm wide and 3.9 to 4.1 µm long. Colonies on Columbia agar with 5% sheep blood (BioMérieux, Marcy-l'Etoile, France) incubated at 37 °C for 2 days were approximatively 0.1 to 0.5 mm in diameter, smooth, circular, and translucent. The temperature range for growth was 28-41.5 °C (optimum 37 °C). The cells grew at pH 6.5 to 8.5 (optimum 7.5) and with an NaCl content ranging from 5 to 10 g/L (optimum 5 g/L). Test results were positive for the fermentation of glycerol, erythritol, D-arabinose, 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-acetyl-glucosamine, amylgdalin, 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 5-ketogluconate. According to the API ZYM system, cells were positive for alkaline phosphatase, leucine arylamidase, acid phosphatase, naphthol-AS-BI-phosphohydrolase, α-galactosidase, β-galactosidase, α-glucosidase, β-glucosidase, and N-acetyl-β-glucosaminidase. The predominant cellular fatty acid was $C_{16:0}$ (64.2%), followed by unsaturated structures $C_{18:1n9}$ (12.5%) and $C_{18:2n6}$ (7.8%). The genome size of strain Marseille-Q4567T was 3,122,925 bp long with a 38.8% G+C content. The type of strain, Marseille-Q4567T (CSUR Q4567), was isolated from a sample of unstimulated whole saliva from a male with stable periodontitis. The sequence data of the 16S rRNA gene of Marseille-Q4567T and the whole genome have been deposited in the GenBank database under accession numbers MW410925 and JAEPJ000000000, respectively.

## 3. 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-Q4567T 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-Q4567T shared the highest 16S rRNA gene sequence similarities (98.23%) with *C. morbi* ATCC 51271T, 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-Q4567T 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-Q4567T 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-Q4567T 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-Q4567T represents a novel species of the genus *Catonella*, for which the name *Catonella massiliensis* sp. nov. is proposed (=CSUR Q4567).

# 4. Materials and Methods

## 4.1. *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™ 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-Q4567T 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-Q4567T 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-Q4567T 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).

## 4.2. Extraction and Genome Sequencing

The genomic DNA (gDNA) of strain Marseille-Q4567T 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-Q4567T 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² cluster density with a cluster passing quality control filters of 95.5%. Within this run, the index representation for strain Marseille-Q4567T 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.

## 4.3. 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 ServerBETA [21] was used to generate a circular map showing a complete view of the genome of strain Marseille-Q4567T. Finally, a mass screening for antimicrobial and virulence genes was performed using ABRIcate (Galaxy Version 1.0.1) [35].

4.4. Phylogenetic Analysis and Genome Comparison

The 16S rRNA gene sequence of strain Marseille-Q4567T (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-Q4567T 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 Q4567T and other closely related species was evaluated using Roary software (Galaxy v 3.13.0) [39].

**Author Contributions:** Conceptualization, B.L.S. and V.M.-C.; methodology, A.A. and M.B.; software, A.A.; validation, B.L.S. and V.M.-C.; formal analysis, A.A.; investigation, A.A.; resources, A.A., B.L.S. and V.M.-C.; data curation, A.A.; writing—original draft preparation, A.A.; writing—review and editing, A.A., M.B., B.L.S., and V.M.-C.; visualization, A.A., B.L.S., and V.M.-C.; supervision, B.L.S. and V.M.-C.; project administration, B.L.S. and V.M.-C.; funding acquisition, B.L.S. and V.M.-C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the French Government under the "Investissements d'avenir" (Investments for the Future) program managed by the Agence Nationale de la Recherche (ANR, French National Agency for Research) (reference: Méditerranée Infection 10-IAHU-03).

**Institutional Review Board Statement:** The study was conducted in accordance with the guidelines of the Declaration of Helsinki and approved by the Comité de Protection des Personnes (C.P.P.) Sud-Ouest et Outre-Mer 1 (n° ID RCB: 2020-A01234-35-CPP 1-20-075 ID 9806).

**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.

**Acknowledgments:** The authors thank Nicholas Armstrong for performing the fatty acid analyses and Anthony Fontanini for the electron microscopy observation.

**Conflicts of Interest:** Boxberger M. received a PhD grant supported by L'Occitane Society. All other authors have no conflict of interest to declare.

**References**

1. Deo, P.N.; Deshmukh, R. Oral microbiome: Unveiling the fundamentals. *J. Oral Maxillofac. Pathol.* **2019**, *23*, 122–128. [CrossRef]

2. Socransky, S.S.; Haffajee, A.D. Periodontal microbial ecology. *Periodontol. 2000* **2005**, *38*, 135–187. [CrossRef]

3. Curtis, M.A.; Diaz, P.I.; Van Dyke, T.E. The role of the microbiota in periodontal disease. *Periodontol. 2000* **2020**, *83*, 14–25. [CrossRef]

4. Nyvad, B.; Takahashi, N. Integrated hypothesis of dental caries and periodontal diseases. *J. Oral Microbiol.* **2020**, *12*, 1710953. [CrossRef] [PubMed]

5. Ye, S.H.; Siddle, K.J.; Park, D.J.; Sabeti, P.C. Benchmarking Metagenomics Tools for Taxonomic Classification. *Cell* **2019**, *178*, 779–794. [CrossRef]
6. Baker, J.L.; Morton, J.T.; Dinis, M.; Alvarez, R.; Tran, N.C.; Knight, R.; Edlund, A. Deep metagenomics examines the oral microbiome during dental caries, revealing novel taxa and co-occurrences with host molecules. *Genome Res.* **2021**, *31*, 64–74. [CrossRef] [PubMed]

7. Qi, Y.; Zang, S.-Q.; Wei, J.; Yu, H.-C.; Yang, Z.; Wu, H.-M.; Kang, Y.; Tao, H.; Yang, M.-F.; Jin, L.; et al. High-throughput sequencing provides insights into oral microbiota dysbiosis in association with inflammatory bowel disease. *Genomics* **2021**, *113*, 664–676. [CrossRef]

8. Chen, T.; Yu, W.-H.; Izard, J.; Baranova, O.V.; Lakshmanan, A.; Dewhirst, F.E. The Human Oral Microbiome Database: A web accessible resource for investigating oral microbe taxonomic and genomic information. *Database* **2010**, *2010*, baq013. [CrossRef] [PubMed]

9. Wade, W.; Thompson, H.; Rybalka, A.; Vartoukian, S. Uncultured Members of the Oral Microbiome. *J. Calif. Dent. Assoc.* **2016**, *44*, 447–456. [PubMed]

10. Vartoukian, S.R.; Palmer, R.M.; Wade, W.G. Strategies for culture of 'unculturable' bacteria: Culturing the unculturable. *FEMS Microbiol. Lett.* **2010**, *309*, 1–7. [CrossRef]

11. Papapanou, P.N.; Sanz, M.; Buduneli, N.; Dietrich, T.; Feres, M.; Fine, D.H.; Flemmig, T.F.; Garcia, R.; Giannobile, W.V.; Graziani, F.; et al. Periodontitis: Consensus report of workgroup 2 of the 2017 World Workshop on the Classification of Periodontal and Peri-Implant Diseases and Conditions. *J. Clin. Periodontol.* **2018**, *45*, S162–S170. [CrossRef]

12. Chapple, I.L.C.; Mealey, B.L.; Van Dyke, T.E.; Bartold, P.M.; Dommisch, H.; Eickholz, P.; Geisinger, M.L.; Genco, R.J.; Glogauer, M.; Goldstein, M.; et al. Periodontal health and gingival diseases and conditions on an intact and a reduced periodontium: Consensus report of workgroup 1 of the 2017 World Workshop on the Classification of Periodontal and Peri-Implant Diseases and Conditions. *J. Periodontol.* **2018**, *89*, S74–S84. [CrossRef] [PubMed]

13. Moore, L.V.; Moore, W.E. *Oribaculum catoniae* gen. nov., sp. nov.; *Catonella morbi* gen. nov., sp. nov.; *Hallella seregens* gen. nov., sp. nov.; *Johnsonella ignava* gen. nov., sp. nov.; and *Dialister pneumosintes* gen. nov., comb. nov., nom. rev., Anaerobic Gram-Negative Bacilli from the Human Gingival Crevice. *Int. J. Syst. Bacteriol.* **1994**, *44*, 187–192. [CrossRef] [PubMed]

14. Wu, X.; Chen, J.; Xu, W.; Zhu, D.; Wang, X.; Chen, Y.; Wu, J.; Cui, C.; Zhang, W.; Yu, L. 16S rDNA analysis of periodontal plaque in chronic obstructive pulmonary disease and periodontitis patients. *J. Oral Microbiol.* **2017**, *9*, 1324725. [CrossRef] [PubMed]

15. Colombo, A.P.V.; Bennet, S.; Cotton, S.L.; Goodson, J.M.; Kent, R.; Haffajee, A.D.; Socransky, S.S.; Hasturk, H.; Van Dyke, T.E.; Dewhirst, F.E.; et al. Impact of Periodontal Therapy on the Subgingival Microbiota of Severe Periodontitis: Comparison Between Good Responders and Individuals With Refractory Periodontitis Using the Human Oral Microbe Identification Microarray. *J Periodontol.* **2012**, *83*, 1279–1287. [CrossRef] [PubMed]

16. Siqueira, J.F.; Rôças, I.N. *Catonella morbi* and *Granulicatella adiacens*: New species in endodontic infections. *Oral Surg. Oral Med. Oral Pathol. Oral Radiol. Endodontol.* **2006**, *102*, 259–264. [CrossRef] [PubMed]

17. Lagkouvardos, I.; Pukall, R.; Abt, B.; Foesel, B.U.; Meier-Kolthoff, J.P.; Kumar, N.; Bresciani, A.; Martínez, I.; Just, S.; Ziegler, C.; et al. The Mouse Intestinal Bacterial Collection (miBC) provides host-specific insight into cultured diversity and functional potential of the gut microbiota. *Nat. Microbiol.* **2016**, *1*, 16131. [CrossRef]

18. Sakamoto, M.; Iino, T.; Ohkuma, M. *Faecalimonas umbilicata* gen. nov., sp. nov., isolated from human faeces, and reclassification of *Eubacterium contortum*, *Eubacterium fissicatena* and *Clostridium oroticum* as *Faecalicatena contorta* gen. nov., comb. nov., *Faecalicatena fissicatena* comb. nov. and *Faecalicatena orotica* comb. nov. *Int. J. Syst. Evol. Microbiol.* **2017**, *67*, 1219–1227. [CrossRef]

19. Ueki, A.; Ohtaki, Y.; Kaku, N.; Ueki, K. Descriptions of *Anaerotaenia torta* gen. nov., sp. nov. and *Anaerocolumna cellulosilytica* gen. nov., sp. nov. isolated from a methanogenic reactor of cattle waste and reclassification of *Clostridium aminovalericicum*, *Clostridium jejuense* and *Clostridium xylanovorans* as *Anaerocolumna* species. *Int. J. Syst. Evol. Microbiol.* **2016**, *66*, 2936–2943. [CrossRef]

20. Koeck, D.E.; Hahnke, S.; Zverlov, V.V. *Herbinix luporum* sp. nov., a thermophilic cellulose-degrading bacterium isolated from a thermophilic biogas reactor. *Int. J. Syst. Evol. Microbiol.* **2016**, *66*, 4132–4137. [CrossRef]

21. Stothard, P.; Wishart, D.S. Circular genome visualization and exploration using CGView. *Bioinformatics* **2004**, *21*, 537–539. [CrossRef]

22. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549. [CrossRef]

23. Seng, P.; Drancourt, M.; Gouriet, F.; La Scola, B.; Fournier, P.; Rolain, J.M.; Raoult, D. Ongoing Revolution in Bacteriology: Routine Identification of Bacteria by Matrix-Assisted Laser Desorption Ionization Time-of-Flight Mass Spectrometry. *Clin. Infect. Dis.* **2009**, *49*, 543–551. [CrossRef]

24. Seng, P.; Abat, C.; Rolain, J.M.; Colson, P.; Lagier, J.-C.; Gouriet, F.; Fournier, P.E.; Drancourt, M.; La Scola, B.; Raoult, D. Identification of Rare Pathogenic Bacteria in a Clinical Microbiology Laboratory: Impact of Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry. *J. Clin. Microbiol.* **2013**, *51*, 2182–2194. [CrossRef]

25. Sasser, M. *Bacterial Identification by Gas Chromatographic Analysis of Fatty Acid Methyl Esters (GC-FAME)*; Microbial ID: Newark, NJ, USA, 2006.

26. Dione, N.; Sankar, S.; Lagier, J.-C.; Khelaifia, S.; Michele, C.; Armstrong, N.; Richez, M.; Abrahão, J.; Raoult, D.; Fournier, P.-E. Genome sequence and description of *Anaerosalibacter massiliensis* sp. nov. *New Microbes New Infect.* **2016**, *10*, 66–76. [CrossRef] [PubMed]

27. Zerbino, D.R.; Birney, E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. *Genome Res.* **2008**, *18*, 821–829. [CrossRef] [PubMed]
28. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. *J. Comput. Biol.* **2012**, *19*, 455–477. [CrossRef] [PubMed]

29. Luo, R.; Liu, B.; Xie, Y.; Li, Z.; Huang, W.; Yuan, J.; He, G.; Chen, Y.; Pan, Q.; Liu, Y.; et al. SOAPdenovo2: An empirically improved memory-efficient short-read de novo assembler. *GigaScience* **2012**, *1*, 18. [CrossRef] [PubMed]

30. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics* **2014**, *30*, 2114–2120. [CrossRef]

31. Hyatt, D.; Chen, G.-L.; LoCascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. *BMC Bioinform.* **2010**, *11*, 119. [CrossRef]

32. Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. *Nucleic Acids Res.* **2016**, *44*, D67–D72. [CrossRef] [PubMed]

33. Seemann, T. Prokka: Rapid Prokaryotic Genome Annotation. *Bioinformatics* **2014**, *30*, 2068–2069. [CrossRef]

34. Cuccuru, G.; Orsini, M.; Pinna, A.; Sbardellati, A.; Soranzo, N.; Travaglione, A.; Uva, P.; Zanetti, G.; Fotia, G. Orione, a web-based framework for NGS analysis in microbiology. *Bioinformatics* **2014**, *30*, 1928–1929. [CrossRef]

35. Seemann, T. ABRicate: Mass Screening of Contigs for Antiobiotic Resistance Genes. Available online: https://github.com/tseemann/abricate (accessed on 10 January 2021).

36. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. *J. Mol. Evol.* **1980**, *16*, 111–120. [CrossRef]

37. Meier-Kolthoff, J.P.; Auch, A.F.; Klenk, H.-P.; Göker, M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. *BMC Bioinformatics* **2013**, *14*, 60. [CrossRef]

38. Lee, I.; Kim, Y.O.; Park, S.-C.; Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. *Int. J. Syst. Evol. Microbiol.* **2016**, *66*, 1100–1103. [CrossRef]

39. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. *Bioinformatics* **2015**, *31*, 3691–3693. [CrossRef]