Bradyrhizobium centrosemae (symbiovar centrosemae) sp. nov., Bradyrhizobium americanum (symbiovar phaseolarum) sp. nov. and a new symbiovar (tropici) of Bradyrhizobium viridifuturi establish symbiosis with Centrosema species native to America

Guidelines for the description of rhizobial symbiovars

Rhizobia are bacteria that form nitrogen-fixing nodules in legume plants. The sets of genes responsible for both nodulation and nitrogen fixation are carried in plasmids or genomic islands that are often mobile. Different strains within a species sometimes have different host specificities, while very similar symbiosis genes may be found in strains of different species. These specificity variants are known as symbiovars, and many of them have been given names, but there are no established guidelines for defining or naming them. Here, we discuss the requirements for guidelines to describe symbiovars, propose a set of guidelines, provide a list of all symbiovars for which descriptions have been published so far, and offer a mechanism to maintain a list in the future.

Novel symbiovars ingae, lysilomae and lysilomaefficiens in bradyrhizobia from tree-legume nodules

Inga vera and Lysiloma tree legumes form nodules with Bradyrhizobium spp. from the japonicum group that represent novel genomospecies, for which we describe here using genome data, symbiovars lysilomae, lysilomaefficiens and ingae. Genes encoding Type three secretion system (TTSS) that could affect host specificity were found in ingae but not in lysilomae nor in lysilomaefficiens symbiovars and uptake hydrogenase hup genes (that affect nitrogen fixation) were observed in bradyrhizobia from the symbiovars ingae and lysilomaefficiens. nolA gene was found in the symbiovar lysilomaefficiens but not in strains from lysilomae. We discuss that multiple genes may dictate symbiosis specificity. Besides, toxin-antitoxin genes were found in the symbiosis islands in bradyrhizobia from symbiovars ingae and lysilomaefficiens. A limit (95%) to define symbiovars with nifH gene sequences was proposed here.

Lupin, a Unique Legume That Is Nodulated by Multiple Microsymbionts: The Role of Horizontal Gene Transfer

Lupin is a high-protein legume crop that grows in a wide range of edaphoclimatic conditions where other crops are not viable. Its unique seed nutrient profile can promote health benefits, and it has been proposed as a phytoremediation plant. Most rhizobia nodulating Lupinus species belong to the genus Bradyrhizobium, comprising strains that are phylogenetically related to B. cytisi, B. hipponenese, B. rifense, B. iriomotense/B. stylosanthis, B. diazoefficiens, B. japonicum, B. canariense/B. lupini, and B. retamae/B. valentinum. Lupins are also nodulated by fast-growing bacteria within the genera Microvirga, Ochrobactrum, Devosia, Phyllobacterium, Agrobacterium, Rhizobium, and Neorhizobium. Phylogenetic analyses of the nod and nif genes, involved in microbial colonization and symbiotic nitrogen fixation, respectively, suggest that fast-growing lupin-nodulating bacteria have acquired their symbiotic genes from rhizobial genera other than Bradyrhizobium. Horizontal transfer represents a key mechanism allowing lupin to form symbioses with bacteria that were previously considered as non-symbiotic or unable to nodulate lupin, which might favor lupin’s adaptation to specific habitats. The characterization of yet-unstudied Lupinus species, including microsymbiont whole genome analyses, will most likely expand and modify the current lupin microsymbiont taxonomy, and provide additional knowledge that might help to further increase lupin’s adaptability to marginal soils and climates.

Bradyrhizobium cenepequi sp. nov., Bradyrhizobium semiaridum sp. nov., Bradyrhizobium hereditatis sp. nov. and Bradyrhizobium australafricanum sp. nov., symbionts of different leguminous plants of Western Australia and South Africa and definition of three novel symbiovars

Bradyrhizobium is a heterogeneous bacterial genus capable of establishing symbiotic associations with a broad range of legume hosts, including species of economic and environmental importance. This study was focused on the taxonomic and symbiovar definition of four strains – CNPSo 4026T, WSM 1704T, WSM 1738T and WSM 4400T – previously isolated from nodules of legumes in Western Australia and South Africa. The 16S rRNA gene phylogenetic tree allocated the strains to the Bradyrhizobium elkanii supergroup. The multilocus sequence analysis (MLSA) with partial sequences of six housekeeping genes – atpD, dnaK, glnII, gyrB, recA and rpoB – did not cluster the strains under study as conspecific to any described Bradyrhizobium species. Average nucleotide identity and digital DNA–DNA hybridization values were calculated for the four strains of this study and the closest species according to the MLSA phylogeny with the highest values being 95.46 and 62.20 %, respectively; therefore, both being lower than the species delineation cut-off values. The nodC and nifH phylogenies included strains WSM 1738T and WSM 4400T in the symbiovars retamae and vignae respectively, and also allowed the definition of three new symbiovars, sv. cenepequi, sv. glycinis, and sv. cajani. Analysis of morphophysiological characterization reinforced the identification of four novel proposed Bradyrhizobium species that are accordingly named as follows: Bradyrhizobium cenepequi sp. nov. (CNPSo 4026T=WSM 4798T=LMG 31653T), isolated from Vigna unguiculata; Bradyrhizobium semiaridum sp. nov. (WSM 1704T=CNPSo 4028T=LMG 31654T), isolated from Tephrosia gardneri; Bradyrhizobium hereditatis sp. nov. (WSM 1738T=CNPSo 4025T=LMG 31652T), isolated from Indigofera sp.; and Bradyrhizobium australafricanum sp. nov. (WSM 4400T=CNPSo 4015T=LMG 31648T) isolated from Glycine sp.

Genotypic and symbiotic diversity studies of rhizobia nodulating Acacia saligna in Tunisia reveal two novel symbiovars within the Rhizobium leguminosarum complex and Bradyrhizobium

Acacia saligna is an invasive alien species that has the ability to establish symbiotic relationships with rhizobia. In the present study, genotypic and symbiotic diversity of native rhizobia associated with A. saligna in Tunisia were studied. A total of 100 bacterial strains were selected and three different ribotypes were identified based on rrs PCR-RFLP analysis. Sequence analyses of rrs and four housekeeping genes (recA, atpD, gyrB and glnII) assigned 30 isolates to four putative new lineages and a single strain to Sinorhizobium meliloti. Thirteen slow-growing isolates representing the most dominant IGS (intergenic spacer) profile clustered distinctly from known rhizobia species within Bradyrhizobium with the closest related species being Bradyrhizobium shewense and Bradyrhizobium niftali, which had 95.17% and 95.1% sequence identity, respectively. Two slow-growing isolates, 1AS28L and 5AS6L, had B. frederekii as their closest species with a sequence identity of 95.2%, an indication that these strains could constitute a new lineage. Strains 1AS14I, 1AS12I and 6AS6 clustered distinctly from known rhizobia species but within the Rhizobium leguminosarum complex (Rlc) with the most closely related species being Rhizobium indicum with 96.3% sequence identity. Similarly, the remaining 11 strains showed 96.9 % and 97.2% similarity values with R. changzhiense and R. indicum, respectively. Based on nodC and nodA phylogenies and cross inoculation tests, these 14 strains of Rlc species clearly diverged from strains of Sinorhizobium and Rlc symbiovars, and formed a new symbiovar for which the name sv. "salignae" is proposed. Bacterial strains isolated in this study that were taxonomically assigned to Bradyrhizobium harbored different symbiotic genes and the data suggested a new symbiovar, for which sv. "cyanophyllae" is proposed. Isolates formed effective nodules on A. saligna.

Bradyrhizobium archetypum sp. nov., Bradyrhizobium australiense sp. nov. and Bradyrhizobium murdochi sp. nov., isolated from nodules of legumes indigenous to Western Australia

The genus Bradyrhizobium is considered as the probable ancestor lineage of all rhizobia, broadly spread in a variety of ecosystems and with remarkable diversity. A polyphasic study was performed to characterize and clarify the taxonomic position of eight bradyrhizobial strains isolated from indigenous legumes to Western Australia. As expected for the genus, the 16S rRNA gene sequences were highly conserved, but the results of multilocus sequence analysis with four housekeeping genes (dnaK, glnII, gyrB and recA) confirmed three new distinct clades including the following strains: (1) WSM 1744^T, WSM 1736 and WSM 1737; (2) WSM 1791^T and WSM 1742; and (3) WSM 1741^T, WSM 1735 and WSM 1790. The highest ANI values of the three groups in relation to the closest type strains were 92.4, 92.3 and 93.3 %, respectively, below the threshold of species circumscription. The digital DNA-DNA hybridization analysis also confirmed new species descriptions, with less than 52 % relatedness with the closest type strains. The phylogeny of the symbiotic gene nodC clustered the eight strains into the symbiovar retamae, together with seven Bradyrhizobium type strains, sharing from 94.2-98.1 % nucleotide identity (NI), and less than 88.7 % NI with other related strains and symbiovars. Morpho-physiological, phylogenetics, genomic and symbiotic traits were determined for the new groups and our data support the description of three new species, Bradyrhizobium archetypum sp. nov., Bradyrhizobium australiense sp. nov. and Bradyrhizobium murdochi sp. nov., with WSM 1744^T (=CNPSo 4013^T=LMG 31646^T), WSM 1791^T (=CNPSo 4014^T=LMG 31647^T) and WSM 1741^T (=CNPSo 4020^T=LMG 31651^T) designated as type strains, respectively.

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).

Efficient Nitrogen-Fixing Bacteria Isolated from Soybean Nodules in the Semi-arid Region of Northeast Brazil are Classified as Bradyrhizobium brasilense (Symbiovar Sojae)

Soybean (Glycine max L.) is an important legume that greatly benefits from inoculation with nitrogen-fixing bacteria. In a previous study, five efficient nitrogen-fixing bacterial strains, isolated from nodules of soybean inoculated with soil from semi-arid region, Northeast Brazil, were identified as a new group within the genus Bradyrhizobium. The taxonomic status of these strains was evaluated in this study. The phylogenetic analysis of the 16S rRNA gene showed the high similarity of the five strains to Bradyrhizobium brasilense UFLA03-321T (100%), B. pachyrhizi PAC48T (100%), B. ripae WR4T (100%), B. elkanii USDA 76T (99.91%), and B. macuxiense BR 10303T (99.91%). However, multilocus sequence analysis of the housekeeping genes atpD, dnaK, gyrB, recA, and rpoB, average nucleotide identity, and digital DNA-DNA hybridization analyses supported the classification of the group as B. brasilense. Some phenotypic characteristics allowed differentiating the five strains and the type strain of B. brasilense from the two neighboring species (B. pachyrhizi PAC48T and B. elkanii USDA 76T). The nodC and nifH genes' analyses showed that these strains belong to symbiovar sojae, together with B. elkanii (USDA 76T) and B. ferriligni (CCBAU 51502T). The present results support the classification of these five strains as Bradyrhizobium brasilense (symbiovar sojae).

The Mimosoid tree Leucaena leucocephala can be nodulated by the symbiovar genistearum of Bradyrhizobium canariense

Leucaena leucocephala is a Mimosoid legume tree indigenous to America that has spread to other continents, although it is not still present in some European countries such as Portugal. Nevertheless, we found that this legume can be nodulated in this country by slow-growing rhizobial strains which were identified as Bradyrhizobium canariense trough the analysis of the core genes recA and glnII. The analysis of the symbiotic gene nodC showed that these strains belong to the symbiovar genistearum, which commonly nodulates Genistoid legumes. Although two strains nodulating L. leucocephala in China and Brazil were classified within the genus Bradyrhizobium, they belong to undescribed species and to the symbiovars glycinearum and tropici, respectively. Therefore, we report here for the first time the ability of L. leucocephala to establish symbiosis with strains of B. canariense sv genistearum confirming the high promiscuity of L. leucocephala, that allows it to establish symbiosis with rhizobia native to different continents increasing its invasiveness potential.

Bradyrhizobium nanningense sp. nov., Bradyrhizobium guangzhouense sp. nov. and Bradyrhizobium zhanjiangense sp. nov., isolated from effective nodules of peanut in Southeast China

Nine slow-growing rhizobia isolated from effective nodules on peanut (Arachis hypogaea) were characterized to clarify the taxonomic status using a polyphasic approach. They were assigned to the genus Bradyrhizobium on the basis of 16S rRNA sequences. MLSA of concatenated glnII-recA-dnaK genes classified them into three species represented by CCBAU 53390^T, CCBAU 51670^T and CCBAU 51778^T, which presented the closest similarity to B. guangxiense CCBAU 53363^T, B. guangdongense CCBAU 51649^T and B. manausense BR 3351^T, B. vignae 7-2^T and B. forestalis INPA 54B^T, respectively. The dDDH (digital DNA-DNA hybridization) and ANI (Average Nucleotide Identity) between the genomes of the three representative strains and type strains for the closest Bradyrhizobium species were less than 42.1% and 91.98%, respectively, below the threshold of species circumscription. Effective nodules could be induced on peanut and Lablab purpureus by all representative strains, while Vigna radiata formed effective nodules only with CCBAU 53390^T and CCBAU 51778^T. Phenotypic characteristics including sole carbon sources and growth features supported the phylogenetic results. Based on the genotypic and phenotypic features, strains CCBAU 53390^T, CCBAU 51670^T and CCBAU 51778^T are designated the type strains of three novel species, for which the names Bradyrhizobium nanningense sp. nov., Bradyrhizobium guangzhouense sp. nov. and Bradyrhizobium zhanjiangense sp. nov. are proposed, respectively.

A Genomotaxonomy View of the Bradyrhizobium Genus

Whole genome analysis of the Bradyrhizobium genus using average nucleotide identity (ANI) and phylogenomics showed the genus to be essentially monophyletic with seven robust groups within this taxon that includes nitrogen-fixing nodule forming bacteria as well as free living strains. Despite the wide genetic diversity of these bacteria no indication was found to suggest that the Bradyrhizobium genus have to split in different taxa. Bradyrhizobia have larger genomes than other genera of the Bradyrhizobiaceae family, probably reflecting their metabolic diversity and different lifestyles. Few plasmids in the sequenced strains were revealed from rep gene analysis and a relatively low proportion of the genome is devoted to mobile genetic elements. Sequence diversity of recA and glnII gene metadata was used to theoretically estimate the number of existing species and to predict how many would exist. There may be many more species than those presently described with predictions of around 800 species in nature. Different arguments are presented suggesting that nodulation might have arose in the ancestral genus Bradyrhizobium.

Widespread Distribution of Highly Adapted Bradyrhizobium Species Nodulating Diverse Legumes in Africa

Bradyrhizobium is one of the most cosmopolitan and diverse bacterial group nodulating a variety of host legumes in Africa, however, the diversity and distribution of bradyrhizobial symbionts nodulating indigenous African legumes are not well understood, though needed for increased food legume production. In this review, we have shown that many African food legumes are nodulated by bradyrhizobia, with greater diversity in Southern Africa compared to other parts of Africa. From a few studies done in Africa, the known bradyrhizobia (i.e., Bradyrhizobium elkanii, B. yuanmingense) along with many novel Bradyrhizobium species are the most dominant in African soils. This could be attributed to the unique edapho-climatic conditions of the contrasting environments in the continent. More studies are needed to identify the many novel bradyrhizobia resident in African soils in order to better understand the biogeography of bradyrhizobia and their potential for inoculant production.

Classification of the inoculant strain of cowpea UFLA03-84 and of other strains from soils of the Amazon region as Bradyrhizobium viridifuturi (symbiovar tropici)

Cowpea (Vigna unguiculata L.) is a legume species that considerably benefits from inoculation with nitrogen fixing bacteria of the genus Bradyrhizobium. One of the strains recommended for inoculation in cowpea in Brazil is UFLA03-84 (Bradyrhizobium sp.). The aim of our study was to define the taxonomic position of the UFLA03-84 strain and of two other strains of Bradyrhizobium (UFLA03-144 and INPA237B), all belonging to the same phylogenetic group and isolated from soils of the Brazilian Amazon. Multilocus sequence analysis (MLSA) of the housekeeping genes atpD, gyrB, recA, and rpoB grouped (with similarity higher than 99%) the three strains with Bradyrhizobium viridifuturi SEMIA 690^T. The analyses of average nucleotide identity and digital DNA-DNA hybridization supported classification of the group as Bradyrhizobium viridifuturi. The three strains exhibited similar behavior in relation to the most of the phenotypic characteristics evaluated. However, some characteristics exhibited variation, indicating phenotypic diversity within the species. Phylogenetic analysis of the nodC and nifH genes showed that the three strains are members of the same symbiovar (tropici) that contains type strains of Bradyrhizobium species coming from tropical soils (SEMIA 690^TB. viridifuturi, CNPSo 1112^TB. tropiciagri, CNPSo 2833^TB. embrapense, and B. brasilense UFLA03-321^T).

Definition of two new symbiovars, sv. lupini and sv. mediterranense, within the genera Bradyrhizobium and Phyllobacterium efficiently nodulating Lupinus micranthus in Tunisia

In this study, a polyphasic approach was used to analyze three representative strains (LmiH4, LmiM2 and LmiT21) from a collection of six previously described strains isolated in Tunisia from root nodules of Lupinus micranthus. The phylogenetic analysis of the concatenated rrs, recA and glnII genes showed that strain LmiH4 had 100% concatenated gene sequence identity with the type strain Bradyrhizobium retamae Ro19T. Similarly, strain LmiM2 shared 100% concatenated gene sequence identity with the species Bradyrhizobium valentinum LmjM3T. However, strain LmiT21 showed an identical concatenated gene sequence with reference strain Phyllobacterium sophorae CCBAU03422T. The recA-glnII concatenated protein-coding genes used produced incongruent phylogenies compared with 16S rDNA phylogeny. The nodC gene analysis showed that the strains were phylogenetically divergent to the Bradyrhizobium symbiovars defined to date, and represented two new symbiovars. Plant infection analysis revealed that the three strains showed moderate host range and symbiotic specificities. Based on their symbiotic characteristics, we propose that the three strains isolated from Lupinus micranthus nodules belong to two new symbiovars, with the first denominated lupini within the two species Bradyrhizobium valentinum (type strain LmiM2) and B. retamae (type strain LmiH4), and the second denominated mediterranense within the species P. sophorae (type strain LmiT21).

Phylogeny and Phylogeography of Rhizobial Symbionts Nodulating Legumes of the Tribe Genisteae

The legume tribe Genisteae comprises 618, predominantly temperate species, showing an amphi-Atlantic distribution that was caused by several long-distance dispersal events. Seven out of the 16 authenticated rhizobial genera can nodulate particular Genisteae species. Bradyrhizobium predominates among rhizobia nodulating Genisteae legumes. Bradyrhizobium strains that infect Genisteae species belong to both the Bradyrhizobium japonicum and Bradyrhizobium elkanii superclades. In symbiotic gene phylogenies, Genisteae bradyrhizobia are scattered among several distinct clades, comprising strains that originate from phylogenetically distant legumes. This indicates that the capacity for nodulation of Genisteae spp. has evolved independently in various symbiotic gene clades, and that it has not been a long-multi-step process. The exception is Bradyrhizobium Clade II, which unlike other clades comprises strains that are specialized in nodulation of Genisteae, but also Loteae spp. Presumably, Clade II represents an example of long-lasting co-evolution of bradyrhizobial symbionts with their legume hosts.

Bradyrhizobium cajani sp. nov. isolated from nodules of Cajanus cajan

Two slow-growing strains, AMBPC1010T and AMBPC1011, were isolated from nodules of Cajanus cajan in the Dominican Republic. 16S rRNA gene analysis placed these strains within the genus Bradyrhizobium, being phylogenetically equidistant to several species of this genus. Analysis of the recA and atpD genes showed that the strains isolated belong to a cluster containing the strains Bradyrhizobium ottawaense OO99T, 'Bradyrhizobium americanum' CMVU44 and Bradyrhizobium daqingense CCBAU 15774T, and presented similarity values lower than 96 % for both genes with respect to the strains nodulating C. cajan. DNA-DNA hybridization analysis showed averages of 36, 40 and 39 % relatedness with respect to the representative strains of Bradyrhizobium ottawaense, 'Bradyrhizobium americanum' and Bradyrhizobium daqingense, respectively. Phenotypic characteristics also differed from those of the most closely related species of the genus Bradyrhizobium. Therefore, based on the data obtained in this study, we propose to classify the strains AMBPC1010T (=LMG 29967T=CECT 9227T) and AMBPC1011 into a novel species named Bradyrhizobium cajani sp. nov.

Phylogenies of symbiotic genes of Bradyrhizobium symbionts of legumes of economic and environmental importance in Brazil support the definition of the new symbiovars pachyrhizi and sojae

Bradyrhizobium comprises most tropical symbiotic nitrogen-fixing strains, but the correlation between symbiotic and core genes with host specificity is still unclear. In this study, the phylogenies of the nodY/K and nifH genes of 45 Bradyrhizobium strains isolated from legumes of economic and environmental importance in Brazil (Arachis hypogaea, Acacia auriculiformis, Glycine max, Lespedeza striata, Lupinus albus, Stylosanthes sp. and Vigna unguiculata) were compared to 16S rRNA gene phylogeny and genetic diversity by rep-PCR. In the 16S rRNA tree, strains were distributed into two superclades-B. japonicum and B. elkanii-with several strains being very similar within each clade. The rep-PCR analysis also revealed high intra-species diversity. Clustering of strains in the nodY/K and nifH trees was identical: 39 strains isolated from soybean grouped with Bradyrhizobium type species symbionts of soybean, whereas five others occupied isolated positions. Only one strain isolated from Stylosanthes sp. showed similar nodY/K and nifH sequences to soybean strains, and it also nodulated soybean. Twenty-one representative strains of the 16S rRNA phylogram were selected and taxonomically classified using a concatenated glnII-recA phylogeny; nodC sequences were also compared and revealed the same clusters as observed in the nodY/K and nifH phylograms. The analyses of symbiotic genes indicated that a large group of strains from the B. elkanii superclade comprised the novel symbiovar sojae, whereas for another group, including B. pachyrhizi, the symbiovar pachyrhizi could be proposed. Other potential new symbiovars were also detected. The co-evolution hypotheses is discussed and it is suggested that nodY/K analysis would be useful for investigating the symbiotic diversity of the genus Bradyrhizobium.