# Plant Nitrate Reductases Regulate Nitric Oxide Production and Nitrogen-Fixing Metabolism During the Medicago truncatula–Sinorhizobium meliloti Symbiosis

Antoine Berger, Alexandre Boscari, Natasha Horta Araújo, Mickaël Maucourt, Mohamed Hanchi, Stéphane Bernillon, Dominique Rolin, Alain Puppo, Renaud Brouquisse

► To cite this version:

Antoine Berger, Alexandre Boscari, Natasha Horta Araújo, Mickaël Maucourt, Mohamed Hanchi, et al.. Plant Nitrate Reductases Regulate Nitric Oxide Production and Nitrogen-Fixing Metabolism During the Medicago truncatula–Sinorhizobium meliloti Symbiosis. Frontiers in Plant Science, 2020, 11, pp.1313. 10.3389/fpls.2020.01313. hal-02983822

HAL Id: hal-02983822

https://hal.inrae.fr/hal-02983822

Submitted on 30 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Plant Nitrate Reductases Regulate Nitric Oxide Production and Nitrogen-Fixing Metabolism During the *Medicago truncatula*–*Sinorhizobium meliloti* Symbiosis

## OPEN ACCESS

**Edited by:**

Juan Carlos Begara-Morales,
University of Jaén, Spain

**Reviewed by:**

Clare Gough,
Centre National de la Recherche
Scientifique (CNRS), France
Mounira Chaki,
University of Jaén, Spain

**\*Correspondence:**

Renaud Brouquisse
renaud.brouquisse@inrae.fr

**Specialty section:**

This article was submitted to
Plant Pathogen Interactions,
a section of the journal
Frontiers in Plant Science

**Received:** 14 February 2020

**Accepted:** 11 August 2020

**Published:** 04 September 2020

**Citation:**

Berger A, Boscari A, Horta Araújo N, Maucourt M, Hanchi M, Bernillon S, Rolin D, Puppo A and Brouquisse R (2020) Plant Nitrate Reductases Regulate Nitric Oxide Production and Nitrogen-Fixing Metabolism During the Medicago truncatula–Sinorhizobium meliloti Symbiosis. *Front. Plant Sci.* 11:1313. doi: 10.3389/fpls.2020.01313

Antoine Berger¹,², Alexandre Boscari¹, Natasha Horta Araújo¹, Mickaël Maucourt³, Mohamed Hanchi¹, Stéphane Bernillon⁴, Dominique Rolin³,⁴, Alain Puppo¹ and Renaud Brouquisse¹*

¹ Institut Sophia Agrobiotech, UMR INRAE 1355, Université Côte d'Azur, CNRS, Sophia Antipolis, France, ² Department of Horticultural Science, University of Florida, Gainesville, FL, United States, ³ Univ. Bordeaux INRAE, UMR Biologie du Fruit et Pathologie, Villenave d'Ornon, France, ⁴ PMB-Metabolome, INRAE, Bordeaux Metabolome Facility, Villenave d'Ornon, France

Nitrate reductase (NR) is the first enzyme of the nitrogen reduction pathway in plants, leading to the production of ammonia. However, in the nitrogen-fixing symbiosis between legumes and rhizobia, atmospheric nitrogen ($N_2$) is directly reduced to ammonia by the bacterial nitrogenase, which questions the role of NR in symbiosis. Next to that, NR is the best-characterized source of nitric oxide (NO) in plants, and NO is known to be produced during the symbiosis. In the present study, we first surveyed the three NR genes (*MtNR1*, *MtNR2*, and *MtNR3*) present in the *Medicago truncatula* genome and addressed their expression, activity, and potential involvement in NO production during the symbiosis between *M. truncatula* and *Sinorhizobium meliloti*. Our results show that *MtNR1* and *MtNR2* gene expression and activity are correlated with NO production throughout the symbiotic process and that *MtNR1* is particularly involved in NO production in mature nodules. Moreover, NRs are involved together with the mitochondrial electron transfer chain in NO production throughout the symbiotic process and energy regeneration in $N_2$-fixing nodules. Using an *in vivo* NMR spectrometric approach, we show that, in mature nodules, NRs participate also in the regulation of energy state, cytosolic pH, carbon and nitrogen metabolism under both normoxia and hypoxia. These data point to the importance of NR activity for the $N_2$-fixing symbiosis and provide a first explanation of its role in this process.

**Keywords:** hypoxia, legumes, *Medicago truncatula*, nitric oxide, nitrogen-fixing symbiosis, nitrate reductase, nodules
# INTRODUCTION

In plants, yeasts, algae, and fungi, nitrate reductase (NR) is a key enzyme of the nitrogen (N) reduction and assimilation pathway. It catalyzes the reduction of nitrate ($NO_3^−$) to nitrite ($NO_2^−$), which is itself reduced to ammonia ($NH_4^+$) by nitrite reductase (NiR), before being assimilated into the amino acids and the nitrogen compounds of the cell (Campbell, 1999). Sixty-five million years ago, the legume family developed a beneficial mutual relationship with soil bacteria, the Rhizobia, which directly reduce atmospheric nitrogen ($N_2$) to ammonia ($NH_4^+$) through the activity of the nitrogenase under nitrogen deficiency situations (Wang et al., 2018). In exchange for $NH_4^+$, plants supply “board and lodging” for the bacteria, providing them with an ecological niche for their development and carbon nutrients for their functioning, in neoformed organs called nodules (Udvardi and Day, 1997; Terpolilli et al., 2012). Thus, in legume nodules, the bacterial nitrogenase substitutes for the NR-NiR pathway to produce $NH_4^+$ and supply the plant with reduced nitrogen. However, many studies have reported high NR expression and activity in symbiotic nodules (see Streeter, 1985a; Streeter, 1985b; Arrese-Igor et al., 1990; Silveira et al., 2001; Kato et al., 2003; Sanchez et al., 2010; Horchani et al., 2011, and references therein) and the question arose—and still arises—of what the NR can be used in the $N_2$-fixing symbiosis, where nitrogen reduction is ensured by nitrogenase.

The involvement of NR activity in NO production has been evidenced in many plant organs and tissues over the past 20 years (Dean and Harper, 1988; Yamasaki and Sakihama, 2000; Rockel et al., 2002; Sakihama et al., 2002; Gupta et al., 2005; Planchet et al., 2005; Kato et al., 2010; Kolbert et al., 2010). Until recently, NR was considered to produce directly NO via the reduction of $NO_2^−$, but another indirect mechanism of NO synthesis involving NR has been proposed (Chamizo-Ampudia et al., 2016; Chamizo-Ampudia et al., 2017). In this mechanism, through its diaphorase activity, NR transfers electron from NAD(P)H to a NO-forming nitrite reductase (NOFNIR) that catalyzes the reduction of $NO_2^−$ to NO. Although the interaction between NR and NOFNIR has been argued in eukaryotic algae, the proof of concept in higher plant is still not demonstrated. NO is a reactive free radical gaseous molecule with a broad spectrum of regulatory functions in plant growth and development, and in response to abiotic and biotic factors (Kolbert et al., 2019). NO is particularly involved in the legume-Rhizobium symbiotic interactions (Hichri et al., 2015; Berger et al., 2019).

During the first hours after inoculation with the symbiotic partner, NO was observed in the roots of *Lotus japonicus*, *Medicago sativa*, and *Medicago truncatula* (Nagata et al., 2008; Fukudome et al., 2016; Hichri et al., 2016). Its production was also detected during the infection process along the infection thread and in the dividing cells of the *M. truncatula* nodule primordium (del Giudice et al., 2011). Similar results were observed in other Medicago species (Pii et al., 2007). In *M. truncatula* mature nodules, NO has been shown to accumulate particularly in the $N_2$-fixing zone (Baudouin et al., 2006; Hichri et al., 2016), and at the onset of nodule senescence a NO

production was reported at the junction of the $N_2$-fixing and senescence zones (Cam et al., 2012). A recent study with *M. truncatula* showed that NO is produced throughout the whole symbiotic process, from infection with *Sinorhizobium meliloti* up to, at least, 8 weeks post-inoculation (wpi), exhibiting production peaks during the first hours of the symbiotic interaction, during early development of the nodule and when the nodule becomes mature (Berger et al., 2020). These observations suggest that NO performs specific signaling and/or metabolic functions during symbiosis. Indeed, two transcriptomic analyses led to the identification of NO-responsive genes either in 4 days post-inoculation (dpi) roots (Boscari et al., 2013) or in developing and mature nodules (Ferrarini et al., 2008). More than 400 plant genes are NO-regulated during the symbiotic process, including genes involved in nodule development normally induced by the symbiont, suggesting that NO participate in signal transduction in the plant-microorganism interaction (Ferrarini et al., 2008; Boscari et al., 2013). On the bacterial partner side, NO has also been shown to regulate a hundred genes, most of them being similarly regulated under microoxic conditions (Bobik et al., 2006; Meilhoc et al., 2010). The biological activity of NO is particularly mediated through redox-dependent protein modifications such as S-nitrosation, tyrosine nitration and metal nitrosylation (Besson-Bard et al., 2008; Hancock, 2019). Several key proteins involved in nodule primary metabolism or stress response were reported to be S-nitrosated, indicating a crucial role of NO in the energy, carbon, and nitrogen metabolism (Puppo et al., 2013). Among these proteins, enzymes such as glutathione peroxidase (Castella et al., 2017), glutamine synthetase (Melo et al., 2011; Seabra and Carvalho, 2015), and leghemoglobins (Mathieu et al., 1998; Navascues et al., 2012; Sainz et al., 2015; Becana et al., 2020) have been shown to be differently regulated by various NO-dependent modifications. Finally, NO has also been shown to play a metabolic function in the maintenance of energy status under hypoxic conditions, such as that prevailing in microoxic nodules (Igamberdiev and Hill, 2004; Igamberdiev and Hill, 2009). Indeed, NO is involved in a respiratory cycle, called Phytoglobin-NO respiration (PNR), allowing the regeneration of ATP under low oxygen concentrations. PNR is divided into four steps including: 1) $NO_3^−$ reduction to $NO_2^−$ by cytosolic NR, 2) $NO_2^−$ transport from the cytosol to the mitochondrial matrix, 3) $NO_2^−$ reduction to NO by the mitochondrial electron transfer chain, and 4) NO diffusion to the cytosol and oxidation to $NO_3^−$ by Phytoglobins. Evidence indicates that PNR potentially functions and participates in the regeneration of ATP in $N_2$-fixing nodules (Horchani et al., 2011; Berger et al., 2020). Thus, the importance and diversity of NO functions in symbiosis, whether through the regulation of gene expression, the modulation of enzyme activity, or its involvement in energy metabolism, require that its production should be timely and tightly regulated during the symbiotic process. Considering that the NR pathway is the main NO production pathway in plants and that its activity is high in nodules, we hypothesized that NR is significantly involved in the production of NO during the symbiotic process.
The aim of this work was to test this hypothesis. To this end, we addressed the expression and activity of the three NR genes present in the *M. truncatula* genome, from the first hours of symbiotic interaction up to 8 wpi, at the onset of nodule senescence. Then, we investigated the impact of NR activity on NO production and energy regeneration during nodule development. Last, we studied its role in carbon and nitrogen metabolism in mature nodules. Based on our data, we discuss the potential roles of NR on NO homeostasis during the symbiotic process and more generally on the N₂-fixing process.

