Azohydromonas caseinilytica sp. nov., a Nitrogen-Fixing Bacterium Isolated From Forest Soil by Using Optimized Culture Method

A bacterial strain, designated strain G-1-1-14^T, was isolated from Kyonggi University forest soil during a study of previously uncultured bacterium. The cells of strain G-1-1-14^T were motile by means of peritrichous flagella, Gram-stain-negative, rod-shaped, and able to grow autotrophically with hydrogen and fix nitrogen. Phylogenetic analysis based on 16S rRNA gene sequence indicated that strain G-1-1-14^T belonged to the genus Azohydromonas. The closest species of strain G-1-1-14^T were Azohydromonas ureilytica UCM-80^T (98.4% sequence similarity), Azohydromonas lata IAM 12599^T (97.5%), Azohydromonas riparia UCM-11^T (97.1%), and Azohydromonas australica IAM 12664^T (97.0%). The genome of strain G-1-1-14^T was 6,654,139 bp long with 5,865 protein-coding genes. The genome consisted of N₂-fixing genes (nifH) and various regulatory genes for CO₂ fixation and H₂ utilization. The principal respiratory quinone was ubiquinone-8, and the major polar lipids were phosphatidylethanolamine, diphosphatidylglycerol, and phosphatidylglycerol. The major fatty acids were summed feature 3 (iso-C_15:0 2-OH and/or C_16:1ω7c), C_16:0, summed feature 8 (C_18:1ω7c and/or C_18:1ω6c), and cyclo-C_17:0. The DNA G + C content was 69.9%. The average nucleotide identity (OrthoANI), in silico DNA–DNA hybridization (dDDH), and conventional DDH relatedness values were below the species demarcation values for novel species. Based on genomic, genetic, phylogenetic, phenotypic, and chemotaxonomic characterizations, strain G-1-1-14^T represents a novel species within the genus Azohydromonas, for which the name Azohydromonas caseinilytica sp. nov. is proposed. The type strain is G-1-1-14^T (= KACC 21615^T = NBRC 114390^T).

Introduction

Metagenomic, next-generation sequencing, and whole-genome sequencing technology help to reclassify and correct the taxonomic position of bacterial species in systematic studies. In recent years, various taxa have been reclassified based on 16S ribosomal RNA (rRNA) gene and whole-genome analysis, and genera are more complete than they previously were. Similarly, the genus Azohydromonas was first proposed by reclassifying the genus [Alcaligenes] (Xie and Yokota, 2005). On the other hand, these technologies revealed that most environmental microorganisms are still “not-yet cultured” on synthetic media and comprise around 99% of microorganisms previously uncultured (Vartoukian et al., 2010; Chaudhary et al., 2019). The simple to complex reasons for failing to cultivate these uncultivated microorganisms include the lack of correct nutrients, use of nutrient-rich media, unsuitable pH, insufficient incubation time, inappropriate temperature, need of specific growth signal, dependence on other microorganisms, or failure to maintain a natural growth environment (Vartoukian et al., 2010; Pulschen et al., 2017; Chaudhary et al., 2019). We have isolated a bacterium, designated strain G-1-1-14^T, from Kyonggi University forest soil using low nutrient and long enrichment time. In addition, the isolated strain showed best match with a previously uncultured bacterium (non-type material sequences) within Megablast (GenBank) and belonged to the genus Azohydromonas.

Currently, the genus Azohydromonas comprise only five species – Azohydromonas lata, Azohydromonas australica, Azohydromonas riparia, and Azohydromonas ureilytica (Xie and Yokota, 2005; Nguyen and Kim, 2017) – with validly published names in addition to “Azohydromonas aeria” (Xue et al., 2020), which has not been validated yet. Members of the genus Azohydromonas are characterized as motile by means of flagella, non-spore-forming, Gram-stain-negative, rod- or coccid-shaped, DNA G + C content between 69 and 72 mol%, ubiquinone-8 (Q-8) as the major respiratory quinone, summed feature 3 (iso-C_15:0 2-OH and/or C_16:1ω7c), C_16:0, and summed feature 8 (C_18:1ω7c and/or C_18:1ω6c) as the principal fatty acids in addition to cells accumulating PHB (poly-β-hydroxybutyrate) granules as storage material (Xie and Yokota, 2005; Nguyen and Kim, 2017). The members are able to fix nitrogen and grow autotrophically with hydrogen (Xie and Yokota, 2005). This study describes the polyphasic study and taxonomic position of strain G-1-1-14^T isolated from forest soil.

Materials and Methods

Isolation and Ecology

Strain G-1-1-14^T was isolated from Kyonggi University forest soil, geographically located at Suwon-si, Gyeonggi-do, South Korea (37°18′5″ N and 127°1′56″ E). The pH of the soil was 6.8. A modified culture method using a six-well Transwell plate (Corning Inc., NY, United States) was used for isolation. Debris-free sieved soil (∼3 g) was kept on the bottom of the Transwell plate and 3 ml of diluted (1/10) Reasoner’s 2A (R2A) broth (MBcell; KisanBio, Seoul, South Korea) was added to the insert. After that, 100 μl of soil suspension (1 g of soil in 9 ml of distilled water) was added to the insert. Then, the Transwell plate was incubated in a shaker at 130 rpm at 28°C for 6 weeks. After 6 weeks, enriched culture was serially diluted and then 100 μl of each dilution was spread on 1/10 R2A agar plates (the agar used in this research was Agar A; Bio Basic, Toronto, ON, Canada) and incubated at 28°C for 4 weeks (Dahal and Kim, 2018c). The short-term maintenance and the long-term preservation of the isolate were done as described previously (Dahal et al., 2018).

