Next Article in Journal
Nonlinear Cellular Mechanical Behavior Adaptation to Substrate Mechanics Identified by Atomic Force Microscope
Previous Article in Journal
Mechanistic Analysis of Physicochemical Cues in Promoting Human Pluripotent Stem Cell Self-Renewal and Metabolism
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 11
10.3390/ijms19113460
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

iTRAQ-Based Analysis of Proteins Co-Regulated by Brassinosteroids and Gibberellins in Rice Embryos during Seed Germination

by
Qian-Feng Li
1,2,
Jin-Dong Wang
1,
Min Xiong
1,
Ke Wei
1,
Peng Zhou
1,
Li-Chun Huang
1,
Chang-Quan Zhang
1,2,
Xiao-Lei Fan
1,2 and
Qiao-Quan Liu
1,2,*
1
Key Laboratory of Crop Genetics and Physiology of Jiangsu Province/Key Laboratory of Plant Functional Genomics of the Ministry of Education/Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding, College of Agriculture, Yangzhou University, Yangzhou 225009, China
2
Co-Innovation Center for Modern Production Technology of Grain Crops of Jiangsu Province/Joint International Research Laboratory of Agriculture and Agri-Product Safety of the Ministry of Education, Yangzhou University, Yangzhou 225009, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(11), 3460; https://doi.org/10.3390/ijms19113460
Submission received: 10 September 2018 / Revised: 27 October 2018 / Accepted: 1 November 2018 / Published: 4 November 2018
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Seed germination, a pivotal process in higher plants, is precisely regulated by various external and internal stimuli, including brassinosteroid (BR) and gibberellin (GA) phytohormones. The molecular mechanisms of crosstalk between BRs and GAs in regulating plant growth are well established. However, whether BRs interact with GAs to coordinate seed germination remains unknown, as do their common downstream targets. In the present study, 45 differentially expressed proteins responding to both BR and GA deficiency were identified using isobaric tags for relative and absolute quantification (iTRAQ) proteomic analysis during seed germination. The results indicate that crosstalk between BRs and GAs participates in seed germination, at least in part, by modulating the same set of responsive proteins. Moreover, most targets exhibited concordant changes in response to BR and GA deficiency, and gene ontology (GO) indicated that most possess catalytic activity and are involved in various metabolic processes. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) analysis was used to construct a regulatory network of downstream proteins mediating BR- and GA-regulated seed germination. The mutation of GRP, one representative target, notably suppressed seed germination. Our findings not only provide critical clues for validating BR–GA crosstalk during rice seed germination, but also help to optimise molecular regulatory networks.

Graphical Abstract

1. Introduction

As sessile organisms, plants need to modulate their growth and development constantly to adapt to the changing environment. A variety of environmental cues, including photoperiod, temperature, and biotic and abiotic stress, are sensed by plants and acted upon by intrinsic hormonal pathways to help determine whether it is appropriate to grow or not [1,2]. At present, the known plant hormone family includes at least nine classes, including auxins, gibberellins (GAs), abscisic acid (ABA), cytokinin (CKs), ethylene, brassinosteroids (BRs), jasmonate (JA), salicylic acid (SA), and strigolactones [3,4]. Although the major functions of each hormone class have been clearly revealed by a series of physiological and genetic studies, the final output that affects plant growth and development is determined by a complex regulatory network mediated by hormonal interactions, and the importance of hormone–hormone interplay is becoming more obvious [5,6,7,8].
Among phytohormones, BRs and GAs are two major growth-promoting hormones with extensive overlapping functions in various plant growth and developmental processes including cell elongation, seed germination, flowering, male fertility and senescence [9]. Interestingly, a great deal of evidence suggests that crosstalk between BRs and GAs regulates plant growth and development. For instance, physiological studies indicate that co-application of BRs and GAs results in a synergistic increase in hypocotyl elongation of light-grown plants [10]. Knock-down of the expression of SPINDLY (SPY), a negative regulator of GA signalling, enhances BR biosynthesis and enlarges lamina joint bending, which represents a BR-related phenotype [11]. Impressive progress has been made in our understanding of the actions of BRs and GAs, especially in their signalling mechanisms and interactions. Two major molecular mechanisms have been established to illustrate the crosstalk between BR and GA pathways. One mechanism, the so-called ‘signalling model’, depends on the physical interaction between BRASSINAZOLE-RESISTANT 1 (BZR1)/BRI1-EMSSUPPRESSOR 1 (BES1) transcription factors and DELLA proteins [12,13,14]. BZR1 and BES1 play critical roles in the BR-signalling pathway by directly regulating a large number of target genes in the Arabidopsis genome [15,16]. The DELLA protein is the master negative regulator of the GA pathway, and undergoes 26S proteasome pathway-mediated degradation in response to bioactive GA3, resulting in the release of plant growth suppression [17,18]. Thus, DELLA directly interacts with BZR1/BES1 and thereby attenuates their stability and activity, and ultimately restrains plant cell elongation [12,13,14]. In the other mechanism, the so-called ‘synthesis model’, BZR1/BES1 binds directly to the promoters of several GA metabolic genes, including GA20ox, GA3ox and GA2ox, and thereby modulates their expression and promotes cell elongation in both Arabidopsis and rice [19,20]. However, how these two different models are integrated together to help adapt to changing internal and external stimuli remains elusive. Moreover, the mechanisms could act differently in different organs or tissues, and whether the two mechanisms are involved in the regulation of other cellular or developmental programs in addition to cell elongation remains unknown. Recently, a mathematical model was also established to predict the influence of crosstalk between BRs and GAs on plant growth and development [21]. Furthermore, there may also be other mechanisms involved in the crosstalk between BRs and GAs. For example, the chromatin remodelling factor PICKLE functions as a critical node that integrates not only BR and GA pathways, but also light signalling to regulate plant growth and development epigenetically [22]. JUNGBRUNNEN1 (JUB1), a NAC (NAM, ATAF1/2, CUC1/2) transcription factor, regulates the metabolism and signalling of BRs and GAs in Arabidopsis [23] and tomato [24]. Recently, Tang et al. (2018) reported that OsmiR396d also participates in the BR–GA interaction module, thereby coordinating plant architecture in rice [25]. More recently, a new BR-responsive miRNA-target module was identified, that is BR could promote the expression of OsGAMYBL2 by suppressing the level of OsmiR159d. In addition, both OsGSK2, a key negative player in BR signaling, and SLR1, a rice DELLA protein, could directly interact with OsGAMYBL2 to modulate its stability and activity, respectively, thus to coordinate the regulation of BR and GA in plant growth and development [26]. Nevertheless, it remains obscure whether these proposed models and characterised integrators are also suitable for other developmental process, such as seed germination. Indeed, recent findings in maize suggest that BRs and GAs do not interact via a single inclusive pathway [27]. Rather, BR and GA signal transduction and downstream responses are affected by the developmental context.
Seed germination is a pivotal process in the life cycle of higher plants, which is regulated by a series of external and internal clues, including plant hormones [28,29]. It has long been acknowledged that BRs and GAs co-regulate seed germination [30,31,32]. A previous RNA-seq analysis of rice embryos disclosed the phytohormone interaction network during seed germination, and six of the top 10 hub genes belong to GA, BR and ABA pathways [33], highlighting their key roles in rice seed germination. However, the underlying molecular mechanism of their interaction is still obscure. Moreover, rice is not only a monocotyledonous model plant, but also an important food crop as the staple food for more than half of the world’s population. Furthermore, although it is noteworthy that rice contains a relatively small embryo and a dominant endosperm, the rice embryo contains most of the genetic information and plays a decisive role during rice seed germination [34]. Thus, studying the interplay between BRs and GAs in the rice embryo during seed germination is important.
Proteomic analysis is a powerful tool for identifying changes in proteins in response to various internal and external stimuli [35]. Following on from the traditional two-dimensional gel electrophoresis (2-DE) technique, isobaric tags for relative and absolute quantification (iTRAQ)-based proteomic analysis is a new quantitative method that can quantitatively analyse protein abundance in eight samples simultaneously with high confidence and repeatability [36]. It has been used to identify multiple biological processes in various plant species including Arabidopsis [37], rice [38], wheat [39], maize [40]. However, this powerful method has not been applied to investigate how BRs and GAs coregulate rice seed germination, and information regarding the interplay between BRs and GAs during seed germination remains limited. Therefore, in the present study, the iTRAQ proteomic approach was performed to identify the common responsive proteins of BRs and GAs in rice embryos during the most important stage for seed germination (Phase II), thus to help reveal the underlying molecular mechanism of their crosstalk in coordinating rice seed germination.

2. Results

2.1. Generation of the Nip (sd1) Near-Isogenic Line (NIL)

