Pseudomonas oligotrophica sp. nov., a Novel Denitrifying Bacterium Possessing Nitrogen Removal Capability Under Low Carbon–Nitrogen Ratio Condition

Pseudomonas is a large and diverse genus within the Gammaproteobacteria known for its important ecological role in the environment. These bacteria exhibit versatile features of which the ability of heterotrophic nitrification and aerobic denitrification can be applied for nitrogen removal from the wastewater. A novel denitrifying bacterium, designated JM10B5aT, was isolated from the pond water for juvenile Litopenaeus vannamei. The phylogenetic, genomic, physiological, and biochemical analyses illustrated that strain JM10B5aT represented a novel species of the genus Pseudomonas, for which the name Pseudomonas oligotrophica sp. nov. was proposed. The effects of carbon sources and C/N ratios on denitrification performance of strain JM10B5aT were investigated. In addition, the results revealed that sodium acetate was selected as the optimum carbon source for denitrification of this strain. Besides, strain JM10B5aT could exhibit complete nitrate removal at the low C/N ratio of 3. Genomic analyses revealed that JM10B5aT possessed the functional genes including napA, narG, nirS, norB, and nosZ, which might participate in the complete denitrification process. Comparative genomic analyses indicated that many genes related to aggregation, utilization of alkylphosphonate and tricarballylate, biosynthesis of cofactors, and vitamins were contained in the genome of strain JM10B5aT. These genomic features were indicative of its adaption to various niches. Moreover, strain JM10B5aT harbored the complete operons required for the biosynthesis of vibrioferrin, a siderophore, which might be conducive to the high denitrification efficiency of denitrifying bacterium at low C/N ratio. Our findings demonstrated that the strain JM10B5aT could be a promising candidate for treating wastewater with a low C/N ratio.

Biological denitrification is an effective and general method for reducing nitrate to nitrogen gas due to its high efficiency, low cost, and environmental friendliness (Pang and Wang, 2021). Heterotrophic denitrification, which utilizes the organic carbon sources as electron donors, has a higher denitrification rate compared with autotrophic denitrification . However, the deficiency of available carbon source is an intractable problem for heterotrophic denitrification in the treatment of wastewater with low C/N ratio. To ensure sufficient electron donors for the denitrification process, the conventional solution is to add external organic carbon (often methanol, ethanol, glucose, and acetic acid) into the wastewater with low C/N ratio . Nevertheless, the addition of external organic carbon may result in the risk of secondary contamination and a high cost (Ling et al., 2021). Therefore, effective and sustainable strategies to enhance heterotrophic denitrification in the treatment of wastewater with a low C/N ratio are needed.
Accumulating approaches have been developed for the nitrate removal from the low C/N ratio wastewater, such as the system coupled with iron-based chemical reduction and autotrophic denitrification , the heterotrophic denitrification system amended with redox-active biochar , the fungal pellets immobilized bacterial bioreactor , and the utilization of biodegradable and inert carriers in the sequencing batch reactors . These approaches are focused on generating more electron donors to enhance nitrate removal capability or providing carriers for microorganisms to create a suitable environment. In addition, some heterotrophic denitrifying bacteria that can achieve efficient denitrification at low C/N ratio have been reported, such as Acinetobacter sp., Comamonas sp., and Pseudomonas sp. (Zhang S. et al., 2020;Chen et al., 2021;Fan et al., 2021). Compared with the abovementioned approaches, these denitrifying bacteria can not only effectively and readily remove nitrate from low C/N ratio wastewater without the utilization of complex engineering but also provide the functional microorganisms possessing excellent nitrate removal ability for the wastewater treatment systems. Among the denitrifying bacteria, Pseudomonas is the dominant bacterial genus in activated sludge and biofilm reactors (Deng et al., 2018;Zhang et al., 2019). Moreover, several Pseudomonas strains possess the denitrification ability and other functions, such as phosphorus removal ability and polyhydroxybutyrate-degrading ability (Di et al., 2019;Li et al., 2021). Therefore, Pseudomonas strains would outcompete other bacteria in the practical application.
In this study, a novel denitrifying bacterium, designated JM10B5a T , was isolated from the pond water for juvenile Litopenaeus vannamei. The phylogenetic, genomic, physiological, and biochemical analyses illustrated that strain JM10B5a T represented a novel species of the genus Pseudomonas, for which the name Pseudomonas oligotrophica sp. nov. was proposed. This strain performed excellent capability for denitrification under the low C/N ratio condition. Genomic information revealed that the functional genes of napA, narG, nirS, norB, and nosZ encoding the enzymatic repertoire for completely denitrification were identified in strain JM10B5a T . Our findings demonstrated that strain JM10B5a T could be a promising candidate for treating wastewater with a low C/N ratio.
P. stutzeri CGMCC 1.1803 T and P. balearica CCUG 44487 T were obtained from China General Microbiological Culture Collection Center (CGMCC) and Culture Collection University of Gothenburg (CCUG), respectively. These two type strains were used as the related strains for phenotypic, chemotaxonomic, and genetic analyses.

Phylogenetic Analysis
Genomic DNA of strains JM10B5a T was extracted from fresh cells, and the 16S rRNA gene sequence was amplified using the universal primers (27F/1492R) as described previously . Sequencing was performed by GENEWIZ, Inc. (Suzhou, China). The BLAST algorithm 2 and EzBioCloud database 3 were used to search for similar sequences (Yoon et al., 2017). Pairwise identities of 16S rRNA gene sequences were calculated using the software DNAMAN version 8. Multiple alignments of the 16S rRNA sequences were performed using the software MAFFT version 7.037 under the L-INS-i iterative refinement (Katoh and Standley, 2014). The phylogenetic tree was reconstructed using the software IQ-TREE version 2.1.2 with the maximum likelihood (ML) method under the TN+F+I+G4 nucleotide substitution model (Felsenstein, 1981;Kalyaanamoorthy et al., 2017;Minh et al., 2020). Support for the inferred ML tree was inferred by the ultrafast bootstrapping with 1,000 replicates (Diep Thi et al., 2018). The visualization and annotation of the resulting phylogenetic tree were performed using the software MEGA version X (Kumar et al., 2018).

Morphological Observations and Analyses of Physiological Characteristics
Cells and colonies of strain JM10B5a T cultivated on the nutrient broth agar medium (NA) at 30 • C for 48 h were observed by a transmission electron microscope (H7650, Hitachi) and naked eyes, respectively. Oxidase activity was determined using oxidase testing strips (HKM), and catalase activity was detected by bubble production after the addition of 3.0% H 2 O 2 (v/v) solution. Growth under an anaerobic environment was determined after 5 days of incubation on NA at 30 • C in an anaerobic pouch (MGC, Mitsubishi). Hydrolysis of casein, starch, tyrosine, and Tweens 20, 40, and 80 was investigated according to the protocols described by Lanyi (1987) and Tindall et al. (2007). Cellular motility was examined using the hanging-drop method (Bernardet et al., 2002).
Additional biochemical characteristics such as the enzyme activities, acid production by fermentation, and utilization of carbon sources were carried out using API ZYM and 20NE kits (bioMérieux) and Biolog GEN III MicroPlate (Biolog) according to the manufacturers' instructions. The API strips and Biolog results were recorded every 24 h after incubation at 30 • C until all reactions were steady. All API and Biolog tests were performed in duplicate under the consistent condition.

Analyses of Chemotaxonomic Characteristics
For fatty acid composition assay, the exponentially growing cells of strain JM10B5a T were harvested. Fatty acids were saponified, methylated, and extracted using the standard protocol of the Sherlock Microbial Identification System (MIDI). The prepared fatty acids were analyzed by gas chromatography (model 7890A; Agilent) using the Microbial Identification software package with the Sherlock MIDI 6.1 system and the Sherlock Aerobic Bacterial Database (TSBA 6.1). Polar lipids were extracted using the chloroform/methanol system and separated by two-dimensional thin-layer chromatography (TLC, Silica gel 60 F254, Merck) (Minnikin et al., 1984). The plates dotted with samples were subjected to two-dimensional development with the first solvent of chloroform-methanol-water (65:25:4, v/v) and the second solvent of chloroform-methanol-acetic acid-water (80:12:15:4, v/v). Total lipids and specific functional groups were detected using the different spray staining reagents on separate TLC plates: 10% ethanolic molybdophosphoric acid, ninhydrin, molybdenum blue, and α-naphthol-sulfuric acid. The quinone was extracted from the freeze-dried cells and determined using high-performance liquid chromatography (Agilent 1200; ODS 250 × 4.6 mm × 5.0 µm; flowing phase, methanol-isopropanol, 2:1; 1.0 ml/min) according to the methods described by Collins and Jones (1981).

Genome Sequencing and Function Analysis
The genomic DNA of strain JM10B5a T was extracted using the HiPure Bacterial DNA Kit (Magen Biotech, Guangzhou) according to the manufacturer's instruction. The genome was sequenced using the Illumina NovaSeq PE150 platform in Shanghai Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Raw reads were filtered and then de novo assembled using the software SPAdes version 3.14.1 under "-careful" mode and with a k-mer value of 127 (Bankevich et al., 2012). Available genomes of the related type strains P. stutzeri CGMCC 1.1803 T and P. balearica CCUG 44487 T were obtained from the GenBank database with the numbers CP002881 and CP007511, respectively. The genome of strain JM10B5a T was in draft status and was quality-checked using the software CheckM version 1.1.2 (Parks et al., 2015). Overall genome relatedness indices (OGRIs) including digital DNA-DNA hybridization values (dDDH) and average nucleotide identity values (ANI) were estimated using the genome-to-genome distance calculator version 2.1 online service with the recommended formula 2 4 and the software FastANI version 1.31, respectively (Meier-Kolthoff et al., 2013;Jain et al., 2018). The UBCG pipeline which could provide an accurate phylogenomic tree for the taxonomic purpose was used Na et al., 2018). The genomes of strain JM10B5a T and the related type strains were annotated using the software Prokka version 1.13 (Seemann, 2014) and Rapid using Subsystem Technology (RAST) version 2.0 with the default parameters (Overbeek et al., 2014).

Assessment of Nitrogen Removal Characteristics
The results of the API 20NE test preliminarily indicated that strain JM10B5a T was capable of reducing nitrate and nitrite. Furthermore, the genome annotation revealed that strain JM10B5a T possessed the functional genes including napA, narG, nirS, norB, and nosZ that participated in the denitrification process. These results suggested that strain JM10B5a T had the potential capability of complete denitrification. To investigate the optimum carbon source of denitrification for strain JM10B5a T , the carbon source was replaced in the DSM-1 (DSM with KNO 3 0.36 g/L as the sole nitrogen source) or DSM-2 (DSM with NaNO 2 0.25 g/L as the sole nitrogen source) by sodium acetate, sodium succinate, sodium citrate, glucose, sucrose, and starch, respectively, with the C/N ratio of 10 (the molar mass of carbon to the molar mass of nitrogen). Furthermore, the initial C/N ratios of 2, 3, 4, and 5 were controlled by changing the addition amount of the optimum carbon source in the DSM-1. Before the experiments, strain JM10B5a T was incubated in NB at 30 • C and 180 rpm for 24 h. The suspension was centrifuged and then washed 3 times with sterile physiological saline to remove the residual medium. The biomass was resuspended in sterile water with the initial OD 600 adjusting to 1.0. The bacterial suspension was inoculated into a conical flask containing sterile medium with the inoculum amount of 4.0% and incubated at 30 • C statically. Samples taken from flasks after 48 h of incubation were used for determining the OD 600 and chemical analyses. The medium without bacterial inoculation was used as control treatment. All the experiments were performed in biological quadruplicate.

Analytical Methods and Statistical Analyses
The concentrations of NO − 2 -N and NO − 3 -N were measured using N-(1-naphthyl) ethylenediamine dihydrochloride and UV spectrophotometry, respectively, according to the standard analytical procedures (Association et al., 2005). All the data were presented as mean and standard error (SE). One-way ANOVA (multiple range test) was performed with Tukey's HSD test (P < 0.05) using SPSS Statistics 19.

16S rRNA Gene Sequence Analyses and Chemotaxonomic Characterization
The almost complete 16S rRNA gene sequence of strain JM10B5a T obtained by amplification (1,375 bp) was included in the complete 16S rRNA gene sequence assembled from genomic sequences (1,537 bp). The sequence comparison showed that strain JM10B5a T fell into the genus Pseudomonas and shared the highest similarity with P. stutzeri CGMCC 1.1803 T (98.0%). All the other type strains showed similarities lower than 98.0% with strain JM10B5a T . Based on phylogenetic analysis using the ML algorithm, strain JM10B5a T was stably located in the genus Pseudomonas and formed a clade with P. balearica CCUG 44487 T (the similarity to JM10B5a T was 97.0%) at the 79.0% bootstrap confidence level (Figure 1). The similarity and phylogenetic analysis based on the 16S rRNA gene sequences indicated that strain JM10B5a T should represent a novel member of the genus Pseudomonas. Therefore, the type strains P. stutzeri CGMCC 1.1803 T and P. balearica CCUG 44487 T were purchased from the culture collection centers and used as references for further comparisons of phenotypic and chemotaxonomic characteristics.
The predominant cellular fatty acids (>10%) of strain JM10B5a T were C 16 : 0 (22.0%), summed feature 3 (C 16 : 1 ω6c FIGURE 1 | ML tree based on 16S rRNA gene sequences of strain JM10B5a T and the related type strains of genus Pseudomonas. Type strain Entomomonas moraniae QZS01 T was used as an outgroup. There were a total of 1,303 positions used to construct the phylogenetic tree. Bootstrap values higher than 60% were shown at branch points. Bar, 0.05 represents the number of substitutions per site. and/or C 16 : 1 ω7c, 21.3%), and summed feature 8 (C 18 : 1 ω6c and/or C 18 : 1 ω7c, 29.1%), which were also detected in other Pseudomonas species (Mamtimin et al., 2021). Comparative fatty acid profiles between strain JM10B5a T and the two related type strains are shown in Table 1. The fatty acid profile of JM10B5a T was similar to that of other related strains, although there were minor quantitative differences observed. The predominant respiratory quinone of strain JM10B5a T was ubiquinone-9 (Q-9) which was consistent with that of other members of the genus Pseudomonas (Zou et al., 2019;Mamtimin et al., 2021). The major polar lipids were phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and diphosphatidylglycerol (DPG), which were consistent with the previously published data for Pseudomonas species (Zou et al., 2019;Li et al., 2020). In addition, minor amounts of two unidentified aminophospholipids (APLs) and three unidentified phospholipids (PLs) were also present (Supplementary Figure 1).

Morphological Observations and Physiological Characterization
The morphological features of strain JM10B5a T were studied on NA medium and formed irregular, filmy sheet, non-transparent, and ecru white colonies after 48 h of incubation at 30 • C (Supplementary Figure 2A). Cells of strain JM10B5a T were observed to be gram-stain-negative, aerobic, rod-shaped (0.6-0.7 × 1.4-1.5 µm), facultative anaerobic, and motile with a single polar flagellum (Supplementary Figure 2B). Growth was determined at 10-45 • C (optimum, 25-30 • C), at pH 5.5-11.0 (optimum, 6.0), and in 0-4.0% (w/v) NaCl (optimum, 0.5-2.0%). The activities of catalase and oxidase were positive, and the further details differentiating strain JM10B5a T from the two related type strains are shown in Table 2. For instance, strain JM10B5a T could not grow in the NB with a NaCl concentration of 5.0 while the two related type strains could grow; the optimum pH of strain JM10B5a T was 6.0, which was inconsistent with that of P. stutzeri CGMCC 1.1803 T and P. balearica CCUG 44487 T of 6.5-7.5 and 8.5-9.5, respectively; the growth of strain JM10B5a T was weaker than that of the two related type strains at 45 • C. Besides, stain JM10B5a T could be distinguished from the reference type strains by positive for adipic acid, trisodium citrate, esterase (C4), and lipase (C14) in the API 20NE and ZYM tests; and positive for Dgalacturonic acid, D-glucuronic acid, gentiobiose, L-arginine, L-aspartic acid, L-galactonic acid lactone, mucic acid, quinic acid, α-hydroxybutyric acid, and γ -aminobutyric acid in the Biolog GNE III MicroPlate system. Importantly, strain JM10B5a T  exhibited the capability for complete denitrification that liberated copious amounts of nitrogen gas from nitrate. The denitrification characteristic of this strain was consistent with that of P. stutzeri and P. balearica strains (Bennasar et al., 1996;Li et al., 2021).

Genome Data and Comparative Genomic Analysis
The draft genome of strain JM10B5a T was assembled into 20 contigs (>500 bp) with a genome size of 3,978,222 bp, and the estimated completeness and contamination were 100 and 0.1%, respectively, indicating that this genome was high quality according to the standards (Parks et al., 2015). Genomes of strains P. stutzeri CGMCC 1.1803 T and P. balearica CCUG 44487 T were obtained from the GenBank database with the numbers CP002881 and CP007511, respectively. Based on the Prokka annotation, a total of 3,744 genes, 3,679 protein-coding sequences (CDS), 60 tRNA genes, and 3 rRNA genes were found in the genome of strain JM10B5a T . The genomic DNA G + C content of JM10B5a T was 67.2%. Compared with the genomes of related type strains, JM10B5a T had ANI values of 78.9-88.1% and dDDH values of 20.8-30.8% (Supplementary Table 1), which were all below 95.0 and 70% cutoff commonly used to define a bacterial species, respectively (Stackebrandt and Goebel, 1994;Richter and Rossello-Mora, 2009). To further determine the taxonomic position of strain JM10B5a T , the phylogenomic tree based on the 92 bacterial core gene sets was reconstructed (Na et al., 2018). As shown in Figure 2, strain JM10B5a T and P. balearica CCUG 44487 T formed a clade with 100% bootstrap value, which was consistent with the phylogenetic tree based on the 16S rRNA gene sequences. Therefore, the comprehensive analyses of OGRIs and phylogenomic tree further indicated that strain JM10B5a T should represent a novel species within the genus of Pseudomonas.
According to the subsystem category distribution of RAST annotation, there were 190, 128, 402, and 336 genes associated with "RNA metabolism, " "DNA metabolism, " "amino acid and derivatives, " and "carbohydrates, " respectively, which might be important to the bacterial growth. In addition, 75 genes were identified in the nitrogen metabolism, including the functional genes related to denitrification and assimilation of nitrite and nitrate (Figures 3A,B). Nitrate reductase catalyzing the reduction in nitrate to nitrite plays an important role in the nitrogen cycle (Kuypers et al., 2018). Previous studies have shown that bacterial dissimilatory nitrate reduction could be catalyzed by two different enzymes: a membrane-bound nitrate reductase (NAR, the catalytic subunit NarG encoded by narG) and periplasmic nitrate reductase (NAP, the catalytic subunit NapA encoded by napA) (Moreno-Vivian et al., 1999). Genes narG and napA were both contained in strain JM10B5a T , which was consistent with Pseudomonas sp. JQ-H3 and Paracoccus denitrificans (Moreno-Vivian et al., 1999;Wang et al., 2019). The dissimilatory reduction in nitrite to nitric oxide can be catalyzed by two unrelated enzymes: a cytochrome cd1 nitrite reductase (cd1-NIR, encoded by nirS) or a Cu-containing nitrite reductase (Cu-NIR, encoded by nirK), which were widespread among bacteria and archaea (Kuypers et al., 2018). Gene nirS rather than nirK was identified in strain JM10B5a T , which was consistent with Pseudomonas sp. JQ-H3, P. bauzanensis DN13-1, and P. stutzeri T13 Zhang M. et al., 2020;Feng et al., 2021). In addition, genes norB encoding for nitric oxide reductase (NorB, the key enzyme for reducing nitric oxide to nitrous oxide) and nosZ encoding for nitrous oxide reductase (NosZ, the key enzyme for reducing nitrous oxide to nitrogen gas) were also identified. Therefore, it was speculated that strain JM10B5a T could perform the complete denitrification pathway: Hence, strain JM10B5a T was a novel denitrifying bacterium.
In addition, comparative genomic analyses showed that strain JM10B5a T encoded 126 and 147 functional genes more than P. balearica CCUG 44487 T and P. stutzeri CGMCC 1.1803 T , respectively, and 67 genes of them were the shared differential genes (Supplementary Table 2). These genes were related to type I secretion system for aggregation, utilization of alkylphosphonate and tricarballylate, histidine ABC transporter, LysR family transcriptional regulator, and many other proteins. Most of these proteins were mainly affiliated with dehydratase, transferase, and biosynthesis of cofactors and vitamins. The genes encoding arginase and galactarate dehydratase made strain JM10B5a T performing the positive for L-arginine and mucic acid ( Table 2). Moreover, strain JM10B5a T was found to harbor the complete operon pvsABCDE and vibrioferrin receptor gene pvuA ( Figure 3C) that are required for the production of a siderophore vibrioferrin. This siderophore was shown to be responsible for vibrioferrin-mediated iron uptake in the terrestrial bacteria Azotobacter vinelandii, Pseudomonas sp., and the marine bacterium Vibrio parahaemolyticus (Yamamoto et al., 1994;Baars et al., 2016;Stanborough et al., 2018;Sood et al., 2019). It was reported that siderophore had a very high affinity to iron with the formation of Fe-siderophore which could facilitate the transport of iron through the membrane directly (Stintzi et al., 2000). Moreover, the denitrification efficiency might be enhanced with the improvement of iron transport from extracellular to intracellular . Therefore, the biosynthesis of vibrioferrin might partly explain that the strain JM10B5a T performed the high denitrification efficiency under the low C/N ratio condition.

Effect of Carbon Source on Denitrification by Strain JM10B5a T
Heterotrophic denitrifying bacteria require organic carbon for cell growth and as the electron donor in the denitrification process (Rajta et al., 2020). Thus, carbon source was considered to be an important factor influencing denitrification. As expected, significant differences were observed using different carbon sources. As shown in Figure 4A, the strain JM10B5a T could grow well with the OD 600 of 0.42 when glucose was used as the sole carbon source. However, it exhibited a higher NO − 3 -N removal efficiency of 99.5% when sodium acetate served as carbon source than that of 58.6% when glucose served as carbon source. As shown in Figure 4B, strain JM10B5a T could grow quite well when glucose served as the carbon source (OD 600 of 0.52), but it performed the relatively high NO − 2 -N removal of 100% when sodium acetate served as the carbon source. These results indicated that sodium acetate was the optimum carbon source for the performance of denitrification for strain JM10B5a T , which was consistent with Bacillus pumilus, Arthrobacter sp., and Streptomyces lusitanus (Elkarrach et al., 2021). Therefore, sodium acetate was employed in the following experiments.

Effect of Carbon/Nitrogen Ratio on Nitrate
Removal by Strain JM10B5a T C/N ratio is a measure of the electron donor to acceptor ratio in the biological denitrification process. This factor can influence the denitrification efficiency and accumulation of NO − 2 -N. To investigate the influence of C/N ratio on denitrification of strain JM10B5a T , the initial concentration of NO − 3 -N was fixed at ∼50.0 mg/L and the different initial C/N ratios (2, 3, 4, and 5) were controlled. As shown in Figure 5, strain JM10B5a T grew poorly with the OD 600 range of 0.10-0.15 at these low C/N ratios. The NO − 3 -N was reduced completely at the C/N ratios of 3-5 which was distinctly higher than that of 84.4% at the C/N ratio of 2. Furthermore, there was no NO − 2 -N accumulation observed at the C/N ratios of 3-5, but the NO − 2 -N concentration of 19.3 mg/L was detected at the C/N ratio of 2. The low NO − 3 -N removal efficiency and NO − 2 -N accumulation at the C/N ratio of 2 were mainly due to the limited carbon source that could not provide sufficient energy for bacterial growth and electron donors for denitrification (Fan et al., 2021). These results were consistent with previous reports that denitrification efficiency could decrease under the extremely low carbon concentration . Numerous studies have suggested that the optimum C/N ratios for most heterotrophic denitrifying bacteria were in the range of 8-15 (Guo et al., 2016;Liu et al., 2018;Zhao et al., 2018). For example, the C/N ratios of 10 and 15 were the most suitable for strains P. stutzeri XL-2 and P. taiwanensis J to achieve efficient nitrate removal, respectively (He et al., 2018;Zhao et al., 2018). In addition, some denitrifying bacteria could achieve complete denitrification at the relatively low C/N ratios. For example, Rout et al. (2017) used Bacillus cereus to remove nitrogen from domestic wastewater with 7.5 as the optimum C/N ratio; Zhang S. et al. (2020) reported that strain Comamonas sp. YSF15 could achieve complete denitrification at C/N of 3. Our results showed that strain JM10B5a T could achieve the complete nitrate removal without accumulation of nitrite at the low carbon-nitrogen ratio of 3 (COD Cr /TN ratio of 2.6), indicating that this strain could be a promising candidate for treating the oligotrophic wastewater.  Description of Pseudomonas oligotrophica sp. nov.

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
MZ, AL, and QY were involved in conceptualization and project administration. MZ involved in data curation, software, visualization, and writing-original draft. MZ and BX designed the formal analysis and methodology. AL, HZ, and BX were involved in funding acquisition. MZ, AL, QY, BX, and HZ investigated the study. HZ collected the resources and supervised the study. MZ and AL validated the study. MZ, QY, and BX were involved in writing-review and editing. All authors contributed to the article and approved the submitted version.