Original Research ARTICLE
Short-Term Treatment with the Urease Inhibitor N-(n-Butyl) Thiophosphoric Triamide (NBPT) Alters Urea Assimilation and Modulates Transcriptional Profiles of Genes Involved in Primary and Secondary Metabolism in Maize Seedlings
- 1Dipartimento di Scienze Agroalimentari, Ambientali e Animali, University of Udine, Udine, Italy
- 2Department of Biotechnology, University of Verona, Verona, Italy
- 3Institute of Plant Biology, University of Zurich, Zurich, Switzerland
To limit nitrogen (N) losses from the soil, it has been suggested to provide urea to crops in conjunction with the urease inhibitor N-(n-butyl) thiophosphoric triamide (NBPT). However, recent studies reported that NBPT affects urea uptake and urease activity in plants. To shed light on these latter aspects, the effects of NBPT were studied analysing transcriptomic and metabolic changes occurring in urea-fed maize seedlings after a short-term exposure to the inhibitor. We provide evidence that NBPT treatment led to a wide reprogramming of plant metabolism. NBPT inhibited the activity of endogenous urease limiting the release and assimilation of ureic-ammonium, with a simultaneous accumulation of urea in plant tissues. Furthermore, NBPT determined changes in the glutamine, glutamate, and asparagine contents. Microarray data indicate that NBPT affects ureic-N assimilation and primary metabolism, such as glycolysis, TCA cycle, and electron transport chain, while activates the phenylalanine/tyrosine-derivative pathway. Moreover, the expression of genes relating to the transport and complexation of divalent metals was strongly modulated by NBPT. Data here presented suggest that when NBPT is provided in conjunction with urea an imbalance between C and N compounds might occur in plant cells. Under this condition, root cells also seem to activate a response to maintain the homeostasis of some micronutrients.
It has been estimated that the global consumption of maize will double before 2050, with this cereal becoming the crop with the greatest production worldwide (Rosegrant et al., 2008).
Due to the low nitrogen (N) use efficiency (< 33%; Raun and Johnson, 1999), large amounts of N fertilizer have to be applied to sustain a high productivity of cultivated maize.
Urea represents the major form of N fertilizer applied worldwide (>50%; Heffer and Prud'hommer, 2014). During the last years a large number of studies demonstrated that roots of higher plants, including crops, possess dedicated transport systems for the acquisition of external urea. When present at low concentration in the soil, urea is taken up by plants through DUR3, a specific transporter present in the plasma membrane of root cells. The DUR3 transporter belongs to the sodium:solute symporter (SSS) family and is the major component of the high affinity urea transport system in plants (Kojima et al., 2007). Several aquaporins have been described as components of a low affinity acquisition system, although their role in urea nutrition remains to be clarified (Gu et al., 2012; Yang et al., 2015; for review see Kojima et al., 2006).
It has been suggested that more than the assimilation, the limiting factor for an efficient use of urea in crop plants is the slow uptake (Wang et al., 2013), which in turn depends on urea persistence in the soil. Once in the soil the stability of urea is affected by the presence of urease, an ubiquitous enzyme released by the microbial population or deriving from the decomposition of organic matter. This enzyme catalyzes the rapid hydrolysis of urea releasing carbon dioxide and ammonia, a N form which might be lost from the soil through volatilization (Dawar et al., 2011; Soares et al., 2012).
In order to delay urea hydrolysis, it has been proposed to apply urea in association with urease inhibitors (Zaman et al., 2008). In co-formulation with urease inhibitors, urea is available as an intact molecule for plant acquisition during a long term, acting essentially as a controlled (or slow) release N-fertilizer (Trenkel, 2010).
Several chemical compounds have been described to act as urease inhibitors (Kiss and Simihaian, 2002). Today the most promising and tested urease inhibitor applied in soils is N-(n-butyl) thiophosphoric triamide (NBPT), often formulated into granular urea fertilizers. This compound is a structural analog of urea (Medina and Radel, 1988) and acts with a mixed inhibition on urease activity (Juan et al., 2009). NBPT forms stable complexes with urease (Watson, 2005) coordinating both nickel (Ni) atoms of the urease active site and binding the oxygen atom of the urea-derived carbamate (Manunza et al., 1999).
Although NBPT is widely used, very few information is available on its interaction with the environment, the effects on soil and water quality, on animal health, and also on plant growth (for review see National Industrial Chemicals Notification and Assessment Scheme, 1997; Arora and Srivastava, 2013). Some studies showed that NBPT can affect soil characteristics, such as pH, temperature, and moisture content (Hendrickson and O'Connor, 1987; Hendrickson and Douglass, 1993; Sigunga et al., 2002; Clough et al., 2004), and lead to phytotoxic symptoms (Krogmeier et al., 1989; Bremner, 1995; Watson, 2000).
In a previous study (Zanin et al., 2015a), it was demonstrated that the presence of NBPT in the nutrient solution limited the capacity of plants to use urea as a N-source, resulting in a limited net-uptake rate of urea and in a reduced accumulation of ureic-N in maize seedlings. At least in part, this effect might be the consequence of NBPT uptake by roots. Indeed, in pea and spinach, NBPT was detected in plant tissues when roots were treated with a high concentration of this molecule (Cruchaga et al., 2011). This evidence highlights that NBPT can cross the root plasma membrane (Cruchaga et al., 2011) although the transport mechanism remains to be defined. Plant metabolism was affected by the accumulation of the urease inhibitor within the cells, as shown by a lower endogenous urease activity. This effect in turn determined an increase of the urea internal pool and a strong reduction of the ammonium levels in plant tissues (Artola et al., 2011; Cruchaga et al., 2011; Zanin et al., 2015a). Also the downstream reactions of the cytosolic urea assimilation were affected by NBPT, as indicated by the low glutamine synthetase activity (Artola et al., 2011) and the low expression of genes coding for glutamine and asparagine synthetases (Zanin et al., 2015a). Indeed, these enzymes might play a key role in the assimilation of ammonium released by the cytosolic hydrolysis of urea (as proposed by Mérigout et al., 2008; for review see Pinton et al., 2016).
Often the treatment with the combination of urea and NBPT is associated with leaf tip scorch and necrotic margins of laminar leaves (Krogmeier et al., 1989; Watson and Miller, 1996; Artola et al., 2011; Cruchaga et al., 2011). Although it is unclear if these visible symptoms are indirectly due to a cytotoxic effect of over-accumulated urea, the inhibitor is known to alter the first steps of ureic-N assimilation (Cruchaga et al., 2011; Zanin et al., 2015a). No information is available on possible further influences of NBPT on primary and secondary plant metabolism.
In the present work, maize seedlings were used for a transcriptomic analysis aiming to provide an overview of metabolic changes occurring in plants when NBPT is applied. Changes in the amino acid profiles as well in the urease activity were also monitored.
Materials and Methods
Plant Material and Growth Conditions
Maize seeds (Zea mays L., inbred line PR33T56, Pioneer Hybrid Italia S.p.A.) were germinated over aerated 0.5 mM CaSO4 solution in a dark growth chamber at 25°C. After 3 days, the seedlings were transferred into an aerated hydroponic system in a controlled climatic conditions: day/night photoperiod, 16/8 h; light intensity, 220 μmol m−2 s−1; temperature (day/night) 25/20°C; relative humidity 70 to 80%. After 2 days (5 day old) plants were transferred to a N-free nutrient solution (Control treatment) containing (μM): KCl 5; CaSO4 500; MgSO4 100; KH2PO4 175; NaFe-EDTA 20; H3BO3 2.5; MnSO4 0.2; ZnSO4 0.2; CuSO4 0.05; Na2MoO4 0.05. When present nitrogen was added to nutrient solution in form of urea [0.5 mM CO(NH2)2] (Urea treatment). The urease inhibitor NBPT [N(n-butyl) thiophosphoric triamide, Apollo Scientific Ltd] was applied in conjunction with urea (Urea+NBPT treatment) at a concentration of 0.897 μM, corresponding to 0.5% w/w of urea, that is the amount commonly used in the commercial formulation of NBPT-urea fertilizer. The pH of solution was adjusted to pH 6.0 with potassium hydroxide (KOH).
Physiological and transcriptional analyses were performed on 5-day-old maize plants exposed up to 24 h to the different N treatments. Nitrogen sources were supplied to nutrient solution after 1 h from the beginning of the light phase (T0 = 0 h of treatment). After 8 and 24 h of treatment, pool of six plants for each sample were analyzed immediately for physiological experiments or stored at −80°C until further processing for molecular works.
Molecular Work: Microarray Analyses and Real-Time RT-PCR Validation
RNA extractions were performed using the Invisorb Spin Plant RNA kit (Stratec Molecular) as reported in the manufacturer's instructions. Maize roots (70 mg) were homogenized in liquid nitrogen and the powder was mixed with 900 μl of DCT solution and dithiothreitol according to the suppliers' instructions. The RNA was evaluated in an agarose/ formaldehyde gel and quantified using a spectrophotometer Nanodrop 2000 instrument (Thermo Scientific). For the microarray analysis, the RNA quality and quantity were determined using a Bioanalyzer Chip RNA 6000 series II (Agilent) and three independent biological replicates were used, for a total of 6 hybridizations. The cDNA synthesis, labeling, hybridization and washing reactions were performed according to the NimbleGen Arrays User's Guide (www.nimblegen.com). Each hybridization was carried out on a NimbleGen microarray (maize chip 12 × 135K Arrays from Roche), representing 59756 transcripts predicted from the B73 reference genome (ftp.maizesequence.org/current/filtered-set/ZmB73_5b_FGS_cdna.fasta.gz). A complete description of the chip is available at the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo) under the series entry (GPL17540). The microarray was scanned using an Axon GenePix 4400 (Molecular Devices) at 532 nm (Cy-3 absorption peak) and GenePix Pro7 software (Molecular Devices) according to the manufacturers' instructions. Images were analyzed using NimbleScan v2.5 software (Roche), which produces Pair Files containing the raw signal intensity data for each probe and Calls Files with normalized expression data- (quantile normalization) derived probe summarization through RMA analysis (Irizarry et al., 2003). Analysis of normalized data (Calls Files) was performed using the open source software of the Bioconductor project (Gentleman et al., 2004) with the statistical R programming language (Ihaka and Gentleman, 1996). Differentially expressed probes were identified by linear model analysis (Smyth, 2005) using the LIMMA package and applying Bayesian correction, adjusted P ≤ 0.05, n = 3, FC ≥ |2.00|. All microarray expression data are available at the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo) under the series entry GSE76828 and GSE53102.
Gene ontology (GO) analysis and GO enrichment was performed using the Singular Enrichment Analysis (SEA) of AgriGO (Du et al., 2010) with Zea mays 6a (http://genome.jgi.doe.gov/; Schnable et al., 2009) as custom reference background. Hypergeometric tests with Yekutieli as multi-test adjusted method were performed using the default parameters to adjust the P-value.
Visualization of transcript expression differences within specific pathways was carried out using MapMan v.3.5.1R2 (Usadel et al., 2009) with the mapping file provided by the MapMan homepage (Zmays_181; http://mapman.gabipd.org/) and custom files for the pathways: urea assimilation and primary metabolism (Supplementary Table 1, Supplementary Figure 1) and phenylpropanoid pathway (Supplementary Table 2, Supplementary Figure 2).
To validate the microarray results, real-time RT–PCR analyses were performed. Total RNA was treated with 1 U mg−1 RNA of DNase I (Sigma Aldrich) and cDNA was synthesized from 1 μg of RNA following the application protocol of the manufacturers [42°C for 1 h with 1 pmol of oligo d(T)23VN (Sigma Aldrich); 15 U of Prime RNase Inhibitor (Eppendorf); 10 U of M-MulV RNase H− (Finnzymes)]. After RNA digestion with 1 U of RNase A (USB) for 1 h at 37°C, gene expression analyses were performed by adding 0.16 μl of the cDNA to the real-time RT-PCR complete mix, FluoCycle™ sybr green (20 μl final volume; Euroclone, Pero, Italy), in a DNA Engine Opticon Real Time PCR Detection system (Biorad). Real-time RT-PCRs for 15 genes identified in the microarray comparisons were performed (Supplementary Table 3). The primers were designed using Primer3 software (Koressaar and Remm, 2007; Untergasser et al., 2012) and they were synthesized by Sigma Aldrich (Supplementary Table 3). The analyses of real-time results were performed using Opticon Monitor 2 software (Biorad) and the qPCR package (version 1.1–8; www.dr-spiess.de/qpcR.html) for the statistical R software (version 2.9.0). Efficiencies of amplification were calculated following the authors' indications (Ritz and Spiess, 2008). Real-time RT-PCR results were validated using two reference genes, on ZmGPDH and ZmTUA. Data were normalized with respect to the transcript level of the reference genes using the 2−ΔΔCT method, where ΔΔCT = (CT, Target − CT, HK)Time x − (CT, Target − CT, HK)Time 0 (Livak and Schmittgen, 2001).
Measurement of Nitrogen, Urea, and Ammonium Content in Maize Tissues
After 24 h of treatment, samples of shoots and roots of maize were dried and their total N content was determined using a Carlo-Erba CHN analyser. For urea and ammonium determination, leaves and roots of maize were sampled and processed as described by Witte et al. (2002). The urea content was quantified using the diacetyl monoxime and thiosemicarbazide reagents and measuring the absorbance at 527 nm. The ammonium quantification was performed using the Barthelot reagent (EN ISO 11732) on a San++ Autoanalyzer (Skalar), the absorbance was determined at 660 nm.
Analyses of Amino Acid Profiles
Plant samples (150 mg) were ground in liquid N2 and homogenized in 600 μl of water:chloroform:methanol (3:5:12 v/v; Hacham et al., 2002). The samples were centrifuged at 18000 g for 5 min. The supernatant was collected and 750 μl water:chloroform (4.5:3) was added and the samples were centrifuged at 18000 g for 5 min. The supernatants were pooled into a fresh tube, completely dried at the speed vacuum and solubilized in 100 μl of water. The samples were transferred to liquid chromatography (LC) vials and analyzed by ultraperformance LC-mass spectrometry (UPLC-MS). Amino acids analysis was performed using an UPLC (Thermo Scientific Dionex UltiMate 3000) coupled to a Bruker Compact Electrospray Ionization-Quadrupole-Time-of-Flight (ESI-Q-TOF; Bruker Daltonics). Liquid chromatography separation was performed on a BEH Amide column (1.7 μm, 2.1 × 150 mm, Waters) as described (Guo et al., 2013). The mobile phase was composed of solvent A (water, 10 mM ammonium formate, 0.15% formic acid) and solvent B (acetonitrile, 2 mM ammonium formate, 0.15% formic acid) with a gradient elution: 0–6 min, 15–20% A; 6−10 min, 20–30% A; 10–12 min, 30–40% A; 12–18 min—equilibration to initial conditions. The flow rate was set up to 0.3 ml min−1 and 5 μl of each sample was injected. Electrospray ionization source was operated in positive mode and parameters were set as follow: gas temperature, 220°C; drying gas, 9 l min−1; nebulizer, 2.2 Bar; capillary voltage, 4500 V; end plate offset, 500 V. The instrument was set to acquire m/z 50–1300. The solvents used were LC-MS/MS grade (Chemie Brunschwig) and all amino acids were purchased from Sigma Aldrich (Switzerland).
All data were analyzed using Data Analysis (version 4.2) and TargetAnalysis (version 1.3; Bruker Daltonics). Absolute amino acid quantification was based on standard curves of all analyzed amino acids between 0.25 and 5 μg. This analysis was performed using QuantAnalysis software (version 2.2; Bruker Daltonics).
Determination of Urease Activity and Total Protein Content
To measure the urease activity, shoots, and roots of maize were sampled and processed as described by Witte et al. (2001). The assay was performed using the urease activity assay kit (MAK120, Sigma Aldrich) as reported in the manufacturer's instructions. Total protein was determined using a commercial Bradford assay (Bio-Rad Laboratories, Hercules, CA) with bovine serum albumin as standard.
Physiological and transcriptional analyses were performed on three independent biological replicates obtained from independent experiments (n = 3); for each sample a pool of six plants was used. Statistical significance was determined by one-way analysis of variances (ANOVA) using Student-Newman-Keuls test (P < 0.05, n = 3). Statistical analyses were performed using SigmaPlot Version 12.0 software.
Statistical analysis of microarray data was performed using linear model analysis (Smyth, 2005) of the LIMMA package after Bayesian correction with Bioconductor software, adjusted P ≤ 0.05, n = 3, FC ≥ |2.00| (for details see the “Microarray Analyses” Section).
Nitrogen, Urea, Ammonium and Protein Contents in Plants
In this set of experiments maize seedlings were treated for 24 h with three different nutritional conditions: nutrient solution without any source of N (Control treatment), nutrient solution containing urea (Urea treatment), nutrient solution containing urea and NBPT (Urea+NBPT treatment). No significant changes in dry weight of shoots and roots was observed among the treatments. However, already 24 h after the application of urea, changes in the N distribution were detectable. In fact plants exposed to urea, either treated with or without NBPT, accumulated higher amounts of N in shoots than in roots, while an opposite behavior was observed in Control plants.
Urea concentration increased both in shoots and in roots when plants were supplied with urea (Urea treatment), while ammonium concentration increased only in roots of these plants. Shoots of plants treated with Urea+NBPT contained significantly less N than Urea-treated plants (Table 1). The addition of NBPT caused a significant increase in urea concentration, this effect being especially evident in shoots. On the other hand, ammonium concentration decreased in root of Urea+NBPT treated plants with respect to Urea-treated plants.
Table 1. Effect of NBPT on the total nitrogen (N), urea, ammonium and protein contents in shoots and roots after 24 h of treatment.
Protein content significantly increased in shoots of plants treated with urea, irrespective of the presence of NBPT, while no significant differences among treatments were observed in root tissues (Table 1).
Urease activity was measured in shoots and roots after 8 and 24 h from the beginning of the treatment. No differences in root enzyme activity was observed after 8 h of plant exposure to the different nutritional conditions (Figure 1). In shoots, only a slight increase in the enzymatic activity was observed when NBPT was added to the nutrient solution (Figure 1). Prolonging the treatment from 8 to 24 h, the ureolitic activity increased, both in shoots and roots of maize plants. This increase was significantly reduced when NBPT was added to urea containing solution, with an inhibition as high as 50% in shoots (Figure 1).
Figure 1. Effect of NBPT on urease activity in shoots and roots of maize at 8 and 24 h of treatment. Five-day-old maize plants were exposed for 24 h to a nutrient solution supplied with 0.5 mM urea in presence or absence of NBPT (Urea+NBPT treatment or Urea treatment, respectively) or exposed to a nutrient solution without addition of any N source (Control treatment). DW, dry weight; FW, fresh weight. Different letters indicate statistically significant differences (Student–Newman–Keuls method ANOVA, n = 3, P < 0.05).
Amino Acid Profiles in Roots and Shoots
In order to better investigate the effect of NBPT on N metabolism, we determined the concentration of six amino acids at two time points (8 and 24 h) during the nutritional treatment. The amino acids were chosen due to their key role in the primary assimilation of ureic-N (glutamine, Gln; glutamate, Glu; asparagine, Asn; aspartate, Asp) and their central function in the secondary metabolism (methionine, Met; tyrosine, Tyr). Among them, the most represented amino acids were Gln, Glu, and Asn, while the concentrations of Asp, Met, and Tyr were lower both in roots and shoots (Figure 2). Eight hours after beginning the exposure to Urea or Urea+NBPT, a decrease in Gln, Glu, and Asn concentrations was observed in shoots. This effect was significant in the presence of NBPT (Figure 2A). Conversely, at the same time point (8 h) the Glu levels increased in roots treated with NBPT, possibly indicating a redistribution of this amino acid within the plant (Figure 2C).
Figure 2. Individual profiles of six amino acids in shoots (A,B) and roots (C,D). Amino acid amounts were measured on the same plants as described for Figure 1 at 8 h (A,C) and 24 h (B,D) of treatment. The values are means + SD of three replicates. Different letters indicate significant differences between treatments (Student–Newman–Keuls method ANOVA, n = 3, P < 0.05). AA, amino acid; FW, fresh weight.
Prolonging the treatment up to 24 h, the levels of amino acids in shoots were comparable in Control, Urea and Urea+NBPT treated plants. Similarly, also in roots the amounts of Glu, Asn, and Asp were not affected by the treatment. On the other hand, the root concentration of Gln decreased when plants were treated with NBPT (Figure 2).
Concerning Asp, Met, and Tyr concentration no significant changes were observed among treatments. The only exception was observed for Tyr in roots of urea treated plants, since slight changes were detected at both time points when compared to Control roots. However, the addition of NBPT in the urea containing solution did not significantly affect the content of this amino acid (Figure 2).
Transcriptomic Changes Induced by NBPT in Roots of Urea-Fed Plants
To evaluate if physiological changes could be linked to modulations of gene expression, transcriptomic analyses were performed in maize roots after exposure for 8 h to Urea+NBPT or Urea treatment. Differentially expressed genes were identified by Linear Models for MicroArray (LIMMA, Smyth, 2005), as the subset of genes displaying at least a 2-fold change in transcript abundance (FC ≥ |2.00|, P ≤ 0.05). The analyses revealed that when NBPT was added to the nutrient solution, 1002 transcripts were significantly upregulated, while 1624 were downregulated, Supplementary Table 4).
According to the terms of biological processes of Gene Ontology (GO) the most representative functional categories in the Urea+NBPT vs. Urea comparison were “metabolic process,” “cellular process,” “biological regulation,” “regulation of biological process,” “establishment of localization,” “localization,” and “response to stimulus” (Figure 3). The GO categories “Metabolic process,” “biological regulation,” and “regulation of biological process” were found to be enriched in the Urea+NBPT vs. Urea comparison with respect to the whole maize transcriptome (used as background reference for Singular Enrichment Analyses, SEA Figure 4). On the other hand “signaling,” “cellular component organization,” “signaling process,” “cell wall organization and biogenesis,” and “cellular component biogenesis” categories were less abundantly represented in the data set Urea+NBPT vs. Urea (< 1%, Figure 3).
Figure 3. Functional category distribution of differentially expressed transcripts in Urea+NBPT vs. Urea, according to the terms of biological process of Gene Ontology (GO).
Figure 4. GO-term enrichment analysis of the significantly regulated genes in Urea+NBPT vs. Urea. Significantly most enriched GO-terms for biological process GOs are plotted according to increasing enrichment of the percentage of genes in the data set of Urea+NBPT vs. Urea (black bars) compared to that of the maize reference set (Background, gray bars). Biological Process GOs—GO:0050896 response to stimulus, GO:0051179 localization, GO:0051234 establishment of localization, GO:0050789 regulation of biological process, GO:0065007 biological regulation, GO:0009987 cellular process, GO:0008152 metabolic process. (SEA analysis using AgriGO toolkit: hypergeometric tests with Yekutieli as multi-test adjusted method, P < 0.05).
The analyses of the genes involved in the metabolic, regulatory and transport pathways might allow improving the comprehension of physiological changes occurring in NBPT-treated plants. Therefore, to identify which pathways were affected by the NBPT treatment, up- and down-regulated transcripts were classified into different groups according to the bin divisions of the MapMan software (Version 3.6.0RC1; Thimm et al., 2004; Figure 5, Supplementary Figure 3).
Figure 5. Transcriptional modulation of genes involved in the cell metabolism in Urea+NBPT vs. Urea roots. Color scale refers to the fold change values of differentially expressed transcripts: red color refers to those transcripts positively regulated by Urea+NBPT treatment, while in green are transcripts negatively regulated by Urea+NBPT treatment.
Some genes involved in the pathway of urea acquisition were found to be modulated by the presence of NBPT (Figure 6). Compared to Urea roots, the treatment with Urea+NBPT upregulated a gene coding for urease (# 202 GRZM2G461569_T03, Table 2). This is an alternative transcript of the urease gene and covers only 20% of the complete coding sequence (primary urease transcript GRZM2G461569_T01). Thus, it appears that this truncated mRNA (# 202) encodes a non-functional urease (for functional annotation see: www.phytozome.jgi.doe.gov).
Figure 6. Transcriptional changes of genes involved in the urea assimilation and primary metabolism in Urea+NBPT vs. Urea roots. PDC, pyruvate decarboxylase; NADH-ME, NADH malic enzyme, MDH, malate dehydrogenase, CA, carbonic anhydrase; Acon, aconitase; GOGAT, glutamate oxoglutarate aminotransferase, GDH, glutamate dehydrogenase; NiR, nitrite transferase; NR, nitrate reductase; NRTs, nitrate transporters; AMTs, ammonium transporters; GS1, glutamine synthetase 1; AsnS, asparagine synthetase; Ure, urease; PIPs, plasma membrane intrinsic proteins; TIPs/NIPs, tonoplast intrinsic protein/ NOD26-like intrinsic protein; DUR3, high affinity urea transporter. Color scale refers to the fold change values of differentially expressed transcripts: red color refers to those transcripts positively regulated by Urea+NBPT treatment, while in green are transcripts negatively regulated by Urea+NBPT treatment.
Table 2. List of modulated transcripts involving in urea acquisition and reported in the Results and Discussion Sections by the comparison of transcriptional profiles of Urea+NBPT roots with profile of Urea roots (Urea+NBPT vs. Urea comparison).
Our results also show that several genes closely related to urea assimilation are downregulated: three transcripts coding for carbonic anhydrases (# 93–95, Table 2), three alternative transcript coding for a glutamine-dependent asparagine synthetase (# 185–187, Table 2) and a transcription factor ZFP16-1 that has been reported to be early responsive to urea nutrition (# 543, Table 2). A similar trend was also observed for some genes coding for aquaporins (nodulin 26-like intrinsic protein, NIP, and tonoplast intrinsic proteins, TIPs, # 1356-1358, Table 2) that might be involved in the low affinity transport of urea. On the other hand, one aquaporin that is putatively localized in the plasma membrane (# 1355 Table 2), as well as transcripts coding for nitrate transport and assimilation (# 183, 1283, 1295, 1338, Table 2), were upregulated by Urea+NBPT (Figure 6).
As regards primary metabolism, a considerable downregulation was observed for genes involved in glycolysis and fermentation (PFK, enolase, G-6P isomerase, pyruvate decarboxylase and a putative phosphoenolpyruvate carboxykinase, GRMZM2G580389, # 77–83, 85, and 1704, Supplementary Table 4) as well as in the TCA cycle (ATP citrate lyase, malate dehydrogenase, carbonic anhydrase, # 90–95, Supplementary Table 4, Figures 5, 6). A similar expression pattern was also observed for genes related to the photosynthetic pathway, including light reactions (# 1–54, Supplementary Table 4), photorespiration (# 55–57, Supplementary Table 4), and Calvin cycle (RubisCo, TPI, FBPase, # 58–60, Supplementary Table 4, Figure 5). This modulation was further supported by a downregulation of some plastidic triose-phosphate translocators (TPTs, # 1318–1319, Supplementary Table 4), indicating a crosslink between Calvin cycle and glycolysis.
Despite glycolysis was negatively affected by NBPT, the upstream reactions for the starch degradation (beta-amylase, heteroglycan glucosidase, # 63–65, Supplementary Table 4) were induced in Urea+NBPT roots. In a similar way, also two genes encoding for key enzymes of TCA cycle (aconitase and NADH-malic enzyme, # 88–89, Supplementary Table 4) and one for the translocation of malate into mitochondria (dicarboxylated transporters, # 1324–1325, Supplementary Table 4) were upregulated (Figure 6).
A consistent modulation was also detected for genes related to the electron transport chain: the expression of six transcripts were upregulated by the NBPT treatment (NADH dehydrogenase, cytochrome c oxidase, and ATP synthase, # 96–97, 101–102, 105–106, Supplementary Table 4), while three genes coding for cytochrome c reductase, a cytochrome c, and a flavoprotein were found to be downregulated (# 98–100, Supplementary Table 4, Figure 6).
Concerning the secondary metabolism, some genes involved in arginine catabolism and polyamine pathway were also modulated. In particular a gene coding for an arginine decarboxylase (ADC1, # 201, Supplementary Table 4) was downregulated by Urea+NBPT, as well as several other genes encoding enzymes for polyamine degradation: polyamine oxidase (PAO1, PAO3) and copper amine oxidase (PAO and CuAO, # 372–373, 437, 1499, and 205, Supplementary Table 4). Along with putrescine, S-adenosyl methionine (SAM) participates as methyl donor to the synthesis of spermidine and spermine, a particular class of polyamines. The treatment with Urea+NBPT led to the upregulation of a gene for the synthesis of SAM (# 188, Supplementary Table 4), but limited the expression of some SAM methyltransferases (# 189, 343–344, 1534, 1568, 2048, Supplementary Table 4), putatively involved in other SAM-consuming pathways. Beside these changes linked to arginine catabolism, results showed that three transcripts encoding proline oxidase (# 198–200, Supplementary Table 4) were upregulated by Urea+NBPT.
Another set of transcripts involved in the secondary metabolism was related to the synthesis of phenylpropanoids (# 193–195, 225, 227, 232–238, 250–258, 289, 412, Table 3, Figure 7). In particular, the NBPT treatment downregulated the expression of three transcripts related to the synthesis of shikimate and arogenate (# 193–195, Table 3), while the downstream reactions for the synthesis of phenylpropanoids were upregulated (# 227, 232–235, 237–238, 250–251, 255–256, 258, 289, 412, Table 3), with the exception of some transcripts involved in the lignin synthesis, that were downregulated (# 225, 236, 252–254, Table 3).
Table 3. List of modulated transcripts involving in the phenylpropanoid synthesis and reported in the Results and Discussion sections by the comparison of transcriptional profiles of Urea+NBPT roots with profile of Urea roots (Urea+NBPT vs. Urea comparison).
Figure 7. Transcriptional modulation of genes involved in the phenylalanine/tyrosine derived secondary metabolite biosynthesis in Urea+NBPT vs. Urea roots. In maize plants, phenylalanine, and tyrosine are precursor for the synthesis of several secondary metabolites. Color scale refers to the fold change values of differentially expressed transcripts: red color refers to those transcripts positively regulated by Urea+NBPT treatment, while in green are transcripts negatively regulated by Urea+NBPT treatment. PAL, phenylalanine ammonia-lyase; C4H, cinnamate-4-hydroxylase; 4CL, 4-coumarate CoA ligase; CAD, cinnamoyl-alcohol dehydrogenase; F5H, ferulate 5-hydroxylase; C3H, coumarate 3-hydroxylase; ALDH, aldehyde dehydrogenase; CCR, cinnamoyl-CoA reductase; HCT, hydroxycinnamoyl-Coenzyme A shikimate/quinate hydroxycinnamoyltransferase; CCoAOMT, caffeoyl/CoA-3-O-metheltransferase; CHS, chalcone synthese; CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; F3′H, flavonoid-3′-hydroxylase; F3GT, flavonoid-3-O-glycosyltransferase; FS, flavone synthase; FOMT, flavonoid O-methyltransferase; FCGT, flavone-C-glycosyltransferase; FLS, flavonol synthese; F3GT, flavonoid-3-O-glycosyltransferase; DFR, dihydroflavonol reductase; ANS, Anthocyanidin synthese; AGT, Flavonoid-O-glycosyltransferase; AAT, anthocyanin acyltransferase; BAN, oxidoreductase|dihydroflavonol reductase like; LAC, laccase. The image is adapted from Tohge et al. (2013).
The presence of NBPT in the root external solution also modulated the expression of genes related to biotic (# 305–332, Supplementary Table 4) and abiotic stresses (# 333–351, Supplementary Table 4), this latter group of transcripts being mainly downregulated by NBPT. As a general response to a stress condition we observed the induction of some genes involved in callose synthesis (glucan synthases, # 71–75, Supplementary material Supplementary Table 4) and coding for cell wall proteins (AGPs and LRRs, # 114–115 and 117–118, Supplementary Table 4). On the other hand some genes related to cell wall degradation (xylosidase and glucanases, # 119–124, Supplementary Table 4) and modification (pectinesterases, # 145–149, Supplementary Table 4) were downregulated by NBPT (Figure 5).
This transcriptional behavior might be closely related to the hormonal response of root cells. Indeed some genes linked to auxin (# 263–274, Supplementary Table 4) and ethylene metabolism (# 283–295, Supplementary Table 4) were mostly negatively modulated in Urea+NBPT treated root cells. A different behavior was observed for genes related to jasmonic acid biosynthesis (# 298–301, Supplementary Table 4), which were upregulated by the treatment with the urease inhibitor.
Finally, a miscellaneous group of transcripts, mainly coding for transporters and transcription factors involved in the acquisition of other nutrients, was found to be modulated by NBPT treatment. Most of them were related to the acquisition and homeostasis of Fe. Indeed two transcripts coding for putative Fe(III)-PS transporters (ZmYSL7, # 1333, 1335; Table 4), one coding for the vacuolar storage of Fe (VIT1, # 219; Table 4), one coding for the plastid protein uptake machinery (TIC21, # 868; Table 4) and some transcripts coding for bHLH transcription factors (ZmPYE-like, # 673–674, 2026 and ZmbHLH29, # 684; Table 4) were modulated by the Urea+NBPT treatment.
Table 4. List of modulated transcripts involving in the acquisition of other nutrients and reported in the Results and Discussion sections by the comparison of transcriptional profiles of Urea+NBPT roots with profile of Urea roots (Urea+NBPT vs. Urea comparison).
Exposure of maize roots to Urea+NBPT influenced also the expression of genes putatively involved in micronutrient complexation: a gene encoding a putative nicotianamine transporter (ZmZIF1, # 2144; Table 4) and a gene mediating the efflux of citrate into the rhizosphere (ZmMATE1, # 1361; Table 4) were negatively affected by Urea+NBPT.
N-(n-butyl) thiophosphoric triamide (NBPT) is the most widely used urease inhibitor (Sanz-Cobena et al., 2012; Abalos et al., 2014). However, only few information is available on the effect of NBPT on N acquisition mechanisms in plants.
An inhibitory effect of NBPT on urea uptake and assimilation was shown in a previous work (Zanin et al., 2015a). In the present work, a transcriptomic analysis was performed in order to highlight the modifications induced by a short-term (8 h) treatment with NBPT in the primary and secondary metabolism of urea-fed maize seedlings.
Urea and Ammonium Assimilation Pathways
It is well-known that the acquisition of external urea is mediated by several proteins localized on the root plasma membrane. A specific transporter, DUR3, and some members belonging to the aquaporin family are considered to be the major components of the high- and low-affinity-transport system in plants, respectively (for review see Kojima et al., 2006).
Confirming previous evidence (Zanin et al., 2015a), transcriptomic analysis revealed that the expression of DUR3 gene was not affected by the NBPT treatment, while a different behavior was observed for several aquaporins. Among the modulated transcripts, one sequence (ZmTIP4;1, # 1356 Table 2) exhibited a high homology to AtTIP4;1 that encodes a tonoplastic aquaporin of arabidopsis permeable to urea and inducible by N deficiency in roots (Liu et al., 2003). When the urease inhibitor NBPT was added to the external solution, root cells accumulated urea (Table 1); this in turn would modulate positively or negatively the expression of different aquaporins (# 1355–1358, Table 2) in order to equilibrate the concentration of urea within different cellular compartments (Figure 6).
In the cytosol, urea should be rapidly hydrolyzed by urease, releasing ammonium and carbonic dioxide but, as reported in wheat and pea plants (Artola et al., 2011; Cruchaga et al., 2011), the treatment with the urease inhibitor NBPT interferes with the assimilation of urea and decreases plant urease activity. Also in maize, the NBPT treatment compromised the urease activity (Figure 1).
After 8 h of treatment, the urea concentration already increased in roots (Zanin et al., 2015a), although the enzyme activity was not yet affect by NBPT (Figure 1). At the same time, microarray analyses revealed that the Urea+NBPT treatment positively modulates the expression of a truncated transcript (5′-smallRNA) putatively coding for a partial and inactive fragment of urease (# 202 Table 2), while the expression of the complete transcript coding for a functional urease was not influenced by NBPT (Zanin et al., 2015a). In response to changes in urea concentration, plant cells might regulate the urease activity through a transcriptional attenuation mechanism, as reported in bacteria. In fact, in Helicobacter pylori, a urease-regulatory mechanisms mediated by 5′-smallRNA downregulates the urease expression through a premature termination of the transcript (Wen et al., 2013).
According to the proposed pathway for ureic-N assimilation (Mérigout et al., 2008; Pinton et al., 2016), ammonium derived from urea catabolism is rapidly assimilated by a cytosolic pathway involving a glutamate dehydrogenase, a glutamine synthase and an asparagine synthetase. In the presence of Urea+NBPT, the root capability to hydrolyze urea (Figure 1) decreased reducing the amount of ammonium available for the primary assimilation (Table 1); this result is confirmed by the low amounts of Gln measured in the roots (Figure 2D). Consistent with these observations, the microarray experiment identified the downregulation of genes encoding asparagine synthetase and carbonic anhydrase (# 185–187 and 93–95, Table 2, Figure 6), two enzymes consuming the ureolitic products (ammonium and carbon dioxide). Furthermore, after 8 h of treatment, the amounts of Glu, Gln, and Asn in shoots were reduced by NBPT (Figure 2A); in roots, Glu content increased while Gln or Asn did not exhibit any significant change (Figure 2C). The root accumulation of Glu supports the idea that NBPT might limit ureic-N assimilation. It is plausible that the NBPT treatment would determine a rapid redistribution of amino acids from shoots to roots.
Other genes involved in N acquisition were modulated when plants were treated with Urea+NBPT. Since the nutrient solution did not contain nitrate and no measurable urea degradation occurred in the external medium, the upregulation of genes coding for putative nitrate transporters and for nitrite reductase (# 1283, 1295, 1338, 183, Table 2, Figure 6) might be activated to compensate for the reduced N assimilation. All these changes involved in the urea assimilation pathway might be triggered by the action of the transcription factor ZFP16-1 (# 543, Table 2), a homolog to ZAT12 in arabidopsis, whose gene expression is known to be responsive to urea nutrition (Mérigout et al., 2008; Zanin et al., 2015b) and to abiotic stresses (Davletova et al., 2005). As an early responsive element (Zanin et al., 2015a), this transcription factor might play a key role in activating the pathway for urea assimilation and in turn the inducible acquisition of urea in plants. The downregulation of ZFP16-1 by NBPT might limit the overall mechanism of ureic-N assimilation and redistribution in plants.
In roots of Urea+NBPT treated plants, many genes involved in primary metabolic pathways, like glycolysis and TCA cycle (# 77–83, 85, 90–95, and 1704, Supplementary Table 4, Figures 5, 6), were downregulated. However, a few transcripts, such as malate/dicarboxylate transporters (for the mitochondrial translocation of malate), a NADPH-malic enzyme and aconitase, were induced by NBPT (# 88–89, 1324–1325, Supplementary Table 4). Aconitase converts citrate to isocitrate, a precursor of α-ketoglutarate, which is a key metabolite for N assimilation. These modulations in transcript levels might indicate the activation of an alternative pathway to sustain the TCA cycle through translocation of cytosolic malate into mitochondria via an antiport with Pi, rather than via a pyruvate/OH− antiport (Taiz and Zeiger, 2002, Figure 6). In this way, the dissipation of the transmembrane pH gradient through the mitochondria inner membrane is avoided. Moreover, in the matrix, malate becomes the substrate of malic enzyme which generates pyruvate and, in turn, participates (as acetyl-CoA) in the TCA cycle. The decarboxylation of malate leads to the consumption of H+ and to the release of a CO2 molecule (Figure 6). This alkalization could be potentially compensated by the transformation of CO2 into bicarbonate, a reaction catalyzed by carbonic anhydrases. However, the transcriptomic analyses revealed that these latter enzymes were downregulated by NBPT (# 93–95, Supplementary Table 4, Figure 6), as a possible response of root cells to preserve the alkalization of mitochondrial matrix.
Thus, in the presence of Urea+NBPT, root cells seem to activate an alternative pathway to limit mitochondrial alkalization of the intermembrane space and, at the same time, to increase the pH of the matrix, as a mechanism to enhance the transmembrane potential suitable for ATP synthesis. In fact there were some positive modulations in the expression of genes coding for mitochondrial proteins involved in ATP synthesis and electron transport (complex I and IV). This response might be a strategy to compensate the wide downregulation occurring in the plastidial electron transport chain (# 1–54, Supplementary Table 4, Figure 6).
Secondary Metabolic Pathway: Changes Triggered by the Urease Inhibitor NBPT
Arginine and Polyamine Synthesis
Urea, either derived from the rhizosphere or from internal N cycling (via arginine catabolism), plays an important role as N metabolite in plants (Witte, 2011). In addition to its role in urea synthesis, arginine is a versatile amino acid and serves as a precursor for the synthesis of polyamines, glutamate, proline and participates also in the biosynthesis of several alkaloids (nicotine, tropane-, and pyrrolidine-alkaloids; Sato and Yamada, 2008).
It is known that the polyamine (PA) homeostasis in plants correlates with several important physiological functions, including the control of the C:N ratio (Mattoo et al., 2006), stress response (Alcázar et al., 2011) and protein regulation (Baron and Stasolla, 2008; Tisi et al., 2011a,b). To understand if the plant response to the urease inhibitor NBPT might have consequences on N cycling, we extracted from the gene list the transcripts encoding proteins that are related to arginine and PAs metabolism. So far two different genes coding for isoforms of arginine decarboxylase (ADC1 and ADC2), the first enzyme for the polyamine synthesis, have been identified in the plant genome. The arginine decarboxylation is a limiting step for PA synthesis; in fact the overexpression of ADC2 enzyme isoform has been shown to promote putrescine accumulation (Alcázar et al., 2005). Our transcriptomic data revealed a downregulation for ADC1 (ADC1 # 201, Supplementary Table 4) but not for ADC2 when plants were treated with NBPT. These two enzymatic isoforms show different expression patterns depending on the nature of the stress and appear to play specific roles in response to different stresses (Alcázar et al., 2010). While the expression of ADC1 is highly responsive to cold conditions, ADC2 is induced by a wide range of stresses, such as drought, wounding/jasmonate treatment, salinity, and potassium deficiency (Urano et al., 2003; Cuevas et al., 2008; Alcázar et al., 2010). Moreover, the arabidopsis homolog of the urea-responsive transcription factor (ZFP16-1 # 543, Table 2, see above) acts as a positive regulator of ADC (Vogel et al., 2005). Therefore, the low expression of ZFP16-1 transcript in Urea+NBPT treated maize roots correlated with the downregulation of ADC1 gene. The effect of NBPT on PAs metabolism is further supported by the inhibition of polyamine degradation (# 205, 372–373, 437, 1499, Supplementary Table 4).
A considerable modulation was observed for the catabolism of proline, an amino acid whose synthesis is metabolically linked with PA synthesis through the common precursor arginine (Balestrasse et al., 2005). Experimental evidence showed that N deficiency caused a reduction in proline levels due to the stimulation of proline catabolism (Sánchez et al., 2002; del Mar Rubio-Wilhelmi et al., 2012). Thus, the upregulation of three transcripts coding for proline oxidases (ZmProOXs) under the Urea+NBPT treatment (# 92–94, Table 2) might be part of a N recycling process induced within the root cells by the inhibitor.
Phenylalanine/Tyrosine Derivatives and Phenylpropanoid Pathway
The modulation of gene expression in response to the Urea+NBPT treatment revealed a group of transcripts involved in the shikimate pathway and the downstream reactions for the synthesis of aromatic amino acids and their derivatives in roots (Table 3, Figure 7). The aromatic amino acids, phenylalanine, and tyrosine (Rösler et al., 1997) are precursors for several secondary metabolites, including phenols (coumarin, flavonoids, lignin), spermidine derivatives/conjugates and others phenylpropanoids (Fellenberg et al., 2009; Tzin and Galili, 2010).
Previous experimental evidence had shown that under abiotic stress, such as N deficiency, the activity and content of enzymes associated with phenylpropanoid synthesis, in particular phenylalanine ammonia-lyase (PAL), increase in plants (Dixon and Paiva, 1995; Kováčik et al., 2007).
Furthermore, limited N-availability induced the expression of several genes associated with phenylpropanoid metabolisms including members of the gene families encoding enzymes such as PAL, 4-coumarate CoA ligase (4CL), and cinnamate-4-hydroxylase (C4H, Wang et al., 2003; Scheible et al., 2004; Fritz et al., 2006; Yang et al., 2011).
The addition of NBPT resulted in a downregulation of the first steps of the shikimate pathway for the synthesis of arogenate (# 193–195, Table 3) while the downstream reactions for the synthesis of cinnamic acid and phenylalanine/tyrosine-derivatives were upregulated (# 227, 232–235, 237,238, 250–251, 255–256, 258, 289, 412, Table 3, Figure 7). This behavior could be part of a N recycling mechanism. While the synthesis of arogenate requires glutamate as ammonium donor for the transamination reaction, the conversion of phenylalanine into cinnamic acid by PAL results in a release of ammonium. In this way, the activation of phenylalanine derivatives biosynthesis might help ensuring the N recycling from phenylalanine and/or tyrosine (Rösler et al., 1997), avoiding the deamination of glutamate.
Lignin synthesis was negatively modulated in Urea+NBPT treated plants (caffeate O-methyltransferase, cinnamoyl-CoA reductase and cinnamyl-alcohol dehydrogenase, # 225, 236, 252-254, Table 3, Figure 7), suggesting a metabolic shift toward the synthesis of other phenolic compounds.
As compared to Urea treated plants, total N and amino acids (as Gln and Asn) contents decreased in Urea+NBPT treated plants, suggesting an imbalance between C and N compounds. This cellular response is typically described under conditions of limiting N availability, where a marked shift from the N-containing molecules to carbon-rich phenols and phenylpropanoids is observed (Fritz et al., 2006; del Mar Rubio-Wilhelmi et al., 2012).
Modulation of Genes Involved in the Acquisition of Other Nutrients
As previously reported, the urease inhibitor coordinates both Ni atoms in the urease active site (Manunza et al., 1999), thus it is plausible that in the root cells NBPT interacts with divalent cations modifying their availability for plants. Moreover, the inhibition of the urease enzyme by NBPT might also indirectly influence Ni homeostasis, with changes in the equilibrium between the active enzyme (urease binding Ni) and the apoenzyme (urease without Ni). Based on this consideration, the levels of the urease cofactor (Ni) might increase in the cytosol thus activating a response for its detoxification. Transcriptomic analyses might suggest that Urea+NBPT treated plants respond to a condition of Ni toxicity.
High concentrations of Ni are linked to a strong inhibition of the photosynthetic apparatus (chlorophyll content, thylakoid membrane, chloroplast grana structure, electron transport chain), occurring both in isolated chloroplasts and whole plant (for review see Chen et al., 2009). Besides the disruption of the photosynthetic apparatus, the toxicity symptoms due to Ni include the inhibition of CO2 assimilation, changes in the proline concentrations (Molas, 1997; Kozlov, 2005; Gajewska et al., 2006), and induced antioxidant enzyme and phenylalanine ammonia-lyase (PAL) activities (Yan et al., 2008). These responses might fit with the changes occurring after Urea+NBPT treatment, since the expression of Rubisco and other genes involved in the Calvin cycle, in proline pathway and some antioxidant and PAL enzymes were modulated by NBPT (# 58–60, 198–200, 1318–1319, 351, 466–472, 235, Supplementary Table 4, Figure 5).
In plants, metal ligands, such as nicotianamine (NA) and carboxylates (citrate and malate), act as intracellular Ni chelators (Devêvre et al., 1996; Ahonen-Jonnarth et al., 2000; Douchkov et al., 2005). The treatment with Urea+NBPT positively modulated the synthesis of the universal methyl donor S-adenosyl-L-methionine (SAM-S # 188, Supplementary Table 4), which along with putrescine (see above) participates to the production of NA (Herbik et al., 1999).
The availability of NA might allow maize plants to cope with the NBPT effect on metals homeostasis, favoring the transport or compartmentalization of divalent metals within the plant. In maize roots we observed a modulation of Yellow-Stripe-Like transporters (# 1335–1336, Table 4), which are homologous of Ni-NA transporters (Gendre et al., 2007).
When NBPT was added to the nutrient solution, the expression of some genes related to the homeostasis of Fe and Zn were also modulated in root cells (Table 4). Also this response might be related to changes in Ni homeostasis, a metal known to compete with Fe and Zn for NA complexation (log stability constants for complexes: Fe(II) 12.1–12.8, Zn(II) 14.6–15.4, Ni(II) 16.1, Benes et al., 1983; for review see Chen et al., 2009; Clemens et al., 2013). In particular, the root exposure to Urea+NBPT led to the upregulation of a gene coding for the vacuolar iron transporter (ZmVIT1, # 219; Table 4) and led to the downregulation of a putative NA transporter (ZmZIF1, # 2144; Table 4). The arabidopsis homologous (AtZIF1, Haydon et al., 2012) is a vacuolar membrane protein, which contributes to the intracellular distribution of NA. This protein would play an important role in Zn and Fe partitioning within the cells by limiting the sequestration of NA in the vacuole (Haydon et al., 2012). Therefore, in Urea+NBPT treated plants a limitation of metal sequestration in the vacuole might occur.
In conclusion, these results indicate that in presence of the urease inhibitor NBPT in the nutrient solution, the ureic-N assimilation is compromised and a general reprogramming of primary and secondary metabolic pathways occurs. In particular, the mechanisms involved in N-recycling appeared to be activated in treated plants (Figure 6). Moreover, root cells show a general response to compensate changes in the internal metal concentrations among different compartments. In particular changes in NA availability might led to unbalanced concentrations of Ni, Fe, and Zn.
LZ and AZ performed the microarray analysis; SV and RD measured the amino acid content; LZ and NT acquired and analyzed the data. LZ, NT, ZV, RP designed and oversaw the research; LZ, NT, RP wrote the article.
This work was supported by the Italian Ministry of University and Research (FIRB grant: RBFR127WJ9); the Department of Agricultural and Environmental Sciences, University of Udine (U.N.I.C.O. Project -2014).
Conflict of Interest Statement
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.
We are grateful to Prof. Enrico Martinoia, University of Zurich, for his helpful revision of the manuscript. Microarray analyses were performed at the Functional Genomics Center (FGC, Department of Biotechnology) at the University of Verona.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00845
Supplementary Figure 1. MapMan Pathway file. Image file for the pathway of urea assimilation and primary metabolism without transcript bins.
Supplementary Figure 2. MapMan Pathway file. Image file for the pathway of phenylalanine/tyrosine derivatives without transcript bins, image adapted from Tohge et al. (2013).
Supplementary Figure 3. Overview of up- (red) and down- (green) modulated transcripts in Urea+NBPT vs. Urea using MapMan-bincode classification. Numbers refer to the MapMan bincode: 1, PS; 2, major CHO metabolism; 3, minor CHO metabolism; 4, glycolysis; 5, fermentation; 6, gluconeogenesis/glyoxylate cycle; 7, OPP; 8, TCA/org transformation; 9, mitochondrial electron transport/ATP synthesis; 10, cell wall; 11, lipid metabolism; 12, N-metabolism; 13, amino acid metabolism; 14, S-assimilation; 15, metal handling; 16, secondary metabolism; 17, hormone metabolism; 18, Co-factor and vitamine metabolism; 19, tetrapyrrole synthesis; 20, stress; 21, redox; 22, polyamine metabolism; 23, nucleotide metabolism; 24, Biodegradation of Xenobiotics; 25, C1-metabolism; 26, miscellaneous; 27, RNA; 28, DNA; 29, protein; 30, signaling; 31, cell; 32, micro RNA, natural antisense; 33, development; 34, transport; 35, not assigned. Color scale refers to the fold change values of differentially expressed transcripts: red color refers to those transcripts positively regulated by Urea+NBPT treatment, while in green are transcripts negatively regulated by Urea+NBPT treatment.
Supplementary Table 1. MapMan Mapping File. Transcript mapping file for the pathway of urea assimilation and primary metabolism.
Supplementary Table 2. MapMan Mapping File. Transcript mapping file for the pathway of phenylalanine/tyrosine derivatives.
Supplementary Table 3. Real-time RT-PCR validation of a set of genes differentially expressed in microarray analysis Urea+NBPT vs. Urea. In the table are shown (starting from left column): Number of transcript according to Tables 2–4, Supplementary Table 4, transcript ID, description, gene symbol, microarray value (fold change and p. value adj.), real-time RT-PCR values (fold change and standard deviation), sequence of forward and reverse primers and primer efficiency.
Supplementary Table 4. Differentially expressed transcripts resulted by the comparison of root transcriptional profiles of Urea+NBPT plants with root transcriptional profile of Urea plants (Urea+NBPT vs. Urea comparison). In the table are shown (starting from left column): transcript ID, fold change (FC) of Urea+NBPT vs. Urea, adjusted P-value of Urea+NBPT vs. Urea, gene model ID, transcript ID, protein ID, Zmays_181 ID, Zmays_181 description, homolog protein of Arabidopsis, and rice (annotation according to Zea mays 6a annotation file, http://genome.jgi.doe.gov), MapMan bincode and name were reported for each transcript.
Abalos, D., Jeffery, S., Sanz-Cobena, A., Guardia, G., and Vallejo, A. (2014). Meta-analysis of the effect of urease and nitrification inhibitors on crop productivity and nitrogen use efficiency. Agric. Ecosyst. Environ. 189, 136–144. doi: 10.1016/j.agee.2014.03.036
Ahonen-Jonnarth, U., Van Hees, P. A. W., Lundstrom, U. S., and Finlay, R. D. (2000). Organic acids produced by mycorrhizal Pinus sylvestris exposed to elevated aluminium and heavy metal concentrations. New Phytol. 146, 557–567. doi: 10.1046/j.1469-8137.2000.00653.x
Alcázar, R., Altabella, T., Marco, F., Bortolotti, C., Reymond, M., Koncz, C., et al. (2010). Polyamines: molecules with regulatory functions in plant abiotic stress tolerance. Planta 231, 1237–1249. doi: 10.1007/s00425-010-1130-0
Alcázar, R., Cuevas, J. C., Planas, J., Zarza, X., Bortolotti, C., Carrasco, P., et al. (2011). Integration of polyamines in the cold acclimation response. Plant Sci. 180, 31–38. doi: 10.1016/j.plantsci.2010.07.022
Alcázar, R., García-Martínez, J. L., Cuevas, J. C., Tiburcio, A. F., and Altabella, T. (2005). Overexpression of ADC2 in arabidopsis induces dwarfism and late-flowering through GA deficiency. Plant J. 43, 425–436. doi: 10.1111/j.1365-313X.2005.02465.x
Artola, E., Cruchaga, S., Ariz, I., Moran, J. F., Garnica, M., Houdusse, F., et al. (2011). Effect of N-(n-butyl) thiophosphoric triamide on urea metabolism and the assimilation of ammonium by Triticum aestivum L. Plant Growth Regul. 63, 73–79. doi: 10.1007/s10725-010-9513-6
Balestrasse, K. B., Gallego, S. M., Benavides, M. P., and Tomaro, M. L. (2005). Polyamines and proline are affected by cadmium stress in nodules and roots of soybean plants. Plant Soil 270, 343–353. doi: 10.1007/s11104-004-1792-0
Clough, T. J., Kelliher, F. M., Sherlock, R. R., and Ford, C. D. (2004). Lime and soil moisture effects on nitrous oxide emissions from a urine patch. Soil Sci. Soc. Am. J. 68, 1600–1609. doi: 10.2136/sssaj2004.1600
Cruchaga, S., Artola, E., Lasa, B., Ariz, I., Irigoyen, I., Moran, J. F., et al. (2011). Short term physiological implications of NBPT application on the N metabolism of Pisum sativum and Spinacea oleracea. J. Plant Physiol. 168, 329–336. doi: 10.1016/j.jplph.2010.07.024
Cuevas, J. C., López-Cobollo, R., Alcázar, R., Zarza, X., Koncz, C., Altabella, T., et al. (2008). Putrescine is involved in arabidopsis freezing tolerance and cold acclimation by regulating abscisic acid levels in response to low temperature. Plant Physiol. 148, 1094–1105. doi: 10.1104/pp.108.122945
Davletova, S., Schlauch, K., Coutu, J., and Mittler, R. (2005). The Zinc-Finger Protein Zat12 plays a central role in reactive oxygen and abiotic stress signaling in arabidopsis. Plant Physiol. 139, 847–856. doi: 10.1104/pp.105.068254
Dawar, K., Zaman, M., Rowarth, J. S., Blennerhassett, J. D., and Turnbull, M. H. (2011). Urea hydrolysis and lateral and vertical movement in the soil: effects of urease inhibitor and irrigation. Biol. Fert. Soils 47, 139–146. doi: 10.1007/s00374-010-0515-3
del Mar Rubio-Wilhelmi, M., Sanchez-Rodriguez, E., Leyva, R., Blasco, B., Romero, L., Blumwald, E., et al. (2012). Response of carbon and nitrogen-rich metabolites to nitrogen deficiency in PSARK::IPT tobacco plants. Plant Physiol. Bioch. 57, 231–237. doi: 10.1016/j.plaphy.2012.06.004
Devêvre, O., Garbaye, J., and Botton, B. (1996). Release of complexing organic acids by rhizosphere fungi as a factor in Norway spruce yellowing in acidic soils. Mycol. Res. 100, 1367–1374. doi: 10.1016/S0953-7562(96)80065-7
Douchkov, D., Gryczka, C., Stephan, U. W., Hell, R., and Baumlein, H. (2005). Ectopic expression of nicotianamine synthase genes results in improved iron accumulation and increased nickel tolerance in transgenic tobacco. Plant Cell Environ. 28, 365–374. doi: 10.1111/j.1365-3040.2005.01273.x
Fellenberg, C., Böttcher, C., and Vogt, T. (2009). Phenylpropanoid polyamine conjugate biosynthesis in Arabidopsis thaliana flower buds. Phytochemistry 70, 1392–1400. doi: 10.1016/j.phytochem.2009.08.010
Fritz, C., Palacio-Rojas, N., Feil, R., and Stitt, M. (2006). Regulation of secondary metabolism by the carbon–nitrogen status in tobacco: nitrate inhibits large sectors of phenylpropanoid metabolism. Plant J. 46, 533–548. doi: 10.1111/j.1365-313X.2006.02715.x
Gajewska, E., Sklodowska, M., Slaba, M., and Mazur, J. (2006). Effect of nickel on antioxidative enzyme activities, proline and chlorophyll contents in wheat shoots. Biol. Plantarum 50, 653–659. doi: 10.1007/s10535-006-0102-5
Gendre, D., Czernic, P., Conéjéro, G., Pianelli, K., Briat, J. F., Lebrun, M., et al. (2007). TcYSL3, a member of the YSL gene family from the hyperaccumulator Thlaspi caerulescens, encodes a nicotianamine -Ni/ Fe transporter. Plant J. 49, 1–15. doi: 10.1111/j.1365-313X.2006.02937.x
Gentleman, R. C., Carey, V. J., Bates, D. M., Bolstad, B., Dettling, M., Dudoit, S., et al. (2004). Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 5:R80. doi: 10.1186/gb-2004-5-10-r80
Gu, R., Chen, X., Zhou, Y., and Yuan, L. (2012). Isolation and characterization of three maize aquaporin genes, ZmNIP2;1, ZmNIP2;4 and ZmTIP4;4 involved in urea transport. BMB Rep. 45, 96–101. doi: 10.5483/BMBRep.2012.45.2.96
Guo, S., Duan, J. A., Qian, D., Tang, Y., Wu, D., Su, S., et al. (2013). Rapid determination of amino acids in fruits of Ziziphus jujuba by hydrophilic interaction ultra-high-performance liquid chromatography coupled with triple-quadrupole mass spectrometry. J. Agr. Food Chem. 61, 2709–2719. doi: 10.1021/jf305497r
Hacham, Y., Avraham, T., and Amir, R. (2002). The N-terminal region of arabidopsis cystathionine γ-synthase plays an important regulatory role in methionine metabolism. Plant Physiol. 128, 454–462. doi: 10.1104/pp.010819
Haydon, M. J., Kawachi, M., Wirtz, M., Hillmer, S., Hell, R., and Krämer, U. (2012). Vacuolar nicotianamine has critical and distinct roles under iron deficiency and for zinc sequestration in arabidopsis. Plant Cell 24, 724–737. doi: 10.1105/tpc.111.095042
Hendrickson, L. L., and Douglass, E. A. (1993). Metabolism of the urease inhibitor N-(n-butyl) thiophosphoric triamide (NBPT). in soils. Soil Biol. Biochem. 25, 1613–1618. doi: 10.1016/0038-0717(93)90017-6
Herbik, A., Koch, G., Mock, H. P., Dushkov, D., Czihal, A., Thielmann, J., et al. (1999). Isolation, characterization and cDNA cloning of nicotianamine synthase from barley - A key enzyme for iron homeostasis in plants. Eur. J. Biochem. 265, 231–239. doi: 10.1046/j.1432-1327.1999.00717.x
Juan, Y. H., Chen, L. J., Wu, Z. J., and Wang, R. (2009). Kinetics of soil urease affected by urease inhibitors at contrasting moisture regimes. Rev. Cienc. Suelo Nutr. Veg. 9, 125–133. doi: 10.4067/s0718-27912009000200004
Kojima, S., Bohner, A., Gassert, B., Yuan, L., and von Wirén, N. (2007). AtDUR3 represents the major transporter for high-affinity urea transport across the plasma membrane of nitrogen-deficient arabidopsis roots. Plant J. 52, 30–40. doi: 10.1111/j.1365-313X.2007.03223.x
Kováčik, J., Klejdus, B., Bačkor, M., and Repčák, M. (2007). Phenylalanine ammonia-lyase activity and phenolic compounds accumulation in nitrogen-deficient Matricaria chamomilla leaf rosettes. Plant Sci. 172, 393–399. doi: 10.1016/j.plantsci.2006.10.001
Kozlov, M. V. (2005). Pollution resistance of mountain birch, Betula pubescens subsp. czerepanovii, near the copper-nickel smelter: natural selection or phenotypic acclimation? Chemosphere 59, 189–197. doi: 10.1152/ajpregu.00259.2015
Krogmeier, M. J., McCarty, G. W., and Bremner, J. M. (1989). Potential phytotoxicity associated with the use of soil urease inhibitors. Proc. Natl. Acad. Sci. U.S.A. 86, 1110–1112. doi: 10.1073/pnas.86.4.1110
Liu, L. H., Ludewig, U., Gassert, B., Frommer, W. B., and von Wirén, N. (2003). Urea transport by nitrogen-regulated tonoplast intrinsic proteins in arabidopsis. Plant Physiol. 133, 1220–1228. doi: 10.1104/pp.103.027409
Manunza, B., Deiana, S., Pintore, M., and Gessa, C. (1999). The binding mechanism of urea, hydroxamic acid and N-(N-butyl) phosphoric triamide to the urease active site. A comparative molecular dynamics study. Soil Biol. Biochem. 31, 789–796. doi: 10.1016/S0038-0717(98)00155-2
Mattoo, A. K., Sobolev, A. P., Neelam, A., Goyal, R. K., Handa, A. K., and Segre, A. L. (2006). Nuclear magnetic resonance spectroscopy-based metabolite profiling of transgenic tomato fruit engineered to accumulate spermidine and spermine reveals enhanced anabolic and nitrogen–carbon interactions. Plant Physiol. 142, 1759–1770. doi: 10.1104/pp.106.084400
Medina, R., and Radel, R. J. (1988). Mechanisms of urease inhibition, in Ammonia Volatilization from Urea Fertilizers. Bull Y-206. National Fertilizer Development Center, eds B. R. Bockand and D. E. Kissel (Muscle Shoals, AL: Tennessee Valley Authority), 137–174.
Mérigout, P., Lelandais, M., Bitton, F., Renou, J. P., Briand, X., Meyer, C., et al. (2008). Physiological and transcriptomic aspects of urea uptake and assimilation in arabidopsis plants. Plant Physiol. 147, 1225–1238. doi: 10.1104/pp.108.119339
Molas, J. (1997). Changes in morphological and anatomical structure of cabbage (Brassica oleracea L.) outer leaves and in ultrastructure of their chloroplasts caused by an in vitro excess of nickel. Photosynthetica 34, 513–522. doi: 10.1023/A:1006805327340
Ritz, C., and Spiess, A. N. (2008). qpcR: an R package for sigmoidal model selection in quantitative real-time polymerase chain reaction analysis. Bioinformatics 24, 1549–1551. doi: 10.1093/bioinformatics/btn227
Rosegrant, M. W., Huang, J., Sinha, A., Ahammad, H., Ringler, C., Zhu, T., et al. (2008). Exploring Alternative Futures for Agricultural Knowledge, Science and Technology (AKST). ACIAR Project Report ADP/2004/045, Washington, DC, IFPRI.
Sánchez, E., Garcia, P. C., López-Lefebre, L. R., Rivero, R. M., Ruiz, J. M., and Romero, L. (2002). Proline metabolism in response to nitrogen deficiency in french bean plants (Phaseolus vulgaris L. cv Strike). Plant Growth Regul. 36, 261–265. doi: 10.1023/A:1016583430792
Sanz-Cobena, A., Sánchez-Martín, L., García-Torres, L., and Vallejo, A. (2012). Gaseous emissions of N2O and NO and NO leaching from urea applied with urease and nitrification inhibitors to a maize (Zea mays) crop. Agric. Ecosyst. Environ. 149, 64–73. doi: 10.1016/j.agee.2011.12.016
Sato, F., and Yamada, Y. (2008). Engineering formation of medicinal compounds in cell cultures, in Advances in Plant Biochemistry and Molecular Biology Volume 1 - Bioengineering and Molecular Biology of Plant Pathways, eds H. J. Bohnert, H. Nguyen and N. G. Lewis (Heidelberg: Elsevier), 311–345. doi: 10.1042/BA20100171
Scheible, W. R., Morcuende, R., Czechowski, T., Fritz, C., Osuna, D., Palacios-Rojas, N., et al. (2004). Genome-wide reprogramming of primary and secondary metabolism, protein synthesis, cellular growth processes, and the regulatory infrastructure of arabidopsis in response to nitrogen. Plant Physiol. 136, 2483–2499. doi: 10.1104/pp.104.047019
Schnable, P. S., Ware, D., Fulton, R. S., Stein, J. C., Wei, F., Pasternak, S., et al. (2009). The B73 maize genome: complexity, diversity, and dynamics. Science 326, 1112–1115. doi: 10.1126/science.1178534
Smyth, G. K. (2005). Limma: linear models for microarray data, in Bioinformatics and Computational Biology Solutions using R and Bioconductor, eds R. Gentleman, V. Carey, S. Dudoit, R. Irizarry, and W. Huber (New York, NY: Springer), 397–420. doi: 10.3892/mmr.2015.4448
Soares, J. R., Cantarella, H., and de Campos Menegale, M. L. (2012). Ammonia volatilization losses from surface-applied urea with urease and nitrification inhibitors. Soil Biol. Biochem. 52, 82–89. doi: 10.1016/j.soilbio.2012.04.019
Thimm, O., Bläsing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., et al. (2004). MapMan: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313X.2004.02016.x
Tisi, A., Angelini, R., and Cona, A. (2011a). Does polyamine catabolism influence root development and xylem differentiation under stress conditions? Plant Signal Behav. 6, 1–4. doi: 10.4161/psb.6.11.17640
Tisi, A., Federico, R., Moreno, S., Lucretti, S., Moschou, P. N., Roubelakis-Angelakis, K. A., et al. (2011b). Perturbation of polyamine catabolism can strongly affect root development and xylem differentiation. Plant Physiol. 157, 200–215. doi: 10.1104/pp.111.173153
Trenkel, M. E. (2010). Slow-and Controlled-Release and Stabilized Fertilizers: an Option for Enhancing Nutrient Use Efficiency in Agriculture. Paris: International Fertilizer Industry Association, (IFA).
Urano, K., Yoshiba, Y., Nanjo, T., Igarashi, Y., Seki, M., Sekiguchi, F., et al. (2003). Characterization of arabidopsis genes involved in biosynthesis of polyamines in abiotic stress responses and developmental stages. Plant Cell Environ. 26, 1917–1926. doi: 10.1046/j.1365-3040.2003.01108.x
Usadel, B., Poree, F., Nagel, A., Lohse, M., Czedik-Eysenberg, A., and Stitt, M. (2009). A guide to using MapMan to visualize and compare Omics data in plants: a case study in the crop species, Maize. Plant Cell Environ. 32, 1211–1229. doi: 10.1111/j.1365-3040.2009.01978.x
Vogel, J. T., Zarka, D. G., Van Buskirk, H. A., Fowler, S. G., and Thomashow, M. F. (2005). Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of arabidopsis. Plant J. 41, 195–211. doi: 10.1111/j.1365-313X.2004.02288.x
Wang, R., Okamoto, M., Xing, X., and Crawford, N. M. (2003). Microarray analysis of the nitrate response in arabidopsis roots and shoots reveals over 1,000 rapidly responding genes and new linkages to glucose, trehalose-6-phosphate, iron, and sulfate metabolism. Plant Physiol.132, 556–567. doi: 10.1104/pp.103.021253
Wang, W. H., Liu, G. W., Cao, F. Q., Cheng, X. Y., Liu, B. W., and Liu, L. H. (2013). Inadequate root uptake may represent a major component limiting rice to use urea as sole nitrogen source for growth. Plant Soil 363, 191–200. doi: 10.1007/s11104-012-1305-5
Watson, C. J., and Miller, H. (1996). Short-term effects of urea amended with the urease inhibitor N-(n-butyl) thiophosphoric triamide on perennial ryegrass. Plant Soil 184, 33–45. doi: 10.1007/BF00029272
Wen, Y., Feng, J., and Sachs, G. (2013). Helicobacter pylori 5′ureB-sRNA, a cis-encoded antisense small RNA, negatively regulates ureAB expression by transcription termination. J. Bacteriol. 195, 444–452. doi: 10.1128/JB.01022-12
Witte, C. P., Isidore, E., Tiller, S. A., Davies, H. V., and Taylor, M. A. (2001). Functional characterisation of urease accessory protein G (ureG). from potato. Plant Mol. Biol. 45, 169–179. doi: 10.1023/A:1006429315851
Witte, C. P., Tiller, S. A., Taylor, M. A., and Davies, H. V. (2002). Leaf urea metabolism in potato. Urease activity profile and patterns of recovery and distribution of 15N after foliar urea application in wild-type and urease antisense transgenics. Plant Physiol. 128, 1129–1136. doi: 10.1104/pp.010506
Yan, R., Gao, S., Yang, W., Cao, M., Wang, S., and Chen, F. (2008). Nickel toxicity induced antioxidant enzyme and phenylalanine ammonia-lyase activities in Jatropha curcas L. cotyledons. Plant Soil Environ. 54, 294–300.
Yang, H., Menz, J., Häussermann, I., Benz, M., Fujiwara, T., and Ludewig, U. (2015). High and low affinity urea root uptake: involvement of NIP5;1. Plant Cell Physiol. 56, 1588–1597. doi: 10.1093/pcp/pcv067
Yang, X. S., Wu, J., Ziegler, T. E., Yang, X., Zayed, A., Rajani, M. S., et al. (2011). Gene expression biomarkers provide sensitive indicators of in planta nitrogen status in maize. Plant Physiol. 157, 1841–1852. doi: 10.1104/pp.111.187898
Zaman, M., Nguyen, M. L., Blennerhassett, J. D., and Quin, B. F. (2008). Reducing NH3, N2O and NO-N losses from a pasture soil with urease or nitrification inhibitors and elemental S-amended nitrogenous fertilizers. Biol. Fertil. Soils 44, 693–705. doi: 10.1007/s00374-007-0252-4
Zanin, L., Tomasi, N., Zamboni, A., Varanini, Z., and Pinton, R. (2015a). The urease inhibitor NBPT negatively affects DUR3-mediated uptake and assimilation of urea in maize roots. Front. Plant Sci. 6:1007. doi: 10.3389/fpls.2015.01007
Zanin, L., Zamboni, A., Monte, R., Tomasi, N., Varanini, Z., Cesco, S., et al. (2015b). Transcriptomic analysis highlights reciprocal interactions of urea and nitrate for nitrogen acquisition by maize roots. Plant Cell Physiol. 56, 532–548. doi: 10.1093/pcp/pcu202
Keywords: gene expression, microarray, nitrogen fertilizers, nitrogen nutrition, root uptake, urea, urea metabolism, Zea mays
Citation: Zanin L, Venuti S, Tomasi N, Zamboni A, De Brito Francisco RM, Varanini Z and Pinton R (2016) Short-Term Treatment with the Urease Inhibitor N-(n-Butyl) Thiophosphoric Triamide (NBPT) Alters Urea Assimilation and Modulates Transcriptional Profiles of Genes Involved in Primary and Secondary Metabolism in Maize Seedlings. Front. Plant Sci. 7:845. doi: 10.3389/fpls.2016.00845
Received: 10 March 2016; Accepted: 30 May 2016;
Published: 22 June 2016.
Edited by:Jan Kofod Schjørring, University of Copenhagen, Denmark
Copyright © 2016 Zanin, Venuti, Tomasi, Zamboni, De Brito Francisco, Varanini and Pinton. 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) or licensor 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.
*Correspondence: Roberto Pinton, firstname.lastname@example.org