The crop and meadow plants selected for this study were ryegrass (Lolium perenne var. Herbus), clover (Trifolium repens L var. Huia), wheat (Triticum aestivum L. var. Cezanne), tomato (Solanum lycopersicum L cv. Agora Hybrid F1), Arabidopsis (A. thaliana Col-0), oilseed rape (Brassica napus cv. Neptune) and Brachypodium (B. distachyon Bd21). Table 1 compiles the total of 40 experimental growth conditions for both leaves and roots.

Plants were germinated and cultured under different environmental conditions (Table 1) in a phytotron (Servicio Fitotrón e Invernadero SGIKer, UPV/EHU) in a common 14/10 h light/darkness photoperiod with 22–24/18°C temperature regime and 60/70–80% relative humidity, respectively. Nitrogen source was ammonium or nitrate, applied as Ca(NO₃)₂ or (NH₄)₂SO₄, respectively, at concentration from 0.5 to 15 mM depending on the plant species and culture medium. Light intensity (from 200–500 μmol photon m⁻² s⁻¹) was selected according to the physiological characteristics of the plant species. For wheat, 700 μmol photon m⁻² s⁻¹ were also supplied as high light intensity. Detailed germination and growth conditions, as well as precise composition of the nutrient solutions, are given for each species in the literature referenced in Table 1.

Ammonium was determined from root and shoot extracts with the phenolhypochlorite-based colorimetric method or by ion-exchange chromatography in the case of ryegrass and clover (Dionex 600 equipment; see references in Table 1). For Arabidopsis, total free amino acids were quantified by the ninhydrin method (Yemm et al., 1955) as Sarasketa et al. (2014). For ryegrass and clover, amino acids were extracted as Belastegui-Macadam et al. (2007) and quantified by HPLC (JEOL Aminotac 500 JLC-500IV); and for wheat, tomato, Brachypodium, and oilseed by capillary electrophoresis (PA-800, Beckman Coulter Inc., United States; Setién et al., 2013).

Soluble protein was extracted from frozen tissues with similar extraction buffers as described in Setién et al. (2013) for wheat, Setién et al. (2014) for tomato, ryegrass, and clover, and in Sarasketa et al. (2016) for Arabidopsis, Brachypodium, and oilseed rape. Soluble protein content was determined by a dye binding protein assay (Bio-Rad Bradford Protein assay) using bovine serum albumin (BSA) as standard (Bradford, 1976).

Enzyme activities were determined in a 96-well microplate reader (BioTek Instruments) from protein extracts. GS activity was measured monitoring the formation of γ-glutamilhydroxamate (γ-GHM) in a semibiosynthetic assay (Setién et al., 2013). GDH, PEPC, ICDH, and NAD-ME activities were measured monitoring the evolution of NAD(P)H monitored at 340 nm. GDH activity was determined in aminating sense (Setién et al., 2013). For ryegrass, wheat, and clover, novel data of carbon enzymes activities (PEPC, ICDH, and NAD-ME) were determined as previously described by Vega-Mas et al. (2015) or Sarasketa et al. (2016).

A total of 80 observations, given for root and leaf, are considered for the analyses comprising the seven species. The physiological and biochemical parameters were analyzed using FactoMineR software of R Commander (Le et al., 2008).¹ A principal component analysis (PCA) was performed using the mean values of the measured parameters (Supplementary Table S1) along the different experiments according to Table 1. The mean value corresponds to several biological replicates (n) for each species, as follows: for ryegrass, clover, oilseed rape, and Brachypodium n = 4 (where each sample is a pool of ten plants); for wheat n = 3 (each sample being a pool of five plants); for tomato n = 9 or n = 3 (for amino acids); and for Arabidopsis n = 3 (where each sample is a pool of three plants). The data were visualized using PCA scores; this representation was used to show an overview of the entire dataset and to identify associations among the traits measured, as well as among the species. Pearson’s correlations were also run, and representations are shown when considered interesting.

In order to outline common physiological and metabolic traits in response to nitrogen nutrition a dataset of seven plant species has been considered in the present study. All species were grown with ammonium or nitrate as N source, although under different environmental factors that are known to affect ammonium tolerance, such as, medium pH, light intensity, or atmospheric CO₂ concentration (Table 1). It must be noted that, at least, four different growth conditions were included for each species, and both organs, root and leaf, were considered in the study. Firstly, a whole-plant PCA was performed with the mean values corresponding to the physiological and metabolic parameters from both organs for the seven species, which included a total of 80 observations (Supplementary Table S1; Supplementary Figure S1). Nine parameters were considered including biomass, internal NH₄⁺ content, total amino acid content (total AA), soluble protein content and enzymatic activities of GS, GDH, PEPC, NAD-ME, and ICDH. The principal components vectors, PC1 and PC2 explained 60.8% of the total variation (Supplementary Figure S1). In this whole-plant multi-species analysis, the first principal component (PC1) accounted for 41.9% of the total variance, without a clear separation of species along this axis. In contrast, PC2, which accounted for 18.9% of the variance, separated wheat from the rest of the species and came out to be more distant from Arabidopsis, suggesting clear differences in the response to nitrogen between these two species. Thus, GS activity and organ biomass were the markers that defined wheat response, while NAD-ME and GDH mostly defined the response of Arabidopsis.

Leaf and roots can show a different response to nitrogen nutrition, since ammonium is normally assimilated in roots meanwhile nitrate is preferentially assimilated in photosynthetic tissues (Yamaya and Kusano, 2014). Therefore, we considered necessary to carry out a second PCA analyzing separately data from roots (40 observations) and leaves (40 observations; Figure 1). Doing so, the variance explained by the first two PCs greatly increased. PC1 and PC2 together explained up to 74.8% (55.1% PC1 and 19.7% PC2; Figures 1C,D) for the root and to 64.4% (37.8% PC1 and 26.6% PC2; Figures 1A,B) for the leaf. Overall, as for the whole-plant multi-species analysis, the PCA for each organ clearly separated some species (Figures 1B,D). In the individuals representation for roots, only wheat and Arabidopsis remained distant to each other and clearly separated from the rest of species, that were situated in the middle and grouped closer among them (Figure 1D). This analysis indicated that GS and root biomass were driving factors in PC2 for wheat root, while soluble protein, PEPC, ICDH, and NAD-ME mainly contributed to Arabidopsis in PC1. For leaves soluble protein content, GDH and NAD-ME were the main positive markers of Arabidopsis and Brachypodium in the PC1 axis (Figures 1A,B), and both species placed separately from the intermediate group formed by tomato, ryegrass, clover, and oilseed rape. As in the root individuals plot, GS and biomass contributed as biomarkers for species segregation along PC2 in the leaf, being positive markers for wheat. Overall, leaf data showed a more clear separation than root data, which suggested that the leaf response to nitrogen nutrition was wider among species. Then, a subsequent PCA was performed for each species including variables related to N-compounds additional to the nine common variables previously considered. Among additional variables, N content, Asn, Gln, Glu, Asp, and Ser + Gly, as well as ratios of Gln/Glu, Asn/Asp, Asn/Gln, and Glu/Asp were included in the analysis for every species except for Arabidopsis, where only the former nine parameters were considered. The loading-plot when analyzing by PCA individually each species confirmed the separation between root and leaves for the seven species (Supplementary Figures S2, S3). Next, to extract information regarding specific responses to nitrogen nutrition we analyzed roots and leaves separately for each species.

The biplot PCAs for root data indicated that the two first PCs explained from a minimum of 64.9% to a maximum of 96.3% of the total variance (Figure 2). Observations of ammonium- and nitrate-fed root clearly separated along PC1 axis (43.3–80.9% of the variance). This revealed that the nitrogen nutrition imposed a different metabolic a biochemical functioning. In all the species several N-compounds (internal NH₄⁺, total amino acids, Gln, Asn, Ser + Gly, and N content) contributed simultaneously to the metabolic profile of ammonium-fed roots, hence, we consider this group of traits as a root “N-cluster” responding to ammonium nutrition (Figure 2). Moreover, GDH behaved as a driving trait for ammonium-grown roots. Root biomass, in contrast, was a common positive marker for nitrate-fed roots; as well as ratios Glu/Asp or Gln/Asn for four out of six species, which can be taken as indicators of the preference of stored amine or amide in such species. GS activity correlated positively with nitrate-fed root of wheat, Brachypodium, but positively with ammonium-fed roots in the case of tomato and ryegrass. Overall, C-enzymes positioned closely to the “N-cluster” in ryegrass, Brachypodium, and oilseed rape, which suggested that C-enzymes would also be markers for ammonium nutrition. However, the general pattern showed a specific trend depending on each species. In the case of wheat, ICDH and NAD-ME segregated from “N-cluster” (Figure 2), as well as NAD-ME in Arabidopsis (Supplementary Figure S3). In tomato, Asp, Glu, and NAD-ME were separated from the main “N-cluster”. Besides, Asp was positioned closely to biomass, common positive marker for nitrate nutrition. And interestingly, clover was the unique species where PEPC acted as positive marker for nitrate-fed root. Correlations between the different parameters are shown in Supplementary Figure S4.

For leaf dataset, similarly to roots, the two first PCs explained a minimum of 73.5% and a maximum of 93.5% of the total variance depending on the species (Figure 3). Leaves corresponding to both N sources positioned segregated along PC1 axis, which explained between 45.5–76.8% of variance. As for roots, a “N-cluster” including main N-compounds (as internal NH₄⁺, total amino acid, Gln and Asn, and soluble protein) could be defined as driver for ammonium nutrition in the leaf. Moreover, the variables in this leaf “N-cluster” showed overall a close correlation values with each other (Supplementary Figure S4). Interestingly, GDH was always closely related to the leaf “N-cluster” and GS also correlated positively to ammonium nutrition in wheat and clover leaves (Figure 3). On the other hand, leaf biomass was a clear positive marker for nitrate-fed leaves, in the case of Brachypodium and oilseed rape; a response also observed for Arabidopsis (Supplementary Figure S3). The behavior of leaf C-enzymes as biochemical marker for nitrogen nutrition was species-dependent, as their position in the PCA representation differed between species. For instance, although PCA did not revealed a clear separation between N nutrition in clover and ryegrass leaves, PEPC and ICDH were traits more positively correlated to nitrate-fed leaves. Nevertheless, for the rest of species, ICDH positively contributed as biochemical marker to leaves under ammonium nutrition (Figure 3; Supplementary Figure S3). In the case of NAD-ME, it was a positive marker for ammonium-fed leaves of clover, tomato and Brachypodium, but it behaved as a positive marker for nitrate-fed leaves of wheat, and oilseed rape.

Overall, the variation of total and soluble protein behaved as the principal contributors of PC1 in both organs (Figures 2, 3). Specially, the higher contribution of total amino acids to PC1 in Brachypodium, oilseed rape, and Arabidopsis roots (Figure 2; Supplementary Figure S3) would reflect their role as better biomarker of ammonium metabolism than the amount of NH₄⁺ accumulated per se. Moreover, the higher correlation of total amino acid contents with soluble protein and C-enzymes (Supplementary Figure S4; Figures 4, 5C), comparing with NH₄⁺, would also support their role as universal markers. Thus, the storage of NH₄⁺ into amino acids and soluble protein could be considered as a metabolic detoxifying strategy, indicative of ammonium tolerance or the capacity to mitigate ammonium toxicity. The fundamental role of amino acid biosynthesis as principal detoxifying mechanism of NH₄⁺ during plant response to ammonium nutrition remains valid for each species, as shown by the high slope of amino acid accumulation when plotted against NH₄⁺ content in both organs (Supplementary Figure S5). Then, in general if NH₄⁺ absorbed by the cell is not assimilated, its excess becomes toxic different ways for cells (Britto and Kronzucker, 2002) and reduces growth for most plant species. In this sense, the accumulation of inorganic N in form of NH₄⁺ in tissues could be rather considered as an indicator of disturbance of NH₄⁺ homeostasis (ion imbalance), which would trigger further cell or whole-plant responses (Sonoda et al., 2003; Ranathunge et al., 2014; Liu and von Wirén, 2017). Interestingly, considering the seven species together, the highly significant correlations of root amino acid content with some enzymes, as ICDH and GDH, but not with NH₄⁺ content would indicate a possible role of ammonium-reduced compounds in the regulation of C and N-metabolism (Figures 4, 5; Supplementary Figure S4).

