Levels of Periplasmic Nitrate Reductase during Denitrification are Lower in Bradyrhizobium japonicum than in Bradyrhizobium diazoefficiens.

Soybean plants host endosymbiotic dinitrogen (N2)-fixing bacteria from the genus Bradyrhizobium. Under oxygen-limiting conditions, Bradyrhizobium diazoefficiens and Bradyrhizobium japonicum perform denitrification by sequentially reducing nitrate (NO3-) to nitrous oxide (N2O) or N2. The anaerobic reduction of NO3- to N2O was previously shown to be lower in B. japonicum than in B. diazoefficiens due to impaired periplasmic nitrate reductase (Nap) activity in B. japonicum. We herein demonstrated that impaired Nap activity in B. japonicum was due to low Nap protein levels, which may be related to a decline in the production of FixP and FixO proteins by the cbb3-type oxidase.

Soybean is a globally important leguminous crop that generally hosts endosymbiotic dinitrogen (N2)-fixing bacteria from the genus Bradyrhizobium (Argaw, 2014). Under oxygen-limiting conditions, some species of soybean-associated Bradyrhizobium may perform denitrification as an alternative respiratory process and sequentially reduce nitrate (NO3–) or nitrite (NO2–) to nitric oxide (NO), nitrous oxide (N2O), and N2 gases. Denitrification is functional in free-living bradyrhizobia at the soybean rhizosphere and in symbiotic bradyrhizobia inside the root nodules (Mesa ; Sameshima-Saito ; Sánchez ; Inaba ).

Bradyrhizobium diazoefficiens (reclassified from Bradyrhizobium japonicum [Delamuta ]; type strain USDA 110) may completely reduce NO3– to N2. These reactions depend on napEDABC, which encodes periplasmic NO3– reductase (Nap); nirK, which encodes copper-containing NO2– reductase (NirK); norCBQD, which encodes c-type NO reductase (Nor); and nosRZDYFLX, which encodes N2O reductase (Nos) (Sánchez ). While B. diazoefficiens is a complete denitrifier, B. japonicum (type strain USDA 6) lacks the nos gene cluster and is, thus, unable to reduce N2O to N2 (Sameshima-Saito ; Itakura ; Kaneko ; Siqueira ).

We previously screened the growth of 11 strains of B. japonicum and 15 strains of B. diazoefficiens and found that anaerobic (≈0% O2) growth by B. japonicum with NO3– as the electron acceptor was significantly lower than that by B. diazoefficiens; however, no significant differences were observed between the growth rates of B. japonicum and B. diazoefficiens strains under microaerobic (≈2% O2) and aerobic (≈12% O2) growth conditions in the presence of NO3– (Siqueira ). The lower growth rate of B. japonicum in anaerobiosis was not explained by the absence of nos, but by markedly impaired Nap activity that had a negative impact on the reduction of NO3– to N2O (Siqueira ). B. japonicum conserved the FixLJ–FixK2 regulatory cascade, which mediates the response of the nap operon to low oxygen and NO3–; the binding sites for FixK/FNR (fumarate and nitrate reductase) regulators upstream of napE were also conserved (Robles ; Bueno ; Siqueira ). Additionally, napA transcript levels were similar in B. japonicum and B. diazoefficiens (Siqueira ). Thus, we proposed that impaired Nap activity in B. japonicum may be due to posttranscriptional events (Siqueira ).

In the present study, we used the type strains B. diazoefficiens USDA 110 (Bd-USDA 110) and B. japonicum USDA 6 (Bj-USDA 6) to expand on our earlier research on the mechanisms responsible for impaired Nap activity in B. japonicum under denitrifying conditions.

Bradyrhizobium cells were precultured at 30°C in HM salt medium (Cole and Elkan, 1973) supplemented with 0.1% L-(+)-arabinose and 0.25% (w/v) yeast extract. HM medium supplemented with trace metals (Sameshima-Saito ) and 10‍ ‍mM KNO3 (HMMN medium) was used in all assays. Escherichia coli cells were grown at 37°C in Luria–Bertani medium (Miller, 1972). The following antibiotics were used for the Bradyrhizobium culture: tetracycline (Tc, 100‍ ‍μg mL–1), spectinomycin (Sp, 100‍ ‍μg mL–1), streptomycin (Sm, 100‍ ‍μg mL–1), and polymyxin B (50‍ ‍μg mL–1). Tc (50‍ ‍μg mL–1), Sp (50‍ ‍μg mL–1), and Sm (50‍ ‍μg mL–1) were used for the E. coli culture.

In growth experiments, precultured cells were inoculated into 34-mL test tubes containing 5‍ ‍mL of HMMN medium. The initial optical density at 660 nm (OD660) was adjusted to 0.01. In the competition experiment, initial cell numbers were adjusted to 3×107‍ ‍cells‍ ‍mL–1 and cells were used at a 1:1 ratio. Aerobic (≈12% O2), microaerobic (≈2% O2), and anaerobic (≈0% O2) treatments were prepared and monitored as previously described (Siqueira ). Cells were grown at 30°C with reciprocal shaking at 300 rpm. Growth was monitored daily by measuring the OD660 of the cultures, and the number of cells was directly counted using a 20-μm-deep hemocytometer (Sunlead Glass) and BX51 Fluorescence Microscope (Olympus).

To construct Bd-USDA 110 tagged with GFP and Bj-USDA 6 tagged with DsRed, pRJPaph-gfp (Ledermann ) and pBjGroEL4::DsRed2 (Hayashi ) plasmids were transferred by triparental mating using pRK2013 as a helper plasmid (Figurski and Helinski, 1979). Single recombination events, indicating the integration of pRJPaph-gfp into the Bd-USDA 110 chromosome or pBjGroEL4::DsRed2 into the Bj-USDA 6 chromosome, were selected based on their acquisition of resistance to Tc or Sp/Sm, respectively. Single recombinant strains were further confirmed using the BX51 Fluorescence Microscope (Olympus).

In immunoblotting and heme-staining analyses, precultured Bradyrhizobium cells were inoculated into 500-mL flasks containing 100‍ ‍mL of HMNN (OD660≈0.03) and incubated under microaerobic or anaerobic conditions at 30°C, with reciprocal shaking at 100 rpm. Soluble and membrane fractions were prepared as described previously (Fischer ; Torres ). Soluble fractions were concentrated using Amicon Ultra-15 centrifugal filter units (Merck). In immunoblotting experiments, rabbit polyclonal antibodies were raised against the synthetic NapA peptide, CPLQKKATGAAKAND (GenScript Japan), which is common between the Nap of Bd-USDA 110 and Bj-USDA 6. Soluble fractions (10‍ ‍μg) were separated by SDS-PAGE (10% polyacrylamide) and transferred to P plus PVDF membranes (ATTO) using a HorizeBLOT 4M-R semi-dry transfer system (ATTO). Membranes were incubated with an anti-NapA antibody (1:10,000) at room temperature for 30‍ ‍min and then with HRP-conjugated anti-rabbit IgG (1:2,500; Promega) for 30‍ ‍min. Immunoreactive bands were detected using the chemiluminescent substrate, PromegaTM ECL Western Blotting Substrate. Regarding heme staining, membrane fractions (10‍ ‍μg) were separated by SDS-PAGE at 4°C (stacking phase at 5% acrylamide and resolving phase at 12% acrylamide), transferred to a 0.45 μm Hybond-P membrane (GE Healthcare), and stained with the PromegaTM ECL Western Blotting Substrate to detect heme-dependent activity by chemiluminescence (Vargas ). Blots were visualized using a Molecular Imager Gel Doc XR+ System with Image Lab software (Bio-Rad). Relative band intensities were calculated using ImageJ software (Schneider ). Experiments were performed with at least two independent biological replicates.

The alignment of the napEDABC genome regions of Bj-USDA 6 (accession number NC_017249) and Bd-USDA 110 (accession number NC_004463) was performed using Genome Matcher version 2.2 (Ohtsubo ). Amino acid alignment was conducted using MEGA version 7.0 and BoxShade version 3.21 (https://www.ch.embnet.org/software/BOX_form.html). A protein structural similarity analysis was performed using TM-align version 20170708 (Zhang and Skolnick, 2005).

The impaired activity of Nap in B. japonicum during anaerobic growth in the presence of NO3– may rely on posttranscriptional effects that affect the structure, function, or level of Nap (Siqueira ). The napEDABC gene cluster of B. diazoefficiens encodes NapE, a membrane protein of 6.6‍ ‍kDa with unknown function; NapD (11.8‍ ‍kDa), a chaperone involved in the maturation of NapA; NapA (94.5‍ ‍kDa), the catalytic subunit that contains a Mobis-MGD active site and [4Fe4S] cluster; and NapB (16.9‍ ‍kDa), the electron transfer subunit, which is a diheme cytochrome c552. NapA and NapB form the heterodimeric subunit NapAB, which is present in the periplasm. NapC (25‍ ‍kDa) is a membrane-anchored, tetraheme, c-type cytochrome that is responsible for the transfer of electrons from the quinol pool to the NapAB complex (Zumft, 1997; Delgado ). A sequence comparison of the napEDABC gene cluster of Bd-USDA 110 and Bj-USDA 6 showed high identity (>85%) at the nucleotide level (Fig. S1). The amino acid sequences of Bj-USDA 6 NapEDABC were conserved, exhibiting 89–97% identity with Bd-USDA 110 (Fig. S1). In NapA of Bj-USDA 6, the twin arginine motif, the residues involved in the binding of the Mobis-MGD cofactor, and the [4Fe4S] cluster were conserved, except in some cases in which there was a substitution with a similar residue. In Bj-USDA 6 NapB, both the cytochrome c-binding sites and motif for translocation via the general secretory pathway were conserved. The structural similarity of the Bj-USDA 6 and Bd-USDA 110 Nap proteins was high, with a TM-score of 0.81, which indicated the same type of folding (Fig. S1). These results suggested that the structure and function of Nap were conserved in Bj-USDA 6.

We compared the amount of Nap in Bj-USDA 6 and Bd-USDA 110 cells incubated under microaerobic and anaerobic conditions in the presence of NO3–. An immunoblotting analysis of soluble fractions revealed a bradyrhizobial NapA protein band at ≈100‍ ‍kDa (Fig. 1A), which was consistent with the expected molecular weight of 94.5‍ ‍kDa (Delgado ). This was further confirmed by the absence of this band in soluble fractions of the B. diazoefficiens ΔnapA mutant (Fig. 1A). The intensity of the NapA band was weaker in Bj-USDA 6 than in Bd-USDA 110, in cells incubated microaerobically (≈0.57 times; Fig. 1A and Table S1), and particularly in cells incubated anaerobically (≈0.33 times; Fig. 1A and Table S1). This result suggested that Bj-USDA 6 produced less NapA than Bd-USDA 110 under low-oxygen conditions (i.e. ≤2% O2) in the presence of NO3–. Additionally, we detected NapC in the membrane fractions of bradyrhizobial cells by heme staining; the staining of heme covalently bound to c-type cytochromes (Vargas ). Bj-USDA 6 and Bd-USDA 110 showed the typical profile of five stained bands previously identified in membrane fractions of Bd-USDA 110 (Fig. 1B) (Preisig ; Mesa ; Delgado ; Bueno ; Torres ). The 32-kDa and 28-kDa bands, almost co-migrating, corresponded to the FixP and FixO proteins, respectively, of cbb3-type high affinity cytochrome oxidase, encoded by the fixNOQP operon (Preisig , 1996); the 25-kDa band corresponded to NapC (Delgado ); the 20-kDa band corresponded to cytochrome CycM (Bott ); and the 16-kDa band corresponded to the NorC subunit of Nor (Mesa ). We confirmed the identities of bands by analyzing soluble fractions of the B. diazoefficiens ΔnapA mutant under microaerobic conditions; bands correspondig to NapC and NorC were absent (Fig. 1B), as described previously (Delgado ). Under anaerobic conditions, the intensity of the NapC band was ≈0.3-fold weaker in Bj-USDA 6 than in Bd-USDA 110 (Fig. 1B and Table S1). Thus, immunoblotting and heme-staining results indicated that Bj-USDA 6 produced lower levels of Nap under low-oxygen conditions in the presence of NO3–. Therefore, impaired Nap activity in B. japonicum during NO3–-dependent anaerobic growth, which was reported previously (Siqueira ), may be caused by a low level of Nap.

Fig. 1.

Detection of Nap in Bradyrhizobium diazoefficiens USDA 110 (Bd) and Bradyrhizobium japonicum USDA 6 (Bj) in HMMN medium under microaerobic and anaerobic conditions. (A) Immunoblotting analysis of NapA in soluble fractions (10‍ ‍μg). (B) Heme staining of membrane fractions (10‍ ‍μg). In panels A and B, the B. diazoefficiens USDA 110 ΔnapA mutant grown in HMMN medium under microaerobic conditions is shown as the control.

Based on the results of heme staining, the amount of NorC was higher in Bj-USDA 6 than in Bd-USDA 110 under microaerobic and anaerobic conditions (≈1.3- and 1.6-fold, respectively) in the presence of NO3– (Fig. 1B and Table S1). This result was supported by a previous finding showing that norB transcript levels were higher in Bj-USDA 6 than in Bd-USDA 110 under anaerobic NO3–-dependent growth (Siqueira ). Heme staining also revealed that besides NapC, FixP and FixO were the only other cytochromes affected in Bj-USDA 6 under low-oxygen conditions; the amounts of FixP and FixO were lower in Bj-USDA 6 than in Bd-USDA 110 by ≈0.45-fold under microaerobic conditions and ≈0.8-fold under anaerobic conditions (Fig. 1B and Table S1). This result suggested that cbb3 oxidase is a critical enzyme for the adaptation process to NO3– respiration under low-oxygen conditions in Bradyrhizobium. Accordingly, the maximal expression of Nap does not appear to occur until the oxygen concentration becomes very low, and this is only observed after cbb3 oxidase has consumed the oxygen present in the growth medium (Bueno ). Since the transcription of nap genes appeared to be unaffected in B. japonicum cells grown anaerobically in the presence of NO3– (Siqueira ), the remaining oxygen in the growth medium may affect the amount of the Nap protein at the posttranscriptional level by inhibiting the translation of nap messenger RNA or inducing the degradation of Nap proteins.

The mechanisms responsible for low levels of cbb3 oxidase in B. japonicum under low-oxygen conditions in the presence of NO3– may be related to B. japonicum being unable to make an effective switch to denitrification in the absence of oxygen, similar to other bacteria (Aida ; Bergaust ; Torres ; Siqueira ). We tested the NO3–-dependent anaerobic growth of bradyrhizobial cells preincubated microaerobically, instead of aerobically as previously reported (Siqueira ). The microaerobic preincubation resulted in a reduced growth rate of Bj-USDA 6, to a similar extent as the aerobic preincubation (Fig. S2), indicating the absence of a significant difference between a rapid (aerobic preincubation) and gradual (microaerobic preincubation) transition from aerobiosis to anaerobiosis. However, further studies are needed to elucidate the mechanisms responsible for the changes observed in cbb3 oxidase levels in B. japonicum under denitrifying conditions.

Our previous findings prompted the hypothesis that B. japonicum may be less competitive than B. diazoefficiens due to energy depletion under anaerobic denitrifying growth (Siqueira ). We examined the competitive growth of Bd-USDA 110 and Bj-USDA 6 tagged with GFP and DsRed proteins, respectively. The tagged strains were tested to exclude the possible effects of tagging on growth rates (Fig. S3). Under microaerobic and anaerobic conditions, in the presence of NO3–, the growth of DsRed-tagged Bj-USDA 6 was lower than that of GFP-tagged Bd-USDA 110, whereas the aerobic growth rates of these two species were similar (Fig. 2). This result indicated that Bd-USDA 110 is more competitive than Bj-USDA 6 at low-oxygen levels in the presence of NO3–. In Bradyrhizobium, growth at low-oxygen levels in the presence of NO3– is supported by both cbb3 oxidase, which is active during free-living microaerobic growth and in N2-fixing bacteroids, and the denitrification pathway (Preisig , 1996; Bueno ). The weaker competitiveness of B. japonicum USDA 6 under low-oxygen conditions in the presence of NO3– may correlate with low levels of Nap (Fig. 1 and Table S1).

Fig. 2.

Average number of cells in co-cultures (1:1) of Bradyrhizobium diazoefficiens USDA 110 tagged with GFP (white bars) and Bradyrhizobium japonicum USDA 6 tagged with DsRed (grey bars) under anaerobic, microaerobic, and aerobic conditions in HMMN medium. Error bars indicate SE. * Values significantly different between USDA 110 and USDA 6 (t-test, P<0.05; n=3).

In conclusion, the present study demonstrated that the level of Nap under low-oxygen conditions in the presence of NO3– was lower in B. japonicum USDA 6 than in B. diazoefficiens USDA 110 and suggests that the capacity to maintain a sufficient quantity of Nap is an advantage under denitrifying conditions. The lower levels of Nap in B. japonicum USDA 6 presumably resulted in a competitive disadvantage against B. diazoefficiens USDA 110 due to energy depletion during denitrifying growth. Ecologically, this advantage of B. diazoefficiens may be an important factor that influences the predominance of B. diazoefficiens over B. japonicum in soils exposed to low-oxygen conditions (Shiina ).

Citation

Siqueira, A. F., Sugawara, M., Arashida, H., Minamisawa, K., and Sánchez, C. (2020) Levels of Periplasmic Nitrate Reductase during Denitrification are Lower in Bradyrhizobium japonicum than in Bradyrhizobium diazoefficiens. Microbes Environ 35: ME19129.

https://doi.org/10.1264/jsme2.ME19129

Supplementary Material