NitrogenFixing Capability
Since bacteria in the soil sample were cultured in nitrogen-free medium on the purpose of isolating nitrogen-fixing strain, strain 39 is possible to have nitrogen-fixing capability. Strains HN-1T isolated from the rhizosphere of rice was detected by acetylene reduction to have nitrogen-fixing capacity (Liu et al. 2019). To confirm genes encoding nitrogenase are contained in their genome, nifH gene was amplified and sequenced. However, the sequencing result of nifH fragment both displayed double peaks, which indicated that there were multiple nifH genes in their genome. Acetylene reduction assays were performed to verify the nitrogenase activity of HN-1T and 39. As shown in Table 1, strains HN-1T and 39 exhibited very high nitrogenase activity compared to other nitrogen-fixing Paenibacillus species, suggesting a high efficiency of the nitrogen fixation process.
Table 1
Nitrogenase activity of strains HN-1T and 39 in comparison with some nitrogen-fixing species of the genus Paenibacillus
Strain
|
Nitrogenase activity
[nmol C2H4 (mg protein h)−1]
|
P. polymyxa DSM 36T
|
1355.1 ± 152.4
|
P. stellifer DSM 14472T
|
6099.5 ± 497.3
|
P. zanthoxyli JH29T
|
6282.4 ± 307.7
|
P. graminis RSA19T
|
4272.9 ± 207.9
|
P. sabinae DSM 17841T
|
7749 ± 371.8
|
P. durus ATCC 35681T
|
5511.5 ± 260.2
|
P. triticisoli BJ-18T
|
743.6 ± 82.9
|
39
|
7160.3 ± 584.6
|
HN-1T
|
6937.2 ± 625.1
|
Molecular Characterization
The newly generated 16S rRNA gene sequence of strain 39 were compared with sequences of bacterial type strains in the Genbank nucleotide database. Phylogenetic trees were inferred using the neighbour-joining and maximum-parsimony methods in the software MEGA7. As no significant differences were found among the phylogenetic trees obtained by the different methods used, only the trees constructed by using the neighbour-joining method are shown. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strains HN-1T and 39 clustered with species of the genus Paenibacillus and formed a monophyletic cluster with P. stellifer DSM 14472T, as the three strains formed a separate phylogenetic branch within the genus Paenibacillus with a high bootstrap resampling value of 100 % (Fig. 1). The strain 39 shares 100% 16S rRNA gene sequence identity with strain HN-1T. The strains HN-1T or 39 showed the highest 16S rRNA gene sequence identity to P. stellifer DSM 14472T (99.5%), followed by P. durus ATCC 35681T (97.4%) and P. sabinae DSM 17841T (97.0%).
Generally, 98.7% sequence identity on the 16S rRNA gene are considered to be within the same species (Kim et al. 2014). However, several reports have been published showing that Paenibacillus species with > 99 % 16S rRNA gene sequence similarity may not belong to the same species (Kamfer et al. 2017; Kim and Cha 2018; Ghio et al. 2019; Guella et al. 2019; Velazquez et al. 2020). Thus, housekeeping genes are now routinely used to complement the 16S rRNA gene analysis for species level determination (da Mota et al. 2004; Holmes et al. 2004; Rodriguez et al. 2019). Due to the low level of discrimination based on 16S rRNA gene between closely related species, the gyrB gene (coding for the b subunit of DNA gyrase) was used as an alternative phylogenetic marker (Wang et al. 2007). The gyrB genes were retrieved from the HN-1T and 39 genomes. The gyrB gene clearly distinguishes HN-1T and 39 from other Paenibacillus species with only 93.04% gene sequence identity to P. stellifer DSM 14472 (Fig. S1). Based on the 95–96% gyrB gene sequence similarity as the interspecies gap (Lee et al. 2008; Liu et al. 2013), strains HN-1T and 39 could be assigned to novel species.
Phenotypic Characteristics
Strains HN-1T and 39 were found to be Gram-positive, facultatively anaerobic, motile and rod-shaped. Colonies of strain HN-1T and 39 grown on LD medium after 72 h of incubation at 30℃ were usually 0.8–1.2 mm in diameter, circular, moist, milky and convex. The transmission electron micrographs of type strain HN-1T showed the presence of peritrichous flagella on cell surface (Fig. 2a). Strain HN-1T produced ellipsoidal spores in swollen sporangia in the terminal region of the cell by scanning electron microscope (Fig. 2b).
In order to determine physiological and biochemical characteristics of HN-1T and 39 in comparison with P. stellifer DSM 14472T and P. sabinae DSM 17841T, a series of tests were carried out following the proposed minimal standards for describing new taxa of facultatively anaerobic, endospore-forming bacteria (Logan et al. 2009). The ability of strains to assimilate different substrates were tested using GEN III microplates by Biolog system (Biolog Microstation TM, CA, USA) (Kiran et al. 2017; Ripa et al. 2019).
The strains HN-1T and 39 grew well in up to 4% NaCl (w/v), however, strain P. stellifer DSM 14472T tolerated only 3% NaCl. The pH range for growth was 5.0–9.0 and the temperature range for growth is 15–42°C. Strains HN-1T and 39 was determined to be negative for the Voges–Proskauer reaction, and positive for the methyl red reaction. Strains HN-1T and 39 were positive for catalase reaction and can produce acid from rhamnose and sorbitol, which differentiated HN-1T and 39 from the most related P. stellifer DSM 14472T. The physiological and biochemical characteristics of HN-1T and 39 are shown in Table 2. In Biolog analysis, significant differences were observed in the metabolization of different substrates on Biolog GEN III microplate by strains HN-1T and 39 from its closest type strain P. stellifer DSM 14472T (Table 3). Strain HN-1T and P. stellifer DSM 14472T differed in the metabolization of D-Fucose, D-Maltose, 3-Methyl glucose, D-Sorbitol, Stachyose, Citric acid, α-Keto-butyric acid, Mucic acid, Methyl pyruvate, Gelatin, Inosine, D-Glucose-6-PO4, Pectin, Aztreonam, Fusidic acid, Nalidixic acid, Vancomycin, Lithium chloride, Sodium bromate, Sodium lactate 1%, Rifamycin sv and Troleandomycin as a sole carbon source. Compared with type strains of closely related Paenibacillus species, strain HN-1T and 39 exhibited nearly identical phenotypic characteristics.
Table 2
Differential phenotypic characteristics between strain HN-1T, 39 and closely related strains.
Characteristic
|
HN-1T
|
39
|
DSM14472T
|
DSM 17841T
|
Growth conditions
|
|
Temperature (℃)
|
15–42
|
15–42
|
15–42
|
15–42
|
pH
|
5–9
|
5–9
|
5–9
|
5–9
|
NaCl
|
0–4%
|
0–4%
|
0–3%
|
0–3%
|
0.001% lysozyme
|
-
|
-
|
-
|
-
|
Oxidase
|
-
|
-
|
-
|
-
|
Catalase
|
+
|
+
|
-
|
+
|
Voges-Proskauer
|
-
|
-
|
-
|
+
|
Methyl red reaction
|
+
|
+
|
+
|
+
|
Nitrate reduction
|
-
|
-
|
-
|
+
|
Hydrolysis of:
|
Starch
|
+
|
+
|
+
|
-
|
Tween 20
|
-
|
-
|
-
|
-
|
Tyrosine
|
-
|
-
|
-
|
-
|
Esculin
|
+
|
+
|
+
|
+
|
Production of:
|
-
|
-
|
-
|
-
|
Indole
|
|
|
|
|
Dihydroxyacetone
|
-
|
-
|
-
|
-
|
Dextrin
|
+
|
+
|
+
|
-
|
Production of acid from:
|
Glucose
|
+
|
+
|
+
|
+
|
Fructose
|
+
|
+
|
+
|
+
|
Rhamnose
|
+
|
+
|
-
|
-
|
Glycerol
|
-
|
-
|
-
|
-
|
Trehalose
|
+
|
+
|
+
|
+
|
Arabinose
|
+
|
+
|
+
|
-
|
Mannitol
|
-
|
-
|
-
|
+
|
Xylose
|
+
|
+
|
+
|
-
|
Sorbitol
|
+
|
+
|
-
|
+
|
Maltose
|
+
|
+
|
+
|
+
|
Starch
|
+
|
+
|
+
|
+
|
Strains: 1. HN-1T; 2. 39; 3. P. stellifer DSM 14472T; 4, P. sabinae DSM 17841T; +. Positive reaction; 2. negative reaction. |
Table 3
Phenotypic fingerprinting of 1. strain HN-1T; 2. strain 39; 3. P. stellifer DSM 14472T; 4. P. sabinae DSM 17841T by oxidation of sole substrates on Biolog GEN III MicroPlates (All data from present study generated under the same conditions).
substrate
|
1
|
2
|
3
|
4
|
|
1
|
2
|
3
|
4
|
Carbohydrates:
|
|
|
|
|
Esters:
|
|
|
|
|
n-Acetyl-β-D-mannosamine
|
-
|
-
|
-
|
-
|
D-Lactic acid methyl ester
|
-
|
-
|
-
|
-
|
n-Acetyl-D-galactosamine
|
-
|
-
|
-
|
-
|
Methyl pyruvate
|
w
|
w
|
-
|
+
|
n-Acetyl-D-glucosamine
|
-
|
-
|
-
|
-
|
Alcohols:
|
|
|
|
|
n-Acetyl neuraminic acid
|
-
|
-
|
-
|
-
|
Glycerol
|
-
|
-
|
-
|
-
|
D-Arabitol
|
-
|
-
|
-
|
-
|
Amides:
|
|
|
|
|
α-D-Glucose
|
+
|
+
|
+
|
+
|
Glucuronamide
|
w
|
w
|
w
|
-
|
α-D-Lactose
|
+
|
+
|
+
|
-
|
Amino acids:
|
|
|
|
|
D-Cellobiose
|
+
|
+
|
+
|
+
|
L-Alanine
|
-
|
-
|
-
|
-
|
D-Fructose
|
+
|
+
|
+
|
+
|
L-Arginine
|
-
|
-
|
-
|
-
|
D-Fucose
|
-
|
-
|
w
|
-
|
L-Aspartic acid
|
-
|
-
|
-
|
-
|
D-Galactose
|
+
|
+
|
+
|
+
|
D-Aspartic acid
|
-
|
-
|
-
|
-
|
D-Gentiobiose
|
+
|
+
|
+
|
+
|
Gelatin
|
w
|
w
|
-
|
-
|
β-Methyl-D-glucoside
|
+
|
+
|
+
|
+
|
Glycyl-l-proline
|
-
|
-
|
-
|
-
|
D-Maltose
|
+
|
+
|
+
|
+
|
L-Glutamic acid
|
-
|
-
|
-
|
-
|
D-Mannitol
|
-
|
-
|
-
|
+
|
L-Histidine
|
-
|
-
|
-
|
-
|
D-Mannose
|
+
|
+
|
w
|
+
|
L-Pyroglutamic acid
|
-
|
-
|
-
|
-
|
D-Melibiose
|
w
|
w
|
+
|
+
|
L-Serine
|
-
|
-
|
-
|
-
|
3-Methyl glucose
|
+
|
+
|
w
|
+
|
D-Serine
|
-
|
-
|
-
|
-
|
D-Raffinose
|
+
|
+
|
+
|
+
|
Nucleosides
|
|
|
|
|
D-Salicin
|
+
|
+
|
+
|
+
|
Inosine
|
w
|
w
|
-
|
-
|
D-Sorbitol
|
+
|
+
|
-
|
+
|
Phosphorylated compounds:
|
|
|
|
|
D-Trehalose
|
+
|
+
|
+
|
+
|
D-Fructose-6-PO4
|
w
|
w
|
w
|
-
|
D-Turanose
|
+
|
+
|
+
|
+
|
D-Glucose-6-PO4
|
w
|
w
|
-
|
-
|
L-Fucose
|
w
|
w
|
w
|
-
|
Polymers:
|
|
|
|
|
L-Rhamnose
|
w
|
w
|
w
|
-
|
Dextrin
|
+
|
+
|
+
|
+
|
Stachyose
|
+
|
+
|
w
|
+
|
Pectin
|
+
|
+
|
w
|
+
|
Sucrose
|
+
|
+
|
+
|
+
|
Tween 40
|
-
|
-
|
-
|
-
|
Carboxylic acid:
|
|
|
|
|
Antibiotics:
|
|
|
|
|
Acetic acid
|
-
|
-
|
-
|
-
|
Aztreonam
|
+
|
+
|
w
|
+
|
Acetoacetic acid
|
w
|
w
|
w
|
+
|
Fusidic acid
|
w
|
w
|
-
|
+
|
γ-Amino-butryric acid
|
-
|
-
|
-
|
-
|
Lincomycin
|
-
|
-
|
-
|
+
|
Bromo-succinic acid
|
-
|
-
|
-
|
-
|
Minocycline
|
-
|
-
|
-
|
+
|
Citric acid
|
w
|
w
|
-
|
-
|
Nalidixic acid
|
+
|
+
|
-
|
+
|
Formic acid
|
-
|
-
|
-
|
-
|
Rifamycin sv
|
-
|
-
|
+
|
+
|
D-Galacturonic acid
|
w
|
w
|
w
|
-
|
Troleandomycin
|
-
|
-
|
-
|
+
|
D-Gluconic acid
|
-
|
-
|
-
|
+
|
Vancomycin
|
+
|
+
|
w
|
+
|
D-Glucuronic acid
|
w
|
w
|
w
|
-
|
pH 5
|
w
|
w
|
w
|
w
|
L-Galactonic acid lactone
|
w
|
w
|
w
|
-
|
pH 6
|
+
|
+
|
+
|
+
|
α-Hydroxy-butyric acid
|
-
|
-
|
-
|
-
|
NaCl 4%
|
+
|
+
|
-
|
w
|
β-Hydroxy-D, L butyric acid
|
-
|
-
|
-
|
-
|
NaCl 8%
|
-
|
-
|
-
|
-
|
p-Hydroxy-phenylacetic acid
|
-
|
-
|
-
|
-
|
Guanidine HCl
|
-
|
-
|
-
|
-
|
α-Keto-butyric acid
|
w
|
w
|
-
|
w
|
Lithium chloride
|
+
|
+
|
w
|
+
|
α-Keto-glutaric acid
|
w
|
w
|
w
|
-
|
Myo-inositol
|
-
|
-
|
-
|
-
|
L-Lactic acid
|
-
|
-
|
-
|
-
|
Niaproof 4
|
-
|
-
|
-
|
-
|
D-Malic acid
|
-
|
-
|
-
|
-
|
Potassium tellurite
|
+
|
+
|
+
|
+
|
L-Malic acid
|
-
|
-
|
-
|
-
|
Sodium bromate
|
+
|
+
|
-
|
w
|
Mucic acid
|
w
|
w
|
-
|
-
|
Sodium butyrate
|
+
|
+
|
+
|
+
|
Propionic acid
|
-
|
-
|
-
|
-
|
Sodium lactate 1%
|
+
|
+
|
w
|
+
|
Quinic acid
|
-
|
-
|
-
|
-
|
Tetrazolium blue
|
-
|
-
|
w
|
w
|
D-Saccharic acid
|
-
|
-
|
-
|
-
|
Tetrazolium violet
|
w
|
w
|
+
|
w
|
+, Positive; −, negative; w, weak reaction |
Chemotaxonomic Characteristics
In order to determine the composition of cellular fatty acid, HN-1T and type strains P. stellifer DSM 14472T and P. sabinae DSM 17841T were incubated in LD medium at 30°C for 2 days. Whole cell fatty acid analysis revealed that anteiso-C15:0, C14: 0, C16: 0, iso-C14: 0, iso-C16: 0 and iso-C15: 0 are present as major (> 5 %) fatty acids, and anteiso-C15: 0, anteiso-C17: 0, iso-C17: 0 and C18:1ω9c are present as minor (< 5 but > 1 %) fatty acids (Table S1). Anteiso-C15:0 is the predominant fatty acid of members of the genus Paenibacillus (Ash et al. 1993), consistent with strain HN-1T (41.14%) being a member of this genus. However, in the closely related type strains P. stellifer DSM 14472T, the fatty acid C16:0 was found to be more abundant than anteiso-C15:0. The major menaquinone of HN-1T was MK-7, in conformity to genus Paenibacillus. The polar lipids detected by two-dimensional TLC are diphosphatidylglycerol (DPG), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), four aminophospholipids (APL) and unidentified glycolipid (Fig. S1), in agreement with the profile of genus Paenibacillus.
Genome sequence and similarity analysis
Genome sequencing was performed to evaluate the genomic relatedness of the strains HN-1T and 39 to its closely related recognized species in the genus Paenibacillus. Genomes of strains HN-1T and 39 were approximately 6.32 and 6.45 Mbp, respectively. The DNA G + C content of the strains HN-1T and 39 were 53.36 and 52.99 mol%, respectively, which is within the range 39–59 mol% reported for the members of genus Paenibacillus (Im et al. 2017). The total number of protein coding genes in HN-1T and 39 were 5631 and 5782, respectively. While, the related strain P. stellifer DSM 14472T had a complete genome of 5.65Mb, comprising 4958 protein coding genes with a DNA G + C content of 55.6 mol% (Genbank: CP009286). More detailed genomic information of strains HN-1T and 39 were shown in Supplementary Table S2 and high-quality draft genomes of strains HN-1T and 39 were deposited in GenBank under accession numbers JAHCMB000000000 and JAHBAZ000000000, respectively.
In this study, digital DNA–DNA hybridization (dDDH) values were calculated using GGDC. The dDDH value of genomes for strain HN-1T and strain 39 was 99.8%. The level of DNA–DNA relatedness was determined to be 42.1% by strain HN-1T with P. stellifer DSM 14472T, which is lower than the threshold value of 70 % for species delineation (Chun et al. 2018).
The average nucleotide identity (ANI) is widely used to define bacterial species (Konstantinidis and Tiedje 2005; Varghese et al. 2015; Chun et al. 2018; Ciufo et al. 2018). The ANI between the whole-genome sequence of HN-1T and that of the most closely related species P. stellifer DSM 14472T was calculated on ezbiocloud (Yoon et al., 2017a). The ANI value of genomes for HN-1T and 39 was 99.93%, but the ANI values of strain HN-1T were less than 91% with the related species P. stellifer DSM 14472T, P. sabinae DSM 17841TT and P. durus 35681T (Table 4), lower than the threshold of 95–96% for differentiating bacterial species (Richter and Rossello´-Mo´ra 2009), suggesting that the new isolate HN-1T represents a distinctive species.
Table 4
Pairwise genome comparisons between the strains HN-1T and 39, and their close relatives
Query Genome
|
Reference genome
|
ANI (%)
|
dDDH (%)
|
HN-1T
|
39
|
99.93
|
99.8
|
HN-1T
|
P. stellifer DSM 14472T (CP009286)
|
90.69
|
42.1
|
HN-1T
|
P. sabinae DSM 17841T (CP004078)
|
76.96
|
22.0
|
HN-1T
|
P. durus ATCC 35681T (CP011114)
|
76.80
|
22.0
|
39
|
P. stellifer DSM 14472T (CP009286)
|
90.83
|
42.2
|
39
|
P. sabinae DSM 17841T (CP004078)
|
77.01
|
22.2
|
39
|
P. durus ATCC 35681T (CP011114)
|
76.78
|
22.2
|
The nitrogen fixation genes of strains HN -1T and 39 were extracted by using Prokka software from the genome sequences (Seemann 2014). The genome of strains HN-1T and 39 contain a compact nif cluster comprising ten genes nifB, nifH, nifD, nifK, nifE, nifN, nifX, orf1, hesA and nifV encoding Mo-nitrogenase, which is unique features of the Paenibacillus nitrogen fixation system. In addition to the nif cluster, the two strains have anfHDGK encoding Fe-nitrogenase and linked to additional copies of nifBENX genes, while the closely related species P. stellifer DSM 14472T contains anfHDGK preceding additional nifV gene. Beyond the nif and anf cluster, there are multiple nifHDK-like genes located at different sites in their genomes. The organization of nif, anf and nif-like genes in type strain HN-1T and the closely related species P. stellifer DSM 14472T was shown in Fig. S3. Previous studies showed that 3 nifH genes of P. sabinae DSM 17841T are functional by complementing K. oxytoca ΔnifH mutant (Hong et al. 2012). Thus, the high nitrogenase activity exhibited by these strains may be due to their additional nif genes.
Paeniacillus durus ATCC 35681T can fix nitrogen even in the presence of nitrate due to the absence of nitrate reductase (Seldin et al. 1984). Whole genome sequence analysis strains HN -1T and 39 revealed that nitrate reductase gene cluster narIJHG were not detected, which suggested these two strains can also fix nitrogen in the nitrate-enriched medium. The draft genome of strains HN -1T and 39 harbor two sets of NAD(P)H-nitrite reductases (nirBD) which are involved in the reduction of nitrite to ammonium in both assimilatory and dissimilatory reduction processes. Additional searches for genes associated with nitric oxide (nirS or nirK) and nitrous oxide reduction (norBC) were performed, but these genes were not detected in their genomes. Therefore, strains HN-1T and 39 may possess dissimilatory nitrate reduction to ammonium pathway, but lack denitrification pathway.
The phenotypic, phylogenetic and genomic data of strains HN-1T and 39 showed that they are different from all other closely related species of genus Paenibacillus. Therefore, we conclude that strain HN-1T or 39 should be recognised as a novel species of the genus Paenibacillus, for which the name Paenibacillus sinensis sp. nov. is proposed.
Description of Paenibacillus sinensis sp. nov.
Paenibacillus sinensis (sin. en’sis. L.gen. n. sinensis of China, where the type strain HN-1T was isolated).
Cells are Gram-positive, facultative anaerobic, rod-shaped (0.4–0.5 µm×2.0–3.2 µm) and motile by means of peritrichous flagella. In slightly swollen sporangia, an ellipsoidal spore is formed and located in terminal position of cells. Colonies on LD medium are circular, convex, cream white, with diameter 1.0–2.0 mm. Nitrogen fixation positive and multiple nifH genes are present. The growth temperature is 15–42°C, optimal at 30°C. The growth pH range is 5.0–9.0, optimal at pH 7.0. NaCl concentration of 0–4% (w/v) is tolerable for growth, optimal at 0–0.2%. Positive tests for catalase, methyl red test, starch and aesculin hydrolysis, but negative for oxidase, Voges–Proskauer reaction, nitrate reduction. The various substrates are assimilated examined using Biolog GEN III microplates: dextrin, D-maltose, D-trehalose, D-cellobiose, D-gentiobiose, sucrose, D-turanose, stachyose, D-raffinose, α-D-lactose, D-melibiose, β-Methyl-D-glucoside, D-salicin, α-D-glucose, D-mannose, D-fructose, D-galactose, 3-methyl glucose, 1% sodium lactate, D-serine, D-sorbitol and pectin were utilized. Strains are resistant to inhibitory chemicals: aztreonam, nalidixic acid, vancomycin, lithium chloride, potassium tellurite, sodium bromate, sodium butyrate, sodium lactate 1% and sensitive to troleandomycin, lincomycin, guanidine HCl, niaproof 4, tetrazolium blue. The major menaquinone is MK-7. The predominant fatty acid is anteiso-C15:0. The major polar lipids are DPG, PE, and PG. The DNA G + C contents for strains HN-1T and 39 are 53.36 and 52.99 mol%, respectively.
The type strain, HN-1T (= CGMCC 1.18902, JCM 34620), was isolated from the rhizosphere soil of rice in Hunan P. R. China. The GenBank (EMBL) accession number for the 16S rRNA gene sequence of strain HN-1T is MF967304 and the GenBank accession number for the draft genome sequence is JAHCMB000000000.