Soybeans inoculated with root zone soils of Canadian native legumes harbour diverse and novel Bradyrhizobium spp. that possess agricultural potential

Diazotrophic community in the sediments of Poyang Lake in response to water level fluctuations

Water level fluctuations (WLFs) are typical characteristic of floodplain lakes and dominant forces regulating the structure and function of lacustrine ecosystems. The sediment diazotrophs play important roles in contributing bioavailable nitrogen to the aquatic environment. However, the relationship between the diazotrophic community and WLFs in floodplain lakes is unknown. In this paper, we carried out a comprehensive investigation on the alpha diversity, abundance, composition and co-occurrence network of the sediment diazotrophs during different water level phases in Poyang Lake. There were no regular variation patterns in the alpha diversity and abundance of the sediment diazotrophs with the water level phase transitions. The relative abundance of some diazotrophic phyla (including Alphaproteobacteria, Deltaproteobacteri, Euryarchaeota, and Firmicutes) and genera (including Geobacter, Deferrisoma, Desulfuromonas, Rivicola, Paraburkholderia, Methylophilus, Methanothrix, Methanobacterium, and Clostridium) was found to change with the water level phase transitions. The results of ANOSIM, PerMANOVA, and DCA at the OTU level showed that the diazotrophic community structure in the low water level phase was significantly different from that in the two high water level phases, while there was no significant difference between the two high water level phases. These results indicated that the diazotrophic community was affected by the declining water level in terms of the composition, while the rising water level contributed to the recoveries of the diazotrophic community. The diazotrophs co-occurrence network was disrupted by the declining water level, but it was strengthened by the rising water level. Moreover, redundancy analysis showed that the variation of the diazotrophic community composition was mostly related to sediment total nitrogen (TN) and total phosphorous (TP). Interestingly, the levels of sediment TN and TP were also found to vary with the water level phase transitions. Therefore, it might be speculated that the WLFs may influence the sediment TN and TP, and in turn influence the diazotrophic community composition. These data can contribute to broadening our understanding of the ecological impacts of WLFs and the nitrogen fixation process in floodplain lakes.

Strains of Bradyrhizobium barranii sp. nov. associated with legumes native to Canada are symbionts of soybeans and belong to different subspecies (subsp. barranii subsp. nov. and subsp. apii subsp. nov.) and symbiovars (sv. glycinearum and sv. septentrionale)

Four bacterial strains isolated from root nodules of soybean plants that had been inoculated with root-zone soils of legumes native to Canada were previously identified as a novel Bradyrhizobium lineage consisting of symbiovars (sv.) glycinearum and septentrionale. Our purpose was to verify the taxonomic status of these strains using phylogenetic, genomic and phenotypic analyses. Multiple phylogenetic analyses including analysis of 51 full-length ribosome protein subunit (rps) gene sequences confirmed placement of the novel strains in a highly supported lineage distinct from named Bradyrhizobium species with B. japonicum USDA 6T as the closest relative. The results of genomic and phylogenomic analyses based on digital DNA–DNA hybridization and genome blast distance phylogeny showed that novel strains in comparisons with type strains of closest relatives were below the established threshold (70 %) for species delineation. Moreover, the novel strains were divided into two subspecies clusters based on the established threshold of 79 %. The genomes of strains 144S4T, 323S2, 1S5 and 38S5T have sizes of 11 399 526, 11 474 152, 10580853 and 10 530 141 bp with DNA G+C contents of 63.1, 63.0, 63.4 and 63.3 mol%, respectively. These strains possess symbiosis islands harbouring key nodulation, nitrogen-fixation and type III secretion system genes as well as abundant insertion sequences and between two and four putative plasmids. Strains 144S4T and 323S2 (sv. glycinearum) are effective with regard to nitrogen fixation in symbiotic association with soybeans whereas strains 1S5 and 38S5T (sv. septentrionale) are ineffective. Data for morphological, physiological and symbiotic characteristics complement the sequence-based results. The data presented here support the description of a new species and two new subspecies for which the names Bradyrhizobium barranii sp. nov. subsp. barranii subsp. nov. (sv. glycinearum) and Bradyrhizobium barranii sp. nov. subsp. apii subsp. nov. (sv. septentrionale) are proposed with strain 144S4T (=LMG 31552T=HAMBI 3722T) as the species type strain and type strain of subsp. barranii subsp. nov., and strain 38S5T (=LMG 31556T=HAMBI 3721T) as the type strain of subsp. apii subsp. nov.

Genetic and phenotypic diversity of microsymbionts nodulating promiscuous soybeans from different agro-climatic conditions

Background

Global food supply is highly dependent on field crop production that is currently severely threatened by changing climate, poor soil quality, abiotic, and biotic stresses. For instance, one of the major challenges to sustainable crop production in most developing countries is limited nitrogen in the soil. Symbiotic nitrogen fixation of legumes such as soybean (Glycine max (L.) Merril) with rhizobia plays a crucial role in supplying nitrogen sufficient to maintain good crop productivity. Characterization of indigenous bradyrhizobia is a prerequisite in the selection and development of effective bioinoculants. In view of this, bradyrhizobia were isolated from soybean nodules in four agro-climatic zones of eastern Kenya (Embu Upper Midland Zone, Embu Lower Midland Zone, Tharaka Upper Midland Zone, and Tharaka Lower Midland Zone) using two soybean varieties (SB8 and SB126). The isolates were characterized using biochemical, morphological, and genotypic approaches. DNA fingerprinting was carried out using 16S rRNA gene and restricted by enzymes HaeIII, Msp1, and EcoRI. 

Results

Thirty-eight (38) bradyrhizobia isolates obtained from the trapping experiments were placed into nine groups based on their morphological and biochemical characteristics. Most (77%) of the isolates had characteristics of fast-grower bradyrhizobia while 23% were slow-growers. Restriction digest revealed significant (p < 0.015) variation within populations and not among the agro-climatic zones based on analysis of molecular variance. Principal coordinate analysis demonstrated sympatric speciation of indigenous bradyrhizobia isolates. Embu Upper Midland Zone bradyrhizobia isolates had the highest polymorphic loci (80%) and highest genetic diversity estimates (H’ = 0.419) compared to other agro-climatic zones.

Conclusion

The high diversity of bradyrhizobia isolates depicts a valuable genetic resource for selecting more effective and competitive strains to improve promiscuous soybean production at a low cost through biological nitrogen fixation.

Bradyrhizobium septentrionale sp. nov. (sv. septentrionale) and Bradyrhizobium quebecense sp. nov. (sv. septentrionale) associated with legumes native to Canada possess rearranged symbiosis genes and numerous insertion sequences

Six bacterial strains isolated from root nodules of soybean plants that had been inoculated with root-zone soil of legumes native to Canada were previously characterized and 1) placed in two novel lineages within the genus Bradyrhizobium and 2) assigned to symbiovar septentrionale. Here we verified the taxonomic status of these strains using genomic and phenotypic analyses. Phylogenetic analyses of five protein encoding partial gene sequences as well as 52 full length ribosome protein subunit gene sequences confirmed placement of the novel strains in two highly supported lineages distinct from named Bradyrhizobium species. The highest average nucleotide identity values of strains representing these two lineages relative to type strains of closest relatives were 90.7 and 92.3% which is well below the threshold value for bacterial species circumscription. The genomes of representative strains 1S1^T, 162S2 and 66S1MB^T have sizes of 10598256, 10733150 and 9032145 bp with DNA G+C contents of 63.5, 63.4 and 63.8 mol%, respectively. These strains possess between one and three plasmids based on copy number of plasmid replication and segregation (repABC) genes. Novel strains also possess numerous insertion sequences, and, relative to reference strain Bradyrhizobium diazoefficiens USDA110^T, exhibit inversion and fragmentation of nodulation (nod) and nitrogen-fixation (nif) gene clusters. Phylogenetic analyses of nodC and nifH gene sequences confirmed placement of novel strains in a distinct lineage corresponding to symbiovar septentrionale. Data for morphological, physiological and symbiotic characteristics complement the sequence-based results. The data presented here support the description of two new species for which the names Bradyrhizobium septentrionale sp. nov. (sv. septentrionale) and Bradyrhizobium quebecense sp. nov. (sv. septentrionale) are proposed, with 1S1^T (=LMG 29930^T=HAMBI 3676^T) and 66S1MB^T (=LMG 31547^T=HAMBI 3720^T) as type strains, respectively.

Soybean Root Nodule and Rhizosphere Microbiome: Distribution of Rhizobial and Nonrhizobial Endophytes

Soybean root nodules are known to contain a high diversity of both rhizobial endophytes and nonrhizobial endophytes (NREs). Nevertheless, the variation of these bacteria among different root nodules within single plants has not been reported. So far, it is unclear whether the selection of NREs among different root nodules within single plants is a random process or is strictly controlled by the host plant to favor a few specific NREs based on their beneficial influence on plant growth. As well, it is also unknown if the relative frequency of NREs within different root nodules is consistent or if it varies based on the location or size of a root nodule. We assessed the microbiomes of 193 individual soybean root nodules from nine plants using high-throughput DNA sequencing. Bradyrhizobium japonicum strains occurred in high abundance in all root nodules despite the presence of other soybean-compatible rhizobia, such as Ensifer, Mesorhizobium, and other species of Bradyrhizobium in soil. Nitrobacter and Tardiphaga were the two nonrhizobial genera that were uniformly detected within almost all root nodules, though they were in low abundance. DNA sequences related to other NREs that have frequently been reported, such as Bacillus, Pseudomonas, Flavobacterium, and Variovorax species, were detected in a few nodules. Unlike for Bradyrhizobium, the low abundance and inconsistent occurrence of previously reported NREs among different root nodules within single plants suggest that these microbes are not preferentially selected as endophytes by host plants and most likely play a limited part in plant growth as endophytes.IMPORTANCE Soybean (Glycine max L.) is a valuable food crop that also contributes significantly to soil nitrogen by developing a symbiotic association with nitrogen-fixing rhizobia. Bacterial endophytes (both rhizobial and nonrhizobial) are considered critical for the growth and resilience of the legume host. In the past, several studies have suggested that the selection of bacterial endophytes within root nodules can be influenced by factors such as soil pH, nutrient availability, host plant genotype, and bacterial diversity in soil. However, the influence of size or location of root nodules on the selection of bacterial endophytes within soybean roots is unknown. It is also unclear whether the selection of nonrhizobial endophytes within different root nodules of a single plant is a random process or is strictly regulated by the host. This information can be useful in identifying potential bacterial species for developing bioinoculants that can enhance plant growth and soil nitrogen.

Bacterial community dynamics across developmental stages of fungal fruiting bodies

Increasing evidence suggest that bacteria form diverse communities in various eukaryotic hosts, including fungi. However, little is known about their succession and the functional potential at different host development stages. Here we examined the effect of fruiting body parts and developmental stages on the structure and potential function of fungus-associated bacterial communities. Using high-throughput sequencing, we characterized bacterial communities and their associated potential functions in fruiting bodies from ten genera belonging to four major mushroom-forming orders and three different developmental stages of a model host species Cantharellus cibarius. Our results demonstrate that bacterial community structure differs between internal and external parts of the fruiting body but not between inner tissues. The structure of the bacterial communities showed significant variation across fruiting body developmental stages. We provide evidence that certain functional groups, such as those related to nitrogen fixation, persist in fruiting bodies during the maturation, but are replaced by putative parasites/pathogens afterwards. These data suggest that bacterial communities inhabiting fungal fruiting bodies may play important roles in their growth and development.

Bradyrhizobium hipponense sp. nov., isolated from Lupinus angustifolius growing in the northern region of Tunisia

Strain aSej3^T was isolated from a root nodule of a Lupinus angustifolius plant growing in Bizerte, Tunisia. 16S rRNA gene analysis placed this strain within the genus Bradyrhizobium. Multilocus sequence analysis (MLSA) including three housekeeping genes (glnII, gyrB and recA) grouped aSej3^T together with Bradyrhizobium rifense CTAW71^T, Bradyrhizobium cytisi CTAW11^T, Bradyrhizobium ganzhouense RITF806^T, Bradyrhizobium lupini USDA 3051^T and Bradyrhizobium canariense BTA-1^T. MLSA with five housekeeping genes (dnaK, glnII, gyrB, recA and rpoB) revealed that this strain shares less than 93.5 % nucleotide identity with other type strains. Genome sequencing and inspection revealed a genome size of 8.83 Mbp with a G+C content of 62.8 mol%. Genome-wide average nucleotide identity and digital DNA-DNA hybridization values were below 87.5 and 36.2 %, respectively, when compared to described Bradyrhizobium species. Strain aSej3^T nodulated L. angustifolius plants under axenic conditions and its nodC gene clustered within the genistearum symbiovar. Altogether, the phylogenetic data and the chemotaxonomic characteristics of this strain support that aSej3^T represents a new species for which we propose the name Bradyrhizobium hipponense sp. nov. with the type strain aSej3^T (=DSM 108913^T=LMG 31020^T).

Strains of Bradyrhizobium cosmicum sp. nov., isolated from contrasting habitats in Japan and Canada possess photosynthesis gene clusters with the hallmark of genomic islands

The taxonomic status of two previously characterized Bradyrhizobium strains (58S1T and S23321) isolated from contrasting habitats in Canada and Japan was verified by genomic and phenotypic analyses. Phylogenetic analyses of five and 27 concatenated protein-encoding core gene sequences placed both strains in a highly supported lineage distinct from named species in the genus Bradyrhizobium with Bradyrhizobium betae as the closest relative. Average nucleotide identity values of genome sequences between the test and reference strains were between 84.5 and 94.2 %, which is below the threshold value for bacterial species circumscription. The complete genomes of strains 58S1T and S23321 consist of single chromosomes of 7.30 and 7.23 Mbp, respectively, and do not have symbiosis islands. The genomes of both strains have a G+C content of 64.3 mol%. Present in the genome of these strains is a photosynthesis gene cluster (PGC) containing key photosynthesis genes. A tRNA gene and its partial tandem duplication were found at the boundaries of the PGC region in both strains, which is likely the hallmark of genomic island insertion. Key nitrogen-fixation genes were detected in the genomes of both strains, but nodulation and type III secretion system genes were not found. Sequence analysis of the nitrogen fixation gene, nifH, placed 58S1T and S23321 in a novel lineage distinct from described Bradyrhizobium species. Data for phenotypic tests, including growth characteristics and carbon source utilization, supported the sequence-based analyses. Based on the data presented here, a novel species with the name Bradyrhizobium cosmicum sp. nov. is proposed with 58S1T (=LMG 31545T=HAMBI 3725T) as the type strain.

Intrinsic role of bacterial secretion systems in phylogenetic niche conservation of Bradyrhizobium spp

Bradyrhizobium is a biologically important bacterial genus. Different Bradyrhizobium strains exhibit distinct niche selection like free living, root nodular and stem nodular. The present in-silico study was undertaken to identify the role of bacterial secretome in the phylogenetic niche conservation (PNC) of Bradyrhizobium sp. Analysis was carried out with the publicly available 19 complete genome assembly and annotation reports. A protocol was developed to screen the secretome related genes using three different database, viz. genome, proteome and gene ortholog. This resulted into 139 orthologs that include type secretion systems (T1SS-T6SS) along with flagella (Flg), type IV pili (T4P) and tight adherence (Tad) systems. Multivariate analysis using bacterial secretome was undertaken to find out the role of these secretion systems in PNC. In free living strains, T3SS, T4SS and T6SS were completely absent. Whereas, in the stem nodulating strains, T3SS and T6SS were absent, but T4SS was found to be present. On the other hand, the T3SS was found to be present only in the root-nodulating strains. The present investigation clearly demonstrated a pattern of PNC based on the distribution of secretion system components. To the best of our knowledge, this is the first report on PNC of Bradyrhizobium using the multivariate analysis of secretome.

Description and complete genome sequences of Bradyrhizobium symbiodeficiens sp. nov., a non-symbiotic bacterium associated with legumes native to Canada

Four bacterial strains isolated from root nodules of soybean plants that had been inoculated with root-zone soil of either Amphicarpaea bracteata (Hog Peanut) or Desmodium canadense (Showy Tick Trefoil) growing in Canada, were previously characterized and placed in a novel lineage within the genus Bradyrhizobium. The taxonomic status of the novel strains was verified by genomic and phenotypic analyses. Phylogenetic analyses of individual and concatenated housekeeping gene sequences (atp D, gln II, rec A, gyr B and rpo B) placed all novel strains in a highly supported lineage distinct from named Bradyrhizobium species. Data for sequence similarities of concatenated housekeeping genes of novel strains relative to type strains of named species were consistent with the phylogenetic data. Average nucleotide identity values of genome sequences (84.5-93.7 %) were below the threshold value of 95-96 % for bacterial species circumscription. Close relatives to the novel strains are Bradyrhizobium amphicarpaeae, Bradyrhizobium ottawaense and Bradyrhizobium shewense. The complete genomes of strains 85S1MBT and 65S1MB consist of single chromosomes of size 7.04 and 7.13 Mbp, respectively. The genomes of both strains have a G+C content of 64.3 mol%. These strains lack a symbiosis island as well as key nodulation, nitrogen-fixation and photosystem genes. Data from various phenotypic tests including growth characteristics and carbon source utilization supported the sequence-based analyses. Based on the data presented here, the four strains represent a novel species for which the name B radyrhizobium symbiodeficiens sp. nov., is proposed, with 85S1MBT (=LMG 29937T=HAMBI 3684T) as the type strain.

Description and complete genome sequence of Bradyrhizobium amphicarpaeae sp. nov., harbouring photosystem and nitrogen-fixation genes

A bacterial strain, designated 39S1MBT, isolated from a root nodule of a soybean plant that had been inoculated with root-zone soil of Amphicarpaea bracteata (hog peanut) growing in Canada, was previously characterized and placed in a novel lineage within the genus Bradyrhizobium. The taxonomic status of strain 39S1MBT was verified by genomic and phenotypic analyses. Phylogenetic analyses of individual and concatenated protein-encoding gene sequences (atpD, glnII, recA, gyrB and rpoB) placed 39S1MBT in a lineage distinct from named species. Data for sequence similarities of concatenated genes relative to type strains of named species supported the phylogenetic data. Average nucleotide identity values of genome sequences (84.5-91.7 %) were well below the threshold value for bacterial species circumscription. Based on these data, Bradyrhizobium ottawaense OO99T and Bradyrhizobium shewense ERR11T are close relatives of 39S1MBT. The complete genome of 39S1MBT consists of a single 7.04 Mbp chromosome without a symbiosis island; G+C content is 64.7 mol%. Present in the genome are key photosystem and nitrogen-fixation genes, but not nodulation and type III secretion system genes. Sequence analysis of the nitrogen fixation gene, nifH, placed 39S1MBT in a novel lineage distinct from named Bradyrhizobium species. Data for phenotypic tests including growth characteristics and carbon source utilization supported the sequence-based analyses. Based on the data presented here, a novel species with the name Bradyrhizobium amphicarpaeae sp. nov. is proposed with 39S1MBT (=LMG 29934T=HAMBI 3680T) as the type strain.

Genome-informed Bradyrhizobium taxonomy: where to from here?

Bradyrhizobium is thought to be the largest and most diverse rhizobial genus, but this is not reflected in the number of described species. Although it was one of the first rhizobial genera recognised, its taxonomy remains complex. Various contemporary studies are showing that genome sequence information may simplify taxonomic decisions. Therefore, the growing availability of genomes for Bradyrhizobium will likely aid in the delineation and characterization of new species. In this study, we addressed two aims: first, we reviewed the availability and quality of available genomic resources for Bradyrhizobium. This was achieved by comparing genome sequences in terms of sequencing technologies used and estimated level of completeness for inclusion in genome-based phylogenetic analyses. Secondly, we utilized these genomes to investigate the taxonomic standing of Bradyrhizobium in light of its diverse lifestyles. Although genome sequences differed in terms of their quality and completeness, our data indicate that the use of these genome sequences is adequate for taxonomic purposes. By using these resources, we inferred a fully resolved, well-supported phylogeny. It separated Bradyrhizobium into seven lineages, three of which corresponded to the so-called supergroups known for the genus. Wide distribution of key lifestyle traits such as nodulation, nitrogen fixation and photosynthesis revealed that these traits have complicated evolutionary histories. We present the first robust Bradyrhizobium species phylogeny based on genome sequence information for investigating the evolution of this important assemblage of bacteria. Furthermore, this study provides the basis for using genome sequence information as a resource to make important taxonomic decisions, particularly at the species and genus levels.

Diversity of bacteria associated with corn roots inoculated with Canadian woodland soils, and description of Pseudomonas aylmerense sp. nov

Bacteria associated with corn roots inoculated with soils collected from the Canadian woodlands were isolated and characterized. Genus-level identification based on 16S rRNA sequence analysis classified the 161 isolates in 19 genera. The majority (64%) of the isolates were affiliated with the genus Pseudomonas. Further analysis of the Pseudomonas isolates based on BLASTn and rpoD-rpoB-gyrB concatenated gene phylogeny revealed three unique clusters that could not be assigned to known species. This study reports the taxonomic description of one of the distinct lineages represented by two strains (S1E40^T and S1E44) with P. lurida LMG 21995^T, P. costantinii LMG 22119^T, P. palleroniana LMG 23076^T, P. simiae CCUG 50988^T and P. extremorientalis LMG 19695^T as the closest taxa. Both strains showed low ANIm (<90%) and genome-based DNA-DNA hybridization (<50%) values, which unequivocally delineated the new strains from the closest relatives. These findings were supported by multilocus sequence analysis (MLSA) and DNA fingerprinting. In addition, growth characteristics and biochemical tests revealed patterns that differed from the related species. Strains S1E40^T and S1E44 are Gram-negative, aerobic, rod-shaped and motile by at least one flagellum; and grew optimally at 30 °C. The predominant polar lipid is phosphatidylethanolamine while the major respiratory quinone is ubiquinone-9. Based on phenotypic and genotypic data presented here, strains S1E40^T and S1E44 represent a novel species for which the name Pseudomonas aylmerense sp. nov. is proposed. The type strain is S1E40^T (= LMG 30784^T = DOAB 703^T = HAMI 3696^T) with a G + C content of 61.6%.