16S rRNA Phylogeny

The genomic DNA of strain G-1-1-14^T was extracted using InstaGene Matrix kit (Bio-Rad, Hercules, CA, United States) according to the manufacturer’s instruction. Amplification of 16S rRNA gene was done by PCR using primers 27F and 1492R (Frank et al., 2008). Sequencing was carried out using a 3770XL DNA analyzer with a BigDye Terminator cycle sequencing kit v.3.1 (Applied Biosystems, CA, United States). Near-complete sequences of the 16S rRNA genes were assembled with SeqMan software (DNASTAR, Inc., Madison, WI, United States). For type material sequences, the closest phylogenetic neighbors were identified by using the EzBioCloud server (Yoon et al., 2017b). In addition, non-type material sequences were compared with the top hits of Megablast (GenBank). All the 16S rRNA gene sequences of the closest phylogenetic members were retrieved from the whole-genome sequences (WGS) and/or GenBank database and aligned using SINA¹ (Pruesse et al., 2012). Phylogenetic trees were reconstructed using MEGA7 (Kumar et al., 2016). Neighbor-joining and maximum-likelihood trees were reconstructed using the Kimura two-parameter model based on 1,000 bootstrap replications and partial deletion with 95% site coverage cutoff (Felsenstein, 1981; Saitou and Nei, 1987). In addition, maximum parsimony tree was inferred with the “Mini-Mini Heuristic” search method.

Genome Features

Whole-genome-based approaches were used for further analysis of the taxonomic status of the novel strain. For whole-genome sequencing, the genomic DNA was extracted by using DNeasy Blood and Tissue kits (Qiagen). Whole-genome shotgun sequencing of strain G-1-1-14^T was performed by Macrogen (Seoul, South Korea) using the Illumina HiSeq platform and assembled by SPAdes (Bankevich et al., 2012). The authenticity of the genome assembly was checked by comparing the 16S rRNA gene sequences using the NCBI Align Sequences Nucleotide BLAST tool (Zhang et al., 2000) and the potential contamination was checked by ContEst16S algorithm (Lee et al., 2017). After analysis, the whole-genome sequence was annotated using NCBI Prokaryotic Genome Annotation Pipeline (PAPG) (Tatusova et al., 2016) and Rapid Annotations using the Subsystems Technology (RAST) server (Aziz et al., 2008). The 92 core genes were extracted from genomes using Prodigal v2.6.3 (Hyatt et al., 2010) and hmmsearch v3.1b2 (Eddy, 2011). The amino acid sequences of 92 core genes were aligned by using MAFFT 7.310 (Katoh and Standley, 2013) and concatenated into a single alignment. The alignment position that had a gap >50% were excluded. Then, the phylogenomic tree was inferred by using FastTree (Price et al., 2010) and RAxML (Stamatakis, 2014) and viewed using MEGA X v10.1 (Kumar et al., 2018). The branch support inference was based on 100 non-parametric bootstrap replicates, and the branch supports of the phylogenomic tree were evaluated using gene support index (GSI). The genome-based relatedness between strain G-1-1-14^T and the phylogenetically related type strains was determined based on average nucleotide identity (ANI) using the OrthoANI algorithm (Yoon et al., 2017a). The digital DNA–DNA hybridization (dDDH) was calculated in silico by the Genome-to-Genome Distance Calculator using the blast method (Meier-Kolthoff et al., 2013). In addition, conventional DDH was measured fluorometrically using photobiotin-labeled DNA probes and microdilution plates as recommended by Ezaki et al. (1989). Moreover, annotation and analysis of the secondary metabolite biosynthesis genes were carried out by using the antiSMASH server (Blin et al., 2019). The COG (Clusters of Orthologous Group) functional categories were assigned by searching against the KEGG (Kyoto Encyclopedia of Genes and Genomes) database (Kanehisa and Goto, 2000). Furthermore, the CRISPR gene and Cas cluster were analyzed using the CRISPRCasFinder online server².

Physiology and Chemotaxonomy

