Arabidopsis Response to Inhibitor of Cytokinin Degradation INCYDE: Modulations of Cytokinin Signaling and Plant Proteome

Cytokinins are multifaceted plant hormones that play crucial roles in plant interactions with the environment. Modulations in cytokinin metabolism and signaling have been successfully used for elevating plant tolerance to biotic and abiotic stressors. Here, we analyzed Arabidopsis thaliana response to INhibitor of CYtokinin DEgradation (INCYDE), a potent inhibitor of cytokinin dehydrogenase. We found that at low nanomolar concentration, the effect of INCYCE on seedling growth and development was not significantly different from that of trans-Zeatin treatment. However, an alteration in the spatial distribution of cytokinin signaling was found at low micromolar concentrations, and proteomics analysis revealed a significant impact on the molecular level. An in-depth proteome analysis of an early (24 h) response and a dose-dependent response after 168 h highlighted the effects on primary and secondary metabolism, including alterations in ribosomal subunits, RNA metabolism, modulations of proteins associated with chromatin, and the flavonoid and phenylpropanoid biosynthetic pathway. The observed attenuation in stress-response mechanisms, including abscisic acid signaling and the metabolism of jasmonates, could explain previously reported positive effects of INCYDE under mild stress conditions.


Introduction
Plants are sessile organisms that, in order to maintain growth and survival under harsh conditions, have evolved unique mechanisms, enabling them to react rapidly to ever-changing ambient conditions. A fundamental part of this mechanism is the plant hormone cytokinin, which is known to function in

Root Growth in Response to tZ and INCYDE Treatment was not Significantly Different
First, to compare the effects of exogenously applied tZ and INCYDE on Arabidopsis physiology, we monitored the primary root growth. Root growth inhibition is a well-known response to cytokinin (e.g., [27]). Seedlings were cultivated as described in Materials and Methods, and after seven days, the length of primary roots was evaluated (Figure 1a). The inhibition of root growth was already observed at the lowest applied concentration (10 nM), and both substances fully retarded seedling growth at 10 µM, which was also accompanied by cotyledon coloring due to the accumulation of anthocyanin (a cytokinin response described in previous studies, e.g., [28]). The INCYDE treatment had a slightly lower effect at 10 nM, but the saturating concentration was lower than that of tZ ( Figure 1b).

INCYDE Treatment Elicited Distinct Spatial Distribution of Cytokinin Signaling
Next, the effect of INCYDE and exogenously supplied tZ on cytokinin signaling was compared. Transgenic Arabidopsis lines bearing ARR5::GUS reporter were grown within the concentration range of INCYDE, and after seven days, the ARR5 promoter activity was visualized by histochemical staining and compared to that of tZ (Figure 2a,b). The observed differences between tZ and INCYDE were Plants 2020, 9,1563 3 of 15 concentration-dependent. The analysis revealed that the growth of seedlings in the presence of INCYDE led to a higher increase in the ARR5 promoter activity in cotyledons (100-500 nM) but a lower cytokinin signaling in the roots (10-100 nM). The most promising response was found for seedlings grown in the presence of 500 nM tZ and INCYDE, with a similar increase in the ARR5 promoter activity in the roots, but a strikingly different pattern in cotyledons ( Figure 2). A similar effect was observed in 14-day-old plants ( Figure S1).

INCYDE Treatment Elicited Distinct Spatial Distribution of Cytokinin Signaling
Next, the effect of INCYDE and exogenously supplied tZ on cytokinin signaling was compared. Transgenic Arabidopsis lines bearing ARR5::GUS reporter were grown within the concentration range of INCYDE, and after seven days, the ARR5 promoter activity was visualized by histochemical staining and compared to that of tZ (Figure 2a,b). The observed differences between tZ and INCYDE were concentration-dependent. The analysis revealed that the growth of seedlings in the presence of INCYDE led to a higher increase in the ARR5 promoter activity in cotyledons (100-500 nM) but a lower cytokinin signaling in the roots (10-100 nM). The most promising response was found for seedlings grown in the presence of 500 nM tZ and INCYDE, with a similar increase in the ARR5 promoter activity in the roots, but a strikingly different pattern in cotyledons ( Figure 2). A similar effect was observed in 14-day-old plants ( Figure S1).

INCYDE Treatment Elicited Distinct Spatial Distribution of Cytokinin Signaling
Next, the effect of INCYDE and exogenously supplied tZ on cytokinin signaling was compared. Transgenic Arabidopsis lines bearing ARR5::GUS reporter were grown within the concentration range of INCYDE, and after seven days, the ARR5 promoter activity was visualized by histochemical staining and compared to that of tZ (Figure 2a,b). The observed differences between tZ and INCYDE were concentration-dependent. The analysis revealed that the growth of seedlings in the presence of INCYDE led to a higher increase in the ARR5 promoter activity in cotyledons (100-500 nM) but a lower cytokinin signaling in the roots (10-100 nM). The most promising response was found for seedlings grown in the presence of 500 nM tZ and INCYDE, with a similar increase in the ARR5 promoter activity in the roots, but a strikingly different pattern in cotyledons ( Figure 2). A similar effect was observed in 14-day-old plants ( Figure S1).   Table S1 for details).

Early INCYDE Response Proteins of Arabidopsis Seedlings Highlight Similarity to Exogenous tZ Treatment
The dose-dependent growth response to INCYDE indicated that the most interesting comparison between tZ and INCYDE treatment was within the 100-1000 nM range ( Figure 1). To provide an insight into the molecular mechanisms behind the observed contrasting response, seven-day-old seedlings were treated for 24 h with 500 nM tZ, INCYDE, or 0.1% (v/v) dimethylsulphoxide (DMSO, mock), as described in Materials and Methods, to determine an early response. In total, 3273 Arabidopsis proteins were identified with reliable quantitative data for more than 2100 of these. Higher biological variability for tZ was observed within the set of four biological replicates (each pooled from at least 20 seedlings), but statistically significant (p < 0.05) separation between tZ, INCYDE and mock-treated samples was apparent (Figure 3a). A detailed analysis revealed 89 and 99 tZ and INCYDE early-response proteins as compared to mock-treatment, respectively, and 69 proteins that showed statistically significant and reproducible differences between INCYDE-and tZ-treated samples (Figure 3b). trans-Zeatin (tZ, blue; n = 16), INCYDE (green; n=20) or DMSO (mock, white; n = 16). Results represent means and standard error, different letters indicate significant differences (Kruskal-Wallis, p < 0.05; see Table S1 for details).

Early INCYDE Response Proteins of Arabidopsis Seedlings Highlight Similarity to Exogenous tZ Treatment
The dose-dependent growth response to INCYDE indicated that the most interesting comparison between tZ and INCYDE treatment was within the 100-1000 nM range ( Figure 1). To provide an insight into the molecular mechanisms behind the observed contrasting response, sevenday-old seedlings were treated for 24 h with 500 nM tZ, INCYDE, or 0.1% (v/v) dimethylsulphoxide (DMSO, mock), as described in Materials and Methods, to determine an early response. In total, 3273 Arabidopsis proteins were identified with reliable quantitative data for more than 2100 of these. Higher biological variability for tZ was observed within the set of four biological replicates (each pooled from at least 20 seedlings), but statistically significant (p < 0.05) separation between tZ, INCYDE and mock-treated samples was apparent (Figure 3a). A detailed analysis revealed 89 and 99 tZ and INCYDE early-response proteins as compared to mock-treatment, respectively, and 69 proteins that showed statistically significant and reproducible differences between INCYDE-and tZtreated samples (Figure 3b).  As illustrated in Figure 4, early INCYDE response proteins were functionally enriched in amino acid and carbohydrate metabolism, but the set encompassed diverse processes of both primary and secondary metabolism. The comparison with tZ-responsive proteins showed only 48 shared proteins, but all with a similar responsiveness. Repressed proteins included an enzyme of anthocyanin metabolism leucoanthocyanidin dioxygenase (AT4G22880), senescence-associated protein SAG24 (AT1G66580) and a protein reportedly involved in the regulation of gravitropic response and auxin transport in roots ROSY1 (AT2G16005, [29]). Significantly accumulated proteins included those required for chloroplast biogenesis and development (chorismate mutase 3, AT3G29200; CPP1, AT5G23040, [30]), enzymes involved in cell wall formation (AT5G15490, AT4G37800; [31]), As illustrated in Figure 4, early INCYDE response proteins were functionally enriched in amino acid and carbohydrate metabolism, but the set encompassed diverse processes of both primary and secondary metabolism. The comparison with tZ-responsive proteins showed only 48 shared proteins, but all with a similar responsiveness. Repressed proteins included an enzyme of anthocyanin metabolism leucoanthocyanidin dioxygenase (AT4G22880), senescence-associated protein SAG24 (AT1G66580) and a protein reportedly involved in the regulation of gravitropic response and auxin transport in roots ROSY1 (AT2G16005, [29]). Significantly accumulated proteins included those required for chloroplast biogenesis and development (chorismate mutase 3, AT3G29200; CPP1, AT5G23040, [30]), enzymes involved in cell wall formation (AT5G15490, AT4G37800; [31]), chloroplastic and cytosolic isoform of glutamine synthetase (AT5G35630, AT1G66200), and extracellular protein with a putative role in circadian rhythm and abiotic stress GER3 (AT5G20630). endocytosis (RABG3c, AT3G16100; AGD8, AT4G17890).