# MATERIALS AND METHODS

## Plants Growth and Inoculation Conditions

*Medicago truncatula* (cv Jemalong A17) were scarified, sterilized, and germinated as in del Giudice et al. (2011). Seedlings were cultivated and inoculated with *Sinorhizobium meliloti* 2011 strain either in Petri dishes as in del Giudice et al. (2011), or in planters as in Horchani et al. (2011). A basic intake of 0.2 mM KNO₃ is provided to plants on Petri dishes and planters. Cultures in Petri dishes were used for short-term experiments up to 14 days post-inoculation (dpi), while those in planters were used for long-term experiments up to 8 weeks post-inoculation (wpi). Roots and/or nodules were harvested at various times of the kinetics. For short-term experiments, 2 cm-long root segments corresponding to the infection zone (del Giudice et al., 2011) were harvested for gene expression and NO production. For long-term experiments only nodules were used.

## NR Sequences Acquisition and Analyses

Protein sequences of NR were obtained from three genomic and protein databases: NCBI (www.blast.ncbi.nlm.nih.gov), Phytozome (www.phytozome.jgi.doe.gov), and Uniprot (www.uniprot.org). Local or multiple alignment search tools (BLAST from NCBI and Phytozome; Water and ClustalW from the European Bioinformatic Institute platform EBI) were used to analyze the sequence's quality, length, and uniqueness. Sequences kept were listed in **Supplementary Table S1**, where a nomenclature code was assigned to each protein in order to simplify data-reading in the phylogenetic tree. The code is made of the name of the species represented, the name allotted to the sequence on the literature or NCBI, and the number of amino acids composing each protein.

## Construction of the Phylogenetic Tree

NR sequences were aligned using Muscle algorithm (Madeira et al., 2019), and their evolutionary history was inferred by Maximum Likelihood (ML) and the JTT matrix-based model (Jones et al., 1992) using MegaX software (Kumar et al., 2018). Initial trees calculated for the heuristic search were acquired automatically by applying Neighbor-Join (NJ) and BioNJ algorithms to a pair-wise distance matrix approximated using the JTT model. The number of distinct/identical residues between each pair of sequences in the multiple alignment was therefore calculated to construct the matrix, and the topology

with a higher log likelihood value was selected. The tree was inferred by Nearest-Neighbor-Interchange (NIN) heuristic method, and to assess the reliability of the inference, the bootstrapping method was applied. The consensus tree deduced from 1,000 bootstrap repetitions was retained, and branches corresponding to divisions reproduced in less than 50% of the bootstrap repetitions were not considered. An outgroup, *Chlamydomonas reinhardtii*'s NR NIT1, was used to root the tree.

## Construction of a Binary Vector for Hairy Root Transformation

For promoter transcriptional fusions, fragments of 1,647, 1,700, and 1,554 bp upstream of the start codon were amplified by PCR using the primers indicated in **Table S2** for MtNR1, MtNR2, and MtNR3 respectively. Each PCR fragment was first cloned into the pDONR207 donor vector and then into the plant expression vector pKGWFS7 (Karimi et al., 2002) using Gateway technology (Invitrogen, http://www.invitrogen.com). For the RNA interference (RNAi) construct, 432- and 441-bp fragments of MtNR1 (MtrunA17Chr3g0115151) and MtNR2 (MtrunA17Chr5g0424491) genes were amplified via polymerase chain reaction (PCR) with specific primers (**Table S2**). PCR products were independently ligated into pGEM-T easy vector (Promega) and subsequently subcloned into pENTR4 vectors in *BamHI*-KpnI restriction sites for MtNR1 and EcoRI and KpnI restriction sites for MtNR2. The pENTR4 vector carrying the MtNR1 or the MtNR2 fragment was recombined with pK7GWIWG5D(II) vector (Horchani et al., 2011) using the LR clonase enzyme mix (catalog no. 11791-019; Invitrogen) to create the RNAi expression vectors. As control, transgenic roots transformed with empty pK7GWIWG5D(II) vector and selected on the base of green fluorescent protein (GFP) marker expression were used. Constructs were checked by sequencing, introduced by electroporation into *Agrobacterium rhizogenes* strain ARqual, and used for *M. truncatula* root transformation as described by Boisson-Dernier et al. (2001).

## Measurement of NO Production

NO detection was performed as in Horchani et al. (2011) using the 4,5-diaminofluorescein probe (DAF-2, Sigma-Aldrich) with the following changes. Either nodules (20–30 mg fresh weight) or root segments (50–100 mg fresh weight) were incubated in the dark at 23°C in 1 ml of detection buffer (10 mM Tris-HCl pH 7.4, 10 mM KCl) in the presence of 10 µM DAF-2. As a control, NO production was measured in the same experimental system through the use of the Cu(II) fluorescein (CuFL) fluorescent probe (Strem Chemicals) instead of DAF-2 in the detection buffer as described in Horchani et al. (2011). Similar results were obtained with both probes. The production of NO was measured with a spectrofluorimeter-luminometer (Xenius, SAFAS, Monaco). Inhibitors are used as described in Horchani et al. (2011). The inhibitors are added to the reaction medium for the determination of NO at the concentration of 1 mM tungstate (Tg), 1 mM allopurinol, 1 mM propyl gallate, and 300 µM potassium cyanide (KCN). NO production was initiated 1 h after the addition of inhibitors. Three independent biological
replicates have been performed with three technical replicates per biological assay.

## Nitrogen-Fixing Capacity Measurement

Nitrogenase activity of the nodules was determined *in vivo* by measuring acetylene reducing activity (ARA, Hardy et al., 1968). Nodulated roots were harvested and incubated at 30°C for 1 h in rubber-capped tubes containing 10% acetylene atmosphere. Ethylene concentrations were determined by gas chromatography (Agilent GC 6890N, Agilent Technologies) equipped with a GS-Alumina separating capillary column. Three independent biological replicates have been performed with five technical replicates per biological assay.

## Phosphorus NMR

For each experiment, 0.9 to 1.1 g fresh weight of 4 wpi-old nodules (around 1,400 to 1,700 nodules) were harvested and incubated at ambient temperature in an aerated perfusion medium containing 1 mM KH₂PO₄, 1 mM MgSO₄·7H₂O, 0.25 mM K₂SO₄, 0.5 mM CaCl₂, 10 mM MES/KOH, pH 6.0, and 25 mM glucose. At the end of the preparation period of approximately 3 h, the nodules were placed between two filters into a 10-mm tightly closed NMR tube, part of a homebuilt perfusion system. The latter, evolved from experimental device described previously (Roby et al., 1987), allows circulation of the perfusion medium controlled in solute composition, temperature, and pH through the living nodule sample. The partial oxygen pressure in the perfusion medium was established by bubbling mixtures of oxygen and nitrogen (either 21:79 or 1:99% O₂:N₂) into the medium reservoir. At various time, effectors were added into the perfusion medium. Each series of *in vivo* biological experiment has been performed at least five times.

³¹P NMR spectra were acquired at 202.47 MHz using a 500.16 MHz NMR spectrometer (Avance III, Bruker). For *in vivo* experiments, ³¹P NMR spectra were recorded for 36 min using a 10 mm ATMA broadband observe probe. A solution of 500 mM Hexamethylphosphoramide (HMPA, ref H3380, Sigma) contained in a concentric capillary provides the chemical shifts and intensity references for the ³¹P NMR spectra. 3,072 scans of 16 K data points were acquired with a 60° pulse angle, a spectral width of 14,204 Hz, acquisition time of 0.58 s and recycle delay of 0.1 s. Preliminary data processing was carried out with TOPSPIN 3.0 software (Bruker Biospin, Karlsruhe, Germany). Each Free Induction Decay (FID) was Fourier transformed (10 Hz line broadening), manually phased and baseline corrected. The resulting spectra were aligned by setting the HMPA signal to 30.73 ppm.

The resonance assignments were based on chemical shifts. ³¹P chemical shifts were determined according to (Rolin et al., 1989). Subcellular pH was estimated by the use of a standard reference curve of pH as a function of chemical shift, which was obtained according to the method of Roberts et al. (1980).

## Measurement of NR Activity

Tissue samples are ground with mortar and pestle in liquid nitrogen. The total proteins are extracted from 100 mg of powder

using the following extraction buffer: 25 mM Tris HCl, pH 8.5, 1 mM EDTA, 20 µM FAD, 0.04% Triton X100, 10 µM NaMO₄, 1 mM DTT, 20 µM L-transepoxysuccinyl-leucylamido-[4-guanidino]butane (E64), 2 mM phenylmethylsulfonyl fluoride (PMSF). The extracts are centrifuged (15,000 g, 15 min). The NR activity is measured by quantifying the NO₂⁻ produced in the reaction mixture containing: the enzyme extract in 0.2 M HEPES, pH7.0, 15 mM KNO₃, and 250 µM of NADH (Miranda et al., 2001). The reaction is stopped after 30 min by boiling the samples for 3 min at 100°C. The nitrite produced is measured using Griess reagent (1% sulphanilamide in 1 M HCl and 0.01% NEDD [N-1-naphthylethylenediamine dihydrochloride] in water) and measured at 540 nm. Soluble proteins are assayed according to the method of Bradford (1976). For each series of experiments, at least three independent biological replicates have been performed with three technical replicates per biological assay.

## RNA Isolation, Reverse Transcription and Gene Expressions

RNAs were isolated from 100 mg of frozen material ground in liquid N₂ using the RNAzol following the manufacturer's recommendations (Sigma-Aldrich). RNA quality was checked, and DNase treatment was carried out before the synthesis by GoScript reverse transcriptase (Promega) of the cDNAs. The RT-qPCR was made with Go-Taq qPCR master Mix kit according to manufacturer's instructions (Promega). RT-qPCR data analysis was carried out using RqPCRBase, an R package working on R computing environment for analysis of quantitative real-time PCR data (Hilliou and Tran, 2013). The expression of the different genes was normalized against two housekeeping genes *Mtc27* (Van de Velde et al., 2006) and *Mta38* (del Giudice et al., 2011). The reference value '1' was attributed to the first time when the cycle threshold (Ct) of the analyzed gene was significantly detectable. The data of comparative expression levels between genes are given on a logarithmic scale expressed as 40 - ΔC<sub>T</sub>, where ΔC<sub>T</sub> is the difference in qRT-PCR threshold cycle number between the respective gene and the reference gene; the number 40 was chosen because PCR run stops after 40 cycles (Bari et al., 2006; Truong et al., 2015; Berger et al., 2020). RT-qPCR analyses were carried out in triplicate using the primers reported in **Table S2**. Three independent biological replicates have been performed.

## Enzymatic and HPLC Metabolite Analyses

Nodule metabolites (amino acids, organic acids, soluble sugars) were extracted by the alcoholic extraction method and resuspended in water as described in Brouquisse et al. (1991). Sucrose was determined enzymatically (Velterop and Vos, 2001) at 340 nm using a microplate spectrophotometer. Succinate and malate were determined by anion exchange HPLC (Dionex) with conductivity detection (Moing et al., 1998). Free amino-acids were analyzed by HPLC using the AccQ.Tag method from Waters (Milford, MA) with fluorescence detection (Moing et al., 1998). Three independent biological replicates have been performed with three technical replicates per biological assay.
# Extraction and Measurement of Nodule Adenine Nucleotides

Adenine nucleotides were extracted and measured as in Horchani et al. (2011). Adenine nucleotides were measured in a Xenius spectrofluorimeter-luminometer using the ATPlite one-step assay system (Perkin-Elmer) according to the manufacturer's instructions. Three independent biological replicates have been performed with three technical replicates per biological assay.

# β-Glucuronidase Detection In Planta

For β-glucuronidase detection, nodulated roots of composite plants were incubated at -20°C for 1 h in a mixture "acetone: phosphate buffer (Na₂HPO₄/NaH₂PO₄ 0.1 M, pH 7.4)", 90%:10% (v/v), then washed twice with the phosphate buffer and incubated 3 h to overnight at room temperature in the dark in phosphate buffer containing potassium ferricyanide (0.5 mM) and X-gluc (0;5 ng·ml⁻¹). Nodulated roots or 80-µm-thick vibroslices, obtained with a HM560V Vibratome (Microm, http://www.microm.de) after embedding plant material in 4.5% low-melting-point agarose, were visualized with a Zeiss Axioplan II microscope (Carl Zeiss, http://www.zeiss.com) using dark-field optics. For each transformation experiment, 40 independent roots or nodules were transformed and individually analyzed for β-glucuronidase activity.

# Statistical Analyses

Statistical analyses were performed using Student t-test or one-way analysis of variance (ANOVA) followed by a Fisher test. Data were considered as significantly different when p < 0.05.

# RESULTS

## NR Gene Expression During the Symbiosis Process

Research in three genomic and protein databases: NCBI (www.blast.ncbi.nlm.nih.gov), Phytozome (www.phytozome.jgi.doe.gov) and Uniprot (www.uniprot.org) revealed that all the plant genomes studied code for one to six NR sequences, with 76% of the species having between one and two NR (Table S1). Among the genomes of the Fabales group in Table S1, only three species seem to have more than two NR sequences: *Glycine max* and *Glycine soja*, with five and four sequences respectively and *M. truncatula* with three sequences. The three *M. truncatula* NR genes (Puppo et al., 2013; Roux et al., 2014), named *MtNR1*, *MtNR2*, and *MtNR3*, respectively encode proteins of 902, 884, and 876 amino acids (Figure S1A). The three NR genes exhibit strong similarities (>70%) to each other and exhibit similar structure with four exons and three introns of different sizes (Figure S1B). Both *MtNR1* and *MtNR3* are closely present in chromosome 3, while *MtNR2* is present in chromosome 5. To study the phylogenetic relationships operating on the NR sequences, a phylogenetic tree using NR from the Fabales family and different model plants was inferred (Figure 1).

Forty seven non-redundant protein sequences from 21 species and belonging to the Solanales, Fabales, Brassicales, Malvales, Vitales, and Poales families were used (Table S1). The most striking result is the separation of the Fabales sequences into three distinct Hypothetical Taxonomic Units (HTU), named HTU1, HTU2, and HTU3 (Figure 1). *MtNR2* is in the largest unit HTU1, whereas *MtNR1* and *MtNR3* are both in HTU3. For species with more than one sequence, we observe a distribution of these sequences over two HTUs.

As reported in the Affymetrix data (Figure S2), both *MtNR1* and *MtNR2* are significantly expressed in the roots and nodules, but *MtNR3* is specifically expressed in the nodules. More precisely, the histochemical detection of GUS activity under the control of NR promoters at 4 dpi shows that the three NR genes are expressed in young developing nodules, i.e. in the controlled area (a block of cells constituted by pericycle, endodermis and dividing cortical cells, Xiao et al., 2014) of nodule primordium (Figures 2A, C, E). In 14 dpi-old nodules, *MtNR1* expression is detected in the whole nodule (Figure 2B). *NR2* expression is also detected in the three zones but also at the periphery of the nodule in the vascular bundle (Figure 2D). In the case of *MtNR3*, GUS staining appears mainly at the level of zones I and II and very slightly in zone III (Figure 2F). These results are consistent with Symbimics data (Roux et al., 2014; Table 1) which report a localization of the expression of the three NR genes close to that which we observe in 14 dpi-old nodules.

Then, we investigated their expression kinetics throughout the symbiotic process. To this end, we used two types of *M. truncatula* cultures: a short-term culture from 0 to 14 dpi and a long-term culture from 0 to 8 wpi. Both *MtNR1* and *MtNR2* are expressed at a significant level in non-inoculated roots (Figures 3 and 4). After inoculation with *Sinorhizobium meliloti*, *MtNR1* expression significantly increases (30- to 45-fold) and exhibits three expression peaks at 10 hpi, 4 dpi, and 5 wpi (Figures 3A, B). Similarly, *MtNR2* expression increases with two peaks at 10 hpi and 4 dpi and then reaches a plateau between 3 and 8 wpi (Figures 3C, D). *MtNR3* expression is detected only from 4 dpi, increases to reach a plateau between 3 and 5 wpi, and then strongly increased (up to 300 times) after 6 wpi at the nodule senescence (Figures 3E, F). The expression level of NR as compared to each other, before inoculation and at four time-points in the symbiosis, is reported in Figure 4. Several features emerged from this analysis. 1) *MtNR2* is the most highly expressed NR gene (30 to 100 times more than *NR1*) throughout the whole process. 2) The three NR genes are significantly expressed in N₂-fixing nodules. 3) Whereas the expression of *MtNR1* decreases at the onset of nodule senescence, the expression of *MtNR3* clearly increases and exceeds that of *MtNR1* in 8 wpi nodules. By comparison with the other enzymes potentially involved in the synthesis of NO such as abscissic aldehyde oxidase (AAO), indole aldehyde oxidase (IAO), NO forming nitrite reductase (NOFNI), sulfite oxidase (SOX), and xanthine dehydrogenase (XDH), it may be noted that the expression of *MtNR1* and *MtNR2* is significantly stronger than that of the other genes (Figure S2).
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Phylogenetic tree of NR protein sequences. 47 nitrate reductase sequences from model plants or belonging to the Fabaceae family were used in this analysis. The tree was inferred by maximum likelihood (ML) with MegaX (Kumar et al., 2018). Percentages displayed next to each branch represent the number of tree replicates in which the associated taxa were assembled together when performing 1,000 bootstraps. The lengths of the branches do not represent a phylogenetic distance.</figcaption></figure>

## NR Activity and NO Level During the Symbiotic Process

To assess the contribution of NR in NO production, the total NR activity was analyzed and compared to the NO production level during the symbiotic process (Figure 5). NR activity in non-inoculated roots is close to 4.2 ± 0.9 nmol.min⁻¹.g⁻¹FW. Following inoculation, NR activity exhibits a significant and reproducible 60–70% decrease within 4 hpi before returning close to its initial value at 10 hpi. After a new decrease at 24 hpi, NR activity strongly increases to 26.5 ± 1.1 nmol.min⁻¹.g⁻¹FW at 4 dpi, strongly decreases between 9 and 14 dpi, and then increases again to peak at 3–5 wpi. After the onset of nodule senescence, NR activity decreases to reach 5.1 ± 0.1 nmol.min⁻¹.g⁻¹ FW at 8 wpi. As far as it is concerned, NO production level shows three transient production peaks, at 10 hpi, 4 dpi, and 3–4 wpi (Figures 5 and S3; Berger et al., 2020). As a whole, the pattern of NR activity follows that of MtNR1 and 2 gene expression (Figure 3) and clearly shows a parallel with the production of NO (Figure 5). It particularly fits with the

expression of MtNR1 and MtNR2 during the three first weeks of the symbiosis and with that of MtNR1 between 3 and 8 wpi, suggesting that NR1 is playing a major role in NO production. To check this hypothesis, a RNAi strategy was used. *M. truncatula* RNAi on the MtNR1, MtNR2, and double MtNR1-2 genes were constructed under the control of the zone III-specific promoter NCR001 (Mergaert et al., 2003; Horchani et al., 2011). Four wpi-old nodules were collected and analyzed for NR activity and NO production. As compared to the control, NR activity is decreased by 47 and 56%, respectively, in the RNAi::NR1 and NR2 nodules (Figure 6A), but the double RNAi::NR1-2 does not make it possible to further reduce the total NR activity. Regarding NO production, it drops by 45% in RNAi::NR1 nodules, but only by 18% in RNAi::NR2 nodules (Figure 6B). The decrease in NO production in the nodules of double RNAi::NR1-2 is of the same order of magnitude as that observed in RNAi::NR1 nodules. These results clearly show that the decrease in NR activity is accompanied by a fall in NO production and that this fall is particularly related to NR1 activity rather than NR2.
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Histochemical localization of *MtNRs* expression in *Medicago truncatula* roots. Localization of GUS activity in transgenic *M. truncatula* roots expressing the gusA reporter gene under the control of a 1.65 Kb *MtNR1* promoter fragment (**A, B**), of a 1.7 kb *MtNR2* promoter (**C, D**) and of a 1.55 kb *pMtNR3* promoter (**E, F**). Whole root segment 4 dpi with *S. meliloti* (**A, C, E**) and longitudinal section of a 2 wpi old nodules (**B, D, F**) were stained for 3, 5, or 16 h with X-gluc for the GUS activity for *pMtNR2*, *pMtNR1*, or *pMtNR3* respectively. Zones I, II, and III of the nodule are represented in red in picture F. p, nodule primordium; Vb, vascular bundles; c, cortex. Scale bars, 50 µm for (**A, C, E**); 100 µm for (**B, D, F**).</figcaption></figure>

## Involvement of NR Activity in NO Production and Energy Regeneration During Nodule Development

The concomitance of the peaks of the NR expression/activity and those of NO production suggests that NRs are involved in the production of NO. Thus, using an inhibitor approach, we investigated the participation of NR and other potential NO sources such as xanthine dehydrogenase and mitochondrial electron transport chain (ETC) in NO production in 10 hpi and 4 dpi roots and in 4 wpi nodules. As reported in **Figure 7**, NO production is 95% inhibited by KCN, used as a negative control. Allopurinol (AP), a specific inhibitor of xanthine dehydrogenase (XDH), moderately inhibits NO production (by 28%) in 10 hpi roots, whereas it is without effect in 4 dpi roots and 4 wpi nodules. In both 10 hpi and 4 dpi roots, propyl gallate (PG), an inhibitor of

