Next Article in Journal
Acinetobacter baumannii IC2 and IC5 Isolates with Co-Existing blaOXA-143-like and blaOXA-72 and Exhibiting Strong Biofilm Formation in a Mexican Hospital
Next Article in Special Issue
Characterization and Genomic Analysis of Fererhizobium litorale gen. nov., sp. nov., Isolated from the Sandy Sediments of the Sea of Japan Seashore
Previous Article in Journal
Feed Regime Slightly Modifies the Bacterial but Not the Fungal Communities in the Intestinal Mucosal Microbiota of Cobia Fish (Rachycentron canadum)
Previous Article in Special Issue
Growing in Saltwater: Biotechnological Potential of Novel Methylotuvimicrobium- and Methylomarinum-like Methanotrophic Bacteria
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome Analysis and Potential Ecological Functions of Members of the Genus Ensifer from Subsurface Environments and Description of Ensifer oleiphilus sp. nov.

by
Alexey P. Ershov
1,
Tamara L. Babich
1,
Denis S. Grouzdev
2,
Diyana S. Sokolova
1,
Ekaterina M. Semenova
1,
Alexander N. Avtukh
3,
Andrey B. Poltaraus
4,
Elena A. Ianutsevich
1 and
Tamara N. Nazina
1,*
1
Winogradsky Institute of Microbiology, Research Center of Biotechnology of the Russian Academy of Sciences, 119071 Moscow, Russia
2
SciBear OU, 10115 Tallinn, Estonia
3
Skryabin Institute of Biochemistry and Physiology of Microorganisms, Russian Academy of Sciences, Pushchino Scientific Center for Biological Research of the Russian Academy of Sciences, 142290 Pushchino, Russia
4
Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Microorganisms 2023, 11(9), 2314; https://doi.org/10.3390/microorganisms11092314
Submission received: 24 July 2023 / Revised: 6 September 2023 / Accepted: 9 September 2023 / Published: 14 September 2023

Abstract

:
The current work deals with genomic analysis, possible ecological functions, and biotechnological potential of two bacterial strains, HO-A22T and SHC 2-14, isolated from unique subsurface environments, the Cheremukhovskoe oil field (Tatarstan, Russia) and nitrate- and radionuclide-contaminated groundwater (Tomsk region, Russia), respectively. New isolates were characterized using polyphasic taxonomy approaches and genomic analysis. The genomes of the strains HO-A22T and SHC 2-14 contain the genes involved in nitrate reduction, hydrocarbon degradation, extracellular polysaccharide synthesis, and heavy metal detoxification, confirming the potential for their application in various environmental biotechnologies. Genomic data were confirmed by cultivation studies. Both strains were found to be neutrophilic, chemoorganotrophic, facultatively anaerobic bacteria, growing at 15–33 °C and 0–1.6% NaCl (w/v). The 16S rRNA gene sequences of the strains were similar to those of the type strains of the genus Ensifer (99.0–100.0%). Nevertheless, genomic characteristics of strain HO-A22T were below the thresholds for species delineation: the calculated average nucleotide identity (ANI) values were 83.7–92.4% (<95%), and digital DNA–DNA hybridization (dDDH) values were within the range of 25.4–45.9% (<70%), which supported our conclusion that HO-A22T (=VKM B-3646T = KCTC 92427T) represented a novel species of the genus Ensifer, with the proposed name Ensifer oleiphilus sp. nov. Strain SHC 2-14 was assigned to the species ‘Ensifer canadensis’, which has not been validly published. This study expanded the knowledge about the phenotypic diversity among members of the genus Ensifer and its potential for the biotechnologies of oil recovery and radionuclide pollution treatment.

1. Introduction

The taxonomic position of the bacterial species of the genera Ensifer and Sinorhizobium has been repeatedly revised over the past two decades [1,2]. These bacteria form a single phylogenetic group within the family Rhizobiaceae, referred to as the Ensifer/Sinorhizobium clade. A number of suggestions have been made to transfer the previously described Ensifer species to the genus Sinorhizobium and vice versa [3,4]. Recent data revealed that despite the high phylogenetic similarity of all detected strains (cpAAI > 86%), some species from the clade differ significantly from the rest in their biology and ecology [5]. This fact has made it possible to separate these genera among themselves, to assign 17 validly described bacterial species to the genus Sinorhizobium and the species E. adhaerens, E. morelensis, E. sesbaniae, and ‘E. canadensis’ to the genus Ensifer [5,6]. Members of the genus Ensifer, originally described by Casida [7], are budding bacteria, able to grow in LB medium and to hydrolyze starch, with an optimal temperature at 27–28 °C; they are generally nonsymbiotic (except for symbiosis of E. sesbaniae with legumes and that of ‘E. canadensis’ with roots of Melilotus albus), slightly acidophilic, and resistant to ampicillin and erythromycin [5,6]. Some members of the genus Ensifer are applicable in biotechnologies of xenobiotics degradation [8,9].
Two bacterial strains, HO-A22T and SHC 2-14, were isolated from an oil field water and radionuclide-contaminated groundwater (Russia), respectively. Microorganisms isolated from petroleum reservoirs are of great potential for application in microbial enhancement of oil recovery (MEOR) technologies. Some of these technologies are based on the production of exopolysaccharides by the reservoir microbiota [10], which leads to the redirection of injected water flows and, as a result, an increase in the coverage of the oil reservoir by flooding. Several bacterial strains from the Ensifer/Sinorhizobium clade are known to release significant amounts of exopolysaccharides into the environment [11,12,13,14], and,, therefore may be useful for MEOR technologies in oil fields. Moreover, injection of nitrate into petroleum reservoirs can activate the growth of denitrifying bacteria, which leads to suppression of the sulfidogenic microbiota and, therefore, decreases the risks of oil field equipment biocorrosion [15,16,17]. Isolates of the genera Ensifer and Sinorhizobium are known to reduce nitrate ions during their growth [18,19,20], making it potentially possible to stimulate bacteria of this group in the formation water in order to control the petroleum reservoir souring.
Bacteria of the genera Ensifer and Sinorhizobium, capable of denitrification, can affect the mobility of variable-valency metals and radionuclides in groundwater contaminated with components of liquid radioactive waste in the areas of the nuclear fuel cycle enterprises. When liquid radioactive waste is buried in surface water bodies and underground layers, toxic components gradually seep into the groundwater [21,22]. Liquid radioactive wastes are enriched in nitrates and organic salts. In the process of denitrification, due to the reduction of nitrates to molecular nitrogen, the ambient redox potential decreases [23], which leads to the transition of radionuclides and heavy metals into a poorly soluble reduced form. This process can be used to create a geochemical barrier for radionuclides in groundwater contaminated with toxic components of liquid radioactive waste.
The microorganisms can also have an effect on radionuclides and heavy metals due to their immobilization in the exopolysaccharide matrix of bacterial biofilms [23,24], in the processes of biosorption and bioaccumulation on the surface or inside the cells [25,26], and in the processes of direct reduction of variable-valency elements to poorly soluble compounds during cellular respiration [22,24,27,28,29]. Thus, the study of microorganisms isolated from radionuclide-contaminated groundwater, as well as from oil reservoirs, has both basic and applied importance.
The aim of our study was to estimate the biotechnological potential of strains HO-A22T and SHC 2-14 belonging to the genus Ensifer based on its genomic and phenotypic characterization and taxonomic description of the strain HO-A22T as Ensifer oleiphilus sp. nov. The study of the strain SHC 2-14 allowed us to assign it to the not validly published species ‘Ensifer canadensis’ and to expand knowledge about the phenotypic diversity of this species. The ecological values of the isolated strains and their potential for the biotechnologies associated with oil recovery and radionuclide pollution treatment were also discussed.

2. Materials and Methods

2.1. Sources of Isolation of Bacterial Strains

Strain HO-A22T was isolated from a sample of injection water represented by a mixture of fresh river water and production water reinjected into the Cheremukhovskoe oil field (Tatarstan, Russia). The water sample was collected in 2016 during a field experiment of a microbial enhanced oil recovery (MEOR) technology. In the terrigenous Cheremukhovskoe heavy oil reservoir, formation water belonged to the chlorine-calcium type. pH of the water sample was 7.86, and its density (at 25 °C) was 997.4 kg·m–3. The injection water contained the following dissolved ions (mg·L–1): hydrocarbonate (244), chloride (114), potassium and sodium (114), sulfate (74), calcium (60), and magnesium (15). The total salinity of the water sample was 620 mg·L–1. The depth of the reservoir was 850–1300 m below sea level, and its average temperature was 22–25 °C. Other physicochemical characteristics of the water sample were described by Nazina et al. [30].
Strain SHC 2-14 was isolated from a groundwater sample collected in 2015 through an observation well from an upper water-bearing horizon at the depth of 10 m near a preserved surface repository for radioactive waste (Tomsk region, Russia) [21]. The temperature of the horizon was 6–8 °C; the groundwater sample had total salinity of 7 g·L–1 and pH 5.9. Redox potential of the groundwater varied from –20 to +150 mV (from +50 to +70 mV on average), so it might be characterized as microaerobic. The sample contained dissolved ions at the following concentrations (mg·L–1): nitrate (4239), calcium (772), sodium (570), hydrocarbonate (170), magnesium (146), carbonate (135), sulfate (133), etc. The presence of significant amounts of strontium, cobalt, cesium, and plutonium in the sample was observed.
Strains HO-A22T and SHC 2-14 were cultivated on the TEG medium containing bacto-trypton (5.0 g), yeast extract (2.5 g), glucose (1.0 g), NaCl (5.0 g), and distilled water (1 L, pH 7.0–7.2). All reagents were purchased from Dia-M (Moscow, Russia) and were chemically pure. Isolation of the strains was performed by serial dilution technique as described before [29]. To obtain pure cultures, the medium was supplemented with bacteriological agar (20 g·L–1), and then single colonies were reinoculated in the liquid broth. The strains have been deposited at the All-Russian Collection of Microorganisms (VKM; Pushchino, Moscow Region, Russia) and at the Korean Collection for Type Cultures (KCTC; Jeongeup-si, Jeollabuk-do, South Korea) under the numbers VKM B-3646T and KCTC 92427T (strain HO-A22T) and VKM B-3385 and KCTC 92426 (strain SHC 2-14). For comparative studies, type strains E. adhaerens AT (=NBRC 100388T = LMG 20216T = ATCC 33212T) and E. morelensis Lc04T (=NBRC 100387T = LMG 21331T = CFN E1007T) were obtained from the Biological Resource Center, NITE (NBRC), Tokyo, Japan.

2.2. DNA Extraction, 16S rRNA Gene Sequencing and Phylogenetic Analysis

The DNA of the strains was extracted from bacterial biomass using Diatom™ DNA Prep 100 (Isogen Laboratory, Moscow, Russia) according to the manufacturer’s recommendations. The purified DNA samples were used as the PCR templates. Pure cultures of bacteria were identified by 16S rRNA gene analysis using the Bacteria-universal primer set 27F–1492R [31]. DNA sequencing was performed according to the Sanger method on an ABI 3730 DNA Analyzer automatic sequencer using the ABI PRISM® BigDye™ Terminator v. 3.1 reagent kit (Applied Biosystems, Waltham, MA, USA).
Initial taxonomic assignment of strains HO-A22T and SHC 2-14 was carried out using EzBioCloud [32]. The phylogenetic position of the isolated strains was determined using 16S rRNA gene sequences of type strains of the family Rhizobiaceae. The 16S rRNA gene sequences were aligned by MUSCLE [33], and the maximum-likelihood tree was performed with a GTR+F+I+G4 model recommended by ModelFinder [34] in IQ-TREE [35]. Branch supports were obtained with 10,000 ultrafast bootstraps [36]. Maximum parsimony and neighbor-joining trees were reconstructed with MPBoot [37] and MEGA11 [38], respectively.

2.3. Genome Sequencing and Analyses

Genomic DNA for the sequencing was extracted from bacterial biomass using the QIAamp DNA Mini Kit (QIAGEN, Germantown, MD, USA). The DNA libraries were constructed with the NEBNext DNA library prep reagent set for Illumina, according to the protocol for the kit. Sequencing of genomic DNA was carried out using the HiSeq 2500 platform (Illumina, Inc., San Diego, CA, USA) with 150-bp paired-end reads. The raw reads were quality checked with FastQC v. 0.11.9 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 8 January 2019), Version 0.12.0 released 1 March 2023), and low-quality reads were trimmed with Trimmomatic v. 0.39 [39], using the default settings for paired-end reads. The quality-filtered reads were then de novo assembled with SPAdes v. 3.13.0 [40] using the default settings. The resulting assembly was quality assessed with QUAST v. 5.0 [41]. The genome coverages were estimated by QualiMap 2 v. 2.2.2 [42] and Bowtie 2 v. 2.3.5.1 [43]. Genome annotation was performed using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP; v. 4.7) [44].
Complete genomic sequences were analyzed to confirm the assignment of new strains at the species and genus levels as recommended by Chun et al. [45]. Phylogenetic analysis of the isolated strains and members of the family Rhizobiaceae was performed based on a concatenated alignment of 120 single-copy genes obtained using GTDB-Tk software version 1.0.2 [46]. The maximum likelihood phylogenetic tree was calculated using IQ-TREE [35], based on the ModelFinder recommendations [34], and the branching support was estimated using UFBoot2 [36]. MPBoot [37] and MEGA11 [38] were used to reconstruct the maximum parsimony and neighbor-joining trees, respectively. DNA–DNA hybridization (dDDH) and average nucleotide identity (ANI) values were determined using Genome-to-Genome Distance Calculator (GGDC) v. 3.0 [47,48] and FastANI v. 1.3 [49], respectively. Average amino acid identity (AAI) was calculated using CompareM v. 0.0.23 with the aai_wf function. POCP values were calculated with the script runPOCP.sh [50] based on the previously described approach [51]. The pangenomic analysis was performed based on a bioinformatic pipeline proposed [52] using Anvi’o version 7.0 [53]. Using the MCL algorithm (Euclidean distance, Ward linkage), genomes were arranged according to the distribution of gene clusters.

2.4. Phenotypic Characterization

TEG medium was used to investigate the morphological, physiological, and biochemical characteristics of the strains, including fatty acids profiles, phospholipids of the cell walls, and quinones. All physiological tests were carried out in triplicate and chemotaxonomic tests in duplicate. Growth at different temperatures (5, 9, 15, 22, 28, 30, 33, 37, and 42 °C), salt tolerance (0, 0.1, 0.3, 0.5, 0.8, 1.2, 1.6, 2.0, 2.5, 3.0, 3.5, and 4.0% (w/v) NaCl) and pH range for growth (5.0–10.5 with increments of ~0.5 pH units) were determined on TEG medium for 3–5 days at aerobic conditions optimal for growth of the strains.
Motility of the studied cells was estimated under an Olympus CX41 system microscope (Tokyo, Japan) with a phase contrast device and a 100× objective with oil immersion. For the scanning electron microscopy, biomass of the isolated strains was grown in a TEG medium with Teflon cubes and then was prepared for the observation using phosphate buffer (pH 7.0), ethyl alcohol, and acetone as described previously [54]. The samples were examined under a scanning electron microscope (Quattro S, “Thermo Fisher Scientific Brno s.r.o.”, Brno-Černovice, Czech Republic) under high vacuum at an accelerating voltage 15 kV.
Catalase activity was determined by the conventional method with hydrogen peroxide. Enzyme activities and substrates consumption by the studied strains were measured using API® ZYM, API® 20E, and API® 50CH tests (bioMérieux SA, Marcy-l’Étoile, France) at 22 °C for 3–5 days. Antibiotics sensitivity was investigated on Petri dishes with susceptibility test discs (HiMedia, Mumbai, Maharashtra, India) with carbenicillin (100 μg), erythromycin (15 μg), gentamicin (10 μg), and kanamycin (30 μg), based on the presence of growth inhibition zones >10 mm around the discs after 5 days of cultivation.
To study utilization of oil alkanes, Adkins mineral medium [55] with pH 7.0–7.2 supplemented with a premade trace elements solution (10 mL·L–1) [56] was used. Sterilized oil from production well 5462 of the Cheremukhovskoe oil field was added at a concentration of 0.5% (v/v). To investigate the denitrifying activity of the strains, Adkins medium with argon as the gas phase was deoxygenated by boiling and supplemented with glucose (5.0 g·L–1), L-cysteine (0.5 g·L–1), and sodium nitrate (0.85 g·L–1).

2.5. Chemotaxonomic Characterization

The fatty acid composition was determined as the percentage of the total ion current peak area using a 7890B+5977B gas chromatograph-mass spectrometer (Agilent Technologies, Santa Clara, CA, USA). The biomass of the strains was dried with methanol and subjected to acidic methanolysis (1.2 M HCl/MeOH, 80 °C, 10 min) as described earlier [57]. The fatty acid composition was analyzed using a Maestro gas chromatograph-mass spectrometer (Interlab, Moscow region, Russia). The analysis of respiratory quinones of strains was performed at the All-Russian Collection of Microorganisms. Isoprenoid quinones were extracted from wet cells, purified according to Collins and Jones [58], and analyzed with an Thermo Finnigan LCQ Advantage MAX mass spectrometer (Thermo Fisher Scientific Inc., Waltham, MA, USA). Two-dimensional thin-layer chromatography on silica gel layers was performed for polar lipids identification [59]. Spraying with ninhydrin, α-naphthol, Dittmer–Lester molybdenum blue reagent, and Dragendorff reagent was used for visualization of aminolipids, glycolipids, phospholipids, and choline, respectively, as described in detail previously [60].

2.6. Analytical Methods

The growth of biomass was monitored by the changes in its optical density at 660 nm using an Ultrospec 2100 pro spectrophotometer (Amersham Biosciences, Amersham, Buckinghamshire, UK). The pH of the cultivation media was measured with a Seven Compact S220 pH meter (Metter Toledo, Greifensee, Switzerland). Formation of N2 or/and N2O by nitrate-reducing bacteria was estimated using a Crystal 5000 gas chromatograph (Chromatec, Yoshkar-Ola, Mari El, Russia) with a HayeSep-N 80/100 column.
Rheological characteristics of the cultivation media, i.e., surface tension (between the medium and air) and interfacial tension (medium/n-hexadecane) were measured with a Surface Tensiomat 21 semi-automatic tensiometer (Cole-Parmer, Vernon Hills, IL, USA) by the ring separation method at 25 °C. To estimate the utilization of oil alkanes, nonpolar compounds of the cultivation media were extracted with n-hexane and fractionated on silica gel columns. The profiles of light aliphatic hydrocarbons were determined using a Crystal 5000.1 gas chromatograph (Chromatec, Yoshkar-Ola, Mari El, Russia) with a ZB-FFAP 15 m capillary column and a flame ionization detector. Evaluation of oil biodegradation by the studied strains and its internal normalization based on the n-alkanes to iso-alkanes ratio were performed as described earlier [61].

2.7. Bacterial Reduction of Radionuclides and Metals

The biomass of the strain SHC 2-14 was grown aerobically on the TEG medium and collected by centrifugation [22]. About 0.125 g·L–1 (dry weight) of the live biomass was incubated for 10 days anaerobically in physiological saline using (metal ion concentrations, mg·L–1) uranyl nitrate (50), sodium pertechnetate (10), potassium bichromate (30), sodium selenate (50), ammonium vanadate (50), or/and sodium molybdate (30) as electron acceptors, and sodium acetate (20 mM) or sodium lactate (30 mM) as electron and carbon donors. The oxidation state of U, Cr, Se, V, and Mo was determined using a Bruker Elexsys E580 pulse EPR spectrometer (Billerica, MA, USA) in 0.6 mL of the liquid phase. Tc(IV) was determined by X-ray photoelectron spectroscopy (XPS) analysis of the solid phase. The samples were centrifuged for 5 min; then the supernatant was decanted, dried at 105 °C, and milled. Milled samples were embedded into indium on an aluminum support plate in form of thin films with mirror-like surfaces. Each sample was measured under a vacuum of 5·10–7 Pa at 24 °C using the low-energy gun for the sample charge compensation. The XPS spectra were obtained using a Kratos Axis Ultra DLD (Kratos Analytical Ltd., Manchester, UK) spectrometer with AlKα monochromatized (1486.6 eV) X-ray source.

2.8. Nucleotide Sequence Accession Numbers

The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA gene sequences of strains HO-A22T and SHC 2-14 are MT495799 and MG051317, respectively. The GenBank/EMBL/DDBJ accession numbers of the genomic assemblies of strains HO-A22T and SHC 2-14 are GCF_013371465.1 (JABWDU000000000.1) and GCF_013141885.1 (SROZ00000000.1), respectively.

3. Results and Discussion

3.1. Morphological Characterization of the Isolated Strains

To study the characteristics of the isolated strains, they were analyzed with morphological, genomic, physiological, and chemotaxonomic methods. Cells of the strains HO-A22T and SHC 2-14 were found to be asporogenous, motile, Gram-stain negative rods, 1.0–1.5 × 0.3–0.5 μm and 1.3–2.2 × 0.6–0.9 μm, respectively (Figure 1). After 5 days of cultivation at 22 °C, the strains formed creamy, smooth, shiny, and round colonies on the medium, up to 8.0 and 6.0 mm in diameter, respectively. The purity of the isolated cultures was confirmed by light microscopy and by 16S rRNA gene analysis.

3.2. Analysis of the 16S rRNA Genes

The 16S rRNA gene sequence similarity between new strains was 99.8%, indicating that they could probably belong to the same species [45]. The 16S rRNA genes of strains HO-A22T and SHC 2-14 showed the highest similarity to those of the members of the genera Ensifer (99.0–100.0%) and Sinorhizobium (97.4–98.8%) of the family Rhizobiaceae. On the phylogenetic tree (Figure 2), the 16S rRNA gene sequences of the strains HO-A22T and SHC 2-14 formed a branch with that of Ensifer morelensis Lc04T [4]. At the time of writing this article, a new species ‘Ensifer canadensis’ has been described, also belonging to the E. morelensis Lc04T branch [6]. The tree had low bootstrap values, which indicated the need for phylogenetic analysis. Nevertheless, the results of phylogenetic analysis of the 16S rRNA genes demonstrated that the new isolates belonged to the Ensifer genus. To clarify the taxonomic position of strains HO-A22T and SHC 2-14, their morphological, physiological, and chemotaxonomic properties were studied, and their genomes were analyzed.

3.3. Phylogenetic Analysis