Among free amino acids, amides (Asn and Gln) perform an important role as N storage and long-distance N transporters thanks to their high N:C ratio, so they are reported as the primary N-reduced compounds that accumulate when plants are fed with NH₄⁺ (Lea et al., 2007; Vega-Mas et al., 2019a,b). In the studied species, these amides represent, together with their respective precursors Glu and Asp, the major pool of amino acids, summing from 30 to 90% of total amino acid content (Supplementary Table S1). However, the nature of stored amides varies depending on species or families, since plants tend to accumulate either Asn or Gln, but also a mixed trend occurs. In monocots (wheat, ryegrass, and Brachypodium) and clover as well, the ammonium nutrition results in the accumulation of huge amounts of free Asn respect to nitrate nutrition (Supplementary Table S1); likewise accompanied by Gln in the case of ryegrass. The high contribution of Asn to PC1 in ammonium conditions would support its physiological role in ammonium nutrition (Figures 2, 3). Indeed, Asn is usually the main amino acid accumulated in some monocot species under ammonium nutrition (Setién et al., 2013; de Souza Miranda et al., 2016; de la Peña et al., 2019). Wheat, for instance, uses free Asn as N storage in the grain (Gao et al., 2016). In clover, the predominance of Asn under ammonium nutrition is in accordance with previous reports in temperate legumes such as pea or alfalfa (Pasqualini et al., 2001; Lea et al., 2007; Ariz et al., 2013). Interestingly, with the exception of Brachypodium, in Asn-accumulating species the equilibrium between amino acids vs. NH₄⁺ accumulation was quite similar in root and leaf (similar slope; Supplementary Figure S5). This would indicate that the assimilatory capacity of inorganic N depends on NH₄⁺ accumulated, which seems to occur in the same extent in photosynthetic and non-photosynthetic tissues within a given species. On the other side, in the case of tomato, Gln is the amino acid that keeps a better correlation with the internal NH₄⁺ and with a high contribution to PC1, both for roots and leaves (Figures 2, 3; Supplementary Figure S4). Serine would also take importance in tomato metabolism, unlike monocots and clover. In oilseed rape roots, the separation of Asn from the principal cluster including reduced N-compounds (Figure 2) would match with the fact that Brassicaceae species are not specialists in accumulating N reduced in the form of Asn, probably because NH₄⁺ is stored in the form of other amino acids, as Gln or Gly + Ser (Supplementary Table S1) or it can be derived to the synthesis of glucosinolates, N-compounds characteristic of this family (Marino et al., 2016; Coleto et al., 2017). Overall, we can propose that those species that tend to behave as “ammonium excluders” and show a tendency to detoxify NH₄⁺ in the form of Asn are more tolerant to ammonium nutrition (wheat, Brachypodium, ryegrass, and clover). This proposal would agree with the fact that Asn is the most abundant amino acid in phloem and xylem of rice (Hou et al., 2019), a species considered as an ammonium-preferring crop (Zhu et al., 2009; Cao et al., 2018). Besides, Asn represents a 20% saving of carbon with respect to Gln and, as suggested initially by Ikeda et al. (2004), the asparagine synthase would be enhanced in Asn-accumulating species under ammonium nutrition similarly to its induction under abiotic stress (Lea et al., 2007).

