Phylogenetic and phylogenomic analyses
The only complete 16S rRNA gene sequence (1,532 bp) extracted from the draft genome shared 99.9% similarity with the nearly complete sequences (1,505 bp) obtained by conventional Sanger sequencing, which confirms their authenticity. Strain WB101T exhibited the highest similarity to Rhodohalobacter barkolensis 15182T (93.5%), followed by Rhodohalobacter halophilus JZ3C29T (93.2%), ‘Rhodohalobacter mucosus’ 8A47T (92.1%), and other validated species of the family Balneolaceae (on the edge of or lower than 90.0%). The phylogenetic tree based on the neighbour-joining algorithm showed that strain WB101T was distinct from R. halophilus JZ3C29T, R. barkolensis 15182T, and ‘R. mucosus’ 8A47T within the genus Rhodohalobacter at a high confidence level (Fig. 1), which was also emphasized by additional MP, ML phylogenetic trees.
The draft genome of strain WB101T was 5,104,032 bp in size, and the DNA G + C content was 42.0 mol%. A total of 138 contigs were assembled from the raw sequences, and the N50 value and coverage depth were 123,641 and 275.0 ×, respectively. The detailed genomic information of strain WB101T and related type strains are listed in Table 1. Strain WB101T had a significantly larger genome than the members of the genus Rhodohalobacter, though they all share a similar number and contents of tRNAs and G + C mol%.
Table 1
Genomic data that differentiates strain WB101T from related type strains.
Genomic data | 1 | 2a | 3b | 4c |
Genome size (bp) | 5,104,032 | 3,908,029 | 3,597,295 | 3,122,722 |
DNA G + C (mol%) | 42.0 | 47.7 | 42.5 | 44.4 |
Number of coding genes | 4,233 | 3,219 | 3,015 | 2,727 |
Number of total pseudo genes | 24 | 9 | 5 | 26 |
Number of tRNAs | 38 | 40 | 39 | 38 |
Number of rRNAs (5S, 16S, 23S) | 3, 1, 1 | 1, 1, 1 | 2, 1, 1 | 1, 2, 1 |
Number of contigs | 138 | 17 | 7 | 92 |
N50 values | 123,641 | 417,940 | 1,807,017 | 84,288 |
Coverage depth | 275.0 × | 198.0 × | 544.0 × | 137.9 × |
1, Strain WB101T; 2, ‘R. mucosus’ 8A47T; 3, R. barkolensis 15182T; 4, R. halophilus JZ3C29T |
Data from: a, Wang et al. (2021); b, Han et al. (2018); c, Xia et al. (2017) |
The OrthoANIu, dDDH, AAI, POCP, and TETRA value between strain WB101T and ‘R. mucosus’ 8A47T, R. barkolensis 15182T, and R. halophilus JZ3C29T are shown in Table 2. The OrthoANIu values between each of these strains were far lower than 95.0–96.0%, the threshold to identify potential novel species (Richter and Rosselló-Móra 2009). Moreover, the ANIb and ANIm values were all below 90.0%. All the dDDH values were below the threshold value (70.0%) for species delineation, demonstrating that each pair formed deep lineages (Li et al. 2010). TETRA values were lower than 0.99, the threshold for novel species (Richter and Rosselló-Móra 2009). For the taxonomic boundaries for genera, all the AAI and POCP values between strain WB101T and each of three related type strains were over 60.0% and 50.0%, respectively, which were argued as genus boundaries (Rodriguez-R and Konstandtinidis 2014; Qin et al. 2014). Additionally, the same taxonomic status shown by phylogenetic trees was also demonstrated by phylogenomic tree based on GTDB (Fig. 2).
Table 2
Phylogenomic value between strain WB101T and related type strains.
Phylogenomic value | 1 | 2 | 3 |
OrthoANIu | 69.2% | 69.5% | 69.4% |
dDDH | 18.5% | 17.6% | 18.1% |
TETRA | 0.81 | 0.93 | 0.92 |
AAI | 68.0% | 68.6% | 68.6% |
POCP | 56.0% | 53.4% | 52.1% |
1, Strain WB101T vs ‘R. mucosus’ 8A47T; 2, Strain WB101T vs R. barkolensis 15182T; 3, Strain WB101T vs R. halophilus JZ3C29T |
The genome of three related type strains from: www.ncbi.nlm.nih.gov |
Genomic analyses
There was no known secondary metabolite cluster identified with high similarity according to antiSMASH. Compared to members of the genus Rhodohalobacter, strain WB101T had a non-ribosomal peptide synthetase cluster (NRPS), which probably suggested an unknown peptide synthesis pathway was encoded by its genome.
Based on the RAST program, the galactosylceramide and sulfatide metabolism pathways with higher copy numbers were encoded by the genome of strain WB101T than three related type strains. The results of analyses by KOALA (KEGG Orthology And Links Annotation) indicated that strain WB101T featured in ASR. Reduction of sulfate (+ 6) occurs within cells of strain WB101T with bifunctional enzyme CysN/CysC (cysNC), sulfate adenylyltransferase subunit 1 and 2 (cysN and cysD), adenylylsulfate kinase (cysC), phosphoadenosine phosphosulfate reductase (cysH), and sulfite reductase (sir). As an intermediate during this process, PAPS (3′-phosphoadenosine-5′-phosphosulfate) is reduced to sulfite. Then, sulfide is transformed to l-cysteine by O-acetyl serine-(thiol)-lyase (Schiff 1979). However, the absence of sulfate permease, encoded by the cysPUWA operon, suggests strain WB101T might not utilise extracellular sulfates (Kushkevych et al. 2020). According to the genomic data from GenBank, the assimilatory sulfate reduction pathways of the genus Rhodohalobacter were incomplete in ‘R. mucosus’ 8A47T and R. barkolensis 15182T, and even absent in R. halophilus JZ3C29T. Conversion among polysulfides can be achieved by sulfhydrogenase subunit gamma, beta, alpha, and delta (hydG, hydB, hydA, and hydD) (Ma et al. 2000). Dissimilatory or assimilatory nitrate reduction was incomplete within the genus Rhodohalobacter and strain WB101T, which was testified by the experiment results. Compared to members of the genus Rhodohalobacter, strain WB101T had a more equipped pathway of carbohydrate metabolism, including pentose and glucuronate interconversions, galactose metabolism, amino sugar, and nucleotide sugar metabolism. The ribokinase (rbsK) allows cells of strain WB101T to transform d-ribose 5-phosphate to d-ribose. Phosphatidylcholine (lecithin) could be synthesized by glycerone-P with glycerol-3-phosphate dehydrogenase (gpsA), acyl phosphate:glycerol-3-phosphate acyltransferase (plsY), 1-acyl-sn-glycerol-3-phosphate acyltransferase (plsC), lysophosphatidate acyltransferase (AGPAT1_2), lysophosphatidic acid acyltransferase / lysophosphatidylinositol acyltransferase (AGPAT3_4), lysophosphatidiate acyltransferase (AGPAT5), lysocardiolipin and lysophospholipid acyltransferase (LCLAT1, AGPAT8), lysophospholipid acyltransferase 1/2 (MBOAT1_2), lysophospholipid acyltransferase (LPT1, ALE1), TAG lipase / sterylester hydrolase / phospholipase A2 / LPA acyltransferase (TGL4), 1-acylglycerol-3-phosphate O-acyltransferase (LOA1), phosphatidate cytidylyltransferase (CDS1, CDS2, cdsA), and phosphatidylcholine synthase (pcs). In summary, compared to species of the genus Rhodohalobacter, strain WB101T had more numerous metabolic pathways, which explained why it had a relatively larger genome and demonstrated its potential application values.
Morphological, physiological, and biochemical characteristics
The cells of strain WB101T were rod-shaped, 0.2–0.5 µm wide and 0.6–3.7 µm long after 3-day-growth in MB (Supplementary Fig. S2). The strain was Gram-stain-negative, oxidase-positive, catalase-positive, and could not motile on the surface of plate medium by flagella or gliding. Optimal growth occurred at 35–37°C (range 25–42°C), 5.0–6.0% (w/v) NaCl (range 1.0–11.0%), and pH 7.5–8.0 (range 7.0–8.5). There formed visible colonies on MA with or without 0.1% (w/v) KNO3 in an anaerobic jar after three weeks. The nitrate reduction was negative, and no H2S was detected by lead acetate papers or TSI.
Casein, Tweens 20, 40, 60, and 80 were hydrolysed by strain WB101T, but DNA, alginate, starch, and carboxymethyl cellulose were not. It was susceptible to penicillin (10µg), ampicillin (10µg), carbenicillin (100µg), erythromycin (15µg), vancomycin (30µg), rifampicin (5µg), chloramphenicol (30µg), norfloxacin (30µg), ofloxacin (5µg), cefotaxime (30µg), clarithromycin (15µg), lincomycin (2µg), ceftriaxone (30µg), and polymyxin B (300µg), while resistant to tobramycin (10µg), tetracycline (30µg), neomycin (30µg), gentamycin (10µg), streptomycin (10µg), and kanamycin (30µg). Remarkably, according to available data of related taxa, including Aliifodinibius salicampi KHM44T (Cho et al. 2017), ‘Aliifodinibius salipaludis’ WN023T (Zhao et al. 2020), Aliifodinibius saliphilus ECH52T (Cho and Whang 2020), Rhodohalobacter barkolensis 15182T (Han et al. 2018), and ‘Rhodohalobacter mucosus’ 8A47T, members in the phylum Balneolota showed the same sensitivity to vancomycin as strain WB101T. However, vancomycin, a kind of glycopeptide antibiotics, is inherently inactive toward Gram-stain-negative bacterium because it cannot cross the bacterial outer membrane (Yarlagadda et al. 2016). Thus, it suggested that cells of these taxa were different from typical Gram-stain-negative bacterium.
According to the API 20NE kits, strain WB101T and three related strains were negative for nitrate reduction and indole production. Compared to R. barkolensis MCCC 1K03442T and R. halophilus JZ3C29T, strain WB101T had a stronger positive reaction for assimilation of capric acid. Strain WB101T showed positive activities for reactions of alkaline phosphatase, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase, cystine arylamidase, trypsin, α-chymotrypsin, acid phosphatase, naphthol-AS-BI-phosphohydrolase, α-galactosidase (weakly), β-galactosidase, β-glucuronidase (weakly), α-glucosidase, β-glucosidase, N-acetyl-β-glucosaminidase, and α-mannosidase (weakly) in API ZYM kits. Positives for β-galactosidase and N-acetyl-β-glucosaminidase differentiated strain WB101T from the three related type strains, and negatives for lipase (C14) and α-fucosidase was consistent in the genus Rhodohalobacter. The results of API 50CH indicated that acids were produced from d-arabinose, l-arabinose (weakly), d-ribose, d-xylose, l-xylose (weakly), d-glucose (weakly), d-fructose (weakly), d-mannose (weakly), l-sorbose, inositol, d-sorbitol, methyl-αd-mannopyranoside, N-acetyl glucosamine, arbutin, esculin ferric citrate, d-lyxose (weakly), d-tagatose, potassium 2-ketogluconate (weakly), and potassium 5-ketogluconate. Additionally, the positives for d-arabinose and d-sorbitol did not present in three related type strains. According to Biolog GEN III, strain WB101T oxidized dextrin, d-maltose, d-trehalose, d-cellobiose, gentiobiose, sucrose, d-turanose, stachyose, d-raffinose, α- d-lactose, d-melibiose, N-acetyl-d-glucosamine, N-acetyl-d-galactosamine, α-d-glucose, d-mannose, d-galactose, d-fucose, l-histidine, pectin, d-galacturonic acid, l-galactonic acid lactone, d-glucuronic acid, glucuronamide, α-keto-glutaric acid, and acetoacetic acid. The oxidation of d-glucuronic acid and glucuronamide were stronger than others.
The further morphological, physiological, and biochemical characteristics that distinguish strain WB101T from related type strains were summarized in Table 3.
Table 3
Characteristics that differentiate strain WB101T from related type strains.
Characteristic data | 1 | 2 | 3 | 4 |
Oxygen requirement | Facultatively anaerobic | Aerobic a | Aerobic b | Facultatively anaerobic c |
Temperature (°C) |
Growth range | 25–42 | 20–42 a | 10–40 b | 20–50 c |
Optimum | 35–37 | 37–40 a | 37 b | 40 c |
NaCl concentration (%, w/v) |
Growth range | 1.0–11.0 | 2.0–14.0 a | 0.5–4.0 b | 2.0–16.0 c |
Optimum | 5.0–6.0 | 4.0–6.0 a | 2.0–3.0 b | 8.0–10.0 c |
pH |
Growth range | 7.0–8.5 | 7.0–9.0 a | 7.0–8.0 b | 7.0–9.0 c |
Optimum | 7.5–8.0 | 7.5–8.0 a | 7.5 b | 7.5–8.5 c |
Hydrolysis of |
Starch | - | + a | + b | + c |
Casein | + | + a | - b | ND c |
Alginate | - | + a | + b | - c |
Tweens 20 | + | + a | - b | - c |
Tweens 40 | + | - a | - b | - c |
Tweens 60 | + | - a | ND b | ND c |
Tweens 80 | + | - a | - b | - c |
Acid production from |
d-Arabinose | + | - | - | - |
Inositol | + | - | - | w |
d-Sorbitol | + | - | - | - |
Enzyme activities |
Oxidase | + | + a | -b | -c |
β-Galactosidase | + | - | - | - |
β-Glucosidase | + | + | - | - |
N-Acetyl-β-glucosaminidase | + | - | - | - |
Oxidation of |
d-Galactose | + | - | - | - |
d-Fucose | + | - | - | - |
Glycerol | - | - | + | + |
Pectin | + | - | - | - |
d-Glucuronic acid | + | - | - | w |
1, Strain WB101T; 2, ‘R. mucosus’ 8A47T; 3, R. barkolensis MCCC 1K03442T; 4, R. halophilus JZ3C29T |
+, positive; w, weakly positive; −, negative; ND, no data available |
The data without particular indications are from this study |
Data from: a, Wang et al. (2021); b, Han et al. (2018); c, Xia et al. (2017) |
Chemotaxonomic analyses
The sole respiratory quinone of strain WB101T was menaquinone 7 (MK-7), which was consistent with species of the genus Rhodohalobacter. Diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), and phosphatidylcholine (PC) were the major polar lipids detected in strain WB101T. In addition, moderate or minor amounts of phosphatidylethanolamine (PE), four unidentified lipids (L1, L2, L3, and L4), and three glycolipids (GL1, GL2, and GL3) were present. Compared to three related type strains, cells of strain WB101T could synthesize lecithin, which was also identified in genomic analyses. The absence of aminolipid (AL1, AL2, and AL3) could differentiate strain WB101T from ‘R. mucosus’ 8A47T and R. barkolensis MCCC 1K03442T. Further details of the polar lipids content of strain WB101T and three related type strains were shown in Supplementary Fig. S3. High content of iso-C15:0 was detected in strain WB101T (53.8%), as well as the member of the genus Rhodohalobacter, which was considered as the major fatty acid (≥ 10.0%). As the minor fatty acids, C12:0 and C16:0 (≥ 1.0%, ≤ 10.0%) appeared in the cells of four strains. Detailed discrepancies between strain WB101T and related type strains were listed in Supplementary Table S1.