The strain HO-A22T genome consists of 6,762,417 bp and comprises 32 scaffolds, an N50 value of 1,120,613 bp, a G+C content of 61.7%, and coverage of 225×. The strain SHC 2-14 genome consists of 8,169,716 bp and comprises 33 scaffolds, an N50 value of 332,742 bp, a G+C content of 61.0%, and coverage of 165×. Annotation led to the identification of 6343 genes, including 6114 protein-coding sequences, 174 pseudogenes, and 55 RNA genes for strain HO-A22T and 7808 genes, including 7411 protein-coding sequences, 340 pseudogenes, and 57 RNA genes for strain SHC 2-14. On the phylogenetic tree, strains HO-A22T and SHC 2-14 formed a branch with the type strains of E. morelensis, E. adhaerens, E. sesbaniae, and ‘E. canadensis’ (Figure 3). This branch was located at a considerable distance from that of the Sinorhizobium type species. The calculated AAI and POCP values between new strains and the analyzed Ensifer strains were 83.7–99.1% and 74.1–88.7%, respectively (Table 1). These values for the strain HO-A22T exceed the thresholds for identifying bacteria as a separate genus (AAI for the family Rhizobiaceae < 76.5% [5] and POCP < 50–60% [51,62,63,64]. These results confirm the data of phylogenetic analysis of the 16S rRNA gene sequences, indicating that the strains HO-A22T and SHC 2-14 belong to the Ensifer genus.
Pangenomic analysis was performed on 18 genomes from the Ensifer/Sinorhizobium clade. The resulting pangenome consisted of 117,350 genes grouped into 22,072 gene clusters (GCs). Among these clusters, 2502 were identified as core clusters present in all analyzed genomes, indicating their significance for the clade (Figure 4).
Core clusters included the genes involved in amino acid transport and metabolism (291 GCs) and translation, ribosomal structure, and biogenesis (211 GCs). The analysis revealed 207 GCs responsible for complete pathways of carbohydrate metabolism, encompassing glycolysis (Embden–Meyerhof pathway), Entner–Doudoroff pathway, pyruvate oxidation, citrate cycle (TCA cycle), pentose phosphate cycle, PRPP biosynthesis, UDP-N-acetyl-D-glucosamine biosynthesis, and the exo cluster for succinoglycan-type exopolysaccharide biosynthesis and production. Additionally, 101 GCs involved in energy metabolism were identified, encompassing genes of NADH:quinone oxidoreductase, succinate dehydrogenase, cytochrome bc1 complex, cytochrome c oxidase, and F-type ATPase. The genomes also contained the genes for assimilatory nitrate and sulfate reduction.
Comparing Ensifer and Sinorhizobium genomes, 313 GCs were found to be shared by all Ensifer genomes but not in Sinorhizobium genomes. These included curA (NADPH-dependent curcumin reductase), cysP (sulfate/thiosulfate transport system substrate-binding protein), gli (D-galactarolactone isomerase), uxaB (tagaturonate reductase), yphE and yphF (sugar transport system genes), adeS (sensor histidine kinase), lplA, lplB, and lplC (putative aldouronate transport system genes), as well as ceuA, ceuC, and ceuD (iron-siderophore transport system genes).
The genome of strain HO-A22T possessed 339 GCs unique to the Ensifer/Sinorhizobium clade, 75 of which were annotated and involved in transcription, signal transduction mechanisms, amino acid metabolism, and inorganic ion transport and metabolism. The pangenome of ‘E. canadensis’ SHC 2-14 and ‘E. canadensis’ T173T consisted of 15,524 genes grouped into 8004 GCs. These genomes harbored 451 unique for the ‘E. canadensis’ clusters, including the genes associated with heavy-metal resistance, such as those of arsenical pump membrane protein ArsB and P-type Cu2+ transporter CopB (EC: 7.2.2.9).
These findings demonstrate the genetic diversity within the Ensifer/Sinorhizobium clade and provide insights into the adaptations of the studied strains to their respective habitats.

3.4. Genomic Analysis and Environmental Implications

The genomes of strains HO-A22T and ‘E. canadensis’ SHC 2-14 both possess numerous carbon metabolism genes, including those for glycolysis (Embden–Meyerhof pathway), gluconeogenesis, TCA cycle, pentose phosphate and Entner–Doudoroff pathways, which are typical for the members of the Ensifer/Sinorhizobium clade [65,66,67,68,69]. The ‘E. canadensis’ SHC 2-14 possessed genes encoding 6-phospho-beta-glucosidase (EC: 3.2.1.86), fructose 1,6-bisphosphate 1-phosphatase (EC: 3.1.3.11), citrate (pro-3S)-lyase beta chain (EC: 4.1.3.6), gluconolactonase (EC: 3.1.1.17), and gluconate dehydratase (EC: 4.2.1.39) that the strain HO-A22T did not possessed.
The mechanisms of reduction of nitrogen compounds by members of the Ensifer/Sinorhizobium group have been most fully studied using the E. (S.) meliloti strain by Torres et al. [70,71]. Strains HO-A22T and ‘E. canadensis’ SHC 2-14 similarly possessed the napADEF, nirKV, and norDEQ genes, i.e., they were able to reduce nitrate to nitrous oxide due to the presence of ferredoxin-type periplasmic nitrate reductase (EC: 1.7.99.4), copper-containing nitrite reductase (NO-forming) (EC: 1.7.2.1), and nitric oxide reductase (EC: 1.7.99.7). At the same time, the genomes of both strains contained the nirD gene encoding nitrite reductase [NAD(P)H] (EC: 1.7.1.4). As a result, they could potentially carry out dissimilatory nitrite reduction to ammonium (DNRA pathway) as well. The ‘E. canadensis’ SHC 2-14 possesses the nosDFRZ genes encoding nitrous-oxide reductase (EC: 1.7.99.6); this corroborates the results of physiological experiments where the strain released molecular nitrogen into the gas phase during its growth [21]. It is important to highlight that the genome of strain SHC 2-14 is characterized by the absence of nifH, nifX, and nifE genes, which are associated with nitrogen fixation. These specific genes were previously identified in the genome of ‘E. canadensis’ T173T and appear to have been acquired through horizontal transfer from Sinorhizobium medicae [6]. The strain HO-A22T lacks these genes, so the terminal product of its denitrification was nitrous oxide. Nevertheless, this strain could potentially be used in biotechnologies to suppress the growth of sulfate-reducing bacteria (and therefore oil field equipment corrosion) by releasing the nitrite ion, an intermediate product of denitrification, if it accumulates in the formation water of an oil field [54].
Members of the Ensifer/Sinorhizobium clade are known for their ability to form extracellular polysaccharides. The eps and exo gene clusters have been identified as playing an essential role in exopolysaccharide production, which is vital for symbiotic rhizobial bacteria, since polysaccharides are directly involved in root nodule occupancy [72]. Bacteria of the genus Ensifer are typically nonsymbiotic species, but they are also capable of releasing extracellular polysaccharides. The epsF and exoFQZ genes encoding glycosyltransferase responsible for exopolysaccharide biosynthesis and the exopolysaccharide production protein were found in the genomes of strains HO-A22T and SHC 2-14, which were therefore able to form these compounds and increase the viscosity of the cultivation media. In addition, the studied strains possessed the genes rhaIMS and rfbCD, responsible for L-rhamnose isomerization (EC: 5.3.1.14) and mutarotation (EC: 5.1.3.32) [72] and dTDP-4-dehydrorhamnose reduction (EC: 1.1.1.133) and epimerization (EC: 5.1.3.13) [73,74], respectively. These processes are known to be the key stages in the synthesis of succinoglycan exopolysaccharides [9,75,76], and, therefore, it becomes possible to suggest that strains HO-A22T and SHC 2-14 produced this type of extracellular polysaccharide. The release of such compounds by strain HO-A22T into the formation water of an oil reservoir may be of practical interest for microbial enhanced oil recovery.
Members of the class Alphaproteobacteria are frequently found in oil-contaminated soils due to their resistance to high concentrations of nonpolar pollutants and the ability to use various oil components for their growth [77]. The genomes of the strains HO-A22T and SHC 2-14 contained a number of genes encoding the enzymes for aerobic degradation of long-chain fatty acids, aliphatic hydrocarbons, benzene, benzoate (BenB—benzoate 1,2-dioxygenase (EC: 1.14.12.10)), phenol, catechol (CatA—catechol 1,2-dioxygenase (EC: 1.13.11.1)), toluene, styrene, and phenanthrene. Various degradation pathways for phenanthrene and other PAHs have been described for the Ensifer/Sinorhizobium group [78,79,80]. The wide range of hydrocarbon compounds that strain HO-A22T could potentially consume attests to its adaptation to the environmental conditions in petroleum reservoirs and the potential of its application in MEOR biotechnologies by stimulating its growth directly in the oil field.
Bacterial strains isolated from radionuclide-contaminated groundwater were found to be tolerant to the high concentrations of heavy metals due to the presence of HMRGs (heavy-metal resistance genes) in their genomes [29,81,82]. The strains HO-A22T and SHC 2-14 were characterized by an abundance of HMRGs associated with copper, cadmium, zinc, lead, cobalt, mercuric, and chromium detoxification. Genomes of the strains contained the petE gene encoding plastocyanin that affected copper concentration and, therefore, provided tolerance to the ion [83]. The strains possessed the cutCE genes responsible for copper homeostasis as well, and that was their presumable alternate mechanism of copper resistance using cytoplasmic copper-binding protein CutE [84,85]. Strains HO-A22T and ‘E. canadensis’ SHC 2-14 contained genes of the CadA enzyme—copper-translocating P-type ATPase (EC: 3.6.3.4), which is also known for cadmium, zinc, and lead detoxification [86]. The copG gene in both genomes encoded a copper amine oxidase, which was responsible for copper resistance as well [87]. Both strains were potentially tolerant of the mercury ions due to two different mechanisms: MerA enzyme is a widespread mercuric ion reductase (EC: 1.16.1.1), which is highly expressed under mercury stress in a number of prokaryotic organisms [88,89,90]; and ZntA protein is a lead-, cadmium-, zinc-, and mercury-transporting ATPase, which is found in phylogenetically closely related bacteria of the Ensifer/Sinorhizobium clade as well [8,91]. Finally, tolerance of the isolated strains to the chromium ions was possible due to the presence of the chrAF genes responsible for chromate transport (ChrA is a chromate efflux transporter) in their genomes [92,93]. The abundance of HMRGs in genomes of the strains HO-A22T and ‘E. canadensis’ SHC 2-14 allowed them to inhabit environments highly polluted with heavy metals, such as radionuclide-contaminated groundwater.
Thus, the genomes of strains HO-A22T and ‘E. canadensis’ SHC 2-14 contained a number of gene cluster activity which was responsible for their experimentally observed physiological characteristics. The isolates possessed the genes for various pathways of carbon metabolism, denitrification and reduction of nitrate to ammonium (DNRA process), synthesis of exopolysaccharides, presumably extracellular succinoglycan-type polysaccharides, degradation of aliphatic and aromatic hydrocarbon components of oil, and detoxification of various heavy metals and/or radionuclides, which allowed the studied strains to grow in the environmental conditions atypical for members of the Ensifer/Sinorhizobium clade and confirmed that they were adapted to the natural habitats from which they were isolated.

3.5. Physiological Characterization

Physiology of the strains HO-A22T and ‘E. canadensis’ SHC 2-14 was investigated using TEG medium, typically at 24 °C because of the low temperature of their natural habitats. The strain HO-A22T grew at 0–3% (w/v) NaCl in the medium (optimum, 0–1.2% NaCl); and ‘E. canadensis’ SHC 2-14 tolerated 0–1.6% (w/v) NaCl with an optimal growth at 0–0.5% NaCl. The strain HO-A22T grew at 15–33 °C (optimum, 22–28 °C); and ‘E. canadensis’ SHC 2-14 tolerated 9–33 °C with an optimum at 28–30 °C. HO-A22T grew at pH 6.0–9.8 with an optimal growth at pH 6.6–8.0; and the strain SHC 2-14 tolerated pH 5.9–9.2 (optimum, pH 6.6–8.6). Growth profiles of strains HO-A22T and SHC 2-14 at various temperatures, NaCl concentrations, and pH are shown in Figure S1. Both strains were catalase-positive. Strain HO-A22T could grow anaerobically reducing nitrite to nitrous oxide (N2O) [94]; ‘E. canadensis’ SHC 2-14 could denitrify as well, but the terminal product of the process was molecular nitrogen (N2). Strains HO-A22T and SHC 2-14 were resistant to carbenicillin (100 μg), erythromycin (15 μg), and kanamycin (30 μg), but sensitive to gentamicin (10 μg).
In the API® ZYM tests with the strains HO-A22T, SHC 2-14, E. adhaerens AT, and E. morelensis Lc04T, all four studied strains were positive or weakly positive for alkaline phosphatase, esterase (C4), esterase lipase (C8), α-chymotrypsin, acid phosphatase, naphthol-AS-BI-phosphohydrolase, ß-galactosidase, α-glucosidase, and ß-glucosidase, but negative for lipase (C14), α-galactosidase, ß-glucuronidase, α-mannosidase, and α-fucosidase. Strains HO-A22T and SHC 2-14 were both positive for N-acetyl-ß-glucosaminidase (which only E. morelensis Lc04T was negative for) and negative for cystine arylamidase and trypsin; HO-A22T was positive for leucine arylamidase and weakly positive for valine arylamidase, whereas ‘E. canadensis’ SHC 2-14 was negative for both of them, unique among the studied strains. Results presented in Table S1 show the similarities in enzyme activities of the studied strains of the genus Ensifer. A number of diagnostic features were not given in the description of the species ‘E. canadensis’, and the results of the study of the strain SHC 2-14 complement the phenotypic properties of this species.
According to the API® 20E results, strains HO-A22T, ‘E. canadensis’ SHC 2-14, E. adhaerens AT, and E. morelensis Lc04T were positive or weakly positive for ß-galactosidase, tryptophane deaminase, rhamnose and arabinose fermentation/oxidation, and NO2 production, but negative for arginine dihydrolase, ornithine decarboxylase, citrate utilization, H2S production, indole production, acetoin production (Voges Proskauer), gelatinase, mannitol, inositol, sorbitol, sucrose, and amygdalin fermentation/oxidation. Strains HO-A22T and ‘E. canadensis’ SHC 2-14 were both negative for lysine decarboxylase and urease (only the strain E. adhaerens AT possessed both of them). HO-A22T was positive for glucose and weakly positive for melibiose fermentation/oxidation, whereas ‘E. canadensis’ SHC 2-14 was negative for them. However, significant physiological differences between the strains were not observed in these experiments (Table S2).
In the API® 50CH tests, all four studied strains showed positive or weakly positive reactions of glycerol, erythritol, D-arabinose, L-arabinose, D-ribose, D-xylose, L-xylose, D-adonitol, D-galactose, D-glucose, D-fructose, D-mannose, L-rhamnose, inositol, D-mannitol, D-sorbitol, aesculin (Fe citrate), D-cellobiose, D-maltose, D-sucrose, D-trehalose, D-turanose, D-lyxose, D-tagatose, D-fucose, L-fucose, D-arabite, and L-arabite utilization, but negative reactions of methyl-βD-xylopyranoside, inulin, starch, glycogen, gluconate K, and 5-ketogluconate K consumption. Both strains HO-A22T and ‘E. canadensis’ SHC 2-14 did not utilize methyl-ɑD-mannopyranoside (only E. morelensis Lc04T weakly consumed it) and D-melicitose (only E. adhaerens AT did that weakly). HO-A22T showed positive or weakly positive reactions of utilization of L-sorbose, dulcitol, N-acetylglucosamine, amygdalin, arbutin, salicin, and 2-ketogluconate K, whereas ‘E. canadensis’ SHC 2-14 showed negative reactions for these substrates. In our study, the ‘E. canadensis’ SHC 2-14 did not consume methyl-ɑD-glucopyranoside, D-lactose, D-melibiose, D-raffinose, xylitol, and gentiobiose; most of these features were indicated as weakly positive in the species description [6]. Results presented in Table 2 and Table S3 show common characteristics of the studied strains and several properties differentiating these strains from other species of the genus Ensifer.
The growth of the strains HO-A22T and ‘E. canadensis’ SHC 2-14 on crude oil was investigated. After 30 days of cultivation at 30 °C, a decrease in the surface and interfacial tension in the cultivation medium by 5–10 mN/m was registered, which indicated the production of surfactants by strains. A decrease in the proportion of short-chain n-alkanes in the structure of oil saturated hydrocarbons degraded by each strain was also demonstrated (Figures S2 and S3), which indicated the consumption of these compounds by the new isolates. That is, both strains HO-A22T and ‘E. canadensis’ SHC 2-14 were able to utilize oil components, which testifies their ability to produce oil-displacing metabolites promising for the enhancement of oil recovery.
The habitat from which ‘E. canadensis’ SHC 2-14 was isolated was characterized by a high content of radionuclides, metals and metalloids. In experiments with the salts of uranium, technetium, iron, chromium, vanadium, selenium, and molybdenum it was shown that the strain reduced more than 98% of V(V) and Se(VI) added to the media in the presence of sodium salts of acetic and lactic acids under anaerobic conditions (Figure S4); high reduction rates of Fe(III), Cr(VI), and U(VI) were noted: 87, 69, and 59% of the amount added into the medium, respectively. The content of Tc(VII) and Mo(VI) in the media decreased slightly. In the control experiment with cells killed by sterilization, there was no change in the concentration of the considered elements in the medium. Reduced radionuclides and heavy metals form extremely poorly soluble compounds and are not transported with groundwater, that is, the obtained results indicated the potential for the biogenic reduction of radionuclides, heavy metals and metalloids in groundwater in the zone of liquid radioactive waste disposal and the distribution of these elements in natural systems.

3.6. Chemotaxonomic Characterization

Cellular fatty acid profile of the strain HO-A22T included C18:1 (66.8%), C16:0 (11.1%), C14:0 3-OH (7.5%), C16:1 (4.8%), C17:1 (4.6%), C18:0 (3.3%), C18:0 3-OH (1.1%), C17:0 (0.4%), and others (0.4%). Cellular fatty acid profile of the ‘E. canadensis’ SHC 2-14 included C18:1 (76.1%), C16:0 (10.5%), C16:1 (7.2%), C14:0 3-OH (3.6%), C18:0 (2.4%), and C17:0 (0.2%). These profiles were compared with the those of the strains E. adhaerens AT and E. morelensis Lc04T (Table S4). C18:1 was the predominant compound in all studied profiles. Other dominant fatty acids were C16:0, C16:1, C18:0, and C14:0 3-OH in all studies strains, except the strain E. adhaerens AT containing the C16:0 3-OH fatty acid instead of the C14:0 3-OH. The fatty acid profile of the strain HO-A22T was close to the one of E. morelensis Lc04T; the strain ‘E. canadensis’ SHC 2-14 lacked the widespread C18:0 3-OH and C17:1 components. Although the strain ‘E. canadensis’ T173T [6] was grown in a TY medium different from the TEG medium in which the four studied strains were grown, the major fatty acids C18:1, C16:0, C18:0, and C16:1 were common to the strain T173T and reference strains of Ensifer species.
The major respiratory quinone for the strains HO-A22T and SHC 2-14 was Q10. The strain HO-A22T did not contain any traces of other menaquinones, whereas the strain SHC 2-14 possessed the minor quinone Q9 (Figure S5). Two studied type strains of the genus Ensifer, as well as the strains HO-A22T and SHC 2-14, contained phosphatidylcholines, phosphatidylethanolamines, phosphatidylglycerols, and diphosphatidylglycerols as major polar lipids of their cell membranes (Figures S6–S10). The ‘E. canadensis’ SHC 2-14 contained the unique polar lipid for this phylogenetic group, labeled as APL4, unidentified aminophospholipid.
Thus, four studied strains of the genus Ensifer had similar chemotaxonomic characteristics. All validly published members of the genus Ensifer were Gram-negative, motile, aerobic, non-spore-forming rods. They could reduce nitrate to nitrite, grow at 1% NaCl, 28 °C, pH 7.0–9.0, produce acid (API® 50CH) from D-arabinose, D-fructose, D-glucose, D-maltose, D-mannose, D-ribose, D-trehalose, D-xylose, and erythritol; none of them grew at 4% NaCl, produced neither H2S nor acid from inulin or starch. Major fatty acids of the type strains were C18:1, C16:0, C16:1, and C18:0. Major polar lipids were phosphatidylcholines, phosphatidylethanolamines, and phosphatidylglycerols. Significant biochemical differences between the strain HO-A22T and validly described species of the genus Ensifer were not observed; and ‘E. canadensis’ SHC 2-14 had slight peculiarities in its cellular fatty acid profile, respiratory quinones, and membranes polar lipids. These results confirmed that these properties of microorganisms have low taxonomic value.

4. Conclusions

Two novel strains, HO-A22T and SHC 2-14, of the genus Ensifer isolated from water injected to the oil field and from radionuclide-contaminated groundwater, were analyzed using polyphasic taxonomy approaches and genomic analysis and compared with respective features of members of the genus Ensifer. Based on the phylogenetic, phenotypic, and chemotaxonomic analyses, it is concluded that strain HO-A22T should be considered as a novel species of the genus Ensifer, for which the name Ensifer oleiphilus sp. nov. is proposed. The protologue of the proposed new species is presented in Table 3. The influence of environmental conditions such as salinity, pH, temperature, the content of crude oil or nitrates, heavy metals and radionuclides on the growth of strains revealed their adaptation to the environmental conditions. The strains HO-A22T and SHC 2-14 can potentially be used in biotechnologies of microbial enhanced oil recovery, suppression of sulfide accumulation in oil reservoirs, and radionuclide pollution treatment due to their abilities to utilize hydrocarbons of crude oil, to produce extracellular polysaccharides, and to reduce nitrates and heavy metals.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11092314/s1, Figure S1: Growth profiles of strains HO-A22T (a–c) and ‘E. canadensis’ SHC 2-14 (d–f) in the TEG medium at various temperatures (a,d), NaCl concentrations (w/v, %) (b,e), and pH (c,f); Figure S2: Chromatograms of the saturated hydrocarbon fraction of sterile oil (a) and oil degraded by the strain HO-A22T (b) and residual content of alkanes (%) in oil degraded by the strain HO-A22T (c). The strain was incubated in a liquid medium with crude oil as a substrate at 30 °C for 30 days; Figure S3: Chromatograms of the saturated hydrocarbon fraction of sterile oil (a) and oil degraded by the strain ‘E. canadensis’ SHC 2-14 (b) and residual content of alkanes (%) in oil degraded by the strain SHC 2-14 (c). The strain was incubated in a liquid medium with crude oil as a substrate at 30 °C for 30 days; Figure S4: Residual content (%) of Se(VI), V(V), Fe(III), Cr(VI), U(VI), Tc(VII), and Mo(VI) in the medium after incubation of the strain ‘E. canadensis’ SHC 2-14 at 22 °C for 10 days. The strain grew in the Adkins mineral medium with sodium acetate (20 mM) or sodium lactate (30 mM); Figure S5: Mass spectra of menaquinones from the strains HO-A22T (a,b) and ‘E. canadensis’ SHC 2-14 (c,d), showing the presence of menaquinone Q10; Figure S6: Two-dimensional thin layer chromatogram of polar lipids from strains HO-A22T (a), ‘E. canadensis’ SHC 2-14 (b), E. adhaerens AT (c), and E. morelensis Lc04T (d). The components were visualized by staining with 5% sulfuric acid in ethanol and heating at 180 °C for 15 min. Abbreviations: PC, phosphatidylcholines; PE, phosphatidylethanolamines; DPG, diphosphatidylglycerols; PG, phosphatidylglycerols; PME, phosphatidylmethylethanolamines; PS, phosphatidylserines; APL1–APL5, aminophospholipids; PL1–PL5, phospholipids; GL1–GL3, glycolipids; GPL1, glycophospholipid; AL1, aminolipid; PA, phosphatidic acids; LPE, lysophosphatidylethanolamines; X, unidentified lipid; Figure S7: Identification of polar lipids from the strain HO-A22T. The components were visualized by molybdenum blue (a); α-naphthol (b); ninhydrin (c); and Dragendorff (d). Abbreviations: AL, aminolipids; GL, glycolipids; PC, phosphatidylcholines; PE, phosphatidylethanolamines; PL, phospholipids; Figure S8: Identification of polar lipids from the strain ‘E. canadensis’ SHC 2-14. The components were visualized by molybdenum blue (a); α-naphthol (b); ninhydrin (c); and Dragendorff (d). Abbreviations: AL, aminolipids; PC, phosphatidylcholines; PL, phospholipids; Figure S9: Identification of polar lipids from the strain E. adhaerens AT. The components were visualized by molybdenum blue (a); α-naphthol (b); ninhydrin (c); and Dragendorff (d). Abbreviations: AL, aminolipids; DPG, diphosphatidylglycerols; GL, glycolipids; LPE, lysophosphatidylethanolamines; PA, phosphatidic acids; PC, phosphatidylcholines; PE, phosphatidylethanolamines; PG, phosphatidylglycerols; PL, phospholipids; PS, phosphatidylserines; Figure S10: Identification of polar lipids from the strain E. morelensis Lc04T. The components were visualized by molybdenum blue (a); α-naphthol (b); ninhydrin (c); and Dragendorff (d). Abbreviations: AL, aminolipids; DPG, diphosphatidylglycerols; LPE, lysophosphatidylethanolamines; PA, phosphatidic acids; PC, phosphatidylcholines; PE, phosphatidylethanolamines; PG, phosphatidylglycerols; PL, phospholipids; PS, phosphatidylserines; Table S1: Enzymatic activities of the strains HO-A22T, ‘E. canadensis’ SHC 2-14, E. adhaerens AT, and E. morelensis Lc04T determined by the API® ZYM test (bioMérieux, France). Designations: +, positive; –, negative; w, weakly positive; Table S2: Comparison of enzymatic activities of the strains HO-A22T, ‘E. canadensis’ SHC 2-14, E. adhaerens AT, and E. morelensis Lc04T, determined by the API® 20E test (bioMérieux, France). Designations: +, positive; –, negative; w, weakly positive; Table S3: Carbohydrate fermentation by the strains HO-A22T, ‘E. canadensis’ SHC 2-14, E. adhaerens AT, and E. morelensis Lc04T, determined by the API® 50CH test (bioMérieux, France). Designations: +, positive; –, negative; w, weakly positive; Table S4: Cellular fatty acid profiles of the strains HO-A22T, ‘E. canadensis’ SHC 2-14, E. adhaerens AT, and E. morelensis Lc04T.

Author Contributions

A.P.E.: Investigation, Data Curation, Writing—Original draft preparation. T.L.B.: Investigation, Validation; Visualization. D.S.G.: Software, Formal analysis; Data Curation, Visualization, Writing—Original draft preparation. D.S.S.: Investigation, Validation. E.M.S.: Investigation, Validation. A.N.A.: Investigation, Validation. A.B.P.: Investigation, Validation. E.A.I.: Investigation, Validation. T.N.N.: Conceptualization; Supervision, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation (grant 21-64-00019); work of TB, DS, and ES was funded by the Ministry of Science and Higher Education of the Russian Federation (grant 122040800164-6).

Data Availability Statement

The 16S rRNA gene sequences of strains HO-A22T and SHC 2-14 have been deposited at the GenBank/EMBL/DDBJ under accession numbers MT495799 and MG051317, respectively. The genomic assemblies of strains HO-A22T and SHC 2-14 have been deposited at the Gen-Bank/EMBL/DDBJ under accession numbers GCF_013371465.1 (JABWDU000000000.1) and GCF_013141885.1 (SROZ00000000.1), respectively.

Acknowledgments

The authors are grateful to N.G. Loiko from the Research Center of Biotechnology of the Russian Academy of Sciences for micrographs of bacteria. SEM studies were carried out at the Shared Research Facility “Electron microscopy in life sciences” at the Moscow State University (Unique Equipment “Three-dimensional electron microscopy and spectroscopy”) and using the equipment purchased under the Moscow State University Development Program.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Willems, A.; Fernández-López, M.; Muñoz-Adelantado, E.; Goris, J.; De Vos, P.; Martínez-Romero, E.; Toro, N.; Gillis, M. Description of new Ensifer strains from nodules and proposal to transfer Ensifer adhaerens Casida 1982 to Sinorhizobium as Sinorhizobium adhaerens comb. nov. Request for an Opinion. Int. J. Syst. Evol. Microbiol. 2003, 53, 1207–1217. [Google Scholar] [CrossRef]
  2. Young, J.M. The genus name Ensifer Casida 1982 takes priority over Sinorhizobium Chen et al. 1988, and Sinorhizobium morelense Wang et al. 2002 is a later synonym of Ensifer adhaerens Casida 1982. Is the combination ‘Sinorhizobium adhaerens’ (Casida 1982) Willems et al. 2003 legitimate? Request for an Opinion. Int. J. Syst. Evol. Microbiol. 2003, 53, 2107–2110. [Google Scholar] [CrossRef]
  3. Young, J.M. Sinorhizobium versus Ensifer: May a taxonomy subcommittee of the ICSP contradict the Judicial Commission? Int. J. Syst. Evol. Microbiol. 2010, 60, 1711–1713. [Google Scholar] [CrossRef]
  4. Wang, Y.C.; Wang, F.; Hou, B.C.; Wang, E.T.; Chen, W.F.; Sui, X.H.; Chen, W.X.; Li, Y.; Zhang, Y.B. Proposal of Ensifer psoraleae sp. nov.; Ensifer sesbaniae sp. nov., Ensifer morelense comb. nov. and Ensifer americanum comb. nov. Syst. Appl. Microbiol. 2013, 36, 467–473. [Google Scholar] [CrossRef]
  5. Kuzmanović, N.; Fagorzi, C.; Mengoni, A.; Lassalle, F.; diCenzo, G.C. Taxonomy of Rhizobiaceae revisited: Proposal of a new framework for genus delimitation. Int. J. Syst. Evol. Microbiol. 2022, 72, 005243. [Google Scholar] [CrossRef]
  6. Bromfield, E.S.P.; Cloutier, S.; Hynes, M.F. Ensifer canadensis sp. nov. strain T173T isolated from Melilotus albus (sweet clover) in Canada possesses recombinant plasmid pT173b harbouring symbiosis and type IV secretion system genes apparently acquired from Ensifer medicae. Front. Microbiol. 2023, 14, 1195755. [Google Scholar] [CrossRef] [PubMed]
  7. Casida, L.E., Jr. Ensifer adhaerens gen. nov., sp. nov.: A bacterial predator of bacteria in soil. Int. J. Syst. Bacteriol. 1982, 32, 339–345. [Google Scholar] [CrossRef]
  8. DiCenzo, G.C.; Debiec, K.; Krzysztoforski, J.; Uhrynowski, W.; Mengoni, A.; Fagorzi, C.; Gorecki, A.; Dziewit, L.; Bajda, T.; Rzepa, G.; et al. Genomic and biotechnological characterization of the heavy-metal resistant, arsenic-oxidizing bacterium Ensifer sp. M14. Genes 2018, 9, 379. [Google Scholar] [CrossRef]
  9. Ma, L.; Chen, S.; Yuan, J.; Yang, P.; Liu, Y.; Stewart, K. Rapid biodegradation of atrazine by Ensifer sp. strain and its degradation genes. Int. Biodeterior. Biodegrad. 2017, 116, 133–140. [Google Scholar] [CrossRef]
  10. Zhao, F.; Guo, C.; Cui, Q.; Hao, Q.; Xiu, J.; Han, S.; Zhang, Y. Exopolysaccharide production by an indigenous isolate Pseudomonas stutzeri XP1 and its application potential in enhanced oil recovery. Carbohydr. Polym. 2018, 199, 375–381. [Google Scholar] [CrossRef]
  11. Alvarez, V.M.; Jurelevicius, D.; Serrato, R.V.; Barreto-Bergter, E.; Seldin, L. Chemical characterization and potential application of exopolysaccharides produced by Ensifer adhaerens JHT2 as a bioemulsifier of edible oils. Int. J. Biol. Macromol. 2018, 114, 18–25. [Google Scholar] [CrossRef]
  12. Lehman, A.P.; Long, S.R. Exopolysaccharides from Sinorhizobium meliloti can protect against H2O2-dependent damage. J. Bacteriol. 2013, 195, 5362–5369. [Google Scholar] [CrossRef] [PubMed]
  13. Quijada, N.M.; Abreu, I.; Reguera, M.; Bonilla, I.; Bolaños, L. Boron nutrition affects growth, adaptation to stressful environments, and exopolysaccharide synthesis of Ensifer meliloti. Rhizosphere 2022, 22, 100534. [Google Scholar] [CrossRef]
  14. Xu, L.; Chen, X.; Li, H.; Hu, F.; Liang, M. Characterization of the biosorption and biodegradation properties of Ensifer adhaerens: A potential agent to remove polychlorinated biphenyls from contaminated water. J. Hazard. Mater. 2016, 302, 314–322. [Google Scholar] [CrossRef] [PubMed]
  15. Agrawal, A.; An, D.; Cavallaro, A.; Voordouw, G. Souring in low-temperature surface facilities of two high-temperature Argentinian oil fields. Appl. Microbiol. Biotechnol. 2014, 98, 8017–8029. [Google Scholar] [CrossRef] [PubMed]
  16. An, B.A.; Shen, Y.; Voordouw, G. Control of sulfide production in high salinity Bakken Shale oil reservoirs by halophilic bacteria reducing nitrate to nitrite. Front. Microbiol. 2017, 8, 1164. [Google Scholar] [CrossRef]
  17. Sokolova, D.S.; Semenova, E.M.; Grouzdev, D.S.; Ershov, A.P.; Bidzhieva, S.K.; Ivanova, A.E.; Babich, T.L.; Sissenbayeva, M.R.; Bisenova, M.A.; Nazina, T.N. Microbial diversity and potential sulfide producers in the Karazhanbas oilfield (Kazakhstan). Microbiology 2020, 89, 459–469. [Google Scholar] [CrossRef]
  18. Torres, M.J.; Avila, S.; Bedmar, E.J.; Delgado, M.J. Overexpression of the periplasmic nitrate reductase supports anaerobic growth by Ensifer meliloti. FEMS Microbiol. Lett. 2018, 365, fny041. [Google Scholar] [CrossRef]
  19. Liu, L.X.; Li, Q.Q.; Zhang, Y.Z.; Hu, Y.; Jiao, J.; Guo, H.J.; Zhang, X.X.; Zhang, B.; Chen, W.X.; Tian, C.F. The nitrate-reduction gene cluster components exert lineage-dependent contributions to optimization of Sinorhizobium symbiosis with soybeans. Environ. Microbiol. 2017, 19, 4926–4938. [Google Scholar] [CrossRef]
  20. Ruiz, B.; Le Scornet, A.; Sauviac, L.; Rémy, A.; Bruand, C.; Meilhoc, E. The nitrate assimilatory pathway in Sinorhizobium meliloti: Contribution to NO production. Front. Microbiol. 2019, 10, 1526. [Google Scholar] [CrossRef]
  21. Safonov, A.V.; Babich, T.L.; Sokolova, D.S.; Grouzdev, D.S.; Tourova, T.P.; Poltaraus, A.B.; Zakharova, E.V.; Merkel, A.Y.; Novikov, A.P.; Nazina, T.N. Microbial community and in situ bioremediation of groundwater by nitrate removal in the zone of a radioactive waste surface repository. Front. Microbiol. 2018, 9, 1985. [Google Scholar] [CrossRef]
  22. Nazina, T.N.; Babich, T.L.; Kostryukova, N.K.; Sokolova, D.S.; Abdullin, R.R.; Tourova, T.P.; Poltaraus, A.B.; Kalmykov, S.N.; Zakharova, E.V.; Myasoedov, B.F.; et al. Microbial diversity and possible activity in nitrate- and radionuclide-contaminated groundwater. In Behavior of Radionuclides in the Environment I. Function of Particles in Aquatic System; Kato, K., Konoplev, A., Kalmykov, S.N., Eds.; Springer Nature Singapore Pte. Ltd.: Singapore, 2020; pp. 35–66. [Google Scholar] [CrossRef]
  23. Safonov, A.V.; Andryushchenko, N.D.; Ivanov, P.V.; Boldyrev, K.A.; Babich, T.L.; German, K.E.; Zakharova, E.V. Biogenic factors of radionuclide immobilization on sandy rocks of upper aquifers. Radiochemistry 2019, 61, 99–108. [Google Scholar] [CrossRef]
  24. Gadd, G.M. Metals, minerals and microbes: Geomicrobiology and bioremediation. Microbiology 2010, 156, 609–643. [Google Scholar] [CrossRef] [PubMed]
  25. Nazina, T.N.; Luk’yanova, E.A.; Zakharova, E.V.; Konstantinova, L.I.; Kalmykov, S.N.; Poltaraus, A.B.; Zubkov, A.A. Microorganisms in a disposal site for liquid radioactive wastes and their influence on radionuclides. Geomicrobiol. J. 2010, 27, 473–486. [Google Scholar] [CrossRef]
  26. Shukla, A.; Parmar, P.; Saraf, M. Radiation, radionuclides and bacteria: An in-perspective review. J. Environ. Radioact. 2017, 180, 27–35. [Google Scholar] [CrossRef] [PubMed]
  27. Lovley, D.R.; Phillips, E.J.P.; Gorby, Y.A.; Landa, E.R. Microbial reduction of uranium. Nature 1991, 350, 413–416. [Google Scholar] [CrossRef]
  28. Lloyd, J.R. Microbial reduction of metals and radionuclides. FEMS Microbiol. Rev. 2003, 27, 411–425. [Google Scholar] [CrossRef]
  29. Nazina, T.; Babich, T.; Kostryukova, N.; Sokolova, D.; Abdullin, R.; Tourova, T.; Kadnikov, V.; Mardanov, A.; Ravin, N.; Grouzdev, D.; et al. Ultramicrobacteria from nitrate- and radionuclide-contaminated groundwater. Sustainability 2020, 12, 1239. [Google Scholar] [CrossRef]
  30. Nazina, T.N.; Sokolova, D.S.; Babich, T.L.; Semenova, E.M.; Ershov, A.P.; Bidzhieva, S.K.; Borzenkov, I.A.; Poltaraus, A.B.; Khisametdinov, M.R.; Tourova, T.P. Microorganisms of low-temperature heavy oil reservoirs (Russia) and their possible application for enhanced oil recovery. Microbiology 2017, 86, 773–785. [Google Scholar] [CrossRef]
  31. Lane, D.J. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematics; Stackebrandt, E., Goodfellow, M., Eds.; John Wiley & Sons: New York, NY, USA, 1991; pp. 115–175. [Google Scholar]
  32. Yoon, S.-H.; Ha, S.-M.; Kwon, S.; Lim, J.; Kim, Y.; Seo, H.; Chun, J. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 2017, 67, 1613–1617. [Google Scholar] [CrossRef]
  33. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  34. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef]
  35. Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  36. Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the ultrafast bootstrap approximation. Mol. Biol. Evol. 2018, 35, 518–522. [Google Scholar] [CrossRef] [PubMed]
  37. Hoang, D.T.; Vinh, L.S.; Flouri, T.; Stamatakis, A.; von Haeseler, A.; Minh, B.Q. MPBoot: Fast phylogenetic maximum parsimony tree inference and bootstrap approximation. BMC Evol. Biol. 2018, 18, 11. [Google Scholar] [CrossRef]
  38. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  39. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  40. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  41. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  42. Okonechnikov, K.; Conesa, A.; García-Alcalde, F. Qualimap 2: Advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics 2016, 32, 292–294. [Google Scholar] [CrossRef]
  43. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  44. 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. [Google Scholar] [CrossRef]
  45. Chun, J.; Oren, A.; Ventosa, A.; Christensen, H.; Arahal, D.R.; da Costa, M.S.; Rooney, A.P.; Yi, H.; Xu, X.W.; De Meyer, S.; et al. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int. J. Syst. Evol. Microbiol. 2018, 68, 461–466. [Google Scholar] [CrossRef] [PubMed]
  46. Chaumeil, P.; Mussig, A.J.; Parks, D.H.; Hugenholtz, P. GTDB-Tk: A toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 2020, 36, 1925–1927. [Google Scholar] [CrossRef] [PubMed]
  47. 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 Bioinf. 2013, 14, 60. [Google Scholar] [CrossRef]
  48. Meier-Kolthoff, J.P.; Carbasse, J.S.; Peinado-Olarte, R.L.; Göker, M. TYGS and LPSN: A database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Res. 2022, 50, D801–D807. [Google Scholar] [CrossRef] [PubMed]
  49. Jain, C.; Rodriguez-R, L.M.; Phillippy, A.M.; Konstantinidis, K.T.; Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 2018, 9, 5114. [Google Scholar] [CrossRef]
  50. Grouzdev, D.S.; Rysina, M.S.; Bryantseva, I.A.; Gorlenko, V.M.; Gaisin, V.A. Draft genome sequences of ‘Candidatus Chloroploca asiatica’ and ‘Candidatus Viridilinea mediisalina’, candidate representatives of the Chloroflexales order: Phylogenetic and taxonomic implications. Stand. Genomic. Sci. 2018, 13, 24. [Google Scholar] [CrossRef]
  51. Qin, Q.-L.; Xie, B.-B.; Zhang, X.-Y.; Chen, X.-L.; Zhou, B.-C.; Zhou, J.; Oren, A.; Zhang, Y.-Z. A proposed genus boundary for the Prokaryotes based on genomic insights. J. Bacteriol. 2014, 196, 2210–2215. [Google Scholar] [CrossRef]
  52. Delmont, T.O.; Eren, E.M. Linking pangenomes and metagenomes: The Prochlorococcus metapangenome. PeerJ 2018, 6, e4320. [Google Scholar] [CrossRef]
  53. Eren, A.M.; Esen, O.C.; Quince, C.; Vineis, J.H.; Morrison, H.G.; Sogin, M.L.; Delmont, T.O. Anvi’o: An advanced analysis and visualization platform for ’omics data. PeerJ 2015, 3, e1319. [Google Scholar] [CrossRef]
  54. Sokolova, D.S.; Semenova, E.M.; Grouzdev, D.S.; Bidzhieva, S.K.; Babich, T.L.; Loiko, N.G.; Ershov, A.P.; Kadnikov, V.V.; Beletsky, A.V.; Mardanov, A.V.; et al. Sulfidogenic microbial communities of the Uzen high-temperature oil field in Kazakhstan. Microorganisms 2021, 9, 1818. [Google Scholar] [CrossRef] [PubMed]
  55. Adkins, J.P.; Cornell, L.A.; Tanner, R.S. Microbial composition of carbonate petroleum reservoir fluids. Geomicrobiol. J. 1992, 10, 87–97. [Google Scholar] [CrossRef]
  56. Pfennig, N.; Lippert, K.D. Über das Vitamin B12-Bedürfnis phototropher Schwefelbakterien. Arch. Mikrobiol. 1966, 55, 245–256. [Google Scholar] [CrossRef]
  57. Bidzhieva, S.K.; Sokolova, D.S.; Grouzdev, D.S.; Kostrikina, N.A.; Poltaraus, A.B.; Tourova, T.P.; Shcherbakova, V.A.; Troshina, O.Y.; Nazina, T.N. Sphaerochaeta halotolerans sp. nov., a novel spherical halotolerant spirochete from a Russian heavy oil reservoir, emended description of the genus Sphaerochaeta, reclassification of Sphaerochaeta coccoides to a new genus Parasphaerochaeta gen. nov. as Parasphaerochaeta coccoides comb. nov. and proposal of Sphaerochaetaceae fam. nov. Int. J. Syst. Evol. Microbiol. 2020, 70, 4748–4759. [Google Scholar] [CrossRef]
  58. Collins, M.D.; Jones, D. Distribution of isoprenoid quinone structural types in bacteria and their taxonomic implication. Microbiol. Rev. 1981, 45, 316–354. [Google Scholar] [CrossRef] [PubMed]
  59. Minnikin, D.E.; O’Donnell, A.G.; Goodfellow, M.; Alderson, G.; Athalye, M.; Schaal, A.; Parlett, J.H. An integrated procedure for the extraction of bacterial isoprenoid quinones and polar lipids. J. Microbiol. Methods 1984, 2, 233–241. [Google Scholar] [CrossRef]
  60. Semenova, E.M.; Grouzdev, D.S.; Sokolova, D.S.; Tourova, T.P.; Poltaraus, A.B.; Potekhina, N.V.; Shishina, P.N.; Bolshakova, M.A.; Avtukh, A.N.; Ianutsevich, E.A.; et al. Physiological and genomic characterization of Actinotalea subterranea sp. nov. from oil-degrading methanogenic enrichment and reclassification of the family Actinotaleaceae. Microorganisms 2022, 10, 378. [Google Scholar] [CrossRef]
  61. Semenova, E.M.; Babich, T.L.; Sokolova, D.S.; Ershov, A.P.; Raievska, Y.I.; Bidzhieva, S.K.; Stepanov, A.L.; Korneykova, M.V.; Myazin, V.A.; Nazina, T.N. Microbial communities of seawater and coastal soil of Russian Arctic region and their potential for bioremediation from hydrocarbon pollutants. Microorganisms 2022, 10, 1490. [Google Scholar] [CrossRef]
  62. Koziaeva, V.; Dziuba, M.; Leão, P.; Uzun, M.; Krutkina, M.; Grouzdev, D. Genome-based metabolic reconstruction of a novel uncultivated freshwater magnetotactic coccus “Ca. Magnetaquicoccus inordinatus” UR-1, and proposal of a candidate family “Ca. Magnetaquicoccaceae”. Front. Microbiol. 2019, 10, 2290. [Google Scholar] [CrossRef]
  63. Li, Y.; Xue, H.; Sang, S.; Lin, C.; Wang, X. Phylogenetic analysis of family Neisseriaceae based on genome sequences and description of Populibacter corticis gen. nov., sp. nov., a member of the family Neisseriaceae, isolated from symptomatic bark of Populus × euramericana canker. PLoS ONE 2017, 12, e0174506. [Google Scholar] [CrossRef] [PubMed]
  64. Orata, F.D.; Meier-Kolthoff, J.P.; Sauvageau, D.; Stein, L.Y. Phylogenomic analysis of the gammaproteobacterial methanotrophs (order Methylococcales) calls for the reclassification of members at the genus and species levels. Front. Microbiol. 2018, 9, 3162. [Google Scholar] [CrossRef] [PubMed]
  65. diCenzo, G.C.; Muhammed, Z.; Østerås, M.; O’Brien, S.A.; Finan, T.M. A key regulator of the glycolytic and gluconeogenic central metabolic pathways in Sinorhizobium meliloti. Genetics 2017, 207, 961–974. [Google Scholar] [CrossRef]
  66. Draghi, W.O.; Del Papa, M.F.; Barsch, A.; Albicoro, F.J.; Lozano, M.J.; Pühler, A.; Niehaus, K.; Lagares, A. A metabolomic approach to characterize the acid-tolerance response in Sinorhizobium meliloti. Metabolomics 2017, 13, 71. [Google Scholar] [CrossRef]
  67. Geddes, B.A.; Oresnik, I.J. Physiology, genetics, and biochemistry of carbon metabolism in the alphaproteobacterium Sinorhizobium meliloti. Can. J. Microbiol. 2014, 60, 491–507. [Google Scholar] [CrossRef]
  68. Hawkins, J.P.; Ordonez, P.A.; Oresnik, I.J. Characterization of mutations that affect the nonoxidative pentose phosphate pathway in Sinorhizobium meliloti. J. Bacteriol. 2018, 200, e00436-17. [Google Scholar] [CrossRef]
  69. Yurgel, S.N.; Qu, Y.; Rice, J.T.; Ajeethan, N.; Zink, E.M.; Brown, J.M.; Purvine, S.; Lipton, M.S.; Kahn, M.L. Specialization in a nitrogen-fixing symbiosis: Proteome differences between Sinorhizobium medicae bacteria and bacteroids. Mol. Plant-Microbe Interact. 2021, 34, 1409–1422. [Google Scholar] [CrossRef]
  70. Torres, M.J.; Rubia, M.I.; Bedmar, E.J.; Delgado, M.J. Denitrification in Sinorhizobium meliloti. Biochem. Soc. Trans. 2011, 39, 1886–1889. [Google Scholar] [CrossRef]
  71. Torres, M.J.; Rubia, M.I.; Coba de la Peña, T.; Pueyo, J.J.; Bedmar, E.J.; Delgado, M.J. Genetic basis for denitrification in Ensifer meliloti. BMC Microbiol. 2014, 14, 142. [Google Scholar] [CrossRef]
  72. Rivers, D.M.; Kim, D.D.; Oresnik, I.J. Inability to catabolize rhamnose by Sinorhizobium meliloti Rm1021 affects competition for nodule occupancy. Microorganisms 2022, 10, 732. [Google Scholar] [CrossRef]
  73. Zivkovic, M.; Miljkovic, M.; Ruas-Madiedo, P.; Strahinic, I.; Tolinacki, M.; Golic, N.; Kojic, M. Exopolysaccharide production and ropy phenotype are determined by two gene clusters in putative probiotic strain Lactobacillus paraplantarum BGCG11. Appl. Environ. Microbiol. 2015, 81, 1387–1396. [Google Scholar] [CrossRef] [PubMed]
  74. Gao, M.; D’Haeze, W.; De Rycke, R.; Wolucka, B.; Holsters, M. Knockout of an azorhizobial dTDP-L-rhamnose synthase affects lipopolysaccharide and extracellular polysaccharide production and disables symbiosis with Sesbania rostrata. Mol. Plant-Microbe Interact. 2001, 14, 857–866. [Google Scholar] [CrossRef]
  75. Jofré, E.; Lagares, A.; Mori, G. Disruption of dTDP-rhamnose biosynthesis modifies lipopolysaccharide core, exopolysaccharide production, and root colonization in Azospirillum brasilense. FEMS Microbiol. Lett. 2004, 231, 267–275. [Google Scholar] [CrossRef] [PubMed]
  76. Zhan, H.J.; Leigh, J.A. Two genes that regulate exopolysaccharide production in Rhizobium meliloti. J. Bacteriol. 1990, 172, 5254–5259. [Google Scholar] [CrossRef]
  77. Pourbabaee, A.A.; Khazaei, M.; Alikhani, H.A.; Emami, S. Root nodulation of alfalfa by Ensifer meliloti in petroleum contaminated soil. Rhizosphere 2021, 17, 100305. [Google Scholar] [CrossRef]
  78. Keum, Y.S.; Seo, J.S.; Hu, Y.; Li, Q.X. Degradation pathways of phenanthrene by Sinorhizobium sp. C4. Appl. Microbiol. Biotechnol. 2006, 71, 935–941. [Google Scholar] [CrossRef] [PubMed]
  79. Muratova, A.; Pozdnyakova, N.; Makarov, O.; Baboshin, M.; Baskunov, B.; Myasoedova, N.; Golovleva, L.; Turkovskaya, O. Degradation of phenanthrene by the rhizobacterium Ensifer meliloti. Biodegradation 2014, 25, 787–795. [Google Scholar] [CrossRef]
  80. Turkovskaya, O.; Muratova, A. Plant–bacterial degradation of polyaromatic hydrocarbons in the rhizosphere. Trends Biotechnol. 2019, 37, 926–930. [Google Scholar] [CrossRef]
  81. Sherpa, M.T.; Najar, I.N.; Das, S.; Thakur, N. Distribution of antibiotic and metal resistance genes in two glaciers of North Sikkim, India. Ecotoxicol. Environ. Saf. 2020, 203, 111037. [Google Scholar] [CrossRef] [PubMed]
  82. Yang, S.; Deng, W.; Liu, S.; Yu, X.; Mustafa, G.R.; Chen, S.; Chen, S.; He, L.; Ao, X.; Yang, Y.; et al. Presence of heavy metal resistance genes in Escherichia coli and Salmonella isolates and analysis of resistance gene structure in E. coli E308. J. Glob. Antimicrob. Resist. 2020, 21, 420–426. [Google Scholar] [CrossRef]
  83. García-Cañas, R.; Giner-Lamia, J.; Florencio, F.J.; López-Maury, L. A protease-mediated mechanism regulates the cytochrome c6/plastocyanin switch in Synechocystis sp. PCC 6803. Proc. Natl. Acad. Sci. USA 2021, 118, e2017898118. [Google Scholar] [CrossRef]
  84. Ge, Q.; Zhu, X.; Cobine, P.A.; De La Fuente, L. The copper-binding protein CutC is involved in copper homeostasis and affects virulence in the xylem-limited pathogen Xylella fastidiosa. Phytopathology 2022, 112, 1620–1629. [Google Scholar] [CrossRef]
  85. Hyre, A.; Casanova-Hampton, K.; Subashchandrabose, S. Copper homeostatic mechanisms and their role in the virulence of Escherichia coli and Salmonella enterica. EcoSal Plus 2021, 9, eESP-0014-2020. [Google Scholar] [CrossRef]
  86. Salam, L.B. Metagenomic insights into the microbial community structure and resistomes of a tropical agricultural soil persistently inundated with pesticide and animal manure use. Folia Microbiol. 2022, 67, 707–719. [Google Scholar] [CrossRef] [PubMed]
  87. Toulouse, C.; Metesch, K.; Pfannstiel, J.; Steuber, J. Metabolic reprogramming of Vibrio cholerae impaired in respiratory NADH oxidation is accompanied by increased copper sensitivity. J. Bacteriol. 2018, 200, e00761-17. [Google Scholar] [CrossRef]
  88. Al-Ansari, M.M. Biodetoxification mercury by using a marine bacterium Marinomonas sp. RS3 and its merA gene expression under mercury stress. Environ. Res. 2022, 205, 112452. [Google Scholar] [CrossRef]
  89. Babich, T.L.; Grouzdev, D.S.; Sokolova, D.S.; Tourova, T.P.; Poltaraus, A.B.; Nazina, T.N. Genome analysis of Pollutimonas subterranea gen. nov.; sp. nov. and Pollutimonas nitritireducens sp. nov.; isolated from nitrate-and radionuclide-contaminated groundwater, and transfer of several Pusillimonas species into three new genera Allopusillimonas, Neopusillimonas, and Mesopusillimonas. Antonie Van Leeuwenhoek 2023, 116, 109–127. [Google Scholar] [CrossRef]
  90. Binish, M.B.; Shini, S.; Sinha, R.K.; Krishnan, K.P.; Mohan, M. Mercuric reductase gene (merA) activity in a mercury tolerant sulphate reducing bacterium isolated from the Kongsfjorden. Arctic. Polar Sci. 2021, 30, 100745. [Google Scholar] [CrossRef]
  91. Lo, K.H.; Lu, C.W.; Liu, F.G.; Kao, C.M.; Chen, S.C. Draft genome sequence of Pseudomonas sp. A46 isolated from mercury-contaminated wastewater. J. Basic Microbiol. 2022, 62, 1193–1201. [Google Scholar] [CrossRef]
  92. Almeida, M.C.; Branco, R.; Morais, P.V. Response to vanadate exposure in Ochrobactrum tritici strains. PLoS ONE 2020, 15, e0229359. [Google Scholar] [CrossRef] [PubMed]
  93. Princy, S.; Subramanian, D.; Natarajan, J.; Prabagaran, S.R. Analysis of chromate transporters in bacterial species for Cr(VI) reduction isolated from tannery effluent contaminated site of Dindigul district, Tamil Nadu, India. Geomicrobiol. J. 2021, 38, 598–606. [Google Scholar] [CrossRef]
  94. Semenova, E.M.; Ershov, A.P.; Sokolova, D.S.; Tourova, T.P.; Nazina, T.N. Diversity and biotechnological potential of nitrate-reducing bacteria from heavy-oil reservoirs (Russia). Microbiology 2020, 89, 685–696. [Google Scholar] [CrossRef]
  95. Tóth, E.; Szuróczki, S.; Kéki, Z.; Bóka, K.; Szili-Kovács, T.; Schumann, P. Gellertiella hungarica gen. nov., sp. nov., a novel bacterium of the family Rhizobiaceae isolated from a spa in Budapest. Int. J. Syst. Evol. Microbiol. 2017, 67, 4565–4571. [Google Scholar] [CrossRef]
  96. Wang, E.T.; Tan, Z.Y.; Willems, A.; Fernández-López, M.; Reinhold-Hurek, B.; Martínez-Romero, E. Sinorhizobium morelense sp. nov.; a Leucaena leucocephala-associated bacterium that is highly resistant to multiple antibiotics. Int. J. Syst. Evol. Microbiol. 2002, 52, 1687–1693. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Cell morphology of the strains HO-A22T (a) and SHC 2-14 (b) grown in TEG medium at 22 °C for 72 h. The samples were examined under a scanning electron microscope (Quattro S, “Thermo Fisher Scientific”, Brno-Černovice, Czech Republic) with accelerating voltage 15 kV. Bars, 10 µm.