the mitochondrial alternative oxidase (AOX), inhibits NO production by 70-90%, indicating that the mitochondria are involved in this reaction. However, PG is ineffective in 4 wpi nodules. It can be noted that PG is also an inhibitor of polyphenol oxidases (Lin et al., 2013), but that the latter are not known to be directly or indirectly involved in the production of NO. The effects of PG on NO production can therefore be attributed to AOX inhibition. Interestingly, tungstate (Tg), a NR inhibitor, inhibits NO production by 88, 92, and 60% respectively in 10 hpi roots, 4 dpi roots, and 4 wpi nodules, and this inhibition is either partially (10 hpi and 4 dpi roots) or totally (4 wpi nodules) relieved by the addition of nitrite, the reaction product of NR (Figure 7). These data suggest that NR and ETC are involved in NO production.

To further investigate the potential involvement of NR and NO in energy metabolism, we analyzed the energy state (i.e. the
<table><caption>TABLE 1 | Access code and Symbimics expression of *Medicago truncatula* NR and NO-producing enzyme genes.</caption><thead><tr><th rowspan="2">Genes</th><th rowspan="2">Code Affymetrix</th><th rowspan="2">Code gene Mt4.0</th><th rowspan="2">Code gene Mt5.0</th><th rowspan="2">Mt20120830-LIPM</th><th colspan="5">DESEQ MEAN</th></tr><tr><th>FI</th><th>FIID</th><th>FIIP</th><th>IZ</th><th>ZIII</th></tr></thead><tbody><tr><td>NR1</td><td>Mtr.42446.1.S1_at</td><td>Medtr3g073180</td><td>MtrunA17Chr3g0115151</td><td>Mt0006_00730</td><td>20.7</td><td>89.1</td><td>21.3</td><td>193.7</td><td>693.5</td></tr><tr><td>NR2</td><td>Mtr.10604.1.S1_at</td><td>Medtr5g059820</td><td>MtrunA17Chr5g0424491</td><td>Mt0008_10301</td><td>156.6</td><td>166.9</td><td>266.7</td><td>712.9</td><td>843</td></tr><tr><td>NR3</td><td>Mtr.31448.1.S1_at</td><td>Medtr3g073150</td><td>MtrunA17Chr3g0115131</td><td>Mt0006_00731</td><td>1.3</td><td>0.5</td><td>1.3</td><td>0</td><td>1.4</td></tr><tr><td>XDH1</td><td>Mtr.23395.1.S1_at</td><td>Medtr2g098030</td><td>MtrunA17Chr2g0328851</td><td>Mt0016_10367</td><td>80</td><td>52.3</td><td>37.3</td><td>39.2</td><td>52.7</td></tr><tr><td>AAO3</td><td>Mtr.29357.1.S1_at</td><td>Medtr5g087390</td><td>NC</td><td>Mt0010_10456</td><td>29.05</td><td>22.75</td><td>10.23</td><td>12.61</td><td>25.35</td></tr><tr><td>IAO3</td><td>Mtr.42638.1.S1_at</td><td>Medtr5g087410</td><td>NC</td><td>Mt0010_10457</td><td>4.95</td><td>9.36</td><td>4.05</td><td>19.45</td><td>62.19</td></tr><tr><td>NOFNiR1</td><td>Mtr.40060.1.S1_at</td><td>Medtr2g035460</td><td>MtrunA17Chr2g0296891</td><td>Mt0101_10060</td><td>43</td><td>29.2</td><td>4.1</td><td>18.8</td><td>125.4</td></tr><tr><td>NOFNiR2</td><td>Mtr.10348.1.S1_at</td><td>Medtr2g035470</td><td>MtrunA17Chr2g0296881</td><td>Mt0101_10060</td><td>43</td><td>29.2</td><td>4.1</td><td>18.8</td><td>125.4</td></tr><tr><td>NOFNiR3</td><td>Mtr.33463.1.S1_s_at</td><td>Medtr2g035480</td><td>MtrunA17Chr2g0296871</td><td>Mt0101_10060</td><td>43</td><td>29.2</td><td>4.1</td><td>18.8</td><td>125.4</td></tr><tr><td>SOX1</td><td>Mtr.28230.1.S1_at</td><td>Medtr6g023975</td><td>MtrunA17Chr6g0459641</td><td>Mt0033_10312</td><td>53</td><td>55.7</td><td>60.2</td><td>137.4</td><td>121.9</td></tr><tr><td>SOX2</td><td>NC</td><td>Medtr7g033410</td><td>MtrunA17Chr7g0227151</td><td>Mt0044_10132</td><td>3.3</td><td>1.5</td><td>2.6</td><td>4.6</td><td>9.2</td></tr></tbody></table>

Affymetrix gene codes were on (https://mtgea.noble.org/v3/), Mt4.0 gene codes were on (https://www.jcvi.org/medicago-truncatula-genome-database), Mt5.0 gene codes were on (https://medicago.toulouse.inra.fr/MtrunA17r5.0-ANR/), Symbimics data were in Roux et al. (2014). NC, no code available; DESEQ MEAN corresponds to the expression value of the RNA-sequencing analysis in the different zones of the nodule; FI, meristematic zone; FIID, distal infection zone; FIIP, proximal infection zone; IZ, interzone; ZIII, N₂-fixing zone. Shades of red highlight the highest expression value for each gene. NR, nitrate reductase; XDH, xanthine dehydrogenase; AAO, abscissic aldehyde oxidase; IAO, indole aldehyde oxidase; NOFNIR, NO forming nitrite reductase; SOX, sulfite oxidase.

ATP/ADP ratio) in 10 hpi and 4 dpi roots and in 4 wpi nodules after a 4 h treatment in the presence of NR effectors (Table 2). In control and NO₂⁻-treated roots and nodules, ATP/ADP ratio is

high, indicating that ATP-regenerating processes are not limited. In the presence of Tg, the ATP/ADP ratio is significantly decreased only in 4 wpi nodules, indicating that the inhibition of NR partially affects the energy state in mature nodules, but not in 10 hpi and 4 dpi roots. This decrease in the ATP/ADP ratio is not observed when 4 wpi nodules are incubated in the presence of both Tg and NO₂⁻, which means that the supply of NO₂⁻ makes it possible to maintain the nodule energy state.

Considered together, these results indicate that 1) both the NR and the mitochondrial ETC are involved in the production of NO, probably via the reduction of nitrate to nitrite by NR and the subsequent reduction of nitrite to NO at the mitochondrial ETC level, and 2) NR activity is linked to energy regeneration processes in mature nodules, but not during the first steps of the symbiosis.

## Involvement of NR in Carbon and Nitrogen Metabolism in Mature Nodules

Based on the above results, it appeared important to investigate more precisely the role of NR activity in the energetic metabolism of mature nodules. To this end, using a perfusion system adapted to NMR spectrometer (Roby et al., 1987), we first followed the effects of *in vivo* transition from normoxia to hypoxia, and *vice versa*, on the energy metabolism of 4 wpi nodules. **Figure 8A** displays typical ³¹P spectra obtained with living *M. truncatula* nodules during normoxia-hypoxia-normoxia transition experiments. The main resonance peaks are attributed to glucose-6-phosphate (G6P), cytosolic and vacuolar inorganic phosphate (Pi-cyt and Pi-vac, respectively), uridine di-phospho-glucose (UDPG), and to α, β and γ ATP. Transition from 21% O₂ (normoxia) to 1% O₂ (hypoxia) leads to an important reduction in ATP content and a significant acidification of cytoplasmic pH from 7.45 to 6.95 after 7 h of hypoxia, as measured by the shift of G6P and Pi-cyt resonances (Figure 8A). Back transition from 1 to 21% O₂ is accompanied by an increase in ATP and a progressive return of cytoplasmic pH to more alkaline values (7.2). In a second series of experiments, nodules were incubated at 21% O₂ in the

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Expression of *Medicago truncatula* NR genes during the symbiotic process. Expression analysis of the NR genes in roots and nodules during the symbiotic process. Short term kinetic 14 dpi (A, C, E), long term kinetic 8 wpi (B, D, F). Expression of NR1 (A, B), NR2 (C, D) and NR3 (E, F). Data are means ± SE of three biological replicates. Each measure was realized in three technical replicates. dpi, day post inoculation; wpi, week post inoculation.</figcaption></figure>
<figure><img src="image_6.png" /></figure>

**FIGURE 4** | Expression of *Medicago truncatula* NR genes at various times of the symbiotic process. The data are drawn from the same experiments as those presented in **Figure 3**. Values followed by different letters are significantly different according to one-way ANOVA analysis followed by a Fisher test (P<0.05). dpi, day post inoculation; hpi, hour post inoculation; wpi, week post inoculation.

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**FIGURE 5** | NR activity and NO production during the *Medicago truncatula* symbiosis. Short term kinetic 14 dpi (A), long term kinetic 8 wpi (B). The nitrate reductase activity is expressed in nmol per min per g of fresh weight. The fluorescence intensity of the NO production was measured using the DAF-2 fluorescent probe. Data are means ± SE of three biological replicates. Each measure was realized in three technical replicates. dpi, day post inoculation; FU, fluorescence unit; wpi, week post inoculation.

presence of 1 mM Tg (Figure 8B). After 5 h of Tg treatment, ATP level decreases by a factor 2.5, and cytoplasmic pH decreases from 7.4 to 6.95, indicating a decrease in energy state and an acidification of the nodules. The addition of 1 mM NO₂⁻ triggers an increase in cytosolic pH to 7.15, indicating a progressive recovery of the cell metabolism (Figure 8B). Thus, the inhibition of NR by Tg and its subsequent relief by NO₂⁻ partially mimics “normoxia-hypoxia-normoxia” transitions and support the involvement of NR in the energy metabolism of nodules.

Then, we investigated the potential involvement of NR in the N₂-fixing metabolism. To this end, mature nodules were incubated for 4 h in either the presence, or the absence, of 1 mM Tg, 10 mM NO₃⁻ and 1 mM NO₂⁻. As reported in Figure 9,

the nitrogenase activity, as measured by its acetylene reducing activity (ARA), is inhibited in the presence of either NO₃⁻ or Tg (or both) but is unaffected by the presence of NO₂⁻, indicating that NR activity supports N₂ fixation. The changes in the content of various metabolites related to carbon and nitrogen metabolism of the nodules are consistent with this hypothesis. Thus, following NR inhibition by Tg (either in the presence or absence of NO₃⁻), sucrose content increased, indicating a slowdown of its consumption by the nodules (Figure 9B). As a consequence, succinate and malate contents decreased (Figures 9C, D), indicating a lower supply of carbon nutrients to bacteroids, and asparagine content decreased (Figure 9E) as a result of the shortage of carbon substrate supply and ARA inhibition. Finally, alanine, which is a marker of hypoxia in
<figure><img src="image_9.png" /></figure>

**FIGURE 6** | NR activity and NO production in NR RNAi transformed *Medicago truncatula* nodules. (A), NR activity is expressed in nmol per min per g of fresh weight. (B). The fluorescence intensity of the NO production was measured using the DAF-2 fluorescent probe. Data are means ± SE of three biological replicates. Each measure was realized in three technical replicates. Asteriks * indicate statistical difference when compared with the control GFP (transgenic roots transformed with empty pK7GWIWG5D(II) vector selected on GFP marker) at P < 0.05 according to Student's t test. FU, fluorescence unit.