The cell morphology of strain G-1-1-14^T, grown on R2A agar for 5 days at 28°C, was examined by transmission electron microscopy (Talos L120C, FEI). Colony morphology was observed by a Zoom Stereo Microscope (SZ61; Olympus, Tokyo, Japan). Gram staining was performed as described previously (Doetsch, 1981). Motility was tested in the R2A medium containing 0.4% (w/v) agar. Oxidase activity was determined using 1% (w/v) tetra-methyl-p-phenylenediamine dihydrochloride. Catalase activity was assessed using 3% (v/v) hydrogen peroxide (H₂O₂). Growth at various temperatures (4–50°C) on R2A agar plates was observed for 10 days. Growth was determined on various media including tryptone soya agar (TSA; Oxoid), R2A agar, nutrient agar (NA; Oxoid), sorbitol MacConkey agar (MA; Oxoid), potato dextrose agar (PDA; Becton), marine agar 2216 (Becton), brain heart infusion (BHI) agar (Oxoid), veal infusion agar (Becton), and Luria–Bertani agar (LBA; Oxoid). The salt tolerance of strain G-1-1-14^T was examined in R2A broth supplemented with NaCl (0–5%, w/v, at 0.5% intervals). The pH range for growth was determined by cultivation at 28°C in R2A broth adjusted to pH 4–12 (at 0.5 pH unit increment) prior to sterilization using citrate/NaH₂PO₄ buffer (for pH 4.0–5.5), phosphate buffer (for pH 6–7.5), Tris buffer (for pH 8–10) and 5 M NaOH (for pH 10.5–12.0) (Dahal and Kim, 2018b). The hydrolysis of Tween 80, Tween 60, and Tween 40 was assessed using the method of Smibert and Krieg (1994). Anaerobic growth was examined on R2A agar at 28°C for 10 days by using the BD GasPak EZ Gas Generating Pouch System. The hydrolysis of starch, chitin, carboxymethyl (CM)-cellulose, tyrosine, and casein was examined as previously described (Dahal and Kim, 2018a). A DNase activity assay was performed with DNase agar (Oxoid). The presence of spores was examined by staining with malachite green. Autotrophic growth on hydrogen and nitrogen fixation was done as recommended by Pedrosa et al. (1980). For the nitrogen fixation test, strain G-1-1-14^T was grown on N-free semi-solid (NFb) medium and N-free agar plate with a bromothymol blue indicator (da Silva Lima et al., 2014). In addition, for autotrophic growth in hydrogen, the culture was supplied with 10% O₂, 10% CO₂, 20% N₂, and 60% H₂ (v/v) in NFb medium. Other physiological tests were performed using API 20NE and API ID 32GN kits (bioMérieux, Marcy-l’Étoile, France). Enzyme activities were observed using an API ZYM kit (bioMérieux) as per the manufacturer’s instructions.

For fatty acid analysis, the cells of strain G-1-1-14^T and the reference strains were harvested from the same culture condition during the late log phase (at 28°C for 4 days on R2A agar plate). The cellular fatty acids were extracted by using the MIDI protocol (Sherlock Microbial Identification System, version 6.0B), analyzed with a gas chromatograph (GC; HP 6890 Series GC System, Hewlett Packard), and identified using the TSBA6 database of the Microbial Identification System (Sasser, 1990). Polar lipids and isoprenoid quinones were extracted from freeze-dried cells according to the procedures described by Minnikin et al. (1984). Appropriate detection reagents were used to identify the spots (Komagata and Suzuki, 1988).

Results and Discussion

Phylogenetic Analysis

The nucleotide sequence of the 16S rRNA gene of strain G-1-1-14^T has been deposited in the GenBank/EMBL/DDBJ database under the accession MN685324. Preliminary comparisons with the 16S rRNA gene sequences in GenBank showed top hits with previously uncultured bacterial clones (Supplementary Figure 1). Among the five closely related uncultured bacterial clones, TSNR003_118, TSSUR003_P21, and TS8 have been isolated from rice paddy soil, whereas bacterial clones BJ201307-105 and SH201208-30 are from rainwater. These results showed that the closest members of strain G-1-1-14^T could be isolated not only from the soil but also from rainwater. In addition, the 16S rRNA gene sequence of strain G-1-1-14^T was analyzed with the EzBioCloud server against the type strain sequences. Strain G-1-1-14^T belonged to the family Alcaligenaceae of the order Burkholderiales and was most closely related to A. ureilytica UCM-80^T (98.4% sequence similarity), A. lata IAM 12599^T (97.5%), A. riparia UCM-11^T (97.1%), and A. australica IAM 12664T (97.0%). Strain G-1-1-14^T was well clustered with the other members of the genus Azohydromonas in the neighbor-joining (NJ), maximum likelihood (ML), and maximum parsimony (MP) trees (Figure 1 and Supplementary Figures 2, 3). In addition, a monophyletic clade formed within the genus Azohydromonas with a strong bootstrap value validated for strain G-1-1-14^T as a novel member of genus Azohydromonas (Figure 1).

FIGURE 1

Figure 1. Maximum likelihood tree based on 16S rRNA gene sequences showing the phylogenetic position of strain G-1-1-14^T among the closely related members of the order Burkholderiales. Filled circles indicate nodes recovered by all three treeing methods (neighbor-joining, maximum likelihood, and maximum parsimony). The numbers at the nodes indicate the percentage of 1,000 bootstrap replicates yielding this topology; only values >50% are shown. Burkholderia cepacia ATCC 25416^T was used as an outgroup. The GenBank accession numbers are given in parentheses. Bar, 0.020 substitution per nucleotide position. Bold values indicates proposed species.

For physiology, biochemical, quinone, and fatty acid analyses, the reference strains A. ureilytica UCM-80^T, A. lata KACC 15149^T, A. riparia UCM-11^T, and A. australica KACC 15148^T were selected and analyzed under identical conditions.

Genome Analysis

The whole-genome shotgun sequence of strain G-1-1-14^T has been deposited at DDBJ/ENA/GenBank under the accession JABBFW000000000. The whole-genome sequence of strain G-1-1-14^T was 6,654,139 bp long with DNA G + C of 69.9%. The genome was assembled in 98 scaffolds with an N50 value of 195,393 bp, 351 subsystems, and genome coverage of 162.0× (Table 1). The genome features of the type species of A. lata and A. australica are also presented in Table 1. The DNA G + C content was calculated based on whole-genome sequences. The OrthoANI values of strain G-1-1-14^T for A. lata NBRC 102462^T, A. australica DSM 1124^T, and “A. aeria” t3-1-3 were 84.6, 86.5, and 85.3%, respectively (Supplementary Table 1). Similarly, the in silico DDH values for A. lata NBRC 102462^T, A. australica DSM 1124^T, and “A. aeria” t3-1-3 were 28.6, 31.5, and 29.7, respectively (Supplementary Table 1). These obtained values were below the threshold ANI value of 95.0–96.0% and dDDH value of <70% used for species delineation (Richter and Rosselló-Móra, 2009; Meier-Kolthoff et al., 2013). The DNA–DNA hybridization of strain G-1-1-14^T with reference strains A. riparia UCM-11^T and A. ureilytica UCM-80^T showed DDH relatedness of 45.8 ± 2.7% and 31.4 ± 2.6%, respectively. The DNA–DNA relatedness between these species clearly showed that strain G-1-1-14^T differs genetically from the Azohydromonas type strains at the species level (Wayne et al., 1987). In addition, the phylogenomic tree constructed using concatenated 92 core genes also proved that strain G-1-1-14^T is a novel member of the genus Azohydromonas (Figure 2). Moreover, due to the unavailability of the genome data of A. ureilytica UCM-80^T and A. riparia UCM-11^T, these strains were not used in the phylogenomic tree reconstruction. However, strain G-1-1-14^T well clustered with the other members of the genus Azohydromonas and distinguished with other clusters (Figure 2). The phylogenomic inference is supported by high bootstrap and GSI values. Furthermore, multilocus sequence analysis (MLSA) trees based on the rpoB and recA genes showed high robustness of strain G-1-1-14^T with the Azohydromonas members (Supplementary Figures 4, 5).

TABLE 1

Table 1. Genome features of strain G-1-1-14^T, Azohydromonas australica DSM 1124^T, and Azohydromonas lata NBRC 102462^T.

FIGURE 2

Figure 2. Phylogenomic tree based on the concatenated alignment of 92 core genes. Bar, 0.050 substitution per nucleotide position. Bootstrap and gene support index (GSI) values are shown in the branch nodes. Values in parentheses are GSI. Bold values indicates proposed species.

The genome of strain G-1-1-14^T contains nitrogen fixation regulatory genes (nif genes) such as nifA, nifB, nifD, nifE, nifH, nifK, nifN, nifO, nifQ, nifT, nifV, nifW, nifX, nifX2, and nifZ (Supplementary Tables 2, 3). Similar nitrogen fixation regulatory genes were also present in the genomes of A. australica DSM 1124^T and A. lata NBRC 102462^T. The presence of these genes showed that strain G-1-1-14^T fixes atmospheric nitrogen. In addition, the genome showed each chain (small and large) of ribulose 1,5-bisphosphate (RuBP) carboxylase, which is an enzyme in the first principal step of carbon fixation. Moreover, the genome consists of the high-affinity carbon uptake protein (hat/hatR), which also regulates CO₂ fixation. The genome consists of the genes (hypF, hypA, hypC, hypB, hypD, and hypE) associated with the [NiFe]-hydrogenase regulatory proteins that are responsible for the uptake of H₂ gas as a source of energy. Furthermore, the genome contains the genes (hoxA, hoxC, hoxO/hyaE, hoxB, hoxT/hybE, and hypC/hoxL) responsible for hydrogen-sensing regulatory systems in addition to the hydrogenase structural genes (hupV, hupU, and hupF) responsible for hydrogenase activity. These results, along with the data of nitrogen fixation and autotrophic growth assessment, proved that strain G-1-1-14^T fixes nitrogen and grows autotrophically in the presence of H₂ gas. Moreover, other members of the genus Azohydromonas have also been reported for autotrophic growth (Xie and Yokota, 2005; Nguyen and Kim, 2017).

Several members of Azohydromonas are well-known nitrogen fixers (Xie and Yokota, 2005). Biological nitrogen fixation (BNF) is a key step in the nitrogen cycle as it transforms atmospheric nitrogen into ammonium (Pereira e Silva et al., 2013). The nitrogen cycle and the related microorganisms play crucial roles in the ecosystem, affecting both the agricultural field and climate sectors (Khanal and Lee, 2020). Strain G-1-1-14^T showed a direct association with the nitrogen cycle as it reported the ability to fix atmospheric nitrogen. In this regard, this strain can play a valuable ecological role that contributes to maintaining the nitrogen cycle in the environment. Furthermore, strain G-1-1-14^T grows autotrophically with the help of H₂ gas. The genome of this strain is also reported to contain the genes hat/hatR, which regulate CO₂ fixation. CO₂ fixation is an important phenomenon in the carbon cycle that reduces the CO₂ emission in the atmosphere. CO₂ is a primary greenhouse gas that is responsible for global warming (Salehizadeh et al., 2020). Any attempt of reducing CO₂ emission in the atmosphere can be a milestone to protecting the natural ecosystem. In addition, strain G-1-1-14^T produces P3(HB) (and regulatory proteins WP_169159239, WP_169160336, and WP_169163299; Supplementary Figure 6), which is also responsible for reducing the CO₂ releases in the atmosphere (Zafar et al., 2014). Species of Azohydromonas have been revealed to synthesize various poly-β-hydroxyalkanoates (PHAs) (Zafar et al., 2014; Sharma et al., 2017). PHAs have been widely applied for the production of biodegradable and biocompatible plastics (Zafar et al., 2014). From this point of view, strain G-1-1-14^T can be considered as a promising bioresource to mitigate plastic-associated problems and environmental issues. Overall, strain G-1-1-14^T can play a key role in the both nitrogen and carbon cycles, suggesting that this strain has a significant potential ecological role in the natural habitat.

The genome of strain G-1-1-14^T consists of 13 putative biosynthetic gene clusters (BGCs) that are responsible for various secondary metabolites including terpene, burkholderic acid, bacteriocin, hserlactone, aryl polyene, lanthipeptide, non-ribosomal peptide synthetase (NRPS), NRPS-like, and N-acetylglutaminylglutamine amide (NAGGN) (Supplementary Table 4). The predicted secondary metabolites may have potential ecological roles, such as hserlactone, which may be related to the communication between fungi and bacteria (Shiner et al., 2005), and aryl polyene, which could protect bacteria from reactive oxidation (Schöner et al., 2016). In addition, RAST analyses revealed the presence of metabolic genes for thiazole/oxazole-modified microcin synthesis, alkaloid biosynthesis, and auxin biosynthesis. Furthermore, DNA metabolism, nitrogen metabolism, carbohydrate metabolism, CO₂ fixation, phosphorus metabolism, and other metabolic as well as physiologically relevant genes were also present in the genome of strain G-1-1-14^T (Supplementary Figure 7). Moreover, 88 genes for secondary metabolite biosynthesis, transport, and catabolism, 1,347 genes with unknown functions, 258 genes for energy production and conversion, and 212 genes for inorganic ion transport and metabolism have been detected in COG functional categories (Figure 3).

FIGURE 3

Figure 3. Clusters of Orthologous Group (COG) functional classification of proteins in the genome of strain G-1-1-14^T.

Phenotypic and Chemotaxonomic Characteristics

The cells of strain G-1-1-14^T (Supplementary Figure 6) were rod-shaped, Gram-stain-negative, aerobic, and motile by means of peritrichous flagella. The cells were able to grow on N₂-free semi-solid medium supplemented with 10% O₂, 10% CO₂, 20% N₂, and 60% H₂ (v/v) (OD₆₀₀ = 0.3). The cells were 2.6–2.9 μm long and 1.2–1.7 μm wide. Colonies on R2A agar were straw white in color, circular, entire, and convex. The differential phenotypic characteristic features of strain G-1-1-14^T are presented in Table 2 with all the type species of the genus Azohydromonas. Having the highest salt tolerance (4.0%); casein hydrolysis; assimilation of 3-hydroxybutyric acid, D-ribose, lactic acid, N-acetyl-glucosamine, propionic acid, salicin, and valeric acid; and the positive activities of trypsin, valine arylamidase, α-chymotrypsin, α-galactosidase, and β-glucosidase clearly distinguished strain G-1-1-14^T from the other type strains (Table 2).

TABLE 2

Table 2. Phenotypic characteristics that differentiate strain G-1-1-14^T from the type strains of the genus Azohydromonas.

The major cellular fatty acids of strain G-1-1-14^T were summed feature 3 (C_16:1ω7c and/or C_18:1ω6c), C_16:0, summed feature 8 (C_18:1ω7c and/or C_18:1ω6c), and cyclo-C_17:0. Although the major fatty acids of the other type strains of the genus Azohydromonas are similar, the presence of minor fatty acids C_14:0 2-OH, C_15:0 3-OH, iso-C_15:0 3-OH, and iso-C_19:0 and the differences of the major and minor fatty acid compositions showed characteristic differences of strain G-1-1-14^T compared to other reference cultures (Table 2). The sole respiratory quinone was Q-8, the same as the major respiratory quinone reported for the other members of the genus Azohydromonas. The main polar lipids for strain G-1-1-14^T were phosphatidylethanolamine (PE), diphosphatidylglycerol (DPG), and phosphatidylglycerol (PG) (Supplementary Figure 8). However, the presence of an unidentified glycolipid (GL), unidentified aminophospholipid (APL), and other minor polar lipids (L1–L4) distinguished strain G-1-1-14^T from the other closely related type strains (Supplementary Figure 8; Nguyen and Kim, 2017).

Conclusion

Based on the above discussed genomic, phylogenetic, phenotypic, and chemotaxonomic characteristic differences, strain G-1-1-14^T represents a novel member in the genus Azohydromonas, for which the name Azohydromonas caseinilytica sp. nov. is proposed.

Description of Azohydromonas caseinilytica sp. nov.

Azohydromonas caseinilytica (ca.se.i.ni.ly’ti.ca. N.L. n. caseinum casein; Gr. adj. lytikos able to dissolve; N.L. fem. adj. lytica dissolving; N.L. fem. adj. caseinilytica casein-dissolving).

The cells (1.2–1.7 μm wide and 2.6–2.9 μm long) are Gram-stain-negative, non-spore-forming, aerobic, rod-shaped, and motile by means of peritrichous flagella. Colonies on R2A agar are straw white colored, entire, convex, and circular. The colony size of the cells is 1–2 mm on R2A agar for 7 days at 28°C. Cells are positive in the catalase and oxidase tests. The cells grow well on R2A, TSA, and veal infusion agar, while no growth occurs on LBA, BHI, PDA, MA, and marine 2216 agar, but poorly grow on NA. The cells are able to grow autotrophically on hydrogen and fix nitrogen. They grow optimally in the absence of NaCl, but tolerate up to 4.0% (w/v). The cells grow at 15–45°C (optimum, 25–35°C) and pH 5.5–11.0 (optimum, pH 6.0–8.0) and produce poly 3-hydroxybutyrate. They are able to hydrolyze casein, tyrosine, Tween 60, Tween 40, and aesculin, but not CM-cellulose, starch, chitin, DNA, Tween 80, and urea. Nitrate is reduced, but nitrite is not. Glucose is not fermented. The type strain is positive for acid phosphatase, alkaline phosphatase, cystine arylamidase, esterase lipase (C8), leucine arylamidase, trypsin, valine arylamidase, α-chymotrypsin, α-galactosidase, and β-glucosidase; weakly positive for esterase (C4) and napthol-AS-BI-phosphohydrolase; and is negative for lipase (C14), N-acetyl-β-glucosaminidase, α-fucosidase, α-glucosidase, α-mannosidase, β-galactosidase, and β-glucuronidase. The following substrates are assimilated: 3-hydroxybutyric acid, D-glucose, D-mannitol, D-ribose, D-saccharose, lactic acid, L-alanine, L-proline, L-rhamnose, L-serine, N-acetyl-glucosamine, potassium gluconate, propionic acid, salicin, and valeric acid. The predominant respiratory quinone is Q-8. The principal cellular fatty acids are summed feature 3 (iso-C_15:0 2-OH and/or C_16:1ω7c), C_16:0, summed feature 8 (C_18:1ω7c and/or C_18:1ω6c), and cyclo-C_17:0. The major polar lipids are phosphatidylethanolamine, diphosphatidylglycerol, and phosphatidylglycerol. The DNA G + C content of the type strain is 69.9%.

The type strain, G-1-1-14^T (= KACC 21615^T = NBRC 114390^T), was isolated from forest soil, geographically located at Suwon-si, Gyeonggi-do, South Korea (37°18′5″ N and 127°1′56″ E). The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA gene sequence and the whole-genome sequence of strain G-1-1-14^T are MN685324 and JABBFW000000000, respectively.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

RD and DC conceived, designed, and conducted all the experiments. D-UK interpreted the data. JK coordinated and supervised the study. RD, DC, and D-UK analyzed all the data and prepared the manuscript. All authors read, discussed, edited, and approved the final draft of the manuscript.

Funding

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT; No. 2019R1F1A1058501).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Ms. Geeta Chhetri (Dongguk University) for help with the phylogenomic tree reconstruction.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.647132/full#supplementary-material

Supplementary Figure 1 | Neighbor–joining tree based on 16S rRNA gene sequences showing the phylogenetic position of strain G-1-1-14^T.

Supplementary Figure 2 | Neighbor–joining tree based on 16S rRNA gene sequences showing the phylogenetic position of strain G-1-1-14^T.

Supplementary Figure 3 | Maximum–parsimony tree based on 16S rRNA gene sequences showing the phylogenetic position of strain G-1-1-14^T.

Supplementary Figure 4 | MLSA (Multilocus sequence analyses) tree based on recA genes showing the phylogenetic position of strain G-1-1-14^T among the closest phylogenetic members.

Supplementary Figure 5 | MLSA (Multilocus sequence analyses) tree based on rpoB genes showing the phylogenetic position of strain G-1-1-14^T among the closest phylogenetic members.

Supplementary Figure 6 | Transmission electron microscopic image of strain G-1-1-14^T grown on R2A agar at 28°C for 7 days.

Supplementary Figure 7 | Subsystem feature counts for strain G-1-1-14^T analyzed by RAST (Rapid Annotation using Subsystem technology).

Supplementary Figure 8 | Thin-layer chromatograms of the polar lipids from strain G-1-1-14^T determined on TLC (thin layer chromatography).

Supplementary Table 1 | Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) of strain G-1-1-14^T with reference strains.

Supplementary Table 2 | The genes and associated functional roles in strain G-1-1-14^T genome revealed by RAST (Rapid Annotation using Subsystem technology) analyses.

Supplementary Table 3 | N₂-fixation regulatory proteins in G-1-1-14^T genome.

Supplementary Table 4 | The distribution of biosynthetic gene clusters in Azohydromonas caseinilytica G-1-1-14^T.

Supplementary Table 5 | Cellular fatty acid profiles (% of totals) of strain G-1-1-14^T and other type strains of the genus Azohydromonas.

Footnotes

1. ^ https://www.arb-silva.de/aligner/
2. ^ https://crisprcas.i2bc.paris-saclay.fr/

References

Aziz, R. K., Bartels, D., Best, A. A., DeJongh, M., Disz, T., Edwards, R. A., et al. (2008). The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genomics 9:75. doi: 10.1186/1471-2164-9-75

PubMed Abstract | CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Blin, K., Shaw, S., Steinke, K., Villebro, R., Ziemert, N., Lee, S. Y., et al. (2019). antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 47, W81–W87.

Google Scholar

Chaudhary, D. K., Khulan, A., and Kim, J. (2019). Development of a novel cultivation technique for uncultured soil bacteria. Sci. Rep. 9:6666.

Google Scholar

da Silva Lima, L., Olivares, F. L., Rodrigues, de Oliveira, R., Vega, M. R. G., et al. (2014). Root exudate profiling of maize seedlings inoculated with Herbaspirillum seropedicae and humic acids. Chem. Biol. Technol. Agric. 1:23.

Google Scholar

Dahal, R. H., and Kim, J. (2018a). Dyadobacter flavus sp. nov. and Dyadobacter terricola sp. nov., two novel members of the family Cytophagaceae isolated from forest soil. Arch. Microbiol. 200, 1067–1074. doi: 10.1007/s00203-018-1521-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahal, R. H., and Kim, J. (2018b). Flavobacterium ureilyticum sp. nov., a novel urea hydrolysing bacterium isolated from stream bank soil. Antonie Van Leeuwenhoek 111, 2131–2139. doi: 10.1007/s10482-018-1105-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahal, R. H., and Kim, J. (2018c). Fluviicola kyonggii sp. nov., a bacterium isolated from forest soil and emended description of the genus Fluviicola. Int. J. Syst. Evol. Microbiol. 68, 1885–1889. doi: 10.1099/ijsem.0.002759

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahal, R. H., Chaudhary, D. K., and Kim, J. (2018). Pinisolibacter ravus gen. nov., sp. nov., isolated from pine forest soil and allocation of the genera Ancalomicrobium and Pinisolibacter to the family Ancalomicrobiaceae fam. nov., and emendation of the genus Ancalomicrobium Staley 1968. Int. J. Syst. Evol. Microbiol. 68, 1955–1962. doi: 10.1099/ijsem.0.002772

PubMed Abstract | CrossRef Full Text | Google Scholar

Doetsch, R. N. (1981). “Determinative Methods of Light Microscopy,” in Manual of Methods for General Bacteriology, eds P. Gerdhardt, R. G. E. Murray, R. N. Costilow, E. W. Nester, W. A. Wood, N. R. Krieg, et al. (Washington, DC: American Society for Microbiology), 21–33.

Google Scholar

Eddy, S. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol. 7:1002195. doi: 10.1371/journal.pcbi.1002195

PubMed Abstract | CrossRef Full Text | Google Scholar

Ezaki, T., Hashimoto, Y., and Yabuuchi, E. (1989). Fluorometric deoxyribonucleic acid-deoxyribonucleic acid hybridization in microdilution wells as an alternative to membrane filter hybridization in which radioisotopes are used to determine genetic relatedness among bacterial strains. Int. J. Syst. Bacteriol. 39, 224–229. doi: 10.1099/00207713-39-3-224

CrossRef Full Text | Google Scholar

Felsenstein, J. (1981). Evolutionary trees from DNA sequences: a maximum likelihood approach. J. Mol. Evol. 17, 368–376. doi: 10.1007/bf01734359

PubMed Abstract | CrossRef Full Text | Google Scholar

Frank, J. A., Reich, C. I., Sharma, S., Weisbaum, J. S., Wilson, B. A., and Olsen, G. J. (2008). Critical evaluation of two primers commonly used for amplification of bacterial 16S rRNA genes. Appl. Environ. Microbiol. 74, 2461–2470. doi: 10.1128/aem.02272-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyatt, D., Chen, G. L., LoCascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi: 10.1186/1471-2105-11-119

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., and Goto, S. (2000). KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30.

Google Scholar

Katoh, K., and Standley, D. M. (2013). MAFFT Multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Khanal, A., and Lee, J. H. (2020). Functional diversity and abundance of nitrogen cycle-related genes in paddy soil. Appl. Biol. Chem. 63:17.

Google Scholar

Komagata, K., and Suzuki, K. (1988). 4 lipid and cell-wall analysis in bacterial systematics. Methods Microbiol. 19, 161–207. doi: 10.1016/s0580-9517(08)70410-0

CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, I., Chalita, M., Ha, S.-M., Na, S.-I., Yoon, S.-H., and Chun, J. (2017). ContEst16S: an algorithm that identifies contaminated prokaryotic genomes using 16S RNA gene sequences. Int. J. Syst. Evol. Microbiol. 67, 2053–2057. doi: 10.1099/ijsem.0.001872

PubMed Abstract | CrossRef Full Text | Google Scholar

Meier-Kolthoff, J. P., Auch, A. F., Klenk, H.-P., and Göker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics 14:60. doi: 10.1186/1471-2105-14-60

PubMed Abstract | CrossRef Full Text | Google Scholar

Minnikin, D. E., O’Donnell, A. G., Goodfellow, M., Alderson, G., Athalye, M., Schaal, A., et al. (1984). An integrated procedure for the extraction of bacterial isoprenoid quinones and polar lipids. J. Microbiol. Methods 2, 233–241. doi: 10.1016/0167-7012(84)90018-6

CrossRef Full Text | Google Scholar

Nguyen, T. M., and Kim, J. (2017). Azohydromonas riparia sp. nov. and Azohydromonas ureilytica sp. nov. isolated from a riverside soil in South Korea. J. Microbiol. 55, 330–336. doi: 10.1007/s12275-017-6519-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedrosa, F. O., Dobereiners, J., and Yates, M. G. (1980). Hydrogen-dependent growth and autotrophic carbon dioxide fixation in Derxia. J. Gen. Microbiol. 119, 547–551. doi: 10.1099/00221287-119-2-547

CrossRef Full Text | Google Scholar

Pereira e Silva, M. C., Schloter-Hai, B., Schloter, M., van Elsas, J. D., and Salles, J. F. (2013). Temporal dynamics of abundance and composition of nitrogen-fixing communities across agricultural Soils. PLoS One 8:e74500. doi: 10.1371/journal.pone.0074500

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2 - Approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490. doi: 10.1371/journal.pone.0009490

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruesse, E., Peplies, J., and Glöckner, F. O. (2012). SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28, 1823–1829. doi: 10.1093/bioinformatics/bts252

PubMed Abstract | CrossRef Full Text | Google Scholar

Pulschen, A. A., Bendia, A. G., Fricker, A. D., Pellizari, V. H., Galante, D., and Rodrigues, F. (2017). Isolation of uncultured bacteria from Antarctica using long incubation periods and low nutritional media. Front. Microbiol. 8:1346.

Google Scholar

Richter, M., and Rosselló-Móra, R. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. 106, 19126–19131. doi: 10.1073/pnas.0906412106

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425.

Google Scholar

Salehizadeh, H., Yan, N., and Farnood, R. (2020). Recent advances in microbial CO2 fixation and conversion to value-added products. Chem. Eng. J. 390:124584. doi: 10.1016/j.cej.2020.124584

CrossRef Full Text | Google Scholar

Sasser, M. (1990). Bacterial identification by gas chromatographic analysis of fatty acid methyl esters (GC-FAME). Newark: MIDI.

Google Scholar

Schöner, T. A., Gassel, S., Osawa, A., Tobias, N. J., Okuno, Y., Sakakibara, Y., et al. (2016). Aryl polyenes, a highly abundant class of bacterial natural products, are functionally related to antioxidative carotenoids. ChemBioChem 17, 247–253. doi: 10.1002/cbic.201500474

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, V., Misra, S., and Kumar Srivastava, A. (2017). Developing a green and sustainable process for enhanced PHB production by Azohydromonas australica. Biocatal. Agric. Biotechnol. 10, 122–129. doi: 10.1016/j.bcab.2017.02.014

CrossRef Full Text | Google Scholar

Shiner, E. K., Rumbaugh, K. P., and Williams, S. C. (2005). Interkingdom signaling: Deciphering the language of acyl homoserine lactones. FEMS Microbiol. Rev. 29, 935–947. doi: 10.1016/j.femsre.2005.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Smibert, R. M., and Krieg, N. R. (1994). “Phenotypic characterization,” in Methods for General and Molecular Bacteriology, eds P. Gerhardt, R. G. E. Murray, W. A. Wood, and N. R. Krieg (Washington DC: American Society for Microbiology), 607–654.

Google Scholar

Stamatakis, A. (2014). RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Tatusova, T., DiCuccio, M., Badretdin, A., Chetvernin, V., Nawrocki, E. P., Zaslavsky, L., et al. (2016). NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44, 6614–6624. doi: 10.1093/nar/gkw569

PubMed Abstract | CrossRef Full Text | Google Scholar

Vartoukian, S. R., Palmer, R. M., and Wade, W. G. (2010). Strategies for culture of “unculturable” bacteria. FEMS Microbiol. Lett. 309, 1–7.

Google Scholar

Wayne, L. G., Brenner, D. J., Colwell, R. R., Grimont, P. A. D., Kandler, O., Krichevsky, M. I., et al. (1987). Report of the ad hoc committee on reconciliation of approaches to bacterial systematics. Int. J. Syst. Evol. Microbiol. 37, 463–464. doi: 10.1099/00207713-37-4-463

CrossRef Full Text | Google Scholar

Xie, C. H., and Yokota, A. (2005). Reclassification of Alcaligenes latus strains IAM 12599T and IAM 12664 and Pseudomonas saccharophila as Azohydromonas lata gen. nov. comb. nov., Azohydromonas australica sp. nov. and Pelomonas saccharophila gen. nov., comb. nov., respectively. Int. J. Syst. Evol. Microbiol. 55, 2419–2425. doi: 10.1099/ijs.0.63733-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, H., Piao, C., Bian, D. R., Guo, M. W., and Li, Y. (2020). Azohydromonas aeria sp. nov., isolated from air. J. Microbiol. 58, 543–549. doi: 10.1007/s12275-020-9423-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, S. H., Ha, S. M., Kwon, S., Lim, J., Kim, Y., Seo, H., et al. (2017b). Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613–1617. doi: 10.1099/ijsem.0.001755

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, S.-H., Ha, S.-M., Lim, J., Kwon, S., and Chun, J. (2017a). A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek 110, 1281–1286. doi: 10.1007/s10482-017-0844-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zafar, M., Kumar, S., Kumar, S., Dhiman, A. K., and Park, H. S. (2014). Maintenance-energy-dependent dynamics of growth and poly(3-Hydroxybutyrate) [P(3HB)] production by Azohydromonas lata MTCC 2311 using simple and renewable carbon substrates. Brazil. J. Chem. Eng. 31, 313–323. doi: 10.1590/0104-6632.20140312s00002434

CrossRef Full Text | Google Scholar

Zhang, Z., Schwartz, S., Wagner, L., and Miller, W. (2000). A greedy algorithm for aligning DNA sequences. J. Comput. Biol. 7, 203–214. doi: 10.1089/10665270050081478

PubMed Abstract | CrossRef Full Text | Google Scholar