Interactive Effects of CO2 Concentration and Water Regime on Stable Isotope Signatures, Nitrogen Assimilation and Growth in Sweet Pepper

Sweet pepper is among the most widely cultivated horticultural crops in the Mediterranean basin, being frequently grown hydroponically under cover in combination with CO2 fertilization and water conditions ranging from optimal to suboptimal. The aim of this study is to develop a simple model, based on the analysis of plant stable isotopes in their natural abundance, gas exchange traits and N concentration, to assess sweet pepper growth. Plants were grown in a growth chamber for near 6 weeks. Two [CO2] (400 and 800 μmol mol−1), three water regimes (control and mild and moderate water stress) and four genotypes were assayed. For each combination of genotype, [CO2] and water regime five plants were evaluated. Water stress applied caused significant decreases in water potential, net assimilation, stomatal conductance, intercellular to atmospheric [CO2], and significant increases in water use efficiency, leaf chlorophyll content and carbon isotope composition, while the relative water content, the osmotic potential and the content of anthocyanins did change not under stress compared to control conditions support this statement. Nevertheless, water regime affects plant growth via nitrogen assimilation, which is associated with the transpiration stream, particularly at high [CO2], while the lower N concentration caused by rising [CO2] is not associated with stomatal closure. The stable isotope composition of carbon, oxygen, and nitrogen (δ13C, δ18O, and δ15N) in plant matter are affected not only by water regime but also by rising [CO2]. Thus, δ18O increased probably as response to decreases in transpiration, while the increase in δ15N may reflect not only a lower stomatal conductance but a higher nitrogen demand in leaves or shifts in nitrogen metabolism associated with decreases in photorespiration. The way that δ13C explains differences in plant growth across water regimes within a given [CO2], seems to be mediated through its direct relationship with N accumulation in leaves. The changes in the profile and amount of amino acids caused by water stress and high [CO2] support this conclusion. However, the results do not support the use of δ18O as an indicator of the effect of water regime on plant growth.