Figure 1. Cell morphology of the strains HO-A22T (a) and SHC 2-14 (b) grown in TEG medium at 22 °C for 72 h. The samples were examined under a scanning electron microscope (Quattro S, “Thermo Fisher Scientific”, Brno-Černovice, Czech Republic) with accelerating voltage 15 kV. Bars, 10 µm.
Microorganisms 11 02314 g001
Figure 2. Maximum-likelihood phylogenetic tree based on 16S rRNA gene sequences (1365 nucleotide sites) reconstructed with evolutionary model GTR+F+I+G4, showing the position of strains HO-A22T and ‘E. canadensis’ SHC 2-14 with closely related members of the family Rhizobiaceae. Grey circles indicate that the corresponding nodes were recovered in the reconstructed tree based on the maximum-parsimony algorithm; white circles indicate that the corresponding nodes were recovered using the neighbor-joining algorithm; black circles indicate that the corresponding nodes were also recovered based on the neighbor-joining and maximum-parsimony algorithms. Bootstrap values (>50%) are listed as percentages at the branching points. Bar, 0.02 substitutions per nucleotide position. The tree was rooted using Rhizobium leguminosarum USDA 2370T as the outgroup. GenBank accession numbers for 16S rRNA genes are indicated in brackets.
Figure 2. Maximum-likelihood phylogenetic tree based on 16S rRNA gene sequences (1365 nucleotide sites) reconstructed with evolutionary model GTR+F+I+G4, showing the position of strains HO-A22T and ‘E. canadensis’ SHC 2-14 with closely related members of the family Rhizobiaceae. Grey circles indicate that the corresponding nodes were recovered in the reconstructed tree based on the maximum-parsimony algorithm; white circles indicate that the corresponding nodes were recovered using the neighbor-joining algorithm; black circles indicate that the corresponding nodes were also recovered based on the neighbor-joining and maximum-parsimony algorithms. Bootstrap values (>50%) are listed as percentages at the branching points. Bar, 0.02 substitutions per nucleotide position. The tree was rooted using Rhizobium leguminosarum USDA 2370T as the outgroup. GenBank accession numbers for 16S rRNA genes are indicated in brackets.
Microorganisms 11 02314 g002
Figure 3. The maximum-likelihood phylogenetic tree derived from concatenated 120 single copy proteins shows the position of strains HO-A22T and ‘E. canadensis’ SHC 2-14 in relation to closely related members of the family Rhizobiaceae. Phylogenetic analysis was performed with the LG+F+I+G4 model using 37,730 amino acid positions. Black circles indicate that the corresponding nodes were also recovered based on the neighbor-joining and maximum-parsimony algorithms. Bar, 0.02 amino acid substitutions per site. Bootstrap values (>50%) are listed as percentages at the branching points. The tree was rooted using Rhizobium leguminosarum USDA 2370T as the outgroup. Accession numbers for the genomic assemblies are indicated in brackets.
Figure 3. The maximum-likelihood phylogenetic tree derived from concatenated 120 single copy proteins shows the position of strains HO-A22T and ‘E. canadensis’ SHC 2-14 in relation to closely related members of the family Rhizobiaceae. Phylogenetic analysis was performed with the LG+F+I+G4 model using 37,730 amino acid positions. Black circles indicate that the corresponding nodes were also recovered based on the neighbor-joining and maximum-parsimony algorithms. Bar, 0.02 amino acid substitutions per site. Bootstrap values (>50%) are listed as percentages at the branching points. The tree was rooted using Rhizobium leguminosarum USDA 2370T as the outgroup. Accession numbers for the genomic assemblies are indicated in brackets.
Microorganisms 11 02314 g003
Figure 4. Pangenome analysis of Ensifer and Sinorhizobium type strains calculated with Anvi’o v. 7.0. Central dendrogram represents relationships between 22,072 gene clusters (117,350 genes) in all analyzed genomes. Dark circular regions represent genes found in those areas for each genome. The phylogenetic tree is reconstructed using the single copy genes. ANI heatmap in blue squares varies between 70 and 100%.
Figure 4. Pangenome analysis of Ensifer and Sinorhizobium type strains calculated with Anvi’o v. 7.0. Central dendrogram represents relationships between 22,072 gene clusters (117,350 genes) in all analyzed genomes. Dark circular regions represent genes found in those areas for each genome. The phylogenetic tree is reconstructed using the single copy genes. ANI heatmap in blue squares varies between 70 and 100%.
Microorganisms 11 02314 g004
Table 1. Comparison of the similarity level (%) of 16S rRNA genes and genomic indexes between the strains HO-A22T and SHC 2-14 and type strains of the genus Ensifer.
Table 1. Comparison of the similarity level (%) of 16S rRNA genes and genomic indexes between the strains HO-A22T and SHC 2-14 and type strains of the genus Ensifer.
GenomeHO-A22TSHC 2-14
16S rRNAdDDHANIAAIPOCP16S rRNAdDDHANIAAIPOCP
HO-A22T10010010010010099.835.888.791.680.6
SHC 2-1499.835.888.791.780.6100100100100100
E. adhaerens AT99.125.484.184.179.799.225.483.883.578.9
E. canadensis T173T99.835.887.991.880.810092.899.199.188.7
E. morelensis Lc04T99.945.992.495.187.499.935.688.791.479.3
E. sesbaniae CCBAU 65729T99.025.683.784.375.699.125.883.783.874.1
ANI analysis of HO-A22T genome was performed with the genomes of Ensifer strains. ANI values for the strain HO-A22T were within the range of 83.7–92.4%, which were below the species threshold, proposed to be 95% [45,49]. The dDDH values against the reference genomes of Ensifer strains were within the range of 25.4–45.9% and were below the 70% threshold to differentiate bacterial species [45,47]. Thus, analysis of genome of the strain HO-A22T testifies to its belonging to a novel species. The ANI and dDDH values between SHC 2-14 and ‘E. canadensis’ T173T were 99.1% and 92.8%, respectively. These high values indicate that strain SHC 2-14 belongs to the ‘E. canadensis’ species.
Table 2. Characteristics differentiating strains HO-A22T and ‘E. canadensis’ SHC 2-14 from members of the genus Ensifer.
Table 2. Characteristics differentiating strains HO-A22T and ‘E. canadensis’ SHC 2-14 from members of the genus Ensifer.
Characteristics E. canadensisE. canadensisE. adhaerensE. morelensisE. sesbaniae
StrainHO-A22TSHC 2-14T173TATLc04TCCBAU 65729T
Catalase++NDw+
Citrate utilizationND+
UreaseND++
Growth at:
  2% NaCl++++
  3% NaClwww
  15 °C+++++
  37 °Cwww+
  pH 5+
  pH 10w+ww