Glutamine synthetase enzyme is the connecting point between inorganic NH₄⁺ and organic N. Cytosolic GS (GS1) is responsible for primary NH₄⁺ assimilation as well as reassimilation during germination or senescence, while chloroplastic GS (GS2) isoform is responsible for the assimilation of NH₄⁺ coming from primary nitrate assimilation and photorespiration in leaves (Thomsen et al., 2014). Cytosolic GS1 is the predominant isoform in roots and its importance to confer tolerance when NH₄⁺ is the sole N source for the plant has been reported (Guan et al., 2016; Sarasketa et al., 2016; Konishi et al., 2017). Interestingly, in wheat, the values for GS activity in both organs were very high and distant from other species, hampering the joint analysis of N nutrition response for the ensemble of the seven species (Supplementary Table S1; Figures 2, 3). Accordingly, to the high amino acid contents in roots (Supplementary Table S1; Supplementary Figure S5) we corroborate that NH₄⁺ absorbed is primarily assimilated in this organ for the species studied, in the same line as several studies previously reported (Li et al., 2014; Liu and von Wirén, 2017). This means that root system manages to assimilate NH₄⁺ despite having 4–10 times lower GS activities in comparison to photosynthetic organ in all the species. The exception was Arabidopsis, whose low GS activity showed similar values in both organs, and this could explain why the synthesis of amino acids seems to plateau in this species (Supplementary Figure S5). In general, there was neither a clear response of root GS activity to nitrogen nutrition nor a correlation with NH₄⁺ or total amino acid content (Supplementary Figure S4). Overall, GS activity seems crucial to equip the leaf with high soluble protein or N, except for oilseed rape and Arabidopsis (Figure 3; Supplementary Figures S3, S4). Some authors have reported that GS is a major checkpoint for controlling the N status (Kichey et al., 2006), and plant species with higher GS activities seem more tolerant to ammonium nutrition (Cruz et al., 2006; Li et al., 2011). Among the several isogenes encoding GS, GLN1;2 or GS1;2 is essential for primary ammonium assimilation, mainly in roots, as revealed by the ammonium sensitivity displayed by mutants in this GS isoform (Yamaya and Kusano, 2014; Guan et al., 2016). The highest GS activity would be concordant with the fact that wheat, ryegrass, clover, and Brachypodium seemed to manage better the ammonium nutrition. In the case of wheat, Brachypodium, and clover, GS in leaf was positively correlated with internal NH₄⁺ contents, along with the most of the variables grouped in the principal cluster (Figure 3; Supplementary Figure S4). Similar results were observed by Kichey et al. (2006) for wheat, with a strong relationship between shoot GS activity and total N content (NH₄⁺ and amino acids), regardless the wheat cultivar examined or the N-feeding conditions. Nevertheless, recently, the ammonium toxicity has been attributed to an excessive assimilation of NH₄⁺ in shoots by the plastidic GS (encoded by GLN2), which seems to promote an acidification in the cell, although the expression of root cytosolic GS (GLN1.2) ameliorated the ammonium toxicity (Hachiya et al., 2019). Since GS enzyme is subjected to a complex regulation (at transcriptional, post-transcriptional, and translational levels) with a negative allosteric behavior in the presence of some metabolites (Guan et al., 2016; Németh et al., 2018), and several mechanisms can contribute to regulate the cell pH in cells (Britto and Kronzucker, 2005), more studies are needed to clarify the role that this key enzyme exerts in plants ammonium tolerance, depending on the organ and the species. Future works will help to confirm if constitutive high GS activity, especially in roots, is a metabolic trait involved in ammonium tolerance.

The duality of GDH functioning has been long-discussed. On one hand, functioning in its aminating sense may help preventing ammonium toxicity under stress conditions (Skopelitis et al., 2006) and works synthesizing glutamate during tomato fruit ripening (Ferraro et al., 2012) or during NH₄⁺ reassimilation (Dubois et al., 2003; Grabowska et al., 2011). On the other hand, GDH may function in its deaminating sense for providing C skeletons, notably under C-limiting scenario (Glevarec et al., 2004; Grabowska et al., 2011; Fontaine et al., 2012). GDH activities are considerably higher in the root than in the leaf for the species included in this study, except Arabidopsis (Supplementary Table S1). Overall, the induction of GDH activity, especially in roots, has been previously reported as a response to ammonium nutrition (Lasa et al., 2002; Cruz et al., 2006; Setién et al., 2013; Sarasketa et al., 2016; de la Peña et al., 2019). In the present work, the positive response of GDH to NH₄⁺, amino acids or Gln/Glu and Asn/Asp ratios (Supplementary Figure S4) for each species, would point out the involvement of this enzyme in the overall detoxifying of excess NH₄⁺ into primary N-compounds. Interestingly, there is a positive correlation of GDH and ICDH with the total content of amino acids in the root-system for the ensemble of the seven species (Figure 5; Supplementary Figure S4). Thus, this high correlation could reflect 2-OG provision by ICDH at optimal rates to GDH for Glu synthesis, giving further support for a high Gln synthesis. Alternatively, GDH could be in charge of regulating the C/N ratio and function in the sense of Glu deamination to fuel back the 2-OG, in what can be considered as an anaplerotic entrance of C to the TCA. This deaminating role would make sense in the context of a severe ammonium toxicity whenever the release of NH₄⁺ from Glu would be compensated with the further assimilation of 2 NH₄⁺ molecules as a four-carbon skeletons Asn, thus saving carbon. Therefore, hypothetically the flexibility of GDH enzyme might lean the C provision toward Glu-Gln (if functioning as aminating) or Asp-Asn (if functioning as deaminating) synthesis.

In leaves, despite the inclusion of GDH, in general terms, in the main cluster of variables in the PCA (Figure 3), it shows a weaker response to internal NH₄⁺, accordingly to the lower GDH activity values compared to the root (Supplementary Table S1; Supplementary Figure S4). The exception is Arabidopsis, where the GDH activity is considerably induced in leaves and maintains a high correlation with not only the total amino acid content but also with NH₄⁺ content. Thus, high GDH activity could compensate the low leaf GS activity in this species considered as “ammonium includer,” and could explain that the total amino acids level in Arabidopsis leaves are in the same range than in the rest of species. It is also worth noting the apparently coordinated regulation of GDH and NAD-ME in the photosynthetic organ for the ensemble of seven species (Figure 5; Supplementary Figure S4). The use of GDH activity as an ammonium-tolerance marker really deserves to be considered in plants, regardless of the sense of the GDH reaction itself. Nevertheless, further experimentation is required for understanding the possible co-regulation of GDH with other TCA enzymes and to shed more light about the specificity toward its substrates in vivo.

The derivation of NH₄⁺ into specific amino acids and soluble proteins will depend not only on N-assimilating enzymes but also on C skeletons-providing enzymes. The common response of the dataset of the whole-plant showed that anaplerotic enzymes PEPC, NAD-ME, and ICDH were acting together to support the C skeletons provision for amino acids and soluble protein in response to N nutrition (Supplementary Figure S1). And the plasticity of these anaplerotic enzymes to match the C demand under different N regimes was further highlighted by the organ- and species-specific responses (Figures 2, 3) and correlations between variables (Figure 5; Supplementary Figure S4).

Phosphoenolpyruvate carboxylase has many roles in plant metabolism, one of its primary roles in roots being refixing the CO₂ released during respiration to feed oxaloacetic acid (OAA) into the TCA cycle or to support the Asp+Asn formation. The high contribution of PEPC to PC1 in roots of Gramineae (wheat, ryegrass, and Brachypodium) and Brassicaceae (oilseed rape and Arabidopsis), close to the “N-cluster” (Figure 2), would note its role in the provision of C skeletons to avoid an overaccumulation of NH₄⁺. The induction of PEPC under ammonium nutrition has been similarly described for other species (Lasa et al., 2002; Britto and Kronzucker, 2005; Ariz et al., 2013; Arias-Baldrich et al., 2017). However, PEPC in general is not well correlated with NH₄⁺ content (Supplementary Figure S4), so this metabolite is not probably the factor that directly regulates PEPC activity. Indeed, our data suggest that downstream ammonium-derived products could be really upregulating PEPC, since in the root this activity keeps a strong correlation with amino acids, with Asn/Asp ratio in wheat and Brachypodium, with Gln and Asn in ryegrass, and with Gln in oilseed rape (R > 0.97; Supplementary Figure S4; Figure 4), main indicators of the N status in these plants. Indeed, Gln and Asn were also reported to upregulate PEPC at transcriptional level (Sugiharto et al., 1992; Suzuki et al., 1994). In this sense, the activation of PEPC in barley was dependent on the newly formed Gln, whenever GS is not inhibited (Díaz et al., 1995). Interestingly, despite the high Asn accumulation in the clover root, PEPC activity would not be participating as C anaplerotic enzyme in provision of OAA, if considering the negative correlation of PEPC with Asn (R = − 0.754; Supplementary Figure S4); as well as if considering its opposite position to the main cluster gathering the N-compounds in the PCA (Figure 2). This result in clover agrees with that obtained for alfalfa roots, where nitrate stimulated PEPC more effectively than ammonium (Pasqualini et al., 2001), but differs from the activation reported for this enzyme in pea roots upon ammonium nutrition (Lasa et al., 2002; Ariz et al., 2013). These contrasting responses suggest a more diverse regulation of anaplerotic PEPC in the provision of oxoacids within the Leguminosae family. Given the lack of PEPC responsiveness does not hamper the Asn formation in clover, alternative enzymes, as malate dehydrogenase, could be performing this role.

