Next Article in Journal
Bioinformatics Analysis and Functional Prediction of Transmembrane Proteins in Entamoeba histolytica
Previous Article in Journal
Genome-Wide Identification and Transcriptional Regulation of Aquaporin Genes in Bread Wheat (Triticum aestivum L.) under Water Stress
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 9
Issue 10
10.3390/genes9100498
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phosphate Deficiency Negatively Affects Early Steps of the Symbiosis between Common Bean and Rhizobia

by
Mariel C. Isidra-Arellano
1,2,†,
María Del Rocio Reyero-Saavedra
1,†,
Maria Del Socorro Sánchez-Correa
1,
Lise Pingault
3,
Sidharth Sen
4,
Trupti Joshi
4,5,6,
Lourdes Girard
7,
Norma A. Castro-Guerrero
8,
David G. Mendoza-Cozatl
8,
Marc Libault
3 and
Oswaldo Valdés-López
1,*
1
Laboratorio de Genómica Funcional de Leguminosas, Facultad de Estudios Superiores Iztacala, Universidad Nacional Autónoma de Mexico, Tlalnepantla 54090, Estado de Mexico, Mexico
2
Posgrado en Ciencias Biológicas, Universidad Nacional Autónoma de Mexico, Coyoacan 04510, Ciudad de Mexico, Mexico
3
Department of Agronomy and Horticulture, University of Nebraska-Lincoln, Beadle Center, Lincoln, NE 68503, USA
4
Informatics Institute, University of Missouri, Columbia, MO 65211, USA
5
Christopher S. Bond Life Sciences Center, University of Missouri, Columbia, MO 65211, USA
6
Department of Health Management and Informatics, School of Medicine, University of Missouri, Columbia, MO 65211, USA
7
Departamento de Biología de Sistemas y Biología Sintética, Centro de Ciencias Genómicas, Universidad Nacional Autónoma de Mexico, Cuernavaca 62210, Morelos, Mexico
8
Division of Plant Sciences, C. S. Bond Life Sciences Center, University of Missouri, Columbia, MO 65211, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2018, 9(10), 498; https://doi.org/10.3390/genes9100498
Submission received: 11 September 2018 / Revised: 5 October 2018 / Accepted: 10 October 2018 / Published: 15 October 2018
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Phosphate (Pi) deficiency reduces nodule formation and development in different legume species including common bean. Despite significant progress in the understanding of the genetic responses underlying the adaptation of nodules to Pi deficiency, it is still unclear whether this nutritional deficiency interferes with the molecular dialogue between legumes and rhizobia. If so, what part of the molecular dialogue is impaired? In this study, we provide evidence demonstrating that Pi deficiency negatively affects critical early molecular and physiological responses that are required for a successful symbiosis between common bean and rhizobia. We demonstrated that the infection thread formation and the expression of PvNSP2, PvNIN, and PvFLOT2, which are genes controlling the nodulation process were significantly reduced in Pi-deficient common bean seedlings. In addition, whole-genome transcriptional analysis revealed that the expression of hormones-related genes is compromised in Pi-deficient seedlings inoculated with rhizobia. Moreover, we showed that regardless of the presence or absence of rhizobia, the expression of PvRIC1 and PvRIC2, two genes participating in the autoregulation of nodule numbers, was higher in Pi-deficient seedlings compared to control seedlings. The data presented in this study provides a mechanistic model to better understand how Pi deficiency impacts the early steps of the symbiosis between common bean and rhizobia.

1. Introduction

Phaseolus vulgaris (common bean) represents the principal source of non-animal proteins and minerals for human consumption in the developing world [1]. Additionally, common bean also provides a variety of secondary metabolites with analgesic and neuroprotective properties [2,3]. Similar to other legumes, common bean is able to symbiotically interacts with atmospheric nitrogen-fixing soil bacteria (rhizobia) and incorporates the fixed-nitrogen (e.g., ammonium) into the food chain [1]. Due to its agronomical, medical and ecological properties, common bean is one of the most important legumes worldwide.
To establish a successful symbiosis with rhizobia, a molecular dialogue between legumes and these soil microbes is required. This molecular dialogue begins with the detection of legume-derived flavonoids by compatible rhizobia [4,5]. In response, rhizobia produce diffusible nodulation factors (NF) which are perceived by legumes via the LysM-domain receptor kinase Nod Factors Receptor 5 (NFR5) and NFR1. These receptors are both located at the root hair plasma membrane similarly to the epidermal LysM receptor eNFR [6,7]. Upon NF perception, a wide variety of molecular responses occurs, which is critical to activate the cellular responses required for rhizobial infection and the formation of the new organ termed nodule [8]. For instance, some of these molecular responses are required to modify the polar growth of the root hair leading to the formation of a curled tip [9,10]. In the center of this curl, an infection chamber is formed where rhizobia are entrapped and multiply to form a micro-colony [11]. Immediately, an infection thread is initiated from the infection chamber and elongates through the nodule primordium located in the root cortex [11,12]. Rhizobia are loaded into the infection thread and subsequently released into the nodule primordium cells where they are differentiated into bacteroids. Then, through the activity of the nitrogenase, the bacteroids fix and assimilate the atmospheric dinitrogen within the plant nodule [13].
Root infection and colonization by rhizobia and nodule formation are controlled by two independent genetic programs, which are tightly coordinated by a core of symbiosis-related genes, which belong to the common symbiosis pathway [8,14]. Among these core genes, the leucine-rich repeat (LRR) receptor-like kinase SYMRK along with the plasma-membrane receptor NFR5 and NFR1 are required to detect and decode the NF signal [15,16]. Downstream to NF perception, different ion channels (e.g., the potassium channel Castor/Pollux and the calcium channels belonging to the CNGC15s family) and nucleoporins (e.g., NENA, NUP83 and NUP85) are activated to initiate and maintain rapid oscillations of the nuclear and perinuclear calcium concentrations (i.e., calcium spiking) [17,18,19,20,21]. These fluxes of calcium are further transduced by a calcium and calmodulin-dependent protein kinase (CCaMK), which, in turn, phosphorylates the transcription factor CYCLOPS [22,23]. CYCLOPS is required to activate other critical transcription factors, including the Nodulation Signaling Pathway 2 (NSP2), Nodule Inception (NIN), Nuclear-Factor YA (NF-YA) and NF-YB as well as Ethylene Responsive Element1 (ERN1) [23,24,25]. The coordinated action of these transcription factors is essential to activate the expression of the genes required for the infection of root hairs by rhizobia [8,14]. Among them, Flotillin2 (FLOT2) is crucial to the formation of the infection thread [26]. In parallel, a second genetic program triggering the formation of the nodule meristem in the roots is activated [8,13]. It has been demonstrated that a delicate balance among the phytohormones auxins, cytokinins, and ethylene is required to trigger formation of the nodule meristem [10,27,28]. Genetic analyses have revealed that mutations in the receptors or transporters of these phytohormones affect nodule development. For example, mutations in the cytokinins receptor MtCRE1/LjLHK1 or in the central regulator of the ethylene signaling EIN2 lead to spontaneous nodulation or over-formation of nodules in the legume models Lotus japonicus and Medicago truncatula, respectively [29,30].
Because symbiotic nitrogen fixation (SNF) is a high-energy demanding process, legumes tightly regulate the number of infections and nodules per plant [31]. This autoregulation of the nodulation (AON) process is locally and systemically activated by the legume. At the systemic level, the plant-derived CLAVATA3/embryo-surrounding region (CLE) peptides are considered as critical regulators of the AON process [31]. In common bean, two AON CLE peptide-encoding genes have been identified, Rhizobia-Induced CLE1 (PvRIC1) and Rhizobia-Induced CLE2 (PvRIC2) [32]. PvRIC1 expression is induced as early as one-day post rhizobia inoculation, whereas PvRIC2 expression is approximately induced five days upon rhizobia inoculation [32]. PvRIC1 and PvRIC2 peptides are synthesized in the root and translocated to the leaves via the xylem. Upon perception of PvRIC1 and PvRIC2 peptides by a Leucine-rich Repeat Receptor-like Kinase (LRR-RLK) called PvNARK, the expression of different genes participating in the negative regulation of the meristematic activity of the nodule primordium cell is activated [31,32]. It has been widely demonstrated that mutations in the PvNARK/GmNARK/LjHAR1/MtSUNN receptor lead to a hypernodulation phenotype along with a reduction in plant growth [32,33,34,35,36].
Most of the arable soils where common bean is produced have low levels of phosphorus (P), an essential macronutrient required for optimal plant growth and development. Common bean, like other plants, exclusively takes up P from the soil as inorganic phosphorus (i.e., orthophosphate, Pi) [37,38]. However, soil Pi rapidly interacts with Fe3+, Ca2+, and Al+3 ions, forming complexes making Pi unavailable for plant uptake [39]. Pi deficiency constrains common bean growth and productivity as well as its ability to form nodules and symbiotically fix atmospheric nitrogen [40]. It has been reported that Pi deficiency decreases common bean nitrogen-fixation efficiency by increasing the permeability of the nodule to oxygen leading to the inactivation of the rhizobia-biosynthesized nitrogenase [41,42]. Furthermore, it has also been shown that Pi deficiency reduces by more than 60% the number of nodules in common bean plants [40].
Different efforts to understand how legumes cope with Pi deficiency while symbiotically interacting with rhizobia have been made [43,44,45,46]. Physiological studies on different legumes have demonstrated that under Pi deficiency conditions, Pi is preferentially relocated from other organs to the nodules [47,48,49,50,51,52]. Furthermore, both transcriptional and metabolic studies in common bean, M. truncatula, and chickpea have revealed a more efficient use of internal Pi sources and the maintenance of Pi homeostasis under limiting conditions [40,43,46]. For instance, it has been reported that the concentrations of different organic acids (e.g., citrate, succinate and pyruvate) were differentially regulated in nodules of both common bean and chickpea plants [40,43,46]. Among them, many are participating in the cation exchange that is required for the Pi uptake. Additionally, it has also been reported that the expression of flavonoid-synthesis related genes decreased in Pi-deficient M. truncatula roots [43].
Nodule efficiency and the molecular responses underlying the adaptation of the nodule to Pi deficiency have been studied. However, the molecular mechanisms leading to a reduction in nodule number under Pi deficiency remain unknown. For instance, it is still unclear whether Pi deficiency interferes with the molecular dialogue between common bean and rhizobia. If so, what part of this molecular dialogue is impaired?
In this study, we report the early molecular and physiological responses triggered by rhizobia in common bean seedlings growing under Pi deficient conditions. A transcriptional analysis by quantitative real-time polymerase chain reaction (qRT-PCR) revealed that the expression of PvNSP2, PvNIN and PvFLOT2, was quickly reduced by 50% in Pi-deficient common bean seedlings. Although curled root hairs were observed regardless of the Pi condition (optimal or deficient conditions), their number was severely reduced in Pi-deficient seedlings. Interestingly, the majority of the rhizobia-induced curled root hairs in Pi-deficient seedlings showed aberrant deformations that have previously been associated with defects in the rhizobia infection process. Additionally, we also observed that the expression of the PvRIC1 and PvRIC2 genes, regardless of the mock or rhizobia treatment, significantly increased in Pi-deficient seedlings. To gain more insights about the transcriptional modifications leading to an impaired communication between common bean and rhizobia under Pi-deficient conditions, we performed a genome-wide transcriptomics analysis. This analysis revealed that the expression of genes associated with early steps of this symbiosis (e.g., signal transduction and hormone balance) was affected in Pi-deficient plants. The data presented in this study pave the way to better understand the genetic mechanisms underlying the establishment of the symbiosis between rhizobia and common bean under Pi-deficient conditions.

2. Material and Methods

2.1. Plant Material

Common bean (P. vulgaris L. cv Negro Jamapa) seeds were produced under greenhouse conditions in the Facultad de Estudios Superiores Iztacala, Universidad Nacional Autónoma de México (UNAM), at Tlalnepantla, Estado de Mexico, Mexico. Seeds were surfaced sterilized by soaking in 70% ethanol for 1 min, followed by treatment for 10 min with 10% bleach. Seeds were subsequently washed 10 times in sterile water. Sterilized seeds were germinated for three days at 25 °C in Petri dishes containing wet and sterile germination paper (Anchor paper CO, St. Paul, MN, USA). Petri dishes containing common bean seeds were kept under dark conditions at 25 °C.