In our previous work, we generated a series of chromosome segment substitutional lines (CSSLs) using the japonica variety Nip as the recipient parent and indica variety 9311 as the donor parent. Among the CSSL lines, the line CSSL120 exhibited a notable semi-dwarf phenotype. After several rounds of backcrossing with Nipponbare and map-based cloning, we finally confirmed that it is the SD1 locus, encoding the GA biosynthesis enzyme GA 20 oxidase (GA20ox), that causes the semi-dwarf phenotype. Alignment of the SD1 coding sequence (CDS) in Nipponbare and 9311 indicated the existence of four single-nucleotide polymorphisms (SNPs), among which one SNP (654) did not change an amino acid, two SNPs (299 and 1019) resulted in an amino acid substitution (Glu to Gly and Gln to Arg, respectively), and the remaining SNP (1026) changed a Tyr residue to a stop codon (Figure S1). Thus, the SD1 allele in 9311 encodes a mutated and defective form of GA20ox, resulting in the GA-deficient endogenous condition. One progeny of the backcross line between CSSL 120 and Nip was genotyped and selected as the NIL of SD1 in the Nip background, designated as Nip (sd1) or sd1, and used in subsequent seed germination and proteomic analysis (Figure S2). In general, sd1 exhibited a semi-dwarf phenotype in terms of plant height (Figure 1A,D), but its other agronomic traits, including grain size, tiller number, length and width of flag leaf, were similar to that of the Nip (SD1) control (Figure 1B,C,E–I). Therefore, seeds of Nip (sd1) were ideal materials for investigating how GA affects rice germination.

2.2. Seed Germination Assays with Either Brassinosteroid (BR) or Gibberellin (GA) Biosynthesis Blocked

Shoot and root length of germinated seeds were measured every 12 h, from 48 to 96 hours after imbibition (HAI). The results indicate that both BRZ treatment and SD1 mutation led to a significant reduction in shoot length during most tested stages (Figure 2A). Regarding root length, a significant change in BR-deficient seeds was observed after 84 HAI, but a change was only evident after 96 HAI for GA-deficient seeds (Figure 2B). These results suggest that both BR and GA play essential roles during rice seed germination.

2.3. Expression Changes in Embryo Proteins in Response to BR or GA Deficiency during Seed Germination

Since proper translation in early stage of germination, especially phase II, is a prerequisite for germination, the embryos of germinated seeds after 36 HAI were used for proteomic analysis. Here, the non-gel-based iTRAQ technique was applied to identify common proteins in the response to both BR and GA deficiency. The embryos were isolated for proteomic assays, and total proteomes from three different origins, Nip/BRZ, Nip/DMSO and sd1/DMSO, were labelled by iTRAQ and quantified to identify differentially expressed proteins. Specifically, protein abundance was compared between Nip/BRZ and Nip/DMSO, and between sd1/DMSO and Nip/DMSO, respectively. Based on MS/MS analysis and comparison of the results, 45 proteins were identified as common targets in responses to both BR and GA deficiency (Figure 3). Among these, expression of 29 proteins was changed consistently in response to both BR and GA deficiency, of which 15 were down-regulated (Table 1) and 14 were up-regulated (Table 2). Among the 29 proteins, five (glutelin B2, prolamin, glycine-rich protein, pectin esterase-like protein, and acyl-desaturase) displayed distinct expression changes (>1.5-fold) in response to both BRZ and sd1 mutation (Figure S3). Meanwhile, 16 proteins exhibited differences in expression pattern in BRZ-treated Nip and mock-treated sd1 rice (Table 3). Thus, the 45 common differentially expressed proteins responsive to both BR and GA deficiency represented candidates for investigation of crosstalk in the coordination of rice seed germination.

2.4. Bioinformaticf Analysis of Common Responsive Proteins Co-Regulated by BR and GA

To gain knowledge of differentially expressed proteins common to both BR and GA in the regulation of seed germination, a number of bioinformatic analyses were performed, including hierarchical clustering, which is widely applied for grouping proteins exhibiting similar expression changes in omics studies. Accordingly, the 45 common proteins were classified into clusters I to IV (Figure 4A). Most proteins in cluster II and III shared concordant expression patterns; protein abundance in cluster II was markedly down-regulated and that in cluster III was up-regulated in both BRZ-treated NIP and the mock-treated sd1 mutant, suggesting that BR and GA deficiency resulted in similar expression changes for many common targets. Using the Clusters of Orthologous Groups (COG) database, the 45 proteins were further divided into different functional categories (Figure 4B), among which the ‘Amino acid transport and metabolism’ category included five proteins, and the ‘Energy production and conversion’ and ‘General function prediction only’ groups each contained four proteins.
The results of molecular function analysis showed that all test proteins could be classified into five functional groups (catalytic activity, binding activity, structural molecule activity, transporter activity, and antioxidant activity). Among these, almost 70% of total proteins belonged to the ‘catalytic activity’ group, and the ‘binding activity’ group was the next most enriched (Figure 5A). Based on GO analysis of biological processes, more than half of test proteins were classified into the ‘metabolic process’ group. The second largest group was the ‘cellular process’ group, containing more than 25% of total proteins. Additionally, a few proteins belonged to five other groups (localisation, cellular component biogenesis, developmental process, biological regulation, and response to stimulus; Figure 5B). When the proteins were classified according to cellular component, the ‘cell part’ group contained the largest number of proteins, followed by the ‘organelle’ and ‘macromolecular complex’ groups (Figure 5C). Finally, the results of GO analysis of the protein class showed that these proteins could be grouped into 10 subfamilies, among which oxidoreductases accounted for the largest proportion (31.2%), followed by hydrolase (15.6%), transferase (9.4%), nucleic acid binding (9.4%), and lyase (9.4%; Figure 5D).
Dissecting the complex network of protein-protein interactions (PPI) could facilitate a better understanding of the molecular mechanism of the crosstalk between BR and GA in seed germination. Totally 19 proteins were included in the constructed network (Figure 6), which was made with a medium confidence cutoff (0.4) using the k-means clustering method. Interestingly, 17 proteins (89.5%) had a high degree of connectivity (more than 2 connections) and were centered on three main node (protein) clusters. Proteins with the highest connectivity are mainly involved in the metabolic process, such as alcohol dehydrogenase, ATP synthase and Type IIA topoisomerase, suggesting the metabolic process is the critical regulation node of seed germination by BR and GA.

2.5. Validation of Several Representative Genes by qRT-PCR

To further validate and complement our proteomic results, qRT-PCR analysis was performed to examine the transcriptions of 15 selected proteins co-regulated by both BR and GA. The selected proteins could be further classified into two different groups; group one displaying consistent protein pattern changes, and group two displaying opposite protein pattern changes (Figure 7A,B). In group one, the abundance of five proteins was elevated in response to both BR and GA deficiency, while three were down-regulated. In general, changes in transcript levels and protein abundance among targets were more consistent in the sd1 mutant than in the BRZ-treated group. Interestingly, the transcription and translation of three targets, histone H1-like protein, pectin esterase-like protein, and acyl-desaturase, were all induced by both BR and GA deficiency (Figure 7B). In group two, although changes in protein abundance were opposite for all seven selected targets in the sd1 mutant and BRZ-treated Nip, changes in transcription of some targets were consistent. For example, transcription of embryo-specific protein, DREPP protein, and DUF248 methyltransferase was significantly increased, while that of a hypothetical protein was notably decreased in both rice materials (Figure 7D). Furthermore, high isoelectric point (pI) alpha-glucosidase and acyl desaturase displayed consistent changes in transcription and protein abundance (Figure 7C,D).

2.6. Generation and Germination Analysis of grp Mutant

The binary vector targeting GRP gene was used for rice transformation. A total of 12 positive transgenic plants were obtained in the T0 generation. After sequencing the target region in all positive transgenic plants, one transgenic plant with homozygous mutation was obtained. 1-bp deletion (T) occurred at the 3th base from the protospacer adjacent motif (PAM) site (Figure 8A). Thus, the homozygous mutation line was designated as grp, whose mature seeds were used for the following germination analysis. In general, germination rate of grp mutant was significantly lower than the wild-type control albeit with or without the treatment of BR and GA (Figure 8B,C). Moreover, the lengths of both shoots and roots of the germinated grp seeds were also shorter than the control (Figure 8D,E), suggesting that GRP mutation has an adverse effect on rice seed germination. Therefore, the identified BR and GA responsive proteins, including GRP, should be involved in BR- and GA-regulated rice seed germination process.

3. Discussion