plant tissues (Gibbs and Greenway, 2003), was found to increase in the presence of Tg (Figure 9F), indicating that NR inhibition mimics a hypoxia situation. The presence of NO₂⁻ reverses the effects of Tg (Figure 9).

# DISCUSSION

## Medicago truncatula NRs Are Involved in NO Production

The main objective of this study was to test the hypothesis that NR regulates NO production and controls N₂-fixing metabolism during the *M. truncatula*-S. meliloti symbiosis.

The Medicago genome encodes three NR genes, MtNR1, MtNR2, and MtNR3 (Puppo et al., 2013; Roux et al., 2014). Whereas most of the species analyzed possess between one and two sequences of NR, *M. truncatula* have three NRs. MtN3 and MtNR1 are both located on chromosome 3 separate from only 20 kb and are grouped in the same HTU on the phylogenetic tree, elements which are in favor of an event of duplication of MtNR1 to give birth to MtNR3. This exhaustive search for NR sequences conducted to identify two potential orthologs to MtNR3 in

species closed to *M. truncatula* such as *Trifolium pratense* and *Cicer arietanum*. However, we did not identify three NR sequences in *T. pratense* and *C. arietanum*. MtNR3 has been shown to be specifically expressed in symbiotic nodules, excluding other organs of the plant (http://mtgea.noble.org/v3/). Its expression from 4 dpi (Figures 3 and 4) and location in the nodule primordium as well as in zones I and II of the nodule (Figure 2, Table 1) suggest that MtNR3 plays a role in cell division and infection processes during nodule organogenesis. Similarly, the sharp increase in its expression beyond 6 wpi (Figure 3F) suggests that NR3 also plays an important role in nodule senescence process that will be worth investigating. However, its low expression level during symbiosis (Figure 4) and its expression pattern, very different from the NO production once (Figures 3 and 5), do not support, a priori, a potential involvement of NR3 in the production of NO, pointing to a specific role of this NR in nodule aging.

MtNR1 and MtNR2 correspond to inducible and/or constitutive forms that have been found in other higher plants such as *A. thaliana* and *G. max* (Santucci et al., 1995). Our results show that three concomitant peaks of MtNR1/NR2 expression and NR activity can be considered over the period analyzed (Figures 3–5): i) during the first hours of the symbiotic interaction at10 hpi; ii) during the early development of the nodule at 4 dpi; and iii) when the nodule becomes mature around 3–4 wpi. The most salient feature of this study is the parallel that can be drawn between the expression of MtNR1 and MtNR2 (Figure 3), the total NR activity (Figure 5), and the production of NO (Figure 5). During the infection phase, at 10 hpi, the increase in MtNR1 and MtNR2 expression is accompanied by a slight increase in NR activity, while the production of NO increases by a factor of 2, suggesting that NRs are probably not the only sources of NO. In contrast, later in the symbiotic process (at 4 dpi and 3–4 wpi), the correlation between gene expression, NR activity, and NO production suggests that NO production is directly related to NR. Moreover, the spatiotemporal expression of MtNR1 and MtNR2 (Figure 2) corresponds to the localization of NO production reported to occur in the nodule primordium (del Giudice et al., 2011) and in the N₂-fixing zone of the mature nodule (Baudouin et al., 2006). While keeping in mind that MtNR2 is much more expressed than MtNR1, these observations strongly suggest that both NR1 and NR2 are involved in the production of NO during the symbiotic process. In a recent study dealing with the role of *M. truncatula* phytoglobins (Phytogbs) during N₂-fixing symbiosis we found that *Phytogb1.1* exhibits an expression pattern similar to that of NO production and is involved in its regulation during the different stages of the symbiosis (Berger et al., 2020). Thus, regardless of the potential mechanisms of post-translational regulation of Phytogbs and NRs activities, our results indicate that NRs probably work in coordination with Phytogb1.1 to regulate the level of NO at each stage of the symbiotic process.

The use of inhibitors of different NO sources (Figure 7) supports the hypothesis that at 10 hpi, 4 dpi, and 4 wpi, NRs are either directly, or indirectly, in combination with the mitochondrial ETC, involved in the production of NO. In fact,
<figure><img src="image_10.png" /><figcaption>**FIGURE 7** | Effects of inhibitors on NO production at various times of the *Medicago truncatula* symbiosis. (A), schematic representation of inhibitor targets. NO production of 10 hpi-old (B) and 4 dpi-old (C) roots and of 4 wpi-old nodules (D) was measured after 4 h in the presence of either 1 mM tungstate (Tg), 1 mM allopurinol (AP), 1 mM propyl gallate (PG), 1 mM nitrite (NO₂⁻), or 300 µM KCN. NO production is expressed as the percentage of the control values. Data are means ± SE of three biological replicates. Each measure was realized in three technical replicates. Values followed by different letters are significantly different according to one-way ANOVA analysis followed by a Fisher test (P < 0.05). dpi, day post inoculation; hpi, hour post inoculation; wpi, week post inoculation. AOX, alternative oxidase; COX, cytochrome oxidase; ETC, electron transfer chain; NR, nitrate reductase; XDH, xanthine dehydrogenase.</figcaption></figure>

<table><caption>**TABLE 2** | Effects of NR effectors on ATP/ADP ratio in *Medicago truncatula* roots and nodules.</caption><thead><tr><th rowspan="2"></th><th colspan="3">ATP/ADP ratio</th></tr><tr><th>10 hpi</th><th>4 dpi</th><th>4 wpi</th></tr></thead><tbody><tr><td>Control</td><td>6.1 ± 1.1</td><td>6.3 ± 1.7</td><td>6.8 ± 1.2</td></tr><tr><td>+1 mM NO<sub>2</sub><sup>-</sup></td><td>5.7 ± 0.9</td><td>4.7 ± 0.4</td><td>6.4 ± 0.7</td></tr><tr><td>+1 mM Tg</td><td>5.2 ± 0.5</td><td>5.6 ± 3.3</td><td>4.4 ± 0.4*</td></tr><tr><td>+1 mM NO<sub>2</sub><sup>-</sup> +1 mM Tg</td><td>4.9 ± 0.6</td><td>4.5 ± 1.0</td><td>5.5 ± 0.7</td></tr></tbody></table>

Either 10 hpi roots, 4 dpi roots, or 4 wpi-old nodules were incubated for 4 h in the presence of either 1 mM nitrite (NO₂⁻), 1 mM tungstate (Tg), or both. Data are the means ± SE of three (10 hpi, 4 dpi) and four (4 wpi) biological replicates. Each measure was realized in three technical replicates. Asterisk * indicates statistical difference with the control at the same time point at P < 0.05 according to Student's t test. dpi, day post inoculation; hpi, hour post inoculation; wpi, week post inoculation.

at these three stages the production of NO is inhibited by Tg. The use of Tg to inhibit NR activity has been successfully used in mature *M. truncatula* nodules to demonstrate the involvement of NR in NO production (Horchani et al., 2011). However, Tg is known to inhibit other molybdoenzymes than NR, including XDH, aldehyde oxidase, sulfite oxidase, and NO forming nitrite reductase (NOFNiR) (Mendel and Hänsch, 2002; Xiong et al., 2012, Chamizo-Ampudia et al., 2016), and the possibility that the Tg-dependent inhibition of NO production is linked to the inhibition of one of these other enzymes cannot be excluded.

Notwithstanding the stability and the specific activity of these enzymes in the root and nodule cells, the weak expression level of their genes remains much lower than those of MtNR1 and MtNR2 (Table 1, Figure S2), and one can reasonably think that NRs remain the major molybdoenzymes in cells. The use of the inhibitor allopurinol (Figure 7) shows that XDH, of which expression remains significant during the symbiotic process (Table 1, Figure S2), is perhaps partially involved in the production of NO during the first hours of symbiosis, but neither during nodule organogenesis nor in the N₂-fixing nodule. Furthermore, the preferential expression of XDH in the infection zone of the nodules (Table 1) reinforces the hypothesis of its involvement in the defense responses which occur at the beginning of the symbiotic process. XDH is a peroxisomal enzyme capable of producing NO⁺ and superoxide anion (O₂⁻) that can complex together to give peroxynitrite, ONOO⁻ (del Rio et al., 2004). ONOO⁻ is known to play a key role in the induction of defense responses, particularly via tyrosine nitration of proteins (Saito et al., 2006), which is consistent with the involvement of NO production in defense responses at 10 hpi.

Otherwise, the partial or total reversion of Tg-dependent inhibition of NO production by nitrite (Figure 7) means that NRs at least produce the nitrite necessary for the production of NO. Indeed, the inhibition of NO production by mitochondrial ETC inhibitors (Figure 7) supports the hypothesis that NR
<figure><img src="image_11.png" /></figure>

**FIGURE 8** | *In vivo* ³¹P-NMR study of metabolic transitions induced in *Medicago truncatula* nodules. (A) Proton-decoupled ³¹P-NMR spectra of nodules in normoxia at 21°C, perfused with a nutritive medium at pH 6.0 with successively 21% O₂ (control), 1% O₂, and then 21% O₂. For each line, duration of the treatment is indicated up to the end of the spectrum. (B) Proton-decoupled ³¹P-NMR spectra of nodules in normoxia at 21°C, perfused with a nutritive medium at pH 6.0 containing successively 21% O₂ (control), 1 mM tungstate (Tg) for 5 h, and then 1 mM nitrite (NO₂⁻) for 2 h. For each line, duration of the treatment is indicated up to the end of the spectrum. Each spectrum series is drawn from a representative experiment of eight (A) and five (B) biological replicates. The values indicated next to the G6P and Pi-cyt resonances correspond to the cytosolic pH values. Exponential apodization and zero filling. Pi, inorganic phosphate; G6P, Glucose-6-phosphate; Pi-cyt, cytoplasmic Pi; Pi-vac, vacuolar Pi; UDPG, UDP-glucose.

activity is indirectly involved in NO production *via* the reduction of nitrate to nitrite and the subsequent nitrite reduction to NO by mitochondrial ETC. This indicates that from the beginning of the symbiotic process, a major part of NO is produced *via* a nitrite reducing pathway.

We should also keep in mind that NO in the nodule is not produced only by the plant partner. Indeed, the bacterial partner was shown to produce from 33 to 90% of NO in *M. truncatula* (Horchani et al., 2011) and soybean (Sanchez et al., 2010) nodules, respectively. If the bacterial denitrification

