- Research article
- Open access
- Published:
Proteomics and transcriptomics analyses of Arabidopsis floral buds uncover important functions of ARABIDOPSIS SKP1-LIKE1
BMC Plant Biology volume 16, Article number: 61 (2016)
Abstract
Background
The ARABIDOPSIS SKP1-LIKE1 (ASK1) protein functions as a subunit of SKP1-CUL1-F-box (SCF) E3 ubiquitin ligases. Previous genetic studies showed that ASK1 plays important roles in Arabidopsis flower development and male meiosis. However, the molecular impact of ASK1-containing SCF E3 ubiquitin ligases (ASK1-E3s) on the floral proteome and transcriptome is unknown.
Results
Here we identified proteins that are potentially regulated by ASK1-E3s by comparing floral bud proteomes of wild-type and the ask1 mutant plants. More than 200 proteins were detected in the ask1 mutant but not in wild-type and >300 were detected at higher levels in the ask1 mutant than in wild-type, but their RNA levels were not significantly different between wild-type and ask1 floral buds as shown by transcriptomics analysis, suggesting that they are likely regulated at the protein level by ASK1-E3s. Integrated analyses of floral proteomics and transcriptomics of ask1 and wild-type uncovered several potential aspects of ASK1-E3 functions, including regulation of transcription regulators, kinases, peptidases, and ribosomal proteins, with implications on possible mechanisms of ASK1-E3 functions in floral development.
Conclusions
Our results suggested that ASK1-E3s play important roles in Arabidopsis protein degradation during flower development. This study opens up new possibilities for further functional studies of these candidate E3 substrates.
Background
The ubiquitin-proteasome system (UPS) plays important roles in targeted protein degradation, thereby regulating a variety of cellular processes [1–3]. Ubiquitination reactions are catalyzed by the sequential actions of E1 ubiquitin activating enzymes, E2 ubiquitin conjugating enzymes, and E3 ubiquitin ligases. Multiple ubiquitin molecules can be attached to the existing ubiquitin moieties on the protein substrates to form polyubiquitin chains and the polyubiquitinated proteins are usually then degraded by the 26S proteasome.
The UPS regulates many processes in plants, including development and biotic/abiotic stress responses [1, 3–5]. This broad spectrum of functions is made possible by the large number of genes encoding components in the UPS. Plants usually contain a few E1 enzymes, tens of E2 enzymes, and hundreds of E3 ligases, which determine substrate specificities. Therefore, the numerous E3 ligases can potentially ubiquitinate many proteins. Moreover, the modular design of multimeric E3 ubiquitin ligases including the SKP1-CUL1-F-box (SCF) complexes greatly expands the likely number of proteins that can be specifically ubiquitinated. The subunits of SCF complexes are encoded by multi-gene families, especially the F-box proteins, which are encoded by hundreds of genes in plants. Thus, the combination of these components can form various SCF complexes to ubiquitinate numerous substrate proteins.
Genetic studies indicate that plant F-box proteins are involved in hormone signaling pathways, self-incompatibility, developmental processes, and others. Among the F-box proteins important for hormone signaling, TRANSPORT INHIBITOR RESPONSE 1 (TIR1) is a receptor of auxin and the SCFTIR1 ubiquitin ligase facilitates the degradation of AUX/IAA proteins, which are repressors of auxin-induced gene expression [6–9]. The F-box protein CORONATINE INSENSITIVE 1 (COI1) has a similar mechanism in regulating jasmonic acid (JA) signaling; COI1 is a receptor of JA and SCFCOI1 destabilizes JAZ proteins, thereby releasing the transcription factor MYC2 for the activation of JA-responsive genes [10–12]. Other signaling pathways for hormones such as ethylene, gibberellic acid (GA), and abscisic acid (ABA) also require components of the UPS [5]. In addition, S-locus F-box proteins (SLFs) function as the pollen-specific determinants of self-incompatibility [13–16]. The F-box protein UNUSUAL FLORAL ORGANS (UFO) is important for normal meristem identity and floral organ development [17–19]. UFO can interact with LEAFY genetically to activate AP3 expression [20–22].
The Arabidopsis homolog of the yeast and human SKP1 genes, the ARABIDOPSIS SKP1-LIKE1 (ASK1), encodes an SCF subunit that bridges Cullin and F-box proteins [23]. It has been shown that ASK1 can interact with F-box proteins UFO [22, 24], COI1 [25], TIR1 [6], and others [24, 26, 27]. Since these F-box proteins have important roles in different pathways, ASK1, as a key component in SCFs, likely has crucial functions in many processes. This was suggested by previous genetic studies of the ask1 mutant, which has defects in male meiosis, floral organ development, and vegetative growth [23, 28–31]. Although a few substrates of SCFs have been identified in Arabidopsis, they are mainly specific to the well-studied F-box proteins described above. A large number of other ASK1-interacting F-box proteins and their substrates remain elusive, as do the biological pathways regulated by E3s containing ASK1.
Most of the known ubiquitin ligase substrates were identified by protein-protein interaction methods, usually when the F-box protein has a known function [10, 32–34]. Recently, mass spectrometry (MS)-based proteomics approaches have been increasingly applied in various areas including differential gene expression, post-translational modifications, disease marker discovery, as well as the identification of ubiquitin ligase substrates either by detection of ubiquitinated proteins [35, 36], or by comparing proteomes of wild-type (WT) and ubiquitin ligase mutants [37]. In this study, we used a proteomics approach, Multidimensional Protein Identification Technology (MudPIT), to identify floral proteins potentially regulated by ASK1 by comparing floral bud proteomes of WT and ask1 mutant plants. Furthermore, we performed comparative transcriptomics analysis of WT and ask1 floral buds to investigate the effect of ASK1 on gene expression. The integrated transcriptomics and proteomics analyses revealed that many proteins are potentially regulated by ASK1-E3s. We discuss several possible ways of how ASK1 might regulate protein stability and further downstream gene expression.
Results and discussion
Transcriptomic analysis of Ler and ask1 floral buds
To determine the effect of the ask1 mutation on the floral transcriptome, WT (Ler) and ask1 floral bud transcriptomes were analyzed using GeneChip Arabidopsis ATH1 Genome Array. The average values from Ler and ask1 microarrays were compared to find genes whose RNA levels differ by at least two fold and Student’s t-test p-value < 0.05. We found that 74 and 42 genes were up-regulated and down-regulated, respectively, in ask1 transcriptome compared with Ler (Additional files 1 and 2). We used agriGO [38] to determine if certain gene categories are over-represented in the up-/down-regulated genes in ask1. We found that genes are enriched in the GO categories of responsive to various stimuli or stresses (Fig. 1). Among the 42 down-regulated genes (including ASK1) in ask1, 19 genes are related to biotic/abiotic signaling pathways (Table 1), including hormone, light/circadian, temperature, salt, and other signaling pathways. Among the 74 up-regulated genes in ask1, 39 genes were annotated to be involved in response to various biotic/ abiotic signals (Table 2). The molecular functions of most of these genes are not well understood except for evidence from transcriptional responses to stimuli (e.g., COLD-REGULATED 15A/15B, DARK INDUCIBLE 10, SENESCENCE 1, etc.) and sequence homology with well-characterized proteins or protein domains (e.g., HAD superfamily acid phosphatase, JUMONJI DOMAIN CONTAINING 5, CONSTANS-LIKE 2, etc.).
Nevertheless, several genes have been functionally characterized, including CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), LATE ELONGATED HYPOCOTYL (LHY), JASMONATE-ZIM-DOMAIN PROTEIN 1 (JAZ1), and JAZ5. CCA1 and LHY encode Myb-like transcription factors that synergistically regulate circadian rhythm of Arabidopsis [39] and thus are important for coordinating internal physiological activities with external environmental cues. JAZ genes are induced by JA through a feedback loop involving JAZ proteins and the G-box-binding MYC2: JAZ proteins bind to and repress the activity of MYC2 in the absence of JA; upon perception of JA, JAZ proteins are degraded after ubiquitination by SCFCOI1 and the released MYC2 can activate transcription of downstream genes, including JAZ genes [40]. According to this feedback regulatory model, it is expected that the ask1 mutation would reduce SCF activities, allowing JAZ proteins to accumulate and repress MYC2 activity and thus reducing the JAZ transcript levels. However, we found that JAZ1 and JAZ5 transcript levels were unexpectedly higher in the ask1 mutant than in WT. This paradox suggests that an uncharacterized mechanism may be involved in modulating the JA signaling pathway. For example, ASK1-containing SCFs might facilitate the removal of a yet unidentified transcription activator that has the ability of inducing the expression of JAZ genes in the absence of JA; when ASK1 is mutated this transcriptional activator is stabilized, thereby inducing the expression of downstream genes including JAZ1 and JAZ5. Further studies are needed to uncover new aspects of these regulatory networks.
We then analyzed possible overrepresentation of cis-elements in the putative promoter regions of these up-/down-regulated genes in the ask1 transcriptome. The frequencies of 6-mer motifs within the 500 bp and 1000 bp putative promoter regions were determined using the Motif Analysis tool from The Arabidopsis Information Resource (TAIR) (Table 3). The G-box (CACGTG) is overrepresented in the putative promoter regions of up-/down-regulated genes, suggesting that corresponding genes might be regulated by G-box-binding transcription factors, which themselves or whose co-factors might be regulated by ASK1-E3 ligases, similar to the JAZ-MYC2 model. Some of these transcription factors or co-factors might be short-lived repressors; when they are stabilized in the absence of ASK1, their target genes are then down-regulated. Others may function as unstable activators, whose stabilization in the absence of ASK1 results in up-regulation of downstream genes. Alternatively, some transcription factors may have dual functions, both activation and repression, as is true for MYC2 [41, 42]. The fact that the genes whose promoters contain these cis-elements are altered in transcription in the ask1 mutant suggests that the protein levels of the corresponding transcription factors were changed in ask1. Another motif, GATAAG (I box), was enriched in the down-regulated genes in ask1. The I box was previously found to be enriched in promoters of light-regulated genes [43] and is required for Arabidopsis rbcS-1A expression [44]. Further experiments are required to test whether the putative cis-elements are functional and to identify cognate transcription factors that connect ASK1-E3 regulation with transcriptional changes.
The enrichment of biotic/abiotic stress related genes in the up-/down-regulated genes in the transcriptome of ask1 mutant floral buds has several possible implications. First, the up-regulation of 39 biotic/abiotic stress related genes in ask1 floral buds (Table 2) suggests the expression of such genes might be tightly constrained to avoid unnecessary expression to ensure continuous and maximal allocation of resources to reproductive organs. In WT floral buds, the expression of these genes may be turned off due to degradation of positive transcriptional regulators by ASK1-E3-mediated ubiquitination, but stresses might block the degradation of such positive regulators. Second, the observation that 19 genes annotated as stress responsive were down-regulated in ask1 floral buds (Table 1) compared with WT floral buds suggests their involvement in normal flower development. Although these genes are annotated as responsive to biotic/abiotic signals, they could be triggered by endogenous signals such as programmed cell death (e.g., tapetum degeneration) and/or controlled dehydration during later stages of anther and pollen development [45]. However, the lack of cell-type-specific transcriptome information makes it difficult to determine the extent to which the transcriptome reprogramming for these developmentally-controlled processes resembles stress responses. In summary, ASK1-E3s might destabilize proteins that are involved in the complex regulations of signaling pathways in normal flower development or in response to external stimuli.
Proteomic analysis of Ler and ask1 floral buds
To probe the effect of ask1 on the floral proteome and to identify potential substrates of ASK1-E3s, we used a label-free proteomic method, MudPIT, to analyze floral bud proteomes of the ask1 mutant and Ler (Fig. 2). Total protein extracts of four Ler and five ask1 floral bud samples were digested in-solution with trypsin without pre-separation to maximize digestion of proteins with different properties (e.g., hydrophobicity and charges) and compartmentalization (cytosol, membrane, nucleus and organelles). MudPIT runs of the four Ler samples (Ler-1 ~ Ler_4) detected 2348, 2258, 1658, and 1400 proteins, respectively, with a false discovery rate (FDR) of <1 %, (Additional file 3). When the four datasets were merged, a total of 3220 non-redundant proteins were detected. MudPIT runs of the five ask1 samples (ask1_1 ~ ask1_5) detected 1780, 1441, 1959, 1007, and 363 proteins, respectively, with FDR < 1 %, (Additional file 3), for a total of 2916 non-redundant proteins. The ask1_5 run detected fewer proteins because the starting protein amount was about 20 % of the others to test whether a smaller amount of input protein extract could lead to different efficiency of protein detection. The test result did not show a huge difference in the detection efficiency when the amount of starting material was changed, i.e., the number of detected proteins is proportional to the starting protein amount. The 363 proteins detected in this test run were included in the total ask1 proteins, but excluded for comparison between individual runs with spectral counting normalization in the following sections.
The stochastic sampling and detection sensitivity of MS-based proteomics results in the partial identification of the whole proteome in each experiment and partially overlapping datasets from different MS runs. Indeed, analysis of our Ler and ask1 samples produced partially overlapping sets of proteins (Fig. 3): 884 (27.5 %) of the 3220 total Ler proteins (FDR < 1 %) were detected by all four MS runs, 684 (17.2 %) proteins by three MS runs, 554 (21.2 %) proteins by two MS runs, and 1096 (34.1 %) proteins detected only once. Among the 2899 ask1 proteins (proteins unique to ask1_5 not included), 568 (19.6 %), 493 (17.0 %), 598 (20.6 %), and 1240 (42.8 %) proteins were detected by four, three, two and one of the four MS runs, respectively. Even the proteins detected once had high confidence (FDR < 1 %) resulting from very stringent MS detection and searching criteria and thus were regarded as detected.
Proteins of some cellular component categories (e.g., membrane proteins) that are usually considered to be difficult to be detected by MS without using detergents were well represented in our Ler and ask1 proteomes (Additional file 4), indicating that our method was able to detect proteins localized to the plasma membrane, organelle membrane, and nuclear envelope. One important improvement to minimize bias was that total protein extracts were digested by trypsin without separation into soluble or insoluble portions. Although insoluble proteins were often thought to be recalcitrant to enzyme digestion, extensive sonication, denaturing treatments, and prolonged tryptic digestion seem to have improved detection efficiency.
Detection of additional floral proteins compared with previous proteomics studies
Comparing our floral bud proteomes with previously published proteomics data, we detected many additional proteins (Fig. 4). WT floral bud proteins from two previous studies [46, 47] were combined into one dataset (named “previous WT”) containing 5461 non-redundant proteins (FDR < 1 %). Compared with the previous WT, additional 516 proteins were only detected in our Ler proteome dataset containing 3220 non-redundant proteins and additional 752 proteins were detected in our Ler + ask1 combined proteome dataset containing 3762 non-redundant proteins. The identification of these additional proteins indicates that proteomics detection has not been saturated. The development of more advanced MS technologies with the capacity of overcoming the huge dynamic range of proteins in biological samples is required to further push the proteomic identification closer to saturation. Until then, we have to keep in mind that there is still much room for improvement to our proteomics data and we need to be cautious when drawing conclusions from the current data.
Proteins detected only in ask1 or higher in ask1
Comparison of our ask1 and Ler proteomes revealed that 542 proteins were found only in ask1 (Fig. 5a). These proteins have accumulated in ask1 but not WT, suggesting that they might be putative ASK1-E3 substrates. However, it remains possible that the accumulation of these proteins might be an indirect effect of the ask1 mutation. Since the limited number of MudPIT runs was not enough to saturate the proteome, we cannot rule out the possibility that some proteins detected only in ask1 could have been detected also in Ler if more MS runs had been done. In order to narrow down putative ASK1-E3 substrates, we combined our Ler floral bud proteome data with the previous WT data [46, 47] yielding a larger WT floral bud proteome dataset consisting of 5977 proteins (Pooled WT). This Pooled WT dataset represents the most comprehensive floral bud proteome currently available obtained by MS methods. The total contribution of our WT (Ler) floral bud proteome dataset to this Pooled WT proteome is about 53.9 % (3220/5977). By comparing ask1 and the Pooled WT proteomes, we found that 236 proteins were only detected in ask1 (ask1-only proteins); these are thus more likely to be ASK1-E3 substrates (Fig. 5b and Additional file 5).
We examined RNA levels of these ask1-only proteins from microarray data to determine whether increased transcription contributed to the accumulation of these proteins. The RNA levels of most of these genes were not significantly different between ask1 and Ler (do not meet the criteria of two-fold cutoff and Student’s t-test p-value < 0.05), except that 11 genes lacked probes on the ATH microarray chip (Additional file 5). Therefore, at least 225 proteins seem to be controlled by ASK1-E3s at the protein level.
Because we used a label-free proteomics method, protein abundance cannot be directly compared between samples. One of the relative quantification approaches, spectral counting, has been demonstrated to show higher reproducibility, larger dynamic range, and stronger linear correlation with relative protein abundance than sequence coverage, peptide number, and ion chromatographic methods [46, 48–51]. Therefore, proteomics datasets of this study were normalized using the spectral counting method as reported [46] and the average values were compared between Ler and ask1 to find 322 proteins with higher abundance (1.5-fold cutoff) in ask1 (ask1-higher; Additional file 6). The previous WT data were not included because they were generated by different proteomics methods and difficult to compare quantitatively with our data. We extracted the RNA expression values from microarray data for these ask1-higher proteins to determine whether their elevated protein levels were due to increased transcript levels (Additional file 6). Only the RNA level of AT2G33380 was 2.2-fold higher in ask1 than Ler. The RNA levels of other genes were not significantly different between ask1 and Ler, except that 12 genes had no probes on the microarray chip, suggesting that at least 309 of ask1-higher proteins are probably regulated at the protein level.
ASK1 regulates abundance of regulatory proteins acting at multiple levels
GO categorization of the ask1-only and ask1-higher proteins shows that some molecular functional categories are overrepresented (Fisher test p-value < 0.05) (Fig. 6). Since regulatory proteins are often affected by the UPS, certain categories are of particular interest, including transcriptional regulators, kinases, and peptidases/proteases. Interestingly, many ribosomal proteins were also found to accumulate in ask1 indicating a possible role of the UPS in translational regulation or extraribosomal functions of ribosomal proteins.
ASK1-E3s affects the levels of transcriptional regulators
The ask1 mutation affected 19 transcription factors and chromatin remodelers, including some with functional information (Table 4). One of them was JAZ9, which might function as a repressor of MYC transcription factors in the JA signaling pathway [52], consistent with the SCFCOI1-dependent ubiquitination of JAZ proteins and subsequent degradation upon JA perception [10–12]. Another protein detected in ask1 was the basic helix-loop-helix (bHLH) transcription factor MYC3, which interacts with JAZ proteins and functions with MYC2 and MYC4 to activate JA responses [53, 54]. MYC3 binds to G-boxes [53], possibly regulating promoters with G-boxes, including some of the genes that exhibited altered expression in ask1, as supported by the enrichment of G-boxes in the promoters of the genes exhibiting altered expression in ask1 as described in the first section (Table 3). The accumulation of JAZ9 and MYC3 in the ask1 proteome raises the possibility that protein stability of both JAZs and MYCs is regulated by the UPS, adding another layer of regulation in the JA signaling pathway. In addition, eight zinc finger transcription factors and several transcription factors of other types (Table 4) also exhibited elevated protein levels in ask1, indicating that destabilization of these transcription factors by ASK1-E3s might be important for floral development.
The ask1 mutation also caused an increased level of a Polycomb group protein, SWINGER (SWN) (Table 4), which interacts with other Polycomb group proteins to repress expression of FLOWERING LOCUS C (FLC) and controls the initiation of endosperm development [55–57]. The elevated SWN protein level is expected to repress FLC expression. Although our ask1 and Ler floral bud transcriptome data did not show significant difference of FLC expression, our unpublished anther transcriptome data showed that averaged microarray values of FLC are 68 and 251 (raw microarray data normalized by R package RMA) in ask1 and Ler anther transcriptomes, respectively, i.e., the FLC level is lower in ask1 anthers than that in Ler anthers (regular Student’s t-test p-value = 0.04). The down-regulation of FLC in the ask1 anther transcriptome suggests that degradation of the SWN protein in developing anthers normally derepress its target genes including FLC. Another affected chromatin remodeling protein is RELATIVE OF EARLY FLOWERING 6 (REF6), which is a histone H3 lysine 27 demethylase [58] and positively regulates flowering and brassinosteroid signaling [59]. Our results suggest that ASK1-E3s may modulate gene transcription by facilitating the degradation of both sequence-specific DNA binding transcription factors and chromatin remodelers.
ASK1-E3 affects ribosomal proteins
Interestingly, many ribosomal proteins were identified only or with higher levels in ask1 (Table 5), indicating that ASK1-E3s may also have a role in translational regulation. Several ribosomal proteins have been genetically studied: NUCLEAR FUSION DEFECTIVE 3 is required for polar nuclei fusion during female gametophyte development [60]; PIGGYBACK1 influences leaf vascular patterning [61]; OLIGOCELLULA7 is involved in ribosome biogenesis and organ size control [62]; and POINTED FIRST LEAF2 plays a role in early leaf development [63]. Ribosomal proteins might be regulated by ASK1-E3s either for ribosome turnover, or for extraribosomal regulatory purposes. It was suggested that ribosomal proteins can be ubiquitinated for selective degradation of ribosomes by autophagy [64, 65]. So, the accumulation of ribosomal proteins in ask1 might result from the failure of ubiquitination by ASK1-E3s and subsequent autophagy. A previous study showed that F-BOX PROTEIN 7 (FBP7), which interacts with ASK1, is required for efficient translation under temperature stress conditions, but the substrate of this F-box protein is not identified [66]. It is possible that F-box proteins, such as FBP7, may regulate translation by ubiquitination of specific ribosomal proteins. Our results suggest that ASK1-E3s might be extensively involved in the regulation of translation.
Alternatively, ribosomal proteins may have extraribosomal functions, which are exemplified by the previous findings that several ribosomal proteins can block the ubiquitination of the tumor suppressor p53 upon ribosomal stress [67]. Therefore, the ribosomal proteins accumulated in ask1 might function as regulatory proteins which themselves may be regulated by ubiquitin-mediated proteolysis. Further studies of genetic and molecular interactions between ribosomal proteins and E3 ubiquitin ligases are required for elucidating the role of UPS in translational regulation.
It is difficult to rule out the possibility that the accumulation of these ribosomal proteins might be a side effect of the ask1 mutation. Since the ask1 is a stable mutant, it is possible that long-term and large-scale disturbance of protein degradation might impose an intracellular stress which in turn affects the ribosome biosynthesis, activity or turnover. Therefore, more research is needed to elucidate the link between protein degradation and ribosomal protein functions.
ASK1-E3s affects regulators of protein activities
A number of kinases accumulated in the ask1 proteome (Table 6), including some that are important for plant responses to various stimuli. Among them, CALCIUM-DEPENDENT PROTEIN KINASE 6 (CPK6) is a positive regulator of salt/drought stress tolerance [68], methyl jasmonate signaling in guard cells [69], and ABA regulation of guard cell ion channels [70]. In addition, LYSM DOMAIN RECEPTOR-LIKE KINASE 1 (LYSM RLK1) is involved in chitin-mediated plant innate immunity [71, 72], and MAP KINASE KINASE 2 (MKK2) regulates cold and salt stress signaling and innate immunity [73–75]. Our results indicate that ASK1-E3s normally destabilize these kinases during normal flower development, possibly to suppress biotic/abiotic stress responses in the absence of stimuli.
Also affected are two cell cycle regulators, cyclin-dependent kinase CELL DIVISION CONTROL 2 (CDC2) and a Cyclin T protein, CYCT1;5. CDC2 is required for male gametogenesis [76]. CYCT1;5 is a subunit of cyclin-dependent kinase C complexes involved in cauliflower mosaic virus infection, plant growth and development [77]. The accumulation of cell cycle regulators in ask1 may affect mitosis and/or meiosis, as suggested by the findings that ask1 mutant plants have reduced cell numbers and defective male meiosis [23, 28–31]. Timely removal of cell cycle regulators is likely to be an important part of ASK1 function in regulating plant development.
We also found that SNF1-RELATED PROTEIN KINASE 1.2 (SnRK1.2)/SNF1 KINASE HOMOLOG 11 (KIN11), which was reported to be degraded during phosphate starvation [78], accumulated in the ask1 proteome. SnRK1.2/KIN11 was also shown to interact with ASK1 [79]. These findings imply that ASK1 might directly recruit SnRK1.2/KIN11 to Cul1 without an F-box protein, resulting in ubiquitination and degradation of SnRK1.2/KIN11. Alternatively, an F-box protein might also interact with ASK1 and SnRK1.2/KIN11, forming an SCF complex for ubiquitinating SnRK1.2/KIN11.
Other kinases are largely unknown, but might have important functions in signal perception and transduction. For example, AT5G43020 and AT3G14350 contain transmembrane domains and could be membrane receptor kinases. In summary, the increased levels of these kinases suggest that ASK1-E3s negatively regulate levels of these protein kinases to control cell cycle, plant immunity, hormone signaling, and other processes.
ASK1-E3 affects regulators of protein stability
The peptidase category is enriched in the ask1-higher proteins (Table 7), indicating that ASK1-E3s may regulate degradation of peptidases, which in turn affect protein processing or turnover. Four peptidases (AT1G53850, AT5G66140, AT1G77440, and AT3G60820) are isoforms of 20S proteasome alpha/beta subunits, indicating that the proteasome core complex may also be regulated by UPS. Two ubiquitin-specific proteases UBIQUITIN-SPECIFIC PROTEASE5 (UBP5) and UBP6 were also detected in ask1-higher proteins, suggesting that deubiquitinases, which antagonize protein ubiquitination, might also be regulated by the UPS. The BRI1 SUPPRESSOR 1 (BRS1), a secreted serine carboxypeptidase, is involved in brassinosteroid signaling possibly by processing some proteins [80]. Other peptidases are largely unknown except information from expression and homology. Peptidases/proteases may normally be subject to negative regulation by ASK1-E3s, thus coupling peptidase-mediated protein processing or degradation with the UPS.
Possible ways that ASK1 regulates gene expression
By integrative analysis of transcriptome and proteome data, we found that ASK1-E3s might regulate gene expression at multiple steps, ranging from transcriptional, translational, to post-translational regulations. ASK1-E3s may destabilize transcription repressors or activators to derepress or inactivate gene transcription, respectively (Fig. 7a). In the absence of ASK1, the accumulation of these transcriptional repressors or activators results in down-regulation or up-regulation of gene transcription, respectively. However, we cannot rule out the possibility that the altered transcriptome and proteome might be indirect consequences of the ask1 mutation.
The proteins accumulated in ask1 might be direct substrates of ASK1-E3s, or stabilized by ASK1-E3 substrates (Fig. 7b). For example, ubiquitin-specific proteases UBP5 and UBP6, which accumulate in the ask1 proteome (Table 7), might be substrates of ASK1-E3s; UBP5 and UBP6 could deubiquitinate and prevent degradation of ubiquitinated proteins, whose protein levels are then increased in ask1. An example in human is the herpesvirus-associated ubiquitin-specific protease (HAUSP), which stabilizes a tumor suppressor p53 by deubiquitination [81]. Ribosomal proteins may share a similar mechanism: accumulation of ribosomal proteins in ask1 may increase protein synthesis; alternatively, if ribosomal proteins have extraribosomal regulatory functions, they may stabilize some proteins in a similar way as those stabilizing p53 in human [67].
In another possible scenario, ASK1-E3s may destabilize some proteolytic enzymes (e.g., E3 ubiquitin ligases or peptidases), which can degrade other proteins (Fig. 7c), forming a double negative regulation cascade. The accumulation of such proteolytic enzymes in ask1 may cause reduced levels of their proteolytic substrates. Proteasome subunits and peptidases that accumulate in ask1 may be involved in degradation of their substrate proteins, which could be detected with lower levels in ask1 proteome (Additional file 7). However, it remains difficult to identify these proteolytic substrates due to lack of functional information of the proteolytic enzymes.
There are probably more proteins regulated by ASK1-E3s than those identified in this study. For example, the substrates of the well-studied F-box proteins, TIR1 and COI1, were not detected except JAZ9 (Table 4). One possible reason is that, due to technical limitations, MS might not have uncovered proteins with low-level and/or spatiotemporally restricted expression (e.g, the putative UFO substrate, LEAFY, which is mainly expressed in the inflorescence meristem and floral meristem [20–22]), and when the substrates of F-box proteins are subject to degradation. Another important reason is the functional redundancies among the 21 ASK family members in Arabidopsis. Since the ASK1 gene is expressed throughout the plant with higher levels in growing organs, its mutation is expected to cause more defects in many plant organs. However, the actual defects are milder than the expected, probably due to the redundancy among the ASK family members. The ASK2 gene is the most closely related gene to ASK1. The single mutant of ask2 is similar to WT plants, but the ask1 ask2 double mutant has severe defects in embryo development and is lethal soon after seed germination [82]. This suggests that the redundancy of ASK1 with ASK2, and possibly other ASK proteins, probably has masked some aspects of the ASK1 function. In other words, some ASK1-E3 substrates might also be ubiquitinated by SCFs containing other ASK proteins (e.g., ASK2-E3s), and thus would be unable to accumulate in the ask1 proteome. This might also explain why most of the well-known substrates of F-box proteins (TIR1 and COI1) were not identified in our MS data.
For example, one new aspect of ASK1 functions was revealed by our reanalysis of female fertility of the ask1 mutant, which was reported to be female fertile in previous studies [23, 28–31]. We loaded excess WT (Ler) pollen onto stigmas of the ask1 mutant, the dysfunctional tapetum 1 (dyt1) mutant (as a control with male sterility and female fertility) [83], and Ler (as a self-pollination control), and finally we counted mature seeds from every silique (Additional file 8). The result clearly showed that the pollinated ask1 pistils yielded significantly fewer seeds (16.0 seeds/silique on average) than Ler (52.5 seeds/silique on average) and dyt1 (52.0 seeds/silique on average) (Student’s t-test p-value < 0.01); while the pollinated dyt1 pistils yielded similar numbers of seeds as Ler (Student’s t-test p-value > 0.05). This finding suggests a previously unrecognized role of ASK1 in female reproductive development in Arabidopsis. Studying the masked aspects of ASK1 functions will need tissue-specific silencing of multiple ASK family members, or tissue-specific ASK1 complementation within the ask1 ask2 double mutant or higher order mutants. In addition, characterization of the ubiquitinated proteome may identify potential substrates of E3 ubiquitin ligases and ubiquitination sites within each protein, providing additional clues about ASK1 function in related processes.
Conclusions
Protein degradation is an integral part of various biological processes. The UPS is of particular interest since it selectively degrades proteins, including many key regulators of many cellular pathways [1–3]. However, searching for specific substrates of E3 ubiquitin ligases has been difficult probably due to rapid degradation of substrate proteins once they have been polyubiquitinated by E3 ubiquitin ligases, relatively weak interaction between E3s and substrates, narrow spatiotemporal window where the E3-substate interaction occurs, and others.
In this study, we have searched for potential E3 substrates by using an Arabidopsis mutant that lacks the functional ASK1 gene encoding a key component of SCF-type E3 ubiquitin ligases and that has developmental defects, particularly in floral organs including petals and anthers [23, 28–31]. We employed a MS-based method, MudPIT, to explore floral bud proteomes and detected 2916 and 3220 proteins in ask1 and WT proteomes, respectively. By comparing the ask1 proteome with a pooled WT floral bud proteome (our WT floral bud proteome combined with two published WT floral bud proteomes), we found 236 proteins that are unique to the ask1 proteome and 322 proteins with higher levels in the ask1 proteome. The accumulation of these proteins in the absence of ASK1-E3s suggests that they may be targeted by ASK1-E3s for degradation in WT. Our transcriptomics analysis of ask1 and WT floral buds showed that the transcripts of genes encoding the proteins accumulated in the ask1 proteome are not significantly affected by the ask1 mutation, suggesting that these proteins are regulated at the protein level and thus are more likely to be candidate substrates of ASK1-E3s. Functional categorization revealed that many of the potential substrates of ASK1-E3s are involved in regulation of transcription, translation, protein phosphorylation, and protein degradation. This indicates a multifaceted role of ASK1 in regulating plant development. Much more work is required to validate these candidate E3 substrates and to investigate their specific molecular functions.
Methods
Plant materials and growth conditions
The Arabidopsis thaliana ecotype Landsberg erecta (Ler) and ask1 mutant within the Ler background [23] were used. Plants were grown on soil (Metro-Mix 360, Sun Gro Horticulture, Bellevue, WA) in a growth room with a temperature of 23 °C and long day conditions (16 h light and 8 h dark). The ask1 mutant plants were selected from the progeny of ASK1/ask1 heterozygous plants by their abnormal phenotypes including reduced plant size compared with WT plants of the same age, reduced number and/or reduced size of petals, sterile anthers, short filaments, and short siliques. Clusters of unopened floral buds from the primary inflorescences (from inflorescence meristem to the biggest unopened bud) of the ask1 mutant and Ler were collected from plants with about 5 open flowers.
Microarray analysis
Ler and ask1 floral bud total RNA was extracted using the NucleoSpin® RNA Plant kit (MACHEREY-NAGEL, Bethlehem, PA). RNA quality analysis was performed on the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA), controlled by the Agilent 2100 Expert software, using the Plant RNA Nano assay following the RNA 6000 Nano kit protocol. Microarray was performed using the GeneChip Arabidopsis ATH1 Genome Array (Affymetrix, Santa Clara, CA) in the Penn State Genomics Core Facility – University Park, PA. Three biological replicates of ask1 and four biological replicates of Ler were performed (Additional file 9).
Data analysis was conducted as previously described with some modifications [84]. Microarray datasets (.CEL files) were normalized by R package RMA and exported as Excel files. Microarray signal values were averaged from biological replicates of each genotype and compared between ask1 and Ler to find differentially expressed genes which show at least 2-fold differences in RNA levels and p-value < 0.05 (regular Student’s t-test). GO categorization was conducted using the Singular Enrichment Analysis (SEA) from agriGO [38]. The Affymetrix ATH1 Genome Array (GPL198) was selected as the background reference which contains 22479 annotated genes. The statistical test was set to Fisher and significance level set to 0.05.
Protein extraction with trichloroacetic acid/acetone method
The protein extraction method was modified from a previous study [85]. Floral buds were ground thoroughly in liquid nitrogen with mortars and pestles and the powder was suspended in -20 °C Acetone with 10 % w/v Trichloroacetic Acid (TCA) and 0.07 % (v/v) β-Mercaptoethanol (1 ml for 0.3 g of tissue powder). After being incubated for 2 h (or overnight) at -20 °C, the protein suspension was centrifuged for 15-20 min at 14,000 rpm. The supernatant was removed and the protein pellet was resuspended and washed with 1 ml of -20 °C Acetone containing 0.07 % (v/v) β-Mercaptoethanol followed by centrifugation for 15-20 min at 14,000 rpm. This washing step was repeated until the pellet was almost white. The protein pellet was vacuum dried for 5-10 min and stored at -20 °C or immediately used for trypsin digestion.
In-solution trypsin digestion of protein extract
About 20-30 mg of crude protein extract from the TCA/Acetone method was resuspended in 1 ml of rehydration buffer [100 mM NH4HCO3, 10 mM Dithiothreitol (DTT), 10 % (v/v) Acetonitrile] and sonicated for 5 times, 20 s each time, duty cycle 40 %, power 3 using a Branson Sonifier S-450A (Branson Ultrasonics, Danbury, CT). Proteins were denatured at 60 °C for 45–60 min and alkylated by 50 mM Iodoacetamide at 37 °C for 30 min in dark. 40 μl of 1 M DTT was added to quench the alkylation reaction. Alkylated proteins were digested by 20 μg of Trypsin Gold, Mass Spectrometry Grade (Promega, Madison, WI) for 16-18 h at 37 °C with moderate shaking. The remaining indigestible debris was removed by centrifugation at 12,000 rpm for 10 min. The supernatant was transferred to a new 1.5 ml tube and centrifuged again to remove residual debris. The supernatant was transferred to a new 1.5 ml tube and was adjust to pH 3.0 with glacial acetic acid. The peptide solution was vacuum dried completely to evaporate off NH4HCO3 and acetonitrile. The pellet was resuspended in 200 μl of H2O and vacuum dried. Three repeats of resuspension and drying were performed in total. Finally the peptides were analyzed in the Proteomics and Mass Spec Core Facility, College of Medicine, Pennsylvania State University, Hershey, PA.
Mass spectrometry analysis/MudPIT
Trypsin-digested peptide samples were analyzed by MudPIT according to the 2D LC-MALDI separation and analysis procedures published previously using a 4800 proteomic analyzer MALDI TOF/TOF tandem system (Applied Biosysems) [86] except several modifications. The ProteinPilot software version 4.2 was used to perform protein identification by searching MS spectra against the protein database which included the Arabidopsis thaliana protein list TAIR10_pep_20101214, 156 common human and lab contaminants (ABSciex_ContaminantDB_20070711), and a reverse “decoy” version of the protein database itself (concatenated Reverse Decoy Database). Proteins with local FDR < 1 % were accepted as detected (Additional files 10, 11, 12, 13, 14, 15, 16, 17, and 18).
Proteomics data analysis
We combined proteins detected in ask1 samples into the ask1 proteome, and combined proteins detected in Ler into the Ler proteome. We first compared our ask1 and Ler proteomes to find proteins that are only detected in ask1 samples. We also obtained previously published proteomics data of wild-type Arabidopsis thaliana floral buds [46, 47] and combined them into a “previous WT” proteome containing 5461 non-redundant proteins (FDR < 1 %). Comparison of our Ler and ask1 proteomes with the previous WT proteome resulted in the finding of additional floral proteins in our data. We combined our Ler floral bud proteome with the previous WT proteome [46, 47] to a “Pooled WT” proteome consisting of 5977 non-redundant proteins. Comparison of our ask1 proteome with the Pooled WT proteome led to the identification of proteins that are considered as “ask1-only” proteins.
Since each sample was analyzed by MudPIT individually without labeling and multiplexing, the abundance of each protein cannot be directly compared across different samples. Instead, the relative abundance of each protein in a sample was normalized using the spectral counting method as previously described [46, 48–51, 87]. The following formula is used to calculate the spectral counting values which represent the normalized relative abundance of proteins:
The “Measured spectra of protein K” is the number of actually detected MS spectra that specifically match to the protein K.
The “Measured spectra of all proteins in dataset” is the sum of the measured spectra of proteins in one sample.
The “Theoretical peptides of all proteins in dataset” is the total number of the in silico tryptic peptides of all proteins detected in one sample. The in silico tryptic digestion was carried out using the tool “digest” from the Galaxy platform (https://usegalaxy.org/). Since trypsin normally does not cut after lysine (K) or arginine (R) residues if it is followed by a Proline (P), we specified these sites as non-cut sites. Partial digestion and fragments containing more than one potential cut site were not included. Peptides containing at least 6 amino acid residues were considered as theoretical peptides.
The “Theoretical peptides of protein K” is the number of tryptic peptides of a protein K that was determined in the above “Theoretical peptides of all proteins in dataset”.
For a protein detected in both ask1 and Ler samples, its spectral counting values were averaged across ask1 and Ler samples, respectively. Then the average spectral counting values of a protein in ask1 and Ler samples were compared. Proteins whose average spectral counting value in ask1 samples is at least 1.5-fold of that in Ler samples were considered as “ask1-higher” proteins.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional file. The raw microarray datasets were deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) with the accession number GSE42841.
Abbreviations
- ASK1:
-
ARABIDOPSIS SKP1-LIKE1
- SCF:
-
SKP1-CUL1-F-box
- ASK1-E3s:
-
ASK1-containing SCF E3 ubiquitin ligases
- UPS:
-
Ubiquitin-proteasome system
- TIR1:
-
TRANSPORT INHIBITOR RESPONSE 1
- COI1:
-
CORONATINE INSENSITIVE 1
- JA:
-
Jasmonic acid
- GA:
-
Gibberellic acid
- ABA:
-
Abscisic acid
- SLFs:
-
S-locus F-box proteins
- UFO:
-
UNUSUAL FLORAL ORGANS
- MS:
-
Mass spectrometry
- MudPIT:
-
Multidimensional Protein Identification Technology
- CCA1:
-
CIRCADIAN CLOCK ASSOCIATED 1
- LHY:
-
LATE ELONGATED HYPOCOTYL
- JAZ1:
-
JASMONATE-ZIM-DOMAIN PROTEIN 1
- TAIR:
-
The Arabidopsis Information Resource
- FDR:
-
False discovery rate
- WT:
-
Wild-type
- bHLH:
-
Basic helix-loop-helix
- SWN:
-
SWINGER
- FLC:
-
FLOWERING LOCUS C
- REF6:
-
RELATIVE OF EARLY FLOWERING 6
- FBP7:
-
F-BOX PROTEIN 7
- CPK6:
-
CALCIUM-DEPENDENT PROTEIN KINASE 6
- LYSM RLK1:
-
LYSM DOMAIN RECEPTOR-LIKE KINASE 1
- MKK2:
-
MAP KINASE KINASE 2
- CDC2:
-
Cell Division Cycle 2
- SnRK1.2:
-
SNF1-RELATED PROTEIN KINASE 1.2
- KIN11:
-
SNF1 KINASE HOMOLOG 11
- UBP:
-
UBIQUITIN-SPECIFIC PROTEASE
- BRS1:
-
BRI1 SUPPRESSOR 1
- HAUSP:
-
Herpesvirus-associated ubiquitin-specific protease
- DYT1:
-
DYSFUNCTIONAL TAPETUM 1
- Ler :
-
Landsberg erecta
- SEA:
-
Singular Enrichment Analysis
- NCBI:
-
National Center for Biotechnology Information
- GEO:
-
Gene Expression Omnibus
- TCA:
-
Trichloroacetic Acid
- DTT:
-
Dithiothreitol
References
Smalle J, Vierstra RD. The ubiquitin 26S proteasome proteolytic pathway. Annu Rev Plant Biol. 2004;55(1):555–90.
Sorokin AV, Kim ER, Ovchinnikov LP. Proteasome system of protein degradation and processing. Biochemistry (Moscow). 2009;74(13):1411–42.
Schwechheimer C, Willige BC, Zourelidou M, Dohmann EM. Examining protein stability and its relevance for plant growth and development. Methods Mol Biol. 2009;479:147–71.
Dreher K, Callis J. Ubiquitin, hormones and biotic stress in plants. Ann Bot. 2007;99(5):787–822.
Stone SL, Callis J. Ubiquitin ligases mediate growth and development by promoting protein death. Curr Opin Plant Biol. 2007;10(6):624–32.
Gray WM, del Pozo JC, Walker L, Hobbie L, Risseeuw E, Banks T, et al. Identification of an SCF ubiquitin-ligase complex required for auxin response in Arabidopsis thaliana. Genes Dev. 1999;13(13):1678–91.
Gray WM, Kepinski S, Rouse D, Leyser O, Estelle M. Auxin regulates SCFTIR1-dependent degradation of AUX/IAA proteins. Nature. 2001;414(6861):271–6.
Kepinski S, Leyser O. The Arabidopsis F-box protein TIR1 is an auxin receptor. Nature. 2005;435(7041):446–51.
Dharmasiri N, Dharmasiri S, Estelle M. The F-box protein TIR1 is an auxin receptor. Nature. 2005;435(7041):441–5.
Thines B, Katsir L, Melotto M, Niu Y, Mandaokar A, Liu G, et al. JAZ repressor proteins are targets of the SCFCOI1 complex during jasmonate signalling. Nature. 2007;448(7154):661–5.
Xu LH, Liu FQ, Lechner E, Genschik P, Crosby WL, Ma H, et al. The SCFCOl1 ubiquitin-ligase complexes are required for jasmonate response in Arabidopsis. Plant Cell. 2002;14(8):1919–35.
Devoto A, Nieto-Rostro M, Xie D, Ellis C, Harmston R, Patrick E, et al. COI1 links jasmonate signalling and fertility to the SCF ubiquitin-ligase complex in Arabidopsis. Plant J. 2002;32(4):457–66.
Hua Z, Meng X, Kao TH. Comparison of Petunia inflata S-Locus F-box protein (Pi SLF) with Pi SLF like proteins reveals its unique function in S-RNase based self-incompatibility. Plant Cell. 2007;19(11):3593–609.
Hua Z, Kao TH. Identification and characterization of components of a putative petunia S-locus F-box-containing E3 ligase complex involved in S-RNase-based self-incompatibility. Plant Cell. 2006;18(10):2531–53.
Fields AM, Wang N, Hua Z, Meng X, Kao TH. Functional characterization of two chimeric proteins between a Petunia inflata S-locus F-box protein, PiSLF2, and a PiSLF-like protein, PiSLFLb-S2. Plant Mol Biol. 2010;74(3):279–92.
Lai Z, Ma W, Han B, Liang L, Zhang Y, Hong G, et al. An F-box gene linked to the self-incompatibility (S) locus of Antirrhinum is expressed specifically in pollen and tapetum. Plant Mol Biol. 2002;50(1):29–42.
Ingram GC, Goodrich J, Wilkinson MD, Simon R, Haughn GW, Coen ES. Parallels between UNUSUAL FLORAL ORGANS and FIMBRIATA, genes controlling flower development in Arabidopsis and Antirrhinum. Plant Cell. 1995;7(9):1501–10.
Wilkinson MD, Haughn GW. UNUSUAL FLORAL ORGANS controls meristem identity and organ primordia fate in Arabidopsis. Plant Cell. 1995;7(9):1485–99.
Levin JZ, Meyerowitz EM. UFO: an Arabidopsis gene involved in both floral meristem and floral organ development. Plant Cell. 1995;7(5):529–48.
Chae E, Tan QK, Hill TA, Irish VF. An Arabidopsis F-box protein acts as a transcriptional co-factor to regulate floral development. Development. 2008;135(7):1235–45.
Lee I, Wolfe DS, Nilsson O, Weigel D. A LEAFY co-regulator encoded by UNUSUAL FLORAL ORGANS. Curr Biol. 1997;7(2):95–104.
Samach A, Klenz JE, Kohalmi SE, Risseeuw E, Haughn GW, Crosby WL. The UNUSUAL FLORAL ORGANS gene of Arabidopsis thaliana is an F-box protein required for normal patterning and growth in the floral meristem. Plant J. 1999;20(4):433–45.
Yang M, Hu Y, Lodhi M, McCombie WR, Ma H. The Arabidopsis SKP1-LIKE1 gene is essential for male meiosis and may control homologue separation. Proc Natl Acad Sci U S A. 1999;96(20):11416–21.
Takahashi N, Kuroda H, Kuromori T, Hirayama T, Seki M, Shinozaki K, et al. Expression and interaction analysis of Arabidopsis Skp1-related genes. Plant Cell Physiol. 2004;45(1):83–91.
Dai L, Xu L, Huang D, Li X, Luo K, Guan C. ASK1 physically interacts with COI1 and is required for male fertility in Arabidopsis. Sci China Ser C: Life Sci. 2002;45(6):631–6.
Woo HR, Chung KM, Park JH, Oh SA, Ahn T, Hong SH, et al. ORE9, an F-box protein that regulates leaf senescence in Arabidopsis. Plant Cell. 2001;13(8):1779–90.
Han L, Mason M, Risseeuw EP, Crosby WL, Somers DE. Formation of an SCFZTL complex is required for proper regulation of circadian timing. Plant J. 2004;40(2):291–301.
Zhao D, Yu Q, Chen M, Ma H. The ASK1 gene regulates B function gene expression in cooperation with UFO and LEAFY in Arabidopsis. Development. 2001;128(14):2735–46.
Zhao D, Yang X, Quan L, Timofejeva L, Rigel NW, Ma H, et al. ASK1, a SKP1 homolog, is required for nuclear reorganization, presynaptic homolog juxtaposition and the proper distribution of cohesin during meiosis in Arabidopsis. Plant Mol Biol. 2006;62(1-2):99–110.
Zhao D, Yang M, Solava J, Ma H. The ASK1 gene regulates development and interacts with the UFO gene to control floral organ identity in Arabidopsis. Dev Genet. 1999;25(3):209–23.
Zhao D, Han T, Risseeuw E, Crosby WL, Ma H. Conservation and divergence of ASK1 and ASK2 gene functions during male meiosis in Arabidopsis thaliana. Plant Mol Biol. 2003;53(1-2):163–73.
Jin J, Ang XL, Shirogane T, Wade Harper J. Identification of substrates for F-box proteins. Methods Enzymol. 2005;399:287–309.
Tan MK, Lim HJ, Harper JW. SCF(FBXO22) regulates histone H3 lysine 9 and 36 methylation levels by targeting histone demethylase KDM4A for ubiquitin-mediated proteasomal degradation. Mol Cell Biol. 2011;31(18):3687–99.
An F, Zhao Q, Ji Y, Li W, Jiang Z, Yu X, et al. Ethylene-induced stabilization of ETHYLENE INSENSITIVE3 and EIN3-LIKE1 is mediated by proteasomal degradation of EIN3 binding F-box 1 and 2 that requires EIN2 in Arabidopsis. Plant Cell. 2010;22(7):2384–401.
Igawa T, Fujiwara M, Takahashi H, Sawasaki T, Endo Y, Seki M, et al. Isolation and identification of ubiquitin-related proteins from Arabidopsis seedlings. J Exp Bot. 2009;60(11):3067–73.
Maor R, Jones A, Nuhse TS, Studholme DJ, Peck SC, Shirasu K. Multidimensional protein identification technology (MudPIT) analysis of ubiquitinated proteins in plants. Mol Cell Proteomics. 2007;6(4):601–10.
Burande CF, Heuze ML, Lamsoul I, Monsarrat B, Uttenweiler-Joseph S, Lutz PG. A label-free quantitative proteomics strategy to identify E3 ubiquitin ligase substrates targeted to proteasome degradation. Mol Cell Proteomics. 2009;8(7):1719–27.
Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38(Web Server issue):W64–70.
Lu SX, Knowles SM, Andronis C, Ong MS, Tobin EM. CIRCADIAN CLOCK ASSOCIATED1 and LATE ELONGATED HYPOCOTYL function synergistically in the circadian clock of Arabidopsis. Plant Physiol. 2009;150(2):834–43.
Chini A, Fonseca S, Fernandez G, Adie B, Chico JM, Lorenzo O, et al. The JAZ family of repressors is the missing link in jasmonate signalling. Nature. 2007;448(7154):666–71.
Yadav V, Mallappa C, Gangappa SN, Bhatia S, Chattopadhyay S. A basic helix-loop-helix transcription factor in Arabidopsis, MYC2, acts as a repressor of blue light-mediated photomorphogenic growth. Plant Cell. 2005;17(7):1953–66.
Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, et al. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell. 2007;19(7):2225–45.
Giuliano G, Pichersky E, Malik VS, Timko MP, Scolnik PA, Cashmore AR. An evolutionarily conserved protein binding sequence upstream of a plant light-regulated gene. Proc Natl Acad Sci U S A. 1988;85(19):7089–93.
Donald RG, Cashmore AR. Mutation of either G box or I box sequences profoundly affects expression from the Arabidopsis rbcS-1A promoter. EMBO J. 1990;9(6):1717–26.
Rogers HJ. Programmed cell death in floral organs: how and why do flowers die? Ann Bot. 2006;97(3):309–15.
Baerenfaller K, Grossmann J, Grobei MA, Hull R, Hirsch-Hoffmann M, Yalovsky S, et al. Genome-scale proteomics reveals Arabidopsis thaliana gene models and proteome dynamics. Science. 2008;320(5878):938–41.
Feng B, Li L, Zhou X, Stanley B, Ma H. Analysis of the Arabidopsis floral proteome: detection of over 2 000 proteins and evidence for posttranslational modifications. J Integr Plant Biol. 2009;51(2):207–23.
Washburn MP, Wolters D, Yates 3rd JR. Large-scale analysis of the yeast proteome by multidimensional protein identification technology. Nat Biotechnol. 2001;19(3):242–7.
Liu H, Sadygov RG, Yates 3rd JR. A model for random sampling and estimation of relative protein abundance in shotgun proteomics. Anal Chem. 2004;76(14):4193–201.
Florens L, Carozza MJ, Swanson SK, Fournier M, Coleman MK, Workman JL, et al. Analyzing chromatin remodeling complexes using shotgun proteomics and normalized spectral abundance factors. Methods. 2006;40(4):303–11.
Paoletti AC, Parmely TJ, Tomomori-Sato C, Sato S, Zhu D, Conaway RC, et al. Quantitative proteomic analysis of distinct mammalian Mediator complexes using normalized spectral abundance factors. Proc Natl Acad Sci U S A. 2006;103(50):18928–33.
Mosblech A, Thurow C, Gatz C, Feussner I, Heilmann I. Jasmonic acid perception by COI1 involves inositol polyphosphates in Arabidopsis thaliana. Plant J. 2011;65(6):949–57.
Fernandez-Calvo P, Chini A, Fernandez-Barbero G, Chico JM, Gimenez-Ibanez S, Geerinck J, et al. The Arabidopsis bHLH transcription factors MYC3 and MYC4 are targets of JAZ repressors and act additively with MYC2 in the activation of jasmonate responses. Plant Cell. 2011;23(2):701–15.
Cheng Z, Sun L, Qi T, Zhang B, Peng W, Liu Y, et al. The bHLH transcription factor MYC3 interacts with the Jasmonate ZIM-domain proteins to mediate jasmonate response in Arabidopsis. Mol Plant. 2011;4(2):279–88.
Wood CC, Robertson M, Tanner G, Peacock WJ, Dennis ES, Helliwell CA. The Arabidopsis thaliana vernalization response requires a polycomb-like protein complex that also includes VERNALIZATION INSENSITIVE 3. Proc Natl Acad Sci U S A. 2006;103(39):14631–6.
Wang D, Tyson MD, Jackson SS, Yadegari R. Partially redundant functions of two SET-domain polycomb-group proteins in controlling initiation of seed development in Arabidopsis. Proc Natl Acad Sci U S A. 2006;103(35):13244–9.
Chanvivattana Y, Bishopp A, Schubert D, Stock C, Moon YH, Sung ZR, et al. Interaction of Polycomb-group proteins controlling flowering in Arabidopsis. Development. 2004;131(21):5263–76.
Lu FL, Cui X, Zhang SB, Jenuwein T, Cao XF. Arabidopsis REF6 is a histone H3 lysine 27 demethylase. Nat Genet. 2011;43(7):715–U144.
Yu X, Li L, Guo M, Chory J, Yin Y. Modulation of brassinosteroid-regulated gene expression by Jumonji domain-containing proteins ELF6 and REF6 in Arabidopsis. Proc Natl Acad Sci U S A. 2008;105(21):7618–23.
Portereiko MF, Sandaklie-Nikolova L, Lloyd A, Dever CA, Otsuga D, Drews GN. NUCLEAR FUSION DEFECTIVE1 encodes the Arabidopsis RPL21M protein and is required for karyogamy during female gametophyte development and fertilization. Plant Physiol. 2006;141(3):957–65.
Pinon V, Etchells JP, Rossignol P, Collier SA, Arroyo JM, Martienssen RA, et al. Three PIGGYBACK genes that specifically influence leaf patterning encode ribosomal proteins. Development. 2008;135(7):1315–24.
Fujikura U, Horiguchi G, Ponce MR, Micol JL, Tsukaya H. Coordination of cell proliferation and cell expansion mediated by ribosome-related processes in the leaves of Arabidopsis thaliana. Plant J. 2009;59(3):499–508.
Ito T, Kim GT, Shinozaki K. Disruption of an Arabidopsis cytoplasmic ribosomal protein S13-homologous gene by transposon-mediated mutagenesis causes aberrant growth and development. Plant J. 2000;22(3):257–64.
Kraft C, Deplazes A, Sohrmann M, Peter M. Mature ribosomes are selectively degraded upon starvation by an autophagy pathway requiring the Ubp3p/Bre5p ubiquitin protease. Nat Cell Biol. 2008;10(5):602–10.
Tagwerker C, Flick K, Cui M, Guerrero C, Dou Y, Auer B, et al. A tandem affinity tag for two-step purification under fully denaturing conditions: application in ubiquitin profiling and protein complex identification combined with in vivo cross-linking. Mol Cell Proteomics. 2006;5(4):737–48.
Calderon-Villalobos LI, Nill C, Marrocco K, Kretsch T, Schwechheimer C. The evolutionarily conserved Arabidopsis thaliana F-box protein AtFBP7 is required for efficient translation during temperature stress. Gene. 2007;392(1-2):106–16.
Zhou X, Hao Q, Liao J, Zhang Q, Lu H. Ribosomal protein S14 unties the MDM2-p53 loop upon ribosomal stress. Oncogene. 2013;32(3):388–96.
Xu J, Tian YS, Peng RH, Xiong AS, Zhu B, Jin XF, et al. AtCPK6, a functionally redundant and positive regulator involved in salt/drought stress tolerance in Arabidopsis. Planta. 2010;231(6):1251–60.
Munemasa S, Hossain MA, Nakamura Y, Mori IC, Murata Y. The Arabidopsis calcium-dependent protein kinase, CPK6, functions as a positive regulator of methyl jasmonate signaling in guard cells. Plant Physiol. 2011;155(1):553–61.
Mori IC, Murata Y, Yang Y, Munemasa S, Wang YF, Andreoli S, et al. CDPKs CPK6 and CPK3 function in ABA regulation of guard cell S-type anion- and Ca(2+)-permeable channels and stomatal closure. PLoS Biol. 2006;4(10), e327.
Gimenez-Ibanez S, Ntoukakis V, Rathjen JP. The LysM receptor kinase CERK1 mediates bacterial perception in Arabidopsis. Plant Signal Behav. 2009;4(6):539–41.
Iizasa E, Mitsutomi M, Nagano Y. Direct binding of a plant LysM receptor-like kinase, LysM RLK1/CERK1, to chitin in vitro. J Biol Chem. 2010;285(5):2996–3004.
Teige M, Scheikl E, Eulgem T, Doczi R, Ichimura K, Shinozaki K, et al. The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis. Mol Cell. 2004;15(1):141–52.
Brader G, Djamei A, Teige M, Palva ET, Hirt H. The MAP kinase kinase MKK2 affects disease resistance in Arabidopsis. Mol Plant Microbe Interact. 2007;20(5):589–96.
Gao M, Liu J, Bi D, Zhang Z, Cheng F, Chen S, et al. MEKK1, MKK1/MKK2 and MPK4 function together in a mitogen-activated protein kinase cascade to regulate innate immunity in plants. Cell Res. 2008;18(12):1190–8.
Iwakawa H, Shinmyo A, Sekine M. Arabidopsis CDKA;1, a cdc2 homologue, controls proliferation of generative cells in male gametogenesis. Plant J. 2006;45(5):819–31.
Cui X, Fan B, Scholz J, Chen Z. Roles of Arabidopsis cyclin-dependent kinase C complexes in cauliflower mosaic virus infection, plant growth, and development. Plant Cell. 2007;19(4):1388–402.
Fragoso S, Espindola L, Paez-Valencia J, Gamboa A, Camacho Y, Martinez-Barajas E, et al. SnRK1 isoforms AKIN10 and AKIN11 are differentially regulated in Arabidopsis plants under phosphate starvation. Plant Physiol. 2009;149(4):1906–16.
Farras R, Ferrando A, Jasik J, Kleinow T, Okresz L, Tiburcio A, et al. SKP1-SnRK protein kinase interactions mediate proteasomal binding of a plant SCF ubiquitin ligase. EMBO J. 2001;20(11):2742–56.
Li J, Lease KA, Tax FE, Walker JC. BRS1, a serine carboxypeptidase, regulates BRI1 signaling in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2001;98(10):5916–21.
Li M, Chen D, Shiloh A, Luo J, Nikolaev AY, Qin J, et al. Deubiquitination of p53 by HAUSP is an important pathway for p53 stabilization. Nature. 2002;416(6881):648–53.
Liu F, Ni W, Griffith ME, Huang Z, Chang C, Peng W, et al. The ASK1 and ASK2 genes are essential for Arabidopsis early development. Plant Cell. 2004;16(1):5–20.
Zhang W, Sun Y, Timofejeva L, Chen C, Grossniklaus U, Ma H. Regulation of Arabidopsis tapetum development and function by DYSFUNCTIONAL TAPETUM1 (DYT1) encoding a putative bHLH transcription factor. Development. 2006;133(16):3085–95.
Feng B, Lu D, Ma X, Peng Y, Sun Y, Ning G, et al. Regulation of the Arabidopsis anther transcriptome by DYT1 for pollen development. Plant J. 2012;72(4):612–24.
Damerval C, De Vienne D, Zivy M, Thiellement H. Technical improvements in two-dimensional electrophoresis increase the level of genetic variation detected in wheat-seedling proteins. Electrophoresis. 1986;7(1):52–4.
Zhao Z, Zhang W, Stanley BA, Assmann SM. Functional proteomics of Arabidopsis thaliana guard cells uncovers new stomatal signaling pathways. Plant Cell. 2008;20(12):3210–26.
Zhu W, Smith JW, Huang CM. Mass spectrometry-based label-free quantitative proteomics. J Biomed Biotechnol. 2010;2010:840518.
Acknowledgements
We thank Anne Stanley (Proteins and Mass Spectrometry Core Research Facility, Penn State College of Medicine, Hershey, PA) for mass spectrometry technical services; and Craig Praul (Director of Expression Analysis, Penn State Genomics Core Facility–University Park, PA) for microarray technical services. This work was supported by a US Department of Energy grant (DE-FG02-02ER15332) to H.M. and by the Department of Biology Huck Institutes of the Life Sciences, Pennsylvania State University (University Park, PA) and Fudan University (Shanghai, China).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
DL and HM conceived and designed the experiments. DL and WN performed the experiments. DL and BAS analyzed the data. DL wrote and revised the manuscript and HM edited the manuscript. All authors read and approved the final manuscript.
Additional files
Additional file 1:
Genes up-regulated in the ask1 mutant microarray compared with L er.
Additional file 2:
Genes down-regulated in the ask1 mutant microarray compared with L er.
Additional file 3:
Proteins detected in individual proteomics samples with FDR < 1 %.
Additional file 4:
Three membrane protein categories are well represented in the detected L er and ask1 proteins (FDR < 1 %).
Additional file 5:
ask1 -only proteins with microarray values.
Additional file 6:
ask1 -higher proteins with spectral counting and microarray values.
Additional file 7:
ask1 -lower proteins with spectral counting and microarray values.
Additional file 8:
The ask1 mutant has reduced female fertility.
Additional file 9:
Normalized microarray data of three biological replicates of ask1 and four biological replicates of L er floral buds.
Additional file 10:
Mass spectrometry data of ask1 _1 sample exported as .zip file.
Additional file 11:
Mass spectrometry data of ask1 _2 sample exported as .zip file.
Additional file 12:
Mass spectrometry data of ask1 _3 sample exported as .zip file.
Additional file 13:
Mass spectrometry data of ask1 _4 sample exported as .zip file.
Additional file 14:
Mass spectrometry data of ask1 _5 sample exported as .zip file.
Additional file 15:
Mass spectrometry data of L er_ 1 sample exported as .zip file.
Additional file 16:
Mass spectrometry data of L er_ 2 sample exported as .zip file.
Additional file 17:
Mass spectrometry data of L er_ 3 sample exported as .zip file.
Additional file 18:
Mass spectrometry data of L er_ 4 sample exported as .zip file.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lu, D., Ni, W., Stanley, B.A. et al. Proteomics and transcriptomics analyses of Arabidopsis floral buds uncover important functions of ARABIDOPSIS SKP1-LIKE1. BMC Plant Biol 16, 61 (2016). https://doi.org/10.1186/s12870-015-0571-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-015-0571-9