Correlation between INCYDE Concentration and Protein Abundance
Next, an INCYDE dose-response was evaluated at a proteome-wide scale. Seven-day-old seedlings cultivated on textile meshes were transferred onto new medium supplemented with 10-1000 nM INCYDE or 0.1% (v/v) DMSO (mock), as described in Materials and Methods. After seven days, shoots were collected for proteome and metabolome analyses. The effect of INCYDE on the growth of fully established seedling was significantly lower than that observed at the early developmental stage (Supplementary Figure S2), but its effect was well-manifested at the proteome level. In total, 29,530 peptide groups were identified, providing sufficient quantitative data for more than 2400 proteins. Independent component analysis ( Figure 5a) showed significant separation of INCYDE-and mock-treated samples, and indicated a high level of similarity between 100 and 1000 nM INCYDE treatments. A detailed comparison of differentially abundant proteins (absolute foldchange 1.4, p < 0.05) confirmed the expected overlap between the INCYDE treatments, and revealed that the highest response was elicited at an INCYDE concentration of 100 nM. The correlation analysis identified 25 and 11 statistically significant positive and negative correlations, respectively (absolute Pearson's correlation coefficients >0.7, p < 0.01). These INCYDE-responsive proteins are involved in Figure 4. Visualization of processes mediated by early INCYDE response proteins. Enrichment analysis based on hypergeometric distribution followed by FDR correction. Two pathways are connected if they share 20% or more proteins. Darker and bigger nodes are more significantly enriched and larger sets, respectively. Analyzed by ShinyGO 0.61 [32].

Correlation between INCYDE Concentration and Protein Abundance
Next, an INCYDE dose-response was evaluated at a proteome-wide scale. Seven-day-old seedlings cultivated on textile meshes were transferred onto new medium supplemented with 10-1000 nM INCYDE or 0.1% (v/v) DMSO (mock), as described in Materials and Methods. After seven days, shoots were collected for proteome and metabolome analyses. The effect of INCYDE on the growth of fully established seedling was significantly lower than that observed at the early developmental stage (Supplementary Figure S2), but its effect was well-manifested at the proteome level. In total, 29,530 peptide groups were identified, providing sufficient quantitative data for more than 2400 proteins. Independent component analysis (Figure 5a) showed significant separation of INCYDEand mock-treated samples, and indicated a high level of similarity between 100 and 1000 nM INCYDE treatments. A detailed comparison of differentially abundant proteins (absolute fold-change 1.4, p < 0.05) confirmed the expected overlap between the INCYDE treatments, and revealed that the highest response was elicited at an INCYDE concentration of 100 nM. The correlation analysis identified 25 and 11 statistically significant positive and negative correlations, respectively (absolute Pearson's correlation coefficients >0.7, p < 0.01). These INCYDE-responsive proteins are involved in signaling, primary and secondary metabolism, including γ-aminobutyric acid biosynthesis, glucosinolate degradation and processes of chloroplast biogenesis and development (Table 1).

INCYDE Response Proteins Modulate Diverse Metabolic Processes
In total, 517 proteins representing an estimated 29% of the total protein extract were found to be INCYDE-responsive (Figure 5b,c). A functional analysis revealed 199 enriched biological processes (Gene Ontology; GO), covering 448 INCYDE-responsive proteins. The strongest response was detected for 100 nM INCYDE treatment, and the GO enrichment ( Figure 6) revealed that this response elicited the highest similarity to annotated response to cytokinin. Other highly enriched categories compared to 10 and 1000 nM treatments included 'response to oxidative stress', 'photosynthesis', 'response to light stimulus' and 'ribosome biogenesis'.

INCYDE Response Proteins Modulate Diverse Metabolic Processes
In total, 517 proteins representing an estimated 29% of the total protein extract were found to be INCYDE-responsive (Figure 5b,c). A functional analysis revealed 199 enriched biological processes (Gene Ontology; GO), covering 448 INCYDE-responsive proteins. The strongest response was detected for 100 nM INCYDE treatment, and the GO enrichment ( Figure 6) revealed that this response elicited the highest similarity to annotated response to cytokinin. Other highly enriched categories compared to 10 and 1000 nM treatments included 'response to oxidative stress', 'photosynthesis', 'response to light stimulus' and 'ribosome biogenesis'. There were only 40 INCYDE-responsive proteins found exclusively for the 10 nM treatment. Besides proteins associated with proteosynthesis (7), photosynthesis (5) and ubiquitin-mediated degradation (6), two INCYDE-responsive proteins are reportedly involved in biotic stress response, namely inhibitor of fungal polygalacturonase PGIP1 (AT5G06860; accumlated) and extracellular lectin, whose expression is induced upon treatment with chitin oligomers (AT3G16530; repressed), and at least three repressed proteins participate in intracellular trafficking (AT5G16880, AT5G05010, AT5G16880). The lowest number of dose-specific proteins was found for 1000 nM treatment ( Figure   Figure 6. Comparison of enriched biological processes in INCYDE-responsive proteomes. For the sake of clarity, categories with at least 15 proteins were included in the final analysis, and only GO terms with the most significant contribution to the separation in PC1 or PC2 are labeled. The separation of individual treatments (represented by boxes) indicates a similarity between 10 and 1000 nM samples.
Next, the proteins found with a similar INCYDE response in all three treatments were analyzed. The dose-independent changes in the amounts of at least ten proteins involved in RNA metabolism, and the accumulation of three proteins involved in chromatin remodeling (H2A. involved in proteosynthesis, ribosome biogenesis, tRNA metabolism and protein folding (Figure 7). The accumulation of phenylpropanoid and flavonoid biosynthetic enzymes was recently reported in barley root proteome response to tZ [33]. Here, phenylalanine ammonia-lyase 1, flavonol synthase and anthocyanidin 3-O-glucosyltransferase were significantly accumulated in response to INCYDE. Cytokinin signaling is intertwined with other phytohormonal networks, and especially that of auxin. It is thus not surprising that the INCYDE treatment resulted in the accumulation of auxin biosynthetic enzyme (AMI1, AT1G08980) and a putative negative regulator of PIN auxin transport (APM1, AT4G33090). All tested concentrations of INCYDE also elicited a depletion of abscisic acid (ABA) signaling transcription factor NFYC4 (AT5G63470), and accumulations of ABA signaling repressor phosphatase 2C (ABI1, AT4G26080) and protein CAR8 (AT1G23140). CAR8 is a putative mediator of ABA receptor interaction with the plasma membrane and could regulate sensitivity to ABA [34]. Several proteins were also connected to photomorphogenesis and plastid biogenesis, including a subunit of COP9 signalosome complex (AT4G14110, decreased), PTAC12 (AT2G34640, accumulated) which is reportedly involved in the initiation of photomorphogenesis [35], epimerase participating in plastid division AT2G21280 and protein CLMP1 (AT1G62390, accumulated), that is required for plastid separation and partitioning during cell division [36]. Last but not least, farnesyl pyrophosphate synthase 1 (AT5G47770) was significantly repressed, which may coincide with the expected repression in the biosynthesis of isopentenyl pyrophosphate, a substrate of cytokinin biosynthetic enzyme isopentenyl transferase.

Cytokinin Dehydrogenase Isoforms Play a Key Role in the Contrasting Response between tZ and INCYDE
The effect of exogenously applied cytokinin strongly depends on its transport and metabolism. It is thus not surprising that the direct tZ effect on the cytokinin signaling was predominantly  [37]; Color coding of proteins is denoted by functional designation given by GO and KEGG pathway enrichments; only the nine most relevant categories are highlighted. The line thickness indicates the strength of data support, the minimum required interaction score is 0.4 (medium confidence).

Cytokinin Dehydrogenase Isoforms Play a Key Role in the Contrasting Response between tZ and INCYDE
The effect of exogenously applied cytokinin strongly depends on its transport and metabolism. It is thus not surprising that the direct tZ effect on the cytokinin signaling was predominantly constrained to the root tissue ( Figure 2). In contrast, the mechanisms of rapid cytokinin inactivation and degradation did not affect INCYDE transport, and the INCYDE induced cytokinin response was determined by the presence of CKX enzymes. Components of cytokinin metabolism and signaling are low abundant proteins, and the quantitative data for CKX were not available in any of the INCYDE response datasets. However, the recently published detailed characterization of Arabidopsis proteome [38] implied that the amounts of cytokinin receptor in root, hypocotyl and cotyledon are comparable, and showed that the observed activation of ARR5::GUS promoter (Figure 2) was similar to the expected profiles of apoplastic isoforms CKX4 and CKX5 (Figure 8). The apoplastic CKX enzymes found in Arabidopsis act mainly on cytokinin free bases and ribosides [16], and these are also the root-to-shoot long-distance signaling forms of cytokinin [39]. It is thus possible that at a given developmental stage, the majority of the INCYDE activity comprises the transport of active cytokinin in plants.
Plants 2020, 9, x FOR PEER REVIEW 10 of 16 2) was similar to the expected profiles of apoplastic isoforms CKX4 and CKX5 (Figure 8). The apoplastic CKX enzymes found in Arabidopsis act mainly on cytokinin free bases and ribosides [16], and these are also the root-to-shoot long-distance signaling forms of cytokinin [39]. It is thus possible that at a given developmental stage, the majority of the INCYDE activity comprises the transport of active cytokinin in plants.

Cytokinin Response Targets Ribosomal Proteins
Plant ribosomes are highly complex structures, with each ribosomal protein being encoded by two to seven paralogs. The ribosome composition reflects external stimuli, and may have a significant impact on plant responses [41]. Previous proteomics reports have indicated that ribosome could be the target of cytokinin signaling [33,42,43]. Here, 122 and 167 ribosomal proteins were identified in seedling and shoot proteome, respectively. First, 60S ribosomal proteins L4-1 and L4-2, previously implicated in response to cytokinin [43], were quantified in both experiments, but there was not any detectable difference compared to the mock-treated plants (Figure 9, L4e). The abundances of only three ribosomal proteins were significantly altered in early INCYDE response. Interestingly, a decrease in protein abundance was observed for 60S ribosomal protein L10-3 involved in UV-B stress response [44], and a previous analysis found expression profile similarities between UV-B receptor and components of cytokinin signaling [2]. The experiment with the dose-dependent response to INCYDE revealed a much higher response in ribosomal proteins. In total, 16 and 18 ribosomal

Cytokinin Response Targets Ribosomal Proteins
Plant ribosomes are highly complex structures, with each ribosomal protein being encoded by two to seven paralogs. The ribosome composition reflects external stimuli, and may have a significant impact on plant responses [41]. Previous proteomics reports have indicated that ribosome could be the target of cytokinin signaling [33,42,43]. Here, 122 and 167 ribosomal proteins were identified in seedling and shoot proteome, respectively. First, 60S ribosomal proteins L4-1 and L4-2, previously implicated in response to cytokinin [43], were quantified in both experiments, but there was not any detectable difference compared to the mock-treated plants (Figure 9, L4e). The abundances of only three ribosomal proteins were significantly altered in early INCYDE response. Interestingly, a decrease in protein abundance was observed for 60S ribosomal protein L10-3 involved in UV-B stress response [44], and a previous analysis found expression profile similarities between UV-B receptor and components of cytokinin signaling [2]. The experiment with the dose-dependent response to INCYDE revealed a much higher response in ribosomal proteins. In total, 16 and 18 ribosomal proteins were accumulated and decreased, respectively. Based on the KEGG (Kyoto Encyclopedia of Genes and Genomes) annotations, these proteins represent 29 ribosomal subunits, but not all of these are major paralogs; the most prominent paralogs form only 50% of the detected isoforms. This indicates that the total ribosome population is not completely altered, or that the alteration is cell-specific, and that this localized distribution was lost in the total shoot protein extracts used in this study. It is noteworthy that the ribosomal subunits did not show mixed-mode response, and that all INCYDE-responsive paralogs for the given subunit always showed the same response.

Positive Effect of INCYDE Could Coincide with Attenuated Stress Response
The response to an exogenous stimuli has been seen to change over time, and it is, thus, not surprising that an early response may be only transient and could significantly differ to that found after prolonged exposure. In the experiment reported here, only 13 early INCYDE response proteins were found to be differentially abundant after 168 h of INCYDE treatment, and only six showed a similar response, including 12-oxo-phytodienoic acid biosynthetic enzyme (AOC2, AT3G25770, repressed) and microsomal prostaglandin E synthase 2 (AT5G42150, repressed). These two enzymes produce substances that trigger plant defense and detoxification response. It has been hypothesized that plants reduce their growth as a primary adaptation response to stress, regardless of the severity of the threat [46], and it is tempting to speculate that the reported INCYDE-promoted growth under suboptimal conditions [23][24][25][26] could coincide with the attenuated stress response. This theory could be supported by an INCYDE-induced attenuation of ABA signaling and repression of at least 21 additional stress-responsive proteins, including UDP-glycosyltransferase 79B2 (AT4G27560, response to cold, drought and salinity; [47]), protein disulfide-isomerase (AT2G32920, response to  [45]. Blue and red gradients indicate subunits with one or more INCYDE-responsive paralogs. The contribution of the paralog to the subunit composition [%] is indicated; white-not detected; grey-no significant difference compared to the mock-treated plants.

Positive Effect of INCYDE Could Coincide with Attenuated Stress Response
The response to an exogenous stimuli has been seen to change over time, and it is, thus, not surprising that an early response may be only transient and could significantly differ to that found after prolonged exposure. In the experiment reported here, only 13 early INCYDE response proteins were found to be differentially abundant after 168 h of INCYDE treatment, and only six showed a similar response, including 12-oxo-phytodienoic acid biosynthetic enzyme (AOC2, AT3G25770, repressed) and microsomal prostaglandin E synthase 2 (AT5G42150, repressed). These two enzymes produce substances that trigger plant defense and detoxification response. It has been hypothesized that plants reduce their growth as a primary adaptation response to stress, regardless of the severity of the threat [46], and it is tempting to speculate that the reported INCYDE-promoted growth under suboptimal conditions [23][24][25][26] could coincide with the attenuated stress response. This theory could be supported by an INCYDE-induced attenuation of ABA signaling and repression of at least 21 additional stress-responsive proteins, including UDP-glycosyltransferase 79B2 (AT4G27560, response to cold, drought and salinity; [47]), protein disulfide-isomerase (AT2G32920, response to endoplasmic reticulum stress; [48]) or MD-2-related lipid-recognition protein 3 (AT5G23820, biotic stress, [49]). The comparison of INCYDE response proteins with the database of previously identified phytohormone-responsive proteins [50] found the highest overlap for cytokinin-responsive proteins (representing 18 and 79 proteins for 24 and 168 h treatments, respectively), but the proteins responsive to stress-associated phytohormones ABA and jasmonic acid represented the second and the third most numerous categories, with 82 ABA-responsive proteins and 55 jasmonic acid-responsive proteins, respectively. The INCYDE response was opposite to that of jasmonic acid or ABA treatment for more than 50% of these proteins, supporting the theory that the INCYDE-repressed degradation of cytokinin inhibits stress perception in plants. Further, this fact indicates a higher level of cytokinin interaction with ABA and jasmonic acid than previously reported.

Plant Material
Seeds of Arabidopsis thaliana ecotype Col-0 and transgenic promoter line ARR5::GUS were used for all experiments. Seeds were surface-sterilized and sown on Petri dishes containing half-strength Murashige and Skoog medium with 1% (w/v) agar with the corresponding concentration of INCYDE or tZ (10 nM, 100 nM, 1 µM, 10 µM) or 0.1% dimethylsulfoxide (DMSO, mock). Seeds were stratified at 4 • C for three days, and cultivated at 21 • C/19 • C day/night temperatures, with a 16 h photoperiod (100 µmol m -2 s -1 photosynthetic photon flux density) for up to 14 days in a growth chamber (AR36LX, Percival). For proteomics experiments, seeds were sown on Uhelon 120T (Silk & Progress, Brněnec, Czech Republic) and cultivated as described for mock. After seven days of cultivation, the Uhelon mesh with the seedlings was transferred onto fresh liquid Murashige and Skoog medium supplemented with 0.1% (v/v) DMSO (mock), tZ or INCYDE in DMSO (final concentration, as for the mock) for five minutes and then transferred onto solidified Murashige and Skoog medium supplemented with the same substance. Plant tissues were harvested after 24 h (whole seedling) and 168 h (shoot tissue) for the analysis of early and prolonged INCYDE response proteins, respectively. The length of the primary root of seven-day-old seedlings was measured using the analysis software ImageJ.

Proteome Analysis
Total protein extracts were prepared as previously described, employing a combination of phenol/acetone/TCA extraction [14]. Portions of samples corresponding to 5 µg of peptide were analyzed by nanoflow reverse-phase liquid chromatography-mass spectrometry using a 15 cm C18 Zorbax column (Agilent), a Dionex Ultimate 3000 RSLC nano-UPLC system and the Orbitrap Fusion Lumos Tribrid Mass Spectrometer (Thermo). Peptides were eluted with up to a 120-min, 4% to 40% acetonitrile gradient. Spectra were acquired using the default settings for peptide identification, employing HCD activation, resolution 60,000 (MS) and 15,000 (MS2), and 60 s dynamic exclusion. The measured spectra were recalibrated and searched against Araport 11 protein database, as described previously [14]. Only proteins with at least two unique peptides were considered for the quantitative analysis. The quantitative differences were determined by Minora, employing precursor ion quantification followed by normalization and background-based t-test for peptide-and protein-based quantitation.

Statistical Analyses
The reported statistical tests were generated and implemented using Instant Clue [53], Rapid Miner (www.rapidminer.com; [54]) and Proteome Discoverer. Significant differences refer to p < 0.05.

Conclusions
Cytokinin metabolism and signaling play important roles in abiotic stress tolerance, and the manipulation of these processes by inhibiting cytokinin degradation could be beneficial for sustainable agriculture. This work provided the first insights into the INCYDE-responsive proteins in Arabidopsis seedlings, and found differences between cytokinin tZ and INCYDE effect at the molecular level. The results showed that the INCYDE inhibition of CKX is different from that of cytokinin accumulation in response to exogenous treatment, and found tissue-specific differences in cytokinin signaling. The inhibitory effect of INCYDE on early seedling growth and development was not found in 14-day-old plantlets, indicating that the response is not only tissue-specific, but also developmentally regulated. Finally, the presented data provide evidence of INCYDE-induced stress response attenuation and a new framework for further detailed investigations of the molecular mechanisms involved in hormonal signaling and stress mitigation.