Ana Zangirolame Gonçalves
Helenice Mercier- Departamento de Botânica, Instituto de Biociências, Universidade de São Paulo (USP), São Paulo, Brazil
Most epiphytes are found in low-nutrient environments with an intermittent water supply. To deal with water limitation, many bromeliads perform crassulacean acid metabolism (CAM), such as Guzmania monostachia, which shifts from C3 to CAM and can recycle CO2 from the respiration while stomata remain closed during daytime and nighttime (CAM-idling mode). Since the absorbing leaf trichomes can be in contact with organic (urea) and inorganic nutrients (NO3−, NH4+) and the urea hydrolysis releases NH4+ and CO2, we hypothesized that urea can integrate the N and C metabolism during periods of severe drought. Under this condition, NH4+ can be assimilated into amino acids through glutamine synthetase (GS), while the CO2 can be pre-fixated by phosphoenolpyruvate carboxylase (PEPC). In this context, we evaluated the foliar transcriptome of G. monostachia to compare the relative gene expression of some genes involved with CAM and the N metabolism when bromeliads were submitted to 7days of drought. We also conducted a controlled experiment with an extended water deficit period (21days) in which bromeliads were cultivated in different N sources (urea, NH4+, and NO3−). Our transcriptome results demonstrated an increment in the expression of genes related to CAM, particularly those involved in the carboxylation metabolism (PEPC1, PPCK, and NAD-MDH), the movement of malate through vacuolar membrane (ALMT9), and the decarboxylation process (PEPCK). Urea stimulated the expression of PEPC1 and ALMT9, while Urease transcripts increased under water deficit. Under this same condition, GS1 gene expression increased, indicating that the NH4+ from urea hydrolysis can be assimilated in the cytosol. We suggest that the link between C and N metabolism occurred through the supply of carbon skeleton (2-oxoglutarate, 2-OG) by the cytosolic isocitrate dehydrogenase since the number of NADP-ICDH transcripts was also higher under drought conditions. These findings indicate that while urea hydrolysis provides NH4+ that can be consumed by glutamine synthetase-cytosolic/glutamate synthase (GS1/GOGAT) cycle, the CO2 can be used by CAM, maintaining photosynthetic efficiency even when most stomata remain closed 24h (CAM-idling) as in the case of a severe water deficit condition. Thus, we suggest that urea could be used by G. monostachia as a strategy to increase its survival under drought, integrating N and C metabolism.
Introduction
Most epiphytic bromeliads inhabit oligotrophic environments with intermittent water and nutrient availability (Benzing, 2008). Periods of severe drought are frequent in the canopy, and it is important to understand how each species copes with harsh environments. One of the principal morphologic/anatomic features of bromeliads to face water and nitrogen (N) scarcity is the formation of a tank by the arrangement of leaves, such as a rosette; this structure is capable of accumulating water and dissolved nutrients, including organic N source, such as urea provided by the excrement of animals (Benzing and Bennett, 2000; Cambuí and Mercier, 2006). Highly efficient absorbing trichomes are distributed along the leaf blade following a gradient in which the basal part has the largest density and also the highest degree of cellular structure specialization (Kleingesinds et al., 2018). Alternatively, the apex has fewer trichomes and a greater number of stomata, thus playing an important role to take up CO2 from the atmosphere. A relevant adaptative physiological characteristic of many epiphytic tank bromeliads to face water stress is the crassulacean acid metabolism (CAM) photosynthesis because this process confers a higher carboxylation efficiency (Luttge, 2004) and a greater Rubisco (Ribulose-1,5-Biphosphate Carboxylase-Oxygenase) efficiency (Matiz et al., 2013; Gonçalves et al., 2020a).
To deal with nutrient scarcity of the canopy, the leaves of epiphytic tank bromeliads acquired the function of the roots of terrestrial plants. The leaf became the principal vegetative organ responsible for plant survival (Zotz, 2013), showing a spatial and functional division due to its ability to uptake nutrients, mainly in the basal region that is in contact with the tank solution, and to perform photosynthesis, mainly in the apical portion (Freschi et al., 2010; Takahashi and Mercier, 2011). Highly efficient uptake and great storage capacity of nutrients is another important physiological feature of the epiphytic tank bromeliads (Zotz and Asshoff, 2010). Effective leaf N uptake has been demonstrated for NH4+ and urea, which are the preferable N sources for Vriesea gigantea and Guzmania monostachia (Inselsbacher et al., 2007; Gonçalves et al., 2020b). Urea, after its intracellular hydrolysis by urease enzyme, provides NH4+ and CO2 to the leaf tissues (Miflin and Lea, 1980). Low and high affinity N transporters have been detected along the leaf blade of G. monostachia. At the base, for example, a significant increase of gene expression of NH4+ (AMT1.2) and urea (DUR3 and PIP1.2) transporters were measured (Gonçalves et al., 2020b). Considering N assimilation, the highest activity of glutamine synthetase (GS) was found at the apex, while urease can be found at the base, indicating the importance of the apical foliar portion for amino acid production (Gonçalves et al., 2020b). Transcriptional analysis of G. monostachia leaf showed that the apex and middle portion had differentially expressed genes related to amino acid metabolism and photosynthesis, while the base showed responses to nutrients, particularly N (Mercier et al., 2019).
In its natural environment, G. monostachia can survive several months with the empty tank (Maxwell et al., 1995). To deal with water limitation, many bromeliad tanks exhibit CAM. This metabolism is a photosynthetic pathway that can be divided into four phases, in which CO2 is obtained overnight when stomata open (phase I) and it may be assimilated by the enzyme phosphoenolpyruvate carboxylase (PEPC), providing oxaloacetic acid (OAA) and subsequently forming malate by the cytosolic malate dehydrogenase (NAD-MDH; Luttge, 2004). Afterwards, malate is transported into the vacuoles through the aluminum-activated malate transporter (ALMT) family (Pereira et al., 2018) where it can be accumulated. At dawn (phase II), both PEPC and Rubisco are able to fix carbon while the stomata are still open (Matiz et al., 2013). In the subsequent light period (phase III), malic acid is remobilized to the cytosol and decarboxylated through NAD-MDH and phosphoenolpyruvate carboxykinase (PEPCK), producing CO2 that is consumed by Rubisco while stomata remain closed (Luttge, 2004). At dusk (phase IV), PEPC activity increases and Rubisco decreases while stomata reopen enabling CO2 assimilation (Luttge, 2004). In G. monostachia, for example, the isoform ALMT9 is involved in the transport of both malate and fumarate to the vacuole (Pereira et al., 2018). During the subsequent light period, stomata remain closed while malic acid is remobilized from the vacuoles to the cytosol where it can be decarboxylated to phosphoenolpyruvate (PEP) and CO2 which is used for refixation in the C3 cycle (Luttge, 2004). Guzmania monostachia, as a facultative CAM species, can shift from C3 to CAM under drought and absorbs CO2 in the dark period (Medina et al., 1977). There is a variant of normal full CAM called CAM-idling in which stomata remain closed during night and day when water availability is unfavorable, like in a dry season. In this mode of CAM, respiratory CO2 is refixed in the dark period and recycled to carbohydrate during the light period (Lüttge, 2006). Guzmania monostachia showed a capacity to express a reversible CAM-idling photosynthesis under drought (Freschi et al., 2010), and it was estimated that the accumulation of nocturnal acid malic can contribute to approximately 50% of the total carbon assimilated under water deficit (Pikart et al., 2018). During drought, atmospheric CO2 absorption by the leaves is restricted due to stomata closure, causing a significant decrease in CO2 availability. Interestingly, the delivery of 13C-labeled urea to the leaves under water deficit enriched the 13C-malate pool of the apical portion of V. gigantea leaves (Matiz et al., 2017), suggesting that this organic N source could be an important source of CO2 for bromeliads under water stress (Matiz et al., 2017).
In this context, our goal was to investigate the connection between the N and C metabolism through urea nutrition during drought periods considering the G. monostachia bromeliad as a plant-model. Since urea hydrolysis releases CO2 and NH4+, we hypothesized that the interaction between drought and urea nutrition can activate the CAM pathway through the carboxylation metabolism (PEPC1, PPCK, and NAD-MDH), the movement of malate through vacuolar membrane (ALMT9), and the decarboxylation process (PEPCK), while NH4+ can be assimilated into glutamine through the activation of the GS/GOGAT cycle. Ultimately, the cytosolic isocitrate desidrogenase (NADP-ICDH) can participate in this process, providing 2-oxoglutarate (2-OG) from isocitrate. This 2-OG is then imported into the chloroplasts, where it can be used as the carbon skeleton required for amino acid synthesis.
Materials and Methods
Organism and Experimental Design
Gene Expression Related to CAM and N Metabolism
This experiment was conducted with adult G. monostachia (about 2.5years) obtained previously by in vitro propagation as described in Rodrigues et al. (2016). Until bromeliads reached this mature-vegetative age, they were kept in a greenhouse of the Department of Botany at the University of São Paulo (São Paulo, Brazil), where they received tap water daily and nutritive solution modified with half of the original concentration of macronutrients (Knudson, 1946) and with micronutrients (Murashige and Skoog, 1962) each, for 15days. These plants were kept in pots (14cm in diameter, 11cm high) containing a mixture of Pinus sp. bark and commercial organic substrate (Tropstrato®). These mature-vegetative bromeliads had around 20 leaves each and the oldest leaves measuring around 20cm length.
Previously to the experiment, bromeliads had their tanks emptied and were transferred to the growth chamber where they were kept for 1month under the following controlled conditions: 12-h photoperiod, photosynthetic photon flux of 250μmolm−2s−1, air temperature of 27/22°C day/night, and relative humidity of 60/70% day/night. During this period, all plants received distilled water daily and were not fertilized. After 30days, bromeliads were divided into two treatments (n=3 replicates): (1) bromeliads were subjected to 7days of water deficit (WD), and (2) bromeliads received distilled water for 7days (W). While the water treatment received distilled water for 7days, bromeliads of the water deficit treatment remained with the tank emptied throughout this period. Subsequently, apical and basal leaf portions of the 8th to 12th innermost nodes were harvested at 7am and 7pm, immediately frozen with liquid nitrogen, and stored at −80°C for RNA extraction.
CAM Modulation Under Different N Sources
In order to evaluate the interaction between CAM and nitrogen metabolism, the experiment was carried out with adult bromeliads in the same controlled conditions of temperature, light, relative humidity, and photoperiod as described above. After 30days, bromeliads received nutrient solutions with a total of 5mM of nitrogen, according to the following treatments (n=8 replicates per treatment): (1) without nitrogen source (control), (2) Ca(NO3)2 (nitrate), (3) (NH4)2SO4 (ammonium), and (4) urea (urea). The nutrient solution was prepared with half of the original concentration of macronutrients (Knudson, 1946; modified without nitrogen sources) and with micronutrients (Murashige and Skoog, 1962). Each experimental bromeliad received 30ml of the nutrient solution inside its tank once every 7days for 30days (totaling four applications).
Twenty-four hours after the last nutrient-solution application, bromeliads had their tanks emptied and were separated into two new treatments for CAM induction (n=4 replicates per treatment): (1) receiving distilled water to keep the tank full of water (W) and (2) water deficit condition, in which bromeliads did not receive distilled water (WD). Bromeliads remained in these conditions for 21days as this environmental condition was shown to be the best to induce the transient shift of the photosynthesis from C3 to CAM-idling in the apical part of the leaves of G. monostachia (Pikart et al., 2020). At the end of the experiment, the apical leaf portion of the 8th to 12th innermost nodes were harvested at 7am in order to analyze the relative water content (RWC) in the leaves (RWC), the activity of the enzyme PEPC, the expression of the CAM-specific phosphoenolpyruvate carboxylase gene (PEPC1), and the tonoplast aluminum-activated malate transporter gene (ALMT9). For the quantification of nocturnal acid accumulation (ΔH+), bromeliad leaves were collected at 7am and 7pm.
RNA Extraction, RNA-Seq de novo Assembly, Annotation, and Analysis
Total RNA was extracted using the PureLink® RNA Mini Kit (Ambion) according to the manufacturer’s protocol. RNA concentration, quality, and integrity were checked as described in Mercier et al. (2019). Samples with integrity number ≥6.5 were pooled to generate libraries using TruSeq RNA Sample Preparation Kit, Set A (Illumina Inc., United States). Paired-end sequences were generated as described by Mercier et al. (2019), and only high-quality sequences (with average PhredScore over 24) were used. Sequence normalization and de novo transcriptome assembly were performed using Trinity (Grabherr et al., 2011; Haas et al., 2013) as described in Mercier et al. (2019), and the reads were deposited at NCBI at the Bioproject ID PRJNA532595. Guzmania monostachia transcripts were annotated using the Viridiplantae (taxa ID 33090, VP) and Monototyledons (taxa ID 4447, MC) database through the BLAST suite (Camacho et al., 2008) and Blast2GO1 (Conesa et al., 2005). The annotated genes related to NR (nitrate reductase), NiR (nitrite reductase), Urease (urease), GS2 (glutamine synthetase-chloroplastic), GS1 (glutamine synthetase-cytosolic), Fd-GOGAT (ferredoxin-glutamate synthase), NADH-GOGAT (NADH-glutamate synthase), NADP-ICDH (NADP-isocitrate dehydrogenase), PEPC1 (phosphoenolpyruvate carboxylase), PPCK1 (phosphoenolpyruvate carboxylase kinase), NAD-MDH (NAD-malate dehydrogenase), ALMT9 (tonoplast aluminum-activated malate transporter), and PEPCK were selected and analyzed through pairwise differential expression (DE) as described by Mercier et al. (2019), and the heatmap was generated using the Morpheus tool2 in order to compare apex and base, samples collected at 7am and 7pm, and water (W) and water deficit (WD) treatments.
Relative Water Content
The RWC of the plants was determined by collecting 10 leaf discs of a known area, from which the values of fresh mass (MF), and, subsequently, the discs remained in distilled water for 24h to obtain the turgid mass (MT). After this period, the discs were dried at 60°C for 48h to determine the dry mass (DM). These results were used in the following equation: RWC (%)=[(MF-MS)/(MT-MS)]*100 (Weatherley, 1950).
Phosphoenolpyruvate Carboxylase Enzyme Activity
The determination of PEPC activity was performed by the in vitro method as described by Nievola et al. (2001) and adapted for G. monostachia (Freschi et al., 2010; Pereira et al., 2013). In this determination, samples of the apical portion of leaves (1g) were collected, macerated in liquid nitrogen and immersed in 5ml of the extraction buffer (pH 8), composed of 200mM Tris–HCl, 10mM MgCl2, 5mM DTT, 1mM EDTA, 0.5% (w/v) BSA, and 10% (v/v) glycerol. The samples were centrifuged for 2min at 15,000g. The supernatant was collected and passed through a SEPHADEX PD10/G-25 column. PEPC activity was analyzed in 2ml of reaction buffer containing 50mM Tris–HCl (pH 8.0), 1mM DTT, 10mM MgCl2, 100mM NaHCO3, 20mM NADH, and 3mM phosphoenolpyruvate (PEP). The reaction was started by adding 200μl of the extract. The consumption of NADH was quantified in a spectrophotometer (340nm) at the initial moment of the reaction and after 4min. The results were expressed in micromoles of NADH consumed per minute per gram of dry mass (μmolmin−1g−1).
Nocturnal Acid Accumulation
The apical leaf portion was analyzed for the degree of CAM expression according to the nocturnal acid accumulation (ΔH+) using the titratable acidity methodology (Matiz et al., 2017). For this, the leaves of the bromeliads were cut into small fragments of 1cm2, and 0.1g was macerated in liquid nitrogen, homogenized in 500μl of MCW solution, and incubated at 60°C for 30min. The MCW solution was prepared in a proportion of 12:5:1 of methanol:chloroform:water, respectively. Subsequently, the tubes received 500μl of distilled water and were centrifuged for 10min at 16,000g. After centrifugation, 1.5ml of the supernatant was collected and completely dried at 60°C. After this procedure, the supernatant was collected and stored at −20°C until the titratable acidity procedure.
For the titratable acidity procedure, a solution of phenolphthalein (10mg in 10ml of ethanol) was diluted in 3ml of boiled distilled water to be pipetted into an Elisa® plate. In the microplate, the following volumes per sample were pipetted in this order: 100μl of boiled distilled water, 100μl of phenolphthalein, 5μl of the respective sample, and 100μl of the curve solution. All samples received all points on the curve in order to calculate their H+ concentrations.
Curve points were obtained as follows: (A) first point with 0.5μmol of H+ (250μl of 0.2M NaOH and 9.75ml of boiled distilled water), (B) second point with 0.4μmol of H+ (200μl of 0.2M NaOH and 9.80ml of boiled distilled water), (C) third point with 0.3μmol of H+ (150μl of 0.2M NaOH and 9.85ml of boiled distilled water), (D) fourth point with 0.2μmol of H+ (100μl of 0.2M NaOH and 9.90ml of boiled distilled water), (E) fifth point with 0.16μmol of H+ (80μl of 0.2M NaOH and 9.92ml of boiled distilled water), (F) sixth point with 0.12μmol of H+ (60μl of 0.2M NaOH and 9.94ml of boiled distilled water), (G) seventh point with 0.08μmol of H+ (40μl of 0.2M NaOH and 9.96ml of boiled distilled water), and (H) eighth point with 0.04μmol of H+ (20μl of 0.2M NaOH and 9.98ml of boiled distilled water).
Gene Expression
The genes encoding the expression of the CAM-specific PEPC enzyme (PEPC1) and the tonoplast aluminum-activated malate transporter (ALMT9) which presented the most discrepant gene expression in the transcriptome of G. monostachia (Mercier et al., 2019) were selected for the primer design using the PerlPrimer software3 (Supplementary Table S2). Among all PEPC isoforms found in the transcriptome, PEPC1 was selected because its genomic sequence was more similar to the genomic sequences of the CAM-specific PEPC enzyme comparing to other CAM species (Beltran and Smith, unpublished data). In addition, ALMT9 isoform was selected because it showed to be involved in malate transport to the vacuole in G. monostachia (Pereira et al., 2018). The reference genes (beta-amylase, BAM; betahydroxysteriod-dehydrogenase decarboxylase, HSDD; Supplementary Table S2) were selected from the G. monostachia transcriptome as described by Gonçalves et al. (2020b).
Total RNA of bromeliad leaves was extracted from 300mg using Trizol® reagent (Invitrogen®) with the PureLink® RNA mini kit (Life Technologies®), according to the manufacturer instructions. RNA concentration, its integrity analysis, and the removal of DNA contamination have been described previously by Gonçalves et al. (2020b). The cDNA was synthesized using the SuperScript IV (Invitrogen®) reaching a final volume of 20μl according to the manufacturer instructions. The PCRs were performed using the StepOnePlusTM® Real-Time PCR (Applied Biosystems®), and products were detected using SYBR® Green PCR Master Mix (Applied Biosystems®). The analyses of gene expressions were performed using BestKeeper (Pfaffl et al., 2004).
Statistical Analyses
The response variables of the RWC, the activity of the PEPC enzyme, the nocturnal acid accumulation (ΔH+), and the expression of the genes PEPC1 and ALMT9 were compared using ANOVA, and Tukey HSD post hoc tests were used for pair-wise comparisons. All statistical analyses were performed with the statistical platform R (R Core Team, 2020).
Results
Gene Expression Related to CAM
Some genes related to CAM showed markedly greater expression in the leaf apex of plants subjected to WD, especially PPCK1 and ALMT9 (Figure 1) when compared to the base. PEPC1 showed a high gene expression throughout the entire leaf; however, its highest expression was observed at 7am under WD (Figure 1). NAD-MDH showed greater number of transcripts at 7am independent of the leaf portion, especially when the plant was subjected to WD (Figure 1). The highest expression of PEPCK was observed at the leaf apex at 7am under WD (Figure 1).
Figure 1. Heatmap indicating log2 (fold-change) of the genes NR (nitrate reductase), NiR (nitrite reductase), Urease (urease), GS2 (glutamine synthetase-chloroplastic), GS1 (glutamine synthetase-cytosolic), Fd-GOGAT (Ferredoxin- glutamine-2-oxoglutarate aminotransferase), NADH-GOGAT (NADH- glutamine-2-oxoglutarate aminotransferase), NADP-ICDH (NADP-isocitrate dehydrogenase), PEPC1 (phosphoenolpyruvate carboxylase), PPCK1 (phosphoenolpyruvate carboxylase kinase), NAD-MDH (NAD-malate dehydrogenase), ALMT9 (tonoplast aluminum-activated malate transporter), and phosphoenolpyruvate carboxykinase (PEPCK) in order to compare the apical and basal leaf portions of Guzmania monostachia submitted to 7days of water deficit (WD) or receiving distilled water (W) for the same period. Samples were collected at 7am and 7pm in order to compare gene expressions in different periods of the day. More details about the methodology are described in the section “Gene Expression Related to CAM and N Metabolism.” Heatmap colors varied from blue to red representing low to high gene expression.
Gene Expression Related to N Metabolism
Both NR and NiR had the highest relative gene expression at 7pm in the apex of the bromeliad leaves subjected to WD (Figure 1). Similarly, the quantity of urease transcripts in the apical leaf portion were the highest at 7pm under drought (Figure 1). GS2 showed greater gene expression at the leaf apex compared to the base, independent of the treatments and time of harvest; however, for the base the GS2 expression was stronger at 7pm, independent of the water availability (Figure 1). Similar to GS2, GS1 showed higher gene expression at the leaf apex, especially at 7am, when plants were subjected to WD (Figure 1). The gene expression of Fd-GOGAT was greater at the leaf apex compared with the base, while at the basal part its expression at 7pm under WD was more intense compared with the other treatments (Figure 1). The NADH-GOGAT showed higher gene expression at 7pm with the apex showing a greater quantity of transcripts than the basal part (Figure 1). The transcripts of NADP-ICDH were more abundant in the apical leaf portion compared to the base, mainly under water deficiency.
CAM Modulation Under Different N Sources
Bromeliads that were subjected to WD and received urea, ammonium or did not receive N showed the lowest RWC in the leaves, while the plants that were fertilized with nitrate had the highest water content (Supplementary Table S1; Figure 2). Hydrated plants that were fertilized with nitrate showed the highest RWC and, after 21days of WD, those bromeliads kept 50% of their water content (Supplementary Table S1; Figure 2). Consistent with these results, we observed that the PEPC enzyme activity and the nocturnal acid accumulation were higher in bromeliads that acquired urea and were submitted to WD, while bromeliads that took up nitrate with WD showed lower PEPC activity and nocturnal acid accumulation (Supplementary Table S1; Figures 3, 4). In addition, bromeliads that received ammonium or did not receive N at all and were subjected to WD presented intermediate PEPC activity and nocturnal acid accumulation (Supplementary Table S1; Figures 3, 4). On the other hand, bromeliads showed the lowest PEPC activities and nocturnal acid accumulation when they received water regardless of nutritional treatment (Supplementary Table S1; Figures 3, 4).
Figure 2. Relative water content in the apical leaf portion of G. monostachia submitted to the following treatments for 1 month: nutrient solution without nitrogen source (control), Ca(NO3)2 (nitrate), (NH4)2SO4 (ammonium), and urea (urea). Twenty-four hours after the last application of the nutrient solutions, bromeliads were separated into two new treatments for 21days: (1) receiving distilled water (water) and (2) water deficit condition in which the bromeliads did not receive distilled water (water deficit). Bars indicate standard error. Letters indicate statistical differences between nutritional treatments, and asterisk (*) indicates differences between water vs. water deficit (ANOVA/Tukey HSD post hoc test, α=0.05). Letters are showing the statistical differences among nutritional treatments also under water treatments. Since the nutritional vs. water treatment was statistically significant (Supplementary Table S1), we represented each nutritional treatment with a single letter.
Figure 3. Activity of the phosphoenolpyruvate carboxylase (PEPC) enzyme in the apical leaf portion of G. monostachia submitted to the following treatments for 1month: nutrient solution without nitrogen source (control), Ca(NO3)2 (nitrate), (NH4)2SO4 (ammonium), and urea (urea). Twenty-four hours after the last application of the nutrient solutions, bromeliads were separated into two new treatments for 21days: (1) receiving distilled water (water) and (2) water deficit condition in which the bromeliads did not receive distilled water (water deficit). Bars indicate standard error. Letters indicate statistical differences between nutritional treatments, and asterisk (*) indicates differences between water vs. water deficit (ANOVA/Tukey HSD post hoc test, α=0.05). Letters are showing the statistical differences among nutritional treatments also under water treatments. Since the nutritional vs. water treatment was statistically significant (Supplementary Table S1), we represented each nutritional treatment with a single letter.
Figure 4. Nocturnal acid accumulation (ΔH+) in the apical leaf portion of G. monostachia submitted to the following treatments for 1month: nutrient solution without nitrogen source (control), Ca(NO3)2 (nitrate), (NH4)2SO4 (ammonium), and urea (urea). Twenty-four hours after the last application of the nutrient solutions, bromeliads were separated into two new treatments for 21days: (1) receiving distilled water (water) and (2) water deficit condition in which the bromeliads did not receive distilled water (water deficit). Bars indicate standard error. Letters indicate statistical differences between nutritional treatments, and asterisk (*) indicates differences between water vs. water deficit (ANOVA/Tukey HSD post hoc test, α=0.05). Letters indicate statistical differences between nutritional treatments, and asterisk (*) indicates differences between water vs. water deficit (ANOVA/Tukey HSD post hoc test, α=0.05). Letters are showing the statistical differences among nutritional treatments also under water treatments. Since the nutritional vs. water treatment was statistically significant (Supplementary Table S1), we represented each nutritional treatment with a single letter.
Concerning the gene expression of PEPC1, a significant increase in its expression was observed in bromeliads that acquired urea under WD, followed by plants that were fertilized with ammonium under WD (Supplementary Table S1; Figure 5A). On the other hand, plants that received nitrate or did not receive N and were subjected to WD showed lower PEPC1 gene expressions comparing to bromeliads with other nutritional treatments (Figure 5A). Additionally, hydrated plants presented the lowest PEPC1 gene expressions regardless the nitrogen they received (Supplementary Table S1; Figure 5A). Consistently, ALMT9 showed its highest expression in bromeliads that acquired urea and were submitted to WD, and the lowest expression was observed in plants that were nourished with nitrate followed by WD (Supplementary Table S1; Figure 5B). Bromeliads that received ammonium or did not have contact with any N sources and were submitted to WD showed intermediate ALMT9 gene expressions (Supplementary Table S1; Figure 5B). Likewise, the gene expressions of ALMT9 were lower in hydrated plants comparing to WD plants (Supplementary Table S1; Figure 5B).
Figure 5. Gene expression (fold change) of the (A) phosphoenolpyruvate carboxylase gene (PEPC1) and (B) the tonoplast aluminum-activated malate transporter gene (ALMT9) in the apical leaf portion of G. monostachia submitted to the following treatments for 1month: nutrient solution without nitrogen source (control), Ca(NO3)2 (nitrate), (NH4)2SO4 (ammonium), and urea (urea). Twenty-four hours after the last application of the nutrient solutions, bromeliads were separated into two new treatments for 21days: (1) receiving distilled water (water) and (2) water deficit condition in which the bromeliads did not receive distilled water (water deficit). Bars indicate standard error. Letters indicate statistical differences between nutritional treatments, and asterisk (*) indicates differences between water vs. water deficit (ANOVA/Tukey HSD post hoc test, α=0.05). Letters indicate statistical differences between nutritional treatments, and asterisk (*) indicates differences between water vs. water deficit (ANOVA/Tukey HSD post hoc test, α=0.05). Letters are showing the statistical differences among nutritional treatments also under water treatments. Since the nutritional vs. water treatment was statistically significant (Supplementary Table S1), we represented each nutritional treatment with a single letter.
Discussion
Our results indicate that drought had a positive effect on the C3-CAM shift at the apical leaf part, increasing the expression of genes related to carboxylation metabolism (PEPC1, PPCK, and NAD-MDH), movement of malate through vacuolar membrane (ALMT9) and decarboxylation process (PEPCK; Figure 6). Together with water deficiency, urea hydrolysis seemed to provide CO2 for malate production and accumulation into vacuoles since we observed that urea stimulated the expression of PEPC1 and ALMT9 genes (Figure 6). At the same time, urea hydrolysis provides NH4+(Miflin and Lea, 1980) that can be assimilated into Gln since our results showed increased GS1 gene expression under water deficiency (Figure 6). Ultimately, NADP-ICDH can supply the carbon skeleton to Fd-GOGAT functioning (Figure 6). Because G. monostachia is an epiphytic bromeliad subjected to the intermittence of water and nutrients in its natural environment, we reinforce that urea nutrition obtained in nature by the symbiotic interaction of bromeliads with anurans may bring advantages since the CO2 can maintain the photosynthetic efficiency even when most stomata remain 24h closed, while NH4+ can maintain amino acid and protein production.
Figure 6. A simplified schematic showing our hypothesis about nitrogen metabolism and its interaction with carbon metabolism in a water stress-induced CAM-idling bromeliad, G. monostachia, submitted to urea nutrition. Urea is taken up at the leaf basal part and after its hydrolysis generates NH4+ and CO2. At the green leaf part and mainly in the apex portion during the night, the CO2 may be assimilated through the crassulacean acid metabolism (CAM) pathway. Urea nutrition associated with water stress positively stimulated PEPC1 and ALMT9 gene expressions and PEPC1 activity, resulting in a higher nocturnal vacuolar acidity. NH4+, in turn, mainly after entering a cell of the apical leaf region may be assimilated during the light–dark period in the cytosol mediated by the enzyme GS1 forming glutamine (Gln). Therefore, glutamine and glutamate (Glu) biosynthesis may occur in separate compartments. The mitochondria may supply the carbon skeleton for Fd-GOGAT functioning. The key C-compound, 2-oxoglutarate (2-OG), may derive from sugar respiration reactions. In this pathway, citrate is used to generate the isocitrate required for 2-OG synthesis, allowing the net Glu production in the chloroplast and then other amino acid formation. Abbreviations for enzymes: PEPC1, phosphoenolpyruvate carboxylase 1; PPCK, phosphoenolpyruvate carboxylase kinase; NAD-MDH, NAD-malate dehydrogenase; PEPCK, phosphoenolpyruvate carboxykinase; Rubisco, ribulose-1,5-bisphosphate carboxylase/oxygenase; GS1, glutamine synthetase 1; Fd-GOGAT, glutamine-2-oxoglutarate aminotransferase; and NADP-ICDH, NADP-isocitrate dehydrogenase. Abbreviations for transporters: ALMT9, tonoplast aluminum-activated malate transporter 9; DUR3, urea transporter; AMT3.1, ammonium transporter; TIP2.1, tonoplast intrinsic protein; PIP1.2, plasma membrane intrinsic protein; AA, amino acid.
C3/CAM-Idling Comparison Through RNA-Seq
Our bromeliads were submitted to drought, and this environmental condition was shown to be the best to induce the transient shift of the photosynthesis from C3 to CAM-idling in the apical part of the leaves of G. monostachia (Freschi et al., 2010; Mercier et al., 2019). However, the molecular mechanisms underlying the stress induction of the transition from C3 to CAM remain poorly understood. Recently, this was observed for the facultative C3-CAM plant Mesembryanthemum crystallinum (Aizoaceae), whose guard cells themselves can shift from C3 to CAM, increasing expression of PEPC1 and PPCK1, which encode key enzymes of CAM photosynthesis (Kong et al., 2020).
Comparison of the differential patterns of mRNA abundance in the apical leaf portion between C3 (W treatment) and CAM-idling (WD treatment) of G. monostachia tissues suggests that the C3-CAM switch may be determined by the greater expressions of some genes involved in the carboxylation metabolism (PEPC1, PPCK, and NAD-MDH), the movement of malate through vacuolar membrane (ALMT9) and the decarboxylation process (PEPCK), independent of harvest time (7am or 7pm). It is interesting to observe that in the basal portion, although there was an increase in the number of PEPC1 transcripts in drought condition, there was practically no expression of PPCK1. This result indicates that the enzyme PEPC1 could not be intensely activated through the phosphorylation process, and, consequently, the CO2 fixation could not happen through the CAM pathway in this part of the leaf. Besides that, the basal part of the leaves has fewer stomata in relation to the middle and apex portions and this morphology may cause some degree of restriction of CO2 influx (Freschi et al., 2010). Previous studies have demonstrated that there is an upregulation of PEPC1 expression (7am, 12pm, and 5pm) in the leaf apex of G. monostachia submitted to 20d of water deficit. However, it was interesting to note that there is also a diurnal modulation of decarboxylation activity of PEPCK with higher gene expression at 7am and an increase in the enzyme activity at noon, avoiding a futile cycling of organic acid during CAM (Pikart et al., 2020).
Considering that the plants were under nutritional restrictions and the base part was not in contact with N, the RNA-seq analysis showed a higher NR and NiR transcript accumulation following water deficit at 7pm in the foliar apex. In addition to the primary role in N uptake and assimilation, the enzyme nitrate reductase has also been described to possess a nitrite-NO reductase activity that uses nitrite to produce nitric oxide (NO; Astier et al., 2018). This signaling molecule could be involved with stomatal closure in response to plant acclimation to drought (Mioto et al., 2014). Detached leaves of G. monostachia exposed for 7days to polyethylene glycol showed that NO increased exclusively in the apical part of the leaves and the stomata remained closed even in the night period (CAM-idling mode; Mioto and Mercier, 2013). In Sedum album, another facultative C3-CAM plant (Crassulaceae), an increase in NR activity enzyme and a corresponding increase in endogenous NO levels were observed under drought treatment. This rise in NO level occurred during C3-CAM transition, correlating well with CAM expression (Habibi, 2020).
Photosynthesis and N metabolism are closely interconnected with the latter being a sink for ATP, reduced power and carbon skeletons produced during photosynthesis (Lea and Ireland, 1999). NH4+, whether converted from nitrate or generated by urea hydrolyzation, is assimilated via the GS/glutamine-2-oxoglutarate aminotransferase (GOGAT) cycle, which can be considered very important to maintain NH4+ at low concentrations in the leaves, helping to counteract ROS formation due to the consumption of reductants, such as Fdred and NADH (Lea and Miflin, 2003). In the context of nutritional deficiency, as the transcriptome of G. monostachia was performed, the base of the leaves presented markedly lower transcript levels of GS2, GS1, Fd or NADH-GOGAT, NADP-ICDH and urease genes. GS1 transcripts increased in the apical part in response to water stress, suggesting some importance in the G. monostachia CAM-idling shift. McNally et al. (1983), studying some C4-type and CAM plants, observed that the foliar GS1 enzyme can contribute 35–80% of total GS-activity. Previous work with epiphytic tank bromeliads, such as V. gigantea and G. monostachia, demonstrated that assimilation of NH4+ took place mainly in the apical leaf part (Takahashi and Mercier, 2011; Gonçalves et al., 2020b). The enzyme GS1 is located in the cytosol and plays an important role in the primary NH4+ assimilation from the environment (Sakakibara et al., 1996). In pine and other conifers NH4+, which is the predominant source of N absorbed for these trees, is assimilated in the cytosol through GS1, and glutamine (Gln) and glutamate (Glu) biosynthesis occur in the chloroplast (Suárez et al., 2002). Tomato GS1/GOGAT cycle was related to stress responsiveness since their expression levels were modified under drought, cold or heat stress treatments (Liu et al., 2016). On the other hand, GS2 and Fd- or NADH-GOGAT transcripts of G. monostachia did not show any increase under water deficiency, and GS2/Fd-GOGAT genes were highly expressed in the apical leaf part. It is well known that the chloroplastidic GS2/Fd-GOGAT cycle plays important roles in the assimilation of N from NO3− in leaf and in plant survival under photorespiratory conditions which lead to the generation of excess NH4+ accumulation (Krapp, 2015), while NADH-GOGAT is crucial for NH4+ assimilation in heterotrophic tissues (Selinski and Scheibe, 2019). For Mesembryanthemum crystalinum, after the transition from C3 to CAM the plants still maintain the capacity to photorespire (Whitehouse et al., 1991). In a previous study using CAM-induced leaves of G. monostachia, it was observed that glycolate oxidase 1 (GLO1) mRNA abundance, which is a gene considered a photorespiration indicator, presented higher expression at noon in CAM-induced plant compared to the dawn and dusk times (Pikart et al., 2020). The shift from C3 to CAM-idling in G. monostachia leaf apex could be related to the increase in the NADP-ICDH gene expression, which encodes the cytosolic enzyme isocitrate dehydrogenase. This enzyme allows a net synthesis of Glu, providing the key organic acid 2-oxoglutarate (2-OG), which is then transported into plastids allowing GOGAT functioning (Hodges, 2002). Since 2-OG could play a significant role modulating the flow of C to N metabolism depending upon environmental conditions (Foyer et al., 2011), it is possible that changes in its synthesis could reflect the decrease of internal CO2 levels caused by day/night stomatal closure, which is typical of the CAM-idling photosynthesis.
Urea Can Positively Modulate C3-CAM Shift Under Water Deficit
Guzmania monostachia, as an example of epiphytic tank bromeliad, is a good model to study the effects of organic nitrogen nutrition on bromeliad development and, particularly, on photosynthetic plasticity when exposed to water deficit (Rodrigues et al., 2014, 2016). In the low-nutrient epiphytic environment, which frequently presents nitrogen and water scarcity, organic compounds such as urea could represent an important nutritional resource that provides both nitrogen and carbon (Matiz et al., 2017). Our investigation reveals for the first time that urea can positively modulate the C3-CAM shift when the leaves were under drought stress. Leaves from urea treatment showed the highest increases in the expressions of PEPC1 and ALMT9 genes comparing with inorganic N sources (NO3− or NH4+) and, consequently, the greatest levels of PEPC activity and nocturnal acid accumulation. Therefore, it seems that the change from C3 to CAM mode was stimulated by this organic N source. On the other hand, plants supplied with NO3− were less vulnerable to water loss under drought conditions while showed lower expressions of PEPC1 and ALMT9 genes. Nitrate is known to increase the endogenous cytokinin level, which may regulate negatively the PEPC activity and decrease CAM in G. monostachia leaves (Pereira et al., 2013, 2018). Furthermore, NO3− uptake could be accompanied by K+ absorption in Arabidopsis thaliana as a charge balance, which can enable osmotic adjustments and may keep the RWC in this plant (Wang et al., 2012; Rubio et al., 2014). Therefore, the maintenance of water content may delay or down-regulate the CAM.
It is well known that abiotic stresses, such as drought, induce stomata closure and, consequently, limit uptake of atmospheric CO2 into the leaf tissue (Chen et al., 2018; Gobara et al., 2020). In this context, urea hydrolysis through foliar urease activity may partially compensate for the CO2 reduction by providing it, in addition to NH4+ (Matiz et al., 2017, 2019). This mechanism may represent an important physiological adaptation to obtain C and N at the same time. For G. monostachia, urea and NH4+ are the preferred N forms of uptake (Gonçalves et al., 2020b). The absorption of urea as a whole molecule took place mainly at the basal region of the leaf through the transporters of low (TIP2.1) and high (DUR 3) affinity and may be stored in the vacuole before the water deficit imposition (Matiz et al., 2019; Gonçalves et al., 2020b). If we consider the CAM-idling photosynthetic mode, whose gas exchange with the atmosphere largely ceases and respiration provides the CO2 (Winter, 2019), this strategy of using urea as a preferential source of N becomes even more relevant. Upper leaf parts of the C3 tank epiphytic bromeliad V. gigantea which close their stomata during the day under drought had a greater increase in urease activity after being supplied with urea at nighttime when compared to daytime. At the same time, a higher PEPC activity was found, indicating that at least part of the CO2 produced by urea hydrolyses was fixed during the night (Matiz et al., 2017). Furthermore, a cytochemical detection of CO2 was performed in the apex leaf portion after urea application during the night, and CO2 deposition was detected in a greater quantity in the cytoplasm near the chloroplasts than in well-watered leaves (Matiz et al., 2017).
It seems that urease activity occurs during both the day and night in G. monostachia leaves since we observed an increase in the expression of urease and GS1 genes under water scarcity. Thus, we consider that the first step of NH4+ assimilation may take place in the cytosol of this bromeliad through GS1 enzyme, giving raise to glutamine (Gln) concentration that could enter the chloroplast. Subsequently, the enzyme GOGAT transfers the amide group of Gln to 2-OG, generating two molecules of glutamate (Glu), one of which may be used to produce other amino acids through the action of transaminases (Forde and Lea, 2007). Recently, in A. thaliana, Gln biosynthesis via GS1 was found to be related to increased tolerance to abiotic stress (Ji et al., 2019). Net NH4+ assimilation requires a source of 2-OG, which is generated through partial respiration of sugars. The provision of carbon skeletons for GOGAT functioning could be supplied when the citrate is exported from the mitochondria and transformed to isocitrate and then into 2-OG by the cytosolic NADP-ICDH enzyme. A number of observations corroborate the idea that cytosolic NADP-ICDH is the major isoform in green leaves and plays an important role in N-assimilation (Foyer et al., 2011). Levels of 2-OG can reflect C/N status and may play a signaling role in the co-ordination of carbon and nitrogen metabolism (Hodges, 2002).
Previous work using G. monostachia detached leaves showed that NH4+ was more effective than nitrate to positively modulate the increase in CAM expression under water stress. The expressions of ALMT and PEPC genes were intensified, and a higher nocturnal malic acid accumulation was observed. On the other hand, NO3− supplied at various concentrations had little effect on the transcript number of ALMT and PEPC genes (Pereira et al., 2018). In the present investigation, we observed that when urea was supplied as the sole N source to the tank of G. monostachia plants, the nocturnal acid malic concentration was around five times greater than that in the presence of NH4+. After urea hydrolysis in the cytoplasm, the availability in CO2 may increase its incorporation by PEPC1, while more oxalacetate acid (OAA) and malate (through the action of NAD-MDH) can be produced, increasing acidity at night. Alternatively, during the day the CO2 generated by urea hydrolysis could be assimilated by ribulose-1,5-biphosphate carboxylase/oxygenase (Rubisco) in chloroplasts. In contrast, NO3− treatment had no effect on C3-CAM shift, curiously showing the higher RWC after 21days of water deficiency. Therefore, urea may have a positive effect on both C and N metabolism of G. monostachia, mainly if we consider a CAM-idling plant with almost no access to CO2 from the atmosphere. For some tank bromeliads, urea can contribute 30% of the total plant N (Romero et al., 2010; Gonçalves et al., 2016). In nature, this organic N source could be excreted by various animals (for instance, anurans) that visit the bromeliad’s tank searching for shelter or a place to breed, and the leaf basal part assumed the role of roots in N absorption and urea hydrolysis (Takahashi and Mercier, 2011).
Conclusion
Considering our biochemical and molecular findings (Figure 6), we conclude that drought had a positive effect on the C3-CAM shift in the apical leaf part, increasing the expression of genes related to carboxylation metabolism (PEPC1, PPCK, and NAD-MDH), movement of malate through vacuolar membrane (ALMT9) and decarboxylation process (PEPCK). Furthermore, our investigation revealed that urea can positively modulate the CAM-idling photosynthesis, which is characterized by the closure of stomata during the day and night. Our results reinforced that, when G. monostachia stay in long periods of drought, urea stimulates the expression of PEPC1 and ALMT9, while Urease transcripts increase under the same conditions. At the same time, urea hydrolysis provides NH4+, which can be assimilated into Gln since our results showed increased GS1 gene expression under water deficiency. We suggest that the link between C and N metabolism in G. monostachia may occur through the supplying of 2-OG by isocitrate dehydrogenase since the number of NADP-ICDH transcripts was also higher under drought. Thus, we suggest that urea nutrition provided in nature most by the symbiotic interaction of bromeliads with anurans may bring advantages to bromeliad survival under severe water stress, providing NH4+, which ultimately can produce amino acids and proteins, while CO2 can maintain photosynthetic efficiency, even when most stomata remain closed for 24h.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, Bioproject ID PRJNA532595.
Author Contributions
AG conceived, designed, and performed the experiments, analyzed the data, and drafted the manuscript. HM conceived and designed the experiments, analyzed the data, and drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) n° 2016/09699-5 to AG, n° 2011/50637-0 and 2018/12667-3 to HM and by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) n° 309504/2014-7 and 303497/2018-1 to HM.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We wish to thank Paulo M. R. Oliveira for technical assistance during gene expression analyses and the construction of Figure 6 of this research, Antonio A. C. M. Neto for technical assistance, and Ítalo V. C. Santos for reference list organization.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.715289/full#supplementary-material
Footnotes
References
Astier, J., Gross, I., and Durner, J. (2018). Nitric oxide production in plants: an update. J. Exp. Bot. 69, 3401–3411. doi: 10.1093/jxb/erx420
Benzing, D. H. (2008). Vascular Epiphytes: General Biology and Related Biota. Cambridge, UK: Cambridge University Press.
Benzing, D. H., and Bennett, B. (2000). Bromeliaceae: Profile of an Adaptive Radiation. Cambridge, UK: Cambridge University Press.
Camacho, C., Madden, T., Coulouris, G., Ma, N., Tao, T., Agarwala, R., et al. (2008). BLAST Command Line Applications User Manual. United States: National Center for Biotechnology Information.
Cambuí, C., and Mercier, H. (2006). Organic nitrogen nutrition in plants. Curr. Top. Plant Biol. 7, 13–17.
Chen, H., Zhang, Q., Lu, Z., and Xu, F. (2018). Accumulation of ammonium and reactive oxygen mediated drought-induced rice growth inhibition by disturbed nitrogen metabolism and photosynthesis. Plant Soil 431, 107–117. doi: 10.1007/s11104-018-3752-0
Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Forde, B. G., and Lea, P. J. (2007). Glutamate in plants: metabolism, regulation, and signalling. J. Exp. Bot. 58, 2339–2358. doi: 10.1093/jxb/erm121
Foyer, C. H., Noctor, G., and Hodges, M. (2011). Respiration and nitrogen assimilation: targeting mitochondria-associated metabolism as a means to enhance nitrogen use efficiency. J. Exp. Bot. 62, 1467–1482. doi: 10.1093/jxb/erq453
Freschi, L., Takahashi, C. A., Cambui, C. A., Semprebom, T. R., Cruz, A. B., Mioto, P. T., et al. (2010). Specific leaf areas of the tank bromeliad Guzmania monostachia perform distinct functions in response to water shortage. J. Plant Physiol. 167, 526–533. doi: 10.1016/j.jplph.2009.10.011
Gobara, B. N. K., Alves, F. R. R., Pikart, F. C., Gonçalves, A. Z., Santos, D. Y. A. C. D., Pinna, G. F. D. A. M. D., et al. (2020). How does a C3 epiphytic tank bromeliad respond to drought? Bot. J. Linn. Soc. 192, 855–867. doi: 10.1093/botlinnean/boz093
Gonçalves, A. Z., Latansio, S., Detmann, K. C., Marabesi, M. A., Neto, A. A. C., Aidar, M. P. M., et al. (2020a). What does the RuBisCO activity tell us about a C3-CAM plant? Plant Physiol. Biochem. 147, 172–180. doi: 10.1016/j.plaphy.2019.12.020
Gonçalves, A. Z., Mercier, H., Oliveira, R. S., and Romero, G. Q. (2016). Trade-off between soluble protein production and nutritional storage in bromeliaceae. Ann. Bot. 118, 1199–1208. doi: 10.1093/aob/mcw174
Gonçalves, A. Z., Oliveira, P. M. R., Neto, A. A. C., and Mercier, H. (2020b). Thinking of the leaf as a whole plant: how does N metabolism occur in a plant with foliar nutrient uptake? Environ. Exp. Bot. 178:104163. doi: 10.1016/j.envexpbot.2020.104163
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29:644. doi: 10.1038/nbt.1883
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Habibi, G. (2020). Comparison of CAM expression, photochemistry and antioxidant responses in Sedum album and Portulaca oleracea under combined str. Physiol. Plant. 170, 550–568. doi: 10.1111/ppl.13187
Hodges, M. (2002). Enzyme redundancy and the importance of 2-oxoglutarate in plant ammonium assimilation. J. Exp. Bot. 53, 905–916. doi: 10.1093/jexbot/53.370.905
Inselsbacher, E., Cambui, C. A., Richter, A., Stange, C. F., Mercier, H., and Wanek, W. (2007). Microbial activities and foliar uptake of nitrogen in the epiphytic bromeliad Vriesea gigantea. New Phytol. 175, 311–320. doi: 10.1111/j.1469-8137.2007.02098.x
Ji, Y., Li, Q., Liu, G., Selvaraj, G., Zheng, Z., Zou, J., et al. (2019). Roles of cytosolic glutamine synthetases in Arabidopsis development and stress responses. Plant Cell Physiol. 60, 657–671. doi: 10.1093/pcp/pcy235
Kleingesinds, C. K., Gobara, B. N. K., Mancilha, D., Rodrigues, M. A., Demarco, D., and Mercier, H. (2018). Impact of tank formation on distribution and cellular organization of trichomes within Guzmania monostachia rosette. Flora 243, 11–18. doi: 10.1016/j.flora.2018.03.013
Knudson, L. (1946). A new nutrient solution for the germination of orchid seed. Amer. Orchid Soc. Bull. 15, 214–217.
Kong, W., Yoo, M.-J., Zhu, D., Noble, J. D., Kelley, T. M., Li, J., et al. (2020). Molecular changes in mesembryanthemum crystallinum guard cells underlying the C3 to CAM transition. Plant Mol. Biol. 103, 653–667. doi: 10.1007/s11103-020-01016-9
Krapp, A. (2015). Plant nitrogen assimilation and its regulation: a complex puzzle with missing pieces. Curr. Opin. Plant Biol. 25, 115–122. doi: 10.1016/j.pbi.2015.05.010
Lea, P. J., and Ireland, R. J. (1999). “Nitrogen metabolism in higher plants.” in Plant Amino Acids: Biochemistry and Biotechnology. ed. Singh, B. K. (Boca Raton, FL: CRC Press), 1–48.
Lea, P. J., and Miflin, B. J. (2003). Glutamate synthase and the synthesis of glutamate in plants. Plant Physiol. Biochem. 41, 555–564. doi: 10.1016/S0981-9428(03)00060-3
Liu, L., Wang, J., Han, Z., Sun, X., Li, H., Zhang, J., et al. (2016). Molecular analyses of tomato GS, GOGAT and GDH gene families and their response to abiotic stresses. Acta Physiol. Plant. 38:229. doi: 10.1007/s11738-016-2251-2
Luttge, U. (2004). Ecophysiology of crassulacean acid metabolism (CAM). Ann. Bot. 93, 629–652. doi: 10.1093/aob/mch087
Lüttge, U. (2006). Photosynthetic flexibility and ecophysiological plasticity: questions and lessons from Clusia, the only CAM tree, in the neotropics. New Phytol. 171, 7–25. doi: 10.1111/j.1469-8137.2006.01755.x
Matiz, A., Mioto, P. T., Aidar, M. P. M., and Mercier, H. (2017). Utilization of urea by leaves of bromeliad Vriesea gigantea under water deficit: much more than a nitrogen source. Biol. Plant. 61, 751–762. doi: 10.1007/s10535-017-0721-z
Matiz, A., Mioto, P. T., Mayorga, A. Y., Freschi, L., and Mercier, H. (2013). “CAM photosynthesis in bromeliads and agaves: what can we learn from these plants.” in Photosynthesis. ed. Dubinsky, Z. (Rijeka, Croatia: Intech), 91–134.
Matiz, A., Mioto, P. T., and Mercier, H. (2019). “Urea in plants: metabolic aspects and ecological implications.” in Progress in Botany. Vol. 81, eds. F. L. Cánovas, U. Leuschner, and M. C. Risueño (Cham: Springer International Publishing), 157–187.
Maxwell, C., Griffiths, H., Borland, A., Young, A., Broadmeadow, M., and Fordham, M. (1995). Short-term photosynthetic responses of the C3-CAM epiphyte Guzmania monostachia var. monostachia to tropical seasonal transitions under field conditions. Funct. Plant Biol. 22, 771–781. doi: 10.1071/PP9950771
McNally, S. F., Hirel, B., Gadal, P., Mann, A. F., and Stewart, G. R. (1983). Glutamine synthetases of higher plants. Plant Physiol. 72, 22–25. doi: 10.1104/pp.72.1.22
Medina, E., Delgado, M., Troughton, J. H., and Medina, J. D. (1977). Physiological ecology of CO2 fixation in bromeliaceae. Flora 166, 137–152. doi: 10.1016/S0367-2530(17)32126-6
Mercier, H., Rodrigues, M. A., Andrade, S. C. D. S., Coutinho, L. L., Gobara, B. N. K., Matiz, A., et al. (2019). Transcriptional foliar profile of the C3-CAM bromeliad Guzmania monostachia. PLoS One 14:e0224429. doi: 10.1371/journal.pone.0224429
Miflin, B. J., and Lea, P. J. (1980). “Ammonia assimilation.” in The Biochemistry of Plants: A Comprehensive Treatise. eds. P. K. Stumpf, E. E. Conn, N. E. Tolbert, D. D. Davis, J. Preiss, and B. J. Miflin et al. (New York: Academic Press), 169–202.
Mioto, P. T., Freschi, L., and Mercier, H. (2014). “Phytohormones and nitric oxide interactions during abiotic stress responses.” in Nitric Oxide in Plants: Metabolism and Role in Stress Physiology. eds. Khan, M. N., Mobin, M., Mohammad, F., and Corpas, F. J. (Cham: Springer), 211–224.
Mioto, P. T., and Mercier, H. (2013). Abscisic acid and nitric oxide signaling in two different portions of detached leaves of Guzmania monostachia with CAM up-regulated by drought. J. Plant Physiol. 170, 996–1002. doi: 10.1016/j.jplph.2013.02.004
Murashige, T., and Skoog, F. (1962). A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant. 15, 473–497. doi: 10.1111/j.1399-3054.1962.tb08052.x
Nievola, C. C., Mercier, H., and Majerowicz, N. (2001). Levels of nitrogen assimilation in bromeliads with different growth habits. J. Plant Nutr. 24, 1387–1398. doi: 10.1081/PLN-100106989
Pereira, P. N., Gaspar, M., Smith, J. A. C., and Mercier, H. (2018). Ammonium intensifies CAM photosynthesis and counteracts drought effects by increasing malate transport and antioxidant capacity in Guzmania monostachia. J. Exp. Bot. 69, 1993–2003. doi: 10.1093/jxb/ery054
Pereira, P. N., Purgatto, E., and Mercier, H. (2013). Spatial division of phosphoenolpyruvate carboxylase and nitrate reductase activity and its regulation by cytokinins in CAM-induced leaves of Guzmania monostachia (Bromeliaceae). J. Plant Physiol. 170, 1067–1074. doi: 10.1016/j.jplph.2013.03.005
Pfaffl, M. W., Tichopad, A., Prgomet, C., and Neuvians, T. P. (2004). Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper – excel-based tool using pair-wise correlations. Biotechnol. Lett. 26, 509–515. doi: 10.1023/B:BILE.0000019559.84305.47
Pikart, F. C., Marabesi, M. A., Mioto, P. T., Gonçalves, A. Z., Matiz, A., Alves, F. R. R., et al. (2018). The contribution of weak CAM to the photosynthetic metabolic activities of a bromeliad species under water deficit. Plant Physiol. Biochem. 123, 297–303. doi: 10.1016/j.plaphy.2017.12.030
Pikart, F. C., Matiz, A., Alves, F. R. R., and Mercier, H. (2020). Diurnal modulation of PEPCK decarboxylation activity impacts photosystem II light-energy use in a drought-induced CAM species. Environ. Exp. Bot. 173:104003. doi: 10.1016/j.envexpbot.2020.104003
R Core Team (2020). R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.
Rodrigues, M. A., Freschi, L., Pereira, P. N., and Mercier, H. (2014). “Interactions between nutrients and crassulacean acid metabolism.” in Progress in Botany. eds. Lüttge, U., Beyschlag, W., and Cushman, J. (Berlin Heidelberg: Springer-Verlag), 167–186.
Rodrigues, M. A., Hamachi, L., Mioto, P. T., Purgatto, E., and Mercier, H. (2016). Implications of leaf ontogeny on drought-induced gradients of CAM expression and ABA levels in rosettes of the epiphytic tank bromeliad Guzmania monostachia. Plant Physiol. Biochem. 108, 400–411. doi: 10.1016/j.plaphy.2016.08.010
Romero, Q. G., Nomura, F., Gonçalves, A. Z., Dias, N. Y. N., Mercier, H., Conforto, E. D. C., et al. (2010). Nitrogen fluxes from treefrogs to tank epiphytic bromeliads: an isotopic and physiological approach. Oecologia 162, 941–949. doi: 10.1007/s00442-009-1533-4
Rubio, F., Fon, M., Ródenas, R., Nieves, C. M., Alemám, F., Rivero, R. M., et al. (2014). A low K+ signal is required for functional high affinity K+ uptake through HAK5 transporters. Physiol. Plant. 152, 558–570. doi: 10.1111/ppl.12205
Sakakibara, H., Shimizu, H., Hase, T., Yamazaki, Y., Takao, T., Shimonishi, Y., et al. (1996). Molecular identification and characterization of cytosolic isoforms of glutamine synthetase in maize roots. J. Biol. Chem. 271, 29561–29568. doi: 10.1074/jbc.271.47.29561
Selinski, J., and Scheibe, R. (2019). Malate valves: old shuttles with new perspectives. Plant Biol. 21:(S1), 21–30. doi: 10.1111/plb.12869
Suárez, M. F., Avila, C., Gallardo, F., Cantón, F. R., García-Gutiérrez, A., Claros, M. G., et al. (2002). Molecular and enzymatic analysis of ammonium assimilation in woody plants. J. Exp. Bot. 53, 891–904. doi: 10.1093/jexbot/53.370.891
Takahashi, C. A., and Mercier, H. (2011). Nitrogen metabolism in leaves of a tank epiphytic bromeliad: characterization of a spatial and functional division. J. Plant Physiol. 168, 1208–1216. doi: 10.1016/j.jplph.2011.01.008
Wang, Y. Y., Hsu, P. K., and Tsay, Y. F. (2012). Uptake, allocation and signaling of nitrate. Trends Plant Sci. 17, 458–467. doi: 10.1016/j.tplants.2012.04.006
Weatherley, P. E. (1950). Studies in the water relations of the cotton plant. I. The field measurement of water deficits in leaves. New Phytol. 49, 81–97. doi: 10.1111/j.1469-8137.1950.tb05146.x
Whitehouse, D. G., Rogerrs, W. J., and Tobin, A. K. (1991). Photorespiratory enzyme activities in C3 and CAM forms of the facultative CAM plant, mesembryanthemum crystallimum L. J. Exp. Bot. 42, 485–492. doi: 10.1093/jxb/42.4.485
Winter, K. (2019). Ecophysiology of constitutive and facultative CAM photosynthesis. J. Exp. Bot. 70, 6495–6508. doi: 10.1093/jxb/erz002
Zotz, G. (2013). The systematic distribution of vascular epiphytes: a critical update. Bot. J. Linn. Soc. 171, 453–481. doi: 10.1111/boj.12010
Keywords: CAM, ICDH, malate, N nutrition, PEPC, PEPCK, PPCK, urea
Citation: Gonçalves AZ and Mercier H (2021) Transcriptomic and Biochemical Analysis Reveal Integrative Pathways Between Carbon and Nitrogen Metabolism in Guzmania monostachia (Bromeliaceae) Under Drought. Front. Plant Sci. 12:715289. doi: 10.3389/fpls.2021.715289
Edited by:
Milton Lima Neto, São Paulo State University, BrazilReviewed by:
Ana Karla M. Lobo, São Paulo State University, BrazilJohan Ceusters, KU Leuven, Belgium
Copyright © 2021 Gonçalves and Mercier. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana Zangirolame Gonçalves, zangirolame@gmail.com