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

HAL
open science

# Leucobacter manosquensis sp. nov.-A Novel Bacterial Species Isolated from Healthy Human Skin

Manon Boxberger, Sibylle Magnien, Angeline Antezack, Clara Rolland, Marine Makoa Meng, Cheikh Ibrahima Lo, Bernard La Scola, Nadim Cassir

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

## To cite this version:

Manon Boxberger, Sibylle Magnien, Angeline Antezack, Clara Rolland, Marine Makoa Meng, et al. "Leucobacter manosquensis sp. nov.-A Novel Bacterial Species Isolated from Healthy Human Skin." Microorganisms, 2023, 11 (10), pp.2535. 10.3390/microorganisms11102535. [https://doi.org/10.3390/microorganisms11102535]

HAL Id: hal-04462954

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

Submitted on 16 Feb 2024

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_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_4.png" /></figure>

Article

# Leucobacter manosquensis sp. nov.—A Novel Bacterial Species Isolated from Healthy Human Skin

Manon Boxberger <fnref n="1" />,<fnref n="2" /> ID, Sibylle Magnien <fnref n="1" />,<fnref n="2" />, Angeline Antezack <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="4" /> ID, Clara Rolland <fnref n="1" />,<fnref n="2" /> ID,
Marine Makoa Meng <fnref n="1" />,<fnref n="2" />, Cheikh Ibrahima Lo <fnref n="1" /> ID, Bernard La Scola <fnref n="1" />,<fnref n="2" /> and Nadim Cassir <fnref n="1" />,<fnref n="2" />,* ID

**Abstract:** Extending our knowledge on human skin microbiota is a challenge to better decipher its role in health and disease. Using the culturomics method, we isolated strain Marseille-Q4368 from the healthy forehead of a 59-year-old woman. We describe here the main characteristics of this bacterium using a taxonogenomic approach. This new bacterial species is Gram-positive, non-motile, and non-spore-forming. Its 16S rRNA sequence exhibited a similarity of 99.59% with *Leucobacter chromiiresistens*, the most closely related species in terms of nomenclature. However, a digital DNA-DNA hybridization analysis between these two species revealed a maximum identity similarity of only 27.5%. We found phenotypical and genomic differences between strain Marseille-Q4368 and its closely related species. These findings underscore the classification of this bacterium as a distinct species. Hence, we propose the name *Leucobacter manosquensis* sp. nov. strain Marseille-Q4368 (=CSUR Q4368 = DSM 112403) for this newly identified bacterial species.

**Keywords:** bacteria; taxonogenomics; culturomics; human skin microbiota; new species; *Leucobacter*

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

check for
updates

**Citation:** Boxberger, M.; Magnien, S.; Antezack, A.; Rolland, C.; Makoa Meng, M.; Lo, C.I.; La Scola, B.; Cassir, N. *Leucobacter manosquensis* sp. nov.—A Novel Bacterial Species Isolated from Healthy Human Skin. *Microorganisms* **2023**, *11*, 2535.
https://doi.org/10.3390/microorganisms11102535

Academic Editor: Roberto Di Marco

Received: 5 September 2023

Revised: 22 September 2023

Accepted: 9 October 2023

Published: 11 October 2023

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

**Copyright:** © 2023 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/).

## 1. Introduction

Pharmaceutical and cosmetic innovation have long been based on the search for innovative compounds, often extracted from plants. A major turning point in the history of medicine, with a strong impact on human health, came with the discovery of antibiotics and their natural producers: microorganisms. Fungal and bacterial microorganisms represent promising sources of new compounds for prophylactic or therapeutic treatments. While the medical applications arising from our knowledge on digestive microbiota are now widely used and recognized, they are not yet as well developed concerning cutaneous microbiota and the potential beneficial impact of pre-probiotics and bacteriotherapy on the skin. Not-yet-cultured microorganisms are certainly an underestimated source of interesting compounds. Extending our knowledge on human skin microbiota is a challenge to better decipher its role in health and disease. Recent studies have shown that human skin harbors a wide diversity of microbes, far beyond commonly isolated bacteria. Many of the bacteria that are part of the skin microbiota are difficult to isolate because of difficulties in reproducing the human skin environment in vitro [1]. Knowledge about the skin microbiota has grown exponentially owing in part to the recent advent of metagenomics and 16S rRNA pyrosequencing as tools for microbiological identification. However, culture methods remain an essential tool for studying the characteristics of microorganisms in vivo.

Since 2012, our laboratory has set up an innovative approach aiming at surpassing the limitations associated with culture-independent techniques. These new strategies, called

<footnote n="1">¹ IHU Méditerranée Infection, 19-21 Boulevard Jean Moulin, 13005 Marseille, France; manon.boxberger@hotmail.fr (M.B.); sibylle.magnien@gmail.com (S.M.)</footnote>

<footnote n="2">² Institut de Recherche Pour le Développement (IRD), Assistance Publique-Hôpitaux de Marseille (AP-HM), MEPHI, Aix-Marseille Université, 19 Boulevard Jean Moulin, 13005 Marseille, France</footnote>

<footnote n="3">³ École de Médecine Dentaire, Faculté des Sciences Médicales et Paramédicales, Aix-Marseille Université, Boulevard Jean Moulin, 13385 Marseille, France</footnote>

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

<footnote>* Correspondence: cassirnadim@gmail.com</footnote>
culturomics and taxonogenomics [2,3], consist of combining several isolation conditions to improve culture sensitivity, and MALDI-TOF mass spectrometry [4] and 16SrRNA gene amplification and sequencing [5] for the identification of a large number of isolated bacteria. In particular, the supplementation of natural compounds and antibiotics in the culture medium allowed us to unhide the minority species or to mimic a natural environment through selection processes [6,7]. This technique has allowed the discovery of more than 700 new bacterial species from various environments since its development [8]. As part of the culturomics project aimed at describing the skin microbiota, we isolated the strain Marseille-Q4368 [1]. The *Leucobacter* genus belongs to the *Microbacteriaceae* family and was first proposed in 1996 by Takeuchi et al., with the first species named *Leucobacter komagatae* [9]. *Leucobacter* takes its name from the morphology of the cells, which are rods and light-colored. Currently, there are 41 known species, of which 33 have been validly published with a correct name (https://lpsn.dsmz.de/genus/leucobacter, accessed on 4 September 2023). Moreover, the species found within the genus have been isolated from different environments and *Leucobacter manosquensis* is the first species of this genus to have been isolated from healthy human skin.

Here, we describe *Leucobacter manosquensis*, a new bacterial species isolated from the healthy skin using a taxonogenomic polyphasic approach. It includes phylogenomic analyses, wall fatty acid composition, and phenotypic characterization.

**2. Materials and Methods**

**2.1. Sample Acquisition and Strain Isolation**

Samples were obtained by swabbing a 10 cm² area of skin from the forehead and hands of volunteer healthy women. The samples were collected using sterile swabs soaked in the Culture Top® transport medium (C-top Ae-Ana, Eurobio, Les Ulis, France). The study was validated by the ethics committee Sud-Est IV (reference ID-RCB: 2019-A01508-49). Each sample mixed with the transport media was serially diluted, and 50 µL of each dilution was incubated under aerobic conditions at 31 °C after being seeded (DS) in Columbia Agar (bioMérieux, Marcy l'Etoile, France). To identify the strains, a Matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry protein analysis was carried out in triplicate using a Microflex spectrometer (Bruker Daltonics, Bremen, Germany). Whole genome sequencing (WGS) was performed as described below. The strain spectra were imported into the MALDI BioTyper software (version 3.0, Bruker, Bremen, Germany) and analyzed using standard pattern matching with the default parameters. Our database (https://www.mediterranee-infection.com/acces-ressources/base-de-donnees/urms-data-base/, accessed on 4 September 2023) was then incremented with the spectra of this new bacterial species.

**2.2. Phenotypic Tests**

Different growth temperatures (20 °C; 31 °C; 37 °C; 45 °C; and 56 °C), atmosphere conditions (anaerobic, aerobic, and microaerophilic) using generator bags (CampyGEN, Oxoid, Columbia, MD, USA), and pH conditions (5; 6.5; 7.5; 8.5) were tested. The biochemical properties of these strains were tested using API ZYM, API 20 NE, and API 50 CH strips (bioMérieux, Marcy l'Etoile, France) in accordance with the manufacturer's instructions. To evaluate the bacterial structure, a colony was collected from agar and immersed into a 2.5% glutaraldehyde fixative solution. The suspension was vortexed, passed ten times through a 21-gauge needle to separate the bacterial colonies, and fixed on an uncoated glass slide using cytocentrifugation. A 1% ammonium molybdate negative stain was applied for 1 min before gently washing the slide with 0.2 µm filtered distilled water. The slide was air-dried and examined using scanning electron microscopy on a TM4000 microscope (Hitachi High-Tech, HHT, Tokyo, Japan) with a 15 kV voltage. A motility test was performed using the semi-solid TTC media.

Cellular fatty acid methyl esters (FAME) analysis was performed using GC/MS. Two samples of each strain were prepared with approximately 110 mg of bacterial biomass
per tube harvested from several culture plates. The fatty acid methyl esters were prepared as described by Sasser (2006) [10]. GC/MS analyses were carried out as described before [11]. Briefly, the fatty acid methyl esters were separated using an Elite 5-MS column and monitored using mass spectrometry (Clarus 500—SQ 8 S, Perkin Elmer, Courtaboeuf, France). A spectral database search was performed using MS Search 2.0 on the Standard Reference Database 1A (NIST, Gaithersburg, USA) and the FAME mass spectral database (Wiley, Chichester, UK). To evaluate the bacterial structure, a colony was collected from agar and immersed into a 2.5% glutaraldehyde fixative solution. The suspension was vortexed, passed ten times through a 21-gauge needle to separate the bacterial colonies, and fixed on an uncoated glass slide using cytocentrifugation. A 1% ammonium molybdate negative stain was applied for 1 min before gently washing the slide with 0.2 µm filtered distilled water. The slide was air-dried and examined using scanning electron microscopy on a TM4000 microscope (Hitachi High-Tech, HHT, Tokyo, Japan) with a 15 kV voltage. A motility test was performed using the semi-solid TTC media as described by Tittsler et al. [12]. The sporulation was evaluated by collecting a colony from agar in 1 mL of filtered phosphate buffer saline (PBS) and using a thermal shock for 20 min at 60 °C. The biochemical properties of the strain Marseille-Q4368 were tested using API ZYM, API 20 NE, and API 50 CH strips (bioMérieux, Marcy L'Etoile, France) in accordance with the manufacturer's instructions.

## 2.3. Genome Sequencing, Annotation, and Genome Comparison

The genomic DNA (gDNA) of the strain Marseille-Q4368 was extracted in two steps: a mechanical treatment was first performed with glass beads being acid-washed (G4649-500 g Sigma, St. Louis, MO, USA) using a FastPrep-24™ 5G grinder (mpBio, Santa Ana, CA, USA) at maximum speed (6.5) for 90 s. Then, after 30 min of lysozyme incubation at 37 °C, the DNA was extracted using the EZ1 biorobot (Qiagen, Hilden, Germany) with the EZ1 DNA tissue kit. The elution volume was of 50 µL. The gDNA was quantified using a Qubit assay with a high-sensitivity kit (Life Technologies, Carlsbad, CA, USA) to 0.2 ng/µL. The genomic DNA was then sequenced using the MiSeq technology (Illumina Inc., San Diego, CA, USA) with the paired end strategy prepared with the Nextera XT DNA sample prep kit (Illumina, San Diego, CA, USA). To prepare the paired end library, dilution was performed to require 1 ng of each genome as input. The "tagmentation" step fragmented and tagged the DNA. Then, limited-cycle PCR amplification (12 cycles) completed the tag adapters and introduced dual-index bar codes. After purification with AMPure XP beads (Beckman Coulter Inc., Fullerton, CA, USA), the libraries were then normalized to specific beads according to the Nextera XT protocol (Illumina). The 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. To improve the quality of the assemblies, an Oxford Nanopore approach was performed on 1D genomic DNA sequencing for the MinIon device using the SQK-LSK109 kit. The library was constructed from 1 µg of genomic DNA without fragmentation and end repair. Adapters were ligated to both ends of the genomic DNA. After purification with the AMPure XP beads (Beckman Coulter Inc., Fullerton, CA, USA), the library was quantified using a Qubit assay with the high-sensitivity kit (Life Technologies, Carlsbad, CA, USA). For bacterial identification, we used the "What's in my pot" (WIMP) workflow (Epi2me, Oxford Nanopore technologies), which uses the RefSeq sequence database from the NCBI (https://www.ncbi.nlm.nih.gov/refseq/, accessed on 4 September 2023).

The NCBI Prokaryotic Genome Annotation Pipeline was used for genome annotation [13]. The genome sequence data were uploaded to the Type (Strain) Genome Server (TYGS), a free bioinformatics platform available at https://tygs.dsmz.de/, for a whole-genome-based taxonomic analysis [14]. Antibiotic resistance genes and the presence of pathogenesis-related proteins were investigated using ABRicate v1.0.1 against ARG-ANNOT [15], EcOH [16], the NCBI Bacterial Antimicrobial Resistance Reference Gene
Database [17], PlasmidFinder [18], ResFinder [19], CARD [20], and VFDB [21] using the Online Galaxy platform [22].

Determination of the closest type strain genomes was carried out in two complementary ways. First, all the user genomes were compared against all the type strain genomes available in the TYGS database using the MASH algorithm [23], and the 10 type strains with the smallest MASH distances chosen per user genome. Second, an additional set of 10 closely related type strains was determined via the 16S rDNA gene sequences. These were extracted from the user genomes using RNAmmer [24]. BLAST was then performed with each sequence [25] against the 16S rDNA gene sequence of each of the 12,983 type strains currently available in the TYGS database. This was used as a proxy to find the best 50 matching type strains for each user genome and to calculate the precise distances using the Genome BLAST Distance Phylogeny approach (GBDP) under the algorithm "coverage" and distance formula d5 [26]. These distances were used to determine the 10 closest type strain genomes for each of the user genomes. All pairwise comparisons among the set of genomes were conducted using GBDP and accurate intergenomic distances were inferred using the algorithm "trimming" and distance formula d5. 100 distance replicates were calculated each. The digital DDH (dDDH) values and confidence intervals were calculated using the recommended settings of the GGDC2. In parallel, the degree of genomic similarity between the strains of interest and closely related species was estimated using the OrthoANI software 0.5.0 with the default parameters [27]. The closest species were determined on the basis of DDH. Trees were inferred using FastME 2.1.6.1 [28] from the GBDP distances calculated using the 16S rDNA gene sequences or a whole genome sequence. The branch lengths were scaled in terms of GBDP distance formula d5. The numbers above the branches are the GBDP pseudo-bootstrap support values >60% from 100 replications, with an average branch support of 84.3%. The tree was rooted at the midpoint and regenerated with the iTOL tool v5 [29].

### 3. Results

#### 3.1. *Identification and Phylogenetic Analysis of Strain Marseille-Q4368*

Following a 24 h incubation under aerobic conditions at 31 °C on Columbia agar supplemented with 5% sheep blood, the bacterial colonies from the Marseille-Q4368T strain were subjected to MALDI-TOF MS analysis for identification. The results obtained were inconclusive for identification at the strain level. Then, a phylogenetic analysis was carried out using the 16S rRNA gene sequence. The 16S rRNA sequence of strain Marseille-Q4368 showed 99.5% of similarity when compared to *Leucobacter chromiiresistens* JG 31, the most closely related species in terms of phylogenetic classification (Figure 1). Considering that the obtained percentage similarity of the 16S rRNA gene exceeds the recommended threshold value (98.65%) for defining the species barrier in prokaryotes, additional investigations encompassing morphological, phenotypic, and genomic analyses have been conducted.

#### 3.2. *Phenotypic Characteristics of Strain Marseille-Q4368*

Strain Marseille-Q4368T was facultative aerobic and grew at 31 °C in 24 h. The growth was achieved at a temperature between 20 °C and 37 °C. The optimal pH was between 6.5 and 8.5. NaCl concentrations ranging from 5 g·L⁻¹ to 15 g·L⁻¹ showed no distinction in growth. The colonies had a translucent gray pigmentation and did not hemolyze. Moreover, the cells were non-motile and did not form spores. This bacterium was Gram-positive (supplementary data—Figure S1). The bacteria were rod-shaped, measuring approximately 1.4 × 0.3 µm. (Figure 2).
<figure><img src="image_7.png" /></figure>

*Leucobacter zeae* DSM 103328
*Leucobacter chironomi* DSM 19883
*Leucobacter muris* DSM 101948
*Leucobacter celer* subsp. astrifaciens DSM 27159
*Leucobacter triazinivorans* JW-1
*Leucobacter massiliensis* 122RC15
*Leucobacter aridicollis* DSM 17380
*Leucobacter dentrificans* KACC 14055
*Leucobacter luti* RF6
*Leucobacter komagatae* DSM 8803
*Leucobacter japonicus* DSM 27158
*Leucobacter musarum* DSM 2710
**Leucobacter manosquensis strain Marseille- Q4368**
*Leucobacter chromiiresistens* JG 31
*Leucobacter chromiiresistans* DSM 22788

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

*Leucobacter dentrificans* KACC 14055
*Leucobacter luti* RF6
*Leucobacter aridicollis* DSM 17380
*Leucobacter komagatae* DSM 8803
*Leucobacter musarum* DSM 2710
*Leucobacter japonicus* DSM 27158
**Leucobacter manosquensis strain Marseille- Q4368**
*Leucobacter chromiiresistens* JG 31
*Leucobacter chromiiresistans* DSM 22788
*Leucobacter zeae* DSM 103328
*Leucobacter chironomi* DSM 19883
*Leucobacter muris* DSM 101948
*Leucobacter massiliensis* 122RC15
*Leucobacter triazinivorans* JW-1
*Leucobacter celer* subsp. astrifaciens DSM 27159

16S

(A)

WGS

(B)

**Figure 1.** (A) 16S rRNA-based phylogenetic tree (B) Whole-genome-based phylogenetic tree highlighting the position of *L. manosquensis* sp nov., strain Marseille-Q4368T relative to other closely related bacterial taxa.
<figure><img src="image_9.png" /><figcaption>**Figure 2.** Scanning electron microscopy of *L. manosquensis* sp. nov., strain Marseille-Q4368T using a TM4000 microscope (Hitachi High-Tech, HHT, Tokyo, Japan). The scale bar represents 10 µm.</figcaption></figure>

For this strain, Marseille-Q4368T, several API galleries were performed, including API ZYM, API 50 CH, and API 20 NE. Positive reactions were observed for the following enzymes: esculin ferric citrate, esterase (C4), esterase lipase (C8), leucine arylamidase, cystine arylamidase, valine arylamidase, naphthol-AS-BI-phosphohydrolase, D-mannitol, D-glucose, D-maltose, malic acid, and α-mannosidase. However, the remaining reactions yielded negative results. Additionally, the catalase test was positive, while the oxidase test was negative. The differential phenotypical characteristics of *Leucobacter manosquensis* strain Marseille-Q4368 and its closest related species are summarized in Table 1. Branched fatty acids were detected in the composition of strain Marseille-Q4368: 12-methyl-tetradecanoic acid (69%), 14-methyl-hexadecanoic acid (16%), and 14-methyl-pentadecanoic acid (7%). Additionally, unsaturated fatty acids were also detected in the cell walls (Table 2).

<table><caption>**Table 1.** Essential phenotypic criteria shared between 1, *Leucobacter manosquensis* strain Marseille-Q4368 and closely phylogenetical species such as 2, *Leucobacter chromiiresistens* JG 31 and 3, *Leucobacter musarum* CBX 152.</caption><thead><tr><th>Properties</th><th>1</th><th>2</th><th>3</th></tr></thead><tbody><tr><td>Cell diameter (µm)</td><td>1.4 × 0.3</td><td>0.3 × 0.7</td><td>NA</td></tr><tr><td>Oxygen requirement</td><td>FA</td><td>Aerobic</td><td>Aerobic</td></tr><tr><td>Gram stain</td><td>+</td><td>+</td><td>+</td></tr><tr><td>Motility</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Endospore formation</td><td>—</td><td>—</td><td>—</td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th>Properties</th><th>1</th><th>2</th><th>3</th></tr></thead><tbody><tr><td>Alkaline phosphatase</td><td>—</td><td>+/-</td><td>NA</td></tr><tr><td>Catalase</td><td>+</td><td>+</td><td>+</td></tr><tr><td>Oxidase</td><td>—</td><td>—</td><td>—</td></tr><tr><td>α-glucosidase</td><td>—</td><td>—</td><td>NA</td></tr><tr><td>β-galactosidase</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Acid from:</td><td></td><td></td><td></td></tr><tr><td>N-acetylglucosamine</td><td>—</td><td>—</td><td>—</td></tr><tr><td>L-arabinose</td><td>—</td><td>+</td><td>—</td></tr><tr><td>D-ribose</td><td>—</td><td>—</td><td>—</td></tr><tr><td>D-mannose</td><td>—</td><td>—</td><td>—</td></tr><tr><td>D-mannitol</td><td>+</td><td>—</td><td>+</td></tr><tr><td>D-glucose</td><td>+</td><td>—</td><td>—</td></tr><tr><td>D-fructose</td><td>—</td><td>—</td><td>—</td></tr><tr><td>D-maltose</td><td>+</td><td>—</td><td>—</td></tr><tr><td>D-lactose</td><td>—</td><td>—</td><td>—</td></tr><tr><td>G + C content (mol%)</td><td>67.34</td><td>NA</td><td>66.77</td></tr><tr><td>Isolation source</td><td>Human healthy skin</td><td>Uncontaminated soil</td><td>Nematode</td></tr></tbody></table>

FA = Facultative anaerobic; + = positive reaction; — = negative reaction; NA = data not available.

<table><caption>**Table 2.** Fatty acid component (%) analysis from *Leucobacter manosquensis* Marseille-Q4368 (Lman), *Leucobacter chromiireducens* JG 31 (Lchr), *Leucobacter aridicollis* L-9T (Lari), and *Leucobacter celer* NAL101 (Lcel).</caption><thead><tr><th>Fatty Acids</th><th>Name</th><th>Lman</th><th>Lchr</th><th>Lari</th><th>Lcel</th></tr></thead><tbody><tr><td>C<sub>15:0</sub> anteiso</td><td>12-methyltetradecanoic acid</td><td>60.7</td><td>55.9</td><td>49</td><td>46.5</td></tr><tr><td>C<sub>17:0</sub> anteiso</td><td>14-methyl-hexadecanoic acid</td><td>23.8</td><td>23.5</td><td>24</td><td>21.1</td></tr><tr><td>C<sub>16:0</sub> iso</td><td>14-methyl-pentadecanoic acid</td><td>13.1</td><td>12.2</td><td>12</td><td>16.0</td></tr><tr><td>C<sub>16:0</sub></td><td>Hexadecanoic acid</td><td>1.4</td><td></td><td>11</td><td></td></tr><tr><td>C<sub>15:0</sub> iso</td><td>13-methyl-tetradecanoic acid</td><td>TR</td><td>0.62</td><td></td><td>9.7</td></tr></tbody></table>

### 3.3. Comparative Genomic Analysis

The genome size of strain Marseille-Q4368 was 3,184,602 bp long with a 67.34% G + C content. The genome assembly of this strain was achieved on three contigs (with 9.44× coverage). Of the 2976 predicted genes, 2870 were protein-coding genes and 59 were RNAs (3 16S rRNAs, 3 5S rRNAs, 3, 23 rRNAs, 47 tRNAs, and 3 ncRNAs) (Figure 3).

The in silico resistome of the strain Marseille-Q4368T and the search for the virulence factors of this strain showed neither resistance nor virulence factors. The distribution of the functional classes of the predicted genes according to the clusters of orthologous groups (COGs) of proteins show that the genome of *Leucobacter manosquensis* showed a coherent structure compared to its closely related species (supplementary data—Figure S2).

Genomic data analysis revealed that the highest OrthoANI value was shared between *L. japonicus* and *L. musarum* at 95.3%. On the other hand, the lowest value was shared between *L. japonicus* and *L. aridicollis* at 75.1%. Consequently, strain Marseille-Q4368 had the highest OrthoANI value of 84.3% with *L. chromiiresistens*, which appears to be the closest phylogenetic species (Table 3). Furthermore, strain Marseille-Q4368 had a dDDH value of 27.5% with *L. chromiiresistens*. However, despite the dDDH value of 62.3% obtained between *L. japonicus* and *L. musarum*, the overall DDH values were below the recommended threshold of 70% to define a new prokaryotic species (Table 3) [31].
<figure><img src="image_10.png" /><figcaption>**Figure 3.** Graphical circular map of the genome from *L. manosquensis* strain Marseille-Q4368T obtained using the CGView Server online tool [30].</figcaption></figure>

<table><caption>**Table 3.** Comparative genomic data including ANI (top right) and dDDH (bottom left) values between *Leucobacter manosquensis* Marseille-Q4368 (Lman) and its closely related species such as *Leucobacter celer* NAL101 (Lcel), *Leucobacter triazinivorans* JW-1 (Ltri), *Leucobacter massiliensis* 122RC15 (Lmas), *Leucobacter zeae* CC-MF41(Lzea), *Leucobacter luti* RF9 (Llut), *Leucobacter chromiiresistens* JG 31 (Lchr), *Leucobacter japonicus* CBX 130 (Ljap), and *Leucobacter musarum* CBX 152 (Lmus).</caption><thead><tr><th></th><th>Lari</th><th>Lcel</th><th>Ltri</th><th>Lmas</th><th>Lzea</th><th>Llut</th><th>Lchr</th><th>Lman</th><th>Ljap</th><th>Lmus</th></tr></thead><tbody><tr><td>Lari</td><td></td><td>76.4%</td><td>76.6%</td><td>76.4%</td><td>76.4%</td><td>76.9%</td><td>76.5%</td><td>76.5%</td><td>75.1%</td><td>75.0%</td></tr><tr><td>Lcel</td><td>21.8%</td><td></td><td>85.7%</td><td>81.3%</td><td>79.7%</td><td>79.1%</td><td>79.5%</td><td>78.2%</td><td>78.2%</td><td>77.9%</td></tr><tr><td>Ltri</td><td>21.9%</td><td>30.1%</td><td></td><td>81.2%</td><td>79.8%</td><td>78.9%</td><td>79.8%</td><td>78.3%</td><td>77.9%</td><td>77.7%</td></tr><tr><td>Lmas</td><td>21.4%</td><td>24.9%</td><td>24.2%</td><td></td><td>78.9%</td><td>78.5%</td><td>79.3%</td><td>77.8%</td><td>77.4%</td><td>77.2%</td></tr><tr><td>Lzea</td><td>21.4%</td><td>23.3%</td><td>22.9%</td><td>22.7%</td><td></td><td>78.4%</td><td>79.4%</td><td>77.8%</td><td>77.6%</td><td>77.5%</td></tr><tr><td>Llut</td><td>21.8%</td><td>23.0%</td><td>22.9%</td><td>22.4%</td><td>22.2%</td><td></td><td>78.8%</td><td>77.4%</td><td>77.0%</td><td>77.0%</td></tr><tr><td>Lchr</td><td>21.7%</td><td>23.6%</td><td>23.3%</td><td>22.8%</td><td>22.8%</td><td>22.4%</td><td></td><td>84.3%</td><td>81.1%</td><td>81.0%</td></tr><tr><td>Lman</td><td>20.9%</td><td>22.6%</td><td>21.8%</td><td>22.0%</td><td>21.5%</td><td>21.5%</td><td>27.5%</td><td></td><td>79.8%</td><td>79.7%</td></tr><tr><td>Ljap</td><td>21.2%</td><td>23.3%</td><td>22.10%</td><td>21.6%</td><td>21.7%</td><td>21.3%</td><td>24.0%</td><td>22.7%</td><td></td><td>95.3%</td></tr><tr><td>Lmus</td><td>21.0%</td><td>23.0%</td><td>22.2%</td><td>21.7%</td><td>21.7%</td><td>21.2%</td><td>23.7%</td><td>22.6%</td><td>62.3%</td><td></td></tr></tbody></table>

In terms of genome size, the strain Marseille-Q4368 (3.18 Mbp) had approximately a similar genomic size compared to all the genomes of the *Leucobacter* genus. However, it possessed a smaller number of proteins compared to the other genomes, except for that of *L. massiliensis* (3.14 Mbp). Additionally, the tRNA molecules were more numerous than the rRNA and other RNA molecules, ranging between 46, 47, and 48 depending on the *Leucobacter* species analyzed. Furthermore, the difference in the G + C percentage between the species was very low, mostly being less than 4%. This reinforces our thesis that the Marseille-Q4368 strain is a new species within the *Leucobacter* genus. More statistical details from the genomes are showed in Table 4.
<table><caption>**Table 4.** Properties from genome analysis.</caption><thead><tr><th>Name</th><th>RefSeq</th><th>Size (Mbp)</th><th>G + C %</th><th>Protein</th><th>rRNA</th><th>tRNA</th><th>Other RNA</th><th>Pseudogene</th><th>Total Genes</th></tr></thead><tbody><tr><td>Leucobacter massiliensis 122RC15</td><td>MWZD00000000</td><td>3.14</td><td>71.0</td><td>2766</td><td>4</td><td>46</td><td>3</td><td>53</td><td>2872</td></tr><tr><td>Leucobacter manosquensis Marseille-Q4368</td><td>JAFEVO00000000</td><td>3.18</td><td>67.3</td><td>2776</td><td>9</td><td>47</td><td>3</td><td>149</td><td>2984</td></tr><tr><td>Leucobacter chromiiresistens JG 31</td><td>AGCW00000000</td><td>3.22</td><td>70.3</td><td>2919</td><td>9</td><td>48</td><td>3</td><td>34</td><td>3013</td></tr><tr><td>Leucobacter musarum subsp. musarum CBX152</td><td>JHBW00000000</td><td>3.44</td><td>66.8</td><td>3096</td><td>4</td><td>47</td><td>3</td><td>20</td><td>3590</td></tr><tr><td>Leucobacter aridicollis L-9</td><td>QYAE00000000</td><td>3.46</td><td>67.3</td><td>3115</td><td>4</td><td>47</td><td>2</td><td>0</td><td>3168</td></tr><tr><td>Leucobacter zeae CC-MF41</td><td>QYAB00000000</td><td>3.47</td><td>70.6</td><td>2914</td><td>4</td><td>47</td><td>3</td><td>137</td><td>3105</td></tr><tr><td>Leucobacter triazinivorans JW-1</td><td>CP035806</td><td>3.48</td><td>69.4</td><td>3048</td><td>6</td><td>46</td><td>3</td><td>38</td><td>3141</td></tr><tr><td>Leucobacter japonicus CBX130</td><td>JHBX00000000</td><td>3.59</td><td>66.8</td><td>3262</td><td>3</td><td>47</td><td>3</td><td>31</td><td>3346</td></tr><tr><td>Leucobacter luti RF6</td><td>QYAG00000000</td><td>3.62</td><td>69.5</td><td>3071</td><td>5</td><td>46</td><td>3</td><td>32</td><td>3156</td></tr><tr><td>Leucobacter celer subsp. astrifaciens CBX151</td><td>JHEI00000000</td><td>4.14</td><td>69.1</td><td>3589</td><td>3</td><td>48</td><td>3</td><td>43</td><td>3686</td></tr></tbody></table>
### 3.4. Description of *Leucobacter manosquensis* sp. nov.

*Leucobacter manosquensis* (ma.nos.quen'sis. N.L. masc. adj. manosquensis, "of Manosque," referring to the place where M&L Laboratories is located). Strain Marseille-Q4368 exhibited facultative aerobic growth, with a growth rate of 31 °C within 24 h. This particular bacterium was Gram-positive and had non-motile, non-spore-forming bacterial cells that were rod-shaped, measuring approximately 1.4 × 0.3 µm. 12-methyl-tetradecanoic acid (69%), 14-methyl-hexadecanoic acid (16%), and 14-methyl-pentadecanoic acid (7%) were the main branched fatty acids retrieved in strain Marseille-Q4368. Esculin ferric citrate, esterase (C4), esterase lipase (C8), leucine arylamidase, cystine arylamidase, valine arylamidase, naphthol-AS-BI-phosphohydrolase, D-mannitol, D-glucose, D-maltose, malic acid, and α-mannosidase were the present enzymes. It tested positive for a catalase test and negative for an oxidase test.

The genome of strain Marseille-Q4368T was 3.18 Mbp with 67.3 mol% of G + C content.

The genome and 16S rRNA sequence were deposited in GenBank under accession numbers MW583464 and JAFEVO000000000, respectively.

The type strain Marseille-Q4368T (=CSUR Q4368 = DSM 112403) was isolated from the forehead of a 59-year-old woman with healthy skin.

### 4. Conclusions

The skin is the largest organ of the human body, and because of its external nature, skin disorders have been the focus of medical research and cosmetic trials since ancient times. Having evolved together, scientific knowledge and cosmetic formulations are now inseparable. It should be noted that cosmetic and pharmaceutical regulations are particularly demanding, and are constantly updated in light of new scientific advances. The skin is colonized by a large number of organisms with which it evolves in symbiosis: the cutaneous microbiota. The number of symbiotic microorganisms is greater than the number of skin cells. Although it is not (yet) possible under actual European regulations to add living microorganisms to a formulation, growing knowledge on this subject already offers clear new prospects. Thus, the discovery of new bacteria could be a starting point for innovating component findings. Leucobacter manosquensis was discovered from a healthy skin swab sample. Both phylogenetic and phenotypic analyses revealed several different characteristics of the strain Marseille-Q4368 when compared to other members of the genus Leucobacter, which belongs to the Microbacteriaceae family and the Actinobacteria phylum and was first proposed in 1996 [9]. Leucobacter manosquensis strain Q4368 is the first species belonging to the Leucobacter genus that has been isolated from healthy human skin. Indeed, the dDDH and OrthoANI percentages of the strain Marseille-Q4368 are below the admitted threshold percentage of <70% and <95%, respectively, to define a new prokaryotic species [31]. Therefore, we propose Marseille-Q4368 as the type strain of a new species, within the genus Leucobacter, named Leubacter manosquensis. Gr. masc. adj. [λευκός] leukos, clear, light; N.L. masc. n. bacter, rod; N.L. masc. n. Leucobacter, colorless rod. Manosquensis, translitt. L. adj., "of Manosque," referring to the place where M&L Laboratories (one of the founders) is located.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11102535/s1; Figure S1: Phenotypic characteristics of Leucobacter manosquensis, (A) Gram staining (B), Colonies (C), Motility test; Figure S2: Distribution of the functional classes of predicted genes according to the clusters of orthologous groups of proteins of Marseille-Q4368 and its closely related species. Table S1: Digital DNA-DNA hybridization values obtained using a sequence comparison of all the studied genomes using the TYGS second value; Table S2: Cellular fatty acid composition (%) of Leucobacter manosquensis. Table S3: Differential phenotypic characteristics of Leucobacter manosquensis strain Marseille-Q4368 and closely related species; Table S4: Description of Leucobacter manosquensis strain Marseille-Q4368.
**Author Contributions:** Conceptualization: N.C. and B.L.S.; methodology: M.B., N.C. and B.L.S.; software: M.B. and C.I.L.; formal analysis: M.B., S.M., A.A., M.M.M., C.R. and C.I.L.; investigation and data curation: M.B., S.M., A.A., C.R. and M.M.M.; writing—original draft preparation: M.B. and N.C.; writing—review and editing: C.I.L. and N.C.; supervision: B.L.S. and N.C.; funding acquisition: B.L.S. and N.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was supported by the French state, managed by the National Research Agency under the "Investissements d'avenir (Investments for the Future)" program under the reference ANR-10-IAHU-03 (Méditerranée Infection), and by the Région Provence-Alpes-Côte-d'Azur and the European funding FEDER PRIMI.

**Data Availability Statement:** Data supporting reported results are available from the corresponding author.

**Acknowledgments:** The authors are indebted to Ludivine Brechard for sequencing the genome and Nicholas Armstrong and Céline Chartier for performing the fatty acid analyses.

**Conflicts of Interest:** The authors declare no conflict of interest. M.B. has a Ph.D. grant due to the collaboration between M&L Laboratories and Aix Marseille University reference PVM:2018-200.

## References

1. Boxberger, M.; Cenizo, V.; Cassir, N.; La Scola, B. Challenges in Exploring and Manipulating the Human Skin Microbiome. *Microbiome* **2021**, *9*, 125. [CrossRef]

2. Lagier, J.-C.; Edouard, S.; Pagnier, I.; Mediannikov, O.; Drancourt, M.; Raoult, D. Current and Past Strategies for Bacterial Culture in Clinical Microbiology. *Clin. Microbiol. Rev.* **2015**, *28*, 208–236. [CrossRef] [PubMed]

3. Lagier, J.-C.; Hugon, P.; Khelaifia, S.; Fournier, P.-E.; La Scola, B.; Raoult, D. The Rebirth of Culture in Microbiology through the Example of Culturomics to Study Human Gut Microbiota. *Clin. Microbiol. Rev.* **2015**, *28*, 237–264. [CrossRef] [PubMed]

4. Seng, P.; Drancourt, M.; Gouriet, F.; La Scola, B.; Fournier, P.-E.; 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. Off. Publ. Infect. Dis. Soc. Am.* **2009**, *49*, 543–551. [CrossRef] [PubMed]

5. Drancourt, M.; Bollet, C.; Carlioz, A.; Martelin, R.; Gayral, J.P.; Raoult, D. 16S Ribosomal DNA Sequence Analysis of a Large Collection of Environmental and Clinical Unidentifiable Bacterial Isolates. *J. Clin. Microbiol.* **2000**, *38*, 3623–3630. [CrossRef] [PubMed]

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

7. Bonnet, M.; Lagier, J.C.; Raoult, D.; Khelaifia, S. Bacterial Culture through Selective and Non-Selective Conditions: The Evolution of Culture Media in Clinical Microbiology. *N. Microbes N. Infect.* **2020**, *34*, 100622. [CrossRef]

8. Tidjani Alou, M.; Naud, S.; Khelaifia, S.; Bonnet, M.; Lagier, J.-C.; Raoult, D. State of the Art in the Culture of the Human Microbiota: New Interests and Strategies. *Clin. Microbiol. Rev.* **2020**, *34*, e00129-19. [CrossRef]

9. Takeuchi, M.; Weiss, N.; Schumann, P.; Yokota, A. Leucobacter Komagatae Gen. Nov., Sp. Nov., a New Aerobic Gram-Positive, Nonsporulating Rod with 2,4-Diaminobutyric Acid in the Cell Wall. *Int. J. Syst. Bacteriol.* **1996**, *46*, 967–971. [CrossRef]

10. Sasser, M. *Bacterial Identification by Gas Chromatographic Analysis of Fatty Acids Methyl Esters (GC-FAME)*; MIDI-Inc.: Newark, DE, USA, 2006.

11. Dione, N.; Sankar, S.A.; 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. *N. Microbes N. Infect.* **2016**, *10*, 66–76. [CrossRef]

12. Tittsler, R.P.; Sandholzer, L.A. The Use of Semi-Solid Agar for the Detection of Bacterial Motility. *J. Bacteriol.* **1936**, *31*, 575–580. [CrossRef] [PubMed]

13. Tatusova, T.; DiCuccio, M.; Badretdin, A.; Chetvernin, V.; Nawrocki, E.P.; Zaslavsky, L.; Lomsadze, A.; Pruitt, K.D.; Borodovsky, M.; Ostell, J. NCBI Prokaryotic Genome Annotation Pipeline. *Nucleic Acids Res.* **2016**, *44*, 6614–6624. [CrossRef] [PubMed]

14. Meier-Kolthoff, J.P.; Göker, M. TYGS Is an Automated High-Throughput Platform for State-of-the-Art Genome-Based Taxonomy. *Nat. Commun.* **2019**, *10*, 2182. [CrossRef] [PubMed]

15. Gupta, S.K.; Padmanabhan, B.R.; Diene, S.M.; Lopez-Rojas, R.; Kempf, M.; Landraud, L.; Rolain, J.-M. ARG-ANNOT, a New Bioinformatic Tool To Discover Antibiotic Resistance Genes in Bacterial Genomes. *Antimicrob. Agents Chemother.* **2014**, *58*, 212–220. [CrossRef] [PubMed]

16. Ingle, D.; Valcanis, M.; Kuzevski, A.; Tauschek, M.; Inouye, M.; Stinear, T.; Levine, M.M.; Robins-Browne, R.M.; Holt, K.E. EcOH In Silico Serotyping of E. Coli from Short Read Data identifies limited novel O-loci but extensive diversity of O:H serotype combinations within and between pathogenic lineages. *Microb. Genom.* **2015**, *2*, e000064.

17. Feldgarden, M.; Brover, V.; Haft, D.H.; Prasad, A.B.; Slotta, D.J.; Tolstoy, I.; Tyson, G.H.; Zhao, S.; Hsu, C.-H.; McDermott, P.F.; et al. Validating the AMRFinder Tool and Resistance Gene Database by Using Antimicrobial Resistance Genotype-Phenotype Correlations in a Collection of Isolates. *Antimicrob. Agents Chemother.* **2019**, *63*, e00483-19. [CrossRef]
18. Carattoli, A.; Zankari, E.; García-Fernández, A.; Voldby Larsen, M.; Lund, O.; Villa, L.; Møller Aarestrup, F.; Hasman, H. In Silico Detection and Typing of Plasmids Using PlasmidFinder and Plasmid Multilocus Sequence Typing. *Antimicrob. Agents Chemother.* **2014**, *58*, 3895–3903. [CrossRef]

19. Bortolaia, V.; Kaas, R.S.; Ruppe, E.; Roberts, M.C.; Schwarz, S.; Cattoir, V.; Philippon, A.; Allesoe, R.L.; Rebelo, A.R.; Florensa, A.F.; et al. ResFinder 4.0 for Predictions of Phenotypes from Genotypes. *J. Antimicrob. Chemother.* **2020**, *75*, 3491–3500. [CrossRef]

20. McArthur, A.G.; Waglechner, N.; Nizam, F.; Yan, A.; Azad, M.A.; Baylay, A.J.; Bhullar, K.; Canova, M.J.; De Pascale, G.; Ejim, L.; et al. The Comprehensive Antibiotic Resistance Database. *Antimicrob. Agents Chemother.* **2013**, *57*, 3348–3357. [CrossRef]

21. Chen, L. VFDB: A Reference Database for Bacterial Virulence Factors. *Nucleic Acids Res.* **2004**, *33*, D325–D328. [CrossRef]

22. Afgan, E.; Baker, D.; van den Beek, M.; Blankenberg, D.; Bouvier, D.; Čech, M.; Chilton, J.; Clements, D.; Coraor, N.; Eberhard, C.; et al. The Galaxy Platform for Accessible, Reproducible and Collaborative Biomedical Analyses: 2016 Update. *Nucleic Acids Res.* **2016**, *44*, W3–W10. [CrossRef] [PubMed]

23. Ondov, B.D.; Treangen, T.J.; Melsted, P.; Mallonee, A.B.; Bergman, N.H.; Koren, S.; Phillippy, A.M. Mash: Fast Genome and Metagenome Distance Estimation Using MinHash. *Genome Biol.* **2016**, *17*, 132. [CrossRef] [PubMed]

24. Lagesen, K.; Hallin, P.; Rødland, E.A.; Stærfeldt, H.-H.; Rognes, T.; Ussery, D.W. RNAmmer: Consistent and Rapid Annotation of Ribosomal RNA Genes. *Nucleic Acids Res.* **2007**, *35*, 3100–3108. [CrossRef] [PubMed]

25. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and Applications. *BMC Bioinform.* **2009**, *10*, 421. [CrossRef] [PubMed]

26. 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 Bioinform.* **2013**, *14*, 60. [CrossRef] [PubMed]

27. Lee, I.; Ouk Kim, Y.; 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]

28. Lefort, V.; Desper, R.; Gascuel, O. FastME 2.0: A Comprehensive, Accurate, and Fast Distance-Based Phylogeny Inference Program: Table 1. *Mol. Biol. Evol.* **2015**, *32*, 2798–2800. [CrossRef]

29. Letunic, I.; Bork, P. Interactive Tree Of Life (ITOL) v5: An Online Tool for Phylogenetic Tree Display and Annotation. *Nucleic Acids Res.* **2021**, *49*, W293–W296. [CrossRef]

30. Grant, J.R.; Stothard, P. The CGView Server: A Comparative Genomics Tool for Circular Genomes. *Nucleic Acids Res.* **2008**, *36*, W181–W184. [CrossRef]

31. Abdallah, R.A.; Beye, M.; Diop, A.; Bakour, S.; Raoult, D.; Fournier, P.-E. The Impact of Culturomics on Taxonomy in Clinical Microbiology. *Antonie Van Leeuwenhoek* **2017**, *110*, 1327–1337. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.