pathway, including the periplasmic nitrate reductase (Nap) and nitrite reductase (Nir), has been described as the main enzymatic source of NO, additional genes encoding putative nitrate and nitrite reductase (called narB and NirN, respectively) have been recently identified that could also participate indirectly in NO synthesis (Ruiz et al., 2019). How the regulatory systems of the plant and the bacterial partners are coordinated to produce NO is one of the main issues to decipher the signaling and metabolic functions of NO at each stage of the symbiotic interaction.
<figure><img src="image_12.png" /><figcaption>**FIGURE 9** | Effects of NR effectors on nitrogenase activity and metabolite contents in *Medicago truncatula* nodules. 4 wpi-old nodules were incubated for 4 h in the presence of either 10 mM nitrate (NO₃⁻), 1 mM nitrite (NO₂⁻), or 1 mM tungstate (Tg). (A), nitrogenase activity (estimated as ARA); (B), sucrose; (C), succinate; (D), malate; (E), asparagine; (F), alanine. Data are means ± SE of three biological replicates. Each measure was realized in three technical replicates. Asterisks * indicate significant difference at P<0.05, when compared with the control (Ctrl) according to Student's *t* test.</figcaption></figure>

## NRs Regulate Energy State and Metabolism in Nitrogen-Fixing Nodules

The combined functioning of NR and mitochondrial ETC to produce NO has been associated with “Phytogb-NO” respiration and energy regeneration in hypoxic organs (Igamberdiev and Hill, 2009), and the role of NR in maintaining the ATP/ADP ratio has already been demonstrated in microoxic N₂-fixing nodules of

*M. truncatula* (Horchani et al., 2011). Our results (Table 2) show that, contrary to what is observed in 4 wpi-old nodules, the inhibition of NRs at the beginning of the symbiotic process (10 hpi and 4 dpi) does not cause a decrease in the ATP/ADP ratio indicating that, at these stages, energy regeneration is not linked to the functioning of Phytogb-NO respiration. These observations were expected since, contrary to mature nodules, roots are normoxic organs, and the cellular energy is supposed to be regenerated by the O₂-dependent mitochondrial respiration. The role of mitochondrial ETC in NO production during the early days of the symbiotic process remains to be clarified.

The role of NR in energy, carbon, and nitrogen metabolism was investigated by using the NR inhibitor Tg. *In vivo* ³¹P-NMR experiments (Figure 8A) showed that hypoxia leads to a fall in ATP and to the acidification of cytoplasmic pH inside the nodules. These data are consistent with earlier observations with maize root tips (Roberts et al., 1984), soybean nodules (Pfeffer et al., 1992) and sycamore maple cells (Gout et al., 2001). In the presence of Tg, the fall in ATP and the cytoplasm acidification resulting from the inhibition of NR and their partial reversion by the addition of nitrite (Figure 8B) confirm that NR activity is fully involved in the energy metabolism of nodules via “Phytogb-NO” respiration, as already observed by Horchani et al. (2011). Logically, the inhibition of energy metabolism triggers the inhibition of the carbon and nitrogen metabolism in nodules, with an accumulation of sucrose, the subsequent decrease in the supply of carbon substrates to bacteroids (succinate, malate), and ultimately, the decrease in nitrogen reduction (ARA) and assimilation (asparagine) activities (Figure 9). It is interesting to note that treatment of nodules with Tg leads to an increase in alanine content (Figure 9F). In hypoxic/anoxic conditions, alanine accumulates in tissues following the inhibition of O₂-dependent respiration and the induction of alanine aminotransferase (Limami et al., 2008, and references inside), and to date alanine is considered as one of the most prevalent and ubiquitous markers of hypoxia in plants (Gibbs and Greenway, 2003; Bailey-Serres and Voesenek, 2008; Limami et al., 2014). The accumulation of alanine, as well as the drop in the energy state (Table 2), is the demonstration that the inhibition of NR activity leads to the inhibition of the mitochondrial respiratory chain, the tricarboxylic acid cycle, and ultimately the carbon and nitrogen metabolism in the same way as the drop in O₂ content in hypoxic tissue.

This conclusion is reinforced by the results obtained in the presence of nitrite. Interestingly, our study shows that over a short period of time, i.e. 4 h, nitrite stimulates ARA, while nitrate inhibits it (Figure 9). The inhibitory effect of nitrate on nitrogen fixation is widely recognized, and the question arose as to whether this effect was potentially linked to nitrite. Indeed, nitrite is known to be a potent inhibitor of nitrogenase *in vitro* (Wong, 1980; Trinchant and Rigaud, 1982). However, several studies have shown that nitrite is not responsible for the *in vivo* inhibition of nitrogenase activity induced by nitrate (Streeter, 1985a; Streeter, 1985b; Becana et al., 1989; Arrese-Igor et al., 1997; Arrese-Igor et al., 1998). The present study confirms this point and shows that over a short period of time, nitrite not only
increases ARA, but even stimulates the energy, carbon, and nitrogen metabolism of nodules (Figures 7–9, Table 2). In a previous *in vivo* ³¹P-NMR study carried out on anoxic maize roots, Libourel et al. (2006) reported that the addition of nitrite makes it possible to limit the drop in cytoplasmic pH following anoxic treatment. Here, we show that 1) similarly to hypoxia, NR inhibition induces a drop in cytoplasmic pH, and 2) this drop in cytoplasmic pH may be reversed by the addition of nitrite (Figure 8), confirming the key role of NR and PNR, in the functioning of microoxic nodules.

Considered together, our data first confirm that NR is a source of nitrite, and indirectly of NO, during the establishment and functioning of N₂-fixing symbiosis. Second, they support the hypothesis that NRs, in close collaboration with Phytogb1.1 (Berger et al., 2020), are strongly involved in the regulation of cellular energy and N₂-fixing metabolism through the functioning of the Phytogb-NO respiration. Finally, they show that NR activity is needed for the N₂-fixing symbiosis and constitutes a first attempt to explain its dual role in nodule functioning: 1) generating NO as a signal for gene regulation and metabolic adaptation, and 2) contributing to the energy supply under the hypoxic conditions prevailing inside the nodule.

## DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the article/ Supplementary Material.

## REFERENCES

Arrese-Igor, C., Garcia-Plazaola, J.II, Hernández, A., and Aparicio-Tejo, P. M. (1990). Effect of low nitrate supply to nodulated lucerne on time course of activites of enzymes involved in inorganic nitrogen metabolism. *Physiol. Plant* 80, 185-190. doi: 10.1111/j.1399-3054.1990.tb04394.x

Arrese-Igor, C., Minchin, F. R., Gordon, A. J., and Nath, A. K. (1997). Possible causes of the physiological decline in soybean nitrogen fixation in the presence of nitrate. *J. Exp. Bot.* 48, 905-913. doi: 10.1093/jxb/48.4.905

Arrese-Igor, C., Gordon, A. J., Minchin, F. R., and Denison, R. F. (1998). Nitrate entry and nitrite formation in the infected region of soybean nodules. *J. Exp. Bot.* 49, 41-48. doi: 10.1093/jxb/49.318.41

Bailey-Serres, J., and Voesenek, L. A. C. J. (2008). Flooding stress: acclimations and genetic diversity. *Annu. Rev. Plant Biol.* 59, 313-339. doi: 10.1146/annurev.arplant.59.032607.092752

Bari, R., Pant, B. D., Stitt, M., and Scheible, W.-R. (2006). PHO2, MicroRNA399, and PHR1 Define a Phosphate-Signaling Pathway in Plants. *Plant Physiol.* 141, 988-999. doi: 10.1104/pp.106.079707

Baudouin, E., Pieuchot, L., Engler, G., Pauly, N., and Puppo, A. (2006). Nitric oxide is formed in Medicago truncatula-Sinorhizobium meliloti functional nodules. *Mol. Plant-Microbe Interact.* 19, 970-975. doi: 10.1094/MPMI-19-0970

Becana, M., Minchin, F. R., and Sprent, J. II (1989). Short-term inhibition of legume N2 fixation by nitrate. *Planta* 180, 40-45. doi: 10.1007/BF02411408

Becana, M., Yruela, I., Sarath, G., Catalán, P., and Hargrove, M. S. (2020). Plant hemoglobins: a journey from unicellular green algae to vascular plants. *New Phytol.* 227, 1618-1635. doi: 10.1111/nph.16444

Berger, A., Boscari, A., Frendo, P., and Brouquisse, R. (2019). Nitric oxide signaling, metabolism and toxicity in nitrogen-fixing symbiosis. *J. Exp. Bot.* 70, 4505-4520. doi: 10.1093/jxb/erz159

Berger, A., Guinand, S., Boscari, A., Puppo, A., and Brouquisse, R. (2020). Medicago truncatula Phytoglobin 1.1 controls symbiotic nodulation and

## AUTHOR CONTRIBUTIONS

ABe, ABo, DR, AP, and RB planned and designed the research. ABe, ABo, MM, MH, NH, SB, DR, and RB performed the experiments. ABe, ABo, DR, AP, and RB analyzed and interpreted the data. ABe, ABo, and RB wrote the manuscript.

## ACKNOWLEDGMENTS