INTRODUCTION
After tomato, sweet pepper is the second largest horticultural product cultivated in the Mediterranean basin in terms of area (del Amor, 2007). Moreover, a significant area is devoted to its cultivation under cover (greenhouses of different categories), and is frequently combined with the application of CO 2 fertilization . Essentially, a high CO 2 concentration ([CO 2 ]) stimulates photosynthesis and may contribute to increasing productivity in greenhouses while palliating environmental problems such as water stress or high temperature. Nevertheless, plant responses to elevated [CO 2 ] may be affected in one or another way as a result of photosynthetic acclimation (Long et al., 2004) or due to differences in growing conditions such as water regime (O'Leary et al., 2015). Thus, some degree of water stress may increase the efficiency of water use, and at the same time positively affect the quality of the final product (López-Marín et al., 2017). However, the interaction of elevated [CO 2 ] with water stress has not been widely studied and there are studies supporting either the positive or the negative effects of elevated [CO 2 ] on water stress tolerance (Medina et al., 2016). Exposure to elevated [CO 2 ] may mitigate the inhibition of photosynthesis under water stress and improve water use efficiency by a positive synergistic effect of both factors on stomatal closure, but can also stimulate oxidative stress and not affect plant growth (Erice et al., 2007(Erice et al., , 2014Bencze et al., 2014;Medina et al., 2016). In addition, plant size may limit the direct physiological effects of elevated [CO 2 ] (Xu et al., 2016). Moreover, some studies indicated that this interaction is highly dependent on the genotypic variability and the severity of water stress (Erice et al., 2014;Medina et al., 2016;Xu et al., 2016).
Therefore, a better understanding of the interactions between high [CO 2 ] and water stress is important for predicting the agricultural consequences of the expected increase in [CO 2 ]. In the case of sweet pepper the genotypic performance and the specific responses to the combinations of [CO 2 ] and water conditions have usually been monitored through photosynthetic and transpirative gas exchange traits (Peñuelas et al., 1995;del Amor et al., 2010). However, the use of methodologies that are able to integrate physiological processes on a larger temporal scale is an alternative that needs exploration. The analysis of the natural abundances of stable isotopes in plant matter may represent an alternative (Dawson et al., 2002;Araus et al., 2008) that is worth exploring.
Carbon isotope composition (δ 13 C), frequently expressed as discrimination ( 13 C) against the surrounding air, provides information on the effect of growing conditions on photosynthetic carbon assimilation (Farquhar et al., 1982;Condon et al., 1990). Plants discriminate against the heavier carbon isotope ( 13 C) during photosynthesis and the extent of this discrimination depends on the ratio of the intercellular vs. the atmospheric [CO 2 ] (C i /C a ) in photosynthetic organs (Farquhar et al., 1982(Farquhar et al., , 1989. Since the C i /C a ratio depends on the balance between the photosynthetic activity and the stomatal conductance (A/g s ) of the plant (Farquhar et al., 1982(Farquhar et al., , 1989Rebetzke et al., 2002;Condon et al., 2004), when analyzed in plant dry matter 13 C becomes a time-integrated indicator of the A/g s ratio and therefore of water use efficiency of the plant (Farquhar and Richards, 1984). Under drought stress, the 13 C (or δ 13 C) is also a good predictor of stomatal conductance  and of water input received by the crop (Araus et al., 2003). For growing conditions where the water regime is not the main environmental variable, it remains challenging to assess whether variation in the carbon isotope signature is the result of changes in intrinsic photosynthetic capacity or stomatal conductance (Scheidegger et al., 2000;Farquhar et al., 2007). Due to this, the analysis of the stable oxygen isotope signature in plants may contribute toward understanding the nature of the changes in δ 13 C ( . The oxygen isotope composition (δ 18 O) in plant tissues can be influenced by three main factors. The first factor is the isotopic composition of the source water taken up by the plant (Roden and Ehleringer, 1999). The second factor is the enrichment in 18 O in the leaves due to evaporation in comparison to source water (Pande et al., 1995). The third factor is the fractionation of oxygen isotopes during biochemical reactions involved in the synthesis of organic matter . Sugars and other metabolites formed in leaves incorporate the leaf water isotopic signal, which is then retained in structural organic compounds, such as cellulose (Barbour, 2007;Gessler et al., 2014). Stomatal conductance plays a crucial role in regulating the water balance of the plant. Providing there is constancy in the δ 18 O of the water taken up by the plant, the δ 18 O of plant matter integrates evaporative conditions throughout the life cycle of the plant, and this is largely unaffected by photosynthesis Farquhar et al., 2007). Therefore, δ 18 O has been proposed for estimating stomatal conductance and transpiration and thus plant water use in different species Sheshshayee et al., 2005;Farquhar et al., 2007;Cabrera-Bosquet et al., 2009a, 2011Cernusak et al., 2009a).
Nitrogen is the most growth-limiting nutrient element for plants . In fact, the nitrogen isotope composition (δ 15 N) in plant matter is an indicator of the effect of growing conditions on the nitrogen metabolism of the plant, even though a complete knowledge of the underlying biochemical mechanisms is lacking (Cernusak et al., 2009b;Tcherkez, 2010;Yousfi et al., 2012). The natural abundance of δ 15 N has been used in sweet pepper to assess the source of nitrogen fertilization (Flores et al., 2007;del Amor and Navarro, 2008). However, to the best of our knowledge, studies on the interactive effect of [CO 2 ] and water regime on the δ 15 N and δ 18 O of the plant are scarce.
This study compared the δ 13 C, δ 18 O, and δ 15 N, together with N concentration, of the leaf dry matter of sweet pepper plants growing under different [CO 2 ] and water regimes. Further the stable isotope signatures of these three elements were correlated with plant growth. Moreover, gas exchange and amino acid profiles were measured on similar leaves. The final aim of this study is to produce a single conceptual model, integrating different key physiological traits, that explains the variability in sweet pepper biomass due to growing conditions and genotypic variability. To that end, plants were grown under hydroponic conditions at a relatively low light intensity and mild to moderate water stress, resembling the growing conditions usually experienced, in a mild Mediterranean climate, by sweet pepper plants within commercial plastic/polycarbonate greenhouses adapted to CO 2 fertilization and which often include a shadow sheet (Dueck et al., 2006;del Amor and Gómez-López, 2009;del Amor et al., 2010;Pérez-Jiménez et al., 2016).

Plant Material and Growth Conditions
The experiment was carried at the Instituto Murciano de Investigación y Desarrollo Agrario y Alimentario (IMIDA), La Alberca, Murcia, Spain. Four sweet pepper (Capsicum annuum L.) cultivars were studied: Tallante (Ta; De Ruiter Vegetable Seeds, Inc.), Coyote (Co; Syngenta Seeds SA), Herminio (He; Syngenta Seeds SA), and Velez (Ve; Enza Zaden BV). These cultivars are commonly used in commercial greenhouses in SE Spain (Almeria and Murcia regions). Seedlings were transplanted to 5-l black containers filled with coconut coir fiber (Pelemix, Alhama de Murcia, Murcia, Spain) and acclimated during 5 days to the new conditions. Then three irrigation treatments were applied for 41 days: the control (an amount of 500 ml of nutrient solution was applied every day), mild stress (same amount every 2 days), and moderate stress (same amount every 3 days). Plants were irrigated with a modified Hoagland solution with the following composition in meq L −1 ; NO − 3 : 12.0; H 2 PO − 4 : 1.0; SO 2− 4 : 3.5; K + : 7.0; Ca 2+ : 4.5; Mg 2+ : 2.0. Irrigation was supplied via pressure-compensating and anti-drain drippers (2 l h −1 ) and fresh nutrient solution was applied to avoid salt accumulation, with a minimum of 35% drainage (del Amor and Gómez-López, 2009).
Plants were grown in a climate chamber designed to mimic usual environmental conditions experienced by pepper plants within greenhouses (del Amor et al., 2010), with fullycontrolled environmental parameters: 70% RH, 16/8 h day/night photoperiod, 28/16 • C and a photosynthetically-active radiation (PAR) of 250 µmol m −2 s −1 provided by a combination of fluorescent lamps (TL-D Master reflex 830 and 840, Philips, the Netherlands) and high-pressure sodium lamps (Son-T Agro, Philips, the Netherlands). Plants were grown at [CO 2 ] of 400 µmol mol −1 (atmospheric [CO 2 ]) and 800 µmol mol −1 (high [CO 2 ]). The [CO 2 ] was regulated by injection of external compressed CO 2 (bottle [CO 2 ] ≥ 99.9%), controlled by an infrared gas analyser (Dräger Politron IR CO 2 , Sweden). For each cultivar within a specific water regime and [CO 2 ], five replications (each consisting in a single plant) were run. Thus, 24 treatments were studied, corresponding to the combination of four cultivars, three irrigation levels, and two [CO 2 ], totalling 120 plants. All the study was done in the same growth chamber. Therefore, the experiment was done consecutively, with the only difference being [CO 2 ]. We carefully verified that the germination and seedling growth conditions were the same (by using a small growth chamber (BINDER KBWF 240, BINDER GmbH, Tuttlingen, Germany), with light, temperature, and RH control. Plants were grown under a randomized complete block design (n = 60).

Photosynthetic and Transpirative Gas Exchange and Chlorophyll Content
At the end of the experimental period, net CO 2 assimilation, stomatal conductance, transpiration, and the C i /C a ratio were measured in the youngest fully-expanded leaf of five plants per treatment, using a CIRAS-2 portable photosynthesis system (PP Systems, Amesbury, MA, USA) with a PLC6 (U) Automatic Universal Leaf Cuvette of 1.7 cm 2 . The cuvette provided light (LED) with a photon flux of 800 µmol m −2 s −1 , 400 or 800 µmol mol −1 [CO 2 ] and a leaf temperature of 22 • C. Water use efficiency (WUE) was determined as the ratio of net CO 2 assimilation to transpiration.
The leaf chlorophyll content on an area basis was determined in the same leaves used for gas exchange with a SPAD-502 (Konica-Minolta Sensing, Japan) portable meter. Three measurements were made on each leaf.

Leaf Water and Osmotic Potentials and Relative Water Content
The leaf water potential ( w ) was measured in the same leaves used for gas exchange determinations, using a Scholander pressure chamber (model 3000, Soil Moisture Equipment Corp., Santa Barbara, USA) as reported elsewhere (Turner, 1988). Measurements were performed after gas exchange determinations and then the leaves were put in Eppendorf tubes with holes at the bottom and rapidly frozen. These tubes were then centrifuged twice, at 4,000 g for 4 min (4 • C), using an Eppendorf centrifuge so that all sap was extracted from the samples. The osmotic potential ( π) of the leaf sap was assessed with a vapor pressure osmometer (Wescor 5500, Logan, Utah, USA) used to measure the osmolality (mmol kg −1 ) of the expressed sap; this was converted to osmotic potential according to the Van't Hoff equation: π (MPa) = −RTC, where R is the gas constant (0.00832 l MPa K −1 mol −1 ), T is the temperature (293 K) and C is the number of moles of solute in 1 kg of water (= 1 l at 293 K). The leaf relative water content (RWC) was measured on same-age leaves as those used for π. Three small disks (2.07 cm 2 ) per leaf from each of the six plants were cut and weighed immediately to obtain fresh mass (FM), and then they were placed for 24 h in the dark in a beaker (30 cm 3 ) filled with distilled water. After this, they were reweighed to obtain turgid fresh mass (TM), and dry mass (DM) after drying at 80 • C for 48 h. The relative water content (RWC), expressed as a percentage, was calculated as RWC = [(FM-DM)/(TM-DM)] ×100%.

Shoot Biomass
Plants were harvested at the end of the experiment (41 days after transplanting), and 120 plants (five plants per treatment) were analyzed. The aerial parts (thereafter referred as shoot biomass, including leaves plus stems and petioles) were dried and the dry weight determined after a minimum of 72 h at 65 • C. The specific leaf area (SLA) was calculated as the ratio between the area and the dry weight of leaf discs of 6.91 cm 2 .
Leaf N-Total Concentration and Stable Carbon, Nitrogen, and Oxygen Isotope Composition Total nitrogen concentration and the stable carbon ( 13 C/ 12 C) and ( 15 N/ 14 N) isotope ratios in the whole pool of shoot leaves were measured using an elemental analyser (Flash 1112 EA, ThermoFinnigan, Germany) coupled with an isotope ratio mass spectrometer (Delta C IRMS, ThermoFinnigan, Germany), operating in continuous mode. Samples of 1 mg and reference materials were weighed into tin capsules, sealed, and then loaded into an automatic sampler (ThermoFinnigan, Germany) prior to EA-IRMS analysis. Measurements were carried out at the CCiT (Centres Científics i Tecnològics) of the University of Barcelona. Nitrogen was expressed as a concentration (percent of dry weight). The 13 C/ 12 C ratios were expressed in δ notation (Coplen, 2008): δ 13 C (‰) = ( 13 C/ 12 C) sample /( 13 C/ 12 C) standard −1, where "sample" refers to plant material and "standard" to international secondary standards of known 13 C/ 12 C ratios (IAEA CH7 polyethylene foil, IAEA CH6 sucrose, and USGS 40 L-glutamic acid) calibrated against Vienna Pee Dee Belemnite calcium carbonate (VPDB) with an analytical precision (SD) of 0.10‰. The same δ notation was used for the 15 N/ 14 N ratio (δ 15 N), but in this case using international secondary standards of known 15 N/ 14 N ratios (IAEA N 1 and IAEA N 2 ammonium sulfate and IAEA NO 3 potassium nitrate) referred to N 2 in air, with an analytical precision of 0.18‰.
For the δ 18 O, the 18 O/ 16 O was determined by an online pyrolysis technique using a Thermo-Chemical Elemental Analyser (TC/EA Thermo Quest Finnigan, Germany) coupled with an IRMS (Delta C Finnigan MAT, Germany). Samples of 1 mg were weighed into silver capsules, sealed and ovendried at 60 • C for not less than 72 h to remove moisture and loaded into an automatic sampler. Results were expressed as δ 18 O values, using two secondary standards (IAEA 601 and IAEA 602) calibrated against Vienna Standard Mean Oceanic Water (VSMOW) (Coplen, 2011); the analytical precision was ≈ 0.25%. Analyses were conducted at Iso-Analytical Limited (Crewe, Cheshire, UK).

Free Amino Acids
The free amino acids were extracted from leaves (frozen at −80 • C): the sap was extracted, after vortexing at 5,000 rpm (10 min, 4 • C), and analyzed by the AccQ·Tag-ultra ultraperformance liquid chromatography (UPLC) method (Waters, UPLC Amino Acid Analysis Solution, 2006). For derivatization, 70 µl of borate buffer was added to 10 µl of the fruit sap and 20 µl of reagent solution. The reaction mixture was mixed instantly and heated at 55 • C for 10 min. After the temperature was lowered, an aliquot of the reaction mixture was used for injection. The column was an Acquity BEH C18 1.7 µm, 2.1-100 mm (Waters), and the wavelengths were set at 266 nm (excitation) and 473 nm (emission). The solvent system consisted of two eluents: (A) AccQ·Tag-ultra eluent A concentrate (5%, v/v) and water (95%, v/v); and (B) AccQ·Tag-ultra eluent B. The following elution gradient procedure was used for the analysis: 0-0.54 min, 99.9% A−0.1% B; 5.74 min, 90.9% A− 9.1% B; 7.74 min, 78.8% A−21.2% B; 8.04 min, 40.4% A−59.6% B; 8.05-8.64 min, 10% A−90% B; 8.73-10 min, 99.9% A−0.1% B. The injection volume was 1 µl, with a flow rate of 0.7 ml min −1 . The temperature of the column was maintained at 55 • C. External standards (Thermo Scientific) were used for the quantification of the amino acids, and Empower 2 (Waters) software for data acquisition and processing.

Statistical Analysis
Data were subjected to factorial ANOVA to test the effects of the growing conditions ([CO 2 ] and water regime), genotype, and their interactions. Mean comparisons were performed using Tukey's honestly significant difference (HSD) test. A bivariate correlation procedure was used to calculate the Pearson correlation coefficients between the different traits measured. Multiple linear regression (stepwise) analysis was used to analyse the criterion included to explain variation in shoot biomass under different growing conditions. Principal component analysis was produced to analyse the interrelationships among the shoot biomass, leaf nitrogen concentration, and chlorophyll content, the stable isotope composition of C, O, and N, and the gas exchange variables. Data were analyzed using IBM SPSS Statistics 24 (SPSS Inc., Chicago, IL, USA). Figures were created using a Sigma-Plot 11.0 program for Windows (Systat Software Inc., Point Richmond, CA, USA). Finally, we performed path analyses (Li, 1975) to quantify the relative contributions of direct and indirect effects of stable isotopes and other key traits on aboveground biomass. This methodology offers the possibility of building associations between variables that is based on prior knowledge. A path analysis determines simple correlations between independent factors (in this case δ 13 C), and regresses them on each intermediary (C i /C a , g s , δ 15 N and N concentration) or dependent factor (shoot biomass) to obtain direct effects in the form of partial regression coefficients (i.e., path coefficients). This model was aimed at understanding biomass responses to genotypic differences across water regimes under different levels of [CO 2 ]. A model with a comparative fit index (CFI) (Arbuckle, 1997) with values > 0.9 was taken as indicative of a good fit. Data were analyzed using the Amos Graphics package (IBM SPSS Amos, USA). A clustered heat map of amino acid profile was built using the R statistics environment (R Development Core Team, Vienna, Austria). Additionally, after generating a correlation matrix of all parameters analyzed in R, we performed a network analysis of significant correlations (Pearson correlation coefficient cut-off of 0.7 and P < 0.05) under ambient and elevated [CO 2 ] using Cytoscape software (Shannon et al., 2003).

Effect of Growing Conditions and Genotype on Shoot Biomass
Compared with control plants water stress negatively affected shoot biomass (SB), plant height (PH), and the leaf water potential ( w ), whereas leaf chlorophyll (LC) content slightly increased and no differences existed in the specific leaf area (SLA), relative water content (RWC), leaf osmotic potential ( π), and anthocyanin content ( Table 1). Increasing ambient [CO 2 ] significantly increased SB, PH, LC, anthocyanin, and π, whereas all the other traits reported in Table 1 decreased. A genotypic effect was significant for all the traits except RWC, anthocyanin, and π.
The interactions between water regime (WR) and [CO 2 ] were significant for all the traits of Table 1 except RWC, anthocyanin, and w , meaning that, except for these two traits, the response to water regime differed depending on [CO 2 ]. Most of the interactions of genotypes with growing conditions were not significant, even when it is worth to mention the significant interaction between genotype and [CO 2 ] for SB.

Effect of Growing Conditions and Genotype on Gas Exchange and Stable Isotope Signatures
Water stress significantly decreased the leaf net CO 2 assimilation (A), stomatal conductance (g s ), the intercellular to ambient CO 2 concentration (C i /C a ), the stable oxygen isotope composition (δ 18 O) and the nitrogen concentration (N), whereas the water use efficiency (WUE), the stable carbon isotope composition (δ 13 C), and the stable nitrogen isotope composition (δ 15 N) increased and no significant differences existed for the transpiration (T) ( Table 2). Rising [CO 2 ] significantly increased A, g s , T, C i /C a , WUE, and δ 18 O, whereas N, δ 13 C and δ 15 N decreased. The genotypic effect was significant for all traits except A and T. The interaction between WR and [CO 2 ] was significant for all the traits except δ 18 O (P = 0.060) and T. The interactions between genotype and [CO 2 ] were significant for A, g s , C i /C a, and WUE, and near significant for T (P =0.056), whereas the interaction between genotype and WR were also significant for g s , C i /C a , and WUE. Except for an interaction between genotype and WR for δ 15 N, no interactions between genotype and growing conditions were found.
Given the significance of the interactions between WR and [CO 2 ] for most of the traits in Tables 1, 2, a subset of traits in these tables was analyzed across water regimes within each [CO 2 ] SB, shoot biomass (g DW); PH, plant height (cm); SLA, specific leaf area (cm 2 /g); LC, leaf chlorophyll content (SPAD units); anthocyanin content (absorbance g −1 dry matter); RWC, leaf relative water content (%); w , leaf water potential (bar); π, leaf osmotic potential (bar).Values within each water regime are the means ± SE of 40 measurements (two levels of CO 2 concentration, four genotypes, and five replicates per genotype), CO 2 concentration values are the means ± SE of 60 measurements (three water regimes, four genotypes, and five replicates per genotype), while genotypic values are the means of 30 ± SE measurements (two CO 2 concentrations, three water regimes, and five replicates per genotype). For each replicate a single plant was used. Means followed by different letters are significantly different (P < 0.05) according to Tukey's honestly significant difference (HSD) test. Analysis of variance for the same variables is shown for the water regime (WR), CO 2 concentration (CO 2 ), genotype (G), and interaction (WR x CO 2 ), (WR x G), (CO 2 x G), (WR x CO 2 x G) effects. The associated percentage of the sum of squares and probabilities (ns, not significant; * P < 0.05; ** P < 0.01 and *** P < 0.001) are shown.
Frontiers in Plant Science | www.frontiersin.org TABLE 2 | Water regime, CO 2 concentration, and genotype effects on gas-exchange parameters, carbon, oxygen, and nitrogen isotope composition, and nitrogen concentration of four sweet pepper genotypes grown under different combinations of CO 2 concentration and water supply.
WATER REGIME A, leaf net CO 2 assimilation (µmol CO 2 m −2 s −1 ); gs, stomatal conductance (mmol CO 2 m −2 s −1 ); C i /Ca, intercellular to ambient CO 2 concentration; T, transpiration rate (mmol H 2 O m −2 s −1 ); WUE, water use efficiency (µmol CO 2 mmol H 2 O −1 ); N, nitrogen concentration (% DW); δ 13 C, stable carbon isotope composition (‰), δ 18 O stable oxygen isotope composition (‰); δ 15 N, stable nitrogen isotope composition (‰). Values within each water regime are the means ± SE of 40 measurements (two levels of CO 2 concentration, four genotypes, and five replicates per genotype), CO 2 concentration values are the means ± SE of 60 measurements (three water regimes, four genotypes, and five replicates per genotype), while genotypic values are the means ± SE of 30 measurements (two CO 2 concentrations, three water regimes, and five replicates per genotype). For each replicate a single plant was used. Means followed by different letters are significantly different (P < 0.05) according to Tukey's honestly significant difference (HSD) test. Analysis of variance for the same variables is shown for the water regime (WR), CO 2 concentration (CO 2 ), genotype (G), and interaction (WR x CO 2 ), (WR x G), (CO 2 x G), (WR x CO 2 x G) effects. The associated percentage of the sum of squares and probabilities (ns, not significant; * P < 0.05; ** P < 0.01 and *** P < 0.001) are shown. (

Relationships of Shoot Biomass with Gas Exchange and Stable Isotopes
The range of A rates measured across the different water conditions at the end of the experiment was only weakly correlated with SB (r = 0.34, P < 0.05) at atmospheric [CO 2 ] and was not correlated at high [CO 2 ] (data not shown). In fact, while A and g s were less affected by water limitation at high compared to atmospheric [CO 2 ], the opposite occurred for leaf growth, where the greatest decrease in biomass occurred at high [CO 2 ].
The single correlations between the signatures of the different isotopes against SB were plotted for each [CO 2 ] level across the three water regimes. δ 13 C correlated negatively against SB at high [CO 2 ], whereas the negative relationship at atmospheric [CO 2 ] did not reach significance ( Figure 1A). By contrast δ 18 O correlated with SB in a weak, albeit significant, positive manner at atmospheric [CO 2 ], whereas the positive relationship did not reach significance at high [CO 2 ] (Figure 1B). δ 15 N correlated negatively with SB at high [CO 2 ],whereas no correlation existed at atmospheric [CO 2 ] (Figure 1C). The nitrogen concentration correlated positively with SB at both [CO 2 ] levels ( Figure 2C) in a stronger manner than any of the three stable isotopes. In addition, the N concentration correlated negatively with δ 13 C at atmospheric [CO 2 ] and in a far stronger way at high [CO 2 ] (Figure 2A). Moreover, δ 13 C correlated with the total shoot nitrogen content (calculated as a SB x N/100) in a weaker manner (r = −0.43 and −0.64, both P < 0.01, for 400 and 800 µmol mol −1 [CO 2 ], respectively) than the nitrogen concentration alone (Figure 2A). In contrast the N concentration correlated negatively with δ 15 N only at high [CO 2 ] ( Figure 2B) and a weak positive correlation existed between nitrogen concentration and δ 18 O at atmospheric [CO 2 ] (r = 0.33, P < 0.05; data not shown). Moreover, g s correlated positively with the N concentration at Means followed by different letters are significantly different (P < 0.05) according to Tukey's honestly significant difference (HSD) test. Analysis of variance for the same variables is shown for the water regime (WR), genotype (G), and interaction (WR x G) effects. The associated percentage of the sum of squares and probabilities (ns, not significant; * P < 0.05; ** P < 0.01 and *** P < 0.001) are shown. Values within each water regime are the means ± SE of 20 measurements (four genotypes and five replicates per genotype). For each replicate a single plant was used.SB, shoot biomass (g DW); A, leaf net CO 2 assimilation (µmol CO 2 m −2 s −1 ); gs, stomatal conductance (mmol CO 2 m −2 s −1 ); C i /Ca, intercellular to ambient CO 2 concentration; N, nitrogen concentration (%, DW); δ 13 C, stable carbon isotope composition (‰), δ 18 O stable oxygen isotope composition (‰); δ 15 N, stable nitrogen isotope composition.
both [CO 2 ] levels, and while it also correlated with SB, this was only at atmospheric [CO 2 ] ( Figure S1). C i /C a correlated with the N concentration and SB in a similar way, but in a somewhat weaker manner than g s . In order to get an overall view of the relationships between shoot biomass and all the different physiological traits, a principal component analysis (PCA) was undertaken that included SB and LC together with the gas exchange traits and stable isotope signatures in Table 2. For all the water regimes and [CO 2 ] combined, the two first components explained 60% of variability. SB was placed nearly opposite to δ 13 C and δ 15 N and surrounded, at a certain distance, by δ 18 O, N and the A, g s , C i /C a , and T gas-exchange parameters, whereas LC was placed further away ( Figure 3A). For the three water regimes at atmospheric [CO 2 ] the two first component axes of the PCA accounted for more than 62% variability. In this case SB, and particularly N, were placed clearly opposite to δ 13 C, LC and δ 15 N, with the first isotope the furthest away and the second the closest to the center (Figure 3B), whereas C i /C a (together with N) surrounded SB. The rest of the gas exchange traits were placed more (g s , and T) or less (A) close to C i /C a , while δ 18 O had the least alignment to SB. In the case of the PCA for the three water regimes at high [CO 2 ], the two first component axes explained around 55% of variability. As in the previous PCAs, δ 13 C and δ 15 N were placed opposite to SB, whereas for the rest of the parameters only N was placed relatively near SB and all the gas exchange traits, together with δ 18 O and LC, were placed on the same side of the representation as SB but far away from it ( Figure 3C).
The relationships between SB with the different traits of Table 2 in combination were assessed through a stepwise regressions analysis (Table 4). At atmospheric [CO 2 ] and the three water regimes combined, the first variable chosen by the model was g s , which on its own explained 30% of the variability in SB. The second variable chosen by the model was N, with the two variables explaining together 44% of the variability in SB. At elevated [CO 2 ] and the three water regimes combined, the first variable chosen by the model was N concentration, which on its own explained 51% of the variability in SB. The second variable chosen by the model was Ci with the two variables explaining 58% of the variability in SB. In control conditions and both [CO 2 ] combined, the first variable chosen by the model was δ 13 C, explaining 34% of the variability in SB. The second variable chosen by the model was g s ; the two variables together explaining 44% of the variability in SB. Under mild stress and both [CO 2 ] combined, the first variable chosen by the model was also δ 13 C, which on its own explained 76% of the variability; the second variable chosen by the model was N concentration, with the two variables together explaining 81% of the variability in SB. Concerning the moderate stress, only T was chosen by the model and it explained merely 16% of the variability in SB.

Path Analysis
Further, we used the traits best correlated with shoot biomass to develop a conceptual model via a path-analysis. Besides SB, the model included δ 13 C and δ 15 N, (as time-integrated indicators of water conditions and nitrogen metabolism, respectively), together with the N concentration and the gas exchange traits, g s and C i /C a (Figure 4). In this case only the models across water regimes within each of the two [CO 2 ] levels were assayed, but excluding the model combining the two [CO 2 ], because the different δ 13 C of the ambient air and the compressed CO 2 . At atmospheric [CO 2 ], both the g s and N concentration had quite    Levels of significance: *P < 0.05 and ***P < 0.001. Abbreviations for variables are as defined in Table 2.
similar direct positive effects on the shoot biomass. However, g s also indirectly affected N concentration through its strong direct effect on C i /C a and δ 13 C. Thus, δ 13 C was negatively related with the N concentration in dry matter. The relationship of δ 13 C on δ 15 N was small, whereas there was no direct relation of δ 15 N on the N concentration. At high [CO 2 ] most of the direct effect on shoot biomass corresponded to the N concentration, whereas the effect of g s was minor and apparently negative. The direct effect of g s on C i /C a was very strong but its direct effect on δ 13 C was minor. However, the effect of δ 13 C on N concentration was stronger than atmospheric [CO 2 ] and included a direct negative effect, together with an indirect effect mediated through changes in δ 15 N.

Amino Acid Profile and Network Analysis
The profile of 17 amino acids were assayed in the leaves of the four sweet pepper genotypes grown under contrasting [CO 2 ] and water regimes and plotted in a hierarchically clustered heat map (Figure 5). The significance of the three main factors and their interactions revealed that water regime, [CO 2 ] and the CO 2 x water stress interaction were the most relevant conditions affecting amino acid contents (all of them except lysine, methionine, and histidine). Although genotypic variability was only significant for three amino acids (proline, glycine, and glutamate), the interaction genotype by [CO 2 ] affected nine amino acids, almost exactly the amino acids that were altered under the CO 2 x water regime. CO 2 enrichment significantly decreased the content of 8 amino acids (serine, asparagine, glutamate, threonine, proline, cysteine, and valine), and tended to decrease the content of other amino acids, and only increased the levels of tyrosine. Water stress, regardless of the severity, decreased the contents of eight amino acid (serine, glycine, asparagine, threonine, alanine, proline, cysteine, and valine). In spite of some genotypic differences the [CO 2 ] x water regime interaction highlighted that, although elevated CO 2 and water stress led to a decrease of amino acid levels, this decrease was less clear under mild compared with the most severe water stress. Based on significant correlations between trait pairs we built a correlation network for ambient and elevated [CO 2 ] treatments, including amino acid content ( Figure S2). They showed similar number of nodes (32 and 31, respectively) and edges (128), although there were some differences with more negative correlations under elevated [CO 2 ]. In both networks amino acid levels were positively correlated between them. However, most amino acids correlated positively with shoot biomass under elevated [CO 2 ], while at ambient [CO 2 ] only cysteine and tyrosine did. In ambient [CO 2 ] network ( Figure  S2A) amino acids levels were also positively associated with SLA, N content, and δ 18 O and negatively with δ 15 N. By contrast, in elevated [CO 2 ] network ( Figure S2B) amino acids levels were associated with more traits: positively with PH, A, and N content and negatively with RWC, w , δ 13 C, T, and δ 15 N.

DISCUSSION
Water stress affected negatively plant growth as compared with control conditions and regardless the [CO 2 ] assayed. Water stress also decreased w but did not have any effect on RWC and π. These kinds of responses do not support osmotic adjustment, which contrasts with the available literature (e.g., Wullschleger and Oosterhuis, 1991;del Amor et al., 2010), and may be due to the mild to moderate nature of the water stress applied. In fact, the relatively small (but significant) changes in A rates, g s , LC, and δ 13 C and the lack of differences for SLA and anthocyanins under stress compared to control conditions support this statement. Moreover, although anthocyanins have protective functions during drought stress, their contribution to osmotic regulation might be low (Manetas, 2006).
The effects of high [CO 2 ] on increasing plant growth and biomass have been widely reported in many plant species including pepper (del Amor et al., 2010), with the positive effect being less evident at the most severe water stress (Peñuelas et al., 1995;Medina et al., 2016). The decrease in SLA has also been reported following increases in [CO 2 ] (Peñuelas et al., 1995;Piñero et al., 2016) and as a consequence of water stress (Xu et al., 2014). These results suggest that the leaf thickness of mesophyll packing increased in response to high [CO 2 ] (Oberbauer et al., 1985). The LC and anthocyanins also increased in response to high [CO 2 ]. Previous results in pepper only exhibited a trend toward higher chlorophyll content following exposure to high [CO 2 ] (Peñuelas et al., 1995). The increase of anthocyanins under elevated [CO 2 ], as it was also observed in Takatani et al. (2014) in response to higher carbon/nitrogen balance observed under these growth conditions. Indeed, the accumulation of anthocyanins is an indicative of nitrogen limitation in the plant according to these authors. The high sucrose levels usually reported under elevated [CO 2 ] might stimulate the expression of MYB75/PAP1 transcription factor that further enhances the FIGURE 4 | Path analyses of four sweet pepper varieties grown under different water regimes and two [CO 2 ]. The conceptual model quantifying the relative strengths of direct and indirect relationships of gas exchange, stable isotope compositions and nitrogen concentration on shoot biomass (SB) is shown in (A). Physiological parameters included in the model are: g s , stomatal conductance; C i /C a , the ratio of the intercellular vs. the atmospheric [CO 2 ]; δ 13 C, carbon isotope composition; δ 15 N, nitrogen isotope composition and N concentration of leaves. The width of arrows is proportional to the path coefficient values. Dashed lines indicate negative relationships. Overall fit statistics for each path model (x 2 and comparative fit index, CFI), the latter useful for small sample sizes (with values >0.9 taken as indicative of a good fit), are shown at the bottom right of each panel. **P < 0.01, ***P < 0.001. production of anthocyanins (Tzin and Galili, 2010). Moreover, the levels of tyrosine, a precursor of anthocyanins, were induced under elevated [CO 2 ] and they could also induce the anthocyanin biosynthesis, as showed in Arabidopsis thaliana . Regarding the genotype x water regime interaction, we observed that in most of the genotypes anthocyanin content increased with water stress as a response to water limitation. However, in the case of Coyote we observed a reduction of anthocyanins with water stress; that could suggest a limitation of the protective functions of anthocyanins in response to water stress in this genotype, which in part accounts for the lowest shoot biomass of this genotype.
In addition, high [CO 2 ] decreased leaf RWC and w , and increased π. A previous study in pepper found a tendency to lower w in response to high [CO 2 ] under well-watered conditions, whereas the RWC increased regardless of the water regime considered, and the π increased but only under water stress (del Amor et al., 2010). However, the high [CO 2 ] concentration assayed in this paper was nearly three times higher than in the present study. The effect of high [CO 2 ] decreasing water status seems fairly surprising, but it may be due to a larger transpiration area compared with plants grown under normal [CO 2 ].
As expected high [CO 2 ] increased A and WUE compared with atmospheric [CO 2 ], but surprisingly did not decrease g s and transpiration under control and mild-stress conditions, and even these rates increased under moderate water stress. A small number of studies have noted similar patterns, including work on sweet pepper (del Amor et al., 2010;del Amor, 2013), other herbaceous crops (Morison, 1998;Medina et al., 2016) and in tree species that were not acclimated (Medlyn et al., 2001). The results of the current study suggest that there was no acclimation of g s to elevated [CO 2 ]. Moreover, the response of photosynthesis to growth in elevated [CO 2 ] is commonly tested by comparing the gas exchange of plants grown at atmospheric-and high-[CO 2 ] at the same [CO 2 ] across both sets of plants (e.g., Drake et al., 1997;Curtis and Wang, 1998).
Effect of [CO 2 ] and Water Regime on δ 13 C, δ 18 O, and δ 15 N The increase in δ 13 C under water stress is in line with earlier work (Farquhar et al., 1982(Farquhar et al., , 1989. However, the effect was more evident at high [CO 2 ], which agrees with the findings of del Amor (2013) for sweet pepper plants that were subjected to different levels of salinity under these two different [CO 2 ]. This may be due to the δ 13 C dilution effect of the CO 2 used to raise the air [CO 2 ] up to 800 ppm (i.e., the industrial CO 2 is a by-product of combusting fossil fuels and therefore its δ 13 C is far more negative than atmospheric CO 2 ). Indeed, δ 13 C decreased around 3‰ for each 100 µmol mol −1 of increase in [CO 2 ], which is in line with previous reports (Schubert and Jahren, 2012;del Amor, 2013). Therefore at high [CO 2 ], water stress slows the increase of new biomass with lower δ 13 C. The effect (pointed out above) of high [CO 2 ] decreasing water status due to a larger transpiration area, compared with plants grown under normal [CO 2 ], may be also involved in the larger range of δ 13 C values as response to water stress. Thus, at high [CO 2 ] the absolute changes in δ 13 C were more in parallel with the SB than with the net assimilation or C i /C a . In that regard, δ 13 C correlated negatively with SB across water regimes at high [CO 2 ], but the negative relationship did not reach significance at atmospheric [CO 2 ]. Significantly, stronger negative relationships between δ 13 C and SB at high vs. atmospheric [CO 2 ] have been reported before for sweet pepper when plants growing across different levels of salinity were compared (del Amor, 2013). Moreover, genotypic effects on δ 13 C were significant at both [CO 2 ]. To the best of our knowledge there are no studies reporting on the genotypic variability in δ 13 C of sweet pepper under different levels of water stress and [CO 2 ].
Water stress slightly decreased δ 18 O but only at high [CO 2 ]. Such decreases in δ 18 O in response to water stress and a lower transpiration are opposite to most of the reports Farquhar et al., 2007;Cabrera-Bosquet et al., 2009a, 2011. The increase in δ 18 O in plant material exposed to high [CO 2 ] has been reported before and may be caused by increases in the δ 18 O of the leaf water (Cooper and Norby, 1994). Although it has been shown that the oxygen isotopic content of atmospheric carbon dioxide has little direct influence on the δ 18 O of either leaf water or cellulose (DeNiro and Epstein, 1979), increasing atmospheric carbon dioxide concentrations may have significant indirect effects on heavy stable isotope enrichment in leaf water. This would result from the expected decreases in transpiration and increases in WUE by plants as [CO 2 ] levels increase (Eamus and Jarvis, 1989). The processes leading to concentration of the heavy isotopes 18 O in leaf water are similar in many respects to evaporation, which alters the isotopic composition of terrestrial surface waters (Cooper and Norby, 1994). Alternative explanations for the increase in δ 18 O at high [CO 2 ] are not supported by our results. Indeed, a higher CO 2 release due to photorespiration at atmospheric compared with high [CO 2 ] levels would cause an increase in δ 18 O at atmospheric [CO 2 ] levels , which does not agree with the increase we found in δ 18 O at high [CO 2 ]. Changes in g s do not seem to be involved because, regardless of the presence or absence of differences in g s within a given water regime between the atmospheric and high [CO 2 ], all water regimes at high [CO 2 ] exhibited higher δ 18 O than the corresponding water regimes at atmospheric [CO 2 ]. Moreover, genotypic variability was only significant at atmospheric [CO 2 ]. In spite of some positive results Cabrera-Bosquet et al., 2009b), the weak performance of δ 18 O in correlating with biomass and yield has been extensively reported Bort et al., 2014;Foulkes et al., 2016;Munjonji et al., 2016). This poor performance appears to be caused by post photosynthetic fractionations of the 18 O signature in the photoassimilates (Sánchez-Bragado et al., 2016). The same reasoning may be extended for the weak and erratic correlations of δ 18 O with gas exchange traits and δ 13 C that we found.
Plants under the strongest water stress exhibited slightly lower δ 15 N values than the other two treatments at atmospheric [CO 2 ], and the opposite occurred at high [CO 2 ]. Decreases in δ 15 N following water limitation Bort et al., 2014) or growing conditions causing water stress, such as salinity (Yousfi et al., 2009(Yousfi et al., , 2012, have been reported before. The effect of rising [CO 2 ] on δ 15 N was clearer than the effect of the water regime; in this case decreasing the isotopic composition, irrespective of the water regime considered. Depletion of the heavier N isotope in plants grown under high [CO 2 ] and water deficit conditions has been reported before in a study with alfalfa (Ariz et al., 2015), whereas another recent study, this time in durum wheat, concluded that elevated [CO 2 ] was the main factor that increased δ 15 N (Medina et al., 2016). The decrease in δ 15 N under elevated [CO 2 ] may reflect decreased g s , but could also be related to a higher nitrogen demand in leaves, as suggested by the decreased in leaf N (Ariz et al., 2015) or shifts in nitrogen metabolism associated with decreases in photorespiration (Tcherkez, 2010).
The photorespiratory nitrogen cycle extending over different plant compartments implicates several reactions related to nitrogen recycling that may have 15 N-kinetic isotope effects (Yu and Woo, 1991;Werner and Schmidt, 2002). Normally, the kinetic isotope effects on these reactions should not become evident when there is a total recycling of the intermediates without input or net production of substrates and products (Werner and Schmidt, 2002). However, plants can excrete gaseous ammonia (Francis et al., 1997;Pearson et al., 1998) as a consequence of photorespiration. If the uptake/loss of NH 3 from plant stomata is rate-limited by the diffusion of NH 3 in air, the transported NH 3 will be depleted in 15 N by around 18‰ relative to the δ 15 N of the source (Farquhar et al., 1980;Tcherkez and Hodges, 2008). As a consequence an increase in δ 15 N may be expected as a response to photorespiration, which may be the case for plants exposed to atmospheric [CO 2 ] compared to enhance [CO 2 ].
The correlation of δ 15 N with SB across water regimes was only significant (and negative) at the high [CO 2 ]. The lack of correlation at atmospheric [CO 2 ] may be due, as for δ 13 C, to the relatively narrow range of variability in SB associated with the response to water regimes at atmospheric [CO 2 ]. Negative relationships between δ 15 N and SB and yield at atmospheric [CO 2 ] have been frequently reported (Robinson et al., 2000;Yousfi et al., 2009;Raimanová and Haberle, 2010;Araus et al., 2013).

Growing Conditions, N Concentration, and Plant Growth
The effects of water stress decreasing the N concentration in leaves have been extensively reported (Shangguan et al., 2000;Yousfi et al., 2012). In accordance with this, the N concentration was positively correlated with the g s (r = 0.36, P < 0.01; r = 0.51, P < 0.001, at atmospheric and high [CO 2 ], respectively), suggesting that the N concentration in leaves depends to some extent on the transpirative stream (Dalla Costa and Gianquinto, 2002). This may be the case for plants growing under hydroponic conditions in particular, where water and nitrogen are provided together through the nutrient solution (Peñuelas et al., 1995) and the transpiration stream largely determines the availability of mineral N in the rhizosphere (Gonzalez-Dugo et al., 2010). In support of that, our results show a general decrease in the content of the different amino acids as response to water stress. Such decrease is nonspecific, as shown by the fact these amino acids belong to different metabolic pathways (Galili et al., 2016). This may highlight that primary substrate(s) for the synthesis of all amino acids were reduced, involving the provision of C skeletons or N. It has been reported for Arabidopsis plants that a reduction in the levels of transpiration, decreased the capacity for nitrogen uptake and the shoot nitrogen concentration of the plant but only when water availability was not restricted (Hepworth et al., 2015).
Increasing the [CO 2 ] had a significant, albeit minor effect, decreasing the N concentration. In fact, the results of an increase in [CO 2 ] decreasing N concentration were only significant at the strongest water stress. A decrease in nitrogen concentration has been widely reported in non-leguminous plants following increases in [CO 2 ] (Jablonski et al., 2002), irrespective of the water regime (Medina et al., 2016). In fact, the predicted growth response to elevated [CO 2 ] is reduced at low N availability (McMurtrie et al., 2008;Vicente et al., 2016b). In our study, water deficit combined with high [CO 2 ] caused the lowest SB and N concentration in the leaves among the six different growing conditions.
As expected at elevated [CO 2 ], the levels of photorespiratory intermediates, glycine and serine decreased (Geiger et al., 1998(Geiger et al., , 1999Yu et al., 2012;Aranjuelo et al., 2013;Noguchi et al., 2015). However, different to previous studies (Fritz et al., 2006;Krapp et al., 2011;Noguchi et al., 2015), the levels of (other) major amino acids (e.g., aspartate, asparagine, glutamate, and alanine) decreased. Such decrease has been reported as response to insufficient nitrogen conditions; specifically under NO − 3 and ambient [CO 2 ] conditions (Noguchi et al., 2015;Vicente et al., 2016b). This pattern suggests that in our system some limitation of N availability may be present. Nevertheless, the relative decrease of glycine as response to inhibition of photorespiration (i.e., at high compared with ambient [CO 2 ]) was much lower than that of serine which is against a low N availability (Sulpice et al., 2013;Noguchi et al., 2015). On the other hand in agreement with previous studies the amount of minor amino acids (e.g., threonine, valine, cysteine, methionine, lysine, leucine, etc.) decrease at elevated [CO 2 ] (Noguchi et al., 2015).
Summarizing the decrease in the amount of several amino acids under elevated CO 2 , reported in our study may have several causes. Some of them, decreased probably due to the inhibition of photorespiration under elevated CO 2 (glycine and serine). The rest of the amino acids could have been reduced due to the inhibition of N assimilation; this has been reported in several plants under elevated CO 2 , such as Arabidopsis and wheat (Aranjuelo et al., 2013;Noguchi et al., 2015;Vicente et al., 2016b). The inhibition of N assimilation under elevated CO 2 is not completely understood, but it could be due, as reported in Vicente et al. (2016b) and in agreement with previous studies (Bloom et al., 2010;Aranjuelo et al., 2013), to (i) lower rates of photorespiration in elevated CO 2 that decrease the availability of reductant in the cytosol, and inhibition of NO − 2 influx into chloroplasts; or to increased demand for NO − 3 to match the increase in net photosynthetic assimilation in elevated CO 2 , leading to a decline of the metabolic NO − 3 pool that restricts induction of nitrate reductase and thus NO − 3 assimilation. Interestingly tyrosine was the only amino acid clearly increasing as response to a high [CO 2 ]. Moreover, under water stress, high [CO 2 ] also increased phenylalanine. Both are aromatic amino acids, which have low N/C ratios, involved as precursors in the synthesis of anthocyanin. In fact an increase in both amino acids has been reported in Arabidopsis leaves at elevated [CO 2 ] in NH + 4 + NO − 3 media (Noguchi et al., 2015) which is our case (i.e., nutrient solution). Cysteine was strongly decreased under elevated [CO 2 ] regardless of water regime. It constitutes the first stable product of the sulfur assimilation, and acts as precursor of the majority of organic sulfur compounds (Hawkesford et al., 2012).
Independently of [CO 2 ], the levels of all amino acids were correlated (Supplementary Material Figure S2). However, the relationships between amino acids with other physiological traits varied depending on [CO 2 ]. It was especially remarkable that amino acid contents correlated positively with SB under high [CO 2 ] and negatively with traits related to water status (RWC, w , δ 13 C, and T), while under atmospheric [CO 2 ] the correlations between amino acids and other traits were scarce. These findings suggest that amino acid pool greatly influences biomass accumulation under elevated [CO 2 ] and its amount is influenced by leaf and plant water status in a highly-dependent manner.
Decreases in transpiration associated with mass flow of soil solution, have been proposed to limit plant N acquisition. In an experiment with cottonwood, where relative humidity and atmospheric [CO 2 ] were manipulated to alter the transpirative stream, N gain was positively correlated across all treatments with root mass, and a significant portion of the remaining variation (44%) was positively related to transpiration per unit root mass (McDonald et al., 2002). Thus, decreases in plant N concentration under water stress are attributable in part to associated decreases in g s and transpiration. However, other mechanisms may be involved for the decrease in plant N under [CO 2 ] enrichment. Thus, in our study no clear differences existed in the rates of g s measured at different [CO 2 ] and even transpiration was increased at high [CO 2 ] relative to atmospheric [CO 2 ]. Besides the potential effect of diminishing transpirative stream, [CO 2 ] enrichment is reported to inhibit the assimilation of nitrate into organic nitrogen compounds (Bloom et al., 2010;Vicente et al., 2016b). This inhibition may be largely responsible for [CO 2 ] acclimation, that is, the decrease in photosynthesis and growth of C 3 plants after long exposures to [CO 2 ] enrichment. Different studies have shown that the effect of elevated [CO 2 ] on reducing N content was related at the transcript level to a down-regulation of genes encoding for Rubisco subunits and N-assimilation enzymes (GS1 and GS2), indicating a co-regulation of primary C and N metabolism (Stitt and Krapp, 1999;Vicente et al., 2015Vicente et al., , 2016aMedina et al., 2016).
As a consequence of the effect of water regime on N concentration, the trait best correlated with total biomass within each of the two [CO 2 ] was the total leaf N concentration (and amino acid contents), although the correlation was stronger at high [CO 2 ]. Overall, the results showed that greater plant growth across water regimes was linked to an increase in shoot N concentration associated to a higher transpiration stream, even when changes in N metabolism appear also involved (Hirel et al., 2007;Medina et al., 2016). A study of Peñuelas et al. (1995) on sweet pepper also found that the effect of [CO 2 ] and water regime was dependent on the accumulated N supply. All these studies have plants grown under hydroponic conditions in common, where irrigation and fertilization are provided together through the nutrient solution under relatively low photosynthetic photon flux density.
The leaf N concentration was negatively correlated with δ 13 C within each [CO 2 ] across water regimes, whereas nitrogen concentration and g s correlated positively (but in a weaker manner). Moreover, δ 13 C correlated with the total shoot nitrogen content in a weaker manner than the nitrogen concentration alone. In fact, both the δ 13 C and nitrogen concentration were expressed on a dry matter basis, which may explain their better correlation. δ 13 C is not just an indicator of water use efficiency (Farquhar et al., 1982(Farquhar et al., , 1989 but it is also strongly negatively affected by the amount of available water (Araus et al., 1997(Araus et al., , 2003 and therefore when analyzed on a dry matter basis it may be considered an indicator on the total (i.e., through time) water used by the plant. Moreover, in the Principal Component Analysis the nitrogen concentration was placed clearly opposite to δ 13 C. In other words, the negative relationship between δ 13 C and nitrogen concentration on a dry leaf basis may be understood as the nitrogen concentration in the leaves being, at least in part, the consequence of the amount of water transpired by the plant. The Path Analysis further supported the direct contribution of the nitrogen concentration in determining total biomass as well as the positive role of an increased water use (assessed through a lower δ 13 C of the dry matter) on the N accumulation. This model worked better under high [CO 2 ]. In sweet pepper growing under hydroponic conditions, strong negative correlations across salinity levels for leaf δ 13 C with both g s and nitrate accumulation have been reported, with the correlations being higher at elevated [CO 2 ] compared with atmospheric [CO 2 ] (del Amor, 2013).

CONCLUSION
The signatures of the three different stable isotopes are significantly affected by water regime, [CO 2 ], and genotypic effects. However, the results do not support the use of δ 18 O as an indicator of the effect of growing conditions and genotypes on plant growth. This study proves that the effect of water regime on sweet pepper growth in a hydroponic system is caused by changes in the amount of nitrogen assimilated, which is associated with the plant's water use. In that sense, the role of δ 13 C in explaining differences in plant growth across water regimes appears mediated via its direct relationship with N accumulation in leaves, particularly at high [CO 2 ]. However, our study does not support stomatal closure as being associated with an elevated [CO 2 ]-induced reduction in N concentration in the shoot.

AUTHOR CONTRIBUTIONS
FdA and JA conceived and designed the experiments. MP, MS, and GO-A contributed to the experimental work. SY and RV contributed to the data analysis and preparation of tables and figures. MS wrote the paper under the supervision of FdA and JA and all three revised the manuscript. All authors read and approved the final manuscript.

ACKNOWLEDGMENTS
This study was supported in part by the European Regional Development Fund (ERDF) 80%-Región de Murcia (FEDER 1420-07) and the 2014-SGR-628 AGAUR project, Generalitat de Catalunya, Spain. MP is the recipient of a pre-doctoral fellowship from INIA. The authors thank Miguel Marín for his technical assistance.