Open Access

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

  • Arthur Fernandes Siqueira1, 2,
  • Ernesto Ormeño-Orrillo3,
  • Rangel Celso Souza4,
  • Elisete Pains Rodrigues5,
  • Luiz Gonzaga Paula Almeida4,
  • Fernando Gomes Barcellos5,
  • Jesiane Stefânia Silva Batista6,
  • Andre Shigueyoshi Nakatani2,
  • Esperanza Martínez-Romero3,
  • Ana Tereza Ribeiro Vasconcelos4 and
  • Mariangela Hungria1, 2Email author
Contributed equally
BMC Genomics201415:420

DOI: 10.1186/1471-2164-15-420

Received: 2 March 2014

Accepted: 20 May 2014

Published: 3 June 2014

Abstract

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.

Keywords

Symbiosis Nodulation Nitrogen fixation Competitiveness Secretion systems Horizontal gene transfer Membrane transporters Surface polysaccharides Secondary metabolism Phytohormone synthesis

Background

Soybean [Glycine max (L.) Merr.] is the most important legume cropped worldwide, and the expansion of this crop since the middle of the 20th century—particularly in tropical areas with soils deficient in nitrogen (N)—is attributable to its high capacity in fixing atmospheric nitrogen (N2) in symbioses with Bradyrhizobium[1]. However, despite the economic and environment importance of this crop, few genomic studies have been reported on elite Bradyrhizobium strains. This represents an important opportunity to understand features correlated with high efficiency of nitrogen fixation with the legume. A pioneering genomic study was performed with B. japonicum strain USDA 110T[2], reclassified now as the type strain of B. diazoefficiens[3] and used in commercial inoculants in the United States, in African and other countries [1]. This research highlighted intriguing information about the symbiosis island and high numbers of transposase-encoding and unknown genes [2]. Since then, the genome of B. japonicum USDA 6T[4] was presented and emphasized evolutionary aspects, possibly because the strain is not outstanding in fixing N2. A few other draft genomes of soybean Bradyrhizobium strains are now available, but their descriptions also highlighted phylogenetic aspects (e.g. [5]). Altogether, the results clearly point out that other genomes, with an emphasis on elite inoculant strains, still have to be sequenced to get a better understanding of the molecular basis of the “perfect symbiosis” involving soybean and Bradyrhizobium[1].

Brazil is the second largest producer of soybean, and is expected to become the largest very soon. Biological fixation of N2 has always been a priority in the Brazilian production system. Research programs have focused both on Bradyrhizobium-strain selection and identification of plant genotypes of superior symbiotic performance with the aim of obtaining high yields without N-fertilizer inputs [1, 6, 7]. Originally, Brazilian soils were free of soybean Bradyrhizobium[8] and commercial foreign inoculants were first used in the 20th century. Searches for locally adapted strains were started immediately and still continue [1, 7, 9]. In the absence of natural biodiversity, these strain-selection programs were, and are, searches for variant genotypes of the introduced Bradyrhizobium strains, superior in N2-fixation capacity and competitiveness obtained mainly by reisolation after long periods of adaptation to local stressful environmental and soil conditions [6, 9, 10]. Two elite strains have emerged from this strategy, B. japonicum CPAC 15 (=SEMIA 5079) and B. diazoefficiens CPAC 7 (=SEMIA 5080) [9, 11]. Interestingly, either due to the adaptation to the edaphic climatic conditions, or to the selection for better symbiotic performance, several morphological, physiological and genetic differences have been reported when variant and parental strains are compared [10, 1216]. In addition, large differences in symbiotic performance, including N2-fixation efficiency and competitiveness in both greenhouse controlled conditions and field experiments have also been detected between parental and variant strains [17, 18].

CPAC 15 (=SEMIA 5079) is a natural variant derived from SEMIA 566, a strain used in inoculants in the late 1960s that belongs to the same serogroup as USDA 123, which is considered as the most competitive serogroup in the United States [1, 6]. CPAC 15 was selected for higher capacity of N2 fixation than the parent, and it has been broadly used in commercial inoculants in Brazil since 1992 [9, 11]. This strain is the most competitive of the four commercial strains used in Brazil, and has been detected in every soil cropped with soybean [1, 6, 19, 20]. In the Midwestern United States, strains belonging to serogroup 123 have been found to occupy 60% to 80% of the nodules formed (e.g. [6, 21, 22]); even higher nodule occupancies have been reported in Brazil [20, 23].

CPAC 7 (=SEMIA 5080) is a natural variant of strain CB 1809 (=SEMIA 586, =3I1b136, =TAL 379; =USDA 136; the latter being a subculture of USDA 122 [24, 25]). CPAC 7 has been selected for higher efficiency of N2 fixation and higher adaptation to tropical soils than the parent [9, 11]. The strain has been used in commercial inoculants in Brazil since 1992 [11]. CPAC 7 is more efficient in fixing N2, but less competitive, than CPAC 15 [6, 17, 20].

Today, Brazil produces over 27 million doses of soybean inoculants per year, including exports to some South American and African countries, the great majority carrying strains CPAC 15 and CPAC 7. In this study, we report the genomes of these two strains, highlighting similarities and differences that may be related to their adaptation to tropical soils, efficiency of N2 fixation and competitiveness. We also compare the two genomes with those of the type strains of B. japonicum and B. diazoefficiens, and some other strains of interest.

Results and discussion

General characteristics of CPAC 15 and CPAC 7 genomes

B. japonicum (Bj) strain CPAC 15 has a 9,582,287-bp genome composed of one circular chromosome with two ribosomal operons and a G + C content of 63.54% (Table 1, Figure 1). Annotation predicted 4,203 ORFs with assigned function (48.6%) and 4,445 hypothetical ORFs (51.4%), according to the criteria established and described in the methods section; in total, 8,648 genes were predicted. A high-quality draft genome with 20-fold coverage and distributed in 13 contigs indicated that B. diazoefficiens (Bd) strain CPAC 7 has a genome size estimated at 9,085,545 bp. No plasmid replication genes were detected, indicating that CPAC 7 does not possess plasmids and that its genome is composed of a single replicon. This result is in agreement with the reported absence of plasmids in Bd strain USDA 122 [26], the parent strain of CPAC 7. In CPAC 7 there was only one ribosomal operon and the G + C content was 63.98% (Table 1, Figure 1). Annotation predicted 4,147 ORFs with assigned functions (50.2%) and 4,084 hypothetical ORFs (49.8%), with 8,231 total predicted genes. General features and statistics of both genomes are presented in Table 1.
Table 1

General information about the genomes of B. japonicum strain CPAC 15 and B. diazoefficiens strain CPAC 7

Characteristic

CPAC 15

CPAC 7

Genome coverage

20-fold

20-fold

Circular chromosome

1

1

Size (bp)

9,582,287

9,085,545

G + C content (%)

63.54

63.98

Number of contigs

1

13

Coding region (% of genome size)

82

83

Average ORF size (bp)

909

927

Number of ORFs with assigned functions

4,203

4,147

Number of hypothetical ORFs

4,445

4,084

rRNA operons

2

1

tRNA genes

51

52

Total number of genes

8,648

8,231

https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-420/MediaObjects/12864_2014_Article_6142_Fig1_HTML.jpg
Figure 1

Representation of the chromosomes of B. japonicum CPAC 15 and B. diazoefficiens CPAC 7. Circles from innermost to outermost indicate: coordinates in Mb, GC skew, GC content, ribosomal operons, tRNA genes, genomic islands (symbiosis islands in black), BLASTN matches to B. japonicum USDA 6T (for CPAC 15) or B. diazoefficiens USDA 110T (for CPAC 7), BLASTN matches to CPAC 7 (for CPAC 15) or CPAC 15 (for CPAC 7), BLASTN matches to USDA 110T (for CPAC 15) or USDA 6T (for CPAC 7). Legend indicates BLASTN percentage identity by color.

Despite many reports of two groups within B. japonicum, it was only recently that a new species, B. diazoefficiens, was described, which includes strains previously classified as B. japonicum Group Ia [3]. Average nucleotide identity (ANI) values are higher within each species than between them [3] and this is clearly observed in the identity percentage of the BLASTN matches when the genomes are compared (three outermost circles in Figure 1). Now we report that the genome of CPAC 7 is smaller than that of Bj CPAC 15 by approximately 500,000 bp. The draft genomes now available for both species (Additional file 1: Table S1) show consistently that strains belonging to B. japonicum have larger genomes than those of B. diazoefficiens.

Considering the ten organisms with the highest genomic similarities determined by the Kyoto Encyclopedia of Genes and Genomes (KEGG) system, 80.9% of the ORFs of CPAC 15 had highest similarity with Bj USDA 6T, 8.2% with Bd USDA 110T, 1.53% with Bradyrhizobium sp. S23321 and less than 0.4% with other organisms. In relation to CPAC 7, 81.3% of its ORFs were highest in similarity with USDA 110T, 5.73% with USDA 6T, 2.32% with S23321 and less than 0.5% for the other organisms (Additional file 1: Table S2). These results are in agreement with known taxonomic affiliations of these strains [3].

Structural genome comparisons

Bj CPAC 15 and Bd CPAC 7 chromosomes were syntenic when compared with each other and with those of Bj USDA 6T and Bd USDA 110T; however, genome rearrangements involving inversions of large genome regions, and also inversions or translocations of small regions were observed (Figure 2). As expected, synteny was higher when genomes of strains from the same species were compared; however, Bj USDA 6T showed a large inversion around the chromosome terminus in comparison to Bj CPAC 15. When B. japonicum and B. diazoefficiens strains were compared, a rearrangement around the chromosome-replication origin (ori) that shifted the position of a region containing the symbiosis island was observed. Although the symbiosis island is located at a different replichore in each species, the orientation relative to the ori was preserved, suggesting that no effects on symbiotic gene expression could be attributed to this rearrangement.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-420/MediaObjects/12864_2014_Article_6142_Fig2_HTML.jpg
Figure 2

Structural genome comparisons between the genomes of B. diazoefficiens strains USDA 110 T and CPAC 7, and B. japonicum strains CPAC 15 and USDA 6 T . Red and blue lanes connecting genomes represent syntenic blocks in direct and reverse orientations, respectively. Positions of the replication origin (ori), replication terminus (ter) and symbiosis island (SI) are indicated.

In Bradyrhizobium, genomic plasticity involving rearrangements was previously reported when species were compared [27, 28]. Here we observed that rearrangements occur also at the strain level. In the genus Rhizobium, rearrangements are common but involve mainly plasmid replicons, whereas chromosomes are more stable [29, 30]. As genome plasticity may be related to adaptation [31], it is not surprising that Bradyrhizobium genomes also displayed rearrangements even when plasmids are relatively uncommon in this genus [26].

Similarities and differences between protein-coding genes of B. japonicum and B. diazoefficiens

Overlaps and differences in genes of Bj CPAC 15, Bj USDA 6T, Bd CPAC 7 and Bd USDA 110T are depicted in a Venn diagram (Figure 3). The common core was composed of 5,770 genes, and included those related to basic cell survival, as ribosomal protein genes, chaperones, DNA-replication and -repair genes, hundreds of ORFs related to energy metabolism, including oxidoreductases, acyltransferases, dehydrogenases, hydrolases, cytochromes, hundreds related to carbon (C) metabolism and transporters of several classes. Other genes included in the common core had roles in nodulation and N2 fixation, and some were hydrogenase genes, genes related to chemotaxis and conjugal transfer, genes related to type I to VI secretion systems, interesting genes such as beta-lactamases, those linked to the metabolism of arsenate and penicillin, many transposases and, as expected, several hundreds encoding hypothetical proteins.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-420/MediaObjects/12864_2014_Article_6142_Fig3_HTML.jpg
Figure 3

Venn diagram depicting overlaps and differences in orthologous genes between B. japonicum strains USDA 6 T and CPAC 15 and B. diazoefficiens strains USDA 110 T and CPAC 7.

Bj CPAC 15 and USDA 6T shared 7,299 genes, whereas Bd CPAC 7 and USDA 110T shared 7,002 genes (Figure 3). In CPAC 15 and USDA 6T, 1,041 and 1,135 genes were exclusive to each strain, respectively. In Bj CPAC 15, 73.1% of the exclusive ORFs were classified as hypothetical (73.1%) and 5.3% as transposases; in addition, there were transporters, ORFs related to secretion systems and transcriptional regulators. In Bj USDA 6T, no functions could be attributed to 93.5% of its unique ORFs and 42 of them (3.7%) encoded transposases. In B. diazoefficiens, 831 and 867 genes were exclusive to CPAC 7 and USDA 110T, respectively; we must remember that these two genomes were smaller than those of B. japonicum (Additional file 1: Table S1). In CPAC 7, 74.6% of the exclusive ORFs were hypothetical and 5.53% encoded transposases, the remaining being represented by several classes of genes including transporters and secretion systems. In USDA 110T, 70.7% of the unique ORFs represented hypothetical genes, there were 88 transposase genes (7.6%), and 2.43% encoded ABC transporters. In summary, considering the whole genome of CPAC 15, 12.1% was represented by hypothetical genes exclusive to this strain. These ORFs ranged from 70 to 3,270 bp, with an average of 194 bp, and they were distributed throughout the genome. Similarly, in CPAC 7 the hypothetical exclusive genes represented 10.1% of the genome, varied in size from 70 to 2,278 bp, with an average of 190 bp, and were also broadly distributed. Altogether, this information highlights our still very poor knowledge of Bradyrhizobium and raises the intriguing question of what features are encoded by this high percentage of hypothetical genes.

The majority of the hypothetical ORFs in Bj CPAC 15 and Bd CPAC 7 had orthologs in at least one of the related genomes, and several had conserved domains and/or had their expression confirmed in other experiments. Some conserved domains allowed assignments in KEGG functional classes. In CPAC 15, the most representative category was that of amino acid metabolism (15.95%), followed by the xenobiotic biodegradation and metabolism (13.24%). In the case of CPAC 7, the xenobiotics category had a higher percentage (18.96%) than in CPAC 15.

Twenty-nine hypothetical proteins of CPAC 15 have been previously reported in its proteomic map [32], and now, with the complete genome, it was possible to localize them in the genome. None was found in the symbiosis island and all except one (BS07054) had orthologs in CPAC 7, USDA 110T and USDA 6T (Additional file 1: Table S3). The proteomic reference map of Bd CPAC 7 also presented a significant number of hypothetical sequences (unpublished) that will now be properly localized in the genome.

Carbohydrate and amino acid metabolism

Predicted genes of Bj CPAC 15 and Bd CPAC 7 fit into 18 KEGG functional classes (Additional file 2: Figure S1). The highest percentages of ORFs were in classes of metabolism of amino acids (14.71 and 15.79% for CPAC 15 and CPAC 7, respectively) and carbohydrates (14.43 and 15.71%, for CPAC 15 and CPAC 7, respectively). The same results were obtained using the COG functional classification (data not shown). Although in RAST a smaller percentage of ORFs was assigned to functional classes, the same two categories were confirmed as the most abundant; in addition, the third class including more ORFs in both strains was that of cofactors, vitamins, prosthetic groups and pigments, followed by fatty acids, lipids and isoprenoids (Additional file 2: Figure S2).

The high number of ORFs of both strains in amino acid and carbohydrate metabolism might be related to their ability to survive in soils of low C and N content. The Brazilian Cerrados is an edaphic type of savanna occupying 207 million hectares; C content in these soils is low, ranging from 0.5 to 2.3 mg/g in the superficial layer and can decrease drastically with agricultural cultivation [33]. In a previous study, we reported that rhizobia adapted to soils of low C content may use a broader range of sources of C to survive than those in soils rich in C [34].

It is noteworthy that, in the amino acid metabolism category, CPAC 7 had a higher percentage of ORFs than did CPAC 15. It is known that a complex amino acid cycle is essential for symbiotic N2 fixation, involving control via amino acid exchange between the symbiotic partners, providing selective pressure for the evolution of mutualism [35, 36]. We may thus hypothesize that the higher efficiency of N2 fixation of CPAC 7 could be related to a larger number of genes involved in the metabolism of amino acids. It has been reported that the synthesis of branched chain amino acids (BCAAs) is switched off in bacteroids of R. leguminosarum, rendering them as functional auxotrophs depending on the host for their supply [37]. In vitro growth tests show that both Bj USDA 6T[3] and CPAC 15 (unpublished) use poorly the BCAAs leucine and isoleucine as sole N sources in comparison to Bd CPAC 7 and USDA 110T. Some genes encoding enzymes of the common and specific pathways for BCAA degradation were missing in both CPAC 15 and USDA 6T, therefore both strains may have to metabolize leucine and isoleucine using other less-efficient pathway. These observations may be linked to the inferior symbiotic N2-fixation efficiency of B. japonicum in comparison to B. diazoefficiens if leucine and isoleucine are as important in the bradyrhizobia-soybean symbioses as they are in the R. leguminosarum-pea (Pisum sativum) interaction [36].

In Rhizobium phaseoli CNPAF 512, mutations in genes of the arginine deiminase pathway used for arginine catabolism under low-oxygen conditions, like those occurring in nodules, diminish N2-fixation efficiency [38]. Interestingly, genes for the arginine deiminase pathway were found in B. diazoefficiens strains but not in B. japonicum, and, again, this may also be related to the superior N2-fixation ability of B. diazoefficiens.

Symbiosis islands

Kaneko et al. [2, 4] proposed that the symbiosis islands of Bd USDA 110T and Bj USDA 6T were split in three regions, named locus A, B and C, after their insertion in the respective genomes. When data from transcriptomic studies of USDA 110T are mapped into these three loci, only locus A and B showed a significant number of genes with altered expression under symbiotic conditions. Locus A is the largest region and includes the 410-kb symbiotic region of USDA 110T originally described by Göttfert et al. [39]. Locus B of USDA 110T, 4.9 kb in size, includes only six genes, five of which are hypothetical; it is conserved in USDA 6T[4] and in Bj CPAC 15, but showed variability in Bd CPAC 7 (data now shown). Since locus A contains all nod, nif and fix genes and seems to be a bona fide symbiotic region we decided to focus our analysis on it.

We demarcated the beginning of the symbiosis island in a DNA-recombinase gene (BS08139 in CPAC 15 and BU02649 in CPAC 7), annotated as hypothetical proteins in the genomes of USDA 110T and USDA 6T. The end of the symbiosis island was clearly defined as a tRNA-valine gene. Interestingly, inside the symbiosis islands there were two tRNA-methionine genes in USDA 6T, USDA 110T and CPAC 7, and one in CPAC 15 (Table 2, Figure 4).
Table 2

Properties of the symbiotic islands of B. japonicum and B. diazoefficiens strains

 

B. japonicum

B. diazoefficiens

 

USDA 6T

CPAC 15

USDA 110T

CPAC 7

Size (bp)

694,648

700,213

681,726

688,007

ORFs (total)

646

569

648

567

Hypothetical

448

247

354

294

Mobile elements

64

90

102

77

tRNA-Met

2

1

2

2

https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-15-420/MediaObjects/12864_2014_Article_6142_Fig4_HTML.jpg
Figure 4

Multiple comparison of the symbiosis islands of B. diazoefficiens strains USDA 110 T and CPAC 7, and B. japonicum strains USDA 6 T and CPAC 15. Percentages of island conservation between strains are indicated at right, as well as the average nucleotide identity (ANI) values of the conserved regions. Green triangles indicate the initial tRNA-Val gene while orange triangles indicate other t-RNA genes located inside the islands.

The sizes and numbers of genes identified in the symbiosis islands of the four strains are shown in Table 2. The largest size of the island of CPAC 15 seems to be related especially to a portion of the island that is unique to this strain, as shown in Figure 4. It is noteworthy that, despite a number of functional studies, a large proportion of predicted genes in the symbiosis island remains as hypothetical, representing 43% in CPAC 15 and 52% in CPAC 7. In addition, we point out that the high percentages of mobile genetic elements, of 16% for CPAC 15 and 14% for CPAC 7, in agreement with similar percentages in USDA 110T and USDA 6T[2, 4], are likely to be implicated in genetic plasticity of the symbiosis island. The symbiosis islands of the four strains encompass all nodulation and N2 fixation genes required for symbioses, a probable non-functional hydrogenase operon, cytochromes for energy supply, ABC-transporters, operons for secretion systems, besides individual genes with roles not yet studied, such as a LuxR–type regulator. All these predicted genes are surrounded by several hypothetical ORFs and mobile genetic elements. The information obtained from the symbiosis island can also be related to some of our previous observations. First, reports of high rates of horizontal gene transfer from the symbiosis island of CPAC 15 to indigenous bradyrhizobia [19] might be explained by the high number of mobile genetic elements surrounding nodulation and N2-fixation genes, despite having a seemingly non-functional T4SS. Second, the potentially high plasticity of the symbiosis island may help to explain the extreme variability in N2-fixing properties among variants of the same strain, especially after adaptation to the soil [10, 1416, 18, 19].

The symbiosis islands of the four soybean strains analyzed were not completely syntenic. There was a backbone of highly conserved regions showing ≥ 98% average nucleotide identity that was interrupted by regions unique to each strain (Figure 4). The conserved backbone regions encompassed only 80% to 88% of the islands when B japonicum strains were compared to each other and with the B. diazoefficiens strains. The most similar islands (>98% as conserved backbone) were those of B. diazoefficiens strains.

In the comparison of CPAC 15 and CPAC 7 it is noteworthy that the first shows higher competitiveness and the second higher efficiency of N2 fixation [1, 6, 9, 11, 13, 1820]; Additional file 1: Table S4, Table S5. One first approach to explain these properties is to look for exclusive ORFs of each strain in the symbiosis islands. The Venn diagram shows that the minimum core for the symbiosis islands is composed of 241 predicted ORFs (Additional file 2: Figure S3), encompassing all nodulation, nif and fix genes, type III and IV secretion system genes, cytochromes and several ABC-transporters.

In relation to the exclusive genes of the symbiosis islands of each strain, in CPAC 15, 52% were represented by hypothetical ORFs, whereas 25% were of mobile genetic elements. In CPAC 7, these fractions were even higher, of 65% for the hypothetical and 29% for the mobile genetic element categories. For USDA 6T and USDA 110T, the percentages of exclusive hypothetical ORFs were of 82% and 48%, while for mobile genetic elements were of 16% and 41%, respectively. We may thus conclude that the ORFs occurring exclusively in each strain and encompassing hypothetical and mobile elements ranged from 77% in CPAC 15 to 98% in USDA 6T. Therefore, several studies still have to be performed to understand the particularities hidden in the symbiosis islands of Bradyrhizobium strains.

In relation to the known genes of the symbiosis island, in both CPAC 15 and CPAC 7 we identified a main core of nodulation genes starting with nolZY, followed by one hypothetical protein and two ABC-transporters, nolA, nodD2D1YABCSUIJ, nolNO, another hypothetical ORF, nodZ, fixR, nifA, fixA, and, in CPAC 15, there were two exclusive hypothetical proteins before fixR; high resemblance was also found in all other nif and fix genes. Similarities of all these genes among the four strains were, in general, of 100%. A slightly lower similarity was found among nodD2 of the four strains, and it would be interesting to develop new studies with this gene, recognized as a repressor of nod-gene expression [40]. The operons were thus in close agreement with the first report of the symbiosis island more than a decade ago [39].

Previously, from the genome draft of CPAC 15 [41] we discussed the role of the two-component regulatory system nodV/nodW, first described by Göttfert et al. [42]. In USDA 110T, NodV responds to an environmental stimulus, as the isoflavonoid genistein and phosphorylates the regulator NodW, which in turn may be required to positively regulate the transcription of one or several unknown genes involved in the nodulation of alternative hosts such as mung bean (Vigna radiata), cowpea (Vigna unguiculata), and siratro (Macroptilium atropurpureum), but not soybean. It has been suggested that NodVW recognizes different host flavonoids, or that, combined with NodD, it may increase Nod-signal synthesis [43]. nodV/nodW were present in the symbiosis island of CPAC 15 (BS 08331/30) and were highly similar in CPAC 7 (BU08157/58), USDA 6T (BJ6T_81070/80) and USDA 110T (bll1715/14). In addition, a second copy of nodW adjacent to the first one was present in all four strains (BS0829 in CPAC 15 and BU08159 in CPAC 7); all genes were neighbored by transposases. Interestingly, two other two-component response regulators were found in the chromosomes of all four strains, resembling nodV/nodW. A second copy of nodV/nodW, called nwsA/nwsB, has been previously reported [44], but the location in the genome was not clear. Here we confirm that there were two copies of nodW and one of nodV outside the symbiosis islands (the first one BS03121/22/23 in CPAC 15 and BU07026/27/28 in CPAC 7, and the second BS06847/46/45 in CPAC 15 and BU07376/77/78 in CPAC 7), also found in the chromosomes of USDA 6T and USDA 110T. Therefore, it remains to be determined if these copies also influence host range.

Hydrogen-uptake genes were also present in the symbiosis islands of CPAC 15 and CPAC 7. Hydrogenase has been long study and associated with higher efficiency of N2 fixation, recycling part of the H2 obligatorily evolved in the reduction of atmospheric N2 by the nitrogenase [45, 46]. Several studies with soybean, including field experiments, have shown increases in plant growth and yield (e.g. [47, 48]). However, hydrogenases probably evolved before the N2-fixation process, when Earth was hydrogen enriched, such that many microorganisms—distributed through several taxa—are now able to use molecular hydrogen as an energy source [49].

Recently, some differences between USDA 110T and USDA 6T in their hydrogenase genes have been pointed out [4]. Now, having the genomes of CPAC 15 and CPAC 7, we can confirm the information available [4] and expand it by showing that differences in hydrogenase expression are observed at the species level. All four strains have genes encoding hydrogenase in the symbiosis island, more specifically, hupA(=S)B(=L)CDFH, hypABFDE. However, some of the genes are duplicated and incomplete, or are pseudogenes, such as hupD. The region is flanked by transposases that are also positioned between genes. Therefore, the operons might not be functional and may have gone through an evolutionary process of gene loss. Only in both B. diazoefficiens strains there are hydrogenase genes in other chromosomal loci, represented by almost identical genes to hupNCUVSLCDFGHIJ, hypKABFCDE and hoxXA with 98% to 100% gene similarity. In CPAC 7, these genes were also in a genomic island as in USDA 110T[4]. On the other hand, in B. japonicum, the only remaining gene in this region was hupN (BJ6T_24840 in USDA 6T and BS01414 in CPAC 15). Indeed, a Hup- phenotype has been reported for CPAC 15, while CPAC 7 showed a Hup+ phenotype [14] (Additional file 1: Table S5). In conclusion, the functional hydrogen-uptake system is located in a genomic island distinct from the symbiosis island. The effects of a functional uptake hydrogenase can thus go much further than the reported benefits in the improvement in the efficiency of N2 fixation.

Other genomic islands

We identified 13 and 16 genomic islands (GIs), including the symbiosis island, in Bj CPAC 15 and Bd CPAC 7, respectively, regarded as originating from horizontal gene transfer (HGT). Several HGT events may be implicated in ecological features, such as adaptation, and thus are worthy of study [50]. The GIs of both strains were scattered throughout the genomes (Figure 1), and a list of their main features is shown in Additional file 1: Table S6. GI 9 of CPAC 15 and GI 8 of CPAC 7 corresponded to “symbiosis” locus C as defined by Kaneko et al. [4].

Almost all the GIs identified contained ORFs related to HGT functions, such as transposases, integrases and resolvases. Some included conjugal-transfer systems. Three GIs of each genome encoded restriction-modification systems that may act as mechanisms of resistance to bacteriophages or foreign DNA. Genes that may confer adaptation capabilities such as type I secretion systems and stress tolerance were also identified (Additional file 1: Table S6). It is worth mentioning that, excluding the symbiosis island, 68% and 59% of the ORFs included in GIs were hypothetical in CPAC 15 and CPAC 7, respectively. Seven and thirteen GIs in CPAC 15 and CPAC 7, respectively, were flanked by a tRNA gene. It is known that, in many cases, GIs are derived from integration events that occurred adjacent to tRNA-coding sequences [51].

CPAC 15 and CPAC 7 GIs corresponded to approximately 20% and 17% of their genomes, respectively. In comparison with CPAC 7, CPAC 15 had 3% more of its genome seemingly originating from HGT. GIs ranged from 0.009 to 0.211 Mb, and one GI of CPAC 15 (GI 1, Additional file 1: Table S6) was 0.345 Mb, larger than the others, except for the symbiosis island.

Within the symbiosis island, which is discussed later, it is noteworthy that, in the genome of CPAC 15, a region of about 75 kb contained several ORFs assigned as hypothetical (BS08270–BS08211). This region is enriched with mobile genetic elements and has two copies of putative HipA protein, in addition to another HipAB locus elsewhere. The HipAB system is associated with the stabilization of low copy-numbers plasmids in cell progeny [52], and is potentially related to successful HGT.

Membrane transporters

Transporters—with an emphasis on ABC-transporters—are critical for bacterial adaptation, with roles in cellular physiology, including the uptake of nutrients, exclusion of cellular residues, energy generation and cellular signaling, among others [53]. Acting in several physiological roles, transporters generate functional diversity, and favor interaction with other organisms in the environment. In the case of our strains, the transporters may be of extreme importance, as the Brazilian Cerrados are very poor in nutrients, with low levels of elements critical for the N2-fixation process, such as phosphorus, molybdenum, and calcium [33].

ABC transporters of Bj CPAC 15 and Bd CPAC 7 represented 7.8% and 8.6% of the entire genome, respectively. The strains shared 63.81% of these genes (data not shown). CPAC 7 had more carbohydrate transporters and CPAC 15 more transporters related to uptake of opine, phosphonate and polyamines.

More genes related to phosphorus transport in CPAC 15 could represent another trait contributing to its high saprophytic capacity and competitiveness in the Cerrados soils. Also interesting was the high number of transporters of opines in CPAC 15. Opines are compounds produced by plant tumors or hairy roots induced by pathogenic Agrobacterium species, and represent sources of C and N for the bacteria. However, it has been reported that rhizobia—probably except for the use of nopaline in B. japonicum—were unable to use opine or nopaline as sources of C and N [54]. It has been suggested that utilization of these molecules could improve competiveness against other strains [55], and thus we may hypothesize that this is another trait that contributes to the superior competitiveness of CPAC 15 (Additional file 1: Table S4, Table S5). However, it is noteworthy that only in CPAC 7 we found all subunits of the transporter of manopine (BU00015 to BU00018), absent in both Bj CPAC 15 and Bj USDA 6T, but present in Bd USDA 110T (blr3544 to blr 3547), and thus they could be specific to B. diazoefficiens. In addition, metabolism of manopine may also be related to tolerance of temperature stress [56], a condition typical of the Cerrados region [9, 33].

Both strains had a transporter for D-methionine (operon metQIN), a general L-amino acid transporter (operon aapJQMP) and a branched-chain amino acid transporter (operon livKHMGF), critical for the biological N2-fixation process. Monosaccharides and oligosaccharides are responsible for various processes related to bacterial growth and survival. In both genomes, ribose transporter proteins RbsB, RbsC, RbsA were found, and in addition, monosacharides may be transported by other proteins, such as GruA (present in both genomes), an analogue of the protein AraG (subunit of L-arabinose ABC transporter). In relation to oligosaccharides transporters, the sn-glycerol-3-phosphate transporter (UgpC, UgpA, UgpE) was present in both strains. A multiple sugar transporter (ChvE, GguB, GguA) was also present in both strains and, interestingly, chvE has been suggested as related to the higher efficiency of N2 fixation in variants of CPAC 15 [57]. Within the group of the ABC-2 transporters, both strains had transporters of antibiotics (yadH, yadG), lipo-oligosaccharide (nodJ, nodI), lipoproteins (lolC, lolE, lolD) and cell-division-related transporters (ftsX, ftsE).

In relation to the transport of other limiting nutrients in the Cerrados soils, the transport of sulphur may occur preferentially by means of the operon cysPTWA, closely similar in both strains, and this operon can also transport molybdate, but with a lower affinity, as shown in studies with Eschericha coli[58]. Other transporters of molybdenum, such as modABC and mobABC were found in both strains, as well as genes related to the transport of sulphonate and nitrate, as such as ssuABC.

Surface polysaccharides

Rhizobial surface polysaccharides, such as capsular- (CPS or KPS), exo- (EPS) and lipopolysaccharides (LPS) play important roles in the symbiotic interaction with host plants [59]. Rhizobia also synthethize cell-surface cyclic beta-glucans that are involved in the symbiosis [59]. The roles of surface polysaccharides include signaling, protection against plant pathogens and nodule organogenesis [59, 60].

In many soybean bradyrhizobia, CPS and EPS have the same composition and structure, and differ only in relation to their association with the cell surface [61, 62]. The EPS of B. japonicum USDA 138 consists of a polymer of a pentameric repeated unit composed of (O-methyl or O-acetyl) galactose, two glucoses, O-acetyl galacturonic acid and mannose [62]. EPSs of similar sugar composition are produced by a range of soybean bradyrhizobia, including Bd strains USDA 110T, USDA 122 and CB 1809, and Bj USDA 123 [62, 63].

Genes related to the biosynthesis of EPS were found in CPAC 7 and CPAC 15 genomes, such as those involved in the generation of sugar precursors (exoB, exoN), assembly (exoF, exoM, exoU), and polymerization or transport of repeating units (exoP, exoT, exoQ). A conserved 11.6 kb cluster grouped genes exoQUM, metA and exoPTB. This group of genes was previously described by Becker et al. [64] in USDA 110T, and we noted that it is also present in USDA 6T. In USDA 110T, mutations in exoPT and exoB cause alterations in EPS composition or in the ratio of low molecular weight versus high molecular weight forms of EPS, and delayed nodulation in soybean. The EPS cluster also included genes for a glycosyltransferase and an acetyltransferase. It is worth noting that the acetyltransferase gene is shorter in USDA 110T in comparison to the other strains, although the functional transferase domain seems to be complete in the predicted polypeptide. A second locus related to EPS biosynthesis has been recently identified by Quelas et al. [65] in USDA 110T. Mutations in lspL and ugdH of this cluster resulted in the production of EPS lacking galacturonic acid and diminished competitiveness for soybean nodulation [65]. This 14.8-kb gene cluster is conserved in CPAC 7, CPAC 15 and USDA 6T. Homologues to genes involved in regulation of EPS biosynthesis, such as exoR, mucR and exoS, were found in all four genomes, but a homologue of the exoZ gene required for acetylation of EPS in S. meliloti was present only in CPAC 15.

Species of Bradyrhizobium are able to produce cyclic glucans with β(1 → 6) and β(1 → 3) glycosidic linkages, containing 10 to 13 glucose residues [66]. Although cyclic β-glucans are mostly studied as secreted polysaccharides, they are predominantly located in the periplasm. These polysaccharides have been implicated in the establishment of the N2-fixing symbiosis in terms of suppression of plant-defense responses by their binding to receptor sites on the host-plant plasma membrane [67]. Other findings, such as the importance of cyclic β-glucans for hypo-osmotic adaptation and motility [68], and their high levels in bacteroids [69] and positive effects on isoflavonoid production in soybean [70], support their function in adaptation of rhizobia to the symbiotic lifestyle. In the four genomes analyzed, genes related to the biosynthesis of Family III glucans [71] had the same organization: ndvBDC. NdvB and NdvC are transferases that synthesize the cyclic glucan and a mutation in any tends to impair the symbiosis with soybean [72]. NdvD is a conserved hypothetical protein, potentially related to the transport of glucan to the periplasm [73]. A locus for Family I glucans biosynthesis composed of mdoGH genes was present only in the genomes of CPAC 15 and USDA 110T. Production of type I glucans has not been reported to date in bradyrhizobia, suggesting that expression of mdoGH genes require special conditions.

Molecules of LPS are major structural components of the outer membrane of Gram-negative bacteria and are typically composed of lipid A, core oligosaccharide and the O-antigen polysaccharide. Implication of LPS in several steps of the symbiotic interaction has been extensively reported [59]. Despite differences in the structure and composition, rhizobial lipid A typically contains secondary or even tertiary acyl chains. Recently, it was demonstrated that Bradyrhizobium strains have (ω-1)-hydroxylated very long chain fatty acids (VLCFAs), accompanied by 3-hydroxydodecanoic and 3-hydroxytetradecanoic acids; this pattern is unique to the slow-growing rhizobia [74]. Rhizobial VLCFAs are worthy of attention because of their importance to the intracellular lifestyle [75]. According to Choma & Komaniecka [4], genes related to the biosynthesis of VLCFAs in USDA 110T probably correspond to bll3807–bll3811. This gene cluster was found to have the same organization in the genomes of CPAC 7, CPAC 15 and USDA 6T. Other lpx genes involved in lipid-A biosynthesis were also conserved between the analyzed strains and were located at equivalent positions in their chromosomes.

A gene cluster (galE-lpcC-rfaF-rfaD) involved in LPS’s core oligosaccharide biosynthesis has been described in B. japonicum 61A101C [76]. This cluster is conserved in CPAC 15 and USDA 6T. The O-antigen is probably the most variable portion of the LPS structure [77]. Ferreira & Hungria [8] showed that CPAC 7 and CPAC 15 have different LPS profiles in SDS-PAGE analysis, and here we show how different these profiles are (Additional file 2: Figure S4). We found strain-specific O-antigen biosynthesis gene clusters. In CPAC 7, a 30.8 kb cluster included genes for the biosynthesis of GDP-L-fucose and GDP-L-rhamnose, and wzm-wzt genes. Strain CPAC 15 had a 47.9 kb gene cluster located at the same relative chromosomal position but with distinct genes including a putative wzx gene. The Wzx flippase is one of the mechanisms for translocation of the O-antigen to the periplasmic space; the other mechanism requires the Wzm/Wzt ABC-type transporter. All analyzed strains belong to different serogroups (unpublished data) which is consistent with the presence of strain-specific O-antigen gene clusters in each strain genome.

Secretion systems

Type I secretion system (T1SS)

Eight genes in Bj CPAC 15, Bj USDA 6T and Bd USDA 110T, and seven in Bd CPAC 7 encoded products with glycine- and aspartate-rich repeats that are characteristic of proteins secreted by the T1SS [78]. The B. japonicum strains shared seven orthologous T1SS-exoprotein genes, whereas the B. diazoefficiens strains shared five. Two of these shared genes were present in all four genomes. Most T1SS exoproteins characterized to date belong to animal pathogens and act as toxins, proteases, lipases or adhesins which promote virulence [78]. Interestingly, the two B. japonicum strains and USDA 110T each had one T1SS gene encoding a putative Zn-dependent metalloprotease, though the products are divergent (42-64% identity). In Rhizobium, T1SS exoproteins have been characterized as bacteriocins promoting competitiveness for nodulation [79], or as factors like NodO, playing a yet-unknown role in nodulation of some hosts [78]. None of the T1SS exoproteins of our strains showed similarity with known bacteriocins or NodO; nevertheless, some of these proteins bear domains related to adhesins and may be involved in competitiveness for root colonization. In S. meliloti 1021, the expression of a secreted peroxidase was highly induced upon oxidative stress [80]. Each B. japonicum strain genome encoded a putative T1SS-secreted peroxidase that might improve fitness under stressful conditions.

Type II secretion system (T2SS)

Genes of the T2SS molecular machinery were found only in Bd USDA 110T and Bj USDA 6T; however, gspC, encoding a core component of the T2SS inner-membrane platform [81] was absent in both, indicating that strains are incapable of type II secretion. In Bj CPAC 15, a gspO gene remnant was located in the region corresponding to the gsp locus, and in Bd CPAC 7 this locus and adjacent genes seem to have been replaced by a polysaccharide-biosynthesis gene cluster. This deletion/replacement of the gsp locus in CPAC 15 and CPAC 7 may have been favored if the ancestor of all of the strains had a nonfunctional T2SS.

Type III secretion system (T3SS)

Several rhizobia, including Bd USDA 110T and Bj USDA 6T possess genes required for biosynthesis of T3SS in their symbiotic compartments [2, 4]. The T3SS acts as a molecular syringe to inject effector proteins into legume host cells, affecting the plant defense systems to promote colonization [82]. USDA 110T mutants in the T3SS show delayed nodulation and diminished N2- fixation capacity with G. max, but are all able to normally nodulate and fix N2 with V. unguiculata[83]. The activities of T3SS effectors may have detrimental effects in some legume hosts [82] even at the cultivar level. The latter effect is exemplified by the ability of T3SS mutants of Bd USDA 122 to nodulate the soybean cultivar Hardee (Rj 2 ), in which nodulation is restricted to certain bradyrhizobial strains such as USDA 122 [84]. The genomes of CPAC 15, CPAC 7, as well as USDA 6T and USDA 110T share a highly conserved T3SS cluster within their symbiosis island. Small differences between these clusters were observed and always involved the presence or absence of genes related to transposases and other mobile genetic elements. These differences are not expected to produce distinct phenotypes related to T3SS activity.

Type IV secretion system (T4SS)

Both Bj CPAC 15 and Bd CPAC 7 genomes encoded genes for T4SS. A T4SS present in CPAC 15 included all mating-pair formation (Mpf) genes required for assembly of the T4SS pilus, a set of DNA-transfer and -replication (Dtr) genes required for preparing the DNA to be transferred, and a gene for a coupling protein that delivers the DNA intermediate to the Mpf pilus [85]. This T4SS of CPAC 15 was part of GI 6 present only in this strain and may participate in conjugative transfer of this island as has been described for other T4SS [86].

Another T4SS was found encoded in the symbiosis island. However, this system seemed to be non-functional as it included only a partial set of Mpf genes, those homologous to trbLFGI, and a partial gene homologous to those encoding TraG-coupling proteins. These T4SS genes were found to be conserved between Bj CPAC 15 and USDA 6T and Bd USDA 110T. The trbLFGI genes were not present in CPAC 7, reinforcing the possibility that this is a relic T4SS system. The symbiosis islands of all strains lacked genes for Dtr components, suggesting that the original T4SS was devoted to export effector proteins, as has been reported for the T4SS located in the Mesorhizobium loti R7A symbiosis island; however, we could not identify proteins similar to known T4SS effectors.

A third T4SS was identified in CPAC 15 within the GI 9. Genes for assembly of a complete Mpf pilus, for DNA processing and for a coupling protein were identified, suggesting that this system is required for conjugative transfer of this locus. GI 9 was conserved in CPAC 7 (as GI 8), USDA 110T and USDA 6T. However, it showed some variability between strains; for example, the T4SS genes were present in USDA 110T but were absent in CPAC 7 and USDA 6T.

Type V secretion system (T5SS)

Genes encoding the autotransporter T5SS subclass were found in all genomes (six in Bj 5079, four in Bj USDA 6T, and five in both B. diazoefficiens strains). Two and three of these genes in the B. japonicum and B. diazoefficiens strains, respectively, code for products possessing the pectin-lyase-type beta-helix fold at their N halves, suggesting that they may act as adhesins binding to alpha-galactose-containing polymers [87]. Putative hemagglutinin-like adhesins secreted by the two partner system T5SS subclasses were also found, two in the B. japonicum strains and one in the B. diazoefficiens strains. Three genes for probable adhesins resembling YadA of Yersinia entercolitica and secreted by the trimeric autotransporter T5SS subclass were identified in the B. diazoefficiens strains. One YadA-like gene was found in USDA 6T, but none in CPAC 15. Proteins secreted by the T5SS are generally regarded as virulence factors promoting host colonization by pathogenic bacteria [88], but their role, if any, in symbiotic relationships, remains to be established.

Type VI secretion system (T6SS)

A 14-gene locus encoding a T6SS was found in each strain genome. A between-strain comparison of these loci showed that they were syntenic and their gene products had identities larger than 91%, indicating that they are orthologous. ImpG/SciC/VasA have been suggested to act as a structural components required for T6SS functioning [89, 90]. In Bd USDA 110T, the gene coding for the ImpG/SciC/VasA component has undergone a mutation (generating an internal stop codon that splits it in two putative genes, bll3589 and bll3590), implying this secretion system is not active in USDA 110T. The T6SS molecular machinery is present in free-living pathogenic and symbiotic bacteria, and may be involved in prokaryote-eukaryote and prokaryote-prokaryote interactions [91]. The function of these secretion systems is not well understood in the context of rhizobia-legume interactions. In R. leguminosarum, the T6SS is required for secretion of some proteins that seem to interfere with symbiosis with pea. The wild-type strain is unable to fix N2, but a mutation in one of the T6SS genes allows N2 fixation in the host [92].

Other possible genetic determinants of ecological traits related to colonization of soil and root

Type IV pili (T4P)

We found two distinct T4P gene clusters shared by the B. diazoefficiens and B. japonicum strain genomes, and two additional clusters present only in the B. diazoefficiens genomes. All of these T4Ps were related to the Cpa and Tad pili found in Caulobacter crescentus and Aggregatibacter actinomycetemcomitans, respectively, and thus they can be classified into subtype IVb. This T4P subtype includes those promoting host colonization of pathogenic bacteria [93], and as-yet uncharacterized T4P of other rhizobia such as S. meliloti 1021, Sinorhizobium sp. NGR234 and Mesorhizobium huakui MAFF303099 [94] which possess 2, 3 and 3 T4P gene clusters, respectively. In B. diazoefficiens, the two additional clusters may favor host colonization.

Flagella

Bd USDA 110T produces two types of flagellum, one called “thin” composed of 33-kD flagellin subunits and a “thick” one with 65-kD subunits [95]. Flagellin genes flaC1C2C3C4 and flaCICII, located in different gene clusters, are required for biosynthesis of thick and thin flagella, respectively [95]. In general, the gene clusters required for biosynthesis of both types in USDA 110T are conserved in Bd CPAC 7 and also in Bj CPAC 15 and USDA 6T; however, differences were observed in the flagellin-coding genes for the thick flagellum. Given that the flagellum encoded in the thick cluster seems to be important for competitiveness [96], it would be interesting to explore if the difference between Bj and Bd strains in correlate with competitiveness.

Quorum sensing

Quorum sensing (QS) is a cell-cell communication mechanism that allows bacteria to regulate gene expression in response to fluctuations in cell-population density and to coordinate group behaviors [9799]. In general, QS in Gram-negative bacteria is mediated by N-acyl-homoserine lactone (AHL)-signaling molecules, which are synthesized by LuxI protein homologs from S-adenosylmethionine (SAM) and acylated acyl-carrier proteins (acyl-ACPs). The LuxR transcriptional regulators bind to cognate HSL autoinducers to activate gene transcription of QS-target genes. QS has recently been linked to several important symbiosis features, including nodulation efficiency, symbiosome development, EPS production and N2 fixation [98, 99]. Approximately 22% (31 of 142) of the B. japonicum and B. elkanii strains tested produced AHLs that induced moderate to elevated levels of β-galactosidase activity with an Agrobacterium tumefaciens biosensor [100]. Both Bj CPAC 15 and Bd CPAC 7 strains each contain a luxI-type synthase, identified as BS04783 and BU02130, respectively, showing over 95% of identity to BjaI protein (blr1063) of USDA 110T and with USDA 6T (BJ6T_10890). Adjacent to these ORFs, there is in both genomes a luxR-type regulator (BS04784 and BU02129), also with high identity (>94%) to BjaR (blr1062) of USDA 110T and with USDA 6T (BJ6T_10880). Neither gene is in any of the genomic islands. According to Lindemann [101], BjaI and BjaR proteins of USDA 110T are responsible for the synthesis of the QS signal, a branched-chain fatty acyl-HSL, isovaleryl-HSL (IV-HSL). Similar molecules may be synthesized by other strains.

In addition to the genes described above, there are isolated LuxR-type regulators in all four strains, CPAC 15 (BS08299), CPAC 7 (BU08021), USDA 110T (blr1880) and USDA 6T (BJ6T_79490) located in the symbiosis island, and a possible role for this gene in the regulation of the symbioses is an exciting subject for investigation.

Iron uptake

Mechanisms that facilitate iron acquisition may confer to bacteria higher survival capacity and competitiveness [102]; even in environments not iron depleted, as is the case for the Cerrados soils [33]. Siderophores are high-affinity low-molecular-mass chelating compounds used to capture iron. In the Bj CPAC 15 and Bd CPAC 7 genomes, we found genes related to the biosynthesis of siderophores, such as citrate and cathecolate siderophore, synthesized via chorismate, in addition to genes related to iron uptake, storage and regulation.

In B. japonicum 61A152, citric acid can act as a siderophore [103], and Lesueur et al. [104] reported that Bradyrhizobium strains excrete citrate when iron-starved. CPAC 15 and CPAC 7, as well as Bd USDA 110T and Bj USDA 6T had three genes encoding citrate synthases. We also found in CPAC 15 and CPAC 7 genes related to the biosynthesis of cathecolate siderophores, previously isolated and characterized in Bradyrhizobium[105]. Nevertheless, we did not find all genes known to be in the biosynthetic pathway, e.g. isochorismate synthase and 2,3-dihydro-2,3-dihydroxybenzoate dehydrogenase were missing, but we found a chorismate synthase and isomerases that might catalyze the conversion of chorismate to isochorismate; in addition, there were several genes encoding dehydrogenase activity, which could replace the 2,3-dihydro-2,3-dihydroxybenzoate dehydrogenase [106]. We also found isochorismatases that were exclusive to CPAC 15 (BS06609, BS07416) and CPAC 7 (BU08273, BU08403), each being absent from the other three strains.

CPAC 7 possessed a gene encoding a receptor for ferric enterobactin (BU07037)—a cathecolate siderophore—which was not identified in CPAC 15. Most soil microorganisms form hydroxamate siderophore receptors [107], and the majority of receptors identified in both soybean strains were also for hydroxamate siderophores. On the other hand, CPAC 15 possessed more exclusive receptors for hydroxamate siderophores than did CPAC 7, five against one, respectively. Some rhizobial strains can utilize siderophores from other organisms by producing appropriate receptors and transporters [108], e.g., USDA 110T and 61A152 are able to utilize exogenous siderophores, such as ferrichrome, rhodotorulate and pyoverdin-types [109]. These differences may permit CPAC 15 to produce and to utilize more siderophores, and to take up more iron for growth and survival, representing another feature that might help to explain its superior competitiveness compared to CPAC 7.

Phytohormone production

Several plant hormones have been reported to positively or negatively regulate nodulation. Whereas auxins and cytokinins (CKs) are positive regulators of nodule organogenesis and development, ethylene generally blocks both infection and nodule-primordia initiation [110, 111]. In addition, we must keep in mind that rhizobia also produce phytohormones during the symbiosis that can change the phytohormone balance in the host [112].

Cytokinin is a key signaling molecule in symbiotic interactions between leguminous plants and rhizobia. Several studies have shown that CK signaling is necessary to induce cortical cell division and nodule organogenesis [110, 111, 113]. CK biosynthesis may take place by the action of the tRNA-dimethylallyltransferase (MiaA; EC 2.5.1.75), product of miaA gene, by a cys-zeatin biosynthesis route that utilizes dimethylallyl diphosphate (DMAPP) and tRNA as substrates. Alternatively, CK synthesis may involve the ipt gene product, an adenylate isopentenyltransferase (IPT; EC 2.5.1.27) that utilizes DMAPP and AMP, ADP or ATP as substrates in the trans-cytokinin pathway [114]. With A. tumefaciens, the gall formation is dependent on IPT activity and other CK biosynthetic proteins, which are responsible for overproduction of this hormone [114]. Like other symbionts, Bradyrhizobium produce CKs both in vitro and during the symbiosis [115117]. In Bj CPAC 15 and Bd CPAC 7, CK production and effects on nodulation have not been studied, but both strains contain miaA (BS01001 and BU03527), showing 99% identity with Bd USDA 110T, 90% with Bj USDA 6T and 76% with Bradyrhizobium sp strain ORS285. In all these strains, miaA is located adjacent to serB and ilvI genes, predicted to encode a phosphoserine phosphatase (EC 3.1.3.3) and an acetolactate synthase large subunit (EC 2.1.1.6), respectively. Therefore, CK production in CPAC 15 and CPAC 7 seems to occur via the cys-zeatin route, since orthologs of the ipt gene were not found. The importance of miaA has been confirmed recently in a study with a mutant of the photosynthetic Bradyrhizobium sp strain ORS285, revealing that tRNA degradation is the major route of CK synthesis and that the deficiency of CK production in the mutant accounts for the delay in nodulation and N2 fixation, as well as with smaller nodules in Aeschynomene[117].

Another phytohormone that positively regulates nodule formation is auxin or indole-3-acetic acid (IAA) [110]. Auxin is involved at different stages of nodule development, including the early stages of cell initiation, stimulation of early cell division, differentiation of nodules and in the systemic regulation of nodule numbers [112]. Because many rhizobia are capable of producing IAA via different pathways, it is also assumed that bacterially produced auxin can alter the balance inside the plant. In addition, rhizobia can indirectly influence auxin homeostasis by interfering with plant-auxin transport [112]. Bradyrhizobium strains synthesize auxin [14, 116, 118120], and effects of IAA produced by Bradyrhizobium on the symbiosis have been reported. Nodules of soybean inoculated with a B. japonicum IAA-overproducing mutant have bacteroids with higher IAA content [121, 122], and positive and negative effects of this increased level on nodulation have been reported [121, 123], indicating that IAA level in nodules are tightly regulated, but that the role of IAA in the symbiosis is not clear yet. In vitro CPAC 15 produces more IAA than CPAC 7 (Additional file 1: Table S5).

In bacteria, at least five pathways have been described for the synthesis of IAA [124126] starting from tryptophan. The indole-3-acetamide (IAM) route is the most studied, particularly in plant pathogens such as A. tumefaciens and Pseudomonas savastanoi[124, 126, 127]. In this pathway, Trp is first converted into IAM by a Trp 2-monooxygenase (TMO), and then to IAA by IAM-hydrolase (IaaH). In A. tumefaciens and Pseudomonas savastanoi, the genes encoding TMO (tms-1 or iaaM) and IaaH proteins (tms-2 or iaaH) constitute an operon and are associated with pathogen virulence and gall formation in plants [124, 126]. The existence of the IAM route in Bradyrhizobium sp. strains has been considered based on the detection of IAM in vitro and also by the determination of IAM-hydrolase activity and sequencing of the bam gene, an ortholog of the iaaH gene [118120]. Three ORFs in the CPAC 15 (BS05569, BS01233, BS02084) and two in the CPAC 7 (BU07366, BU07366) genomes showed high similarity to bam and iaaH genes. BS05569 and BU07366 are clustered with ORFs involved in cysteine and methionine metabolism (metA and metY) and with ORFs encoding a polysaccharide deacetylase and a 1-aminocyclopropane-1-carboxylate deaminase (acdS) related to ethylene biosynthesis. BS01233 and BU07366 are clustered with a glutathione peroxidase, whereas BS02084, exclusive to CPAC 15, is positioned together with an indolepyruvate ferredoxin oxidoreductase (IOR, EC 1.2.7.8). The other enzyme of the IAM route, the Trp 2-monooxygenase (TMO), was not found in the vicinity of any ORF encoding IAM-hydrolase protein and blast searches with tms-1 and iaaH orthologs against CPAC 15 and CPAC 7 genomes were not successful. Previously, attempts to amplify the TMO gene and to determine TMO activity of Bradyrhizobium strains failed [120]. In addition, the genome analysis of CPAC 15 and CPAC 7 supports the proposal that, in Bradyrhizobium, the synthesis of IAM is independent of TMO activity and related to action of nitrile hydratase (Nhase; EC 4.2.1.84), which converts indol-3-acetonitrile (IAN) into IAM [120].

The genes for cobalt-containing nitrile hydratase (Nhase) α and β subunits (nthA and nthB) were located in the genomes of CPAC 15 (BS03457 and BS03458) and CPAC 7 (BU04142 and BU04143) downstream of an ORF encoding a Nhase activator protein (BS 03459 or BU04144, respectively), which is required for full activity of their respective NHase. Another enzyme involved in IAA biosynthesis via indol-3-acetonitrile (IAN) is nitrilase (Nit; EC 3.5.5.1), which converts IAN directly into IAA. CPAC 7 contained a nitrilase gene (BU00160) with high similarity to that of USDA 110T (blr3397), reported to be an aliphatic nitrilase with high affinity for hydrocinnamonitrile and with low activity towards indol-3-acetonitrile [128]. In CPAC 15, three genes showing similarity to nitrilase (BS06879, BS00934, BS02760) were found; interestingly, BS02760 is adjacent to an ORF predicted to encode an acetaldoxime dehydratase, which catalyzes IAN production via intermediate indole-3-acetaldoxime (IAOx).

It is noteworthy that Boddey & Hungria [14] reported that the capacity for synthesis of IAA in vitro by CPAC 15 is greater (31.5 μM) in comparison to CPAC 7, USDA 110T and USDA 123 (13.12, 6.85 and 4.88 μM, respectively), what could be explained by the higher number of ORFs detected in CPAC 15. This feature should be explored more thoroughly, as the use of Bradyrhizobium as plant-growth promoting rhizobacteria (PGPR) of non-legumes has been suggested [129]. The delivery of phytohormones straight to plants is suggested by the findings that Bradyrhizobium is found as an endophyte of non-legumes, such as wild rice (Oryza breviligulata) in Senegal and Guinea [130] and sugarcane (Saccharum spp.) in Brazil [131]. Therefore, the study of elite Bradyrhizobium strains, such as CPAC 15, may have important implications beyond the symbiosis with soybean. Contributions of Bradyrhizobium to the N-nutrition of non-legumes is also worth investigating given recent findings of bradyrhizobial nifH transcripts inside sugarcane [132].

A transcriptomic study with USDA 110T incubated with 1 mM IAA resulted in 1,323 genes differentially expressed, most related to responses to heat, cold, oxidative, osmotic and desiccation stresses and in EPS biosynthesis [133]. Therefore, again, IAA may play several roles.

Secondary metabolism

Although secondary metabolites are not essential for growth and reproduction of organisms, they represent natural molecules of adaptation [134], playing several roles, including development of the symbioses, responses to predators, as well as antibiotics and effectors in ecological competitiveness [134136]. Several secondary metabolic pathways were identified in both Bj CPAC 15 and Bd CPAC 7, with the largest numbers being related to the biosynthesis of butirosin and neomycin (aminoglycoside family). Although several ORFs related to the isoquinoline alkaloid pathway were found, the majority was repeated or encoded the same protein. These genes are responsible for subpathways for galanthamine, colchicine and kreysigine synthesis. In both strains, within the pathway for stilbenoid, diarylheptanoid and gingerol biosynthesis, there was a subpathway for the production of 6-gingerol and curcumin, two compounds showing antitumor activity [137]. Other rhizobial molecules have been shown to possess antitumor activity [138].

ORFs related to the hexosyltransferase family (EC 2.4.1.-), the “transferring groups other than amino-acyl groups” family (EC 2.3.1.-) and the methyltransferase family (EC 2.1.1.-) were present in high number in the anthocyanin biosynthesis pathway in both strains. These compounds have antibiotic, antiviral and antifungal activities [139] and may contribute to strain survival in soil.

Both strains have genes related to the conversion of DIBOA-glucoside in TRIBOA-glucoside and DIMBOA-glucoside, and most ORFs for benzoxazinoid pathway show that this subpathway may be functional. Usually benzoxazinones are released by plants and have antifungal and antibacterial activities [140] and being metabolized by strains CPAC 15 and CPAC 7 may represent an advantage in terms of survival.

In relation to the biosynthesis of puromycin, novobiocin and streptomycin, we identified genes in both genomes, but the biosynthesis pathways were incomplete in both strains. In the puromycin pathway (pur cluster), we found pur10, pur4, and pur5 in both strains, and a putative pur7 in CPAC 15, suggesting that it is capable of metabolizing ATP into 3'-amino-3'-deoxy-AMP, which can then be used in other secondary biosynthetic processes. Similarly for streptomycin-biosynthesis, algC (phosphomannomutase/phosphoglucomutase) and glk (glucokinase) were present, capable of converting D-glucose-1-phosphate and D-glucose into D-glucose-6-phosphate, respectively, which can then be used for other purposes including biosynthesis of butirosin and neomycin.

The high number of ORFs related to secondary metabolites in both strains might be related to the ecology of these bacteria, in particular to the adaptation to the Cerrados soils. In addition, in almost all categories of secondary metabolites CPAC 15 had higher numbers of ORFs than CPAC 7 (Additional file 2: Figure S5), possibly related to the superior saprophytic capacity and competitiveness of the former.

In a second approach, both genomes were analyzed to find genes or gene clusters that encode proteins involved in secondary metabolite biosynthesis using the antiSMASH software [141] (Additional file 1: Table S7). Both CPAC 15 and CPAC 7 genomes encoded a two-gene cluster for bacteriocin synthesis. Bacteriocins are important in competitiveness of rhizobial strains. It has been reported that bacteriocin-producing strains have higher capacity for competitiveness and occupancy of nodules [142]. In addition, both strains share gene clusters responsible for the production of terpenes, chains of isoprene units, and peptide secondary metabolites. Only CPAC 15 encoded a type I polyketide synthase, a multi-domain enzyme that may be responsible for the synthesis of a yet unknown compound.

Conclusion

The soybean-Bradyrhizobium symbiosis is considered as one of the most efficient in fixing N2 and probably the highest in global economic importance. Gaining a better understanding of the microsymbiont is thus important as it can open new frontiers for the improvement of the symbiosis not only with soybean, but also with other legumes. Here we described the genomes of two inoculant elite strains, Bj CPAC 15 and Bd CPAC 7—applied yearly to millions of hectares in Brazil, other countries of South America and, more recently, in African countries—as well as the genome comparison with B. japonicum and B. diazoefficiens type strains. Most of the genomes were syntenic. However, genome rearrangements were observed both between and within species, conferring high genomic plasticity, which may generate environmental adaptability and contribute to differences in symbiotic performance. Our still-poor knowledge of these important bacteria was emphasized with the observation that it was not possible to attribute any known function to about 50% of the predicted genes, including genes of the symbiosis island. The high proportion of ORFs related to the metabolism of amino acids and carbohydrates of CPAC 7 and CPAC 15 might be related to the broad adaptability of both strains to tropical soils with low N and C contents. CPAC 15, outstanding in saprophytic capacity and competitiveness had more genes in all categories of secondary metabolites, transporters of nutrients, iron-acquisition and IAA metabolism. In CPAC 7, characterized by higher efficiency of N2 fixation, more ORFs were related to the metabolism of amino acids, particularly in the metabolism of leucine, isoleucine and arginine; and a functional hydrogenase operon. Other differences between the species were detected, including size of the genomes, number of ribosomal operons, genes related to hydrogenase activity, LPS and type IV pili. Interesting genes were predicted in all strains, including those related to type I, III, IV, V and VI secretion systems, EPS, quorum sensing and genes related to the synthesis of phytohormones, among others. Although the backbone of the symbiotic island was conserved among the strains, the conserved region was interrupted by regions unique to each strain and many hypothetical ORFs and mobile elements that may contribute to differences in symbiotic performance.

Methods

Background information about CPAC 15 and CPAC 7

B. japonicum (Bj) CPAC 15 (=SEMIA 5079, = DF 24, =CNPSo 7) shows high saprophytic capacity (i.e. capable of surviving for long periods in soils even under nutrient-limiting and environmental stressing conditions) and high competitiveness in comparison with other soybean Bradyrhizobium strains, whereas B. diazoefficiens (Bd) CPAC 7 (=SEMIA 5080, =CNPSo 6) has higher efficiency of N2 fixation but lower competitiveness and saprophytic capacity than CPAC 15 [1, 6, 9, 11, 13, 1820]. Both strains are well adapted to tropical environmental stressful conditions, including high temperature and low soil moisture and fertility, with an emphasis on soil acidity and low phosphorus content. Some morphological, physiological, symbiotic and genetic differences previously reported between the two strains are summarized in Additional file 1: Table S5.

Bacterial strains and growth conditions

Bj strain CPAC 15 and Bd strain CPAC 7 were obtained from the “Diazotrophic and Plant Growth Promoting Bacteria Culture Collection of Embrapa Soja” (WFCC Collection # 1213, WDCC Collection # 1054) located in Londrina, Paraná, Brazil. CPAC refers to Embrapa Cerrados, Planaltina, Federal District, Brazil. Bacterial growth conditions and DNA extraction were performed as described before [41].

Sequencing, assembly and gap closure

Library generation and genome sequencing of both strains were performed at “Darcy Fontoura de Almeida” Computational Genomics Unity (UGCDFA) of the National Laboratory of Scientific Computation (LNCC) (Petrópolis, Rio de Janeiro, Brazil). The genomes were sequenced using a whole-genome shotgun strategy, with a combination of Roche 454 GS-FLX shotgun and 3 kb-insert paired-end libraries. Libraries were prepared following GS FLX Titanium series protocols. 454 sequence reads were assembled using both Newbler 2.6 (454 Life Sciences, Roche Diagnostics Corporation, Branford, CT) and Celera (WGS, version 7.0) assemblers. The two assemblies were aligned to each other with Cross Match (Phred/Phrap/Consed package) since both yielded results that can be complementary and effective to close gaps. Sanger sequences previously obtained for strain CPAC 15 at Embrapa Soja (Brazil) [41] were also joined to the pyrosequence data. For closing gaps, a primer walking strategy was used, PCR products were Sanger sequenced at Embrapa Soja. All the consensus sequences of each contig and the sequence that closed each gap were aligned and joined using Consed (version 20.0). Information for scaffolding and gap closure was also obtained by mapping each assembly against the reference genomes of B. diazoefficiens USDA 110T (BA000040.2) and B. japonicum USDA 6T (AP012206.1).

Annotation

Annotation and analysis of the sequences were carried out using the System for Automated Bacterial Integrated Annotation (SABIA) [143]. An automatic functional annotation was performed using the KEGG database according to the following criteria: i) ORFs with a BLASTP hit on KEGG database with a minimum 50% similarity, 60% query coverage and 80% subject coverage were assigned “valid”. The first three hits were analyzed and the product was imported from KEGG Orthology (KO) if there was one associated with the hit, or from KEGG GENES definition if no KO was associated with the first three hits; ii) ORFs that had (1) no BlastP hits found on the NCBI-nr, KEGG, UniProtKB/Swiss-Prot, TCDB and Interpro databases, or (2) the first three BlastP hits product on KEGG containing the keyword hypothetical were assigned “hypothetical”. A manual annotation by comparison with the UniProt/Swiss-Prot, KEGG, NCBI-nr and InterPro databases was performed for ORFs that did not fit the above criteria. In addition, annotation was also performed using the Rapid Annotation using Subsystem Technology (RAST) server [144]. Data for CPAC 15 and CPAC 7 were submitted to the GenBank database and were assigned Bioprojects numbers PRJNA20963 and PRJNA47329, respectively. Accession numbers for the genomes at the GenBank are CPAC 7 (ADOU00000000) and CPAC 15 (CP007569).

Comparative genomics and other bioinformatics analyses

Genome alignments of CPAC 15, CPAC 7, USDA 110T and USDA 6T were performed with Mummer [145] and Mauve [146]. The genomes were also analyzed by the Bidirectional Best-Hits (BBH) clustering method [147] that compares the genome of each strain against each of the other genomes, using the BLAST program [148] to identify pairs of corresponding genes (clusters) and to recognize the best “hit” in other genomes. The parameters applied were 60% coverage, 60% of similarity and an E-value of <10-5. Venn diagrams were built with GeneVenn and VENNY [149, 150] to compare the list of shared and unique genes.

For transporters, a more detailed analyses was performed with the Transporter Classification Database—TCDB database [151]. Another database used was PFAM [152] to find domains related to the ABC 1, 2 and 3 transporters. Once the information was found, a search in the in silico proteome of CPAC 15 and CPAC 7 was performed with HMMER [153] to search for similarities in the protein sequences using an E-value <10-4.

Genomic islands were predicted using IslandViewer [154]. Search for T4SS effector proteins were performed using the SecReT4 web server [155]. Other software used for specific analyses were cited in the text.

Availability of supporting data

Supporting data are included as additional files. In addition, detailed information about the genomes are available at the homepages:

http://ligeirinha.lncc.br/bj5079-final-bin/general_info.cgi for B. japonicum CPAC 15; http://ligeirinha.lncc.br/bj5080-final-bin/general_info.cgi for B. diazoefficiens CPAC 7.

Notes

Declarations

Acknowledgments

The study was partially supported by CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazil), Project Repensa (562008/2010-1). Authors acknowledge Ligia Maria O. Chueire for help in several steps of this work and to Dr. Allan R.J. Eaglesham for English review and suggestions on the manuscript. A.F.S. acknowledges an MSc fellowship from CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil), J.S.S.B is a postdoc fellow from CNPq and A.S.N. is postdoc fellow from Fundação Araucária/CAPES. M.H., F.G.B. and A.T.R.V. are also research fellows from CNPq. Embrapa Soja and Centro de Ciencias Genómicas had a bilateral collaboration project by CNPq/CONACYT (490048/2009-9).

Authors’ Affiliations

(1)
Department Biochemistry and Biotechnology, Universidade Estadual de Londrina (UEL)
(2)
Embrapa Soja, C.P. 231
(3)
Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México
(4)
Laboratório Nacional de Computação Científica
(5)
UEL, Depto. General Biology
(6)
Department Structural, Molecular and Genetic Biology, Universidade Estadual de Ponta Grossa (UEPG)

References

  1. Hungria M, Mendes IC, de Bruijn F: Nitrogen fixation with soybean: the perfect symbiosis?. Biological Nitrogen Fixation. 2014, New Jersey: Wiley Publisher, Hoboken, in pressGoogle Scholar
  2. Kaneko T, Nakamura Y, Sato S, Minamisawa K, Uchiumi T, Sasamoto S, Watanabe A, Idesawa K, Iriguchi M, Kawashima K, Kohara M, Matsumoto M, Shimpo S, Tsuruoka H, Wada T, Yamada M, Tabata S: Complete genomic sequence of nitrogen-fixing symbiotic bacterium Bradyrhizobium japonicum USDA110. DNA Res. 2002, 9 (6): 189-197. 10.1093/dnares/9.6.189.PubMedGoogle Scholar
  3. Delamuta JR, Ribeiro RA, Ormeño-Orrillo E, Melo IS, Martínez-Romero E: Polyphasic evidence supporting the reclassification of Bradyrhizobium japonicum group Ia strains as Bradyrhizobium diazoefficiens sp. nov. Int J Syst Evol Microbiol. 2013, 63 (9): 3342-3351.PubMedGoogle Scholar
  4. Kaneko T, Maita H, Hirakawa H, Uchiike N, Minamisawa K, Watanabe A, Sato S: Complete genome sequence of the soybean symbiont Bradyrhizobium japonicum strain USDA6T. Genes. 2011, 2 (4): 763-787. 10.3390/genes2040763.PubMed CentralPubMedGoogle Scholar
  5. Tian CF, Zhou YJ, Zhang YM, Li QQ, Zhang YZ, Li DF, Wang S, Wang J, Gilbert LB, Li YR, Chen WX: Comparative genomics of rhizobia nodulating soybean suggests extensive recruitment of lineage-specific genes in adaptations. Proc Natl Acad Sci U S A. 2012, 109 (22): 8629-8634. 10.1073/pnas.1120436109.PubMed CentralPubMedGoogle Scholar
  6. Hungria M, Campo RJ, Mendes IC, Graham PH: Contribution of biological nitrogen fixation to the N nutrition of grain crops in the tropics: the success of soybean (Glycine max L. Merr.) in South America. Nitrogen Nutrition and Sustainable Plant Productivity. Edited by: Singh RP, Shankar N, Jaiwal PK. 2006, Houston: Studium Press, LLC, 43-93.Google Scholar
  7. Hungria M, Franchini JC, Campo RJ, Crispino CC, Moraes JZ, Sibaldelli RNR, Mendes IC, Arihara J: Nitrogen nutrition of soybean in Brazil: contributions of biological N-2 fixation and N fertilizer to grain yield. Can J Plant Sci. 2006, 86: 927-939. 10.4141/P05-098.Google Scholar
  8. Ferreira MC, Hungria M: Recovery of soybean inoculants strains from uncropped soils in Brazil. Field Crops Res. 2002, 79 (2–3): 139-152.Google Scholar
  9. Hungria M, Vargas MAT: Environmental factors affecting N2 fixation in grain legumes in the tropics, with an emphasis on Brazil. Field Crops Res. 2000, 65 (2–3): 151-164.Google Scholar
  10. Torres AR, Kaschuk G, Saridakis GP, Hungria M: Genetic variability in Bradyrhizobium japonicum strains nodulating soybean [Glycine max (L.) Merrill]. World J Microbiol Biotechnol. 2012, 28: 1831-1835. 10.1007/s11274-011-0964-3.PubMedGoogle Scholar
  11. Peres JAR, Mendes LC, Suhet AR, Vargas MAT: Eficiência e competitividade de estirpes de rizóbio para soja em solos de cerrado. Rev Bras Ciênc Solo. 1993, 17: 357-363.Google Scholar
  12. Hungria M, Nishi CYM, Cohn J, Stacey G: Comparison between parental and variant soybean Bradyrhizobium strains with regard to the production of lipo-chitin nodulation signals, early stages of root infection, nodule occupancy, and N2 fixation rates. Plant Soil. 1996, 186: 331-341. 10.1007/BF02415528.Google Scholar
  13. Nishi CYM, Boddey LH, Vargas MAT, Hungria M: Morphological, physiological and genetic characterization of two new Bradyrhizobium strains recently recommended as Brazilian commercial inoculants for soybean. Symbiosis. 1996, 20 (2): 147-162.Google Scholar
  14. Boddey LH, Hungria M: Phenotypic grouping of Brazilian Bradyrhizobium strains which nodulate soybean. Biol Fertil Soils. 1997, 25: 407-415. 10.1007/s003740050333.Google Scholar
  15. Santos MA, Vargas MAT, Hungria M: Characterization of soybean Bradyrhizobium strains adapted to the Brazilian savannas. FEMS Microbiol Ecol. 1999, 30 (3): 261-272. 10.1111/j.1574-6941.1999.tb00654.x.PubMedGoogle Scholar
  16. Barcellos FG, Menna P, da Silva Batista JS, Hungria M: Evidence of horizontal transfer of symbiotic genes from a Bradyrhizobium japonicum inoculant strain to indigenous diazotrophs Sinorhizobium (Ensifer) fredii and Bradyrhizobium elkanii in a Brazilian savannah soil. Appl Environ Microbiol. 2007, 73 (8): 2635-2643. 10.1128/AEM.01823-06.PubMed CentralPubMedGoogle Scholar
  17. Nishi CYM, Hungria M: Effects of soybean [Glycine max (L) Merrill] inoculation in a soil with established population of Bradyrhizobium with strains SEMIA 566, 586, 587, 5019, 5079 and 5080. Pesq Agropec Bras. 1996, 31: 359-368.Google Scholar
  18. Hungria M, Boddey LH, Santos MA, Vargas MAT: Nitrogen fixation capacity and nodule occupancy by Bradyrhizobium japonicum and B. elkanii strains. Biol Fertil Soils. 1998, 27: 393-399. 10.1007/s003740050449.Google Scholar
  19. Batista JSS, Hungria M, Barcellos FG, Ferreira MC, Mendes IC: Variability in Bradyrhizobium japonicum and B. elkanii seven years after introduction of both the exotic microsymbiont and the soybean host in a Cerrados soil. Microbial Ecol. 2007, 53: 270-284. 10.1007/s00248-006-9149-2.Google Scholar
  20. Mendes I, Hungria M, Vargas M: Establishment of bradyrhizobium japonicum and B elkanii strains in a Brazilian Cerrado oxisol. Biol Fertil Soils. 2004, 40 (1): 28-35. 10.1007/s00374-004-0739-1.Google Scholar
  21. Cregan PB, Keyser HH, Sadowsky MJ: Soybean genotype restricting nodulation of a previously unrestricted serocluster 123 bradyrhizobia. Crop Sci. 1989, 29 (2): 307-312. 10.2135/cropsci1989.0011183X002900020015x.Google Scholar
  22. Weber DF, Keyser HH, Uratsu SL: Serological distribution of Bradyrhizobium japonicum from United States soybean production areas. Agron J. 1989, 81: 786-789. 10.2134/agronj1989.00021962008100050018x.Google Scholar
  23. Vargas MAT, Mendes IC, Suhet AR, Peres JRR: Serological distribution of Bradyrhizobium japonicum from Brazilian “Cerrados” areas under soybean cultivation. Rev Microbiol. 1993, 24 (4): 239-243.Google Scholar
  24. Van Berkum P, Fuhrmann JJ: Evolutionary relationships among the soybean bradyrhizobia reconstructed from 16S rRNA gene and internally transcribed spacer region sequence divergence. Int J Syst Evol Microbiol. 2000, 50 (6): 2165-2172. 10.1099/00207713-50-6-2165.PubMedGoogle Scholar
  25. Keyser HH, Griffin RF: BeltsviIle Rhizobium Culture Collection Catalog. 1987, Beltsville, MD, USA: USDA-ARSGoogle Scholar
  26. Cytryn EJ, Jitacksorn S, Giraud E, Sadowsky MJ: Insights learned from pBTAi1, a 229-kb accessory plasmid from Bradyrhizobium sp. strain BTAi1 and prevalence of accessory plasmids in other Bradyrhizobium sp. strains. ISME J. 2008, 2 (2): 158-170. 10.1038/ismej.2007.105.PubMedGoogle Scholar
  27. Giraud E, Moulin L, Vallenet D, Barbe V, Cytryn E, Avarre JC, Jaubert M, Simon D, Cartieaux F, Prin Y, Bena G, Hannibal L, Fardoux J, Kojadinovic M, Vuillet L, Lajus A, Cruveiller S, Rouy Z, Mangenot S, Segurens B, Dossat C, Franck WL, Chang WS, Saunders E, Bruce D, Richardson P, Normand P, Dreyfus B, Pignol D, Stacey G, et al: Legumes symbioses: absence of nod genes in photosynthetic bradyrhizobia. Science. 2007, 316 (5829): 1307-1312. 10.1126/science.1139548.PubMedGoogle Scholar
  28. Okubo T, Fukushima S, Itakura M, Oshima K, Longtonglang A, Teaumroong N, Mitsui H, Hattori M, Hattori R, Hattori T, Minamisawa K: Genome analysis suggests that the soil oligotrophic bacterium Agromonas oligotrophica (Bradyrhizobium oligotrophicum) is a nitrogen-fixing symbiont of Aeschynomene indica. Appl Environ Microbiol. 2013, 79 (8): 2542-2551. 10.1128/AEM.00009-13.PubMed CentralPubMedGoogle Scholar
  29. Orozco-Mosqueda MC, Altamirano-Hernandez J, Farias-Rodriguez R, Valencia-Cantero E, Santoyo G: Homologous recombination and dynamics of rhizobial genomes. Res Microbiol. 2009, 160 (10): 733-741. 10.1016/j.resmic.2009.09.011.Google Scholar
  30. Althabegoiti MJ, Ormeño-Orrillo E, Lozano L, Torres Tejerizo G, Rogel M, Mora J, Martínez-Romero E: Characterization of Rhizobium grahamii extrachromosomal replicons and their transfer among rhizobia. BMC Microbiol. 2014, 14 (1): 6-10.1186/1471-2180-14-6.PubMed CentralPubMedGoogle Scholar
  31. Dobrindt U, Hacker J: Whole genome plasticity in pathogenic bacteria. Curr Opin Microbiol. 2001, 4 (5): 550-557. 10.1016/S1369-5274(00)00250-2.PubMedGoogle Scholar
  32. Batista JS, Torres AR, Hungria M: Towards a two-dimensional proteomic reference map of Bradyrhizobium japonicum CPAC 15: spotlighting “hypothetical proteins”. Proteomics. 2010, 10 (17): 3176-3189.33. 10.1002/pmic.201000092.PubMedGoogle Scholar
  33. Goedert WJ: Solos dos Cerrados: Tecnologias e Estratégias de Manejo. 1985, Brasília-DF: Embrapa-CPACGoogle Scholar
  34. Hungria M, Chueire LMO, Coca RG, Megías M: Preliminary characterization of fast growing rhizobial strains isolated from soyabean nodules in Brazil. Soil Biol Biochem. 2001, 33 (10): 1349-1361. 10.1016/S0038-0717(01)00040-2.Google Scholar
  35. Prell J, Poole P: Metabolic changes of rhizobia in legume nodules. Trends Microbiol. 2006, 14 (4): 161-168. 10.1016/j.tim.2006.02.005.PubMedGoogle Scholar
  36. Lodwig EM, Hosie AH, Bourdes A, Findlay K, Allaway D, Karunakaran R, Downie JA, Poole PS: Amino-acid cycling drives nitrogen fixation in the legume-Rhizobium symbiosis. Nature. 2003, 422 (6933): 722-726. 10.1038/nature01527.PubMedGoogle Scholar
  37. Prell J, White JP, Bourdes A, Bunnewell S, Bongaerts RJ, Poole PS: Legumes regulate Rhizobium bacteroid development and persistence by the supply of branched-chain amino acids. Proc Natl Acad Sci U S A. 2009, 106 (30): 12477-12482. 10.1073/pnas.0903653106.PubMed CentralPubMedGoogle Scholar
  38. D’hooghe I, VanderWauven C, Michiels J, Tricot C, De Wilde P, Vanderleyden J, Stalon V: The arginine deiminase pathway in Rhizobium etli: DNA sequence analysis and functional study of the arcABC genes. J Bacteriol. 1997, 179: 7403-7409.PubMed CentralPubMedGoogle Scholar
  39. Göttfert M, Rothlisberger S, Kundig C, Beck C, Marty R, Hennecke H: Potential symbiosis-specific genes uncovered by sequencing a 410-kilobase DNA region of the Bradyrhizobium japonicum chromosome. J Bacteriol. 2001, 183: 1405-1412. 10.1128/JB.183.4.1405-1412.2001.PubMed CentralPubMedGoogle Scholar
  40. Garcia M, Dunlap J, Loh J, Stacey G: Phenotypic characterization and regulation of the nolA gene of Bradyrhizobium japonicum. Mol Plant Microbe Interact. 1996, 9: 625-636. 10.1094/MPMI-9-0625.PubMedGoogle Scholar
  41. Godoy LP, Vasconcelos ATR, Chueire LMO, Souza RC, Nicolás MF, Barcellos FG, Hungria M: Genomic panorama of Bradyrhizobium japonicum CPAC 15, a commercial inoculant strain largely established in Brazilian soils and belonging to the same serogroup as USDA 123. Soil Biol Biochem. 2008, 40: 2742-2753.Google Scholar
  42. Göttfert M, Grob P, Hennecke H: Proposed regulatory pathway encoded by the nodV and nodW genes, determinants of host specificity in Bradyrhizobium japonicum. Proc Natl Acad Sci U S A. 1990, 87: 2680-2684. 10.1073/pnas.87.7.2680.PubMed CentralPubMedGoogle Scholar
  43. Loh J, Lohar DP, Andersen B, Stacey G: A two-component regulator mediates population-density-dependent expression of the Bradyrhizobium japonicum nodulation genes. J Bacteriol. 2002, 184: 1759-1766. 10.1128/JB.184.6.1759-1766.2002.PubMed CentralPubMedGoogle Scholar
  44. Grob P, Michel P, Hennecke H, Gottfert M: A novel response-regulator is able to suppress the nodulation defect of a Bradyrhizobium japonicum nodW mutant. Mol Gen Genet. 1993, 241: 531-541.PubMedGoogle Scholar
  45. Schubert KR, Evans HJ: Hydrogen evolution: a major factor affecting the efficiency of nitrogen fixation in nodulated symbionts. Proc Natl Acad Sci U S A. 1976, 73 (4): 1207-1211. 10.1073/pnas.73.4.1207.PubMed CentralPubMedGoogle Scholar
  46. Neves MCP, Hungria M: The physiology of nitrogen fixation in tropical grain legumes. CRC Crit Rev Plant Sci. 1987, 6 (3): 267-321. 10.1080/07352688709382252.Google Scholar
  47. Hungria M, Neves MCP, Dobereiner J: Relative efficiency, ureide transport and harvest index in soybeans inoculated with isogenic HUP mutants of Bradyrhizobium japonicum. Biology Fertil Soils. 1989, 7: 325-329.Google Scholar
  48. Sl A, Maier RJ, Hanus FJ, Russell SA, Emerich DW, Evans HJ: Hydrogenase in Rhizobium japonicum increases nitrogen fixation by nodulated soybeans. Science. 1979, 203 (4386): 1255-1257. 10.1126/science.203.4386.1255.Google Scholar
  49. Vignais PM, Billoud B: Occurrence, classification, and biological function of hydrogenases: an overview. Chem Rev. 2007, 107 (10): 4206-4272. 10.1021/cr050196r.PubMedGoogle Scholar
  50. Wiedenbeck J, Cohan FM: Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches. FEMS Microbiol Rev. 2011, 35 (5): 957-976. 10.1111/j.1574-6976.2011.00292.x.PubMedGoogle Scholar
  51. Ramsay JP, Sullivan JT, Stuart GS, Lamont IL, Ronson CW: Excision and transfer of the Mesorhizobium loti R7A symbiosis island requires an integrase IntS, a novel recombination directionality factor RdfS, and a putative relaxase RlxS. Mol Microbiol. 2006, 62 (3): 723-734. 10.1111/j.1365-2958.2006.05396.x.PubMedGoogle Scholar
  52. Dombrecht B, Vanderleyden J, Michiels J: Stable RK2-derived cloning vectors for the analysis of gene expression and gene function in Gram-negative bacteria. Mol Plant Microbe Interact. 2001, 14 (3): 426-430. 10.1094/MPMI.2001.14.3.426.PubMedGoogle Scholar
  53. Davidson AL, Dassa E, Orelle C, Chen J: Structure, function, and evolution of bacterial ATP-binding cassette systems. Microbiol Mol Biol Rev. 2008, 72 (2): 317-364. 10.1128/MMBR.00031-07.PubMed CentralPubMedGoogle Scholar
  54. Lippincott JA, Beiderbeck R, Lippincott BB: Utilization of octopine and nopaline by Agrobacterium. J Bacteriol. 1973, 116: 378-383.PubMed CentralPubMedGoogle Scholar
  55. Oger P, Mansouri H, Dessaux Y: Effect of crop rotation and soil cover on alteration of the soil microflora generated by the culture of transgenic plants producing opines. Mol Ecol. 2000, 9: 881-890. 10.1046/j.1365-294x.2000.00940.x.PubMedGoogle Scholar
  56. Gong ZZ, Lee H, Xiong LM, Jagendorf A, Stevenson B, Zhu JK: RNA helicase-like protein as an early regulator of transcription factors for plant chilling and freezing tolerance. Proc Natl Acad Sci U S A. 2002, 99: 11507-11512. 10.1073/pnas.172399299.PubMed CentralPubMedGoogle Scholar
  57. Barcellos F, Batista J, Menna P, Hungria M: Genetic differences between Bradyrhizobium japonicum variant strains contrasting in N2-fixation efficiency revealed by representational difference analysis. Arch Microbiol. 2009, 191 (2): 113-122. 10.1007/s00203-008-0432-0.PubMedGoogle Scholar
  58. Sirko A, Hryniewicz M, Hulanicka D, Bock A: Sulfate and thiosulfate transport in Escherichia coli K-12: nucleotide sequence and expression of the cysTWAM gene cluster. J Bacteriol. 1990, 172 (6): 3351-3357.PubMed CentralPubMedGoogle Scholar
  59. Fraysse N, Couderc F, Poinsot V: Surface polysaccharide involvement in establishing the rhizobium-legume symbiosis. Eur J Biochem. 2003, 270 (7): 1365-1380. 10.1046/j.1432-1033.2003.03492.x.PubMedGoogle Scholar
  60. Masson-Boivin C, Giraud E, Perret X, Batut J: Establishing nitrogen-fixing symbiosis with legumes: how many rhizobium recipes?. Trends Microbiol. 2009, 17 (10): 458-466. 10.1016/j.tim.2009.07.004.PubMedGoogle Scholar
  61. Mort AJ, Bauer WD: Composition of the capsular and extracellular polysaccharides of Rhizobium japonicum: Changes with culture age and correlations with binding of soybean seed lectin to the bacteria. Plant Physiol. 1980, 66 (1): 158-163. 10.1104/pp.66.1.158.PubMed CentralPubMedGoogle Scholar
  62. Mort AJ, Bauer WD: Application of two new methods for cleavage of polysaccharides into specific oligosaccharide fragments. Structure of the capsular and extracellular polysaccharides of Rhizobium japonicum that bind soybean lectin. J Biol Chem. 1982, 257 (4): 1870-1875.PubMedGoogle Scholar
  63. Dudman WF: The extracellular polysaccharides of Rhizobium japonicum: compositional studies. Carbohydr Res. 1976, 46 (1): 97-110. 10.1016/S0008-6215(00)83534-X.Google Scholar
  64. Becker BU, Kosch K, Parniske M, Muller P: Exopolysaccharide (EPS) synthesis in Bradyrhizobium japonicum: sequence, operon structure and mutational analysis of an exo gene cluster. Mol Gen Genet. 1998, 259 (2): 161-171. 10.1007/s004380050801.PubMedGoogle Scholar
  65. Quelas JI, Mongiardini EJ, Casabuono A, Lopez-Garcia SL, Althabegoiti MJ, Covelli JM, Perez-Gimenez J, Couto A, Lodeiro AR: Lack of galactose or galacturonic acid in Bradyrhizobium japonicum USDA 110 exopolysaccharide leads to different symbiotic responses in soybean. Mol Plant Microbe Interact. 2010, 23 (12): 1592-1604. 10.1094/MPMI-05-10-0122.PubMedGoogle Scholar
  66. McIntosh M, Stone BA, Stanisich VA: Curdlan and other bacterial (1→3)-beta-D-glucans. Appl Microbiol Biotechnol. 2005, 68 (2): 163-173. 10.1007/s00253-005-1959-5.PubMedGoogle Scholar
  67. Mithofer A: Suppression of plant defence in rhizobia-legume symbiosis. Trends Plant Sci. 2002, 7 (10): 440-444. 10.1016/S1360-1385(02)02336-1.PubMedGoogle Scholar
  68. Gay-Fraret J, Ardissone S, Kambara K, Broughton WJ, Deakin WJ, Le Quere A: Cyclic-beta-glucans of Rhizobium (Sinorhizobium) sp. strain NGR234 are required for hypo-osmotic adaptation, motility, and efficient symbiosis with host plants. FEMS Microbiol Lett. 2012, 333 (1): 28-36. 10.1111/j.1574-6968.2012.02595.x.PubMedGoogle Scholar
  69. Gore RS, Miller KJ: Cyclic [beta]-1,6–1,3 glucans are synthesized by Bradyrhizobium japonicum bacteroids within soybean (Glycine max) root nodules. Plant Physiol. 1993, 102 (1): 191-194. 10.1104/pp.102.1.191.PubMed CentralPubMedGoogle Scholar
  70. Miller KJ, Hadley JA, Gustine DL: Cyclic [beta]-1,6-1,3-glucans of Bradyrhizobium japonicum USDA 110 elicit isoflavonoid production in the soybean (Glycine max) host. Plant Physiol. 1994, 104 (3): 917-923.PubMed CentralPubMedGoogle Scholar
  71. Bohin JP: Osmoregulated periplasmic glucans in Proteobacteria. FEMS Microbiol Lett. 2000, 186 (1): 11-19. 10.1111/j.1574-6968.2000.tb09075.x.PubMedGoogle Scholar
  72. Bhagwat AA, Mithofer A, Pfeffer PE, Kraus C, Spickers N, Hotchkiss A, Ebel J, Keister DL: Further studies of the role of cyclic beta-glucans in symbiosis. An NdvC mutant of Bradyrhizobium japonicum synthesizes cyclodecakis-(1→3)-beta-glucosyl. Plant Physiol. 1999, 119 (3): 1057-1064. 10.1104/pp.119.3.1057.PubMed CentralPubMedGoogle Scholar
  73. Chen R, Bhagwat AA, Yaklich R, Keister DL: Characterization of ndvD, the third gene involved in the synthesis of cyclic β-(13), (16)-D-glucans in Bradyrhizobium japonicum. Can J Microbiol. 2002, 48 (11): 1008-1016. 10.1139/w02-099.PubMedGoogle Scholar
  74. Choma A, Komaniecka I: Straight and branched (omega-1)-hydroxylated very long chain fatty acids are components of Bradyrhizobium lipid A. Acta Biochim Pol. 2011, 58 (1): 51-57.PubMedGoogle Scholar
  75. Haag AF, Arnold MFF, Myka KK, Kerscher B, Dall’Angelo S, Zanda M, Mergaert P, Ferguson GP: Molecular insights into bacteroid development during Rhizobium–legume symbiosis. FEMS Microbiol Rev. 2013, 37 (3): 364-383. 10.1111/1574-6976.12003.PubMedGoogle Scholar
  76. Chang W-S, Park K-M, Koh S-C, So J-S: Characterization of the Bradyrhizobium japonicum galE gene: its impact on lipopolysaccharide profile and nodulation of soybean. FEMS Microbiol Lett. 2008, 280 (2): 242-249. 10.1111/j.1574-6968.2008.01066.x.PubMedGoogle Scholar
  77. Ormeño-Orrillo E: Lipopolisacáridos de Rhizobiaceae: estructura y biosíntesis. Rev Latinoam Microbiol. 2005, 47 (3–4): 165-175.PubMedGoogle Scholar
  78. Linhartová I, Bumba L, Masin J, Basler M, Osicka R, Kamanová J, Prochazkova K, Adkins I, Hejnova-Holubova J, Sadilkova L, Morová J, Sebo P: RTX proteins: a highly diverse family secreted by a common mechanism. FEMS Microbiol Rev. 2010, 34 (6): 1076-1112.PubMed CentralPubMedGoogle Scholar
  79. Oresnik IJ, Twelker S, Hynes MF: Cloning and characterization of a Rhizobium leguminosarum gene encoding a bacteriocin with similarities to RTX toxins. Appl Environ Microbiol. 1999, 65 (7): 2833-2840.PubMed CentralPubMedGoogle Scholar
  80. Barloy-Hubler F, Chéron A, Hellégouarch A, Galibert F: Smc01944, a secreted peroxidase induced by oxidative stresses in Sinorhizobium meliloti 1021. Microbiology. 2004, 150 (3): 657-664. 10.1099/mic.0.26764-0.PubMedGoogle Scholar
  81. Douzi B, Filloux A, Voulhoux R: On the path to uncover the bacterial type II secretion system. Philos Trans R Soc Lond B Biol Sci. 2012, 367 (1592): 1059-1072. 10.1098/rstb.2011.0204.PubMed CentralPubMedGoogle Scholar
  82. Deakin WJ, Broughton WJ: Symbiotic use of pathogenic strategies: rhizobial protein secretion systems. Nat Rev Micro. 2009, 7 (4): 312-320.Google Scholar
  83. Krause A, Doerfel A, Göttfert M: Mutational and transcriptional analysis of the type III secretion system of Bradyrhizobium japonicum. Mol Plant Microbe Interact. 2002, 15 (12): 1228-1235. 10.1094/MPMI.2002.15.12.1228.PubMedGoogle Scholar
  84. Tsukui T, Eda S, Kaneko T, Sato S, Okazaki S, Kakizaki-Chiba K, Itakura M, Mitsui H, Yamashita A, Terasawa K, Minamisawa K: The Type III secretion system of Bradyrhizobium japonicum USDA122 mediates symbiotic incompatibility with Rj2 soybean plants. Appl Environ Microbiol. 2013, 79 (3): 1048-1051. 10.1128/AEM.03297-12.PubMed CentralPubMedGoogle Scholar
  85. Juhas M, Crook DW, Hood DW: Type IV secretion systems: tools of bacterial horizontal gene transfer and virulence. Cell Microbiol. 2008, 10 (12): 2377-2386. 10.1111/j.1462-5822.2008.01187.x.PubMed CentralPubMedGoogle Scholar
  86. Juhas M, Crook DW, Dimopoulou ID, Lunter G, Harding RM, Ferguson DJ, Hood DW: Novel type IV secretion system involved in propagation of genomic islands. J Bacteriol. 2007, 189 (3): 761-771. 10.1128/JB.01327-06.PubMed CentralPubMedGoogle Scholar
  87. Jenkins J, Mayans O, Pickersgill R: Structure and evolution of parallel beta-helix proteins. J Struct Biol. 1998, 122 (1–2): 236-246.PubMedGoogle Scholar
  88. Thanassi DG, Stathopoulos C, Karkal A, Li H: Protein secretion in the absence of ATP: the autotransporter, two-partner secretion and chaperone/usher pathways of gram-negative bacteria (review). Mol Membr Biol. 2005, 22 (1–2): 63-72.PubMedGoogle Scholar
  89. Schell MA, Ulrich RL, Ribot WJ, Brueggemann EE, Hines HB, Chen D, Lipscomb L, Kim HS, Mrázek J, Nierman WC, DeShazer D: Type VI secretion is a major virulence determinant in Burkholderia mallei. Mol Microbiol. 2007, 64 (6): 1466-1485. 10.1111/j.1365-2958.2007.05734.x.PubMedGoogle Scholar
  90. Zheng J, Leung KY: Dissection of a type VI secretion system in Edwardsiella tarda. Mol Microbiol. 2007, 66 (5): 1192-1206. 10.1111/j.1365-2958.2007.05993.x.PubMedGoogle Scholar
  91. Records AR: The type VI secretion system: a multipurpose delivery system with a phage-like machinery. Mol Plant Microbe Interact. 2011, 24 (7): 751-757. 10.1094/MPMI-11-10-0262.PubMedGoogle Scholar
  92. Bladergroen MR, Badelt K, Spaink HP: Infection-blocking genes of a symbiotic Rhizobium leguminosarum strain that are involved in temperature-dependent protein secretion. Mol Plant Microbe Interact. 2003, 16 (1): 53-64. 10.1094/MPMI.2003.16.1.53.PubMedGoogle Scholar
  93. Pelicic V: Type IV pili: e pluribus unum?. Mol Microbiol. 2008, 68 (4): 827-837. 10.1111/j.1365-2958.2008.06197.x.PubMedGoogle Scholar
  94. Tomich M, Planet PJ, Figurski DH: The tad locus: postcards from the widespread colonization island. Nat Rev Microbiol. 2007, 5 (5): 363-375. 10.1038/nrmicro1636.PubMedGoogle Scholar
  95. Kanbe M, Yagasaki J, Zehner S, Goettfert M, Aizawa SI: Characterization of two sets of subpolar flagella in Bradyrhizobium japonicum. J Bacteriol. 2007, 189: 1083-1089. 10.1128/JB.01405-06.PubMed CentralPubMedGoogle Scholar
  96. Althabegoiti MJ, Covelli JM, Perez-Gimenez J, Quelas JI, Mongiardini EJ, Lopez MF, Lopez-Garcia SL, Lodeiro AR: Analysis of the role of the two flagella of Bradyrhizobium japonicum in competition for nodulation of soybean. FEMS Microbiol Lett. 2011, 319: 133-139. 10.1111/j.1574-6968.2011.02280.x.PubMedGoogle Scholar
  97. Miller MB, Bassler BL: Quorum sensing in bacteria. Annu Rev Microbiol. 2001, 55: 165-199. 10.1146/annurev.micro.55.1.165.PubMedGoogle Scholar
  98. Gonzalez JE, Marketon MM: Quorum sensing in nitrogen-fixing rhizobia. Microbiol Mol Biol Rev. 2003, 67: 574-592. 10.1128/MMBR.67.4.574-592.2003.PubMed CentralPubMedGoogle Scholar
  99. Sanchez-Contreras M, Bauer WD, Gao MS, Robinson JB, Downie JA: Quorum-sensing regulation in rhizobia and its role in symbiotic interactions with legumes. Philos Trans R Soc Lond B Biol Sci. 2007, 362: 1149-1163. 10.1098/rstb.2007.2041.PubMed CentralPubMedGoogle Scholar
  100. Pongsilp N, Triplett EW, Sadowsky MJ: Detection of homoserine lactone-like quorum sensing molecules in Bradyrhizobium strains. Curr Microbiol. 2005, 51: 250-254. 10.1007/s00284-005-4550-5.PubMedGoogle Scholar
  101. Lindemann A, Pessi G, Schaefer AL, Mattmann ME, Christensen QH, Kessler A, Hennecke H, Blackwell HE, Greenberg EP, Harwood CS: Isovaleryl-homoserine lactone, an unusual branched-chain quorum-sensing signal from the soybean symbiont Bradyrhizobium japonicum. Proc Natl Acad Sci U S A. 2011, 108: 16765-16770. 10.1073/pnas.1114125108.PubMed CentralPubMedGoogle Scholar
  102. Winkelmann G: Microbial siderophore-mediated transport. Biochem Soc Trans. 2002, 30 (4): 691-696.PubMedGoogle Scholar
  103. Guerinot ML, Meidl EJ, Plessner O: Citrate as a siderophore in Bradyrhizobium japonicum. J Bacteriol. 1990, 172: 3298-3303.PubMed CentralPubMedGoogle Scholar
  104. Lesueur D, Diem HG, Meyer JM: Iron requirement and siderophore production in Bradyrhizobium strains isolated from Acacia mangium. J Appl Microbiol. 1993, 74: 675-682.Google Scholar
  105. Nambiar PTC, Sivaramakrishnan S: Detection and assay of siderophores in cowpea rhizobia (Bradyrhizobium) using radioactive Fe. Lett Appl Microbiol. 1987, 4: 37-40. 10.1111/j.1472-765X.1987.tb01577.x.Google Scholar
  106. Datta B, Chakrabartty P: Siderophore biosynthesis genes of Rhizobium sp. isolated from Cicer arietinum L. Biotech. 2013, 1-11. doi:10.1007/s13205-013-0164-yGoogle Scholar
  107. Bossier P, Hofte M, Verstraete W: Ecological significance of siderophores in soil. Adv Microb Ecol. 1988, 10: 385-414. 10.1007/978-1-4684-5409-3_9.Google Scholar
  108. Geetha SJ, Joshi SJ: Engineering rhizobial bioinoculants: a strategy to improve iron nutrition. Sci World J. 2013, doi:10.1155/2013/315890Google Scholar
  109. Plessner O, Klapatch T, Guerinot ML: Siderophore utilization by Bradyrhizobium japonicum. Appl Environ Microbiol. 1993, 59: 1688-1690.PubMed CentralPubMedGoogle Scholar
  110. Ding Y, Oldroyd GED: Positioning the nodule, the hormone dictum. Plant Signal Behav. 2009, 4 (2): 89-93. 10.4161/psb.4.2.7693.PubMed CentralPubMedGoogle Scholar
  111. Ryu H, Cho H, Choi D, Hwang I: Plant hormonal regulation of nitrogen-fixing nodule organogenesis. Mol Cells. 2012, 34: 117-126. 10.1007/s10059-012-0131-1.PubMed CentralPubMedGoogle Scholar
  112. Malthesius U: Auxin: at the root of nodule development?. J Funct Plant Biol. 2008, 35 (8): 651-668. 10.1071/FP08177.Google Scholar
  113. Frugier F, Kosuta S, Murray JD, Crespi M, Szczyglowski K: Cytokinin: secret agent of symbiosis. Trends Plant Sci. 2008, 13: 115-120. 10.1016/j.tplants.2008.01.003.PubMedGoogle Scholar
  114. Frebort I, Kowalska M, Hluska T, Frebortova J, Galuszka P: Evolution of cytokinin biosynthesis and degradation. J Exp Bot. 2011, 62: 2431-2452. 10.1093/jxb/err004.PubMedGoogle Scholar
  115. Sturtevant DB, Taller BJ: Cytokinin production by Bradyrhizobium japonicum. Plant Physiol. 1989, 89: 1247-1252. 10.1104/pp.89.4.1247.PubMed CentralPubMedGoogle Scholar
  116. Boiero L, Perrig D, Masciarelli O, Penna C, Cassán F, Luna V: Phytohormone production by three strains of Bradyrhizobium japonicum and possible physiological and technological implications. Appl Microbiol Biotechnol. 2007, 74 (4): 874-880. 10.1007/s00253-006-0731-9.PubMedGoogle Scholar
  117. Podlesakova K, Fardoux J, Patrel D, Bonaldi K, Novak O, Strnad M, Giraud E, Spichal L, Nouwen N: Rhizobial synthesized cytokinins contribute to but are not essential for the symbiotic interaction between photosynthetic bradyrhizobia and Aeschynomene legumes. Mol Plant Microbe Interact. 2013, 26: 1232-1238. 10.1094/MPMI-03-13-0076-R.PubMedGoogle Scholar
  118. Sekine M, Ichikawa T, Kuga N, Kobayashi M, Sakurai A, Syono K: Detection of the IAA biosynthetic-pathway from tryptophan via indole-3-acetamide in Bradyrhizobium spp. Plant Cell Physiol. 1988, 29: 867-874.Google Scholar
  119. Sekine M, Watanabe K, Syono K: Molecular cloning of a gene for indole-3-acetamide hydrolase from Bradyrhizobium japonicum. J Bacteriol. 1989, 171: 1718-1724.PubMed CentralPubMedGoogle Scholar
  120. Vega-Hernández MC, Leon-Barrios M, Perez-Galdona R: Indole-3-acetic acid production from indole-3-acetonitrile in Bradyrhizobium. Soil Biol Biochem. 2002, 34: 665-668. 10.1016/S0038-0717(01)00229-2.Google Scholar
  121. Hunter WJ: Influence of 5-methyltryptophan-resistant Bradyrhizobium japonicum on soybean root nodule indole-3-acetic-acid content. Appl Environ Microbiol. 1987, 53: 1051-1055.PubMed CentralPubMedGoogle Scholar
  122. Hunter WJ: Indole-3-acetic-acid production by bacteroids from soybean root-nodules. Physiol Plant. 1989, 76: 31-36. 10.1111/j.1399-3054.1989.tb05448.x.Google Scholar
  123. Kaneshiro T, Kwolek WF: Stimulated nodulation of soybeans by Rhizobium japonicum mutant (B-14075) that catabolizes the conversion of tryptophan to indol-3yl-acetic acid. Plant Sci. 1985, 42 (3): 141-146.124. 10.1016/0168-9452(85)90119-0.Google Scholar
  124. Spaepen S, Vanderleyden J: Auxin and plant-microbe interactions. Cold Spring Harb Perspect Biol. 2011, 3: a001438-PubMed CentralPubMedGoogle Scholar
  125. Spaepen S, Vanderleyden J, Remans R: Indole-3-acetic acid in microbial and microorganism-plant signaling. FEMS Microbiol Rev. 2007, 31 (4): 425-448. 10.1111/j.1574-6976.2007.00072.x.PubMedGoogle Scholar
  126. Patten CL, Blakney AJC, Coulson TJD: Activity, distribution and function of indole-3-acetic acid biosynthetic pathways in bacteria. Crit Rev Microbiol. 2013, 39: 395-415. 10.3109/1040841X.2012.716819.PubMedGoogle Scholar
  127. Lehmann T, Hoffmann M, Hentrich M, Pollmann S: Indole-3-acetamide-dependent auxin biosynthesis: a widely distributed way of indole-3-acetic acid production?. Eur J Cell Biol. 2010, 89: 895-905. 10.1016/j.ejcb.2010.06.021.PubMedGoogle Scholar
  128. Zhu D, Mukherjee C, Yang Y, Rios BE, Gallagher DT, Smith NN, Biehl ER, Hua L: A new nitrilase from Bradyrhizobium japonicum USDA 110: Gene cloning, biochemical characterization and substrate specificity. J Biotechnol. 2008, 133 (3): 327-333. 10.1016/j.jbiotec.2007.10.001.PubMedGoogle Scholar
  129. Antoun H, Beauchamp CJ, Goussard N, Chabot R, Lalande R: Potential of Rhizobium and Bradyrhizobium species as plant growth promoting rhizobacteria on non-legumes: effect on radishes (Raphanus sativus L.). Plant Soil. 1998, 204: 57-67. 10.1023/A:1004326910584.Google Scholar
  130. Chaintreuil C, Giraud E, Prin Y, Lorquin J, Ba A, Gillis M, de Lajudie P, Dreyfus B: Photosynthetic bradyrhizobia are natural endophytes of the African wild rice Oryza breviligulata. Appl Environ Microbiol. 2000, 66 (12): 5437-5447. 10.1128/AEM.66.12.5437-5447.2000.PubMed CentralPubMedGoogle Scholar
  131. Rouws LFM, Leite J, de Matos GF, Zilli JE, Coelho MRR, Xavier GR, Fischer D, Hartmann A, Reis VM, Baldani JI: Endophytic Bradyrhizobium spp. isolates from sugarcane obtained through different culture strategies. Environ Microbiol Rep. 2013, n/a-n/a. doi:10.1111/1758-2229.12122Google Scholar
  132. Fischer D, Pfitzner B, Schmid M, Simoes-Araújo JL, Reis VM, Pereira W, Ormeño-Orrillo E, Hai B, Hofmann A, Schloter M, Martinez-Romero E, Baldani JI, Hartmann A: Molecular characterisation of the diazotrophic bacterial community in uninoculated and inoculated field-grown sugarcane (Saccharum sp.). Plant Soil. 2012, 356 (1–2): 83-99.Google Scholar
  133. Donati AJ, Lee HI, Leveau JHJ, Chang WS: Effects of indole-3-acetic acid on the transcriptional activities and stress tolerance of Bradyrhizobium japonicum. PLoS One. 2013, 8 (10): e76559-10.1371/journal.pone.0076559. 134PubMed CentralPubMedGoogle Scholar
  134. O’Brien J, Wright GD: An ecological perspective of microbial secondary metabolism. Curr Opin Biotechnol. 2011, 22: 552-558. 10.1016/j.copbio.2011.03.010.PubMedGoogle Scholar
  135. Davies J: Secondary metabolites - their function and evolution - introduction. Ciba Found Symp. 1992, 171: 1-2.Google Scholar
  136. Demain AL: Microbial natural products: alive and well in 1998. Nature Biotechnol. 1998, 16 (1): 3-4. 10.1038/nbt0198-3.Google Scholar
  137. Park K-K, Chun K-S, Lee J-M, Lee SS, Surh Y-J: Inhibitory effects of [6]-gingerol, a major pungent principle of ginger, on phorbol ester-induced inflammation, epidermal ornithine decarboxylase activity and skin tumor promotion in ICR mice. Cancer Lett. 1998, 129 (2): 139-144. 10.1016/S0304-3835(98)00081-0.PubMedGoogle Scholar
  138. Zhao LQ, Chen YL, Ren S, Han Y, Cheng HB: Studies on the chemical structure and antitumor activity of an exopolysaccharide from Rhizobium sp N613. Carbohydr Res. 2010, 345: 637-643. 10.1016/j.carres.2009.11.017.PubMedGoogle Scholar
  139. Cushnie TPT, Lamb AJ: Antimicrobial activity of flavonoids. Int J Antimicrob Agents. 2005, 26 (5): 343-356. 10.1016/j.ijantimicag.2005.09.002.PubMedGoogle Scholar
  140. Niemeyer HM: Hydroxamic acids (4-hydroxy-1,4-benzoxazin-3-ones), defence chemicals in the gramineae. Phytochemistry. 1988, 27 (11): 3349-3358. 10.1016/0031-9422(88)80731-3.Google Scholar
  141. Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R: antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011, 39 (suppl. 2): W339-W346.PubMed CentralPubMedGoogle Scholar
  142. Dowling DN, Broughton WJ: Competition for nodulation of legumes. Annu Rev Microbiol. 1986, 40: 131-157. 10.1146/annurev.mi.40.100186.001023.PubMedGoogle Scholar
  143. Almeida LG, Paixao R, Souza RC, Costa GC, Barrientos FJ, Santos MT, Almeida DF, Vasconcelos AT: A system for automated bacterial (genome) integrated annotation–SABIA. Bioinformatics. 2004, 20 (16): 2832-2833. 10.1093/bioinformatics/bth273.PubMedGoogle Scholar
  144. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008, 9: 75-10.1186/1471-2164-9-75.PubMed CentralPubMedGoogle Scholar
  145. Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5 (2): R12-10.1186/gb-2004-5-2-r12.PubMed CentralPubMedGoogle Scholar
  146. Darling AE, Mau B, Perna NT: progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010, 5 (6): e11147-10.1371/journal.pone.0011147.PubMed CentralPubMedGoogle Scholar
  147. Overbeek R, Fonstein M, D’Souza M, Pusch GD, Maltsev N: The use of gene clusters to infer functional coupling. Proc Natl Acad Sci U S A. 1999, 96: 2896-2901. 10.1073/pnas.96.6.2896.PubMed CentralPubMedGoogle Scholar
  148. Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedGoogle Scholar
  149. Oliveros JC: VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007, http://bioinfogp.cnb.csic.es/tools/venny/index.html,Google Scholar
  150. Pirooznia M, Nagarajan V, Youping D: GeneVenn - a web application for comparing gene lists using Venn diagrams. Bioinformation. 2007, 1 (10): 420-422. 10.6026/97320630001420.PubMed CentralPubMedGoogle Scholar
  151. Saier MH, Tran CV, Barabote RD: TCDB: the Transporter Classification Database for membrane transport protein analyses and information. Nucleic Acids Res. 2006, 34 (Database issue): D181-D186.PubMed CentralPubMedGoogle Scholar
  152. Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer ELL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38 (Database issue): D211-D222.PubMed CentralPubMedGoogle Scholar
  153. Eddy SR: Accelerated profile HMM searches. PLoS Comp Biol. 2011, 7 (10): e1002195-10.1371/journal.pcbi.1002195.Google Scholar
  154. Langille MGI, Brinkman FSL: IslandViewer: an integrated interface for computational identification and visualization of genomic islands. Bioinformatics. 2009, 25 (5): 664-665. 10.1093/bioinformatics/btp030.PubMed CentralPubMedGoogle Scholar
  155. Bi D, Liu L, Tai C, Deng Z, Rajakumar K, Ou H-Y: SecReT4: a web-based bacterial type IV secretion system resource. Nucleic Acids Res. 2013, 41 (D1): D660-D665. 10.1093/nar/gks1248.PubMed CentralPubMedGoogle Scholar

Copyright

© Siqueira et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.