Introduction {#sec1}
============

Bacteria constitute an important and highly diversified taxonomic group within the life tree of living organisms. Decoding the bacterial diversity underlying their normal and pathogenic functions is fundamental [@bib1]. A high-throughput bacterial culture approach based on diversified culture conditions, known as culturomics, was designed to isolate as yet uncultured species to unveil human microbial diversity, and to complement 16S rRNA metagenomics [@bib2], [@bib3], [@bib4]. Furthermore, a new taxonomic strategy named taxono-genomics was developed to include the analysis of complete genome sequences in combination with phenotypic characteristics [@bib5]. Herein, we report a short description of strain Marseille-P2802, isolated from the urine of a man treated for diabetic nephropathy.

Isolation and growth conditions {#sec2}
===============================

We isolated from the urine of a 59-year-old man treated with chronic haemodialysis for diabetic nephropathy, a potential new bacterial strain that could not be identified by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The screening was performed on a Microflex LT spectrometer (Bruker Daltonics, Bremen, Germany), as previously described [@bib6]. Spectra obtained ([Fig. 1](#fig1){ref-type="fig"}) were imported and analysed using the [Biotyper]{.smallcaps} 3.0 software against the Bruker database, which was continually incremented with local MEPHI database. The strain was isolated from a human urine sample, after 7 days growth on 5% sheep blood-enriched Columbia agar (bioMérieux, Marcy l'Etoile, France) at 37°C in an anaerobic atmosphere (anaeroGEN; Oxoid, Dardilly, France).Fig. 1MALDI-TOF MS reference spectrum of *Varibaculum massiliense* sp. nov. The reference spectrum was generated by comparison of spectra from 12 individual colonies.Fig. 1

Phenotypic characteristics {#sec3}
==========================

Microcolonies with mean diameter 0.5 mm were entire edged, translucent, greyish and glistering. Cells were a Gram positive rod-shaped bacterium and were slightly curved, non-motile and non-spore-forming ([Fig. 2](#fig2){ref-type="fig"}). The Strain Marseille-P2802 exhibited neither catalase nor oxidase activities <http://www.mediterranee-infection.com/article.php?larub=280&titre=umrs-database> [@bib7]. A comparative study of the biochemical characteristics of this strain with other closely related *Varibaculum* species is presented in [Table 1](#tbl1){ref-type="table"}. The biochemical characteristics of the Marseille-P2802^T^ strain obtained using the API ZYM and 20A strips (bioMérieux) are presented in [Table 2](#tbl2){ref-type="table"}. The major fatty acid found for this strain was hexadecanoic acid (47%). The 18 carbons, mostly unsaturated structures, were also abundant and represented almost 50% of the total composition: 18:1n9 (22%), 18:2n6 (12%), 18:0 (12%) and 18:1n7 (4%) ([Table 3](#tbl3){ref-type="table"}).Fig. 2Scanning electron microscopy (SEM) of stained *Varibaculum massiliense* sp. nov. A colony was collected from agar and immersed in a 2.5% glutaraldehyde fixative solution. Then a drop of the suspension was directly deposited on a poly-[l]{.smallcaps}-lysine-coated microscope slide for 5 minutes and treated with 1% phosphotungstic acid aqueous solution (pH 2.0) for 2 minutes to increase SEM image contrast. The slide was gently washed in water; air-dried and examined in a tabletop SEM (Tecnai G20). Scales and acquisition settings are shown in the figure.Fig. 2Table 1Differential characteristics of *Varibaculum massiliense* sp. nov., *Varibaculum timonense*[@bib15], *Varibaculum cambriensis*[@bib16] and *Varibaculum anthropi*[@bib17]Table 1Properties*V. massilienseV. timonenseV. cambriensisV. anthropi*Cell diameter (μm)0.5--0.60.4--0.5NANAOxygen requirementanaerobicanaerobicanaerobicanaerobicGram stain+−++Salt requirement−−−−Motility−−−−Endospore formation−−−−Alkaline phosphatase++−−Catalase−−−−Indole−+−−Urease++−+β-galactosidase−−−−*N*-acetyl-glucosamine−−−−Arabinose−−−−Lipase (C8)+++NATrypsin++−NAMannose−−−−Mannitol−+−−Glucose−−++Maltose++++SourceUrine sampleStool sampleHuman sourcesClinical sample[^1]Table 2Phenotypic characterization of *Varibaculum massiliense* sp. nov., based on analytical profile index (API) testsTable 2TestsCharacteristicsResultsAPI ZYMAlkaline phosphatase+Esterase (C4)−Esterase lipase (C8)+Lipase (C14)−Leucine arylamidase+Valine arylamidase−Cystine arylamidase−Trypsin+α-chymotrypsin−Acid phosphatase+Naphthalo-AS-BI-phosphohydrolase+α-galactosidase−β-galactosidase−β-glucuronidase−α-glucosidase+β-glucosidase−*N*-acetyl-β-glucosaminidase−α-mannosidase−α-fucosidase−API 20AIndole producttion−Urease+Glucose−Mannitol−Lactose+Sucrose+Maltose+Salicin−Xylose+Arabinose−Gelatin−Esculin−Glycerol+Cellobiose−Mannose−Melezitose+Raffinose+Sorbitol−Rhamnose−Trehalose−Table 3Cellular fatty acid composition (%)Table 3Fatty acidsNameMean relative % [a](#tbl3fna){ref-type="table-fn"}16:00hexadecanoic acid46.5 ± 0.618:1n99-octadecenoic acid21.9 ± 0.518:2n69,12-octadecadienoic acid11.8 ± 0.418:00octadecanoic acid11.6 ± 0.218:1n711-octadecenoic acid3.6 ± 0.214:00tetradecanoic acid1.3 ± 0.217:00heptadecanoic acid1.2 ± 0.315:00pentadecanoic acidTR17:0014-methyl-hexadecanoic acidTR17:0015-methyl-hexadecanoic acidTR16:1n79-hexadecenoic acidTR17:1n710-heptadecenoic acidTR[^2]

Strain identification {#sec4}
=====================

In order to classify this bacterium, the 16S rRNA gene was amplified using the primer pair fD1 and rP2 (Eurogentec, Angers, France) and sequenced using the Big Dye® Terminator v1.1 Cycle Sequencing Kit and 3500xLGenetic Analyzer capillary sequencer (Thermofisher, Saint-Aubin, France), as previously described [@bib8]. The 16S rRNA nucleotide sequence was assembled and corrected using [CodonCode Aligner]{.smallcaps} software (<http://www.codoncode.com>). Strain Marseille-P2802 exhibited 98.6% 16S rRNA similarity with *Varibaculum cambriense* strain CCUG 44998 (GenBank accession number NR_042127), the phylogenetically closest species with standing in nomenclature ([Fig. 3](#fig3){ref-type="fig"}). We consequently proposed to classify strain Marseille-P2802 as a new species within the genus *Varibaculum* belonging to the family *Actinomycetaceae* within the phylum Actinobacteria.Fig. 3Phylogenetic tree highlighting the position of *Varibaculum massiliense* sp. nov. with regard to other closely related species. GenBank accession numbers of 16S rRNA are indicated in parentheses. Sequences were aligned using MUSCLE with default parameters, phylogenetic inferences were obtained using the maximum likelihood method and the MEGA 7 software. Bootstrap values obtained by repeating the analysis 1000 times to generate a majority consensus tree are indicated at the nodes. The scale bar indicates a 2% nucleotide sequence divergence.Fig. 3

Genome sequencing {#sec5}
=================

Genomic DNA was extracted using the EZ1 biorobot with the EZ1 DNA tissue kit (Qiagen, Hilden, Germany) and then sequenced on a MiSeq sequencer (Illumina Inc, San Diego, CA, USA) with the Nextera Mate Pair sample prep kit and Nextera XT Paired End (Illumina), as previously described [@bib9]. The assembly was performed using a pipeline containing several softwares ([Velvet]{.smallcaps} [@bib10], [Spades]{.smallcaps} [@bib11] and [Soap Denovo]{.smallcaps} [@bib12]), and trimmed data ([MiSeq]{.smallcaps} and [Trimmomatic]{.smallcaps} [@bib13] softwares) or untrimmed data (only [MiSeq]{.smallcaps} software). [GapCloser]{.smallcaps} was used to reduce assembly gaps. Scaffolds \<800 bp and scaffolds with a depth value lower than 25% of the mean depth were removed. The best assembly was selected by using different criteria (number of scaffolds, N50, number of N). The genome of strain Marseille-P2802 was 2.14 Mb long with a 52.3 mol% of G+C content. The degree of genomic similarity of strain Marseille-P2802 with closely related species was estimated using the [OrthoANI]{.smallcaps} software [@bib14]. [OrthoANI]{.smallcaps} values among closely related species ([Fig. 4](#fig4){ref-type="fig"}) ranged from 85.95% between *Varibaculum cambriense* strain DSM 15806 and *V. cambriense* strain DORA_20Q618 to 61.29% between *V. cambriense* strain DSM 15806 and *Varibaculum timonense* Marseille-P3369. When *Varibaculum massiliense* was compared with these closely related species, values ranged from 62.01% with *V. timonensis* to 81.26% with *V. cambriense* strain DORA_20Q618.Fig. 4Heatmap generated with [OrthoANI]{.smallcaps} values calculated using the OAT software between *Varibaculum massiliense* sp. nov. and other closely related species with standing in nomenclature.Fig. 4

Conclusion {#sec6}
==========

On the basis of unique phenotypic features, including the MALDI-TOF spectrum, a 16S rRNA sequence divergence \>1.3% and an [OrthoANI]{.smallcaps} value \<95% with the phylogenetically closest species with standing in nomenclature, we formally proposed strain Marseille-P2802 as the type strain of *Varibaculum massiliense* sp. nov., a new species within the genus *Varibaculum*.

Description of *Varibaculum massiliense* strain Marseille-P2802 sp. nov. {#sec6.1}
------------------------------------------------------------------------

Marseille-P2802 is the type strain of *Varibaculum massiliense* sp. nov. (mas.si.li.en\'sis, L. fem. adj., from *Massilia*, the Latin name of Marseille, where the strain was first cultivated). The strain grows strictly under anaerobic conditions at 37°C. The potential pathogenicity of the type strain Marseille-P2802 (= CSUR P2802 = DSM 103074) is unknown. This strain has a genome length of 2.14 Mb with a 52.3% G+C content. The 16S rRNA gene sequence and whole-genome shotgun sequence of Marseille-P2802 were deposited in GenBank under accession numbers LT576396 and FNWI00000000, respectively.

Nucleotide sequence accession number {#sec6.2}
------------------------------------

The 16S rRNA gene and genome sequences were deposited in GenBank under accession number LT576396 and FNWI00000000, respectively.

Deposit in culture collections {#sec6.3}
------------------------------

Strain Marseille-P2802 was deposited in two different strain collections under number (= CSURP2802 = DSM 103074).

Conflict of interest {#appsec1}
====================

None declared.

Funding sources {#appsec2}
===============

This study was supported by the Institut Hospitalo-Universitaire Méditerranée Infection, the French National Research Agency under the programme *Investissements d\'avenir*, reference ANR-10-IAHU-03, the Région Provence Alpes Côte d'Azur and European funding FEDER PRIMI.

The authors thank Aurelia Caputo for submitting the genomic sequence to GenBank.

[^1]: **+**, positive result; **−**, negative result; NA, data not available.

[^2]: Mean peak area percentage; TR, trace amounts \<1%.
