MLSA of the five concatenated HKGs
Of all 223 Proteus strains collected in this study, phylogenetic tree of concatenated 5-gene divided them into eleven clusters (Fig. 1), representing of thirteen species. Among them, ten clusters contained one type strain of each. However, cluster 5 was comprised of three type strains, i.e, Proteus genospecies ATCC51470T, P. cibarius JCM 30699T and P. terrae LMG 28659T.
Expectedly, among the 223 Proteus strains, P. mirabilis (cluster 1) is the largest cluster (n = 157, 70.4%) distinctly separated from others, and there are three subclusters within this cluster. Cluster 5 is the second largest among Proteus strains (n = 33, 14.8%), followed by P. vulgaris (cluster 2) (n = 14, 6.3%) and P. penneri (cluster 7) (n = 6, 2.7%).
Identification of Proteus species by phylogenetic analysis of five individual genes
Phylogenetic trees based on five individual HKGs were also constructed (Fig. 2). Phylogenetic trees of the five HKGs (dnaJ, mdh, pyrC, recA and rpoD) can be divided into eleven clusters, representative thirteen species, and corresponding to that of the concatenated tree. Meanwhile, phylogenetic trees of four individual HKGs (dnaJ, mdh, pyrC and rpoD) were the same as that of the concatenated tree, both in numbers of species (cluster) and strain numbers within each species (cluster). There is one inconsistence between trees of recA and concatenated 5-gene: recA identified four strains as uncluster, whereas the four strains were identified by concatenated 5-gene and other four HKGs as species 6 (Fig. 2).
Inter- And Intra-species Distance Of HKGs
The inter- and intra-species distances of HKGs were summarized in boxplot of the concatenated 5-gene (Fig. 3). All inter-species indicated the obviously different distances than that of intra-species. Among the inter-species boxplot, two species, P. mirabilis, and P. hauseri indicated compacted distance ranges (both standard deviation, SD = 0.004), whereas the remaining nine species shared dispersive distance ranges (SD ranges from 0.024 to 0.065); On the other hand, among intra-species boxplot, P. hauseri possessed compacted distance range (SD = 0.000) compared to that of five species (SD range from 0.012 to 0.058). Meanwhile, boxplots of the five individual genes (Figure S1) indicated same trends of intra and inter-species distance as that of the concatenated 5-gene, although there are small part overlapping in species 5 and 6 of pyrC. The detailed genetic distance and median values of individual genes and the concatenated 5-genes were summarized in Table S1.
Web-based DNA-DNA Hybridizations Among Species
To confirm the correctness of strains among the eleven species, we used web-based DDH, like is DDH and ANI to detect their similarity values. Among the eleven species defined in this study, the is DDH and ANI values of the type /representative strains were 23.5–51.4% and 80.8–94.4% (Table 2), less than the proposed cutoff level for species delineation, i.e., 70% and 95%, respectively. Notably, among the three subclusters within cluster 5 (Fig. 1), either among the three published type strains (Proteus genospecies ATCC51470T, P. cibarius JCM30699T and P. terrae LMG28659T) or representative strain (CA142267) among the three subclusters, their is DDH and ANI values were more than the proposed cutoff level for species delineation. The results indicated strains within the cluster 5 actually belonged to same species.
Table 2
is DDH relatedness and ANI values among the eleven species.
Strains # | a | b | c | d | e | f | g | h | i | j | k | l | m | n |
Species | P. mirabilis | P. penneri | P. vulgaris | P.hauseri | Proteus genospecies 4 | Proteus genospecies 5 | P. cibarius | P. terrae | CA142267 | Proteus genospecies 6 | P. columbae | P. alimentorum | P. faecis | P. cibi |
DDH/ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI | DDH | ANI |
Strains# | a | 100.0 | 100.0 | | | | | | | | | | | | | | | | | | | | | | | | | | |
b | 24.7 | 82.1 | 100.0 | 100.0 | | | | | | | | | | | | | | | | | | | | | | | | |
c | 24.4 | 81.7 | 37.5 | 89.4 | 100.0 | 100.0 | | | | | | | | | | | | | | | | | | | | | | |
d | 23.5 | 80.8 | 25.8 | 83.3 | 25.1 | 82.6 | 100.0 | 100.0 | | | | | | | | | | | | | | | | | | | | |
e | 24.6 | 81.8 | 41.5 | 90.9 | 44.5 | 91.6 | 25.8 | 83.2 | 100.0 | 100.0 | | | | | | | | | | | | | | | | | | |
f | 24.3 | 81.8 | 44.6 | 91.7 | 37.7 | 89.6 | 25.8 | 83.1 | 40.0 | 90.2 | 100.0 | 100.0 | | | | | | | | | | | | | | | | |
g | 25.1 | 82.6 | 44.5 | 91.7 | 37.5 | 89.6 | 25.7 | 83.2 | 39.9 | 90.2 | 72.3 | 96.9 | 100.0 | 100.0 | | | | | | | | | | | | | | |
h | 24.1 | 81.5 | 44.6 | 91.8 | 37.5 | 89.5 | 25.7 | 83.3 | 39.9 | 90.2 | 93.1 | 99.2 | 71.2 | 96.7 | 100.0 | 100.0 | | | | | | | | | | | | |
i | 24.4 | 82.0 | 44.9 | 91.8 | 37.6 | 89.7 | 25.8 | 83.1 | 39.8 | 90.3 | 74.4 | 97.1 | 72.3 | 96.8 | 72.6 | 96.8 | 100.0 | 100.0 | | | | | | | | | | |
j | 25.0 | 82.4 | 44.8 | 91.7 | 36.0 | 89.0 | 26.1 | 83.4 | 37.8 | 89.5 | 47.3 | 92.3 | 47.1 | 92.5 | 47.3 | 92.4 | 48.3 | 92.8 | 100.0 | 100.0 | | | | | | | | |
k | 24.8 | 81.9 | 46.3 | 92.1 | 36.8 | 89.1 | 26.0 | 83.4 | 37.9 | 89.6 | 48.9 | 92.7 | 48.3 | 92.6 | 48.8 | 93.0 | 49.8 | 93.0 | 57.1 | 94.4 | 100 | 100 | | | | | | |
l | 24.9 | 82.0 | 46.3 | 92.1 | 36.4 | 89.0 | 25.9 | 83.3 | 37.8 | 89.6 | 48.7 | 92.8 | 48.1 | 92.7 | 48.5 | 92.7 | 49.2 | 92.3 | 52.4 | 93.7 | 53.9 | 93.9 | 100 | 100 | | | | |
m | 24.1 | 81.6 | 52.7 | 93.7 | 35.8 | 88.7 | 25.7 | 83.1 | 36.5 | 89.1 | 45.9 | 91.9 | 45.5 | 91.9 | 45.9 | 92.0 | 46.6 | 92.1 | 47.1 | 92.2 | 48.6 | 92.7 | 48.5 | 92.6 | 100 | 100 | | |
n | 24.6 | 81.9 | 45.3 | 92.0 | 35.5 | 88.3 | 25.8 | 83.0 | 36.5 | 89.0 | 45.8 | 92.0 | 45.5 | 92.0 | 45.9 | 92.0 | 46.6 | 92.3 | 49.2 | 92.0 | 51.4 | 93.3 | 50.5 | 93 | 49.5 | 92.9 | 100 | 100 |
# Strain: a, P. mirabilis ATCC 29906T(GenBank accession no. ACLE00000000.1); b, P. penneri ATCC 33519T (PHFJ00000000); c, P. vulgaris KCTC 2579T (PHNN000000000); d, P.hauseri JCM 1668T (PGWU00000000); e, Proteus genospecies 4 ATCC 51469T (PENV00000000); f, Proteus genospecies 5 ATCC 51470T (PENU00000000); g, P. cibarius JCM 30699T (PGWT00000000); h, P. terrae LMG 28659T (PENS00000000); i, CA142267; j, Proteus genospecies 6 ATCC 51471T (PENT00000000); k, P. columbae 08MAS2615T (NGVR00000000); l, P. alimentorum 08MAS0041T (NBVR00000000); m, P. faecis TJ1636T (PENZ00000000); n, P. cibi FJ2001126-3T (PENW00000000). Results were percentages based on Formula 2, calculate distances and DDH estimates with GGDC 2; ANI values were estimated using the web-based service ANI calculator (http://www.ezbiocloud.net/tools/ani). The grey shadow marked strains with is DDH > 70% and ANI > 95%, respectively, indicating they belong to the same species. |
Reclassification of Proteus genospecies 5 and P. cibarius to P. terrae
Since either MLSA of the five concatenated HKGs or phylogenetic analysis of five individual genes indicated that three type strains, i.e, Proteus genospecies ATCC51470T, P. cibarius JCM 30699T and P. terrae LMG 28659T were fell into one cluster (cluster 5 in Fig. 1), further web-based DNA-DNA hybridizations, like is DDH and ANI confirmed that, among the three subclusters within cluster 5, either among the three type strains or representative strain (CA142267) among the three subclusters, their is DDH and ANI values were higher than the proposed cutoff level for species delineation (70% for is DDH and 95% for ANI, Table 2). The genomic analysis provided evidences that strains within the cluster 5 actually belonged to same species.
Further phenotypic characteristics were detected among type strains of Proteus genospecies 5, P. cibarius and P. terrae, and slight distinctive properties were observed (Table 3). Only minor differences were obtained between type strains of the three species, this including growth in optimum temperature, growth range in NaCl and pH, the utilization of Dnase, lipase and citric acid, and DNA G + C content. Combined with the genetic, genomic and phenotypic characteristics, three species, P. terrae reported by Behrendt et al. 2015, P. cibarius reported by Hyun et al. 2016 and Proteus genospecies 5 reported by O’Hara et al. 2000 should be regarded as the heterotypic synonyms of Proteus terrae reported by Behrendt et al. 2015.
Table 3
Distinctive phenotypic characteristics among type strains P. terrae, P. cibarius and Proteus genospecies 5#.
Characteristic | P. terrae | P. cibarius | Proteus genomospecies 5 |
Growth in optimum temperature (°C) | 37 | 35 | 37 |
Growth range in NaCl (%,w/v) | 0–15 | 0–12 | 0–15 |
Growth range in pH | 4–9 | 4–9 | 4–9 |
Dnase (25 °C) (3 days) | + | + | + |
Lipase (olive oil) (7 days) | - | + | - |
CIT | - | - | + |
DNA G + C content (mol %) | 37.9 | 37.8 | 37.8 |
#Species strain: P. terrae LMG 28659T; P. cibarius JCM 30699T; Proteus genospecies 5 ATCC 51470T. |