2.2. Bacterial Strains and Culture Conditions

Rhizobium tropici CIAT 899 strain was used to inoculate common bean seedlings. R. tropici cells were grown for 2 days at 30 °C on PY medium (5 mg/L peptone; 3 mg/L yeast extract) supplemented with 0.7 M CaCl2 and 20 μg/mL nalidixic acid (Nal). After two days, R. tropici cells were harvested and resuspended in sterile water at O.D600 = 0.3. 1 mL of this bacterial suspension was used to individually inoculate common bean seedlings.
Escherichia coli DH5α strain carrying the plasmid pFAJ1700::pLacZ::uidA was used as a plasmid donor for triparental mating with R. tropici. E. coli cells were grown in Luria-Bertani (LB) broth supplemented with 10 μg/mL tetracyclin (Tc) overnight at 37 °C.
E. coli pRK2013 strain was used as a helper for triparental mating with R. tropici. pRK2013 cells were grown in LB broth supplemented with 30 μg/mL kanamycin overnight at 37 °C.
R. tropici::uidA strain was obtained by triparental mating using the E. coli DH5α strain carrying the plasmid pFAJ1700::placZ:uidA, containing the uidA gene under the lacZ promoter as a donor, and the pRK2013 as conjugation helper. R. tropici derivatives carrying the plasmid were selected in PY medium as Nal and Tc resistant transjugants. R. tropici::uidA cells were grown in PY medium supplemented with 0.7 M CaCl2, 10 μg/mL Tc, and 20 μg/mL Nal for 2 days at 30 °C. After two days, cells were harvested and resuspended in sterile water at OD600nm = 0.3. 1 mL of this bacterial suspension was used to individually inoculate common bean seedlings.

2.3. Treatment

Three-day-old common bean seedlings were transferred into a nitrogen-free Fähraeus medium [53] plates supplemented with 1 mM PO4 (optimal-Pi conditions) or 5 μM PO4 (low-Pi conditions). Three days after transplanting, seedlings were inoculated with autoclaved-water (Mock) or with R. tropici CIAT899. Inoculated seedlings were kept under dark conditions at room temperature. At 1, 24 and 48 h post-inoculation, the zone of the roots susceptible to rhizobium infection (hereafter referred as susceptible zone) was isolated, flash frozen in liquid nitrogen, and stored at −80 °C until used for total RNA extraction. The susceptible zone encompassed a region of ~1 cm and is analogous to the distal region of the elongation zone and the entire differentiated zone of the root [54]. Four biological replicates, each one containing susceptible zones from four different seedlings, were independently generated.

2.4. Quantification of Soluble Phosphate

Soluble Pi content was determined in roots from common bean seedlings grown for three days in optimal- or low-Pi conditions using the colorimetric assay reported in [55]. Briefly, roots were harvested, rinsed with ultra-pure water, weighed, immediately homogenized in 10 N trichloroacetic acid and fractioned by centrifugation. A soluble fraction was incubated with the staining solution (350 mM FeSO4, 16% (NH4)6Mo7O24) for 10 min at room temperature. Optical density was determined at 660 nm. For each experimental condition, five replicates, each one containing four roots from different seedlings, were analyzed.

2.5. Root Hair Deformation Analysis

Three-day-old common bean seedlings grown for three days under optimal-Pi or low-Pi conditions were inoculated with 1 mL of R. tropici suspension (OD600nm = 0.3). Forty-eight hours after inoculation, roots were collected and stained with methylene blue to maximize contrast, and then observed with a bright-field microscope. A total of 20 independent biological replicates were generated, each one including 10 plants.

2.6. Quantification of Number of Infection Threads

To quantify the number of infection threads per root, common bean seedlings grown for three days under optimal-Pi or low-Pi conditions were inoculated with 1 mL of R. tropici::uidA suspension (OD600nm = 0.3). Three days post-inoculation, susceptible zones were collected, rinsed twice with autoclaved ultra-pure water and immersed in uidA staining solution (0.05% 5-bromo-4-chloro-3-indolxyl-β-d-glucuronic acid, 100 mM sodium phosphate buffer (pH7), 0.5 mM potassium ferrocyanide, 0.5 mM potassium ferricyanide, 10 mM Na2EDTA and 0.1% Triton X-100) and incubated for one hour at 37 °C. After rinsing in phosphate buffer, the number of infection threads formed at the susceptible zone was quantified under a bright-field microscopy. Additionally, to evaluate the rhizobia adhesion to the common bean roots, the intensity of the blue precipitate (as a product of the activity of the R. tropici::uidA) formed in the susceptible zone was densitrometrically quantified. To this end, pictures of susceptible zones showing blue precipitate were captured using a SZX10 stereomicroscope (Olympus, Center Valley, PA, USA) equipped with an Olympus UC50 camera (Olympus). The intensity of the blue precipitate in each susceptible zone was quantified as number of pixels by using the software CellSens Standard (Olympus). For these experiments, four biological replicates, each one with 10 roots from different seedlings, were included.

2.7. Gene Expression Analysis

To analyze the expression of PvNSP2, PvNIN, PvFLOT2, PvRIC1 and PvRIC2 genes, total RNAs were extracted from 0.5 g of rhizobia-inoculated or mock-inoculated susceptible zones from common bean seedlings growing under optimal-Pi or low-Pi conditions using ZR Plant RNA Miniprep kit (Zymo Research, Irvine, CA, USA) according to manufacturer’s instructions. Genomic DNA (gDNA) was removed from purified RNA using DNaseI RNase-free (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. 1 μg of gDNA-free total RNA was used to synthesize cDNA using Thermo Scientific RevertAid Reverse Transcriptase (Thermo Fisher Scientific, USA) according to manufacturer’s instructions. cDNA samples were used to analyze the expression of the aforementioned genes by qRT-PCR in a Step-One qPCR thermocycler (Applied Biosystems, Foster City, CA, USA). The housekeeping gene PvActin was used to normalize gene expression levels [56]. The expression level of different genes was calculated according to the equation E = nPeff(−ΔCt) where Peff is the primer set efficiency calculated using LinRegPCR [57] and ΔCt is the difference of the cycle threshold (Ct) values of the housekeeping gene and a given gene. qRT-PCR primers were designed by using the online available software Primer 3 v.0.4.0 [58] and using the following parameters: 80–120 pb as product size, 19–23 nt as primer size, 60 °C as primer temperature melting, and 40–60% of primer GC content. To design primers, the sequence from the 3’UTR region of each gene was obtained from the P. vulgaris genome Version 2.1 [59,60]. The nucleotide sequences of the qRT-PCR primers used in this study are provided in Table S1. For this experiment, four biological replicates, each one containing roots from four different seedlings, were included.

2.8. Preparation of Messenger RNA-Seq Library and High-Throughput Sequencing

Total RNA was isolated from 0.5 g of 48 h rhizobia-inoculated or mock-inoculated roots from common bean plants growing under optimal-Pi or low-Pi conditions, as described earlier in this paper. Stranded messenger RNA-seq (mRNA-seq) libraries were generated from 1.5 μg of gDNA-free total RNA from each experimental condition and prepared using TruSeq RNA Sample Prep kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Briefly, sample concentration was determined by Qubit flourometer (Invitrogen, Carlsbad, CA, USA) using the Qubit HS RNA assay kit, and the RNA integrity was checked using the Fragment Analyzer automated electrophoresis system (Advanced Analytical, Ankeny, IA, USA). Poly-A containing mRNAs were purified from total RNAs and subsequently fragmented. cDNAs were generated from fragmented RNAs, and the index-containing adapters were ligated to their ends. The amplified cDNAs were purified by the addition of Axyprep Mag PCR Clean-up beads (Axygen, Corning, NY, USA). The quality of each purified RNA-seq library was evaluated using the Fragment Analyzer automated electrophoresis system, and quantified with the Qubit flourometer using the Qubit HS dsDNA assay kit.
In total, 12-stranded mRNA-seq [two Pi conditions (sufficient and deficient), two treatments (mock or rhizobia), and three independent biological replicates] were high-throughput sequenced in 75 bp single-end at the University of Missouri DNA Core Facility by using an Illumina NextSeq500 sequencer (Illumina, USA) according to the manufacturer’s instructions. Image analysis and base calling were performed using the Illumina pipeline (http://www.illumina.com).

2.9. Mapping and Processing mRNA-Seq Reads

mRNA-seq analysis was conducted using a custom in-house developed informatics pipeline which first performs a quality check on the raw sequencing FASTQ files using the tool FastQC [61] and trim-galore tool [62] to remove low-quality reads and trim adaptor sequences. We then indexed the P. vulgaris Version 2.1 [59,60] reference genome using Bowtie2 short read aligner [63]. The trimmed mRNAseq reads were then aligned, allowing two mismatches, to this index reference genome using TopHat, which also reports splice junction sites [64].

2.10. Quantification and Identification of Differentially Regulated Genes

The next step in our pipeline was to quantify the expression level of each transcript/gene by measuring the level of mRNA-seq alignments using the tool Cufflinks [65]. The quantified expression level was represented by Fragments per Kilobase of transcript per Million mapped reads (FPKM) values. The differential expression levels of quantified genes were then calculated using the tool Cuffdiff [65] which takes the observed log fold change of a gene’s expression from control and treatment samples and then reports if the gene is significantly differentially expressed (q-values). The thresholds for Cuffdiff were set at 2-fold change and q-value ≤ 0.05 before calling genes to be differentially expressed, which are usually annotated with their unique identifiers defined by the gene models annotation and gene coordinates files. Data obtained from this differentially expression analysis was independently validated by qRT-PCR.

2.11. Gene Functional Classification

The biological relevance of the differentially regulated genes was assessed by a gene function enrichment analysis using the method Singular Enrichment Analysis (SEA) available in the web-based tool AgriGo [66,67], followed by a summarization and visualization of statistically significant non-redundant gene ontology (GO) terms by using the web-based tool REVIGO [68,69]. Briefly, GO terms enriched in each set of genes were compared to the P. vulgaris Version 2.1 gene reference background [59,60]. p-values for the GO terms were obtained through Fisher’s exact test, and a q-value was computed to produce lists of significant GO terms with an estimated False discovery rate (FDR) of 5%. Enriched GO terms with q < 0.05 were further summarized and visualized on REVIGO. Additionally, MapMan gene functional classification was used [70,71]. For MapMan analysis, the common bean mapping file, available at MapMan website [72], was used.

2.12. Statistical Analyses

All the statistical analyses were conducted using R software 3.0.1. The specific tests performed are indicated in the legend of the corresponding figure.

3. Results and Discussion

3.1. Rhizobia-Induced Root Hair Deformation and Infection Thread Formation is Affected in Pi-Deficient Common Bean Seedlings

It has been demonstrated that the reduction in nodule number is one of the negative effects of Pi deficiency in legumes [40]. However, it is unclear whether this nutritional stress interferes with the molecular dialogue between the common bean and rhizobia; therefore, the precise stages of this molecular dialogue that are potentially affected by Pi deficiency remain unknown. To better characterize the effect of the Pi deficiency in the molecular dialogue between the common bean and rhizobia, we first developed a system which allows us to track early molecular and physiological responses to rhizobia under Pi limiting conditions. This system allows the growth of common bean seedlings under optimal- (1 mM) or low- (5 μM) Pi conditions, as well as under sterile and controlled environmental conditions (Figure S1). By using this system, we were able to obtain Pi-deficient common bean seedlings, which led to a 50% reduction in plant Pi contents (Figure S1). Once we confirmed that our experimental system yields Pi-deficient common bean seedlings, we proceeded to assess early responses to rhizobia in common bean seedlings growing under optimal- or low-Pi conditions.
One of the early physiological responses triggered by rhizobia is the root hair deformation, which is required for rhizobia entrapping and the formation of the infection chamber [9,10]. To test whether Pi deficiency affects the rhizobia-induced root hair deformation and infection thread formation, common bean seedlings growing for three days under optimal- or low-Pi conditions were inoculated with R. tropici CIAT899::uidA. Forty-eight hours post inoculation, 95% (n = 190) of the seedlings growing under optimal-Pi conditions showed rhizobia-induced root hair deformation (Figure 1a,c). In contrast, only 50% (n = 200) of the Pi-deficient seedlings displayed deformed root hairs (Figure 1b,c). Interestingly, most of the deformations observed in Pi-deficient common bean seedlings were aberrant (Figure 1d). For instance, we observed root hairs with swelled tip, root hairs with more than one tip outgrowth, and root hairs with spatula-like deformation (Figure 1e). Additionally, we observed that Pi deficiency significantly affected rhizobia attachment, which was reflected by a reduction in the blue precipitate (as a product of the activity of the R. tropici::uidA strain) at the susceptible zone (Figure 2a,b). Accordingly, we observed that Pi-deficient common bean seedlings developed 60% less infection thread than the seedlings growing under optimal-Pi condition and inoculated with rhizobia (Figure 2c–e). Similar aberrant rhizobia-induced root hair deformations have been observed in M. truncatula ern1/ern2, M. truncatula dmi1 and in L. japonicus scarn mutant plants [9,53,73,74]. This defect has been associated with the inhibition in the polar growth of the root hair cells, which affects the rhizobia infection process including the infection thread formation [9,73,74]. Additionally, an early study reported that Pi deficiency affects the rhizobia attachment to Medicago roots [75]. Similarly, our data suggest that Pi deficiency negatively affects the polar elongation of the root hairs, likely leading to inefficient attachment and entrapping of rhizobia.

3.2. Pi Deficiency Decreases the Expression of Genes Required for Root Hair Curling and Rhizobia Infection

The coordinated action of different transcription factors is essential for a successful symbiosis with rhizobia. Upon NF perception, the transcription factors NSP2 and NSP1 form a DNA binding complex [24,76]. This complex positively regulates the expression of the transcription factor NIN, which is required to activate the expression of different genes involved in rhizobia infection and colonization processes [77]. It has been demonstrated that mutations in the NSP2 gene significantly affect the rhizobia-induced root hair deformation and rhizobia infection process, particularly in the formation of the infection chamber and infection thread [74,78]. Likewise, early reports have demonstrated that M. truncatula nin mutants undergo excessive rhizobia-induced root hair curling but are impaired in infection and nodule formation [77].
Pi deficiency reduces both rhizobia-induced root hair deformation and infection thread formation in common bean seedlings suggesting that the expression of genes participating in the root hair curling and rhizobia infection is compromised under this nutritional condition. To test this hypothesis, we evaluated the expression of PvNSP2 (Phvul.009G122700), PvNIN (Phvul009G115800), and PvFLOT2 (Phvul.009G090700) in common bean seedlings growing under optimal- and low-Pi conditions and inoculated with rhizobia. Our analysis revealed that the expression of these three symbiosis-related genes was induced upon one-, 24, and 48 h post-inoculation with R. tropici in common bean seedlings growing under optimal Pi conditions (Figure 3). In contrast, we observed a 50% reduction in the expression of these three symbiotic genes in Pi-deficient common bean seedlings inoculated with rhizobia (Figure 3). Interestingly, at 48 h post-inoculation with rhizobia, PvNIN was the only gene that showed a substantial induction (in average by 5 fold-change) in Pi-deficient seedlings in response to rhizobia (Figure 3c).
To explore the mechanism behind the reduction in the fold-change expression of these three symbiosis-related genes in response to Pi-deficiency, we analyzed their expression levels in response to mock and rhizobia treatment (Figure 3b,d,f). By comparing the expression levels of PvNSP2 and PvNIN from both optimal-Pi and Pi-deficient seedlings, we did not find significant differences at one and 24 h post-mock treatment (Figure 3b,d). Similarly, we did not find significant differences in the expression levels of PvFLOT2 across the entire time-course when comparing optimal- and low-Pi seedlings post-mock treatment (Figure 3f). At 48 h, we observed an increase in the expression levels of PvNSP2 and a diminution in the expression levels of PvNIN in mock-inoculated Pi-deficient seedlings. By comparing the expression levels obtained from mock- and rhizobia-inoculated Pi-deficient seedlings, we observed a significant increase in the expression levels of PvNSP2 and PvNIN upon 24 and 48 h of inoculation with rhizobia (Figure 3b,d). Altogether, these transcriptome analyses suggest that the defects in root hair deformation and infection thread formation observed in common bean seedlings grown under Pi-deficient conditions are the consequences of an inefficient activation of critical transcription factors like PvNIN and PvNSP2, and symbiotic genes like PvFLOT2. The fact that substantial induction of the expression of PvNIN in Pi-deficient seedlings was detected upon 48 h post-inoculation with rhizobia only also indicates that Pi deficiency delays the activation of critical symbiosis-related genes.

3.3. Pi Deficiency Increases the Expression of the PvRIC1 and PvRIC2 Genes in Common Bean Seedlings

It has been reported that the expression of PvRIC1 and PvRIC2 is induced one- and five-days post-inoculation with rhizobia, respectively [32]. Based on the fact that Pi deficiency triggers a reduction in the number of nodules in common bean [40] and that the peptides PvRIC1 and PvRIC2 actively participate in the AON process, we hypothesized that the expression of PvRIC1 and PvRIC2 might be affected in Pi-deficient common bean seedlings. To test this hypothesis, we quantified the expression of PvRIC1 (Phvul.005G096901) and PvRIC2 (Phvul011G135900) in common bean seedlings grown under optimal- or low-Pi conditions, and under mock and rhizobia-treated conditions. This analysis revealed that the expression of PvRIC1 is induced as early as one-hour post inoculation with rhizobia in both control and Pi-deficient common bean seedlings and have its maximum induction level at 48 h post-inoculation with rhizobia (Figure 4a). In contrast, we were not able to detect a consistent and reproducible induction of PvRIC2 at 1 and 24 h post-inoculation with rhizobia. However, we observed a clear induction of the expression of PvRIC2 in response to rhizobia in both control and Pi-deficient common bean seedlings after 48 h post-inoculation (Figure 4c). Although the induction of the expression of these two AON-related genes was observed in both control and Pi-deficient seedlings, we observed that the expression levels of both PvRIC1 and PvRIC2 genes were always higher in Pi-deficient seedlings than in control seedlings (Figure 4b,d). Additionally, by comparing the expression levels of mock-inoculated seedlings, we observed that Pi-deficient seedlings always showed higher expression levels similar to those observed in control-seedlings inoculated with rhizobia (Figure 4b,d). Altogether, these data indicate that Pi deficiency induces the expression of PvRIC1 and PvRIC2 in the absence of rhizobia, and that the expression of these AON-related genes is highly induced in response to rhizobia in Pi-deficient common bean seedlings. Under this scenario, it is likely that: (1) Pi-deficient common bean seedlings are preconditioned to reduce the number of rhizobia infections. This hypothesis can be supported by the fact that the overexpression of RIC1 and RIC2 inhibits nodulation in soybean [79]. (2) Higher expression of PvRIC1 and PvRIC2 in response to Pi deficiency might be critical for common bean to properly adapt to Pi deprivation. This hypothesis is supported by the fact that other CLE peptides (e.g., CLE14) [80] play a critical role in the plant adaptation to Pi deficiency [81,82].

3.4. mRNA-Seq Analysis

To obtain a better insight into the genetic responses of common bean plants leading to a reduced nodulation under Pi-deficient conditions, we conducted an mRNA-seq analysis. For this genome-wide transcriptional analysis, based on our transcriptional data obtained by qRT-PCR, we selected 48 h post-inoculation to capture more transcriptional responses to rhizobia.
A total of 12 cDNA libraries were generated from control (optimal Pi: OP) and Pi-deficient roots (low Pi: LP) inoculated 48 h with mock (OPM or LPM, respectively) or rhizobia (OPR or LPR, respectively). Three independent biological replicates were generated for each condition. These libraries were sequenced using the Illumina NextSeq500 platform. After filtering low-quality reads, a total of 176,365,724 single-end reads (75 bp) were aligned to the P. vulgaris genome reference sequence Version 2.1 [59,60] using Bowtie and Tophat software and allowing two mismatches [64]. Of these, 164,037,947 were uniquely mapped to the common bean genome and were used for further analysis. The full RNA-seq datasets generated from the 12 cDNA libraries used in this study were deposited in the Gene Expression Omnibus [83] under the accession number GSE118968.

3.5. Global Transcriptional Responses of Pi-Deficient Common Bean Seedlings Interacting with Rhizobia

In order to identify additional genes with a potential role in the regulation of the early stages of the symbiosis between common bean and rhizobia under low-Pi conditions, the mRNA-seq data sets were analyzed using Cuffdiff [84] with an additional cutoff of 2-fold in pairwise comparisons (e.g., OPR/OPM). In total, 867 (511 up-regulated and 356 down-regulated) and 383 (285 up-regulated and 98 down-regulated) differentially regulated genes were identified when comparing OPR/OPM and LPR/LPM, respectively (Tables S2 and S3). Among them, PvNSP2 (Phvul.009G122700) and PvFLOT2 (Phvul.009G090700) were also classified as differentially regulated by mRNA-seq approach (Tables S2 and S3). However, we did not detect the expression of PvNIN (Phvul009G115800) in our mRNA-seq analysis, which we explain by a difference in the sensitivity of mRNA-seq and qRT-PCR approaches.
To confirm these mRNA-seq results, the expression of 25 randomly selected genes was analyzed via qRT-PCR and using independent biological material from the one used for mRNA-seq analysis (Figure S2). The expression pattern obtained by qRT-PCR showed the same trend observed by mRNA-seq in response to rhizobia. We found some differences in the fold-change computed by the two independent methods, which we explain by the relative sensitivity and data processing (e.g., algorithms used for data normalization) of each method, along with technical and biological aspects (e.g., efficiency/specificity of qRT-PCR and inherent variation in the responses of common bean to both Pi availability and rhizobia). In spite of these differences, the fold-changes computed for each analyzed gene by qRT-PCR were statistically different between both Pi conditions.
Upon confirmation of the reliability of our mRNA-seq datasets, we performed additional pairwise comparisons to identify genes that might explain the reduction in the nodule number triggered by Pi deficiency. To this end, we compared the differentially regulated genes found in the pairwise comparisons OPR/OPM and LPR/LPM to identify commonly and uniquely regulated genes under these Pi conditions and interacting with rhizobia. This comparative analysis revealed that 234 (180 up-regulated and 54 down-regulated) genes were commonly regulated in both Pi conditions while interacting with rhizobia (Figure 5a,b and Table S4). In contrast, 633 (331 up-regulated and 302 down-regulated) and 149 (105 up-regulated and 44 down-regulated) genes were uniquely regulated in Pi-optimal and Pi-deficient seedlings, respectively (Figure 5a,b and Tables S5 and S6).
Further hierarchical clustering analyses on the 234 common regulated genes revealed a clear variation in the expression levels computed for these genes (Figure 5c,d). For instance, we identified two clusters showing obvious differences when comparing the expression levels computed for both Pi-optimal and Pi-deficient seedlings inoculated with rhizobia (Figure 5c). The first identified cluster (cluster I) contains genes whose expression levels were higher in rhizobia-inoculated Pi-deficient seedlings than in rhizobia-inoculated optimal-Pi seedlings, whereas the second cluster (cluster II) contains genes whose expression levels were higher in optimal-Pi seedlings than in Pi-deficient seedlings, both of them inoculated with rhizobia (Figure 5c). To extend this analysis and identify genes differentially regulated between Pi-optimal and Pi-deficient seedlings while interacting with rhizobia, we performed a pair-wise comparison of the expression levels between rhizobia-inoculated optimal-Pi and rhizobia-inoculated Pi-deficient conditions. This comparison led us to the identification of 131 differentially expressed genes (77 showing higher expression in Pi-deficient seedlings and 54 showing higher expression in control seedlings) (Figure S3 and Table S7). The observed variation in the expression of these 131 genes was on average five-fold. For instance, Phvul.001G226900 and Phvul.008G244400, which encode for a Wall-associated receptor kinase galacturonan-binding protein and Leucine-Rich Repeat Protein Kinase-Like protein, respectively, showed in average by a four-fold increase in control common bean seedlings. In contrast, the genes Phvul.001G228500 and Phvul.011G074400, which encode for a U-Box domain-containing protein and Fasciclin-Like Arabinogalactan Protein11, respectively, showed on average by a six-fold increase in Pi-deficient common seedlings. Collectively, this transcriptional data clearly indicate that Pi deficiency affects the activity of different rhizobia-regulated genes by either inactivating, decreasing or increasing their expression in common bean.

3.6. Expression of Hormone- and Signal-Transduction-Related Genes Is Compromised in Rhizobia Inoculated Pi-Deficient Common Bean Plants

To gain more insights into the molecular functions of the differentially regulated genes identified in this study, we performed a gene-enrichment analysis and a functional classification in MapMan to identify the molecular mechanisms leading to an inefficient symbiosis between common bean and rhizobia under Pi-deficient conditions. This analysis indicated that Regulation of transcription, Hormones, Regulation/Signal transduction, Transport, Biotic interactions and Protein modification/degradation were the most overrepresented functional categories in both control and Pi-deficient common bean seedlings interacting with rhizobia (Figure 6).
Phytohormones play a critical role in different stages of the symbiosis between legumes and rhizobia [85]. For instances, it has been demonstrated that both auxin and cytokinin positively regulate the expression of genes that conforms the common symbiosis pathway, as well as the progression of the symbiotic signaling events from the epidermis to the cortex [27,86]. Recently, it was demonstrated that the auxin accumulation in the root hairs upon rhizobia recognition is essential for the infection thread formation [10]. Furthermore, it has also been demonstrated that the correct perception of cytokinin is crucial for the polar auxin transport and the subsequent auxin accumulation in the rhizobia-infected root hairs [87].
Based on the fact that auxin and cytokinin are crucial regulators of both early and late (e.g., nodule development) stages of the symbiosis between legume and rhizobia [88], and the fact that we found that the Hormones category was one of the overrepresented functional categories, we sought in our transcriptomic data to determinate whether the expression of genes related to these two hormones was affected by Pi deficiency. This analysis revealed that 12 auxin-related genes (six of them involved in auxin biosynthesis, one gene involved in auxin homeostasis and five auxin-responsive genes) and eight cytokinin-related genes (five of them involved in cytokinin homeostasis, two involved in cytokinin biosynthesis, and one involved in cytokinin signaling) were differentially regulated in optimal-Pi common bean seedlings and inoculated with rhizobia (Figure 7a,b). In contrast, only four auxin-related genes (two of them involved in auxin biosynthesis, one involved in auxin homeostasis, and one auxin-responsive gene) and one cytokinin-related (cytokinin homeostasis) gene were differentially regulated in response to rhizobia when common bean seedlings were grown under Pi-deficient condition (Figure 7). The fact that Pi deficiency affects the expression of auxin- and cytokinin-related genes in response to rhizobia, lead us to propose that the defects in early molecular (e.g., expression of PvNIN and PvFLOT2) and physiological (e.g., root hair deformation) events observed in Pi-deficient common bean seedlings might be partially explained by defects in the auxin-cytokinin balance.
In addition to auxin- and cytokinin-related genes, we also found that genes involved in the biosynthesis of the phytohormones ethylene, abscisic acid (ABA)- and jasmonate, as well as in the signal transduction related to these three phytohormones, were differentially regulated in response to rhizobia inoculation. Although we observed that common bean seedlings growing under optimal-Pi conditions and inoculated with rhizobia showed more differentially regulated genes (six ABA-related genes, seven ethylene-related genes, and six jasmonate-related genes) related to these three hormones than Pi-deficient seedlings inoculated with rhizobia, the majority of them were down-regulated (Figure 8). In contrast, all ethylene-, ABA- and jasmonate-related genes identified in common bean seedlings grown under Pi-deficient conditions and inoculated with rhizobia were up-regulated (Figure 8).
Ethylene, jasmonate, and ABA are essential to activate a series of molecular and physiological responses allowing plants to cope with Pi deficiency [89]. However, growing evidence indicates that these three phytohormones negatively control the establishment of the root nodule symbiosis [90,91,92]. For instances, it has been demonstrated that the accumulation of ethylene or jasmonate inactivates the expression of early symbiosis-related genes (e.g., ENOD11), inhibits the initiation and maintenance of the calcium spiking, decreases the number of rhizobia-induced root hair deformation and infection thread, and reduces nodule formation in M. truncatula plants [90,91]. Recently, it was demonstrated that the negative impact of ethylene on the rhizobia infection process and nodule formation is under the control of the regulator EIN2 [28]. Hence, our transcriptional data suggest that the activation of genes related to the biosynthesis and signal transduction of these three phytohormones is required in order to allow common bean to properly adapt to Pi deficiency. However, the activation of the ethylene-, jasmonate- and ABA signaling pathways compromise the activation of early molecular and physiological events of the root nodule symbiosis, as we demonstrated in this study.
The proper recognition of rhizobia by the host legume is crucial for a successful symbiosis. To this end, legumes rely on the LysM-domain receptor kinases NFR5 [93], NFR1 [94] and eNFR [7], as well as in the LRR receptor-like kinase SYMRK, to properly detect and decode the NF signal [15,16]. Recently, it was demonstrated that a fourth LysM receptor kinase receptor, EPR3, is required for the perception of compatible rhizobial exopolysaccharides and promotes a successful infection and colonization [95]. Additionally, based on the fact that a gene encoding cell wall-associated receptor kinases (WAKs) is significantly up-regulated upon one-hour treatment with NF in M. truncatula roots, it has been proposed that WAKs might play a critical role in the early steps of this symbiosis, primarily in the cell wall-cytoplasm signaling [96]. In our transcriptomic analysis, we also found different types of protein receptors, including LysM receptor kinases, LRR receptor kinases, and WAKs, that were differentially regulated in both control and Pi-deficient seedlings inoculated with rhizobia. However, we observed that Pi-deficient seedlings inoculated with rhizobia showed less differentially regulated genes encoding for protein receptors in response to rhizobia (Figure S4). These data suggest that the misregulation of the expression of additional protein receptor might contribute to an inefficient communication between common bean and rhizobia under Pi-limiting conditions.

4. Conclusions

The data presented in this study lead us to conclude that Pi deficiency negatively affects early events of the symbiosis between legumes and rhizobia. The negative effects in the symbiosis between the common bean and rhizobia under Pi deficiency might be triggered by an imbalance between auxin and cytokinin, as well as by the activation of the ethylene, jasmonate and ABA-related signaling pathway, which are negative regulators of early molecular events of this symbiosis (Figure 9). Finally, the reduction in the number of nodules observed under Pi-deficient conditions can be also explained by a constitutive activation of the AON process (Figure 9); however, further experimentation is needed to confirm this hypothesis.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/10/498/s1, Figure S1: In vitro system to assess early responses to rhizobia on Pi-deficient common bean seedlings. Seedlings were grown under OP (a) or LP (b) conditions. Seedlings growing under LP conditions show a 50% decrease in their Pi content (c). Asterisk indicates a significant difference according to One-way ANOVA (p-value < 0.001), Figure S2: qRT-PCR validation of the mRNA-seq datasets. Data showed are the fold-change obtained from the pairwise comparison between the expression levels obtained from seedlings grown under optimal-Pi (OP) or low-Pi (LP) inoculated with rhizobia (R; OPR or LPR) or mock (M; OPM or LPM), Figure S3: Genes showing significantly higher expression levels in Pi-deficient common seedlings inoculated with rhizobia. Data showed are the expression levels obtained from common bean seedlings grown under optimal-Pi (OP) or low-Pi (LP) conditions and inoculated with rhizobia. Panels on the right site show the functional categories regulated under the described experimental conditions. TF: transcription factors; PM: protein modification; PD: protein degradation; IAA: auxin; ABA: abscisic acid; Eth: ethylene; CK: cytokinins; JA: jasmonate; SA: salicylic acid; GA: gibberellic acid; RK: receptor kinases; CaR: calcium regulation; GP: G-protein; Pi: phosphoinositides; C and Nutr: C and nutrients; MK: MAP kinases; U: unspecified; Thx: thioredoxin; Asc/Glut: arcorbic/gluath; Glutx: glutaredoxin; Px: periredoxin; Dism/Cat: dismutase/catalase, Figure S4: Expression pattern of protein receptors- and secondary metabolism-related genes. Data showed are the fold-change obtained from the pairwise comparison between the expression levels obtained from seedlings grown under optimal-Pi (OP) or low-Pi (LP) inoculated with rhizobia (R; OPR or LPR) or mock (M; OPM or LPM). Fold-changes were false-colored; genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in green color, Table S1: Primer sequences used for quantitative real-time polymerase chain reaction (qRT-PCR) analysis, Table S2: Rhizobia-induced genes under optimal-Pi conditions, Table S3: Rhizobia-induced genes under low-Pi conditions, Table S4: Commonly regulated genes in both optimal- and low-Pi conditions, Table S5: Genes uniquely regulated under optimal-Pi conditions, Table S6: Genes uniquely regulated under low-Pi conditions, Table S7: Comparison of the expression levels of rhizobia-induced genes under OP or LP conditions.

Author Contributions

O.V.-L. and M.L. proposed and designed the study. M.C.I.-A. and L.G. generated the R. tropici::uidA strain. M.C.I.-A., M.d.R.R.-S. and M.d.S.S.-C. analyzed the rhizobia-induced root hair deformation. M.C.I.-A. analyzed the infection thread formation. M.C.I.-A. and M.d.R.R.-S., analyzed the expression symbiosis-related genes. M.C.I.-A., M.d.R.R.-S., M.d.S.S.-C., D.G.M.-C., and N.A.C.-G. generated biological material for mRNAseq analysis. S.S., T.J., L.P., O.V.-L., and M.C.I.-A. analyzed the mRNAseq data. O.V.-L., M.L., M.C.I.-A. wrote the manuscript. All authors read and approved the manuscript.

Funding

This research was supported by a grant from Programa de Apoyo a Proyectos de Investigación e Inovación Tecnológica (PAPIIT grant #IN213017) to OVL. This work was also supported by the United States National Sciences Foundation grant numbers IOS-1339194 and IOS-1453613 to ML, and IOS-1734145 to DMC. Mariel C. Isidra-Arellano is a Doctoral student from Programa de Doctorado en Ciencias Biológicas, Universidad Nacional Autónoma de México (UNAM) and receives a fellowship from Consejo Nacional de Ciencia y Tecnologia (CONACyT, CVU: 919676). Maria del Rocio Reyero-Saavdra is a Doctoral student from Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM) and receives a fellowship from CONACyT (347027/239879).

Acknowledgments

We thank the Joint Genome Institute for prepublication access to the Phaseolus vulgaris V2.1 genome sequence. High-throughput services were performed at the University of Missouri DNA Core Facility. We would like to thank Jose L. Reyes (Instituto de Biotecnología UNAM), Georgina Hernandez (Centro de Ciencias Genómicas, UNAM) and Kiwamu Tanaka (Washington State University, Pullman, USA) for constructive discussion.

Conflicts of Interest

Authors have no conflict of interests to declare.

References

  1. Castro-Guerrero, N.A.; Isidra-Arellano, M.C.; Mendoza-Cozatl, D.M.; Valdés-López, O. Common bean: A legume model on the rise for unraveling adaptations to iron, zinc and phosphate deficiencies. Front. Plant Sci. 2016, 7, 600. [Google Scholar] [CrossRef] [PubMed]
  2. Blair, M.W.; Izquierdo, P.; Astudillo, C.; Grusak, M.A. A legume biofortification quandary: Variability and genetic control of seed coat micronutrient accumulation in common beans. Front. Plant Sci. 2013, 4, 275. [Google Scholar] [CrossRef] [PubMed]
  3. Jha, A.B.; Ashokkumar, K.; Diapari, M.; Ambrose, S.J.; Zhang, H.; Tar’an, B.; Bett, K.E.; Vandenberg, A.; Warkentin, T.D.; Purves, R.W. Genetic diversity of folate profiles in seeds of common bean, lentil, chickpea and pea. J. Food Compos. Anal. 2015, 42, 134–140. [Google Scholar] [CrossRef]
  4. Peters, N.K.; Frost, J.N.; Long, S.R. A plant flavonone, loteolin, induces expression of Rhizobium meliloti genes. Sciences 1996, 233, 977–980. [Google Scholar] [CrossRef]
  5. Liu, C.W.; Murray, J. The role of flavonoids in nodulation host-range specificity and update. Plants 2016, 5, 33. [Google Scholar] [CrossRef] [PubMed]
  6. Broghanmer, A.; Krusell, L.; Blaise, M.; Sauer, J.; Sullivan, J.T.; Maolanon, N.; Vinther, M.; Lorentzen, A.; Madsen, E.B.; Jensen, K.J.; et al. Legume receptors perceive the rhizobial lipochitin oligosaccharide signal molecules by direct binding. Proc. Natl. Acad. Sci. USA 2012, 109, 13859–13864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Murakami, E.; Cheng, J.; Gysel, K.; Bozsoki, Z.; Kawaharada, Y.; Hjuler, C.T.; Sørensen, K.K.; Tao, K.; Kelly, S.; Venice, F.; et al. Epidermal LysM receptor ensure robust symbiotic signalling in Lotus japonicus. eLife 2018, 7, e33506. [Google Scholar] [CrossRef] [PubMed]
  8. Venkateshwaran, M.; Volkening, J.D.; Sussman, M.R.; Ané, J.M. Symbiosis and the social network of higher plants. Curr. Opin. Plant Biol. 2013, 16, 118–127. [Google Scholar] [CrossRef] [PubMed]
  9. Lei, M.J.; Wang, Q.; Li, X.; Chen, A.; Luo, L.; Xie, Y.; Li, G.; Luo, D.; Mysore, K.S.; Wen, J.; et al. The small GTPase ROP10 of Medicago truncatula is required for both tip growth of root hairs and Nod Factor-induced root hair deformation. Plant Cell 2015, 27, 806–822. [Google Scholar] [CrossRef] [PubMed]
  10. Nadzieja, M.; Kelly, S.; Stougaard, J.; Reid, D. Epidermal auxin biosynthesis facilitates rhizobial infection in Lotus japonicus. Plant J. 2018, 95, 101–111. [Google Scholar] [CrossRef] [PubMed]
  11. Fournier, J.; Timmers, A.C.J.; Sieberer, B.J.; Jauneau, A.; Chabaud, M.; Barker, D.G. Mechanism of infection thread elongation in root hairs of Medicago truncatula and dynamic interplay with associated rhizobial colonization. Plant Physiol. 2008, 148, 1985–1995. [Google Scholar] [CrossRef] [PubMed]
  12. Fournier, J.; Teillet, A.; Chabaud, M.; Ivanov, S.; Genre, A.; Limpens, E.; de Carvalho-Nievel, F.; Barker, D. Remodeling of the infection chamber before infection thread formation reveals a two-step mechanism for rhizobial entry into the host legume root hair. Plant Physiol. 2015, 167, 1233–1242. [Google Scholar] [CrossRef] [PubMed]
  13. Xiao, T.T.; Schilderink, S.; Moling, S.; Deinum, E.E.; Kondorosi, E.; Franssen, H.; Kulikova, O.; Niebel, A.; Bisseling, T. Fate map of Medicago truncatula root nodules. Development 2014, 141, 3517–3528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Genre, A.; Russo, G. Does a common pathway transduce symbiotic signals in plant-microbe interactions? Front. Plant Sci. 2016, 7, 96. [Google Scholar] [CrossRef] [PubMed]
  15. Bersoult, A.; Camut, S.; Perhald, A.; Kereszt, A.; Kiss, G.B.; Cullimore, J.V. Expression of the Medicago truncatula DM12 gene suggests roles of the symbiotic nodulation receptor kinase in nodules and during early nodule development. Mol. Plant Microbe Interact. 2005, 18, 869–876. [Google Scholar] [CrossRef] [PubMed]
  16. Hogg, B.V.; Cullimore, J.V.; Ranjeva, R.; Bono, J.J. The DMI1 and DMI2 early symbiotic genes of Medicago truncatula are required for a high-affinity nodulation factor-binding site associate to a particulate fraction of roots. Plant Physiol. 2006, 140, 365–373. [Google Scholar] [CrossRef] [PubMed]
  17. Imaizumi-Anraku, H.; Takeda, N.; Charpentier, M.; Perry, J.; Miwa, H.; Umehara, Y.; Kouchi, H.; Murakami, Y.; Mulder, L.; Vickers, K.; et al. Plastid proteins crucial for symbiotic fungal and bacterial entry into plant host. Nature 2005, 433, 527–531. [Google Scholar] [CrossRef] [PubMed]
  18. Kanamori, N.; Madsen, L.H.; Radatoiu, S.; Frantescu, M.; Quistgaard, E.M.; Miwa, H.; Downie, J.A.; James, E.K.; Felle, H.H.; Haaning, L.L.; et al. A nucleoporin is required for induction of Ca2+ spiking in legume nodule development and essential for rhizobial and fungal symbiosis. Proc. Natl. Acad. Sci. USA 2006, 103, 359–364. [Google Scholar] [CrossRef] [PubMed]
  19. Saito, K.; Yoshikawa, M.; Yano, K.; Miwa, H.; Uchida, H.; Asamizu, E.; Sato, S.; Tabata, S.; Imaizumi-Anraku, H.; Kouchi, H.; et al. NUCLEOPORIN85 is required for calcium spiking, fungal and bacterial symbioses, and seed production in Lotus japonicus. Plant Cell 2007, 19, 610–624. [Google Scholar] [CrossRef] [PubMed]
  20. Groth, M.; Takeda, N.; Perry, J.; Uchida, H.; Dräxl, S.; Brachmann, A.; Sato, S.; Tabata, S.; Kawaguchi, M.; Wang, T.L.; et al. NENA, a Lotus japonicus homolog of Sec13, is required for rhizodermal infection by arbuscular mycorrhiza fungi and rhizobia but dispensable for cortical endosymbiotic development. Plant Cell 2010, 22, 2509–2526. [Google Scholar] [CrossRef] [PubMed]
  21. Charpentier, M.; Sun, J.; Vaz Martins, T.; Radhakrishnan, G.V.; Findlay, K.; Soumpourou, E.; Thouin, J.; Véry, A.A.; Sanders, D.; Morris, R.J.; et al. Nuclear-localized cyclic nucleotide-gated channels mediate symbiotic calcium oscillations. Science 2016, 352, 1102–1105. [Google Scholar] [CrossRef] [PubMed]
  22. Miller, J.B.; Pratap, A.; Miyahara, A.; Zhou, L.; Borneman, S.; Morris, R.J.; Oldroyd, G.E. Calcium/Calmodulin-dependent protein kinase is negatively and positively regulated by calcium, providing a mechanism for decoding calcium responses during symbiosis signaling. Plant Cell 2013, 25, 5053–5066. [Google Scholar] [CrossRef] [PubMed]
  23. Singh, S.; Katzer, K.; Lambert, J.; Cerri, M.; Parniske, M. CYCLOPs, a DNA-binding transcriptional activator orchestrates symbiotic root nodule development. Cell Host Microbe 2014, 15, 139–152. [Google Scholar] [CrossRef] [PubMed]
  24. Hirsch, S.; Kim, J.; Muñoz, A.; Heckmann, A.B.; Downie, J.A.; Oldroyd, G.E.D. GRAS proteins form a DNA binding complex to induce gene expression during nodulation signaling in Medicago truncatula. Plant Cell 2009, 21, 545–557. [Google Scholar] [CrossRef] [PubMed]
  25. Soyano, T.; Kouchi, H.; Hirota, A.; Hayashi, M. Nodule Inception directly targets NF-Y sub-unit genes to regulate essential processes of root nodule development in Lotus japonicus. PLoS Genet. 2013, 9, e1003352. [Google Scholar] [CrossRef] [PubMed]
  26. Haney, C.H.; Long, S.R. Plant flotillins are required for infection by nitrogen-fixing bacteria. Proc. Natl. Acad. Sci. USA 2010, 107, 478–483. [Google Scholar] [CrossRef] [PubMed]
  27. van Zeijl, A.; Op den Camp, R.H.M.; Deinum, E.E.; Charnikhova, T.; Franssen, H.; Op den Camp, H.J.M.; Bouwmeester, H.; Kohlen, W.; Bisseling, T.; Geurts, R. Rhizobium lipo-chitooligosaccharide signaling triggers accumulation of cytokinins in Medicago truncatula roots. Mol. Plant 2015, 8, 1213–1226. [Google Scholar] [CrossRef] [PubMed]
  28. Reid, D.; Liu, H.; Kelly, S.; Kawaharada, Y.; Mun, T.; Andersen, S.U.; Desbrosses, G.; Stougaard, J. Dynamics of Ethylene production in response to compatible Nod Factor. Plant Physiol. 2018, 176, 1764–1772. [Google Scholar] [CrossRef] [PubMed]
  29. Penmetsa, V.R.; Cook, D.R. A legume ethylene-insensitive mutant hyperinfected by its rhizobial symbiont. Science 1997, 275, 527–530. [Google Scholar] [CrossRef] [PubMed]
  30. Trichine, L.; Sandal, N.; Madsen, L.H.; Radutoiu, S.; Albrektsen, A.S.; Sato, S.; Asamizu, E.; Tabata, S.; Stougaard, J. A gain-of-function mutation in a cytokinin receptor triggers spontaneous root nodule organogenesis. Science 2007, 315, 104–107. [Google Scholar] [CrossRef] [PubMed]
  31. Ferguson, B.J.; Mens, C.; Hastwell, A.H.; Zhang, M.; Su, H.; Jones, C.H.; Chu, X.; Gresshoff, P.M. Legume nodulation: The host controls the party. Plant Cell Environ. 2018, 1, 11. [Google Scholar] [CrossRef] [PubMed]
  32. Ferguson, B.J.; Dongxue, L.; Hastwell, A.H.; Reid, D.E.; Li, Y.; Jackson, S.A.; Gresshoff, P. The soybean (Glycine max) nodulation-suppressive CLE peptide, GmRIC1, functions interspecifically in common white bean (Phaseolus vulgaris), but not in a suppernodulating line mutated in the receptor PvNARK. Plant Biotechnol. J. 2014, 12, 1085–1097. [Google Scholar] [CrossRef] [PubMed]
  33. Krusell, L.; Madsen, L.H.; Sato, S.; Aubert, G.; Genua, A.; Szczyglowski, K.; Duc, G.; Kaneko, T.; Tabata, S.; de Bruijn, F.; et al. Shoot control of root development and nodulation is mediated by a receptor-like kinase. Nature 2002, 420, 422–426. [Google Scholar] [CrossRef] [PubMed]
  34. Nishimura, R.; Hayashi, M.; Wu, G.J.; Kouchi, H.; Imaizumi-Anraku, H.; Murakami, Y.; Kawasaki, S.; Akao, S.; Ohmori, M.; Nagasawa, M.; et al. HAR1 mediates systemic regulation of symbiotic organ development. Nature 2002, 420, 426–429. [Google Scholar] [CrossRef] [PubMed]
  35. Searle, I.R.; Men, A.E.; Laniya, T.S.; Buzas, D.M.; Iturbe-Ormaetxe, I.; Carroll, B.J.; Gresshoff, P.M. Long-distance signaling in nodulation directed by a CLAVATA1-like receptor kinase. Science 2003, 299, 109–112. [Google Scholar] [CrossRef] [PubMed]
  36. Schnabel, E.; Journet, E.P.; de Carvalho-Niebel, F.; Duc, G.; Frugoli, J. The Medicago truncatula SUNN gene encodes a CLV1-like leucine-rich repeat receptor kinase that regulates nodule number and root length. Plant Mol. Biol. 2005, 58, 809–822. [Google Scholar] [CrossRef] [PubMed]
  37. Broughton, W.J.; Hernandez, G.; Blair, M.; Beebe, S.; Gepts, P.; Vanderleyden, J. Beans (Phaseolus spp.)—Model food legumes. Plant Soil 2003, 252, 55–128. [Google Scholar] [CrossRef]
  38. Valdés-López, O.; Hernández, G. Transcriptional regulation and signaling in phosphorus starvation: What about legumes? J. Integr. Plant Biol. 2008, 50, 1213–1222. [Google Scholar] [CrossRef] [PubMed]
  39. Vance, C.P.; Uhde-Stone, C.; Allan, D.L. Phosphorus acquisition and use: Critical adaptations by plants for securing a nonrenewable resource. New Phytol. 2003, 157, 423–447. [Google Scholar] [CrossRef]
  40. Hernández, G.; Valdés-López, O.; Ramírez, M.; Goffard, N.; Weiller, G.; Aparicio-Fabre, R.; Fuentes, S.I.; Erban, A.; Kopka, J.; Udvardi, M.K.; et al. Global changes in the transcript and metabolic profiles during symbiotic nitrogen fixation in phosphorus-stressed common bean plants. Plant Physiol. 2009, 151, 1221–1238. [Google Scholar] [CrossRef] [PubMed]
  41. Hernandez, G.; Devron, J.J. Influence of oxygen and acetylene during in situ open-flow assays of nitrogenase activity (C2H2 reduction) in Phaseolus vulgaris root nodule. J. Plant Physiol. 1991, 138, 587–590. [Google Scholar] [CrossRef]
  42. Vadez, V.; Rodier, F.; Payre, H.; Devon, J.J. Nodule conductance to O2 and nitrogenase-linked respiration in bean genotypes varying in the tolerance of N2 fixation to P deficiency. Plant Physiol. Biochem. 1996, 34, 871–878. [Google Scholar]
  43. Cabeza, R.A.; Liese, R.; Lingner, A.; von Steieglitz, I.; Neumann, J.; Salinas-Riester, G.; Pommerenke, C.; Dittert, K.; Schulze, J. RNA-seq transcriptome profiling reveals that Medicago trucatula nodules acclimate N2 fixation before emerging P deficiency reaches the nodules. J. Exp. Bot. 2014, 65, 6035–6048. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, Z.; Cui, Q.; Liang, C.; Sun, L.; Tian, J.; Liao, H. Identification of differentially expressed proteins in soybean nodules under phosphorus deficiency through proteomics analysis. Proteomics 2011, 11, 4648–4659. [Google Scholar] [CrossRef] [PubMed]
  45. Suleiman, S.; Ha, C.V.; Schulze, J.; Tran, L.S.P. Growth and nodulation of symbiotic Medicago truncatula at different levels of phosphorus availability. J. Exp. Bot. 2013, 64, 2701–2712. [Google Scholar] [CrossRef] [PubMed]
  46. Nasr Esfahani, M.; Inoue, K.; Chu, H.D.; Nguyen, K.H.; Van Ha, C.; Watanabe, Y.; Burritt, D.J.; Herrera-Estrella, L.; Mochida, K.; Tran, L.P. Comparative transcriptome analysis of nodules of two Mesorhizobium-chickpea associations with differential symbiotic efficiency under phosphate deficiency. Plant J. 2017, 91, 911–926. [Google Scholar] [CrossRef] [PubMed]
  47. Sa, T.; Israel, D.W. Energy status and functioning of phosphorus-deficient soybean nodules. Plant Physiol. 1991, 97, 928–935. [Google Scholar] [CrossRef] [PubMed]
  48. Al-Niemi, T.S.; Kahn, M.L.; McDermott, T.R. P metabolism in the bean-Rhizobium tropicisymbiosis. Plant Physiol. 1997, 113, 1233–1242. [Google Scholar] [CrossRef] [PubMed]
  49. Vadez, V.; Beck, D.P.; Lasso, J.H.; Devron, J.J. Utilization of the acetylene reduction assay to screen for tolerance of symbiotic N2 fixation to limiting P nutrition in common bean. Physiol. Plant 1997, 99, 227–232. [Google Scholar] [CrossRef]
  50. Tang, C.; Hinsinger, P.; Drevon, J.J.; Jaillard, B. Phosphorus deficiency impairs early nodule functioning and enhances proton release in roots of Medicago truncatula L. Ann. Bot. 2001, 88, 131–138. [Google Scholar] [CrossRef]
  51. Hogh-Jensen, H.; Schjoerring, J.K.; Soussana, J.F. The influence of phosphorus deficiency on growth and nitrogen fixation of white clover plants. Ann. Bot. 2002, 90, 745–753. [Google Scholar] [CrossRef] [PubMed]
  52. Schulze, J.; Temple, G.; Temple, S.; Beschow, H.; Vance, C.P. Nitrogen fixation by white lupin under phosphorus deficiency. Ann. Bot. 2006, 98, 731–740. [Google Scholar] [CrossRef] [PubMed]
  53. Catoira, R.; Galera, C.; de Billy, F.; Penmetsa, R.V.; Journet, E.P.; Maillet, F.; Rosenberg, C.; Cood, D.; Gough, C.; Dénarié, J. Four genes of Medicago truncatula controlling components of a nod factor transduction pathway. Plant Cell 2000, 12, 1647–1665. [Google Scholar] [CrossRef] [PubMed]
  54. Dalla Via, V.; Narduzzi, C.; Aguilar, O.M.; Zanetti, M.E.; Blanco, A. Changes in the common bean transcriptome in response to secreted and surface signal molecules of Rhizobium etli. Plant Physiol. 2015, 169, 1356–1370. [Google Scholar] [CrossRef] [PubMed]
  55. Taussky, H.H.; Shorr, E. A microcolorimetric method for determination of inorganic phosphorus. J. Biol. Chem. 1953, 202, 675–685. [Google Scholar] [PubMed]
  56. Reyero-Saavedra, M.D.R.; Qiao, Z.; Sánchez-Correa, M.D.S.; Díaz-Pineda, M.E.; Reyes, J.L.; Covarrubias, A.A.; Libault, M.; Valdés-López, O. Gene silencing of Argonaute5 negatively affects the establishment of the legume-rhizobia. Genes 2017, 8, 352. [Google Scholar] [CrossRef] [PubMed]
  57. Ramakers, C.; Ruijter, M.J.; Deprez, H.R.; Moorman, A.F. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 2003, 339, 62–66. [Google Scholar] [CrossRef]
  58. Primer3. Available online: http://bioinfo.ut.ee/primer3-0.4.0/primer3/ (accessed on 10 July 2017).
  59. Schmutz, J.; McClean, P.E.; Mamidi, S.; Wu, G.A.; Cannon, S.B.; Grimwood, J.; Chavarro, C.; Torres-Torres, M.; Geffroy, V.; Moghaddam, S.M.; et al. A reference genome for common bean and geno-wide analysis of dual domestications. Nat. Genet. 2014, 46, 707–713. [Google Scholar] [CrossRef] [PubMed]
  60. Phytozome. Available online: https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Pvulgaris (accessed on 8 June 2017).
  61. Barbraham Bionformatics. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 8 June 2017).
  62. Babraham Bioinformatics. Available online: http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ (accessed on 8 June 2017).
  63. Langmead, B.; Salzberg, S.L. Fast Gapped-Read Alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  64. Trapnell, C.; Pachter, L.; Sazlberg, S.L. TopHat: Discovering splice junctions with RNA-seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  65. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, C.; Salzberg, S.L.; Rinn, J.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protocols 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Du, Z.; Zhou, X.; Liang, Y.; Zhang, Z.; Su, Z. agriGO: A GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38, W64–W70. [Google Scholar] [CrossRef] [PubMed]
  67. AGRIGO. Available online: http://systemsbiology.cau.edu.cn/agriGOv2/ (accessed on 20 August 2017).
  68. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. REVIGO summarizes and visualizes long list of gene ontology terms. PLoS ONE 2011, 6, e21800. [Google Scholar] [CrossRef] [PubMed]
  69. REVIGO. Available online: http://revigo.irb.hr (accessed on 20 August 2017).
  70. Thimm, O.; Bläsing, O.; Gibon, Y.; Nagel, A.; Meyer, S.; Krüger, P.; Selbig, J.; Müller, L.A.; Rhee, S.Y.; Stitt, M. MAPMAN: A user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37, 914–939. [Google Scholar] [CrossRef] [PubMed]
  71. Usadel, B.; Poree, F.; Nagel, A.; Lohse, M.; Czedik-Eysenberg, A.; Stitt, M. A guide to using MapMan to visualize and compare Omics data in plants: A case study in the crop species, maize. Plant Cell Environ. 2009, 32, 1211–1229. [Google Scholar] [CrossRef] [PubMed]
  72. MapMan. Available online: https://mapman.gabipd.org (accessed on 15 April 2018).
  73. Qiu, L.; Lin, J.S.; Xu, J.; Sato, S.; Parniske, M.; Wang, J.; Xie, F. SCARN a novel class of SCAR protein that is required for root-hair infection during legume nodulation. PLoS Genet. 2015, 11, e1005623. [Google Scholar] [CrossRef] [PubMed]
  74. Cerri, M.R.; Frances, L.; Laloum, T.; Auriac, M.C.; Niebel, A.; Oldroyd, G.E.; Barke, D.G.; Fournier, J.; de Carvalho-Niebel, F. Medicago truncatula ERN transcription factors: Regulatory interplay with NSP1/NSP2 GRAS factors and expression dynamics throughout rhizobia infection. Plant Physiol. 2012, 160, 2155–2172. [Google Scholar] [CrossRef] [PubMed]
  75. Howieson, J.G.; Robson, A.D.; Ewing, M.A. External phosphate and calcium concentrations and pH, but not the products of rhizobial nodulation genes, affects the attachment of Rhizobium meliloti to roots of annual medics. Soil Biol. Biochem. 1993, 25, 567–573. [Google Scholar] [CrossRef]
  76. Smit, P.; Raedts, J.; Portyanko, V.; Debellé, F.; Gough, C.; Bisseling, T.; Geurts, R. NSP1 of the GRAS protein family is essential for rhizobial Nod factor-induced transcription. Science 2005, 308, 1789–1791. [Google Scholar] [CrossRef] [PubMed]
  77. Marsh, J.F.; Rakocevic, A.; Mitra, R.M.; Brocard, L.; Sun, J.; Eschstruth, A.; Long, S.; Schultze, M.; Ratet, P.; Oldroyd, G.E.D. Medicago truncatula NIN is essential for rhizobial-independent nodule organogenesis induced by autoactive calcium/calmodulin-dependent protein kinase. Plant Physiol. 2007, 144, 324–335. [Google Scholar] [CrossRef] [PubMed]
  78. Oldroyd, G.E.; Long, S.R. Identification and characterization of nodulation-signalling pathway 2, a gene of Medicago truncatula involved in Nod Factor signaling. Plant Physiol. 2003, 131, 1027–1032. [Google Scholar] [CrossRef] [PubMed]
  79. Reid, D.E.; Ferguson, B.J.; Gresshoff, P.M. Inoculation- and nitrate-induced CLE peptides of soybean control NARK-dependent Nodule formation. Mol. Plant Microbe Interact. 2011, 24, 606–618. [Google Scholar] [CrossRef] [PubMed]
  80. Funayama-Noguchi, S.; Noguchi, K.; Yoshida, C.; Kawaguchi, M. Two CLE genes are induced by phosphate in roots of Lotus japonicus. J. Plant Res. 2011, 124, 155–163. [Google Scholar] [CrossRef] [PubMed]
  81. de Bang, T.C.; Lundquist, P.K.; Dai, X.; Boschiero, C.; Zhuang, Z.; Pant, P.; Torres-Jerez, I.; Roy, S.; Nogales, J.; Veerappan, V.; et al. Genome-Wide identification of Medicago peptides involved in macronutrient responses and nodulation. Plant Physiol. 2017, 175, 1669–1689. [Google Scholar] [CrossRef] [PubMed]
  82. Gutierrez-Alanis, D.; Yong, L.; Jiménez, P.; Alatorre-Cobos, F.; Oropeza-Aburto, A.; Mora-Macías, J.; Sánchez-Rodríguez, F.; Cruz-Ramirez, A.; Herrera-Estrella, L. Phosphate Starvation-Dependent Iron Mobilization Induces CLE14 Expression to Trigger Root Meristem Differentiation through CLV2/PEPR2 Signaling. Dev. Cell 2017, 41, 555–570. [Google Scholar] [CrossRef] [PubMed]
  83. GEO. Available online: https://www.ncbi.nlm.nih.gov/geo/ (accessed on 23 August 2018).
  84. Trapnell, C.; Hendrickson, D.G.; Sauvageau, M.; Goff, L.; Rinn, J.L.; Pachter, L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat. Biotechnol. 2013, 31, 46–53. [Google Scholar] [CrossRef] [PubMed]
  85. Ferguson, B.J.; Mathesius, U. Phytohormone regulation of legume-rhizobia interactions. J. Chem. Ecol. 2014, 40, 770–790. [Google Scholar] [CrossRef] [PubMed]
  86. Hirsch, A.M.; Bhuvaneswari, T.V.; Torrey, J.G.; Bisseling, T. Early nodulin genes are induced in alfalfa roots outgrowths elicited by auxin transport inhibitors. Proc. Natl. Acad. Sci. USA 1989, 86, 1244–1248. [Google Scholar] [CrossRef] [PubMed]
  87. Ng, J.L.; Hassan, S.; Truong, T.T.; Hocart, C.H.; Laffont, C.; Frugier, F.; Mathesius, U. Flavonoids and auxin transport inhibitors rescue symbiotic nodulation in Medicago truncatula cytokinin perception cre1. Plant Cell 2015, 27, 2210–2226. [Google Scholar] [CrossRef] [PubMed]
  88. Liu, H.; Zhang, C.; Yang, J.; Yu, N.; Wang, E. Hormone modulation of legume-rhizobial symbiosis. J. Integr. Plant Biol. 2018, 60, 632–648. [Google Scholar] [CrossRef] [PubMed]
  89. Rouached, H.; Arpat, A.B.; Poirier, Y. Regulation of phosphate starvation responses in plants: Signaling players and cross-talks. Mol. Plant 2010, 3, 288–299. [Google Scholar] [CrossRef] [PubMed]
  90. Oldroyd, G.E.D.; Engstrom, E.M.; Long, S.R. Ethylene inhibits the Nod factor signaling transduction pathway of Medicago truncatula. Plant Cell 2001, 13, 1835–1849. [Google Scholar] [CrossRef] [PubMed]
  91. Sun, J.; Cardoza, V.; Mitchell, D.M.; Brigth, L.; Oldroyd, G.; Harris, J.M. Crosstalk between jasmonic acid, ethylene and Nod factor signaling allows integrations of diverse inputs for regulation of nodulation. Plant J. 2006, 46, 961–970. [Google Scholar] [CrossRef] [PubMed]
  92. Ding, Y.; Kalo, P.; Yendrek, C.; Sun, J.; Liang, Y.; Marsh, J.F.; Harris, J.M.; Oldroyd, G.E. Abscisic acid coordinates Nod factors and cytokinins signaling during the regulation of nodulation in Medicago truncatula. Plant Cell 2008, 20, 2681–2695. [Google Scholar] [CrossRef] [PubMed]
  93. Madsen, E.B.; Madsen, L.H.; Radutoiu, S.; Olbryt, M.; Rakwalska, M.; Szczyglowski, K.; Sato, S.; Kaneko, T.; Tabata, S.; Sandal, N.; et al. A receptor kinase gene of the LysM type is involved in legume perception of rhizobial signals. Nature 2003, 425, 637–640. [Google Scholar] [CrossRef] [PubMed]
  94. Radutoiu, S.; Madsen, L.H.; Madsen, E.B.; Felle, H.H.; Umehara, Y.; Sato, S.; Nakamura, Y.; Tabata, S.; Sandal, N.; Stougaard, J. Plant recognition of symbiotic bacteria requires two LysM receptor-like kinases. Nature 2003, 425, 585–592. [Google Scholar] [CrossRef] [PubMed]
  95. Kwaharada, Y.; Kelly, S.; Nielsen, M.W.; Hjuler, C.T.; Gysel, K.; Muszynski, A.; Carlson, R.W.; Thygesen, M.B.; Sandal, N.; Asmussen, M.H.; et al. Receptor-mediated exopolysaccharide perception controls bacterial infection. Nature 2015, 523, 308–312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Rose, C.M.; Venkateshwaran, M.; Volkening, J.D.; Grimsrud, P.A.; Maeda, J.; Park, K.; Howes-Podoll, M.; den Os, D.; Yeun, L.H.; Westphall, M.S.; et al. Rapid phosphoproteomic and transcriptomic changes in the rhizobia-legume symbiosis. Mol. Cell. Proteomics 2012, 11, 724–744. [Google Scholar] [CrossRef] [PubMed]
  97. Vankateshwaran, M.; Jayaramar, D.; Chabaud, M.; Genre, A.; Ballon, A.J.; Maeda, J.; Forshey, K.; den Os, D.; Kwiecien, N.W.; Coon, J.J.; et al. A role for the mevalonate pathway in early symbiotic signaling. Proc. Natl. Acad. Sci. USA 2015, 112, 9781–9786. [Google Scholar] [CrossRef] [PubMed]
Genes 09 00498 g001 550
Figure 1. Pi deficiency negatively affects the rhizobia-induced root hairs deformation in common bean. Rhizobia-induced root hairs in common bean seedling growing under optimal-phosphate (Pi) (OP) (a) or low-Pi (LP) conditions (b). Scale bars in panel represent 100 μm. Arrows indicate some examples of the rhizobia-induced root hairs observed in both OP or LP conditions (c) Number of rhizobia-induced root hair deformation in common bean seedlings growing under OP conditions or LP conditions. (d) Number of aberrant rhizobia-induced root hair deformations in common bean seedlings growing under OP conditions or LP conditions. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). Asterisk indicates a significant difference according to one-way analysis of variance (ANOVA) (p-value < 0.001). Data showed was obtained from 20 biological replicates, each one containing 10 roots from different common bean seedlings. (e) Types of aberrant rhizobia-induced root hair deformation observed in common bean seedlings growing under LP conditions. Scale bars in panel represent 50 μm.
Figure 1. Pi deficiency negatively affects the rhizobia-induced root hairs deformation in common bean. Rhizobia-induced root hairs in common bean seedling growing under optimal-phosphate (Pi) (OP) (a) or low-Pi (LP) conditions (b). Scale bars in panel represent 100 μm. Arrows indicate some examples of the rhizobia-induced root hairs observed in both OP or LP conditions (c) Number of rhizobia-induced root hair deformation in common bean seedlings growing under OP conditions or LP conditions. (d) Number of aberrant rhizobia-induced root hair deformations in common bean seedlings growing under OP conditions or LP conditions. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). Asterisk indicates a significant difference according to one-way analysis of variance (ANOVA) (p-value < 0.001). Data showed was obtained from 20 biological replicates, each one containing 10 roots from different common bean seedlings. (e) Types of aberrant rhizobia-induced root hair deformation observed in common bean seedlings growing under LP conditions. Scale bars in panel represent 50 μm.
Genes 09 00498 g001
Genes 09 00498 g002a 550Genes 09 00498 g002b 550
Figure 2. Pi deficiency reduces rhizobia attachment and infection thread formation in common bean. (a–b) Rhizobia attachment in roots from common bean seedlings growing under optimal-Pi (OP) or low-Pi conditions (LP). Seedlings were inoculated with R. tropici::uidA and stained by the β-glucorinidase activity. Scale bars in panel represent 2 mm. Infection thread formation in common bean seedlings growing under OP (c) or LP conditions (d). Arrows indicates some examples of the infection thread observed in OP or LP conditions. Scale bars in panel represent 100 μm. (e) Number of infection thread (IT) per susceptible zone. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). Asterisk indicates a significant difference according to one-way ANOVA (p-value < 0.001). Data showed was obtained from 20 biological replicates, each one containing ten roots from different common bean seedlings.
Figure 2. Pi deficiency reduces rhizobia attachment and infection thread formation in common bean. (a–b) Rhizobia attachment in roots from common bean seedlings growing under optimal-Pi (OP) or low-Pi conditions (LP). Seedlings were inoculated with R. tropici::uidA and stained by the β-glucorinidase activity. Scale bars in panel represent 2 mm. Infection thread formation in common bean seedlings growing under OP (c) or LP conditions (d). Arrows indicates some examples of the infection thread observed in OP or LP conditions. Scale bars in panel represent 100 μm. (e) Number of infection thread (IT) per susceptible zone. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). Asterisk indicates a significant difference according to one-way ANOVA (p-value < 0.001). Data showed was obtained from 20 biological replicates, each one containing ten roots from different common bean seedlings.
Genes 09 00498 g002aGenes 09 00498 g002b
Genes 09 00498 g003 550
Figure 3. Pi deficiency reduces the activation of early symbiotic genes. Rhizobia-triggered expression of PvNSP2 (a), PvNIN (c) and PvFLOT2 (e) genes in common bean seedlings growing under OP or LP conditions and inoculated with mock or rhizobia. Data showed are the fold-changes (e.g., OPR/OPM) obtained from four independent replicates, each one containing susceptible zones from four different plants. Dashed line shows basal expression (e.g., no induction triggered by rhizobia). Expression levels of PvNSP2 (b), PvNIN (d) and PvFLOT2 (f) genes in common bean seedlings grown under OP (dark-gray) or LP (light-gray) conditions and inoculated with mock (−) or rhizobia (+). Data showed was obtained from four independent replicates, each one containing susceptible zones from four different plants. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). One-way ANOVA, followed by a Tukey honest significant difference (HSD) test was performed (p-value < 0.01). Statistical classes sharing a letter are not significantly different.
Figure 3. Pi deficiency reduces the activation of early symbiotic genes. Rhizobia-triggered expression of PvNSP2 (a), PvNIN (c) and PvFLOT2 (e) genes in common bean seedlings growing under OP or LP conditions and inoculated with mock or rhizobia. Data showed are the fold-changes (e.g., OPR/OPM) obtained from four independent replicates, each one containing susceptible zones from four different plants. Dashed line shows basal expression (e.g., no induction triggered by rhizobia). Expression levels of PvNSP2 (b), PvNIN (d) and PvFLOT2 (f) genes in common bean seedlings grown under OP (dark-gray) or LP (light-gray) conditions and inoculated with mock (−) or rhizobia (+). Data showed was obtained from four independent replicates, each one containing susceptible zones from four different plants. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). One-way ANOVA, followed by a Tukey honest significant difference (HSD) test was performed (p-value < 0.01). Statistical classes sharing a letter are not significantly different.
Genes 09 00498 g003
Genes 09 00498 g004 550
Figure 4. Pi deficiency increases the expression of PvRIC1 and PvRIC2 even in the absence of rhizobia. Rhizobia-triggered expression of PvRIC1 (a), and PvRIC2 (c) genes in common bean seedlings growing under OP or LP conditions and inoculated with mock or rhizobia. Data showed are the fold-changes (e.g., OPR/OPM) obtained from four independent replicates, each one containing susceptible zones from four different plants. Dashed line shows basal expression (e.g., no induction triggered by rhizobia). Expression levels of PvRIC1 (b) and PvRIC2 (d) at various times (1, 24 and 48 h) post inoculation with mock (−) or rhizobia (+) in common bean seedlings growing under OP (dark-gray) or LP (light-gray) conditions. Despite several attempts, we were not able to detect reproducible and consistent expression (non-detected: ND) of PvRIC2 at one and 24 h post-inoculation with mock or rhizobia. Data showed was obtained from four independent replicates, each one containing susceptible zones from four different plants. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). One-way ANOVA, followed by a Tukey HSD test was performed (p-value < 0.01). Statistical classes sharing a letter are not significantly different.
Figure 4. Pi deficiency increases the expression of PvRIC1 and PvRIC2 even in the absence of rhizobia. Rhizobia-triggered expression of PvRIC1 (a), and PvRIC2 (c) genes in common bean seedlings growing under OP or LP conditions and inoculated with mock or rhizobia. Data showed are the fold-changes (e.g., OPR/OPM) obtained from four independent replicates, each one containing susceptible zones from four different plants. Dashed line shows basal expression (e.g., no induction triggered by rhizobia). Expression levels of PvRIC1 (b) and PvRIC2 (d) at various times (1, 24 and 48 h) post inoculation with mock (−) or rhizobia (+) in common bean seedlings growing under OP (dark-gray) or LP (light-gray) conditions. Despite several attempts, we were not able to detect reproducible and consistent expression (non-detected: ND) of PvRIC2 at one and 24 h post-inoculation with mock or rhizobia. Data showed was obtained from four independent replicates, each one containing susceptible zones from four different plants. Box plots represent first and third quartile (horizontal box side), minimum and maximum (outside whiskers). One-way ANOVA, followed by a Tukey HSD test was performed (p-value < 0.01). Statistical classes sharing a letter are not significantly different.
Genes 09 00498 g004
Genes 09 00498 g005 550
Figure 5. Expression of rhizobia-responsive genes is affected in Pi-deficient common bean seedlings. Number of over-lapping and non-overlapping up-regulated (a) or down-regulated (b) rhizobia-induced genes among common bean seedlings grown under OP or LP conditions and inoculated with rhizobia (R) or mock (M) during 48 h. Differentially regulated genes in each Pi conditions were identified by Cuffdiff at false discovery rate (FDR) < 0.05, with an additional cutoff of two-fold in a pairwise comparison between treatments (e.g., OPrhizobia/OPmock). Over- and non-overlapping genes were identified after a pairwise comparison between the differentially regulated genes identified in each Pi conditions. A heat map showing the expression levels of up-reguled genes (c) or down-regulated genes (d) from common bean seedlings growing under OP or LP conditions and inoculated with mock or rhizobia. Expression levels were false colored. Genes showing higher expression levels are shown in different scales of red color, whereas genes showing lower expression levels are showed in different scales of green color.
Figure 5. Expression of rhizobia-responsive genes is affected in Pi-deficient common bean seedlings. Number of over-lapping and non-overlapping up-regulated (a) or down-regulated (b) rhizobia-induced genes among common bean seedlings grown under OP or LP conditions and inoculated with rhizobia (R) or mock (M) during 48 h. Differentially regulated genes in each Pi conditions were identified by Cuffdiff at false discovery rate (FDR) < 0.05, with an additional cutoff of two-fold in a pairwise comparison between treatments (e.g., OPrhizobia/OPmock). Over- and non-overlapping genes were identified after a pairwise comparison between the differentially regulated genes identified in each Pi conditions. A heat map showing the expression levels of up-reguled genes (c) or down-regulated genes (d) from common bean seedlings growing under OP or LP conditions and inoculated with mock or rhizobia. Expression levels were false colored. Genes showing higher expression levels are shown in different scales of red color, whereas genes showing lower expression levels are showed in different scales of green color.
Genes 09 00498 g005
Genes 09 00498 g006 550
Figure 6. Functional categories affected in Pi-deficient common seedlings inoculated with rhizobia. MapMan regulation overview map showing differences in transcript levels between common bean seedlings growing under OP (a) or LP (b) conditions and inoculated for 48 h with rhizobia. Fold-changes were false-colored; genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in green color. The complete sets of genes submitted to MapMan analysis are shown in Tables S2 and S3. Cd and Cc: cell division and cell cycle; DR: DNA repair; DS: DNA synthesis; CO: cell organization; Vt and Pt: vesicle transport and protein targeting; Sb and A: stress biotic and abiotic; RS: RNA synthesis; RT: regulation of transcription; RP: RNA processing; Ps and AAa: protein synthesis and amino acid activation; Dev: development; Horm: hormones; Reg: regulation; PM: protein modification; PD: protein degradation; Unknown: unclassified no ontology and unknown; EF: enzyme families; MH: metal-handling; Transp: transport.
Figure 6. Functional categories affected in Pi-deficient common seedlings inoculated with rhizobia. MapMan regulation overview map showing differences in transcript levels between common bean seedlings growing under OP (a) or LP (b) conditions and inoculated for 48 h with rhizobia. Fold-changes were false-colored; genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in green color. The complete sets of genes submitted to MapMan analysis are shown in Tables S2 and S3. Cd and Cc: cell division and cell cycle; DR: DNA repair; DS: DNA synthesis; CO: cell organization; Vt and Pt: vesicle transport and protein targeting; Sb and A: stress biotic and abiotic; RS: RNA synthesis; RT: regulation of transcription; RP: RNA processing; Ps and AAa: protein synthesis and amino acid activation; Dev: development; Horm: hormones; Reg: regulation; PM: protein modification; PD: protein degradation; Unknown: unclassified no ontology and unknown; EF: enzyme families; MH: metal-handling; Transp: transport.
Genes 09 00498 g006
Genes 09 00498 g007 550
Figure 7. Pi deficiency affects the expression of rhizobia-induced auxin- and cytokinin-related genes. Heat-map showing the fold-change of genes involved in the biosynthesis (*), homeostasis (#), signal transduction (..) of auxin (a) or cytokinin (b) in common bean seedlings growing under OP or LP conditions and inoculated for 48 h with rhizobia. Genes regulated by auxins or cytokinin are indicated with (.). Fold-changes were false-colored. Genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in blue color. Genes involved in biosynthesis.
Figure 7. Pi deficiency affects the expression of rhizobia-induced auxin- and cytokinin-related genes. Heat-map showing the fold-change of genes involved in the biosynthesis (*), homeostasis (#), signal transduction (..) of auxin (a) or cytokinin (b) in common bean seedlings growing under OP or LP conditions and inoculated for 48 h with rhizobia. Genes regulated by auxins or cytokinin are indicated with (.). Fold-changes were false-colored. Genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in blue color. Genes involved in biosynthesis.
Genes 09 00498 g007
Genes 09 00498 g008 550
Figure 8. Pi deficiency induce the expression of ethylene-, jasmonate-, and abscisic acid-related genes. Heat-map showing the fold-change of genes involved in the biosynthesis (*), homeostasis (#), signal transduction (..) of ethylene (a), jasmonate (b) or abscisic acid (c) in common bean seedlings growing under OP or LP conditions and inoculated for 48 h with rhizobia. Genes regulated by ethylene, jasmonate or abscisic acid are indicated with (.). Fold-changes were false-colored. Genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in blue color.
Figure 8. Pi deficiency induce the expression of ethylene-, jasmonate-, and abscisic acid-related genes. Heat-map showing the fold-change of genes involved in the biosynthesis (*), homeostasis (#), signal transduction (..) of ethylene (a), jasmonate (b) or abscisic acid (c) in common bean seedlings growing under OP or LP conditions and inoculated for 48 h with rhizobia. Genes regulated by ethylene, jasmonate or abscisic acid are indicated with (.). Fold-changes were false-colored. Genes showing higher expression are shown in different scales of red color, whereas genes showing lower expression were colored in blue color.
Genes 09 00498 g008
Genes 09 00498 g009 550
Figure 9. Proposed model of the effects of Pi deficiency in the early steps of the symbiosis between common bean and rhizobia. Higher induction, lower induction and down-regulation of the expression of different rhizobia-induced genes is indicated in red, orange, and green, respectively. Upon NF perception, a variety of molecular responses occurs. For example, the accumulation of mevalonate at the cytoplasm is required for the activation of the Ca2+ spiking [97]. These fluxes of calcium are further transduced by CCaMK, which in turn activates the transcription factor CYCLOPS. In turn, CYCLOPS is required to activate the transcription factors NSP2/NSP1 and NIN. The coordinated action of these transcription factors activates the expression of different genes (e.g., FLOT2) required for the infection of the root by rhizobia. There is evidence indicating that these molecular events are positively regulated by a delicate balance among the phytohormones auxins (IAA), cytokinins (CK) and ethylene (Eth). These molecular events occur during the first twenty-four hours of interaction between legumes and rhizobia. At this stage, the expression of RIC1 and RIC2 genes, which participate in the AON process, is not fully activated. The coordinated action of these molecular events leads to a successful symbiosis. Based on our data, Pi deficiency might affects the phytohormone balance by activating the expression of genes participating in the biosynthesis of abscisic acid (ABA), ethylene (Eth), and jasmonate (JA). These three hormones negatively affect the Ca2+ spiking [90,91], leading to a reduction in the expression of symbiosis-related genes (e.g., NSP2, NIN, and FLOT2) required for the rhizobia infection process. Additionally, Pi deficiency induces the expression of RIC1 and RIC2, even in the absence of rhizobia. All these defects in the early stages of this symbiosis lead to a reduction in the number of nodules in common bean.
Figure 9. Proposed model of the effects of Pi deficiency in the early steps of the symbiosis between common bean and rhizobia. Higher induction, lower induction and down-regulation of the expression of different rhizobia-induced genes is indicated in red, orange, and green, respectively. Upon NF perception, a variety of molecular responses occurs. For example, the accumulation of mevalonate at the cytoplasm is required for the activation of the Ca2+ spiking [97]. These fluxes of calcium are further transduced by CCaMK, which in turn activates the transcription factor CYCLOPS. In turn, CYCLOPS is required to activate the transcription factors NSP2/NSP1 and NIN. The coordinated action of these transcription factors activates the expression of different genes (e.g., FLOT2) required for the infection of the root by rhizobia. There is evidence indicating that these molecular events are positively regulated by a delicate balance among the phytohormones auxins (IAA), cytokinins (CK) and ethylene (Eth). These molecular events occur during the first twenty-four hours of interaction between legumes and rhizobia. At this stage, the expression of RIC1 and RIC2 genes, which participate in the AON process, is not fully activated. The coordinated action of these molecular events leads to a successful symbiosis. Based on our data, Pi deficiency might affects the phytohormone balance by activating the expression of genes participating in the biosynthesis of abscisic acid (ABA), ethylene (Eth), and jasmonate (JA). These three hormones negatively affect the Ca2+ spiking [90,91], leading to a reduction in the expression of symbiosis-related genes (e.g., NSP2, NIN, and FLOT2) required for the rhizobia infection process. Additionally, Pi deficiency induces the expression of RIC1 and RIC2, even in the absence of rhizobia. All these defects in the early stages of this symbiosis lead to a reduction in the number of nodules in common bean.
Genes 09 00498 g009

Share and Cite

MDPI and ACS Style

Isidra-Arellano, M.C.; Reyero-Saavedra, M.D.R.; Sánchez-Correa, M.D.S.; Pingault, L.; Sen, S.; Joshi, T.; Girard, L.; Castro-Guerrero, N.A.; Mendoza-Cozatl, D.G.; Libault, M.; et al. Phosphate Deficiency Negatively Affects Early Steps of the Symbiosis between Common Bean and Rhizobia. Genes 2018, 9, 498. https://doi.org/10.3390/genes9100498

AMA Style

Isidra-Arellano MC, Reyero-Saavedra MDR, Sánchez-Correa MDS, Pingault L, Sen S, Joshi T, Girard L, Castro-Guerrero NA, Mendoza-Cozatl DG, Libault M, et al. Phosphate Deficiency Negatively Affects Early Steps of the Symbiosis between Common Bean and Rhizobia. Genes. 2018; 9(10):498. https://doi.org/10.3390/genes9100498

Chicago/Turabian Style

Isidra-Arellano, Mariel C., María Del Rocio Reyero-Saavedra, Maria Del Socorro Sánchez-Correa, Lise Pingault, Sidharth Sen, Trupti Joshi, Lourdes Girard, Norma A. Castro-Guerrero, David G. Mendoza-Cozatl, Marc Libault, and et al. 2018. "Phosphate Deficiency Negatively Affects Early Steps of the Symbiosis between Common Bean and Rhizobia" Genes 9, no. 10: 498. https://doi.org/10.3390/genes9100498

APA Style

Isidra-Arellano, M. C., Reyero-Saavedra, M. D. R., Sánchez-Correa, M. D. S., Pingault, L., Sen, S., Joshi, T., Girard, L., Castro-Guerrero, N. A., Mendoza-Cozatl, D. G., Libault, M., & Valdés-López, O. (2018). Phosphate Deficiency Negatively Affects Early Steps of the Symbiosis between Common Bean and Rhizobia. Genes, 9(10), 498. https://doi.org/10.3390/genes9100498

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop