# Colibacter massiliensis gen. nov. sp. nov., a novel Gram-stain-positive anaerobic diplococcal bacterium, isolated from the human left colon

Hussein Anani, Rita Abou Abdallah, May Khoder, Anthony Fontanini, Morgane Mailhe, Davide Ricaboni, Didier Raoult, Pierre-Edouard Fournier

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

## To cite this version:

Hussein Anani, Rita Abou Abdallah, May Khoder, Anthony Fontanini, Morgane Mailhe, et al.. Colibacter massiliensis gen. nov. sp. nov., a novel Gram-stain-positive anaerobic diplococcal bacterium, isolated from the human left colon. Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-53791-1. hal-02458106

HAL Id: hal-02458106

https://amu.hal.science/hal-02458106v1

Submitted on 5 Jul 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# *Colibacter massiliensis* gen. nov. sp. nov., a novel Gram-stain-positive anaerobic diplococcal bacterium, isolated from the human left colon

Hussein Anani<fnref n="1" />,², Rita Abou Abdallah<fnref n="1" />,², May Khoder²,³, Anthony Fontanini²,³, Morgane Mailhe²,³, Davide Ricaboni²,³, Didier Raoult²,³,⁴ & Pierre-Edouard Fournier <fnref n="1" />,²*

The gut microbiota is considered to play a key role in human health. As a consequence, deciphering its microbial diversity is mandatory. A polyphasic taxonogenomic strategy based on the combination of phenotypic and genomic analyses was used to characterize a new bacterium, strain Marseille-P2911. This strain was isolated from a left colon sample of a 60-year old man who underwent a colonoscopy for an etiological investigation of iron-deficiency anemia in Marseille, France. On the basis of 16S rRNA sequence comparison, the closest phylogenetic neighbor was *Anaeroglobus geminatus* (94.59% 16S rRNA gene sequence similarity) within the family Veillonellaceae. Cells were anaerobic, Gram-stain-positive, non-spore-forming, catalase/oxidase negative cocci grouped in pairs. The bacterium was able to grow at 37 °C after 2 days of incubation. Strain Marseille-P2911 exhibited a genome size of 1,715,864-bp with a 50.2% G + C content, and digital DNA-DNA hybridization (dDDH) and OrthoANI values with A. geminatus of only 19.1 ± 4.5% and 74.42%, respectively. The latter value being lower than the threshold for genus delineation (80.5%), we propose the creation of the new genus *Colibacter* gen. nov., with strain Marseille-P2911T (=DSM 103304 = CSUR P2911) being the type strain of the new species *Colibacter massiliensis* gen. nov., sp. nov.

With 10¹<fnref n="1" /> to 10¹² cells per gram<fnref n="1" />, the commensal microbiota that resides in the colon has become a focus of interest that is attracting the attention of scientists. Although this complex flora plays a role in homeostasis, it has been demonstrated to participate in triggering intestinal and extra-intestinal diseases in sensitive people². The microbial diversity of the colon microbiota may result from a co-evolution with its host³ but may be affected by environmental conditions⁴. In order to identify all bacteria (including uncultured and fastidious) present in the colon, we use the culturomics strategy based on diversified culture conditions (temperature, media, and atmosphere), and rapid bacterial identification using matrix-assisted desorption ionization-time of flight mass spectrometry (MALDI-TOF MS)⁵⁻⁷. For putative new taxa, we use the taxonogenomic method that combines phenotypic characteristics and whole genome sequencing analysis to describe new bacterial species⁸⁻¹⁰. In 2016, we isolated the new bacterial strain Marseille-P2911 (=CSUR P2911 = DSM 103304), from a left colon sample of a 60-year-old patient who underwent a colonoscopy for the etiological investigation of iron-deficiency anemia in Marseille, France¹¹. This bacterium was identified by matrix-assisted desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) using a Microflex spectrometer¹² (Bruker Daltonics, Bremen, Germany). The strain was predicted to be affiliated with members of the family of Veillonellaceae but distinct form species with a validly published name. In the present study, we aimed at comparing strain Marseille-P2911 to its closely related phylogenetic neighbors, and at proposing the creation of the new genus *Colibacter massiliensis* gen. nov., sp. nov.

<footnote n="1">¹Aix Marseille Univ, Institut de Recherche pour le Développement (IRD), Service de Santé des Armées, AP-HM, UMR Vecteurs Infections Tropicales et Méditeranéennes (VITROME), Marseille, France. ²Institut Hospitalo-Universitaire Méditerranée-Infection, Marseille, France. ³Aix-Marseille Université, Institut de Recherche pour le Développement (IRD), UMR Microbes Evolution Phylogeny and Infections (MEPHI), Marseille, France. ⁴Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia. *email: pierre-edouard.fournier@univ-amu.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Phylogenetic tree highlighting the position of *Colibacter massiliensis* (in red) relative to other closely related bacterial taxa. Strains in bold are the eight closely related bacterial type strains with available genomes. The sequence of *Escherichia coli* ATCC 11775 served as outgroup. Genbank accession numbers for each 16s rRNA gene sequences are presented in parentheses. Sequences were here aligned by Muscle v3.8.31 with default parameters and finally 1247 nucleotides were considered to infer the phylogenetic relationship, using the Maximum Likelihood method with 1000 bootstrap replicates, within MEGA 7 software version 7.0.</figcaption></figure>