This work was supported by the “Institut National de la Recherche en Agronomie”, the “Centre National de la Recherche Scientifique”, the Côte d’Azur University, and the French Government (National Research Agency, ANR) through the LABEX SIGNALIFE program (reference # ABR-11-LABX-0028-01) and the STAYPINK project (ANR-15-CE20-0005). NMR and HPLC analyses were performed at the Bordeaux Metabolome Facility that is partly supported by the national infrastructure for metabolomics and fluxomics MetaboHUB (ANR-11-INBS-0010) grant. Some of the results presented in this manuscript are part of the PhD thesis of ABe (Berger, 2019).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.01313/full#supplementary-material

nitrogen fixation via the regulation of nitric oxide level. *New Phytol.* 227, 84-98. doi: 10.1111/nph.16462

Berger, A. (2019). *Contribution des phytoglobines et des nitrate réductase à la régulation de l'oxyde nitrique et de la fixation de l'azote dans la symbiose Medicago truncatula / Sinorhizobium meliloti*. *PhD thesis* (Nice (France): Université Côte d'Azur).

Besson-Bard, A., Pugin, A., and Wendehenne, D. (2008). New insights into nitric oxide signaling in plants. *Annu. Rev. Plant Biol.* 59, 21-39. doi: 10.1146/annurev.arplant.59.032607.092830

Bobik, C., Meilhoc, E., and Batut, J. (2006). FixJ: A major regulator of the oxygen limitation response and late symbiotic functions of Sinorhizobium meliloti. *J. Bacteriol.* 188, 4890-4902. doi: 10.1128/JB.00251-06

Boisson-Dernier, A., Chabaud, M., Garcia, F., Bécard, G., Rosenberg, C., and Barker, D. G. (2001). Agrobacterium rhizogenes-transformed roots of Medicago truncatula for the study of nitrogen-fixing and endomycorrhizal symbiotic associations. *Mol. Plant Microbe Interact.* 14, 695-700. doi: 10.1094/MPMI.2001.14.6.695

Boscari, A., del Giudice, J., Ferrarini, A., Venturini, L., Zaffini, A.-L., Delledonne, M., et al. (2013). Expression Dynamics of the Medicago truncatula Transcriptome during the Symbiotic Interaction with Sinorhizobium meliloti: Which Role for Nitric Oxide? *Plant Physiol.* 161, 425-439. doi: 10.1104/pp.112.208538

Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. *Anal. Biochem.* 72, 248-254. doi: 10.1016/0003-2697(76)90527-3

Brouquisse, R., James, F., Raymond, P., and Pradet, A. (1991). Study of glucose starvation in excised maize root tips. *Plant Physiol.* 96, 619 LP - 626. doi: 10.1104/pp.96.2.619

Cam, Y., Pierre, O., Boncompagni, E., Hérouart, D., Meilhoc, E., and Bruand, C. (2012). Nitric oxide (NO): A key player in the senescence of Medicago truncatula root nodules. *New Phytol.* 196, 548-560. doi: 10.1111/j.1469-8137.2012.04282.x
Campbell, W. H. (1999). Nitrate reductase structure, function and regulation: Bridging the Gap between Biochemistry and Physiology. *Annu. Rev. Plant Physiol. Plant Mol. Biol.* 50, 277–303. doi: 10.1146/annurev.arplant.50.1.277

Castella, C., Mirtziou, I., Seassau, A., Boscari, A., Montrichard, F., Papadopoulou, K., et al. (2017). Post-translational modifications of Medicago truncatula glutathione peroxidase 1 induced by nitric oxide. *Nitric. Oxide* 68, 125–136. doi: 10.1016/j.niox.2017.02.004

Chamizo-Ampudia, A., Sanz-Luque, E., Llamas, Á, Ocaña-Calahorro, F., Mariscal, V., Carreras, A., et al. (2016). A dual system formed by the ARC and NR molybdoenzymes mediates nitrite-dependent NO production in Chlamydomonas. *Plant Cell Environ.* 39, 2097–2107. doi: 10.1111/pce.12739

Chamizo-Ampudia, A., Sanz-Luque, E., Llamas, A., Galvan, A., and Fernandez, E. (2017). Nitrate Reductase Regulates Plant Nitric Oxide Homeostasis. *Trends Plant Sci.* 22, 163–174. doi: 10.1016/j.tplants.2016.12.001

Dean, J. V., and Harper, J. E. (1988). The Conversion of Nitrite to Nitrogen Oxide(s) by the Constitutive NAD(P)H-Nitrate Reductase Enzyme from Soybean. *Plant Physiol.* 88, 389 LP – 395. doi: 10.1104/pp.88.2.389

del Giudice, J., Cam, Y., Damiani, I., Fung-Chat, F., Meilhoc, E., Bruand, C., et al. (2011). Nitric oxide is required for an optimal establishment of the medicago truncatula-sinorhizobium meliloti symbiosis. *New Phytol.* 191, 405–417. doi: 10.1111/j.1469-8137.2011.03693.x

del Rio, L. A., Javier Corpas, F., and Barroso, J. B. (2004). Nitric oxide and nitric oxide synthase activity in plants. *Phytochemistry* 65, 783–792. doi: 10.1016/j.phytochem.2004.02.001

Ferrarini, A., De Stefano, M., Baudouin, E., Pucciariello, C., Polverari, A., Puppo, A., et al. (2008). Expression of Medicago truncatula genes responsive to nitric oxide in pathogenic and symbiotic conditions. *Mol. Plant Microbe Interact.* 21, 781–790. doi: 10.1094/MPMI-21-6-0781

Fukudome, M., Calvo-Begueria, L., Kado, T., Osuki, K.-I., Rubio, M. C., Murakami, E.-I., et al. (2016). Hemoglobin LjGlb1-1 is involved in nodulation and regulates the level of nitric oxide in the Lotus japonicus-Mesorhizobium loti symbiosis. *J. Exp. Bot.* 67, 5275–5283. doi: 10.1093/jxb/erw290

Gibbs, J., and Greenway, H. (2003). Review: Mechanisms of anoxia tolerance in plants. I. Growth, survival and anaerobic catabolism. *Funct. Plant Biol.* 30, 1–47. doi: 10.1071/PP98095

Gout, E., Boisson, A.-M., Aubert, S., Douce, R., and Bligny, R. (2001). Origin of the Cytoplasmic pH Changes during Anaerobic Stress in Higher Plant Cells. Carbon-13 and Phosphorous-31 Nuclear Magnetic Resonance Studies. *Plant Physiol.* 125, 912 LP – 925. doi: 10.1104/pp.125.2.912

Gupta, K. J., Stoimenova, M., and Kaiser, W. M. (2005). In higher plants, only root mitochondria, but not leaf mitochondria reduce nitrite to NO, in vitro and in situ. *J. Exp. Bot.* 56, 2601–2609. doi: 10.1093/jxb/eri252

Hancock, J. T. (2019). Considerations of the importance of redox state for reactive nitrogen species action. *J. Exp. Bot.* 70, 4323–4331. doi: 10.1093/jxb/erz067

Hardy, R. W., Holsten, R. D., Jackson, E. K., and Burns, R. C. (1968). The acetylene-ethylene assay for N₂ fixation: laboratory and field evaluation. *Plant Physiol.* 43, 1185–1207. doi: 10.1104/pp.43.8.1185

Hichri, I., Boscari, A., Castella, C., Rovere, M., Puppo, A., and Brouquisse, R. (2015). Nitric oxide: A multifaceted regulator of the nitrogen-fixing symbiosis. *J. Exp. Bot.* 66, 2877–2887. doi: 10.1093/jxb/erv051

Hichri, I., Boscari, A., Meilhoc, E., Català, M., Barreno, E., Bruand, C., et al. (2016). “Nitric oxide: a multitask player in plant-microorganism symbioses,” in *Gasotransmitters in Plants: The Rise of a New Paradigm in Cell Signaling*. Eds. L. Lamattina and C. Garcia-Mata (Switzerland: Springer International Publishing), 239–268. doi: 10.1007/978-3-319-40713-5_12

Hilliou, F., and Tran, T. (2013). RqPCRAnalysis: Analysis of quantitative real-time PCR data. *Bioinforma. 2013* - Proc. *Int. Conf. Bioinforma. Model. Methods Algorithms*. (Science and Technology Publications, Lda), 202–211. doi: 10.5220/0004312002020211

Horchani, F., Prévot, M., Boscari, A., Evangelisti, E., Meilhoc, E., Bruand, C., et al. (2011). Both plant and bacterial nitrate reductases contribute to nitric oxide production in Medicago truncatula nitrogen-fixing nodules. *Plant Physiol.* 155, 1023–1036. doi: 10.1104/pp.110.166140

Igamberdiev, A. U., and Hill, R. D. (2004). Nitrate, NO and haemoglobin in plant adaptation to hypoxia: An alternative to classic fermentation pathways. *J. Exp. Bot.* 55, 2473–2482. doi: 10.1093/ixb/erh272

Igamberdiev, A. U., and Hill, R. D. (2009). Plant mitochondrial function during anaerobiosis. *Ann. Bot.* 103, 259–268. doi: 10.1093/aob/mcn100

Jones, D. T., Taylor, W. R., and Thornton, J. M. (1992). The rapid generation of mutation data matrices. *Comput. Appl. Biosci.* 8, 275282. doi: 10.1093/bioinformatics/8.3.275

Karimi, M., Inze, D., and Depicker, A. (2002). GATEWAY vectors for Agrobacterium-mediated plant transformation. *Trends Plant Sci.* 7, 193–195. doi: 10.1016/S1360-1385(02)02251-3

Kato, K., Okamura, Y., Kanahama, K., and Kanayama, Y. (2003). Nitrate-independent expression of plant nitrate reductase in Lotus japonicus root nodules. *J. Exp. Bot.* 54, 1685–1690. doi: 10.1093/jxb/erg189

Kato, K., Kanahama, K., and Kanayama, Y. (2010). Involvement of nitric oxide in the inhibition of nitrogenase activity by nitrate in Lotus root nodules. *J. Plant Physiol.* 167, 238–241. doi: 10.1016/j.jplph.2009.08.006

Kolbert, Z., Barroso, J. B., Brouquisse, R., et al. (2019). A forty year journey: The generation and roles of NO in plants. *Nitric. Oxide* 93, 53–70. doi: 10.1016/j.niox.2019.09.006

Kolbert, Z., Ortega, L., and Erdei, L. (2010). Involvement of nitrate reductase (NR) in osmotic stress-induced NO generation of Arabidopsis thaliana L. roots. *J. Plant Physiol.* 167, 77–80. doi: 10.1016/j.jplph.2009.08.013

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: Molecular evolutionary genetics analysis across computing platforms. *Mol. Biol. Evol.* 35 (6), 15471549. doi: 10.1093/molbev/msy096

Libourel, I. G. L., Bethke, P. C., De Michele, R., and Jones, R. L. (2006). Nitric oxide gas stimulates germination of dormant Arabidopsis seeds: use of a flow-through apparatus for delivery of nitric oxide. *Planta* 223, 813–820. doi: 10.1007/s00425-005-0117-8

Limami, A. M., Glévarec, G., Ricoult, C., Cliquet, J.-B., and Planchet, E. (2008). Concerted modulation of alanine and glutamate metabolism in young Medicago truncatula seedlings under hypoxic stress. *J. Exp. Bot.* 59, 2325–2335. doi: 10.1093/jxb/ern102

Limami, A. M., Diab, H., and Lothier, J. (2014). Nitrogen metabolism in plants under low oxygen stress. *Planta* 239, 531–541. doi: 10.1007/s00425-013-2015-9

Lin, Y.-F., Hu, Y.-H., Lin, H.-T., Liu, X., Chen, Y.-H., Zhang, S., et al. (2013). Inhibitory Effects of Propyl Gallate on Tyrosinase and Its Application in Controlling Pericarp Browning of Harvested Longan Fruits. *J. Agric. Food Chem.* 61, 2889–2895. doi: 10.1021/jf305481h

Madeira, F., Park, Y. M., Lee, J., Buso, N., Gur, T., Madhusoodanan, N., et al. (2019). The EMBL-EBI search and sequence analysis tools APIs in 2019. *Nucleic Acids Res.* 47 (W1), W636W641. doi: 10.1093/nar/gkz268

Mathieu, C., Sophie, M., Frendo, P., Puppo, A., and Davies, M. (1998). Direct detections of radicals in intact soybean nodules: presence of nitric oxide-leghemoglobin complexes. *Free Radic. Biol. Med.* 24, 1242–1249. doi: 10.1016/S0891-5849(97)00440-1

Meilhoc, E., Cam, Y., Skapski, A., and Bruand, C. (2010). The response to nitric oxide of the nitrogen-fixing symbiont Sinorhizobium meliloti. *Mol. Plant Microbe Interact.* 23, 748–759. doi: 10.1094/MPMI-23-6-0748

Melo, P. M., Silva, L. S., Ribeiro, I., Seabra, A. R., and Carvalho, H. G. (2011). Glutamine Synthetase Is a Molecular Target of Nitric Oxide in Root Nodules of Medicago truncatula and Is Regulated by Tyrosine Nitration. *Plant Physiol.* 157, 1505–1517. doi: 10.1104/pp.111.186056

Mendel, R. R., and Hänsch, R. (2002). Molybdoenzymes and molybdenum cofactor in plants. *J. Exp. Bot.* 53, 1689–1698. doi: 10.1093/jxb/erf038

Mergaert, P., Nikovics, K., Kelemen, Z., Maunoury, N., Vaubert, D., Kondorosi, A., et al. (2003). A Novel Family in Medicago truncatula Consisting of More Than 300 Nodule-Specific Genes Coding for Small, Secreted Polypeptides with Conserved Cysteine Motifs. *Plant Physiol.* 132, 161–173. doi: 10.1104/pp.102.018192

Miranda, K. M., Espey, M. G., and Wink, D. A. (2001). A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. *Nitric. Oxide Biol. Chem.* 5, 62–71. doi: 10.1006/niox.2000.0319

Moen, A., Svanella, L., Rolin, D., Gaudillère, M., Gaudillère, J., and Monet, R. (1998). Compositional changes during the fruit development of two peach cultivars differing in juice acidity. *J. Amer. Soc Hortic. Sci.* 123, 770–775. doi: 10.21273/JASHS.123.5.770

Nagata, M., Murakami, E., Shimoda, Y., Shimoda-Sasakura, F., Kucho, K., Suzuki, A., et al. (2008). Expression of a class 1 hemoglobin gene and production of nitric oxide in response to symbiotic and pathogenic bacteria in Lotus japonicus. *Mol. Plant-Microbe Interact.* 21, 1175–1183. doi: 10.1094/MPMI-21-9-1175
Navascues, J., Perez-Rontome, C., Gay, M., Marcos, M., Yang, F., Walker, F., et al. (2012). Leghemoglobin green derivatives with nitrated hemes evidence production of highly reactive nitrogen species during aging of legume nodules. *Proc. Natl. Acad. Sci.* 109, 2660–2665. doi: 10.1073/pnas.1116559109

Pfeffer, P., Rolin, D., Kumosinski, T., MacFall, J., and Schmidt, J. (1992). ³¹P relaxation responses associated with N2/O2 diffusion in soybean nodule cortical cells and excised cortical tissue. *Plant Physiol.* 100, 1682–1690. doi: 10.1104/pp.100.4.1682

Pii, Y., Crimi, M., Cremonese, G., Spena, A., and Pandolfini, T. (2007). Auxin and nitric oxide control indeterminate nodule formation. *BMC Plant Biol.* 7, 21. doi: 10.1186/1471-2297-7-21

Planchet, E., Jagadis Gupta, K., Sonoda, M., and Kaiser, W. M. (2005). Nitric oxide emission from tobacco leaves and cell suspensions: rate limiting factors and evidence for the involvement of mitochondrial electron transport. *Plant J.* 41, 732–743. doi: 10.1111/j.1365-313X.2005.02335.x

Puppo, A., Pauly, N., Boscari, A., Mandon, K., and Brouquisse, R. (2013). Hydrogen Peroxide and Nitric Oxide: Key Regulators of the Legume—*Rhizobium* and Mycorrhizal Symbioses. *Antioxid. Redox Signal.* 18, 2202–2219. doi: 10.1089/ars.2012.5136

Roberts, J. K. M., Ray, P. M., Wade-Jardetzky, N., and Jardetzky, O. (1980). Estimation of cytoplasmic and vacuolar pH in higher plant cells by 31P NMR. *Nature* 283, 870–872. doi: 10.1038/283870a0

Roberts, J. K., Callis, J., Wemmer, D., Walbot, V., and Jardetzky, O. (1984). Mechanisms of cytoplasmic pH regulation in hypoxic maize root tips and its role in survival under hypoxia. *Proc. Natl. Acad. Sci.* 81, 3379 LP – 3383. doi: 10.1073/pnas.81.11.3379

Roby, C., Martin, J. B., Bligny, R., and Douce, R. (1987). Biochemical changes during sucrose deprivation in higher plant cells. Phosphorus-31 nuclear magnetic resonance studies. *J. Biol. Chem.* 262, 5000–5007.

Rockel, P., Strube, F., Rockel, A., Wildt, J., and Kaiser, W. M. (2002). Regulation of nitric oxide (NO) production by plant nitrate reductase in vivo and in vitro. *J. Exp. Bot.* 53, 103–110. doi: 10.1093/jexbot/53.366.103

Rolin, D. B., Boswell, R. T., Sloger, C., Tu, S.-I., and Pfeffer, P. E. (1989). In Vivo ³¹P NMR Spectroscopic Studies of Soybean Bradyrhizobium Symbiosis. *Plant Physiol.* 89, 1238–1246. doi: 10.1104/pp.89.4.1238

Roux, B., Rodde, N., Jardinaud, M. F., Timmers, T., Sauviac, L., Cottret, L., et al. (2014). An integrated analysis of plant and bacterial gene expression in symbiotic root nodules using laser-capture microdissection coupled to RNA sequencing. *Plant J.* 77, 817–837. doi: 10.1111/tpj.12442

Ruiz, B., Le Scornet, A., Sauviac, L., Rémy, A., Bruand, C., and Meilhoc, E. (2019). The Nitrate Assimilatory Pathway in Sinorhizobium meliloti: Contribution to NO Production. *Front. Microbiol.* 10:1526. doi: 10.3389/fmicb.2019.01526

Sainz, M., Calvo-Begueria, L., Pérez-Rontomé, C., Wienkoop, S., Abian, J., Staudinger, C., et al. (2015). Leghemoglobin is nitrated in functional legume nodules in a tyrosine residue within the heme cavity by a nitrite/peroxide-dependent mechanism. *Plant J.* 81, 723–735. doi: 10.1111/tpj.12762

Saito, S., Yamamoto-Katou, A., Yoshioka, H., Doke, N., and Kawakita, K. (2006). Peroxynitrite Generation and Tyrosine Nitration in Defense Responses in Tobacco BY-2 Cells. *Plant Cell Physiol.* 47, 689–697. doi: 10.1093/pcp/pcj038

Sakihama, Y., Nakamura, S., and Yamasaki, H. (2002). Nitric Oxide Production Mediated by Nitrate Reductase in the Green Alga Chlamydomonas reinhardtii: an Alternative NO Production Pathway in Photosynthetic Organisms. *Plant Cell Physiol.* 43, 290–297. doi: 10.1093/pcp/pcf034

Sanchez, C., Gates, A. J., Meakin, G. E., Uchiumi, T., Girard, L., Richardson, D. J., et al. (2010). Production of nitric oxide and nitrosylleghemoglobin complexes in soybean nodules in response to flooding. *Mol. Plant-Microbe Interact.* 23, 702–711. doi: 10.1094/MPMI-23-5-0702

Santucci, D. M., Haas, B., and Smarrelli, J. (1995). Regulation of the inducible soybean nitrate reductase isoform in mutants lacking constitutive isoform(s). *Biochim. Biophys. Acta - Protein Struct. Mol. Enzymol.* 1247, 46–50. doi: 10.1016/0167-4838(94)00200-Z

Seabra, A. R., and Carvalho, H. G. (2015). Glutamine synthetase in Medicago truncatula, unveiling new secrets of a very old enzyme. *Front. Plant Sci.* 6, 578. doi: 10.3389/fpls.2015.00578

Silveira, J. A. G., Matos, J. C. S., Cecatto, V. M., Viegas, R. A., and Oliveira, J. T. A. (2001). Nitrate reductase activity, distribution, and response to nitrate in two contrasting Phaseolus species inoculated with Rhizobium spp. *Environ. Exp. Bot.* 46, 37–46. doi: 10.1016/S0098-8472(01)00082-X

Streeter, J. G. (1985a). Nitrate Inhibition of Legume Nodule Growth and Activity: II. Short Term Studies with High Nitrate Supply. *Plant Physiol.* 77, 325–328. doi: 10.1104/pp.77.2.325

Streeter, J. G. (1985b). Nitrate inhibition of legume nodule growth and activity: I. Long term studies with a continuous supply of nitrate. *Plant Physiol.* 77, 321–324. doi: 10.1104/pp.77.2.321

Terpolilli, J. J., Hood, G. A., Poole, P. S.R. K. B. T.-A (2012). Chapter 5 - What determines the efficiency of N2-fixing rhizobium-legume symbioses? *Adv. Microb. Physiol.* 60, 325–389.

Trinchant, J. C., and Rigaud, J. (1982). Nitrite and nitric oxide as inhibitors of nitrogenase from soybean bacteroids. *Appl. Environ. Microbiol.* 44, 1385–1388. doi: 10.1128/AEM.44.6.1385-1388.1982

Truong, H. N., Thelineau, E., Bonneau, L., Fournier, C., Potin, S., Balzergue, S., et al. (2015). The Medicago truncatula hypermycorrhizal B9 mutant displays an altered response to phosphate and is more susceptible to *Aphanomyces euteiches*. *Plant Cell Environ.* 38, 73–88. doi: 10.1111/pce.12370

Udvardi, M. K., and Day, D. A. (1997). Metabolite transport across symbiotic membranes of legume nodules. *Annu. Rev. Plant Physiol. Plant Mol. Biol.* 48, 493–523. doi: 10.1146/annurev.arplant.48.1.493

Van de Velde, W., Guerra, J. C. P., De Keyser, A., De Rycke, R., Rombauts, S., Maunoury, N., et al. (2006). Aging in legume symbiosis. A molecular view on nodule senescence in Medicago truncatula. *Plant Physiol.* 141, 711–720. doi: 10.1104/pp.106.078691

Velterop, J. S., and Vos, F. (2001). A rapid and inexpensive microplate assay for the enzymatic determination of glucose, fructose, sucrose, L-malate and citrate in tomato (Lycopersicon esculentum) extracts and in orange juice. *Phytochem. Anal.* 12, 299–304. doi: 10.1002/pca.598

Wang, Q., Liu, J., and Zhu, H. (2018). Genetic and Molecular Mechanisms Underlying Symbiotic Specificity in Legume-Rhizobium Interactions. *Front. Plant Sci.* 9:313. doi: 10.3389/fpls.2018.000313

Wong, P. P. (1980). Nitrate and Carbohydrate Effects on Nodulation and Nitrogen Fixation (Acetylene Reduction) Activity of Lentil (Lens esculenta Moench). *Plant Physiol.* 66, 78–81. doi: 10.1104/pp.66.1.78

Xiao, T. T., Schilderink, S., Moling, S., Deinum, E. E., Kondorosi, E., Franssen, H., et al. (2014). Fate map of Medicago truncatula root nodules. *Development* 141, 3517–3528. doi: 10.1242/dev.110775

Xiong, J., Fu, G., Yang, Y., Zhu, C., and Tao, L. (2012). Tungstate: is it really a specific nitrate reductase inhibitor in plant nitric oxide research? *J. Exp. Bot.* 63, 33–41. doi: 10.1093/jxb/err268

Yamasaki, H., and Sakihama, Y. (2000). Simultaneous production of nitric oxide and peroxynitrite by plant nitrate reductase: in vitro evidence for the NR-dependent formation of active nitrogen species. *FEBS Lett.* 468, 89–92. doi: 10.1016/S0014-5793(00)01203-5

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Berger, Boscari, Horta Araújo, Maucourt, Hanchi, Bernillon, Rolin, Puppo and Brouquisse. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.