BRs and GAs engage in direct interactions and coordinate a series of plant growth and development events, including seed germination. Nevertheless, it remains unclear how BR and GA coordinate seed germination and what are the common downstream targets that mediate seed germination. In the present study, we identified BR- and GA-responsive proteins in rice embryos during the early stages of seed germination using an iTRAQ proteomic approach and identified downstream targets of BRs and GAs that are intimately involved in the regulation of seed germination. In general, more differentially expressed proteins were identified in mutant sd1 than that in the BRZ-treated NIP samples. It is possible that GA deficiency caused a more dramatic expression change of downstream regulated proteins, which is consistent with the fact that GA is a major determinant of seed germination. Another possibility is that although BRZ treatment could mimic the BR-deficient conditions, it still could not absolutely exhaust the endogenous BRs. Therefore, a relatively mild effect was caused by pharmaceutical treatment, which might result in a smaller number of changed proteins. A Venn diagram illustrated that a total of 45 common responsive proteins were identified, of which ~64% exhibited coincident changes in expression patterns under BR- or GA-deficient germination conditions, suggesting that crosstalk between BR and GA pathways controls seed germination, at least in part, by modulating the same set of key proteins. Interestingly, GO analysis indicated that most of the common targets were proteins with catalytic activity that participate in various metabolic processes, consistent with the fact that a series of germination events occur in phase II, the most important stage for seed germination, including metabolism reactivation, coleoptile elongation, and mobilisation of reserves [41].
Among the common targets, five proteins (glutelin B2, prolamin, glycine-rich protein, pectin esterase, and acyl-desaturase) shared remarkably consistent expression changes (>1.5-fold or <0.67-fold) in both BR- and GA-deficient conditions (Figure S3). For instance, expression of both glutelin B2 and prolamin, two storage proteins in rice seeds, declined significantly when endogenous BR or GA biosynthesis was attenuated, suggesting the decreased accumulation of storage proteins in embryos might contribute partly to the defect in seed germination. This is consistent with previous reports that suppression of seed storage proteins by RNAi, including both glutelins and prolamins, delayed rice seed germination [42]. A recent proteomic study also indicated that the abundance of several glutelin storage proteins was down-regulated during rice seed germination under salt stress [43]. Moreover, abundance of the other three proteins, including glycine-rich protein, pectin esterase and acyl desaturase, was enhanced. Other studies have emphasised a potential correlation between glycine-rich proteins and seed germination. For example, the quantitative trait locus (QTL) qLTG3-1 controlling low-temperature tolerance at the germination stage was cloned and found to encode a hybrid glycine-rich protein (HyGRP). Strong expression of qLTG3-1 in embryos enhances rice seed germination under low-temperature conditions by weakening the tissue covering the embryo [44]. To validate our proteomic result and further dissect the function of GRP in seed germination, a grp mutant was generated using the clustered regularly interspaced short palindromic repeats-associated protein 9 CRISPR/Cas9 system. Germination analysis showed that GRP mutation remarkably suppressed rice seed germination, which could not be fully recovered by GA or BL treatment (Figure 8B–E), suggesting the critical roles of GRP in phytohormone-regulated seed germination. Further diversification analysis of the HyGRP gene family in five monocots, namely rice, barley, brachypodium, maize and sorghum, indicated that the HyGRP gene family was formed in the ancestral genome before the divergence of these species. The conserved collinearity of chromosomal regions around qLTG3-1 and its orthologs implies that their functions in seed germination are conserved [45]. In addition, studies indicated that mutation of PME48, a member of the pectin esterase family, would cause a significant delay in imbibition and germination of pollen grains in Arabidopsis [46]. On the other hand, enhancing the activity of pectin esterase by direct protein-protein interaction could also improve seed germination performance [47]. However, the function of acyl desaturase in seed germination have not yet been reported yet and require further investigation.
In addition, several other proteins identified in the present proteomic study also played key roles in seed germination. For example, protein abundance of a plant SKP1-like family protein, component of the SCF (SKP-cullin-F-box protein) E3 ligase complexes that ubiquitinate target proteins, is increased in response to both BR and GA deficiency (Table 2). In Arabidopsis, abiotic stress could induce the expression of Arabidopsis SKP1-like protein 13 (ASK13); in turn, ASK13 positively regulates seed germination and seedling growth under abiotic stress [48]. Another study also demonstrated that transcript abundance of GsSKP21 was induced under the treatment of alkali and salt stresses in Glycine soja. Overexpression of GsSKP21 in Arabidopsis results in decreased ABA sensitivity and promoted seed germination [49]. Moreover, the accumulation of glutamate decarboxylase (GAD) and alcohol dehydrogenase (ADH) are also promoted in both sd1 mutant and BRZ-treated Nip (Table 2). GAD functions in catalyzing the unidirectional decarboxylation of glutamate to form γ-aminobutyric acid (GABA). Overexpression of GAD could promote the amounts of GABA in Arabidopsis seeds [50]. GABA homeostasis is involved in the regulation of seed germination in both soybean [51] and fava bean [52], and pollen germination in Picea wilsonii [53]. As to ADH, its deficient rice mutant rad had lower amounts of the α-amylases, resulting in less starch degradation and lower glucose concentrations, subsequent a delayed germination in partially oxygenated water [54]. A rice genotype Khao Hlan On (KHO), which is favorable for direct seedling, accumulates more alcohol dehydrogenase ADH1 and exhibits faster germination than the control under the anaerobic germination condition [55]. Interestingly, the content of Peptidase S10, another identified protein in this study, was decreased in response to both BR and GA deficiency (Table 1). It was reported that Site-2 Protease (S2P), a membrane-associated transcription factor peptidase, desensitizes ABA signaling during seed germination through regulating the expression of negative regulators of ABA signaling, and therefore promotes seed germination [56].
Furthermore, the potential protein association network constructed by STRING analysis highlighted the fact that proteins involved in the metabolic process of seed germination are regulated by BR and GA. For example, alcohol dehydrogenase, with the highest degree of connectivity in the regulatory network, was induced in response to both BR and GA deficiency. A number of reports have suggested that various dehydrogenases are involved in regulating seed germination in sweet sorghum [57], Arabidopsis [58] and the legume Medicago truncatula [59]. Topoisomerase, another key factor in the centre of the protein association network, can enhance seed germination under salt stress conditions in both Arabidopsis [60] and tobacco [61]. The qRT-PCR approach is widely applied to validate proteomic data at the transcriptional level. However, changes in transcript level and protein abundance are not always consistent. Herein, 15 differentially expressed proteins were selected for validation, and only six displayed concordant changes at both transcriptional and translational levels (Figure 7A–D). Such inconsistencies are often not unusual. For instance, transcriptome and proteome analyses were integrated to study the potential mechanism underlying Ogura cytoplasmic male sterility (OguCMS) in cabbage [62]. Using the same set of anther material for RNA-seq and iTRAQ analysis, 1323 differentially expressed genes and 833 differentially abundant proteins were identified. However, only 92 targets overlapped in the transcriptomic and proteomic data, and even more strikingly, only 22 out of the 92 targets shared consistent expression patterns. Similar results have also been reported for both plants and animals, implying a weak correlation between gene transcription and protein abundance [63,64,65,66,67]. Such discrepancies in transcriptomic and proteomic results may be due to the different regulatory mechanisms for mRNAs and proteins. Moreover, genes with distinct functions can be controlled by direct or delayed translation modules, hence delayed translation may also partly contribute to these weak correlations.

4. Materials and Methods

4.1. Plant Materials

In this study, two different rice materials were used for proteomic experiment; japonica cultivar Nipponbare (Nip) that has a functional SD1 (SEMIDWARF1) gene (Nip (SD1)), and a near-isogenic line (NIL) of SD1 locus generated by introduction of a mutated sd1 allele from indica cultivar 9311 into the Nip genetic background, designated Nip (sd1). SD1, the gene of the Green Revolution, also known as OsGA20ox2, encodes a key gibberellin biosynthesis enzyme. Sd1, the null allele of SD1, produces a semi-dwarf phenotype and was selected during the rice Green Revolution [68,69]. In addition, another japonica cultivar Zhonghua 11 (ZH11) was used for Agrobacterium-mediated transformation.

4.2. Plant Growth and Seed Germination Assay

Both Nip (SD1) and Nip (sd1) were grown in the same paddy field at the Wenhui Road Campus, Yangzhou University under identical climatic conditions, and mature seeds were collected for germination assays. Since the seeds of BR-deficient plants were abnormal, and homozygous mutants are often sterile [70], we used brassinazole (BRZ; a commonly used inhibitor of BR biosynthesis) to treat Nip seeds to mimic the BR-deficient germination conditions. Because BRZ was dissolved in dimethyl sulphoxide (DMSO), Nip seeds were imbibed in Milli-Q water supplemented with BRZ or DMSO solution in darkness, while seeds of Nip (sd1) were also imbibed in DMSO mock solution. In brief, 100 rice seeds were manually dehulled, sterilised with 70% ethanol, and rinsed twice with sterile water. Seeds were then imbibed with 1 µM BRZor mock solutions at 26 °C in darkness. Embryos of seeds germinated for 36 h were separated and quickly frozen in liquid nitrogen for proteomic analysis. The lengths of shoots and roots of 30 germinated seeds were measured after germination for 60, 72, 84 and 96 h. As to the seed germination assay of grp mutant, the seeds were imbibed in Milli-Q water supplemented with 1 µM BL, 10 µM GA or the mock solution, respectively. Three biological replicates were performed for all experiments.

4.3. Protein Extraction

Dithiothreitol (DTT; 10 mM) was added to embryo samples and grinding was performed in liquid nitrogen. The suspension was sonicated for 15 min, centrifuged for 15 min (30,000× g, 4 °C), and the supernatant was mixed thoroughly with a 5× volume of chilled acetone containing 10% (v/v) trichloroacetic acid (TCA) and incubated at ×20 °C overnight. After centrifugation, the precipitate was washed three times with chilled acetone, air-dried, and dissolved in lysis buffer (7 M urea, 2 M thiourea, 4% NP40, 20 mM Tris-HCl, pH 8.0–8.5). The suspension was sonicated and centrifuged as above, the supernatant was transferred to a new tube, and 10 mM DTT was added and incubated at 56 °C for 1 h. Next, 55 mM iodoacetamide (IAM) was added and incubated for 1 h in darkness, and the supernatant was mixed thoroughly with chilled acetone (2 h, −20 °C). After centrifugation, the pellet was air-dried and dissolved in 500 µL of 0.5 M triethylammonium bicarbonate (TEAB; Applied Biosystems, Milan, Italy), sonicated for 15 min (200 W), and the supernatant was collected for proteomic assays.

4.4. iTRAQ Labelling and Liquid Chromatography-Electrospray Ionisation-Tandem Mass Spectroscopy (LC-ESI-MS/MS)

In brief, 100 µg of total proteins was digested with Gold-grade Trypsin, dried, reconstituted, and labelled with 8-plex iTRAQ reagent (Applied Biosystems, Foster City, CA, USA). Labelled peptide mixtures were pooled into 20 fractions, desalted, and vacuum-dried. Each fraction was separated and data were acquired using a TripleTOF 5600 System (AB SCIEX, Concord, ON, Canada). A detailed description of iTRAQ labelling and LC-ESI-MS/MS analysis is described in our previous publication [71].

4.5. Proteomic Data Analysis

The collected raw data were processed using Proteome Discoverer 1.2 (PD 1.2, Thermo Scientific, San Jose, CA, USA), and the Mascot search engine was used for protein identification. Only peptides with a 95% confidence interval were regarded as positively identified, at least one unique peptide was required for every confident protein identification, and for quantitation, at least two unique peptides were required. Quantitative protein ratios were weighted and normalised by median ratio in Mascot. Only ratios with a fold change >1.2 and p-values < 0.05 were considered differentially abundant.

4.6. Bioinformatic Analysis

Cluster software (version 3.0) was used for hierarchical clustering analysis of protein expression patterns [72], and functional annotation of proteins was conducted using the Blast2GO program to search against the non-redundant protein database (NR; NCBI). Gene ontology (GO) was used for gene function classification and description of genes or gene product attributes. Using the PANTHER database, proteins were analysed with four sets of ontologies: molecular function, biological process, cellular component, and protein class [73]. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) version 10.0 was then used for interactome analysis of identified proteins [74]. Using various prediction methods, including experiments, co-expression, databases, neighbourhood, co-occurrence and text-mining, information on potential protein–protein interactions was retrieved with a confidence cut-off of 0.4.

4.7. Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Analysis

Total RNA was isolated from the collected embryos of germinated seeds using an RNeasy Plant Mini Kit (Qiagen, Hilden, Germany), and RNA was digested with DNase I (Qiagen) and reverse-transcribed with a SuperScript first-strand synthesis system (Invitrogen, Van Allen Way Carlsbad, CA, USA). After first-strand cDNA synthesis, quantitative real-time polymerase chain reaction (qRT-PCR) was performed using SYBR Premix Ex Taq (TaKaRa, Dalian, Liaoning, China) and an ABI PRISMTM 7700 sequence detector system (Applied Biosystems). The relative transcript abundance was estimated using the 2−ΔΔCt method, and ubiquitin conjugase (UBC) was selected as the reference gene. Primer information is provided in Table S1 in Supplementary Materials.

4.8. Plasmid Construction and Rice Transformation

The target site of glycine-rich protein (GRP) was designed for knock out of the gene using CRISPR/Cas9 system. The sgRNA-Cas9 expressing plasmid was constructed following the method described by Wang et al. (2015) [75]. Subsequently, the construct was transformed into Agrobacterium tumefaciens EHA105 and then introduced into primary calli derived from mature seeds of japonica rice ZH11 via Agrobacterium-mediated transformation [76]. The target sequence of GRP and primers for sequencing are provided in Table S1.

5. Conclusions

Understanding the mechanisms of seed germination is beneficial for breeding elite rice varieties with high seed vigor. Our proteome analysis identified a number of common responsive proteins of both BRs and GAs during rice seed germination, most of which were involved in the metabolic process and exhibited concordant changes in response to BR and GA deficiency. GO and STRING analyses revealed the functions and possible roles of these proteins in mediating BR- and GA-regulated seed germination. Further functional analysis of one representative protein GRP demonstrated its critical role in normal seed germination. In brief, our study not only provides critical clues for validating BR–GA crosstalk models during rice seed germination, but also helps to optimise the underlying molecular regulatory networks.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/19/11/3460/s1, Figure S1: Alignment of SD1/OsGA20ox2 CDS from Nipponbare (Nip; top) and 9311 (bottom). Four SNPs were detected (shown in red), Figure S2: Genotyping the sd1 mutant, Figure S3: Five selected responsive proteins displaying consistent prominent changes in protein abundance (fold change ratios ≥1.5 or ≤0.67) in response to both BR and GA deficiency, Table S1: Primers and SG sequence used in this study.

Author Contributions

Conceptualization and Funding acquisition, Q.-F.L. and Q.-Q.L.; Investigation, Q.-F.L., J.-D.W., M.X., K.W.; Methodology, J.-D.W., M.X., P.Z., and L.-C.H.; Resources, C.-Q.Z.; Validation, Q.-F.L. and X.-L.F.; Writing—original draft, Q.-F.L.; Writing—review and editing, Q.-Q.L.

Funding

This research was supported by the National Natural Science Foundation of China [31561143008, 31601275, 31771745], the Jiangsu Provincial Department of Education of China [17KJA210001], Jiangsu PAPD, “333” and Innovative and Entrepreneurial Talent projects.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vert, G.; Chory, J. Crosstalk in cellular signaling: Background noise or the real thing? Dev. Cell 2011, 21, 985–991. [Google Scholar] [CrossRef] [PubMed]
  2. Lau, O.S.; Deng, X.W. Plant hormone signaling lightens up: Integrators of light and hormones. Curr. Opin. Plant Biol. 2010, 13, 571–577. [Google Scholar] [CrossRef] [PubMed]
  3. Santner, A.; Calderon-Villalobos, L.I.; Estelle, M. Plant hormones are versatile chemical regulators of plant growth. Nat. Chem. Biol. 2009, 5, 301–307. [Google Scholar] [CrossRef] [PubMed]
  4. Santner, A.; Estelle, M. Recent advances and emerging trends in plant hormone signalling. Nature 2009, 459, 1071–1078. [Google Scholar] [CrossRef] [PubMed]
  5. Vanstraelen, M.; Benkova, E. Hormonal interactions in the regulation of plant development. Annu. Rev. Cell Dev. Biol. 2012, 28, 463–487. [Google Scholar] [CrossRef] [PubMed]
  6. Depuydt, S.; Hardtke, C.S. Hormone signalling crosstalk in plant growth regulation. Curr. Biol. 2011, 21, R365–R373. [Google Scholar] [CrossRef] [PubMed]
  7. Howe, G.A.; Major, I.T.; Koo, A.J. Modularity in jasmonate signaling for multistress resilience. Annu. Rev. Plant Biol. 2018, 69, 387–415. [Google Scholar] [CrossRef] [PubMed]
  8. Verma, V.; Ravindran, P.; Kumar, P.P. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016, 16, 86. [Google Scholar] [CrossRef] [PubMed]
  9. Li, Q.F.; He, J.X. Mechanisms of signaling crosstalk between brassinosteroids and gibberellins. Plant Signal. Behav. 2013, 8, e24686. [Google Scholar] [CrossRef] [PubMed]
  10. Tanaka, K.; Nakamura, Y.; Asami, T.; Yoshida, S.; Matsuo, T.; Okamoto, S. Physiological roles of brassinosteroids in early growth of Arabidopsis: Brassinosteroids have a synergistic relationship with gibberellin as well as auxin in light-grown hypocotyl elongation. J. Plant Growth Regul. 2003, 22, 259–271. [Google Scholar] [CrossRef]
  11. Shimada, A.; Ueguchi-Tanaka, M.; Sakamoto, T.; Fujioka, S.; Takatsuto, S.; Yoshida, S.; Sazuka, T.; Ashikari, M.; Matsuoka, M. The rice SPINDLY gene functions as a negative regulator of gibberellin signaling by controlling the suppressive function of the DELLA protein, SLR1, and modulating brassinosteroid synthesis. Plant J. 2006, 48, 390–402. [Google Scholar] [CrossRef] [PubMed]
  12. Bai, M.Y.; Shang, J.X.; Oh, E.; Fan, M.; Bai, Y.; Zentella, R.; Sun, T.P.; Wang, Z.Y. Brassinosteroid, gibberellin and phytochrome impinge on a common transcription module in Arabidopsis. Nat. Cell Biol. 2012, 14, 810–817. [Google Scholar] [CrossRef] [PubMed]
  13. Gallego-Bartolome, J.; Minguet, E.G.; Grau-Enguix, F.; Abbas, M.; Locascio, A.; Thomas, S.G.; Alabadi, D.; Blazquez, M.A. Molecular mechanism for the interaction between gibberellin and brassinosteroid signaling pathways in Arabidopsis. Proc. Natl. Acad. Sci. USA 2012, 109, 13446–13451. [Google Scholar] [CrossRef] [PubMed]
  14. Li, Q.F.; Wang, C.; Jiang, L.; Li, S.; Sun, S.S.; He, J.X. An interaction between BZR1 and DELLAs mediates direct signaling crosstalk between brassinosteroids and gibberellins in Arabidopsis. Sci Signal. 2012, 5, ra72. [Google Scholar] [CrossRef] [PubMed]
  15. Sun, Y.; Fan, X.Y.; Cao, D.M.; Tang, W.; He, K.; Zhu, J.Y.; He, J.X.; Bai, M.Y.; Zhu, S.; Oh, E.; et al. Integration of brassinosteroid signal transduction with the transcription network for plant growth regulation in Arabidopsis. Dev. Cell 2010, 19, 765–777. [Google Scholar] [CrossRef] [PubMed]
  16. Yu, X.; Li, L.; Zola, J.; Aluru, M.; Ye, H.; Foudree, A.; Guo, H.; Anderson, S.; Aluru, S.; Liu, P.; et al. A brassinosteroid transcriptional network revealed by genome-wide identification of BESI target genes in Arabidopsis thaliana. Plant J. 2011, 65, 634–646. [Google Scholar] [CrossRef] [PubMed]
  17. Sun, T.P. The molecular mechanism and evolution of the GA–GID1–DELLA signaling module in plants. Curr. Biol. 2011, 21, R338–R345. [Google Scholar] [CrossRef] [PubMed]
  18. Gao, X.H.; Xiao, S.L.; Yao, Q.F.; Wang, Y.J.; Fu, X.D. An updated GA signaling ‘relief of repression’ regulatory model. Mol. Plant 2011, 4, 601–606. [Google Scholar] [CrossRef] [PubMed]
  19. Tong, H.; Xiao, Y.; Liu, D.; Gao, S.; Liu, L.; Yin, Y.; Jin, Y.; Qian, Q.; Chu, C. Brassinosteroid regulates cell elongation by modulating gibberellin metabolism in rice. Plant Cell 2014, 26, 4376–4393. [Google Scholar] [CrossRef] [PubMed]
  20. Unterholzner, S.J.; Rozhon, W.; Papacek, M. Brassinosteroids are master regulators of gibberellin biosynthesis in Arabidopsis. Plant Cell 2015, 27, 2261–2272. [Google Scholar] [CrossRef] [PubMed]
  21. Allen, H.R.; Ptashnyk, M. Mathematical modelling and analysis of the brassinosteroid and gibberellin signalling pathways and their interactions. J. Theor. Biol. 2017, 432, 109–131. [Google Scholar] [CrossRef] [PubMed]
  22. Zhang, D.; Jing, Y.; Jiang, Z.; Lin, R. The chromatin-remodeling factor PICKLE integrates brassinosteroid and gibberellin signaling during skotomorphogenic growth in Arabidopsis. Plant Cell 2014, 26, 2472–2485. [Google Scholar] [CrossRef] [PubMed]
  23. Shahnejat-Bushehri, S.; Tarkowska, D.; Sakuraba, Y.; Balazadeh, S. Arabidopsis NAC transcription factor JUB1 regulates GA/BR metabolism and signalling. Nat. Plants 2016, 2, 16013. [Google Scholar] [CrossRef] [PubMed]
  24. Shahnejat-Bushehri, S.; Allu, A.D.; Mehterov, N.; Thirumalaikumar, V.P.; Alseekh, S.; Fernie, A.R.; Mueller-Roeber, B.; Balazadeh, S. Arabidopsis NAC transcription factor JUNGBRUNNEN1 exerts conserved control over gibberellin and brassinosteroid metabolism and signaling genes in tomato. Front. Plant Sci. 2017, 8, 214. [Google Scholar] [CrossRef] [PubMed]
  25. Tang, Y.; Liu, H.; Guo, S.; Wang, B.; Li, Z.; Chong, K.; Xu, Y. OsmiR396d affects gibberellin and brassinosteroid signaling to regulate plant architecture in rice. Plant Physiol. 2018, 176, 946–959. [Google Scholar] [CrossRef] [PubMed]
  26. Gao, J.; Chen, H.; Yang, H.; He, Y.; Tian, Z.; Li, J. A brassinosteroid responsive miRNA-target module regulates gibberellin biosynthesis and plant development. New Phytol. 2018, 220, 488–501. [Google Scholar] [CrossRef] [PubMed]
  27. Best, N.B.; Hartwig, T. nana plant2 encodes a maize ortholog of the Arabidopsis brassinosteroid biosynthesis gene Dwarf1, identifying developmental interactions between brassinosteroids and gibberellins. Plant Physiol. 2016, 171, 2633–2647. [Google Scholar] [PubMed]
  28. Wojtyla, Ł.; Lechowska, K.; Kubala, S.; Garnczarska, M. Different modes of hydrogen peroxide action during seed germination. Front. Plant Sci. 2016, 7, 66. [Google Scholar] [CrossRef] [PubMed]
  29. Vishal, B.; Kumar, P.P. Regulation of seed germination and abiotic stresses by gibberellins and abscisic acid. Front. Plant Sci. 2018, 9, 838. [Google Scholar] [CrossRef] [PubMed]
  30. Chen, J.G.; Pandey, S.; Huang, J.; Alonso, J.M.; Ecker, J.R.; Assmann, S.M.; Jones, A.M. GCR1 can act independently of heterotrimeric G-protein in response to brassinosteroids and gibberellins in Arabidopsis seed germination. Plant Physiol. 2004, 135, 907–915. [Google Scholar] [CrossRef] [PubMed]
  31. Steber, C.M.; McCourt, P. A role for brassinosteroids in germination in Arabidopsis. Plant Physiol. 2001, 125, 763–769. [Google Scholar] [CrossRef] [PubMed]
  32. Shu, K.; Liu, X.D.; Xie, Q.; He, Z.H. Two Faces of One Seed: Hormonal Regulation of Dormancy and Germination. Mol. Plant 2016, 9, 34–45. [Google Scholar] [CrossRef] [PubMed]
  33. Wei, T.; He, Z.; Tan, X.; Liu, X.; Yuan, X.; Luo, Y.; Hu, S. An integrated RNA-Seq and network study reveals a complex regulation process of rice embryo during seed germination. Biochem. Biophys. Res. Commun. 2015, 464, 176–181. [Google Scholar] [CrossRef] [PubMed]
  34. Han, C.; Zhen, S.; Zhu, G.; Bian, Y.; Yan, Y. Comparative metabolome analysis of wheat embryo and endosperm reveals the dynamic changes of metabolites during seed germination. Plant Physiol. Biochem. 2017, 115, 320–327. [Google Scholar] [CrossRef] [PubMed]
  35. Domzalska, L.; Kedracka-Krok, S.; Jankowska, U.; Grzyb, M.; Sobczak, M.; Rybczynski, J.J.; Mikula, A. Proteomic analysis of stipe explants reveals differentially expressed proteins involved in early direct somatic embryogenesis of the tree fern Cyathea delgadii Sternb. Plant Sci. 2017, 258, 61–76. [Google Scholar] [CrossRef] [PubMed]
  36. Evans, C.; Noirel, J.; Ow, S.Y.; Salim, M.; Pereira-Medrano, A.G.; Couto, N.; Pandhal, J.; Smith, D.; Pham, T.K.; Karunakaran, E.; et al. An insight into iTRAQ: Where do we stand now? Anal. Bioanal. Chem. 2012, 404, 1011–1127. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, F.X.; Luo, Y.M.; Ye, Z.Q.; Cao, X.; Liang, J.N.; Wang, Q.; Wu, Y.; Wu, J.H.; Wang, H.Y.; Zhang, M.; et al. iTRAQ-based proteomics analysis of autophagy-mediated immune responses against the vascular fungal pathogen Verticillium dahliae in Arabidopsis. Autophagy 2018, 14, 598–618. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, J.; Islam, F.; Li, L.; Long, M.; Yang, C.; Jin, X.; Ali, B.; Mao, B. Complementary RNA-sequencing based transcriptomics and iTRAQ proteomics reveal the mechanism of the alleviation of quinclorac stress by salicylic acid in Oryza sativa ssp. japonica. Int. J. Mol. Sci. 2017, 18, 1975. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, S.; Chen, J.; Hou, F.; Feng, Y.; Zhang, R. iTRAQ-based quantitative proteomic analysis reveals the lateral meristem developmental mechanism for branched spike development in tetraploid wheat (Triticum turgidum L.). BMC Genom. 2018, 19, 228. [Google Scholar] [CrossRef] [PubMed]
  40. Bu, T.T.; Shen, J.; Chao, Q.; Shen, Z.; Yan, Z.; Zheng, H.Y.; Wang, B.C. Dynamic N-glycoproteome analysis of maize seedling leaves during de-etiolation using Concanavalin A lectin affinity chromatography and a nano-LC-MS/MS-based iTRAQ approach. Plant Cell Rep. 2017, 36, 1943–1958. [Google Scholar] [CrossRef] [PubMed]
  41. Han, C.; Yang, P. Studies on the molecular mechanisms of seed germination. Proteomics 2015, 15, 1671–1679. [Google Scholar] [CrossRef] [PubMed]
  42. Cho, K.; Lee, H.J.; Jo, Y.M.; Lim, S.H.; Rakwal, R.; Lee, J.Y.; Kim, Y.M. RNA interference-mediated simultaneous suppression of seed storage proteins in rice grains. Front. Plant Sci. 2016, 7, 1624. [Google Scholar] [CrossRef] [PubMed]
  43. Xu, E.; Chen, M.; He, H.; Zhan, C.; Cheng, Y.; Zhang, H.; Wang, Z. Proteomic analysis reveals proteins involved in seed imbibition under salt stress in rice. Front. Plant Sci. 2016, 7, 2006. [Google Scholar] [CrossRef] [PubMed]
  44. Fujino, K.; Sekiguchi, H.; Matsuda, Y.; Sugimoto, K.; Ono, K.; Yano, M. Molecular identification of a major quantitative trait locus, qLTG3-1, controlling low-temperature germinability in rice. Proc. Natl. Acad. Sci. USA 2008, 105, 12623–12628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Fujino, K.; Obara, M.; Sato, K. Diversification of the plant-specific hybrid glycine-rich protein (HyGRP) genes in cereals. Front. Plant Sci. 2014, 5, 489. [Google Scholar] [CrossRef] [PubMed]
  46. Leroux, C.; Bouton, S.; Kiefer-Meyer, M.C.; Fabrice, T.N.; Mareck, A.; Guénin, S.; Fournet, F.; Ringli, C.; Pelloux, J.; Driouich, A.; et al. PECTIN METHYLESTERASE48 is involved in Arabidopsis pollen grain germination. Plant Physiol. 2015, 167, 367–380. [Google Scholar] [CrossRef] [PubMed]
  47. Zúñiga-Sánchez, E.; Soriano, D.; Martínez-Barajas, E.; Orozco-Segovia, A.; Gamboa-deBuen, A. BIIDXI, the At4g32460 DUF642 gene, is involved in pectin methyl esterase regulation during Arabidopsis thaliana seed germination and plant development. BMC Plant Biol. 2014, 14, 338. [Google Scholar] [CrossRef] [PubMed]
  48. Rao, V.; Petla, B.P.; Verma, P.; Salvi, P.; Kamble, N.U.; Ghosh, S.; Kaur, H.; Saxena, S.C.; Majee, M. Arabidopsis SKP1-like protein13 (ASK13) positively regulates seed germination and seedling growth under abiotic stress. J. Exp. Bot. 2018, 69, 3899–3915. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, A.; Yu, Y.; Duan, X.; Sun, X.; Duanmu, H.; Zhu, Y. GsSKP21, a Glycine soja S-phase kinase-associated protein, mediates the regulation of plant alkaline tolerance and ABA sensitivity. Plant Mol. Biol. 2015, 87, 111–124. [Google Scholar] [CrossRef] [PubMed]
  50. Fait, A.; Nesi, A.N.; Angelovici, R.; Lehmann, M.; Pham, P.A.; Song, L.; Haslam, R.P.; Napier, J.A.; Galili, G.; Fernie, A.R. Targeted enhancement of glutamate-to-γ-aminobutyrate conversion in Arabidopsis seeds affects carbon-nitrogen balance and storage reserves in a development-dependent manner. Plant Physiol. 2011, 157, 1026–1042. [Google Scholar] [CrossRef] [PubMed]
  51. Matsuyama, A.; Yoshimura, K.; Shimizu, C.; Murano, Y.; Takeuchi, H.; Ishimoto, M. Characterization of glutamate decarboxylase mediating gamma-amino butyric acid increase in the early germination stage of soybean (Glycine max [L.] Merr). J. Biosci. Bioeng. 2009, 107, 538–543. [Google Scholar] [CrossRef] [PubMed]
  52. Yang, R.; Hui, Q.; Gu, Z. Effects of ABA and CaCl2 on GABA accumulation in fava bean germinating under hypoxia-NaCl stress. Biosci. Biotechnol. Biochem. 2016, 80, 540–546. [Google Scholar] [CrossRef] [PubMed]
  53. Ling, Y.; Chen, T.; Jing, Y.; Fan, L.; Wan, Y.; Lin, J. γ-Aminobutyric acid (GABA) homeostasis regulates pollen germination and polarized growth in Picea wilsonii. Planta 2013, 238, 831–843. [Google Scholar] [CrossRef] [PubMed]
  54. Takahashi, H.; Greenway, H.; Matsumura, H.; Tsutsumi, N.; Nakazono, M. Rice alcohol dehydrogenase 1 promotes survival and has a major impact on carbohydrate metabolism in the embryo and endosperm when seeds are germinated in partially oxygenated water. Ann. Bot. 2014, 113, 851–859. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Miro, B.; Longkumer, T.; Entila, F.D.; Kohli, A.; Ismail, A.M. Rice seed germination underwater: Morpho-physiological responses and the bases of differential expression of alcoholic fermentation enzymes. Front. Plant Sci. 2017, 8, 1857. [Google Scholar] [CrossRef] [PubMed]
  56. Zhou, S.F.; Sun, L.; Valdés, A.E.; Engström, P.; Song, Z.T.; Lu, S.J.; Liu, J.X. Membrane-associated transcription factor peptidase, site-2 protease, antagonizes ABA signaling in Arabidopsis. New Phytol. 2015, 208, 188–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Guo, Y.; Song, Y.; Zheng, H.; Zhang, Y.; Guo, J.; Sui, N. NADP-malate dehydrogenase of sweet sorghum improves salt tolerance of Arabidopsis thaliana. J. Agric. Chem. 2018, 66, 5992–6002. [Google Scholar] [CrossRef] [PubMed]
  58. Sew, Y.S.; Stroher, E. Loss of Mitochondrial Malate Dehydrogenase Activity Alters Seed Metabolism Impairing Seed Maturation and Post-Germination Growth in Arabidopsis. Plant Physiol. 2016, 171, 849–863. [Google Scholar] [PubMed]
  59. Glevarec, G.; Bouton, S.; Jaspard, E.; Riou, M.T.; Cliquet, J.B.; Suzuki, A.; Limami, A.M. Respective roles of the glutamine synthetase/glutamate synthase cycle and glutamate dehydrogenase in ammonium and amino acid metabolism during germination and post-germinative growth in the model legume Medicago truncatula. Planta 2004, 219, 286–297. [Google Scholar] [CrossRef] [PubMed]
  60. Jain, M.; Tyagi, A.K.; Khurana, J.P. Constitutive expression of a meiotic recombination protein gene homolog, OsTOP6A1, from rice confers abiotic stress tolerance in transgenic Arabidopsis plants. Plant Cell Rep. 2008, 27, 767–778. [Google Scholar] [CrossRef] [PubMed]
  61. John, R.; Ganeshan, U.; Singh, B.N.; Kaul, T.; Reddy, M.K.; Sopory, S.K.; Rajam, M.V. Over-expression of topoisomerase II enhances salt stress tolerance in tobacco. Front. Plant Sci. 2016, 7, 1280. [Google Scholar] [CrossRef] [PubMed]
  62. Xing, M.; Sun, C.; Li, H.; Hu, S.; Lei, L.; Kang, J. Integrated analysis of transcriptome and proteome changes related to the Ogura cytoplasmic male sterility in cabbage. PLoS ONE 2018, 13, e0193462. [Google Scholar] [CrossRef] [PubMed]
  63. Zhu, H.G.; Cheng, W.H.; Tian, W.G.; Li, Y.J.; Liu, F.; Xue, F.; Zhu, Q.H.; Sun, Y.Q.; Sun, J. iTRAQ-based comparative proteomic analysis provides insights into somatic embryogenesis in Gossypium hirsutum L. Plant Mol. Biol. 2018, 96, 89–102. [Google Scholar] [CrossRef] [PubMed]
  64. Menotta, M.; Orazi, S.; Gioacchini, A.M.; Spapperi, C.; Ricci, A.; Chessa, L.; Magnani, M. Proteomics and transcriptomics analyses of ataxia telangiectasia cells treated with Dexamethasone. PLoS ONE 2018, 13, e0195388. [Google Scholar] [CrossRef] [PubMed]
  65. Ma, Q.; Liu, X.; Feng, W.; Liu, S.; Zhuang, Z. Analyses of the molecular mechanisms associated with salinity adaption of Trachidermus fasciatus through combined iTRAQ-based proteomics and RNA sequencing-based transcriptomics. Prog. Biophys. Mol. Biol. 2018, 136, 40–53. [Google Scholar] [CrossRef] [PubMed]
  66. Li, M.; Wang, K.; Li, S.; Yang, P. Exploration of rice pistil responses during early post-pollination through a combined proteomic and transcriptomic analysis. J. Proteom. 2016, 131, 214–226. [Google Scholar] [CrossRef] [PubMed]
  67. Peng, X.; Qin, Z.; Zhang, G.; Guo, Y.; Huang, J. Integration of the proteome and transcriptome reveals multiple levels of gene regulation in the rice dl2 mutant. Front. Plant Sci. 2015, 6, 351. [Google Scholar] [CrossRef] [PubMed]
  68. Sasaki, A.; Ashikari, M.; Ueguchi-Tanaka, M.; Itoh, H.; Nishimura, A.; Swapan, D.; Ishiyama, K.; Saito, T.; Kobayashi, M.; Khush, G.S.; et al. Green revolution: A mutant gibberellin-synthesis gene in rice. Nature 2002, 416, 701–702. [Google Scholar] [CrossRef] [PubMed]
  69. Kuroha, T.; Nagai, K.; Gamuyao, R.; Wang, D.R.; Furuta, T.; Nakamori, M.; Kitaoka, T.; Adachi, K.; Minami, A.; Mori, Y.; et al. Ethylene-gibberellin signaling underlies adaptation of rice to periodic flooding. Science 2018, 361, 181–186. [Google Scholar] [CrossRef] [PubMed]
  70. Hong, Z.; Ueguchi-Tanaka, M.; Fujioka, S.; Takatsuto, S.; Yoshida, S.; Hasegawa, Y.; Ashikari, M.; Kitano, H.; Matsuoka, M. The Rice brassinosteroid-deficient dwarf2 mutant, defective in the rice homolog of Arabidopsis DIMINUTO/DWARF1, is rescued by the endogenously accumulated alternative bioactive brassinosteroid, dolichosterone. Plant Cell 2005, 17, 2243–2254. [Google Scholar] [CrossRef] [PubMed]
  71. Li, Q.F.; Xiong, M.; Xu, P.; Huang, L.C.; Zhang, C.Q.; Liu, Q.Q. Dissection of brassinosteroid-regulated proteins in rice embryos during germination by quantitative proteomics. Sci. Rep. 2016, 6, 34583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Eisen, M.B.; Spellman, P.T.; Brown, P.O.; Botstein, D. Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. USA 1998, 95, 14863–14868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Mi, H.; Poudel, S.; Muruganujan, A.; Casagrande, J.T.; Thomas, P.D. PANTHER version 10: Expanded protein families and functions, and analysis tools. Nucleic Acids Res. 2016, 44, D336–D342. [Google Scholar] [CrossRef] [PubMed]
  74. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein-protein interaction networks, integrated over the tree of life. Nucleic. Acids. Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef] [PubMed]
  75. Wang, C.; Shen, L.; Fu, Y.; Yan, C.; Wang, K. A simple CRISPR/Cas9 system for multiplex genome editing in rice. J. Genet. Genom. 2015, 42, 703–706. [Google Scholar] [CrossRef] [PubMed]
  76. Li, Q.F.; Yu, J.W.; Lu, J.; Fei, H.Y.; Luo, M.; Cao, B.W.; Huang, L.C.; Zhang, C.Q.; Liu, Q.Q. Seed-specific expression of OsDWF4, a rate-limiting gene involved in brassinosteroids biosynthesis, improves both grain yield and quality in rice. J. Agric. Food Chem. 2018, 66, 3759–3772. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 03460 g001 550
Figure 1. Phenotype of near-isogenic line (NIL) NIP (sd1). (A) Gross morphology of NIP (SD1) and NIP (sd1) NILs grown in the paddy field. (B) Seed length of (top) NIP (SD1) and (bottom) NIP (sd1) NILs. (C) Seed width of (top) NIP (SD1) and (bottom) NIP (sd1) NILs. Scale bar = 1 cm. (D) Plant height, (E) tiller number, (F) flag leaf length, (G) flag leaf width, (H) grain length, and (I) grain width of NIP (SD1) and NIP (sd1) NILs.
Figure 1. Phenotype of near-isogenic line (NIL) NIP (sd1). (A) Gross morphology of NIP (SD1) and NIP (sd1) NILs grown in the paddy field. (B) Seed length of (top) NIP (SD1) and (bottom) NIP (sd1) NILs. (C) Seed width of (top) NIP (SD1) and (bottom) NIP (sd1) NILs. Scale bar = 1 cm. (D) Plant height, (E) tiller number, (F) flag leaf length, (G) flag leaf width, (H) grain length, and (I) grain width of NIP (SD1) and NIP (sd1) NILs.
Ijms 19 03460 g001
Ijms 19 03460 g002 550
Figure 2. Both brassinosteroid (BR) and gibberellin (GA) deficiency suppress rice shoot and root elongation during seed germination. (A) Shoot length of germinated seeds. (B) Root length of germinated seeds. Asterisks indicate significant differences between NIP/brassinazole (BRZ) or sd1/dimethyl sulphoxide (DMSO) and NIP/DMSO by Student’s t tests (* p < 0.05, ** p < 0.01).
Figure 2. Both brassinosteroid (BR) and gibberellin (GA) deficiency suppress rice shoot and root elongation during seed germination. (A) Shoot length of germinated seeds. (B) Root length of germinated seeds. Asterisks indicate significant differences between NIP/brassinazole (BRZ) or sd1/dimethyl sulphoxide (DMSO) and NIP/DMSO by Student’s t tests (* p < 0.05, ** p < 0.01).
Ijms 19 03460 g002
Ijms 19 03460 g003 550
Figure 3. Venn diagram illustrating the number of common responsive proteins identified from embryos of germinated seeds in both BRZ-treated NIP and the sd1 mutant.
Figure 3. Venn diagram illustrating the number of common responsive proteins identified from embryos of germinated seeds in both BRZ-treated NIP and the sd1 mutant.
Ijms 19 03460 g003
Ijms 19 03460 g004 550
Figure 4. Hierarchical cluster analysis and Clusters of Orthologous Groups (COG) classification of differentially abundant proteins in the embryos of BRZ-treated NIP and the sd1 mutant. (A) Hierarchical cluster analysis of common proteins in response to both BRZ treatment and OsGA20ox2 mutation. (B) COG classification of all common proteins.
Figure 4. Hierarchical cluster analysis and Clusters of Orthologous Groups (COG) classification of differentially abundant proteins in the embryos of BRZ-treated NIP and the sd1 mutant. (A) Hierarchical cluster analysis of common proteins in response to both BRZ treatment and OsGA20ox2 mutation. (B) COG classification of all common proteins.
Ijms 19 03460 g004
Ijms 19 03460 g005 550
Figure 5. Functional analyses of proteins responsive to both BR and GA deficiency among different gene ontology (GO) categories in embryos of germinated seeds. Common responsive proteins are divided into molecular function (A), biological process (B), cellular component (C) and protein class (D).
Figure 5. Functional analyses of proteins responsive to both BR and GA deficiency among different gene ontology (GO) categories in embryos of germinated seeds. Common responsive proteins are divided into molecular function (A), biological process (B), cellular component (C) and protein class (D).
Ijms 19 03460 g005
Ijms 19 03460 g006 550
Figure 6. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) construction of a potential protein-protein interaction network responsive to both BR and GA deficiency during seed germination. The network, made with a medium confidence cutoff (0.4) using the k-means clustering method, includes five clusters presented as different colours. Line thickness indicates the strength of data support.
Figure 6. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) construction of a potential protein-protein interaction network responsive to both BR and GA deficiency during seed germination. The network, made with a medium confidence cutoff (0.4) using the k-means clustering method, includes five clusters presented as different colours. Line thickness indicates the strength of data support.
Ijms 19 03460 g006
Ijms 19 03460 g007 550
Figure 7. Validation of proteomic data by quantitative real-time polymerase chain reaction (qRT-PCR). (A) Eight selected common responsive proteins displaying consistent changes in protein abundance in response to both BR and GA deficiency. (B) Changes in transcription for the eight selected genes. (C) Seven selected common responsive proteins displaying opposite changes in protein abundance in response to BR and GA deficiency. (D) Changes in transcription of the seven selected genes. Ubiquitin-Conjugating Enzyme (UBC) was used as an internal control. Asterisks indicate significant differences according to Student’s t-tests (** p < 0.01).
Figure 7. Validation of proteomic data by quantitative real-time polymerase chain reaction (qRT-PCR). (A) Eight selected common responsive proteins displaying consistent changes in protein abundance in response to both BR and GA deficiency. (B) Changes in transcription for the eight selected genes. (C) Seven selected common responsive proteins displaying opposite changes in protein abundance in response to BR and GA deficiency. (D) Changes in transcription of the seven selected genes. Ubiquitin-Conjugating Enzyme (UBC) was used as an internal control. Asterisks indicate significant differences according to Student’s t-tests (** p < 0.01).
Ijms 19 03460 g007
Ijms 19 03460 g008 550
Figure 8. Germination analysis of grp mutant. (A) Schematic diagram of the target in GRP gene. PAM, the protospacer adjacent motif. (B) Rice seeds were dehulled and imbibed in Milli-Q water supplemented with GA, BL or DMSO solution in darkness for 84 h. The scale bar indicates 0.5 cm. (C) Rice seed germination rates after 84 h of imbibition. (D) Shoot length of germinated seeds after 84 h of imbibition. (E) Root length of germinated seeds after 84 h of imbibition. Asterisks indicate significant differences according to Student’s t-tests (* p < 0.05, ** p < 0.01).
Figure 8. Germination analysis of grp mutant. (A) Schematic diagram of the target in GRP gene. PAM, the protospacer adjacent motif. (B) Rice seeds were dehulled and imbibed in Milli-Q water supplemented with GA, BL or DMSO solution in darkness for 84 h. The scale bar indicates 0.5 cm. (C) Rice seed germination rates after 84 h of imbibition. (D) Shoot length of germinated seeds after 84 h of imbibition. (E) Root length of germinated seeds after 84 h of imbibition. Asterisks indicate significant differences according to Student’s t-tests (* p < 0.05, ** p < 0.01).
Ijms 19 03460 g008
Table
Table 1. Proteins down-regulated in embryos of germinated rice seeds in response to both BR and GA deficiency by iTRAQ analysis.
Table 1. Proteins down-regulated in embryos of germinated rice seeds in response to both BR and GA deficiency by iTRAQ analysis.
Gene IDNIP (BRZ/Mock) Fold Change asd1/NIP Fold Change bProtein Score cUnique Peptide dSequence Coverage (%) eDescription
Os01g06528000.6930.56813338.9Protein of unknown function DUF231
Os02g02488000.2080.586144337.1Similar to glutelin type-B2 precursor
Os02g02490000.4060.7768041022.7Glutelin
Os03g02316000.7990.71224139.4Similar to branched-chain-amino-acid aminotransferase 3
Os03g02407000.8260.731133110.1Similar to Erwinia induced protein 2
Os03g03379000.8310.818371513.5Similar to predicted protein
Os03g08120000.5870.70330526.6DNA topoisomerase type IIA
Os04g01657000.8040.68531239.5Cysteine synthase
Os05g03291000.5550.661917342.7Prolamin
Os06g01122000.8120.781274315.4Purine and other phosphorylases
Os07g02143000.4740.7631815228.3Seed allergenic protein RAG2 precursor
Os07g05296000.650.71333625.1Similar to thiazole biosynthetic enzyme 1-1
Os08g05304000.4650.757128.8Moco-containing protein
Os09g04842000.8080.6152271758Hypothetical protein
Os11g02136000.6820.785180511.6Peptidase S10
a NIP (BRZ/Mock) fold change indicates fold change of protein abundance between BRZ treated samples and the mock-treated control. b sd1/NIP fold change indicates fold change of protein abundance between sd1 mutant and the NIP control. c Protein score is based on combined mass spectroscopy (MS) and tandem mass spectroscopy (MS/MS) spectra. The proteins that had a statistically significant (p < 0.05) protein score of 70 or more were considered successfully identified. d Unique peptide indicates the number of identified unique peptides for confident protein identification. e The protein with a sequence coverage larger than 5% was considered as credible.
Table
Table 2. Proteins up-regulated in embryos of germinated rice seeds in response to both BR and GA deficiency by iTRAQ analysis.
Table 2. Proteins up-regulated in embryos of germinated rice seeds in response to both BR and GA deficiency by iTRAQ analysis.
Gene IDNIP(BRZ/Mock)sd1/NIPProtein Unique Sequence Description
Fold Change aFold Change bScore cPeptide dCoverage (%) e
Os01g03584001.2591.2292745143.4Similar to 40S ribosomal protein S4
Os01g08808001.5851.506891516.8Similar to acyl-[acyl-carrier-protein] desaturase
Os02g05501001.381.238419110.8Similar to vacuolar ATP synthase 16 kDa proteolipid subunit
Os02g05870001.5292.233569220.8Similar to glycine-rich protein
Os03g01924001.2581.382257326.8GRIM-19 family protein
Os03g07500001.2881.5792653125.9Similar to ethylene-responsive protein
Os03g07742001.4841.261170418.4Similar to NADH-ubiquinone oxidoreductase subunit 8
Os03g07990001.9971.316126218.4Similar to histone H1
Os04g02496001.3711.272151223.9Rhodanese-like domain containing protein
Os05g05126001.2211.428365215.4X8 domain-containing protein
Os07g01194001.5541.554142411.6Similar to pectin esterase-like protein
Os08g04658001.3041.357952833.4Similar to glutamate decarboxylase
Os09g05395001.3671.402784336.6Similar to SKP1-like protein 1A
Os11g02105001.251.4271940639.8Similar to alcohol dehydrogenase
a NIP (BRZ/Mock) fold change indicates fold change of protein abundance between BRZ treated samples and the mock-treated control. b sd1/NIP fold change indicates fold change of protein abundance between sd1 mutant and the NIP control. c Protein score is based on combined mass spectroscopy (MS) and tandem mass spectroscopy (MS/MS) spectra. The proteins that had a statistically significant (p < 0.05) protein score of 70 or more were considered successfully identified. d Unique peptide indicates the number of identified unique peptides for confident protein identification. e The protein with a sequence coverage larger than 5% was considered as credible.
Table
Table 3. Proteins inconsistently regulated in embryos of germinated rice seeds in response to BR and GA deficiency by iTRAQ analysis.
Table 3. Proteins inconsistently regulated in embryos of germinated rice seeds in response to BR and GA deficiency by iTRAQ analysis.
Gene IDNIP (BRZ/Mock)sd1/NIPProtein Unique Sequence Description
Fold Change aFold Change bScore cPeptide dCoverage (%) e
Os01g02105000.7621.3661133645Similar to SOUL-like protein
Os01g02330001.2370.611391534.3DREPP plasma membrane polypeptide family protein
Os01g02947001.2390.694707833.7Haem peroxidase, plant/fungal/bacterial family protein
Os02g02093000.691.33378113.7Hypothetical conserved gene
Os03g03411001.2340.767506428.3Similar to 60S ribosomal protein L18
Os03g03791001.2520.81117636.9Protein of unknown function DUF248
Os03g07004001.2970.4081056337.3Similar to LOX4
Os03g08429000.7681.2810471343.9Similar to steroleosin-B
Os04g01653000.8251.2717428.9Conserved hypothetical protein
Os04g03908000.8181.24425201654.8NAD(P)-binding domain containing protein
Os04g04972001.2750.73516949.3Cellulase precursor
Os04g05465000.8071.31836213.5Similar to oleosin
Os05g02685000.4511.478131512Similar to serine carboxypeptidase 2
Os05g04744001.2530.775197210.9Prenylated rab acceptor PRA1 family protein
Os06g06757000.4631.229831418Similar to high pI alpha-glucosidase
Os11g05824000.6261.38570225.1Similar to embryo-specific protein
a NIP (BRZ/Mock) fold change indicates fold change of protein abundance between BRZ treated samples and the mock-treated control. b sd1/NIP fold change indicates fold change of protein abundance between sd1 mutant and the NIP control. c Protein score is based on combined mass spectroscopy (MS) and tandem mass spectroscopy (MS/MS) spectra. The proteins that had a statistically significant (p < 0.05) protein score of 70 or more were considered successfully identified. d Unique peptide indicates the number of identified unique peptides for confident protein identification. e The protein with a sequence coverage larger than 5% was considered as credible.

Share and Cite

MDPI and ACS Style

Li, Q.-F.; Wang, J.-D.; Xiong, M.; Wei, K.; Zhou, P.; Huang, L.-C.; Zhang, C.-Q.; Fan, X.-L.; Liu, Q.-Q. iTRAQ-Based Analysis of Proteins Co-Regulated by Brassinosteroids and Gibberellins in Rice Embryos during Seed Germination. Int. J. Mol. Sci. 2018, 19, 3460. https://doi.org/10.3390/ijms19113460

AMA Style

Li Q-F, Wang J-D, Xiong M, Wei K, Zhou P, Huang L-C, Zhang C-Q, Fan X-L, Liu Q-Q. iTRAQ-Based Analysis of Proteins Co-Regulated by Brassinosteroids and Gibberellins in Rice Embryos during Seed Germination. International Journal of Molecular Sciences. 2018; 19(11):3460. https://doi.org/10.3390/ijms19113460

Chicago/Turabian Style

Li, Qian-Feng, Jin-Dong Wang, Min Xiong, Ke Wei, Peng Zhou, Li-Chun Huang, Chang-Quan Zhang, Xiao-Lei Fan, and Qiao-Quan Liu. 2018. "iTRAQ-Based Analysis of Proteins Co-Regulated by Brassinosteroids and Gibberellins in Rice Embryos during Seed Germination" International Journal of Molecular Sciences 19, no. 11: 3460. https://doi.org/10.3390/ijms19113460

APA Style

Li, Q. -F., Wang, J. -D., Xiong, M., Wei, K., Zhou, P., Huang, L. -C., Zhang, C. -Q., Fan, X. -L., & Liu, Q. -Q. (2018). iTRAQ-Based Analysis of Proteins Co-Regulated by Brassinosteroids and Gibberellins in Rice Embryos during Seed Germination. International Journal of Molecular Sciences, 19(11), 3460. https://doi.org/10.3390/ijms19113460

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop