Comparative genomic analysis of Bradyrhizobium strains with natural variability in the efficiency of nitrogen fixation, competitiveness, and adaptation to stressful edaphoclimatic conditions

Bradyrhizobium is known for fixing atmospheric nitrogen in symbiosis with agronomically important crops. This study focused on two groups of strains, each containing eight natural variants of the parental strains, Bradyrhizobium japonicum SEMIA 586 (=CNPSo 17) or Bradyrhizobium diazoefficiens SEMIA 566 (=CNPSo 10). CNPSo 17 and CNPSo 10 were used as commercial inoculants for soybean crops in Brazil at the beginning of the crop expansion in the southern region in the 1960s-1970s. Variants derived from these parental strains were obtained in the late 1980s through a strain selection program aimed at identifying elite strains adapted to a new cropping frontier in the central-western Cerrado region, with a higher capacity of biological nitrogen fixation (BNF) and competitiveness. Here, we aimed to detect genetic variations possibly related to BNF, competitiveness for nodule occupancy, and adaptation to the stressful conditions of the Brazilian Cerrado soils. High-quality genome assemblies were produced for all strains. The core genome phylogeny revealed that strains of each group are closely related, as confirmed by high average nucleotide identity values. However, variants accumulated divergences resulting from horizontal gene transfer, genomic rearrangements, and nucleotide polymorphisms. The B. japonicum group presented a larger pangenome and a higher number of nucleotide polymorphisms than the B. diazoefficiens group, possibly due to its longer adaptation time to the Cerrado soil. Interestingly, five strains of the B. japonicum group carry two plasmids. The genetic variability found in both groups is discussed considering the observed differences in their BNF capacity, competitiveness for nodule occupancy, and environmental adaptation.IMPORTANCEToday, Brazil is a global leader in the study and use of biological nitrogen fixation with soybean crops. As Brazilian soils are naturally void of soybean-compatible bradyrhizobia, strain selection programs were established, starting with foreign isolates. Selection searched for adaptation to the local edaphoclimatic conditions, higher efficiency of nitrogen fixation, and strong competitiveness for nodule occupancy. We analyzed the genomes of two parental strains of Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens and eight variant strains derived from each parental strain. We detected two plasmids in five strains and several genetic differences that might be related to adaptation to the stressful conditions of the soils of the Brazilian Cerrado biome. We also detected genetic variations in specific regions that may impact symbiotic nitrogen fixation. Our analysis contributes to new insights into the evolution of Bradyrhizobium, and some of the identified differences may be applied as genetic markers to assist strain selection programs.

Genetic variation in symbiotic islands of natural variant strains of soybean Bradyrhizobium japonicum and Bradyrhizobium diazoefficiens differing in competitiveness and in the efficiency of nitrogen fixation

Soybean is the most important legume cropped worldwide and can highly benefit from the biological nitrogen fixation (BNF) process. Brazil is recognized for its leadership in the use of inoculants and two strains, Bradyrhizobium japonicum CPAC 15 (=SEMIA 5079) and Bradyrhizobium diazoefficiens CPAC 7 (=SEMIA 5080) compose the majority of the 70 million doses of soybean inoculants commercialized yearly in the country. We studied a collection of natural variants of these two strains, differing in properties of competitiveness and efficiency of BNF. We sequenced the genomes of the parental strain SEMIA 566 of B. japonicum, of three natural variants of this strain (S 204, S 340 and S 370), and compared with another variant of this group, strain CPAC 15. We also sequenced the genome of the parental strain SEMIA 586 of B. diazoefficiens, of three natural variants of this strain (CPAC 390, CPAC 392 and CPAC 394) and compared with the genome of another natural variant, strain CPAC 7. As the main genes responsible for nodulation (nod, noe, nol) and BNF (nif, fix) in soybean Bradyrhizobium are located in symbiotic islands, our objective was to identify genetic variations located in this region, including single nucleotide polymorphisms (SNPs) and insertions and deletions (indels), that could be potentially related to their different symbiotic phenotypes. We detected 44 genetic variations in the B. japonicum strains and three in B. diazoefficiens. As the B. japonicum strains have gone through a longer period of adaptation to the soil, the higher number of genetic variations could be explained by survival strategies under the harsh environmental conditions of the Brazilian Cerrado biome. Genetic variations were detected in genes enconding proteins such as a dephospho-CoA kinase, related to the CoA biosynthesis; a glucosamine-fructose-6-phosphate aminotransferase, key regulator of the hexosamine biosynthetic pathway; a LysR family transcriptional regulator related to nodulation genes; and NifE and NifS proteins, directly related to the BNF process. We suggest potential genetic variations related to differences in the symbiotic phenotypes.

Multidisciplinary approaches for studying rhizobium–legume symbioses

The rhizobium–legume symbiosis is a major source of fixed nitrogen (ammonia) in the biosphere. The potential for this process to increase agricultural yield while reducing the reliance on nitrogen-based fertilizers has generated interest in understanding and manipulating this process. For decades, rhizobium research has benefited from the use of leading techniques from a very broad set of fields, including population genetics, molecular genetics, genomics, and systems biology. In this review, we summarize many of the research strategies that have been employed in the study of rhizobia and the unique knowledge gained from these diverse tools, with a focus on genome- and systems-level approaches. We then describe ongoing synthetic biology approaches aimed at improving existing symbioses or engineering completely new symbiotic interactions. The review concludes with our perspective of the future directions and challenges of the field, with an emphasis on how the application of a multidisciplinary approach and the development of new methods will be necessary to ensure successful biotechnological manipulation of the symbiosis.

The Bacteroid Periplasm in Soybean Nodules Is an Interkingdom Symbiotic Space

The functional role of the periplasm of nitrogen-fixing bacteroids has not been determined. Proteins were isolated from the periplasm and cytoplasm of Bradyrhizobium diazoefficiens bacteroids and were analyzed using liquid chromatography tandem mass spectrometry proteomics. Identification of bacteroid periplasmic proteins was aided by periplasm prediction programs. Approximately 40% of all the proteins identified as periplasmic in the B. diazoefficiens genome were found expressed in the bacteroid form of the bacteria, indicating the periplasm is a metabolically active symbiotic space. The bacteroid periplasm possesses many fatty acid metabolic enzymes, which was in contrast to the bacteroid cytoplasm. Amino acid analysis of the periplasm revealed an abundance of phosphoserine, phosphoethanolamine, and glycine, which are metabolites of phospholipid metabolism. These results suggest the periplasm is a unique space and not a continuum with the peribacteroid space. A number of plant proteins were found in the periplasm fraction, which suggested contamination. However, antibodies to two of the identified plant proteins, histone H2A and lipoxygenase, yielded immunogold labeling that demonstrated the plant proteins were specifically targeted to the bacteroids. This suggests that the periplasm is an interkingdom symbiotic space containing proteins from both the bacteroid and the plant.

Proteomic analysis of free-living Bradyrhizobium diazoefficiens: highlighting potential determinants of a successful symbiosis

BACKGROUND

Strain CPAC 7 (=SEMIA 5080) was recently reclassified into the new species Bradyrhizobium diazoefficiens; due to its outstanding efficiency in fixing nitrogen, it has been used in commercial inoculants for application to crops of soybean [Glycine max (L.) Merr.] in Brazil and other South American countries. Although the efficiency of B. diazoefficiens inoculant strains is well recognized, few data on their protein expression are available.

RESULTS

We provided a two-dimensional proteomic reference map of CPAC 7 obtained under free-living conditions, with the successful identification of 115 spots, representing 95 different proteins. The results highlighted the expression of molecular determinants potentially related to symbiosis establishment (e.g. inositol monophosphatase, IMPase), fixation of atmospheric nitrogen (N2) (e.g. NifH) and defenses against stresses (e.g. chaperones). By using bioinformatic tools, it was possible to attribute probable functions to ten hypothetical proteins. For another ten proteins classified as "NO related COG" group, we analyzed by RT-qPCR the relative expression of their coding-genes in response to the nodulation-gene inducer genistein. Six of these genes were up-regulated, including blr0227, which may be related to polyhydroxybutyrate (PHB) biosynthesis and competitiveness for nodulation.

CONCLUSIONS

The proteomic map contributed to the identification of several proteins of B. diazoefficiens under free-living conditions and our approach-combining bioinformatics and gene-expression assays-resulted in new information about unknown genes that might play important roles in the establishment of the symbiosis with soybean.

Comparative genomics of Bradyrhizobium japonicum CPAC 15 and Bradyrhizobium diazoefficiens CPAC 7: elite model strains for understanding symbiotic performance with soybean

BACKGROUND

The soybean-Bradyrhizobium symbiosis can be highly efficient in fixing nitrogen, but few genomic sequences of elite inoculant strains are available. Here we contribute with information on the genomes of two commercial strains that are broadly applied to soybean crops in the tropics. B. japonicum CPAC 15 (=SEMIA 5079) is outstanding in its saprophytic capacity and competitiveness, whereas B. diazoefficiens CPAC 7 (=SEMIA 5080) is known for its high efficiency in fixing nitrogen. Both are well adapted to tropical soils. The genomes of CPAC 15 and CPAC 7 were compared to each other and also to those of B. japonicum USDA 6T and B. diazoefficiens USDA 110T.

RESULTS

Differences in genome size were found between species, with B. japonicum having larger genomes than B. diazoefficiens. Although most of the four genomes were syntenic, genome rearrangements within and between species were observed, including events in the symbiosis island. In addition to the symbiotic region, several genomic islands were identified. Altogether, these features must confer high genomic plasticity that might explain adaptation and differences in symbiotic performance. It was not possible to attribute known functions to half of the predicted genes. About 10% of the genomes was composed of exclusive genes of each strain, but up to 98% of them were of unknown function or coded for mobile genetic elements. In CPAC 15, more genes were associated with secondary metabolites, nutrient transport, iron-acquisition and IAA metabolism, potentially correlated with higher saprophytic capacity and competitiveness than seen with CPAC 7. In CPAC 7, more genes were related to the metabolism of amino acids and hydrogen uptake, potentially correlated with higher efficiency of nitrogen fixation than seen with CPAC 15.

CONCLUSIONS

Several differences and similarities detected between the two elite soybean-inoculant strains and between the two species of Bradyrhizobium provide new insights into adaptation to tropical soils, efficiency of N2 fixation, nodulation and competitiveness.

Proteomic Analysis of Soybean [Glycine max (L.) Merrill] Roots Inoculated with Bradyrhizobium japonicum Strain CPAC 15

This research intended to analyze the expression pattern of proteins in roots of the Brazilian soybean cultivar Conquista when inoculated with Bradyrhizobium japonicum CPAC 15, a strain broadly used in commercial inoculants in Brazil. At ten days after bacterial inoculation, whole-cell proteins were extracted from roots and separated by 2-D gel electrophoresis. Comparative analysis revealed significant changes in the intensity of 37 spots due to the inoculation (17 up-regulated and 20 down-regulated proteins), identified by MALDI-TOF/TOF-TOF. Identified proteins were associated with COG functional categories of information storage and processing, cellular processes and signaling, metabolism, and also in the “poorly characterized” and “not in COG” categories. Among the up-regulated proteins, we identified sucrose synthase (nodulin-100), β-tubulin, rubisco activase, glutathione-S-transferase, a putative heat-shock 70-kDa protein, pyridine nucleotide-disulphideoxidoreductase and a putative transposase. Proteomic analysis allowed for the identification of some putative symbiotic functions and confirmed the main biological processes triggered in the nitrogen-fixing symbiosis with soybean.

Proteomic profiling of Rhizobium tropici PRF 81: identification of conserved and specific responses to heat stress

Background

Rhizobium tropici strain PRF 81 (= SEMIA 4080) has been used in commercial inoculants for application to common-bean crops in Brazil since 1998, due to its high efficiency in fixing nitrogen, competitiveness against indigenous rhizobial populations and capacity to adapt to stressful tropical conditions, representing a key alternative to application of N-fertilizers. The objective of our study was to obtain an overview of adaptive responses to heat stress of strain PRF 81, by analyzing differentially expressed proteins when the bacterium is grown at 28°C and 35°C.

Results

Two-dimensional gel electrophoresis (2DE) revealed up-regulation of fifty-nine spots that were identified by MALDI-TOF/TOF-TOF. Differentially expressed proteins were associated with the functional COG categories of metabolism, cellular processes and signaling, information storage and processing. Among the up-regulated proteins, we found some related to conserved heat responses, such as molecular chaperones DnaK and GroEL, and other related proteins, such as translation factors EF-Tu, EF-G, EF-Ts and IF2. Interestingly, several oxidative stress-responsive proteins were also up-regulated, and these results reveal the diversity of adaptation mechanisms presented by this thermotolerant strain, suggesting a cross-talk between heat and oxidative stresses.

Conclusions

Our data provide valuable protein-expression information relevant to the ongoing genome sequencing of strain PRF 81, and contributes to our still-poor knowledge of the molecular determinants of the thermotolerance exhibited by R. tropici species.

The genetics of symbiotic nitrogen fixation: comparative genomics of 14 rhizobia strains by resolution of protein clusters

The symbiotic relationship between legumes and nitrogen fixing bacteria is critical for agriculture, as it may have profound impacts on lowering costs for farmers, on land sustainability, on soil quality, and on mitigation of greenhouse gas emissions. However, despite the importance of the symbioses to the global nitrogen cycling balance, very few rhizobial genomes have been sequenced so far, although there are some ongoing efforts in sequencing elite strains. In this study, the genomes of fourteen selected strains of the order Rhizobiales, all previously fully sequenced and annotated, were compared to assess differences between the strains and to investigate the feasibility of defining a core ‘symbiome’—the essential genes required by all rhizobia for nodulation and nitrogen fixation. Comparison of these whole genomes has revealed valuable information, such as several events of lateral gene transfer, particularly in the symbiotic plasmids and genomic islands that have contributed to a better understanding of the evolution of contrasting symbioses. Unique genes were also identified, as well as omissions of symbiotic genes that were expected to be found. Protein comparisons have also allowed the identification of a variety of similarities and differences in several groups of genes, including those involved in nodulation, nitrogen fixation, production of exopolysaccharides, Type I to Type VI secretion systems, among others, and identifying some key genes that could be related to host specificity and/or a better saprophytic ability. However, while several significant differences in the type and number of proteins were observed, the evidence presented suggests no simple core symbiome exists. A more abstract systems biology concept of nitrogen fixing symbiosis may be required. The results have also highlighted that comparative genomics represents a valuable tool for capturing specificities and generalities of each genome.

Proteomics reveals differential expression of proteins related to a variety of metabolic pathways by genistein-induced Bradyrhizobium japonicum strains

The rhizobia-legume symbiosis requires a coordinated molecular interaction between the symbionts, initiated by seed and root exudation of several compounds, mainly flavonoids, that trigger the expression of nodulation genes in the bacteria. Since the role of flavonoids seems to be broader than the induction of nodulation genes, we aimed at characterizing genistein-induced proteins of Bradyrhizobium japonicum CPAC 15 (=SEMIA 5079), used in commercial soybean inoculants in Brazil, and of two genetically related strains grown in vitro. Whole-cell proteins were extracted both from induced (1 μM genistein) and from non-induced cultures of the three strains, and separated by two-dimensional electrophoresis. Spot profiles were compared between the two conditions and selected spots were excised and identified by mass spectrometry. Forty-seven proteins were significantly induced by genistein, including several hypothetical proteins, the cytoplasmic flagellar component FliG, periplasmic ABC transporters, a protein related to biosynthesis of exopolysaccharides (ExoN), and proteins involved in redox-state maintenance. Noteworthy was the induction of the PhyR-σ(EcfG) regulon, recently demonstrated to be involved in the symbiotic efficiency of, and general stress response in B. japonicum. Our results confirm that the role of flavonoids, such as genistein, can go far beyond the expression of nodulation-related proteins in B. japonicum.

Proteomics Analysis of the Effects of Cyanate on Chromobacterium violaceum Metabolism

Chromobacterium violaceum is a gram-negative betaproteobacterium that has been isolated from various Brazilian ecosystems. Its genome contains the cyn operon, which gives it the ability to metabolize highly toxic cyanate into ammonium and carbon dioxide. We used a proteomics approach to investigate the effects of cyanate on the metabolism of this bacterium. The proteome of cells grown with and without cyanate was compared on 2-D gels. Differential spots were digested and identified by mass spectrometry. The bacterium was able to grow at concentrations of up to 1 mM cyanate. Eighteen spots were differentially expressed in the presence of cyanate, of which 16 were downregulated and only two were upregulated. An additional 12 spots were detected only in extracts of cells unexposed to cyanate, and one was expressed only by the exposed cells. Fourteen spots were identified, corresponding to 13 different proteins. We conclude that cyanate promotes expression of enzymes that combat oxidative stress and represses enzymes of the citric acid cycle, strongly affecting the energetic metabolism of the cell. Other proteins that were under-expressed in bacteria exposed to cyanate are involved in amino-acid metabolism or are hypothetical proteins, demonstrating that cyanate also affects expression of genes that are not part of the cyn operon.