Acid production (API® 50CH) from:
  D-fucose++++ND
  D-galactose+w++ND
  D-lactose+w+++
  D-melicitoseNDw+
  D-turanose+w++ND
  Dulcitol+NDw
  Gentiobiose+ww+ND
  L-rhamnose+w+++
  L-sorbose+NDw+
  Salicin+++
Major fatty acidsC18:1, C16:0, C14:0 3-OH, C16:1, C17:1C18:1, C16:0, C16:1, C14:0 3-OH, C18:0C18:1, C16:0, C12:0 aldehyde, C18:0, C18:0 cyclo ω8cC18:1, C16:0, C16:1, C14:0 3-OH, C18:0, C18:0 3-OHC18:1, C16:0, C16:0 3-OH, C18:0, C16:1C18:1, C12:0 aldehyde, C16:0, C16:1, C18:0
Major polar lipidsPC, PE, DPG, PGPC, DPG, PE,
APL4, PG
NDPE, PC, PG, DPGPE, DPG, PC, PGPC, PE, PG
G+C content, %61.761.061.062.8 a61.7 b60.4
Isolation sourceInjection water of an oil fieldRadionuclide-contaminated groundwaterRoot nodule of Melilotus albusSoil aRoot nodule of Leucaena leucocephala bRoot nodule of Sesbania cannabina
All strains were Gram-negative, motile, aerobic, non-spore-forming rods. Positive results for all strains were obtained for: growth at pH 7.0–9.0, 28 °C, 1% NaCl; acid production (API® 50CH) from D-fructose, D-glucose, D-mannose, and D-trehalose. Negative results for all strains were obtained for growth at 4% NaCl. Data for ‘E. canadensis’ T173T are from Bromfield et al. [6], and those of E. sesbaniae CCBAU 65729T are from Wang et al. [4]. Data for all other strains are from this study, except as labeled: a, Data from Tóth et al. [95]; b, Data from Wang et al. [96]. Designations: +, positive; –, negative; w, weakly positive; APL4, unidentified aminophospholipid; DPG, diphosphatidylglycerols; PC, phosphatidylcholines; PE, phosphatidylethanolamines; PG, phosphatidylglycerols; ND, not determined.
Table 3. Protologue description of Ensifer oleiphilus sp. nov.
Table 3. Protologue description of Ensifer oleiphilus sp. nov.
ParameterEnsifer oleiphilus sp. nov.
Genus nameEnsifer
Species nameEnsifer oleiphilus
Species statussp. nov.
Species etymologyo.le.i.phi’lus. L. neut. n. oleum oil; N.L. masc. adj. philus (from Gr. masc. adj. philos) loving; N.L. masc. adj. oleiphilus loving oil, referring to the ability of type strain to utilize hydrocarbons of crude oil
Description of the new taxon and diagnostic traitsCells are Gram-negative, motile, aerobic, non-spore-forming rods. Colonies up to 8.0 in diameter are creamy, smooth, shiny, and round after 5 days of cultivation in the TEG medium at 22 °C. Growth occurs on TEG, PCA, and LB media, in the presence of 0–3.0% (w/v) NaCl (optimum, 0–1.2% NaCl), at pH 6.0–9.8 (optimum, pH 6.6–8.0) and at 15–33 °C (optimum, 22–28 °C). Catalase-positive, but negative for urease and citrate utilization. Chemoorganoheterotrophic, facultatively anaerobic. Nitrate is reduced to the nitrous oxide. Acid is produced from D-fucose, D-galactose, D-lactose, D-turanose, dulcitol, gentiobiose, L-rhamnose, L-sorbose, and salicin, but not from D-melicitose. The predominant fatty acids are C18:1, C16:0, C14:0 3-OH, C16:1, C17:1, and C18:0. Major menaquinone is Q10. The major polar lipids are phosphatidylcholine, phosphatidylethanolamine, diphosphatidylglycerol, and phosphatidylglycerol. The genome size of the type strain is 6.76 Mb with a G+C value of 61.7%. The strain type, HO-A22T (=VKM B-3646T = KCTC 92427T), was isolated from a denitrifying enrichment obtained from an injection water of the Cheremukhovskoe oil field (Tatarstan, Russia). The GenBank/EMBL/DDBJ accession number for the 16S rRNA gene sequence is MT495799 and the genomic assembly accession number is GCF_013371465.1.
Country of originRussian Federation
Region of originNurlat region, Tatarstan
Date of isolation2018
Source of isolationInjection water of an oil field
Sampling date2016
Latitude55°0′ N
Longitude51°1′ E
Depth (meters below sea level)850–1300
Number of strains in study1
Information related to the Nagoya ProtocolNot applicable
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.

Share and Cite

MDPI and ACS Style

Ershov, A.P.; Babich, T.L.; Grouzdev, D.S.; Sokolova, D.S.; Semenova, E.M.; Avtukh, A.N.; Poltaraus, A.B.; Ianutsevich, E.A.; Nazina, T.N. Genome Analysis and Potential Ecological Functions of Members of the Genus Ensifer from Subsurface Environments and Description of Ensifer oleiphilus sp. nov. Microorganisms 2023, 11, 2314. https://doi.org/10.3390/microorganisms11092314

AMA Style

Ershov AP, Babich TL, Grouzdev DS, Sokolova DS, Semenova EM, Avtukh AN, Poltaraus AB, Ianutsevich EA, Nazina TN. Genome Analysis and Potential Ecological Functions of Members of the Genus Ensifer from Subsurface Environments and Description of Ensifer oleiphilus sp. nov. Microorganisms. 2023; 11(9):2314. https://doi.org/10.3390/microorganisms11092314

Chicago/Turabian Style

Ershov, Alexey P., Tamara L. Babich, Denis S. Grouzdev, Diyana S. Sokolova, Ekaterina M. Semenova, Alexander N. Avtukh, Andrey B. Poltaraus, Elena A. Ianutsevich, and Tamara N. Nazina. 2023. "Genome Analysis and Potential Ecological Functions of Members of the Genus Ensifer from Subsurface Environments and Description of Ensifer oleiphilus sp. nov." Microorganisms 11, no. 9: 2314. https://doi.org/10.3390/microorganisms11092314

APA Style

Ershov, A. P., Babich, T. L., Grouzdev, D. S., Sokolova, D. S., Semenova, E. M., Avtukh, A. N., Poltaraus, A. B., Ianutsevich, E. A., & Nazina, T. N. (2023). Genome Analysis and Potential Ecological Functions of Members of the Genus Ensifer from Subsurface Environments and Description of Ensifer oleiphilus sp. nov. Microorganisms, 11(9), 2314. https://doi.org/10.3390/microorganisms11092314

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop