Contrasting Expression of Rhizobial Phytase in Nodules of Two Soybean Cultivars Grown Under Low Phosphorus Availability

Phosphorus deficiency can be a major limitation to legume growth when plant nitrogen nutrition depends on symbiotic nitrogen fixation. One possible approach to overcome this constraint is the selection of plant and rhizobial genotypes capable of metabolizing complex forms of phosphorus in the nodules. The aim of this research was to study the rhizobial phytase transcript abundance in nodules of two soybean cultivars (Glycine max (L.) Merr.) grown under two different phosphorus conditions in hydroaeroponic conditions. An in situ RT-PCR of a rhizobial phytase was performed in microtome sections of soybean nodules of two cultivars growing under phosphorus sufficiency vs. phosphorus deficiency. The results showed that the plant cultivar may influence the level of transcript abundance of the bacterial phytase and in consequence affect the phosphorus use efficiency of nitrogen-dependent Bradyrhizobium spp.-soybean symbioses. Thus, the selection of a good combination of plant and rhizobial genotypes should be a priority when breeding for phosphorus deficiency is performed.


INTRODUCTION
Phosphorus (P) is the second limiting element for plant growth after nitrogen (N), thus being an essential nutrient for agricultural production. Total soil P, in both inorganic and organic forms, is abundant in most soils (Richardson, 2001). However, only 1-5% of the total P is dissolved in the soil solution as orthophosphate ions, HPO 2− 4 or H 2 PO − 4 , and absorbed by plants (Richardson and Simpson, 2011;Khan et al., 2013). Despite being often overlooked, soil organic P (Po) comprises 30-50% of soil total P (Hinsinger et al., 2015), whether as part of the biomass or associated to the soil organic matter (Maougal et al., 2014a). Phytates, in particular, may represent up to 80% of the total Po (Quiquampoix and Mousain, 2005). However, organic forms of P are unavailable to plants unless mineralization takes place (Lazali et al., 2013). Phytases are enzymes responsible for the mineralization of phytates. They are classified according to their hydrolytic mechanism, the substrate specificity, their protein structure and optimal pH: histidine acid phosphatase (HAP) (EC 3.1.3.8), β-propeller phytase (BPP) (EC 3.1.3.8), cystein phytase (CPhy) and purple acid phytase (PAP) (EC 3.1.3.2) (Mullaney and Ullah, 2003;Jorquera et al., 2008a). There are a lot of rhizospheric bacterial genera that can mineralize phytate such as Enterobacter, Pseudomonas, Burkholderia, Bacillus and Klebsiella, among others (Yoon et al., 1996;Richardson and Hadobas, 1997;Richardson, 2001;Hill et al., 2007;Jorquera et al., 2008b;Maougal et al., 2014a;Ramesh et al., 2014). Several plant species have the ability to increase the activity of phytases in the rhizosphere in response to P deficiency (Li et al., 1997). P has a crucial role for symbiotic nitrogen fixation (SNF) and the subsequent ammonium assimilation, which take place in legume nodules; both processes are highly energy-consuming and depend on the energy status of nodules (Araújo et al., 2008;Lazali et al., 2016). P deficiency comprises a limitation for SNF in legumes (Graham et al., 2003), resulting in a decrease in SNF and, consequently, in slow plant growth (Israel, 1987). Previous reports showed that mutations on three Sinorhizobium meliloti P transport systems, OrfA-Pit, PstSCAB, and phoCDET, failure to obtain sufficient P for growth during the infection process, forming "empty nodules, " which contain very few infected cells, in the alfalfa (Medicago sativa) symbiont and fail to fix N 2 (Bardin et al., 1996;Yuan et al., 2006). Atlhough, the expression of one of these transports could be necessary and sufficient for SNF in alfalfa (Yuan et al., 2006). The content of P is greater in nodules than in other plant organs under conditions of either P sufficiency or deficiency (Tang et al., 2001;Schulze et al., 2006). Legumes allocate different type of P sources in nodules, and leaves start to deplete before nodules, when they grow under P deficiency (Sulieman et al., 2010). It has been reported for different legumes, such as M. sativa, Medicago truncatula, soybean (Glycine max), common bean (Phaseolus vulgaris), pea (Pisum sativum), and Sophora flavescens, that nodules become a P sink under P deficient growing conditions (Sa and Israel, 1991;Kouas et al., 2005;Sulieman et al., 2013a;Hu et al., 2018).
The present study analyzes this last approach, aiming to test if one of the two strains used to inoculate Uruguayan soybean (SEMIA 587) could mineralize the phytate present in soybean nodules, increasing P content and therefore the efficiency of the SNF. In Uruguay, SNF has been upheld as a long-term public policy since the 1960's (Altier et al., 2013). As a result of research and outreach policies, 100% of the farmers have adopted the use of inoculation technology (Lindrström et al., 2010). In Uruguay, there are two officially recommended strains for soybean inoculation, Bradyrhizobium elkanii strain SEMIA 587 and Bradyrhizobium elkanii strain SEMIA 5019 (syn. 29w) (Olivera et al., 2016). Specific references of these two strains can be found in Hungria et al. (1998). First, we tested the ability of the strain SEMIA 587 to mineralize phytate in vitro and afterwards, we studied the rhizobial phytase transcript abundance in soybean nodules for a better understanding of the role of this enzyme in the acquisition and cycling of phytate in soybean. We also addressed whether the transcript abundance levels significantly vary between two cultivars.

Rhizobial HAP Phytase Gene Primers Design
The design of the rhizobial HAP phytase primers was performed online at the National Center of Biotechnology Information (NCBI, http://blast.ncbi.nlm.nih.gov/Blast.cgi). The highest hit by BLASTp to rhizobiaceae group of a known amino acid sequence of the Escherichia coli HAP phytase (GenBank accession number: AAN28334.1) was obtained with a conserved hypothetical HAP-type phytase of Azorhizobium caulinodans (GenBank accession number: 5692057 AZC_4319). The gene sequence of this hypothetical HAP-type phytase was obtained from A. caulinodans full genome (GenBank accession number: AP009384.1). Thereafter, the following pair of primers were designed with the primer-BLAST tool from this HAP gene-sequence: forward primer (HAPfor), 5 ′ CAGTTCACGCCAAAGATGCC3 ′ ; reverse primer (HAPrev), 5 ′ CGCGTATGGTCCATCCTGAA3 ′ . Reverse primer containing three consecutive mismatched bases near its 3 ′ terminus can hybridize to the target RNA at a temperature (in this study: 42 • C) which is not high enough for DNA denaturation. Furthermore, due to the lower affinity between this mismatched primer and genomic DNA (gDNA), it is a very inefficient primer for gDNA amplification at high annealing temperatures (in this study: 62 • C). Consequently, the mismatched primer can be used to selectively amplify the newly synthesized complementary DNA (cDNA), even in the presence of high quantities of gDNA (Koo and Jaykus, 2000).

In vitro RT-PCR of Rhizobial Phytase Transcripts
HAP designed primer pair along with a primer pair that amplifies BBP described by Farhat et al. (2008) (BPPfor, 5 ′ GAT GCAGCTGATGATCCTGCG3 ′ ; BPPrev, 5 ′ ATTTTCTCCGTC CTGTGCGAC3 ′ ), were used to amplified B. elkanii SEMIA 587 phytases expressed in presence of phytate as unique P source. B. elkanii SEMIA 587 RNA and gDNA were extracted with the RNeasy Plus and DNeasy Plant Mini Kits (QIAGEN), respectively. RNA was extracted from B. elkanii SEMIA 587 cultivated in Angle medium supplemented with Na-phytate and KH 2 PO 4 .
A touchdown PCR was performed with HAP and BPP primers and with cDNA and gDNA as templates, to test if at different temperatures it was able to amplify a phytase. PCR reactions were carried out in 25 µL volumes containing 1 µL target, 10 X Buffer (Invitrogen), 50 mM MgCl 2 , 10 mM of each dNTP, 10 µM forward primer (HAP or BPP), 20 µM reverse primer (HAP or BPP), 0.1% bovine serum albumin (BSA), 5 U Taq Polymerase (Invitrogen), under the following thermal conditions: 95 • C for 2 min, 30 cycles of 95 • C for 30 s, 70 • C for 45 s, and gradually reduced 1 • C till reaching 60 • C, 72 • C for 45 s, with a final extension at 72 • C for 3 min.
When PCR products were obtained in the touchdown PCR, a gradient PCR was performed with cDNA and gDNA as templates, to obtain the optimal hybridization temperature of primers. PCR reactions were carried out in 25 µL volumes containing 1 µL target, 10 X Buffer (Invitrogen), 50 mM MgCl 2 , 10 mM of each dNTP, 10 µM forward primer, 20 µM reverse primer, 0.1% BSA, 5 U Taq Polymerase (Invitrogen), under the following thermal conditions: 94 • C for 3 min, 35 cycles of 94 • C for 30 s, 64-59 • C for 45 s, 72 • C for 45 s, with a final extension at 72 • C for 3 min. Consequently, the optimal hybridization temperature of primers HAP was 62 • C.
All amplification products were checked by electrophoresis in 1% agarose gel.

Hydroaeroponic Culture of Soybean
Two commercial soybean cultivars, namely N5909 (Nidera) and SR532 (Santa Rosa), were used in this study. SR532 and N5909 have a short and long production cycle, respectively. B. elkanii SEMIA 587 was used to inoculate these two cultivars. Before sowing, the soybean seeds were surface sterilized, first with 95% ethanol for 3 min and afterwards, with 3% NaOCl for 3 min more. Then, seeds were rinsed five times with sterile distilled water. Later, seeds were imbibed in sterile distilled water for 4 h, and then transferred for germination on sterile 15 g L −1 agarwater plates (Beyhaut et al., 2006). After germination, inoculation was performed by 30 min soaking 10-day-old seedlings in a suspension of B. elkanii SEMIA 587 containing ∼10 8 cfu mL −1 . The inoculum was prepared by growing the strain SEMIA 587 in yeast extract mannitol broth during four days at 28 • C at 120 rpm.
Twenty inoculated seedlings were transferred into each container, with base dimensions 0.2 × 0.4 m and height 0.2 m, refilled weekly with 40 L of Vadez et al. (1996) nutrient solution. In order to avoid initial N deficiency and ensure optimal nodulation, urea was supplied at 1 mmol plant −1 during the initial 15 days of growth. Thereafter, the plants were grown in N-free solution. They were grown under both deficient and sufficient P supply in hydroaeroponic conditions with a permanent flow of 400 mL min −1 of compressed air to ensure the oxygenation of the culture solution (Hernandez and Drevon, 1991). P was supplied weekly in the form of KH 2 PO 4 , at 75 vs. 250 µmol plant −1 as P deficiency vs. P sufficiency, respectively (Vadez et al., 1996). The pH was adjusted to 6.8 with 0.2 g L −1 CaCO 3 (Hernandez and Drevon, 1991). The experiment was carried out in a greenhouse, 28/20 • C under 16/8 h day/night cycle, with an additional light supply of 400 µmol photons m −2 s −1 and 70% relative humidity during the day.

Shoot, Root and Nodule Biomass
Forty two days after transfer (DAT), soybean plants were harvested and shoot, roots and nodules were separated. Dry weights (DW) were determined after drying at 70 • C for three days. Shoot accumulated P was determined by the vanadomolybdate method in Laboratorio Oriental (http://labo.com.uy/).
The efficiency in use of the rhizobial symbiosis (EURS) was estimated by the slope of the regression model of shoot biomass as a function of nodule biomass (y = ax + b), where a corresponds to the EURS and b corresponds to the plant biomass production without nodules (Zaman-Allah et al., 2007).

In situ RT-PCR of Rhizobial Phytase Transcripts
The in situ RT-PCR was made according to Lazali et al. (2013) and Maougal et al. (2014b). Three nodules of 3 mm diameter of three plants of each cultivar and P treatment were carefully detached from roots at 42 DAT. They were thoroughly washed with diethyl pyrocarbonate (DEPC) treated water, then fixed in 4% paraformaldehyde, 45% ethanol and 5% acetic acid, kept for 2 h under vacuum, and stored overnight at 4 • C. Fixed nodules were extensively rinsed with two baths of DEPC treated water during 5 min each, two baths of 1 X phosphate buffered saline (PBS, 5 mM Na 2 HPO 4 , 300 mM NaCl, pH 7.5) during 10 min each, one bath of 1 X PBS plus 0.2% glycine during 10 min and one last bath with 1 X PBS during 10 min.
Thereafter, the fixed nodules were dehydrated in baths of increasing concentrations of ethanol and butanol: one bath of 50% ethanol during 30 min; two baths of 70% ethanol during 30 min and 1 h; one bath of 70% of ethanol overnight at 4 • C; two baths with 100% ethanol for 1 h; and in the end, at least, three baths with 100% butanol during one week. Then, the nodules were included in paraffin and cut into transversal 12 µm thick sections using a microtome (Micro-cut H1200 Vibrating Microtome). The resulting sections were placed on slides treated FIGURE 1 | B. elkanii SEMIA 587 growth curves in Angle medium at pH 6 (A) or at pH 7 (B) supplemented with 1.5 mM KH 2 PO 4 (light gray) or with 6 mM Na-phytate (dark gray). Data represent average and standard deviation of three replicates. (C) RT-PCR or PCR products of B. elkanii SEMIA 587 grown in Angle medium supplemented with Na-phytate or with KH 2 PO 4 as P source in 1% agarose gel. Lane 1, GeneRuler 1 kb Plus DNA Ladder (Thermo Scientific); Lane 2, RT-PCR performed with HAP primers and cDNA as template obtained from SEMIA 587 cultivated with phytate as P source; Lane 3, RT-PCR performed with HAP primers (and BSA) and cDNA as template obtained from SEMIA 587 cultivated with phytate as P source; Lane 4, RT-PCR performed with BPP primers and cDNA as template obtained from SEMIA 587 cultivated with phytate as P source; Lane 5, PCR performed with HAP primers (and BSA) and gDNA as template obtained from SEMIA 587 cultivated with phytate as P source; Lane 6, PCR performed with BPP primers and gDNA as template obtained from SEMIA 587 cultivated with phytate as P source; Lane 7, RT-PCR performed with HAP primers (and BSA) and cDNA as template obtained from SEMIA 587 cultivated with KH 2 PO 4 as P source; Lane 8, RT-PCR performed with BPP primers and cDNA as template obtained from SEMIA 587 cultivated with KH 2 PO 4 as P source.
with silane to increase cell adhesion. Thereafter, the sections were deparaffinated through three baths of Safesolv (Q Path R ) for 10 min each, one bath of Safesolv during 5 min; three baths of 100% ethanol during 10 min each; one bath of 70% ethanol during 5 min; one bath of 50% ethanol during 5 min. This was continued by section rehydration with one bath of 1 X PBS during 10 min, one bath of 1 X PBS plus glycine 0.2% during 2 min, one bath of 1 X PBS during 10 min and one bath of 1X PBS for 5 min at 65 • C.
The first cDNA strand was synthesized from total RNA of nodular sections, that were incubated at 42 • C during 1 h in 100 µL reverse transcriptase mix containing 5 X Buffer (Promega), 10 mM of each dNTP, 10 µM HAPrev, 0.1% BSA, and 200 U M-MLV reverse transcriptase H-(Promega).
Negative controls (NRT) were prepared by omitting the reverse transcriptase. Afterwards, the reverse transcriptase mix was removed, by adding two times 1 mL of 1 X PBS and incubating the slides two times more in 50 mL of 1 X PBS during 10 min. PCR reactions were carried out in 100 µL volumes containing 10 X Buffer (Invitrogen), 50 mM MgCl 2 , 10 mM of each dNTP, 1 mM digoxigenin-11-uridine-triphosphate (Dig-11-dUTP), 10 µM HAPfor, 20 µM HAPrev, and 5 U Taq polymerase (Invitrogen), under the following thermal conditions: 35 cycles of 94 • C for 30 s, 62 • C for 45 s, and 72 • C for 45 s, with a final extension at 72 • C for 3 min.
For detection of the amplified cDNA in the nodules, the PCR mix was removed, and samples were washed in 100 µl of blocking solution under gentle agitation in the dark FIGURE 2 | (A) Shoot dry weight, (B) nodules dry weight, and (C) root dry weight of G. max cultivars N5909 and SR532 inoculated with B. elkanii SEMIA 587 under P sufficiency (250 µmol P plant −1 week −1 ) (white) vs. P deficiency (75 µmol P plant −1 week −1 ) (gray). Data represent average and standard deviation of ten replicates harvested at 42 DAT. Treatments with a common letter are not significantly different (p > 0.05). at 37 • C. The blocking solution was replaced by 100 µl of alkaline phosphatase conjugated anti-dioxygenin Fab fragment (Roche Diagnostics) diluted 1:1000 in 2% BSA. The samples were incubated at room temperature for 90 min, then washed three times to remove excess antibody. Detection of alkaline phosphatase was carried out using an ELF-97 R endogenous phosphatase detection kit (Molecular Probes). The ELF substrate was diluted 1:40 in the alkaline detection buffer (Molecular Probes), vigorously shaken, and then filtered through a 0.22 µm filter (Millex R -GV, Millipore) to remove any aggregates of the substrate that may have formed during storage. Samples were incubated in 20 µL ELF substrate-buffer solution for 20 min in the dark and transferred to washing buffer (1 X PBS, 25 mM EDTA, and 5 mM levamisole, pH 8.0). Three washings of 1 min were performed. Samples were mounted for observations with a BX61 R microscope (Olympus) equipped with an epifluorescence condenser, a Hoechst/DAPI filter set configured with an excitation filter of 360-370 nm, a dichroic mirror of 400 nm and an emission of 420 nm, and a gray View II R camera (ORCA AG). Image analysis was performed using ImageJ software (Schneider et al., 2012) as an image analysis program. The intensity of the fluorescent signal emitted by the rhizobial phytase transcript was measured as number of green pixels per image.

Statistical Analysis
Shoot, nodules and root DW, shoot accumulated P and intensity of the fluorescent signal emitted by the rhizobial phytase transcript were analyzed through analysis of variance (ANOVA) and tested for differences between cultivars and P growing conditions using a post-hoc Tukey's HSD test (p < 0.05) in R3.1.3 software (https://www.r-project.org/) using the "agricolae" package (de Mendiburu, 2020). The relationship between nodule and shoot biomass was tested by regression analysis. FIGURE 4 | Shoot accumulated phosphorus in G. max cultivars N5909 and SR532 inoculated with B. elkanii SEMIA 587 under P sufficiency (250 µmol P plant −1 week −1 ) (white) vs. P deficiency (75 µmol P plant −1 week −1 ) (gray). Data represent average and standard deviation of ten replicates harvested at 42 DAT. Treatments with a common letter are not significantly different (p > 0.05).

RESULTS
Phytate Mineralization by B. elkanii SEMIA 587 B. elkanii SEMIA 587 showed similar growth curves when it was cultivated in Angle medium supplemented with Na-phytate or KH 2 PO 4 , as P sources, at both pH evaluated (6 and 7) (Figures 1A,B). Therefore, this strain was able to mineralize phytate and to use it as a P source.
Two primers pairs, that amplify two types of phytases, HAP and BPP, were tested through a retro-transcription and a subsequent touchdown PCR in B. elkanii SEMIA 587. An amplification product was obtained when SEMIA 587 was cultivated with Na-phytate as P source and when HAP primers were used, but not when BPP primers were used ( Figure 1C). When SEMIA 587 was cultivated with KH 2 PO 4 as P source, no amplification products were obtained either when HAP or BPP primers were used ( Figure 1C). Additionally, we demonstrated that HAP gene was not amplified from gDNA with these HAP primers (Figure 1C), due to the lower affinity between the mismatched primer and the gDNA at the used annealing temperature (62 • C).

Efficiency in Use of the Rhizobial Symbiosis for Plant Growth
Cultivar N5909 cultivated under P sufficiency showed more shoot, nodules and root DW than under P deficiency (p = 0.0001, p = 0.0007, p = 0.0300, respectively), whereas no significant differences of shoot, nodules and root DW under P sufficiency and P deficiency in cultivar SR532 were found (Figure 2).
The EURS was estimated as the regression slope of shoot biomass as a function of nodule biomass, whenever the correlation between both parameters was significant (up to R 2 = 0.719). Under P sufficiency, SR532 was the most efficient, and showed similar EURS under P deficiency, whereas N5909 increased its EURS under P deficiency (Figure 3).
Cultivar N5909 accumulated more P in its shoots under P sufficiency than under P deficiency (p = 0.0008), whereas cultivar SR532 did not show any significant difference in shoot accumulated P between P treatments (Figure 4).

Rhizobial Phytase Transcript Abundance in Soybean Nodules
Using the In situ RT-PCR technique we were able to show the transcript abundance of the rhizobial phytase in nodular sections of soybean cultivars N5909 and SR532 inoculated with B. elkanii SEMIA 587 cultivated under P sufficiency and P deficiency (Figures 5A-D). The NRT nodules not subjected to reverse transcription did not display any fluorescent signal (Figures 5E,F), confirming that fluorescence was not due to methodological artifacts. The rhizobial phytase was transcripted in all bacteroidal symbiosomes in the infected zone, with a high transcript abundance rate since the fluorescent signal, as assessment of number of transcripts, was intense for all the treatments (Figures 5A-D).
In the case of cultivar SR532, the fluorescent signal was more intense under P deficiency in comparison to P sufficiency (p < 0.0001), whereas in the case of cultivar N5909, there were no significant differences between the intensity of the fluorescent signal under P deficiency or sufficiency. The intensity of the fluorescent signal in SR532 under P deficiency was equal to the intensity of this signal in N5909 under both P conditions (Figure 6). The ANOVA showed a significant interaction between the factors P supply and cultivar for the variable intensity of fluorescent signal (p = 0.0010).

DISCUSSION
The present study shows for the first time that a rhizobial strain could mineralize phytate and use it as a P source and that we could locate and observe the distribution of the transcripts of a rhizobial phytase within nodules of N 2 -fixing soybean plants, thanks to the in situ RT-PCR methodology carried out. This method has proven suitable to quantify the level of transcripts in prokaryote (Maougal et al., 2014b) and in eukaryotes cells (van Aarle et al., 2007;Bargaz et al., 2012;Lazali et al., 2013), and in the case of the present study, we were able to quantify the expression of rhizobial phytase transcripts with this method. The rhizobial phytase transcripts were localized only in the bacteroidal cells. It has been shown that P in G. max nodules was higher in the infected zone in comparison to the cortical region (revised by Hu et al., 2018). Thus, the high level of transcripts in the infection zone may respond to bacteroidal requirement for metabolism and survival, as well as bacteria multiplication (Bargaz et al., 2012;Lazali et al., 2013).
Legumes nodules are plant tissues with high concentration of P and when plants are under P-depletion, they follow a strategy to maintain SNF and viable leaf tissue as long as possible, expressing several P cycling genes in nodules (Cabeza et al., 2014). Plant phosphatase transcripts have been previously shown as affected by plant available P (Araújo et al., 2008;Bargaz et al., 2012;Lazali et al., 2013). Previously, the transcript abundance and activity of the PAP phytase of common beans was shown to increase when plants were grown under P deficiency, independently of the plant FIGURE 6 | Intensity of the fluorescent signal emitted by the rhizobial phytase transcripts in nodules sections of G. max cultivars N5909 and SR532 inoculated with B. elkanii SEMIA 587 under P sufficiency (250 µmol P plant −1 week −1 ) (white) vs. P deficiency (75 µmol P plant −1 week −1 ) (gray). Data represent average and standard deviation of three nodules of three plants. The y-axis intensity of fluorescence is expressed in number of green pixels and must by 10 8 . Treatments with a common letter are not significantly different (p > 0.05).
genotype (Araújo et al., 2008;Lazali et al., 2013). In addition, Bargaz et al. (2012) demonstrated that the transcript abundance and activity of the plant phosphoenol pyruvate phosphatase (PEPase) increased under P deficiency for several common bean genotypes. Contrary, the two soybean cultivars performed contrastingly when grown under low-and high-P supply. While cultivar N5909 did not show differences in intensity of the fluorescent signal under low-and high-P, cultivar SR532 showed a high transcript abundance of the rhizobial phytase in the nodules when grown under P deficiency. Moreover, the transcript abundance of the rhizobial phytase in the nodules of SR532 under P sufficiency seemed to be repressed compared to the transcript abundance in nodules of N5909 under the same conditions. These facts suggest that the transcript abundance rate of the rhizobial phytase depends on the plant genotype and P level, although further testing of other soybean genotypes is required. Additionally, the high transcript abundance level of the rhizobial phytase in the nodules of the cultivar SR532 under P deficiency, suggests that this enzyme contributes to the adaptation of this cultivar to P deficiency. Cultivar SR532 under P deficiency shows a correlation of high bacterial phytase transcripts with high EURS, suggesting a high regulation between EURS and the P content. Overall, our results suggest that the high tolerance of SR532 to P deficiency, as compared to N5909, is associated with better capacity to maintain SNF under low P supply based on the ability to increase the rhizobial phytase transcript abundance under such condition. Phytases might play a major role for internal plant metabolism rather than for obtaining P from the soil phytate (Tang et al., 2006). Phytate seems to act as a phosphate storage (Raboy, 2003) and rhizobial phytases can mineralize the phytate present in the soybean nodules, increasing the concentration of orthophosphates ions necessary for the SNF (Lazali et al., 2013).
Although further studies are required to identify new rhizobial strains with phytase activity and to test their interaction with various soybean genotypes, the differential transcript abundance shown in this research suggests that this phenomenon can contribute to the external and internal use efficiency of phytate for SNF. Sulieman et al. (2013b) found that different strains of S. meliloti showed different symbiotic efficiency with M. truncatula growing under different P levels. Thus, the identification and selection of an efficient combination of rhizobia and soybean genotype might play an important role in soybean breeding and production, aiming to reduce synthetic P fertilizer inputs and to improve sustainability in soils with low available P. We propose to consider the ability of a strain to mineralize phytate at the moment of evaluate new N 2 fixing rhizobial strains, and to evaluate how the plant genotype influences this process, looking for combinations of rhizobia and plant that enhances P metabolism, thus SNF, under low available P soils, to select the best adapted symbiosis to P deficient soils.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
VC performed the experimental work, analyzed data, and wrote the manuscript. LA and CT coordinated the experimental work. NA, EB, and J-JD conceived, planned the study, and helped to write the manuscript. All authors read and approved the final manuscript.