Development and application of a multilocus sequence analysis method for the identification of genotypes within genus Bradyrhizobium and for establishing nodule occupancy of soybean (Glycine max L. Merr)

Host-imposed control mechanisms in legume-rhizobia symbiosis

Legumes are ecologically and economically important plants that contribute to nutrient cycling and agricultural sustainability, features tied to their intimate symbiosis with nitrogen-fixing rhizobia. Rhizobia vary dramatically in quality, ranging from highly growth-promoting to non-beneficial; therefore, legumes must optimize their symbiosis with rhizobia through host mechanisms that select for beneficial rhizobia and limit losses to non-beneficial strains. In this Perspective, we examine the considerable scientific progress made in decoding host control over rhizobia, empirically examining both molecular and cellular mechanisms and their effects on rhizobia symbiosis and its benefits. We consider pre-infection controls, which require the production and detection of precise molecular signals by the legume to attract and select for compatible rhizobia strains. We also discuss post-infection mechanisms that leverage the nodule-level and cell-level compartmentalization of symbionts to enable host control over rhizobia development and proliferation in planta. These layers of host control each contribute to legume fitness by directing host resources towards a narrowing subset of more-beneficial rhizobia.

Genetic variation and genetic complexity of nodule occupancy in soybean inoculated with USDA110 and USDA123 rhizobium strains

BACKGROUND

Symbiotic nitrogen fixation differs among Bradyrhizobium japonicum strains. Soybean inoculated with USDA123 has a lower yield than strains known to have high nitrogen fixation efficiency, such as USDA110. In the main soybean-producing area in the Midwest of the United States, USDA123 has a high nodule incidence in field-grown soybean and is competitive but inefficient in nitrogen fixation. In this study, a high-throughput system was developed to characterize nodule number among 1,321 Glycine max and 69 Glycine soja accessions single inoculated with USDA110 and USDA123.

RESULTS

Seventy-three G. max accessions with significantly different nodule number of USDA110 and USDA123 were identified. After double inoculating 35 of the 73 accessions, it was observed that PI189939, PI317335, PI324187B, PI548461, PI562373, and PI628961 were occupied by USDA110 and double-strain nodules but not by USDA123 nodules alone. PI567624 was only occupied by USDA110 nodules, and PI507429 restricted all strains. Analysis showed that 35 loci were associated with nodule number in G. max when inoculated with strain USDA110 and 35 loci with USDA123. Twenty-three loci were identified in G. soja when inoculated with strain USDA110 and 34 with USDA123. Only four loci were common across two treatments, and each locus could only explain 0.8 to 1.5% of phenotypic variation.

CONCLUSIONS

High-throughput phenotyping systems to characterize nodule number and occupancy were developed, and soybean germplasm restricting rhizobium strain USDA123 but preferring USDA110 was identified. The larger number of minor effects and a small few common loci controlling the nodule number indicated trait genetic complexity and strain-dependent nodulation restriction. The information from the present study will add to the development of cultivars that limit USDA123, thereby increasing nitrogen fixation efficiency and productivity.

Phylogenetic distribution and evolutionary dynamics of nod and T3SS genes in the genus Bradyrhizobium

Bradyrhizobium are abundant soil bacteria and the major symbiont of legumes. The recent availability of Bradyrhizobium genome sequences provides a large source of information for analysis of symbiotic traits. In this study, we investigated the evolutionary dynamics of the nodulation genes (nod) and their relationship with the genes encoding type III secretion systems (T3SS) and their effectors among bradyrhizobia. Based on the comparative analysis of 146 Bradyrhizobium genome sequences, we identified six different types of T3SS gene clusters. The two predominant cluster types are designated RhcIa and RhcIb and both belong to the RhcI-T3SS family previously described in other rhizobia. They are found in 92/146 strains, most of them also containing nod genes. RhcIa and RhcIb gene clusters differ in the genes they carry: while the translocon-encoding gene nopX is systematically found in strains containing RhcIb, the nopE and nopH genes are specifically conserved in strains containing RhcIa, suggesting that these last two genes might functionally substitute nopX and play a role related to effector translocation. Phylogenetic analysis suggests that bradyrhizobia simultaneously gained nod and RhcI-T3SS gene clusters via horizontal transfer or subsequent vertical inheritance of a symbiotic island containing both. Sequence similarity searches for known Nop effector proteins in bradyrhizobial proteomes revealed the absence of a so-called core effectome, i.e. that no effector is conserved among all Bradyrhizobium strains. However, NopM and SUMO proteases were found to be the main effector families, being represented in the majority of the genus. This study indicates that bradyrhizobial T3SSs might play a more significant symbiotic role than previously thought and provides new candidates among T3SS structural proteins and effectors for future functional investigations.

Identification of soybean Bradyrhizobium strains used in commercial inoculants in Brazil by MALDI-TOF mass spectrometry

Biological nitrogen fixation (BNF) with the soybean crop probably represents the major sustainable technology worldwide, saving billions of dollars in N fertilizers and decreasing water pollution and the emission of greenhouse gases. Accordingly, the identification of strains occupying nodules under field conditions represents a critical step in studies that are aimed at guaranteeing increased BNF contribution. Current methods of identification are mostly based on serology, or on DNA profiles. However, the production of antibodies is restricted to few laboratories, and to obtain DNA profiles of hundreds of isolates is costly and time-consuming. Conversely, the matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) MS technique might represent a golden opportunity for replacing serological and DNA-based methods. However, MALDI-TOF databases of environmental microorganisms are still limited, and, most importantly, there are concerns about the discrimination of protein profiles at the strain level. In this study, we investigated four soybean rhizobial strains carried in commercial inoculants used in over 35 million hectares in Brazil and also in other countries of South America and Africa. A supplementary MALDI-TOF database with the protein profiles of these rhizobial strains was built and allowed the identification of unique profiles statistically supported by multivariate analysis and neural networks. To test this new database, the nodule occupancy by Bradyrhizobium strains in symbiosis with soybean was characterized in a field experiment and the results were compared with serotyping of bacteria by immuno-agglutination. The results obtained by both techniques were highly correlated and confirmed the viability of using the MALDI-TOF MS technique to effectively distinguish bacteria at the strain level.

Elucidation of the Genome of Bradyrhizobium sp. Strain USDA 3456, a Historic Agricultural Diazotroph from Cowpea (Vigna unguiculata)

Bradyrhizobium sp. strain USDA 3456 is a historic strain from the United States Department of Agriculture (USDA) Agricultural Research Service (ARS) National Rhizobium Germplasm Collection isolated from Vigna unguiculata (cowpea) in 1966. Strain USDA 3456 has been utilized in global agricultural applications, including improving soil nitrogen fertility. The draft genome sequence here provides a genetic reference of a novel diazotroph.

Bradyrhizobium sp. strain USDA 3456 is a historic strain from the United States Department of Agriculture (USDA) Agricultural Research Service (ARS) National Rhizobium Germplasm Collection isolated from Vigna unguiculata (cowpea) in 1966. Strain USDA 3456 has been utilized in global agricultural applications, including improving soil nitrogen fertility. The draft genome sequence here provides a genetic reference of a novel diazotroph.

Phylogenomic Analyses of Bradyrhizobium Reveal Uneven Distribution of the Lateral and Subpolar Flagellar Systems, Which Extends to Rhizobiales

Dual flagellar systems have been described in several bacterial genera, but the extent of their prevalence has not been fully explored. Bradyrhizobium diazoefficiens USDA 110^T possesses two flagellar systems, the subpolar and the lateral flagella. The lateral flagellum of Bradyrhizobium displays no obvious role, since its performance is explained by cooperation with the subpolar flagellum. In contrast, the lateral flagellum is the only type of flagella present in the related Rhizobiaceae family. In this work, we have analyzed the phylogeny of the Bradyrhizobium genus by means of Genome-to-Genome Blast Distance Phylogeny (GBDP) and Average Nucleotide Identity (ANI) comparisons of 128 genomes and divided it into 13 phylogenomic groups. While all the Bradyrhizobium genomes encode the subpolar flagellum, none of them encodes only the lateral flagellum. The simultaneous presence of both flagella is exclusive of the B. japonicum phylogenomic group. Additionally, 292 Rhizobiales order genomes were analyzed and both flagellar systems are present together in only nine genera. Phylogenetic analysis of 150 representative Rhizobiales genomes revealed an uneven distribution of these flagellar systems. While genomes within and close to the Rhizobiaceae family only possess the lateral flagellum, the subpolar flagellum is exclusive of more early-diverging families, where certain genera also present both flagella.

A Proteomic Network for Symbiotic Nitrogen Fixation Efficiency in Bradyrhizobium elkanii

Rhizobia colonize legumes and reduce N₂ to NH₃ in root nodules. The current model is that symbiotic rhizobia bacteroids avoid assimilating this NH₃. Instead, host legume cells form glutamine from NH₃, and the nitrogen is returned to the bacteroid as dicarboxylates, peptides, and amino acids. In soybean cells surrounding bacteroids, glutamine also is converted to ureides. One problem for soybean cultivation is inefficiency in symbiotic N₂ fixation, the biochemical basis of which is unknown. Here, the proteomes of bacteroids of Bradyrhizobium elkanii USDA76 isolated from N₂ fixation-efficient Peking and -inefficient Williams 82 soybean nodules were analyzed by mass spectrometry. Nearly half of the encoded bacterial proteins were quantified. Efficient bacteroids produced greater amounts of enzymes to form Nod factors and had increased amounts of signaling proteins, transporters, and enzymes needed to generate ATP to power nitrogenase and to acquire resources. Parallel investigation of nodule proteins revealed that Peking had no significantly greater accumulation of enzymes needed to assimilate NH₃ than Williams 82. Instead, efficient bacteroids had increased amounts of enzymes to produce amino acids, including glutamine, and to form ureide precursors. These results support a model for efficient symbiotic N₂ fixation in soybean where the bacteroid assimilates NH₃ for itself.

Bradyrhizobium sacchari sp. nov., a legume nodulating bacterium isolated from sugarcane roots

Members of the genus Bradyrhizobium are well-known as nitrogen-fixing microsymbionts of a wide variety of leguminous species, but they have also been found in different environments, notably as endophytes in non-legumes such as sugarcane. This study presents a detailed polyphasic characterization of four Bradyrhizobium strains (type strain BR 10280T), previously isolated from roots of sugarcane in Brazil. 16S rRNA sequence analysis, multilocus sequence analysis (MLSA) and analysis of the 16S-23S rRNA internal transcribed spacer showed that these strains form a novel clade close to, but different from B. huanghuaihaiense strain CCBAU 23303T. Average nucleotide identity (ANI) analyses confirmed that BR 10280T represents a novel species. Phylogenetic analysis based on nodC gene sequences also placed the strains close to CCBAU 23303T, but different from this latter strain, the sugarcane strains did not nodulate soybean, although they effectively nodulated Vigna unguiculata, Cajanus cajan and Macroptilium atropurpureum. Physiological traits are in agreement with the placement of the strains in the genus Bradyrhizobium as a novel species for which the name Bradyrhizobium sacchari sp. nov. is proposed.

African origin of Bradyrhizobium populations nodulating Bambara groundnut (Vigna subterranea L. Verdc) in Ghanaian and South African soils

Flavonoids secreted by legumes play a major role as signal molecules for attracting compatible rhizobia. The aim of this study was to assess and understand the diversity of microsymbionts nodulating Bambara groundnut (Vigna subterranea L. Verdc.) landraces of different seedcoat colours using restriction fragment length polymorphism and phylogenetic analysis. Seedcoat pigmentation of landraces had effect on the diversity of microsymbionts of Bambara groundnut. Even when planted together in one hole, nodulating bradyrhizobia clustered differently. For example, 16S rDNA-RFLP typing of rhizobial samples TUTVSBLM.I, TUTVSCRM.I and TUTVSRDM.I originating respectively from Black, Cream and Red landraces that were co-planted in the same hole at Manga in the Sudano-sahelian savanna, as well as TUTVSCRK.I and TUTVSRDK.I respectively from Cream and Red landraces co-planted at Kpalisogu in the Guinea savanna, revealed different 16S rDNA- RFLP types. Phylogenetic analysis of 16S rDNA, glnII, recA and atpD sequences showed that Vigna subterranea was nodulated specifically by a diverse group of Bradyrhizobium species (e.g. Bradyrhizobium vignae, and a novel group of Bradyrhizobium spp.) in soils from Ghana and South Africa. The recA gene phylogeny showed incongruency with the other housekeeping genes, indicating the possibility of lateral gene transfer and/or recombination events. The grouping of isolates according to symbiotic gene (nifH and nodD) phylogenies revealed inter- and intra-specific symbiotic plasmid transfer and different evolutionary history. The results also showed that a cropping history and physico-chemical environment of soils increased bradyrhizobial diversity in Ghana and South Africa.

Swimming performance of Bradyrhizobium diazoefficiens is an emergent property of its two flagellar systems

Many bacterial species use flagella for self-propulsion in aqueous media. In the soil, which is a complex and structured environment, water is found in microscopic channels where viscosity and water potential depend on the composition of the soil solution and the degree of soil water saturation. Therefore, the motility of soil bacteria might have special requirements. An important soil bacterial genus is Bradyrhizobium, with species that possess one flagellar system and others with two different flagellar systems. Among the latter is B. diazoefficiens, which may express its subpolar and lateral flagella simultaneously in liquid medium, although its swimming behaviour was not described yet. These two flagellar systems were observed here as functionally integrated in a swimming performance that emerged as an epistatic interaction between those appendages. In addition, each flagellum seemed engaged in a particular task that might be required for swimming oriented toward chemoattractants near the soil inner surfaces at viscosities that may occur after the loss of soil gravitational water. Because the possession of two flagellar systems is not general in Bradyrhizobium or in related genera that coexist in the same environment, there may be an adaptive tradeoff between energetic costs and ecological benefits among these different species.

MALDI-TOF mass spectrometry as a tool for differentiation of Bradyrhizobium species: application to the identification of Lupinus nodulating strains

Genus Bradyrhizobium includes slow growing bacteria able to nodulate different legumes as well as species isolated from plant tumours. The slow growth presented by the members of this genus and the phylogenetic closeness of most of its species difficults their identification. In the present work we applied for the first time Matrix-Assisted Laser Desorption Ionization-Time-of-Flight Mass Spectrometry (MALDI-TOF MS) to the analysis of Bradyrhizobium species after the extension of MALDI Biotyper 2.0 database with the currently valid species of this genus. With this methodology it was possible to identify strains belonging to phylogenetically closely related species of genus Bradyrhizobium allowing the discrimination among species with rrs gene identities higher than 99%. The application of MALDI-TOF MS to strains isolated from nodules of different Lupinus species in diverse geographical locations allowed their correct identification when comparing with the results of rrs gene and ITS analyses. The nodulation of Lupinus gredensis, an endemic species of the west of Spain, by B. canariense supports the European origin of this species.

Pathogen typing in the genomics era: MLST and the future of molecular epidemiology

Multi-locus sequence typing (MLST) is a high-resolution genetic typing approach to identify species and strains of pathogens impacting human health, agriculture (animals and plants), and biosafety. In this review, we outline the general concepts behind MLST, molecular approaches for obtaining MLST data, analytical approaches for MLST data, and the contributions MLST studies have made in a wide variety of areas. We then look at the future of MLST and their relative strengths and weaknesses with respect to whole genome sequence typing approaches that are moving into the research arena at an ever-increasing pace. Throughout the paper, we provide exemplar references of these various aspects of MLST. The literature is simply too vast to make this review comprehensive, nevertheless, we have attempted to include enough references in a variety of key areas to introduce the reader to the broad applications and complications of MLST data.

Microevolution of symbiotic Bradyrhizobium populations associated with soybeans in east North America

Microevolution and origins of Bradyrhizobium populations associated with soybeans at two field sites (A and B, 280 km apart in Canada) with contrasting histories of inoculation was investigated using probabilistic analyses of six core (housekeeping) gene sequences. These analyses supported division of 220 isolates in five lineages corresponding either to B. japonicum groups 1 and 1a or to one of three novel lineages within the genus Bradyrhizobium. None of the isolates from site A and about 20% from site B (the only site with a recent inoculation history) were attributed to inoculation sources. The data suggest that most isolates were of indigenous origin based on sequence analysis of 148 isolates of soybean-nodulating bacteria from native legumes (Amphicarpaea bracteata and Desmodium canadense). Isolates from D. canadense clustered with B. japonicum group 1, whereas those from A. bracteata were placed in two novel lineages encountered at soybean field sites. One of these novel lineages predominated at soybean sites and exhibited a significant clonal expansion likely reflecting selection by the plant host. Homologous recombination events detected in the 35 sequence types from soybean sites had an effect on genetic diversification that was approximately equal to mutation. Interlineage transfer of core genes was infrequent and mostly attributable to gyrB that had a history of frequent recombination. Symbiotic gene sequences (nodC and nifH) of isolates from soybean sites and native legumes clustered in two lineages corresponding to B. japonicum and B. elkani with the inheritance of these genes appearing predominantly by vertical transmission. The data suggest that soybean-nodulating bacteria associated with native legumes represent a novel source of ecologically adapted bacteria for soybean inoculation.