The presence of NAD-ME activity in mitochondria outfits the TCA cycle with an alternative pathway that allows the decarboxylation of malate rendering pyruvate that can be incorporated to TCA cycle (Wedding, 1989). Non-photosynthetic malic enzymes are key to govern malate-pyruvate balance, metabolites with multiple functions in the cell. Recently, the role of malic enzymes has also been reviewed in plant development and stress responses. For instance, malic enzymes can provide reductant power in the form of NAD(P)H, which are essential to remove reactive oxygen species under stressful scenarios (Sun et al., 2019). The proximate position of NAD-ME to the main “N-cluster” in roots, except for tomato and clover (Figure 2), suggests the coordinated regulation of this activity with other anaplerotic enzymes, PEPC and ICDH, to fuel the synthesis of oxoacids. Concerning mitochondrial isoform NAD-ME, a few papers described this enzyme as showing a positive response in the context of ammonium nutrition in spinach, pea, and tomato roots depending on the C-status of the plant (Kandlbinder et al., 1997; Lasa et al., 2002; Vega-Mas et al., 2015). Thus, the participation of NAD-ME seems to be claimed depending on the C-status of the plants, which is determined by the environmental growth conditions, as it occurs in tomato growing under elevated CO₂, in such a way that under non-limiting C conditions the role of NAD-ME as anaplerotic enzyme may not be essential (Supplementary Table S1).

The enhanced TCA cycle, either by provision of OAA or pyruvate would ultimately promote the fuelling of C through ICDH activity to form 2-OG, the key organic acid regulating the coordination of C and N metabolism (Hodges et al., 2003). The cytosolic NADP-ICDH isoforms represent mostly (90–95%) the total ICDH activity (Chen, 1998; Sulpice et al., 2010; Leterrier et al., 2012). In ryegrass, tomato, Brachypodium, and oilseed rape, ICDH is well correlated with PC1, which leads to consider this enzyme as biomarker of ammonium nutrition, specifically in roots, and working in coordination with NAD-ME (Figure 3). Meanwhile, in wheat and clover, ICDH is more related to the PC2, but with dissimilar tendency in both species. In wheat, ICDH would support the high GS activity, while in clover it is opposed to GS and would be coordinately regulated with NAD-ME to provide oxoacids for amino acids synthesis (Figures 3, 5).

The role of photosynthetic organ to assimilate NH₄⁺ excess was revealed more critical in those species behaving as “ammonium includers.” Overall, the storage of NH₄⁺ into soluble protein in the leaf would be supported by PEPC response for the set of the seven species (Figure 5; Supplementary Figure S4). The positive relation of ICDH to the NH₄⁺ accumulated and Asp or Glu formation (Supplementary Figure S4) in Arabidopsis, a species that can be considered among “ammonium includers”, and in oilseed rape, the weakest “ammonium excluder” among the species studied, would point it out as a key C-enzyme in an attempt to assimilate the NH₄⁺ in the leaf in these species. The assimilation of NH₄⁺ reaching photosynthetic tissues in “ammonium excluders” would be also favored by an enhanced ICDH, surely to endorse the high leaf GS activity characteristic of these species (Supplementary Table S1; Figure 5). With the exception of Arabidopsis, the fact that ICDH activity highly correlated with GS activity (Figure 5) in the multi-species representation would indicate this C enzyme is well-coordinated with N assimilation, probably to sustain the 2-OG production for the enhanced synthesis of Glu + Gln.

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.