Strain identification
MALDI-TOF MS failed to identify strain Marseille-P3249. Therefore, 16S rRNA gene sequencing was performed and using a blast comparison against the NCBI nucleotide database, strain Marseille-P3249 exhibited a 98.3% sequence similarity with Gemella bergeri strain 617 − 93, being the phylogenetically closest species with standing in nomenclature (Fig. 1) (18). Thus, and according to Kim et al., this strain may be classified within a new bacterial species within the Gemella genus as it exhibits more than 1.35% sequence divergence with its phylogenetically closest species with a validly published name (19). A gel view performed to compare the mass spectra of strain Marseille-P3249 and its phylogenetically-close species confirmed the novelty of this strain with its unique peak profile (Fig. 2).
General Characteristics Of Strain Marseille-p3249
Cells from strain Marseille-P3249 were Gram-positive cocci. Colonies grew in both aerobic and anaerobic atmospheres at temperatures ranging between 25°C and 37°C in optimally at 37°C in aerobic conditions. This strain grew at a pH range between 6 and 8.5 and NaCl concentrations below 50 g/L. In aerobic conditions at 37°C strain Marseille-P3249 formed colonies after 24hrs on COS agar of 0.5 to 1.2 mm in diameter. Cells had an average diameter of 0.78 µm (Table 1, Fig. 3). Cells were not motile and non-spore forming.
Table 1
Differential characteristics of Gemella massiliensis strain Marseille-P3249T (GMA), Gemella assaccharolytica EU427463 (GAS) (22), Gemella cuniculi AJ251987 (GCU) (23), Gemella morbillorum L14327 (GMO) (24), Gemella bergeri Y13365 (GBE) (18) and Gemella sanguinis Y13364 (GSA) (25).
Properties
|
GMA
|
GAS
|
GCU
|
GMO
|
GBE
|
GSA
|
Cell diameter (µm)
|
0.78
|
0.5
|
Na
|
0.3–0.8
|
Na
|
Na
|
Oxygen requirement
|
Fa
|
Fa
|
Fa
|
Fa
|
Fa
|
Fa
|
Gram stain
|
+
|
V
|
+
|
+
|
+
|
+
|
Endospore formation
|
-
|
-
|
Na
|
Na
|
-
|
-
|
Production of
|
|
|
|
|
|
|
Alkaline phosphatase
|
-
|
-
|
+
|
Na
|
-
|
+
|
Catalase
|
-
|
-
|
-
|
-
|
-
|
-
|
Urease
|
-
|
-
|
-
|
Na
|
-
|
-
|
β-galactosidase
|
-
|
Na
|
Na
|
Na
|
-
|
-
|
N-acetyl-β-glucosamine
|
-
|
Na
|
-
|
Na
|
Na
|
-
|
L-Arabinose
|
-
|
-
|
-
|
Na
|
-
|
-
|
D-Ribose
|
-
|
-
|
-
|
Na
|
-
|
-
|
D-Mannose
|
-
|
-
|
Na
|
+
|
Na
|
Na
|
D-Mannitol
|
-
|
-
|
+
|
+
|
-
|
+
|
D-glucose
|
-
|
-
|
+
|
+
|
+
|
+
|
D-fructose
|
+
|
+
|
Na
|
-
|
-
|
Na
|
D-maltose
|
-
|
-
|
V
|
+
|
W
|
+
|
D-lactose
|
-
|
-
|
-
|
-
|
-
|
V
|
G + C content (mol%)
|
30.5
|
26.7
|
28.9
|
30.8
|
30.3
|
29.7
|
Habitat
|
Sputum sample
|
Clinical specimen
|
Abcess of a rabit
|
Clincal specimen
|
Clincal specimen
|
Clincal specimen
|
Fa = Facultative anaerobic; Na = data not available; V = Variable; W = weakly positive |
Using an API 50CH strip (bioMérieux), this strain was able to metabolize D-fructose, amygdaline and L-sorbose. As for API ZIM (bioMérieux), positive reactions were observed for esterase (C4), esterase lipase (C8), leucine arylamidase, phosphatase acid and naphtol phosphohydrolase. Using API 20A (bioMérieux), the strain showed a positive reaction for esculin ferric citrate only. Negative reactions were obtained with alkaline phosphatase, lipase (C14), valine arylamidase, cystine arylamidase, trypsin, α-chymotrypsin, α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, N-acetyl-β-glucosaminidase, α-mannosidase and α-fucosidase.
The major fatty acids were Hexadecanoic acid (34 %), 9-Octadecenoic acid (28 %), Octadecanoic acid (15 %) and 9,12-Octadecadienoic acid (13 %). A wide variety of other fatty acids were described with low abundances. Four of them, rarely detected as cellular fatty acids, were composed of longer aliphatic chains (C20 and C22) (Table 2).
Table 2
Cellular fatty acid composition (%).
Fatty acids
|
Name
|
Mean relative % (a)
|
16:00
|
Hexadecanoic acid
|
34.1 ± 0.3
|
18:1n9
|
9-Octadecenoic acid
|
27.6 ± 0.2
|
18:00
|
Octadecanoic acid
|
14.8 ± 0.1
|
18:2n6
|
9,12-Octadecadienoic acid
|
12.5 ± 0.4
|
18:1n7
|
11-Octadecenoic acid
|
2.3 ± 0.1
|
18:1n5
|
13-Octadecenoic acid
|
2.1 ± 0.1
|
14:00
|
Tetradecanoic acid
|
1.2 ± 0.1
|
17:00
|
Heptadecanoic acid
|
TR
|
15:00
|
Pentadecanoic acid
|
TR
|
15:0 anteiso
|
12-methyl-tetradecanoic acid
|
TR
|
16:1n7
|
9-Hexadecenoic acid
|
TR
|
(a) Mean peak area percentage; TR = trace amounts < 1 % |
Strain Marseille P3249 exhibited MICs (µg/mL) of 0.012, 0.016, 0.016, 0.016, 0.016, 0.19, > 6, 0.125, 0.03, 0.64, 0.032, 0.75, > 256 and > 256 for benzylpenicillin, amoxicillin, cefotaxime, ceftriaxone, imipenem, erythromycin, daptomycin, amikacin, rifampicin, minocycline, teicoplanin, vancomycin, metronidazole and colistin, respectively.
Genome Characteristics Of Strain Marseille-p3249
The genome was 1,804,813-bp long with a 30.5 mol% G + C content (Fig. 4). It is composed of 7 scaffolds (composed of 8 contigs). Of the 1,727 predicted genes, 1,677 were protein-coding genes and 50 were RNAs (5 genes were 5S rRNA, 2 genes were 16S rRNA, 2 genes were 23S rRNA, and 41 genes were tRNA genes). A total of 1 276 genes (76.09%) were assigned a putative function (by cogs or by NR blast). Twenty-six genes were classified as ORFans (1.55%). The remaining genes were annotated as hypothetical proteins (304 genes (18,13%)) (Table 3). The distribution of genes into COG functional categories is detailed in Table 4.
Table 3
Nucleotide gene content and gene count levels of strain Marseille-P3249T
Infromation
|
Value
|
%
|
Size (bp)
|
1,804,813
|
100
|
Number of GC
|
551,183
|
30.54
|
Number total of genes
|
1,727
|
100
|
Number total of protein genes
|
1,677
|
97.10
|
Number total of RNA Genes
|
50
|
2.89
|
Number total of TRNA Genes
|
41
|
2.37
|
Number total of RNA (5S, 16S, 23S) Genes
|
9
|
0.52
|
Coding sequence size
|
1,547,868
|
85.76
|
Coding sequence gene protein size
|
1,535,274
|
85.06
|
Coding sequence tRNA gene size
|
3,178
|
0.17
|
Coding sequence (5S, 16S, 23S) gene size
|
9,416
|
0.52
|
Number of protein coding gene
|
1,677
|
100
|
Number of protein associated to cogs
|
1,136
|
67.74
|
Number of protein NOT associated to cogs
|
541
|
32.25
|
Number of protein associated to orfan
|
26
|
1.55
|
Number of gene associated to resistance genes
|
1
|
0.05
|
Number of genes associated to virulence
|
369
|
22.00
|
Table 4
Number of genes associated with the 25 general COG functional categories
Code
|
Value
|
% of total
|
Description
|
[J]
|
175
|
10.44
|
Translation
|
[A]
|
0
|
0.00
|
Rna processing and modification
|
[K]
|
75
|
4.47
|
Transcription
|
[L]
|
64
|
3.82
|
Replication, recombination and repair
|
[B]
|
0
|
0.00
|
Chromatin structure and dynamics
|
[D]
|
18
|
1.07
|
Cell cycle control, mitosis and meiosis
|
[Y]
|
0
|
0.00
|
Nuclear structure
|
[V]
|
43
|
2.56
|
Defense mechanisms
|
[T]
|
37
|
2.21
|
Signal transduction mechanisms
|
[M]
|
48
|
2.86
|
Cell wall/membrane biogenesis
|
[N]
|
9
|
0.54
|
Cell motility
|
[Z]
|
0
|
0.00
|
Cytoskeleton
|
[W]
|
2
|
0.12
|
Extracellular structures
|
[U]
|
20
|
1.19
|
Intracellular trafficking and secretion
|
[O]
|
46
|
2.74
|
Posttanslational modification, protein turnover,chaperones
|
[X]
|
63
|
3.76
|
Mobilome: prophages, transposons
|
[C]
|
63
|
3.76
|
Energy production and conversion
|
[G]
|
74
|
4.41
|
Carbohydrate transport and metabolism
|
[E]
|
100
|
5.96
|
Amino acid transport and metabolism
|
[F]
|
61
|
3.64
|
Nucleotide transport and metabolism
|
[H]
|
58
|
3.46
|
Coenzyme transport and metabolism
|
[I]
|
43
|
2.56
|
Lipid transport and metabolism
|
[P]
|
68
|
4.05
|
Inorganic ion transport and metabolism
|
[Q]
|
9
|
0.54
|
Secondary metabolites biosynthesis, transport and catabolism
|
[R]
|
96
|
5.72
|
General function prediction only
|
[S]
|
71
|
4.23
|
Function unknown
|
|
541
|
32.26
|
Not in COGs
|
Genome Comparison
The draft genome sequence of strain Marseille-P3249T was larger than those of Gemella cuniculi DSM 15828T, Gemella sanguinis ATCC 700632T and Gemella haemolysans ATCC 10379T (1.86, 1.90 and 1.91 Mb respectively), but smaller than those of Gemella asaccharolytica WAL 1945JT, Gemella bergeri ATCC 700627T and Gemella morbillorum NCTC11323T (1.28, 1.60 and 1.75 Mb respectively) (Table 5).
Table 5
Genome information of the species involved in the genomic comparative analyses.
Species
|
Size (Mb)
|
GC (%)
|
Gene Content
|
Gemella asaccharolytica WAL 1945JT
|
1.28
|
26.6
|
1,251
|
Gemella bergeri ATCC 700627T
|
1.60
|
30.3
|
1,524
|
Gemella morbillorum NCTC11323T
|
1.75
|
30.7
|
1,622
|
Gemella massiliensis Marseille-P3249T
|
1.80
|
30.5
|
1,677
|
Gemella cuniculi DSM 15828T
|
1.86
|
28.9
|
1,687
|
Gemella haemolysans ATCC 10379T
|
1.91
|
30.8
|
1,710
|
Gemella sanguinis ATCC 700632T
|
1.90
|
29.6
|
1,861
|
Additionally, the G + C content of strain Marseille-P3249T is smaller than those of G. asaccharolytica WAL 1945JT, G. cuniculi DSM 15828T, G. sanguinis ATCC 700632T and G. bergeri ATCC 700627T (26.6, 28.9, 29.6 and 30.3%, respectively), but larger than those of G. morbillorum NCTC11323T and G. haemolysans ATCC 10379T (30.7 and 30.8 %, respectively).
The gene content of strain Marseille-P3249T was larger than those of G. asaccharolytica WAL 1945JT, G. bergeri ATCC 700627T and G. morbillorum NCTC11323T (1,251, 1,524 and 1,622 respectively), but smaller than those of G. cuniculi DSM 15828T, G. haemolysans ATCC 10379T and G. sanguinis ATCC 700632T (1,687, 1,710 and 1,861, respectively).
Strain Marseille-P3249 shared the highest number of orthologous proteins with G. cuniculi (1039). Furthermore, this bacterium shared 1031, 1032, 1054, and 778 orthologous proteins with G. haemolysans, G. morbillorum, G. sanguinis and G. asaccharolytica, respectively. Strain Marseille-P3249 exhibited OrthoANI values of 76.5, 76.8, 75.9, 75.4, 94.8 and 70.3% with G. morbillorum, G. sanguinis, G. haemolysans G. cuniculi, G. bergeri and G. asaccharolytica, respectively (Fig. 5).
Strain Marseille-P3249 exhibited dDDH values of 21.3 ± 4.7%, 22.6 ± 4.7%, 21.7 ± 4.7%, 22.1 ± 4.7%, 21.9 ± 4.7% and 59.7 ± 5.6% with G. asaccharolytica, G. cuniculi, G. haemolysans, G. morbillorum, G. sanguini and G. bergeri (Table 6). These results confirm the novelty of the isolated strain, since 70% is the recommended dDDH threshold to delimitate a new bacterial species (20, 21).
Table 6
Digital DNA-DNA hybridization values (%) obtained by strain Marseille-P3249T with other closely-related species using the GGDC formula 2 software (dDDH estimates based on identities / HSP length).
|
GMA
|
GAS
|
GCU
|
GHA
|
GMO
|
GSA
|
GBE
|
GMA
|
100%
|
21.30 ± 4.7%
|
22.60 ±4.7%
|
21.70 ±4.7%
|
22.10 ± 4.7%
|
21.90 ±4.7%
|
59.70 ±5.6%
|
GAS
|
|
100%
|
23.40 ±4.7%
|
23.20 ±4.7%
|
22.40 ± 4.7%
|
21.60 ±4.6%
|
21.00 ±4.7%
|
GCU
|
|
|
100%
|
21.80 ±4.7%
|
22.00 ± 4.7%
|
22.10 ±4.7%
|
22.70 ±4.7%
|
GHA
|
|
|
|
100%
|
22.90 ± 4.7%
|
23.50 ±4.8%
|
22.10 ± 4.7%
|
GMO
|
|
|
|
|
100%
|
23.00 ±4.8%
|
21.90 ±4.7%
|
GSA
|
|
|
|
|
|
100%
|
21.90 ±4.8%
|
GBE
|
|
|
|
|
|
|
100%
|
Abbreviations : GMA, Gemella massiliensis Marseille-P3249; GAS, Gemella asaccharolytica strain KA00071; GCU, Gemella cuniculi DSM 15828; GHA, Gemella haemolysans strain NCTC10459; GMO, Gemella morbillorum strain NCTC11323; GSA, Gemella sanguinis strain FDAARGOS 742; GBE, Gemella bergeri ATCC 700627. |
Description of Gemella massiliensis sp. nov.
We propose strain Marseille-P3249 is the type strain of the new species Gemella massiliensis sp. nov. (mas.il.i.en’sis, L. gen. fem., adj., massiliensis, pertaining to Massilia, the Latin name of the city of Marseille, where this bacterium was discovered). Strain Marseille-P3249 is a facultative anaerobic bacterium but grows optimally at 37°C under aerobic conditions. Using a 50 CH strip, this strain exhibits positive reactions for D-fructose, amygdaline and L-sorbose. Positive reactions are also observed for esterase (C4), esterase lipase (C8), leucine arylamidase, phosphatase acid and naphtol phosphohydrolase using an API strip. In addition, using an API 20A (bioMérieux), positive reactions are observed for esculin ferric citrate only.
The major fatty acids are Hexadecanoic acid (34 %), 9-Octadecenoic acid (28 %), Octadecanoic acid (15 %) and 9,12-Octadecadienoic acid (13 %).
The genome is 1,804,813 bp long with 30.5 mol% G + C content. The 16S rRNA and whole genome sequences of G. massiliensis sp. nov., were deposited in EMBL-EBI under accession numbers LT628479 and FQLS00000000, respectively. The type strain Marseille-P3249T (= CSUR P3249 = DSMZ 103940) was isolated from the sputum sample of a healthy French man.