## Results

**Strain identification and classification.** Strain Marseille-P2911 was isolated from the left colon liquid sample of a 60-year-old man who underwent a colonoscopy for the etiological investigation of an iron-deficiency anemia. The patient provided signed informed consent, and the study was approved by the ethics committee of the Institut Fédératif de Recherche IFR48 under number 2016-010. Strain Marseille-P2911 could not be identified by our systematic MALDI-TOF MS screening as the score was lower than 1.8, suggesting that the corresponding species was not in the database (Fig. S1). Moreover, strain Marseille-P2911 exhibited a 94.59% 16S rRNA sequence similarity with *Anaeroglobus geminatus* strain AIP 313.00 (GenBank accession no. AF338413), the phylogenetically-closest species with standing in nomenclature (Fig. 1). As this value is lower than the 95% threshold proposed by Stackebrandt and Ebers for defining a new genus, strain Marseille-P2911 was considered as a representative of a putatively new genus within the family Veillonellaceae in the phylum Firmicutes.

**Phenotypic characteristics.** Growth was observed on 5% sheep blood-enriched Columbia agar (BioMérieux) at 37 °C after 2 days of incubation. Colonies from strain Marseille-P2911 showed neither pigmentation nor haemolysis. They were circular with a diameter of 0.1 mm, and transparent. Bacterial cells were Gram-positive, non-motile diplococci with a diameter of 0.4 to 0.6 µm, as determined by transmission electron microscopy (Fig. 2). Strain Marseille-P2911 grew only in anaerobic conditions. The sporulation test (20 minutes at 80°C) was negative. In addition, this bacterium had neither oxidase nor catalase activities. Biochemical
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Transmission electron microscopy of *Colibacter massiliensis* using a Technai G2 electron microscope operating at 60 Kv (the scale bar represents 500 nm).</figcaption></figure>

properties from strain Marseille-P2911, as determined by the use of API ZYM, API 20E and API 50 CH strips (bioMérieux), are detailed below in the description of *Colibacter massiliensis* gen. nov., sp. nov. paragraph.

According to the French Microbiology Society¹³, susceptibility tests showed that strain Marseille-P2911 was resistant to penicillin, piperacillin/tazobactam, ticarcillin-clavulanic acid, cefoxitin, ceftriaxone, ciprofloxacin, fosfomycin, trimethoprim-sulfamethoxazole, metronidazole, rifampin, amikacin and teicoplanin.

By comparison with closely related taxa, strain Marseille-P2911 differed in a combination of negative oxidase activity and positive sucrose/raffinose metabolism (Table 1). The most abundant fatty acids were saturated ones (83%), including dodecanoic acid that represented 25% of all fatty acids extracted from the bacterial cell wall. However, a specific 14:0 3-OH and several branched fatty acids were also described (Table S1).

**Genome sequencing information and genome properties.** The genome size of strain Marseille-P2911 was 1,715,864 bp long with a 50.2% G + C content. It was assembled into 2 scaffolds. Of the 1,655 predicted genes, 1,567 were protein-coding genes and 62 were RNAs (one complete rRNA operon, three additional 5S rRNAs and 49 tRNA genes). A total of 1,350 genes (81.57%) were assigned a putative function (by COGs) and 305 genes (18.43%) were annotated as hypothetical proteins. The genome properties and distribution of genes into COGs functional categories are detailed in Table S2 and Fig. 3. Genes putatively gained by hypothetical lateral gene transfer (LGT) were classified according to the bacterial families of origin (Fig. 4). Although we cannot rule out the possibility that some of the transfers may be from as yet unidentified taxa, most hypothetical lateral gene transfer (LGT)-acquired genes were obtained from members of the Veillonellaceae (90.6%) and Selenomonadaceae (1.6%) families.

**Comparison with closely related bacterial strains.** The genome of strain Marseille-P2911 was compared to the available genomes of eight closely related bacterial type strains. This comparison revealed that the genome size of our strain (1.72 Mb) was larger than that of *Dialister micraerophilus* (1.28 Mb), but smaller than those of *Anaeroglobus geminatus*, *Megasphaera cerevisiae*, *M. elsdenii*, *M. hexanoica*, *M. massiliensis*, *M. micronuciformis* and *M. paucivorans* (1.8, 3.24, 2.5, 2.88, 2.74, 1.77 and 2.91 Mb, respectively) (Table S3). The G + C content of strain Marseille-P2911 (50.2 mol %) was equal to that of *M. massiliensis*, but greater than those of all compared species (Table S3) except *Megasphaera elsdenii* (52.8%). The gene content of strain Marseille-P2911 (1,655) was similar to that of *D. micraerophilus* but smaller than those of other compared genomes (Table S3). The distribution of genes into COG categories was similar in all nine compared genomes (Fig. 5). Strain Marseille-P2911 shared 1169, 1161, 1129, 1109, 1083, 1081, 1071 and 666 orthologous genes with *M. micronuciformis*, *M. cerevisiae*, *M. hexanoica*, *M. massiliensis*, *A. geminatus*, *M. elsdenii*, *M. paucivorans* and *D. micraerophilus*, respectively (Table 2). Moreover, the MAGI analysis showed that AGIOS (Average Genomic Identity of Orthologous gene Sequences) values ranged from 47.76% between *M. elsdenii* and *D. micraerophilus*, to 50.93% between *M. paucivorans* and *D. micraerophilus*, among studied bacterial strains. Regarding strain Marseille-P2911, AGIOS values ranged from 48.86% with *D. micraerophilus* to 50.78% with *A. geminatus* (Table 2). Using dDDH analysis, strain Marseille-P2911 exhibited values ranging from 18.1% with *M. massiliensis* to 27.1% with *D. micraerophilus* (Table S4). These values are lower than the 70% threshold used for delineating prokaryotic species, thus confirming that this strain represents a new species. Finally, strain Marseille-P2911 exhibited average nucleotide identity (ANI) values ranging from 63.37% with *D. micraerophilus* to 74.42% with
<table><caption>**Table 1.** Compared phenotypic characteristics of studied species. Differential characteristics of *Colibacter massiliensis*¹¹ (CM), *Anaeroglobus geminatus*³⁷ (AG), *Dialister micraerophilus*³⁸ (DM), *Megasphaera paucivorans*³⁹ (MP), *Megasphaera micronuciformis*⁴⁰ (MMi), *Megasphaera cerevisiae*⁴¹,⁴² (MC), *Megasphaera elsdenii*⁴²,⁴³ (ME), *Megasphaera massiliensis*⁴⁴ (MMa) and *Megasphaera hexanoica*⁴⁵ (MH). +positive reaction, –negative reaction, ±variable, w week reaction, na data not available.</caption><thead><tr><th>Characteristics</th><th>CM</th><th>AG</th><th>DM</th><th>MP</th><th>MMi</th><th>MC</th><th>ME</th><th>MMa</th><th>MH</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><td>-</td></tr><tr><td>Production of Catalase</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Oxidase</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>Nitrate reductase</td><td>na</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>na</td><td>-</td></tr><tr><td>Gelatin hydrolysis</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>-</td></tr><tr><td>Utilisation of L-Arabinose</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>±</td><td>-</td><td>+</td><td>-</td></tr><tr><td>D-Galactose</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>w</td></tr><tr><td>D-Fructose</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>D-Mannose</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>±</td></tr><tr><td>D-Rhamnose</td><td>na</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>D-Mannitol</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>-</td></tr><tr><td>Salicin</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>D-Cellobiose</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>D-Maltose</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>+</td><td>-</td></tr><tr><td>D-Lactose</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>D-Trehalose</td><td>+</td><td>-</td><td>-</td><td>na</td><td>na</td><td>na</td><td>na</td><td>+</td><td>-</td></tr><tr><td>Inulin</td><td>-</td><td>+</td><td>na</td><td>-</td><td>+</td><td>-</td><td>-</td><td>na</td><td>-</td></tr><tr><td>Starch</td><td>-</td><td>-</td><td>-</td><td>na</td><td>na</td><td>na</td><td>na</td><td>na</td><td>-</td></tr><tr><td>Glycogen</td><td>-</td><td>+</td><td>na</td><td>-</td><td>+</td><td>-</td><td>-</td><td>na</td><td>-</td></tr><tr><td>D-Glucose</td><td>+</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>-</td><td>-</td></tr><tr><td>D-Xylose</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>Sucrose</td><td>+</td><td>-</td><td>-</td><td>na</td><td>-</td><td>na</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

A. geminatus. An ANI value lower than 80.5% suggesting that two strains belong to distinct genera, we considered that strain Marseille-P2911 was representative of a new genus (Fig. 6). Consequently, based on the presented phenotypic and genomic data, we propose the creation of the new genus *Colibacter* gen. nov., with strain Marseille-P2911T being the type strain of the novel species *Colibacter massiliensis* gen. nov., sp. nov.

Discussion

The gastrointestinal tract harbors a complex microbial microflora whose dynamic composition is important for health¹⁴. Here we aimed at describing a new bacterial species to enrich the knowledge on the human microbiome, using the culturomics and taxonogenomic strategies⁷,¹⁰,¹⁵.

The phylogenetic and phenotypic analysis of the new strain Marseille-P2911 revealed several distinct traits when compared to other members of the family Veillonellaceae¹⁶, suggesting that it could be classified in a new species of a new genus. The family Veillonellaceae is currently made of six genera of Gram-negative bacteria, including *Veillonella* (twelve species), *Megasphaera* and *Dialister* (five species each), *Allisonella*, *Anaeroglobus* and *Negativicoccus* (one species each)¹⁷.

As strain Marseille-P2911, many members of the family Veillonellaceae were detected in humans. For example, *Allisonella histaminiformans* was previously isolated from the vagina and A. geminatus was isolated from the gastrointestinal tract, stool and skin¹⁷. Moreover, these two species were demonstrated to be pathogenic in several diseases such as community-acquired pneumonia and advanced caries¹⁷.

The genomic content of strain Marseille-P2911 (dDDH, orthoANI, and AGIOS values) comforted its new species status. We observed a significant similarity to genes from the family Veillonellaceae (90.6%) (Fig. 4). A small rate of hypothetical lateral gene transfer (9.4%) from other bacterial families was observed, notably several toxin/antitoxin system-related genes putatively acquired from the families Eubacteriaceae, Lachnospiraceae, Lactobacillaceae, Selenomonadaceae and Streptococcaceae). So, we formally propose the creation of the new genus and species *Colibacter massiliensis* gen. nov., sp. nov., within the family Veillonellaceae. The type strain, Marseille-P2911T, was deposited in the DSMZ and CSUR collections under accession numbers DSM 103304 and CSUR P2911, respectively. The 16S rRNA and genome sequences are available in GenBank under accession numbers LT576403 and FMIY00000000, respectively. The Digital Protologue TaxoNumbers (http://imedea.uibcsic.es/dprotologue/index.php) of C. massiliensis gen. nov., sp. nov. is GA00103.

**Description of Colibacter gen. nov.** *Colibacter* (Co.li.bac'ter, N.L. masc. n. *colibacter*, composed of colon, the organ from which the strain was first isolated, and bacter, a rod. *Colibacter*, a bacterium from the colon). Cells are diplococci, anaerobic, Gram-positive, nonmotile, asporogenous, catalase and oxidase negative.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Graphical circular map of the genome. From outside to the center: genes on the forward strand colored by Clusters of Orthologous Groups of proteins (COG) categories (only genes assigned to COG). Genes on the reverse strand colored by COG categories (only genes assigned to COG). RNA genes (tRNAs green, rRNA red). G + C content and G + C skew.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Network showing the origin of genes in *C. massiliensis* according to the bacterial families.</figcaption></figure>

**Description of Colibacter massiliensis gen. nov., sp. nov.** *Colibacter massiliensis* (mas.si.li.ensis; L. fem. adj. *massiliensis*, for Massilia, the Roman name of Marseille, where the type strain was first isolated).

In addition to the description features of the genus, cells have a diameter varying from 0.4 to 0.6 µm. Colonies grown on 5% sheep blood-enriched Columbia agar (bioMérieux) are circular and transparent after 2 days of incubation in anaerobic atmosphere. Growth occurs at 37 °C. Cells grow anaerobically only. Using an API ZYM strip, a positive reaction was observed for alkaline phosphatase, esterase (C4), lipase esterase (C8), leucine arylamidase, valine arylamidase, cystine arylamidase, α-chymotrypsine, phosphatase acid, naphtol-AS-BI-phosphohydrolase, β-galactosidase, α-glucosidase, β-glucosidase and N-acetyl-β-glucosaminidase, but negative reactions were
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Distribution of functional classes of predicted genes according to the clusters of orthologous groups of proteins of *Colibacter massiliensis* among other bacterial taxa type strain.</figcaption></figure>

<table><caption>**Table 2.** Numbers of orthologous proteins shared between genomes (upper right) and AGIOS values (lower left). The numbers of proteins per genome are indicated in bold. *Colibacter massiliensis* (CM), *Anaeroglobus geminatus* (AG), *Dialister micraerophilus* (DM), *Megasphaera paucivorans* (MP), *Megasphaera micronuciformis* (MMi), *Megasphaera cerevisiae* (MC), *Megasphaera elsdenii* (ME), *Megasphaera massiliensis* (MMa) and *Megasphaera hexanoica* (MH).</caption><thead><tr><th></th><th>MP</th><th>MMi</th><th>MMa</th><th>MH</th><th>AG</th><th>CM</th><th>ME</th><th>DM</th><th>MC</th></tr></thead><tbody><tr><td>MP</td><td><strong>2717</strong></td><td>1089</td><td>1309</td><td>1340</td><td>972</td><td>1071</td><td>1300</td><td>672</td><td>1418</td></tr><tr><td>MMi</td><td>50.45</td><td><strong>1745</strong></td><td>1154</td><td>1162</td><td>1090</td><td>1169</td><td>1124</td><td>665</td><td>1201</td></tr><tr><td>MMa</td><td>49.55</td><td>50.37</td><td><strong>2536</strong></td><td>1405</td><td>1029</td><td>1109</td><td>1441</td><td>675</td><td>1444</td></tr><tr><td>MH</td><td>50.26</td><td>50.18</td><td>50.6</td><td><strong>2785</strong></td><td>1019</td><td>1129</td><td>1367</td><td>692</td><td>1442</td></tr><tr><td>AG</td><td>49.85</td><td>50.59</td><td>50.39</td><td>50.58</td><td><strong>1969</strong></td><td>1083</td><td>1003</td><td>594</td><td>1075</td></tr><tr><td>CM</td><td>49.81</td><td>50.42</td><td>50.46</td><td>50.31</td><td>50.78</td><td><strong>1661</strong></td><td>1081</td><td>666</td><td>1161</td></tr><tr><td>ME</td><td>48.97</td><td>49.83</td><td>50.78</td><td>50.42</td><td>50.4</td><td>50.38</td><td><strong>2333</strong></td><td>683</td><td>1434</td></tr><tr><td>DM</td><td>50.93</td><td>50.1</td><td>48.48</td><td>49.22</td><td>49.07</td><td>48.86</td><td>47.76</td><td><strong>1245</strong></td><td>683</td></tr><tr><td>MC</td><td>50.61</td><td>50.47</td><td>50.08</td><td>50.3</td><td>50.41</td><td>50.08</td><td>49.79</td><td>50.18</td><td><strong>3086</strong></td></tr></tbody></table>

obtained for lipase, trypsin, α-galactosidase, β-glucuronidase, α-mannosidase and α-fucosidase. Using an API 20E strip, a positive reaction were obtained for β-galactosidase, arginine dihydrolase, urea and gelatin hydrol-ysis, glucose fermentation, mannitol, inositol, sorbitol, rhamnose, sucrose, melibiose, amygdalin and arabinose. Using an API 50 CH strip, a positive result was shown for D-ribose, D-galactose, D-glucose, D-fructose, D-mannose, D-mannitol, D-sorbitol, N-acetylglucosamine, amygdalin, arbutin, esculin iron citrate, salicin, D-celibiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-trehalose, D-melezitose, D-raffinose and D-tagatose. In contrast, negative reactions were observed for glycerol, methyl-α-D-glucopyranoside, D-turanose, potassium 5-Ketogluconate, erythritol, D-arabinose, L-arabinose, D-xylose, L-xylose, D-adonitol, methyl-β-D-xylopyranoside, L-sorbose, L-rhamnose, dulcitol, inositol, methyl-α-D-mannopyranoside, inulin, starch, glycogen, xylitol, gentiobiose, D-lyxose, D-fucose, L-fucose, D-arabinol, L-arabinol, potassium gluconate and potassium 2-ketogluconate. The major fatty acid found for this strain was dodecanoic acid. The genome is 1,715,864 bp long with 50.2% G + C content.

The type strain, Marseille-P2911<sup>T</sup>, isolated from the left colon of a patient, was deposited in the DSMZ and CSUR collections under accession numbers DSM 103304 and CSUR P2911, respectively. The 16S rRNA and genome sequences are available in GenBank under accession numbers LT576403 and FMIY00000000, respectively.

## Materials and Methods

**Strain isolation and phenotypic tests.** The left colon liquid sample of a 60-year-old-man, who underwent a colonoscopy for an aetiological investigation of iron-deficiency anemia, was initially collected in La Timone Hospital in Marseille, France. From this sample, strain Marseille-P2911 was isolated after 3 days of pre-incubation in an anaerobic blood culture bottle (VersaTREK REDOX 2, Thermo Scientific, Villebon sur Yvette, France) supplemented with 5 mL of 0.2 µm-filtered rumen. This enriched liquid medium was then inoculated on 5% sheep blood-enriched Columbia agar (BioMérieux, Marcy l'Etoile, France) followed by an incubation at 37°C in anaerobic atmosphere (AnaeroGEN Compact, Oxoid, Thermo Scientific, Dardilly, France). MALDI-TOF mass spectrometry (MS) protein analysis was carried out using a Microflex spectrometer<sup>18</sup> (Bruker Daltonics,
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Heatmap generated with OrthoANI values calculated using the OAT software between *Colibacter massiliensis* and other closely related taxa with standing in nomenclature.</figcaption></figure>

Bremen, Germany). Strain Marseille-P2911 spectra were imported into the MALDI BioTyper software (version 2.0, Bruker) and analysed by standard pattern matching (with default parameter settings). Interpretation of the scores was performed as previously described¹⁹.

Moreover, the 16S rRNA gene was sequenced using the fD1-rP2 primer pairs as previously described²⁰, using a 3130-XL sequencer (Applied Biosciences, Saint Aubin, France). A phylogenetic tree was obtained using the Maximum Likelihood method and Kimura 2-parameter within the MEGA 7 software²¹. Several growth temperatures (20, 28, 37, 45 and 55°C) on 5% sheep blood-enriched Columbia agar medium (BioMérieux, Marcy l'Etoile, France) were tested. Growth of strain Marseille-P2911 was tested under different atmospheres (anaerobic, aerobic and microaerophilic (CampyGEN, Oxoid). API ZYM and API 50CH strips (BioMérieux) were used to evaluate the biochemical properties of the strain test according to the manufacturer's instructions. All API experiments were performed under anaerobic conditions. Using API 50CH, API 20E and API ZYM strips, strain Marseille-P2911 was incubated for 48, 24 and 4 hours, respectively. The standard disc method was applied for antimicrobial susceptibility testing according to the French Microbiology Society¹³. Finally, cellular fatty acid methyl ester (FAME) analysis was performed by GC/MS. Two samples were prepared with approximately 18 mg of bacterial biomass per tube harvested from several culture plates. Briefly, fatty acid methyl esters were separated using an Elite 5-MS column and monitored by mass spectrometry (Clarus 500 - SQ 8S, Perkin Elmer, Courtaboeuf, France)²². GC/MS analyses were carried out as described before²³. Spectral database search was performed using MS Search 2.0 operated with the Standard Reference Database 1 A (NIST, Gaithersburg, USA) and the FAMEs mass spectral database (Wiley, Chichester, UK). For transmission electronic microscopy, detection formvar-coated grids were dropped onto a 40 µL bacterial suspension before incubation at 37°C for 30 minutes. Then, the grids were incubated on 1% ammonium molybdate for 10 seconds, dried on blotting paper and finally observed using a Tecnai G20 transmission electron microscope (FEI, Limeil-Brevannes, France) at an operating voltage of 60 Kv. All methods were performed in accordance with the relevant guidelines and regulations.

**Extraction and genome sequencing.** After a pretreatment by lysozyme incubation at 37°C for 2 hours, DNA of strain Marseille-P2911 was extracted using an EZ1 biorobot (Qiagen) with the EZ1 DNA Tissue kit. The elution volume was 50 µL. Genomic DNA (gDNA) was quantified by a Qubit assay with the high sensitivity kit (Life technologies, Carlsbad, CA, USA) at 92.6 ng/µL. Genomic DNA was sequenced on a MiSeq sequencer (Illumina Inc, San Diego, CA, USA) with the Mate Pair strategy. The gDNA was barcoded in order to be mixed with 11 others projects with the Nextera Mate Pair sample prep kit (Illumina). The Mate Pair library was prepared with 1.5 µg of genomic DNA using the Nextera Mate Pair Illumina guide. The gDNA sample was simultaneously fragmented and tagged with a Mate Pair junction adapter. The pattern of the fragmentation was validated on an Agilent 2100 BioAnalyzer (Agilent Technologies Inc, Santa Clara, CA, USA) with a DNA 7500 labchip. DNA
fragments ranged in size from 1.5 kb up to 11 kb with an optimal size at 8.4 kb. No size selection was performed and 600 ng of tagmented fragments were circularized. The circularized DNA was mechanically sheared to small fragments with an optimal at 706 bp on the Covaris device S2 in microtubes (Covaris, Woburn, MA, USA). The library profile was visualized on a High Sensitivity Bioanalyzer LabChip (Agilent Technologies Inc, Santa Clara, CA, USA) and the final concentration library was measured at 13.256 nmol/l. The libraries were normalized at 2 nM after a denaturation step and dilution at 15 pM, loaded onto the reagent cartridge and then onto the instrument along with the flow cell. Automated cluster generation and sequencing run were performed in a single 39-hours run in a 2 × 251-bp format. Total information of 8.3 gb was obtained from a 910 K/mm² cluster density with a cluster passing quality control filters of 92.8% (16,316,000 passing filter paired reads). Within this run, the index representation for strain Marseille-P2911 was determined to 8.4%. The 688,244 paired reads were quality-checked using FastQC, trimmed using Trimmomatic version 0.36.6²⁴ and assembled in two scaffolds using the SPAdes version 3.5.0 software²⁵. The option “careful” was used in order to reduce the number of mismatches and short indels. Default parameters were applied for k values, i.e., k-mer values of 127, 99, 77, 55, 33, and 21. SSPACE²⁶ and GapFiller²⁷ were used to combine contigs, using default parameters²⁸,²⁹.

**Genome annotation and genome comparison.** The genome was annotated as previously described¹⁹. In addition, we used the Genome-to Genome Distance Calculator (GGDC) web server available at (http://ggdc.dsmz.de) to estimate the overall similarity among the compared genomes and to replace the wet-lab DNA-DNA hybridization (DDH) by a digital DDH (dDDH)³⁰,³¹. Average nucleotide identity analysis was also estimated using the orthoANI³² and MAGI³³ softwares. Antibiotic resistance genes (ARG) were searched using the ARG-ANNOT database and Bio-Edit interface³⁴. Assembled sequences were searched against the ARG-ANNOT database under moderately stringent conditions (e-value of 10⁻⁵) for the in silico ARG prediction. These putative ARGs were further confirmed through a BLAST search against non-redundant (nr) database in GenBank.

The presence of pathogenesis-related proteins was investigated using PathogeneFinder 1.1³⁵. Finally, predicted protein sequences of strain Marseille-P2911 were used as queries to search the NCBI GenBank non-redundant protein sequence database. These results were formatted to generate network of protein sequences using the Cytoscape tool³⁶.

Received: 6 February 2019; Accepted: 6 November 2019;

Published online: 20 November 2019

**References**

1. Berg, R. The indigenous gastrointestinal microflora. *Trends in Microbiology* **4**, 430–435 (1996).

2. Shanahan, F. The colonic microbiota in health and disease: *Current Opinion in Gastroenterology* **29**, 49–54 (2013).

3. Ley, R. E., Peterson, D. A. & Gordon, J. I. Ecological and Evolutionary Forces Shaping Microbial Diversity in the Human Intestine. *Cell* **124**, 837–848 (2006).

4. Chow, J., Lee, S. M., Shen, Y., Khosravi, A. & Mazmanian, S. K. Host-Bacterial Symbiosis in Health and Disease. In *Advances in Immunology* **107**, 243–274 (Elsevier, 2010).

5. Lagier, J.-C. et al. Microbial culturomics: paradigm shift in the human gut microbiome study. *Clinical Microbiology and Infection* **18**, 1185–1193 (2012).

6. Lagier, J.-C. et al. The Rebirth of Culture in Microbiology through the Example of Culturomics To Study Human Gut Microbiota. *Clinical Microbiology Reviews* **28**, 237–264 (2015).

7. Lagier, J.-C. et al. Culture of previously uncultured members of the human gut microbiota by culturomics. *Nature Microbiology* **1** (2016).

8. Ramasamy, D. et al. A polyphasic strategy incorporating genomic data for the taxonomic description of novel bacterial species. *International Journal of Systematic and Evolutionary Microbiology* **64**, 384–391 (2014).

9. Fournier, P.-E. & Drancourt, M. New Microbes New Infections promotes modern prokaryotic taxonomy: a new section "TaxonoGenomics: new genomes of microorganisms in humans". *New Microbes and New Infections* **7**, 48–49 (2015).

10. Abdallah, R. A. et al. The impact of culturomics on taxonomy in clinical microbiology. *Antonie van Leeuwenhoek* **110**, 1327–1337 (2017).

11. Mailhe, M. et al. 'Colibacter massiliensis' gen. nov., sp. nov., a new bacterial species isolated from human left colon. *New Microbes and New Infections* **17**, 23–24 (2017).

12. Seng, P. et al. Ongoing Revolution in Bacteriology: Routine Identification of Bacteria by Matrix-Assisted Laser Desorption Ionization Time-of-Flight Mass Spectrometry. *Clinical Infectious Diseases* **49**, 543–551 (2009).

13. Leclercq, R. et al. EUCAST expert rules in antimicrobial susceptibility testing. *Clinical Microbiology and Infection* **19**, 141–160 (2013).

14. Cremer, J., Arnoldini, M. & Hwa, T. Effect of water flow and chemical environment on microbiota growth and composition in the human colon. *Proceedings of the National Academy of Sciences* **114**, 6438–6443 (2017).

15. Fournier, P.-E., Lagier, J.-C., Dubourg, G. & Raault, D. From culturomics to taxonomogenomics: A need to change the taxonomy of prokaryotes in clinical microbiology. *Anaerobe* **36**, 73–78 (2015).

16. Rogosa, M. Transfer of Peptostreptococcus elsdenii Gutierrez et al. to a New Genus, Megasphaera CM. elsdenii (Gutierrez et al.) comb. nov. **3** (2018).

17. Marchandin, H. & Jumas-Bilak, Estelle. The Family Veillonellaceae. The Prokaryotes: Firmicutes and Tenericutes. 433–453 (2013).

18. Seng, P. et al. Identification of Rare Pathogenic Bacteria in a Clinical Microbiology Laboratory: Impact of Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry. *Journal of Clinical Microbiology* **51**, 2182–2194 (2013).

19. Hadjadj, L. et al. Non contiguous-finished genome sequence and description of Microbacterium gorillae sp. nov. *Standards in Genomic Sciences* **11** (2016).

20. Drancourt, M. et al. 16S Ribosomal DNA Sequence Analysis of a Large Collection of Environmental and Clinical Unidentifiable Bacterial Isolates. *J. CLIN. MICROBIOL.* **38**, 8 (2000).

21. Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Molecular Biology and Evolution* **33**, 1870–1874 (2016).

22. Sasser, M. Bacterial Identification by Gas Chromatographic Analysis of Fatty Acids Methyl Esters (GC-FAME). **6** (2006).

23. Dione, N. et al. Genome sequence and description of Anaerosalibacter massiliensis sp. nov. *New Microbes and New Infections* **10**, 66–76 (2016).

24. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* **30**, 2114–2120 (2014).
25. Bankevich, A. et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. *Journal of Computational Biology* **19**, 455–477 (2012).

26. Boetzer, M. et al. Scaffolding pre-assembled contigs using SSPACE. *Bioinformatics* **27**, 578–579 (2011).

27. Boetzer, M. & Pirovano, W. Toward almost closed genomes with GapFiller. *Genome B* (2012).

28. Anani, H. et al. Draft genome and description of Merdibacter massiliensis gen.nov., sp. nov., a new bacterium genus isolated from the human ileum. *Scientific Reports* **9** (2019).

29. Anani, H., Khodor, M., Raoult, D. & Fournier, P.-E. Whole-Genome Sequence of French Clinical Olivibacter jilunii Strain P8502. *Microbiology Resource Announcements* **8** (2019).

30. Auch, A. F. Digital DNA-DNA hybridization. *Standards in Genomic Sciences* **18** (2010).

31. 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* **14**, 60 (2013).

32. Lee, I., Ouk Kim, Y., Park, S.-C. & Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. *International Journal of Systematic and Evolutionary Microbiology* **66**, 1100–1103 (2016).

33. Weber, T. et al. antiSMASH 3.0—a comprehensive resource for the genome mining of biosynthetic gene clusters. *Nucleic Acids Research* **43**, W237–W243 (2015).

34. Gupta, S. K. et al. ARG-annot, a New Bioinformatic Tool To Discover Antibiotic Resistance Genes in Bacterial Genomes. *Antimicrobial Agents and Chemotherapy* **58**, 212–220 (2014).

35. Cosentino, S., Voldby Larsen, M., Møller Aarestrup, F. & Lund, O. PathogenFinder - Distinguishing Friend from Foe Using Bacterial Whole Genome Sequence Data. *PLoS ONE* **8**, e77302 (2013).

36. Smoot, M. E., Ono, K., Ruscheinski, J., Wang, P.-L. & Ideker, T. Cytoscape 2.8: new features for data integration and network visualization. *Bioinformatics* **27**, 431–432 (2011).

37. Carlier, J. P. Anaeroglobus geminatus gen. nov., sp. nov., a novel member of the family Veillonellaceae. *International Journal of Systematic and Evolutionary Microbiology* **52**, 983–986 (2002).

38. Jumas-Bilak, E. Dialister micraerophilus sp. nov. and Dialister propionicaciens sp. nov., isolated from human clinical samples. *International Journal of Systematic and Evolutionary Microbiology* **55**, 2471–2478 (2005).

39. Juvonen, R. Megasphaera paucivorans sp. nov., Megasphaera sueciensis sp. nov. and Pectinatus haikarae sp. nov., isolated from brewery samples, and emended description of the genus Pectinatus. *International Journal of Systematic and Evolutionary Microbiology* **56**, 695–702 (2006).

40. Marchandin, H. Phylogenetic analysis of some Sporomusa sub-branch members isolated from human clinical specimens: description of Megasphaera micronuciformis sp. nov. *International Journal of Systematic and Evolutionary Microbiology* **53**, 547–553 (2003).

41. Engelmann, U. & Weiss, N. Megasphaera cerevisiae sp. nov.: A New Gram-negative Obligately Anaerobic Coccus Isolated from Spoiled Beer. *Systematic and Applied Microbiology* **6**, 287–290 (1985).

42. Shetty, S. A., Marathe, N. P., Lanjekar, V., Ranade, D. & Shouche, Y. S. Comparative Genome Analysis of Megasphaera sp. Reveals Niche Specialization and Its Potential Role in the Human Gut. *PLoS ONE* **8**, e79353 (2013).

43. Gutierrez, J., Davis, R. E., Lindahl, I. L. & Warwick, E. J. Bacterial changes in the rumen during the onset of feed-lot bloat of cattle and characteristics of *Peptostreptococcus elsdenii* n. sp. *Applied Microbiology* **7**, 16–22 (1959).

44. Padmanabhan, R. et al. Non-contiguous finished genome sequence and description of Megasphaera massiliensis sp. nov. *Standards in Genomic Sciences* **8**, 525–538 (2013).

45. Jeon, B. S., Kim, S. & Sang, B.-I. Megasphaera hexanoica sp. nov., a medium-chain carboxylic acid-producing bacterium isolated from a cow rumen. *International Journal of Systematic and Evolutionary Microbiology* **67**, 2114–2120 (2017).

## Acknowledgements

The study was supported by the Méditerranée Infection foundation, the National Research Agency under the program "Investissements d'avenir", reference ANR-10IAHU-03 and by Région Provence Alpes Côte d'Azur and European funding FEDER PRIMI.

## Author contributions

H.A. performed the genomic analysis and phenotypic tests, R.A. and M.K. participated in the genomic analysis, A.F. did the electronic microscopy, D.Ri. and M.M. isolated the bacterium, D.Ra., P.E.F. and H.A. wrote the paper and designed the study. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-53791-1.

Correspondence and requests for materials should be addressed to P.-E.F.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019