Skip to main content
  • Research article
  • Open access
  • Published:

Investigation of cell wall proteins of C. sinensis leaves by combining cell wall proteomics and N-glycoproteomics

Abstract

Background

C. sinensis is an important economic crop with fluoride over-accumulation in its leaves, which poses a serious threat to human health due to its leaf consumption as tea. Recently, our study has indicated that cell wall proteins (CWPs) probably play a vital role in fluoride accumulation/detoxification in C. sinensis. However, there has been a lack in CWP identification and characterization up to now. This study is aimed to characterize cell wall proteome of C. sinensis leaves and to develop more CWPs related to stress response. A strategy of combined cell wall proteomics and N-glycoproteomics was employed to investigate CWPs. CWPs were extracted by sequential salt buffers, while N-glycoproteins were enriched by hydrophilic interaction chromatography method using C. sinensis leaves as a material. Afterwards all the proteins were subjected to UPLC-MS/MS analysis.

Results

A total of 501 CWPs and 195 CWPs were identified respectively by cell wall proteomics and N-glycoproteomics profiling with 118 CWPs in common. Notably, N-glycoproteomics is a feasible method for CWP identification, and it can enhance CWP coverage. Among identified CWPs, proteins acting on cell wall polysaccharides constitute the largest functional class, most of which might be involved in cell wall structure remodeling. The second largest functional class mainly encompass various proteases related to CWP turnover and maturation. Oxidoreductases represent the third largest functional class, most of which (especially Class III peroxidases) participate in defense response. As expected, identified CWPs are mainly related to plant cell wall formation and defense response.

Conclusion

This was the first large-scale investigation of CWPs in C. sinensis through cell wall proteomics and N-glycoproteomics. Our results not only provide a database for further research on CWPs, but also an insight into cell wall formation and defense response in C. sinensis.

Background

Plant cell walls are a primary subcellular structure and are located in the outside of the cells. They offer the skeletal framework to tissues and play essential roles in protection, cell-to-cell adhesion and communication. Cell walls are mainly composed of complex polysaccharidic networks of celluloses, hemicelluloses, and pectins with a small proportion of cell wall proteins (CWPs), lignins, and lipids [1]. Among them, CWPs constitute around 10% of cell wall dry weight [2,3,4], but play important roles in various kinds of biological processes including cell wall metabolism, cell wall composition and modification, cell enlargement, signal transduction, biotic and abiotic stress response and other physiological processes [5,6,7,8,9,10].

In view of the importance of CWP function, the identification and characterization of CWPs have been widely studied in some plant species such as Arabidopsis [11,12,13,14,15,16,17,18,19], B. distachyon [20,21,22], flax [23, 24], sugarcane [10, 25, 26], rice [27,28,29], and others in recent decades by cell wall proteomics strategy using destructive and non-destructive extraction methods. These studies have greatly contributed to a broader knowledge of CWPs. However, to our knowledge, there is still a lack of understanding about CWPs due to difficult extraction and high contamination of intracellular proteins. N-Glycosylation is a common form of eukaryotic protein post-translational modification, and most plant proteins are N-glycosylated through the conventional endoplasmic reticulum (ER)-golgi apparatus (GA) secretory pathway [30, 31]. Consequently, N-glycosylation of plant CWPs is particularly prevalent and extensive. Conversely, large-scale and detailed characterization of N-glycoproteins has a great potential to increase our understanding of CWPs, and therefore N-glycoproteomics can be employed to investigate CWPs [32,33,34,35,36]. Based on these findings, this study combined N-glycoproteomics and cell wall proteomics to investigate CWPs.

C. sinensis is an important woody economic crop cultivated widely from tropical to temperate regions, its leaves are usually used for making tea. It is reported that the leaves of C. sinensis can accumulate much higher level of fluoride (F) than those of most other plants without exhibiting any toxicity symptoms under normal soil conditions [37,38,39,40], suggesting that there may be a special mechanism responsible for F accumulation/detoxification. Previous research has shown that cell wall immobilization and vacuolar compartmentation contribute to F accumulation/detoxification [40, 41], and recently we have found the important roles of CWPs in F accumulation/detoxification by a comparative proteomics analysis [42]. However, CWP identification and characterization have rarely been studied in C. sinensis up to now.

Herein, to broaden the knowledge of CWPs and provide a base for revealing the molecular mechanisms underlying F accumulation/detoxification-related CWPs, cell wall proteomics and N-glycoproteomics profiling of C. sinensis leaves was performed. In this study, CaCl2, EGTA, and LiCl were used sequentially to extract CWPs, and hydrophilic interaction chromatography (HILIC) was also employed to enrich N-glycoproteins. The peptides of obtained proteins were analyzed by ultrahigh performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS). Afterwards, all the identified proteins were subjected to multiple bioinformatics analysis. All in all, 578 CWPs were identified by combining cell wall proteomics and N-glycoproteomics. This study makes first attempt to investigate cell wall proteome and N-glycoproteome in C. sinensis. Our results will expand the understanding of CWPs and reveal the mechanism related to plant growth, development, and defense response.

Results

Identification of CWPs

To identify more CWPs, a combining strategy of cell wall proteomics and glycoproteomics was employed in this study as shown in Fig. 1. After UPLC-MS/MS analysis and database search, a total of 3618 target CWPs (TCWPs) and 262 N-glycoproteins were identified from C. sinensis leaves (Additional files 1, 2: Table S1, S2). To pick out CWPs, all identified proteins (3880) were subjected firstly to WallProtDB database retrieval. Among them, 627 TCWPs and 187 N-glycoproteins were identified as potential CWPs.

Fig. 1
figure 1

Experiment workflow. The extraction, precipitation, digestion, fractionation, and MS/MS and data analyses of TCWPs were operated according to blue arrow instruction. Likewise, those of glycoproteins were operated according to red arrow instruction. The figure is under copyright, each researcher can use and adapt it by citing our paper

According to previous reports, only those proteins (i) having a predicted signal peptide (SP), (ii) lacking ER retention signal (KDEL or HDEL) and (iii) no more than one transmembrane domain (TMD) were defined as CWPs [11, 12, 25]. To obtain CWPs as many as possible and enhance CWP coverage, all identified proteins were subjected to multiple bioinformatics analyses including SP, TMD, ER retention signal, and subcellular localization. Based on the above-mentioned definition of CWPs and the report by Day et al. [24], a total of 501 TCWPs and 195 N-glycoproteins were identified as CWPs. Among them, 484 TCWPs and 187 N-glycoproteins were also retrieved in WallProtDB database, whereas 17 TCWPs and 8 N-glycoproteins were absent and thus determined firstly as CWPs (Table 1; Additional file 3: Table S3). As for the remaining proteins, 38 TCWPs were designated as plasma membrane proteins (Additional file 4: Table S4), others including 3079 TCWPs and 67 N-glycoproteins was defined as intracellular proteins. Taken together, 501 CWPs were identified by cell wall proteomics and 195 CWPs were identified by N-glycoproteomics, respectively, 118 CWPs were in common through both approaches (Table 1; Fig. 2A; Additional file 3: Table S3).

Table 1 578 CWPs identified from C. sinensis leaves
Fig. 2
figure 2

Identification (A) and functional classification (B) of CWPs identified from C. sinensis leaves

Functional classification of CWPs

To better understand the biological functions of CWPs, CWPs were categorized on the basis of their functional domains proposed by Jamet et al. [6]. A total of 578 CWPs (501 + 195–118) were divided into nine groups (Fig. 2B). Among them, proteins acting on polysaccharides (PACs; 147) were the largest functional class, occupying 25.4% of total CWPs. Proteases (Ps; 94) were the second largest class, accounting for 16.3% of total CWPs. Oxido-reductases (ORs; 62) were the third largest class, occupying 10.7% of identified CWPs, followed by proteins involved in signaling (PSs; 56), proteins related to lipid metabolism (PLMs; 43), and proteins with interaction domains (PIDs;34), accounting for 9.7, 7.4 and 5.9% of the identified CWPs, respectively. Structural proteins (SPs; 5, 0.9%) had lowest abundance, only containing 5 members. The remaining CWPs related to various functions were categorized as miscellaneous proteins (MPs; 75, 13.0%), and CWPs with previously uncharacterized domains were referred to as proteins of unknown function (PUFs; 62, 10.7%).

CWP comparison between C. sinensis and other two species

As expected, the functional distribution of CWPs identified from C. sinensis leaves was in good concordance with that from A. thaliana rosettes and B. distachyon leaves (Additional file 5: Fig.S1) with PACs, Ps, and ORs representing top three functional classes. Notably, the percentage of PSs in C. sinensis (9.7%) leaves was obviously higher than that in A. thaliana rosettes (3.7%) and B. distachyon (4.0%) leaves, respectively [12, 20]. Such a difference may be attributed to the longer lifecycle of the woody evergreen leaf in C. sinensis.

Main representative of functional classes

In PAC class, GHs are the major representative (Table 1). In this study, a total of 110 GHs were identified accounting for 74.8% in PAC class, which fell into 23 families including GH1, GH3, GH5, GH9, GH10, GH13, GH16, GH17, GH18, GH19, GH20, GH27, GH28, GH29, GH31, GH32, GH35, GH37, GH38, GH51, GH65, GH79, and GH127 according to CAZy nomenclature based on sequence homology (Fig. 3). As expected, the most representative families were GH3 and GH17, as previously documented [12, 20, 26]. Moreover, GH1, GH5, GH16, GH18, GH19, GH27, GH28, GH31, GH35, and GH38 were also well representative families with at least 5 members in each family (Fig. 3). In addition, less representative CWPs acting on polysaccharides were also identified, including carbohydrate esterase [11 pectinesterases (known as pectin methylesterases (PMEs)) and 3 pectinesterase inhibitors (PMEIs)], 4 glycosyl transferases (GTs, including GT2, GT31, GT48, and GT68), 6 expansins, 4 PNGase A, 3 pectin acetylesterases (PAEs), 2 pectate lyases (PLs), and 4 carbohydrate acylation (trichome birefringence-like proteins).

Fig. 3
figure 3

Glycoside hydrolases identified from C. sinensis leaves

In Ps class, Asp proteases (28), Ser carboxypeptidases (28), Ser proteases (19), and Cys proteases (14) are main families, occupying 94.7% of Ps. ORs functional class mainly comprised of class III peroxidase (PODs, 29), multicopper oxidases (13), BBE (berberine bridge enzyme) (S)-reticulin (6), and laccases (5). Other CWPs related to redox processes were identified including monocopper oxidase-like proteins (SKU5 and SKS1), blue copper proteins, and ascorbate oxidases.

PSs class mainly contained fasciclin-like arabinogalactan proteins (FLAs, 9) and receptor-like protein kinases (RLKs) superfamily proteins (38). Among them, RLKs comprised 21 LRR-RLKs, 6 cysteine-rich recetor-like protein kinases, 3 S-locus receptor kinase subfamily proteins, 2 wall-associated receptor kinases, and 6 lectin receptor kinase subfamily proteins. PLMs class mainly consisted of lipid-transfer proteins (LTPs, 10) and GDSL esterase/lipases (GDSLs, 16). As for SPs class, only five CWPs were identified in this study including 3 leucine-rich repeat extensin-like proteins (LRR-EXTs), 1 non-classical arabinogalactan protein 31-like (AGP), and 1 hydroxyproline-rich glycoprotein (HPRG). Identified MPs mainly encompasses purple acid phosphatases (PAPs, 17), blue copper binding proteins (BCPs, 9), dirigent proteins (DIRs, 8), germin-like proteins (GLPs, 5), thaumatins (7), and proteins with a cupin domain (5).

Discussion

Identification and functional classification of identified CWPs

Totally, 3618 TCWPs were identified in C. sinensis leaves by sequential salt extractions and UPLC-MS/MS. Among them, 627 TCWPs were homologs of the CWPs in WallProtDB database, whereas 501 TCWPs were in good accordance with the CWPs defined by multiple bioinformatics analyses. In except for firstly defined CWPs, there is an identification difference between WallProtDB database and bioinformatics analyses, which might ascribe to a low homology of CWPs between C. sinensis and other plant species indexed in WallProtDB.

Finally, 501 TCWPs and 3079 TCWPs were designated as CWPs and intracellular proteins, indicating that TCWPs were subjected to the contamination during TCWPs preparation. Similar high contamination of intracellular proteins was also detected in sugarcane [25] and rice [29], accounting for 81.6 and 80.5%, respectively. So far there have been rare cell wall proteomics studies, therefore CWP extraction remains to be improved. In spite of the high contamination of intracellular proteins, this study adopted cell wall proteomics to enhance the CWP coverage in C. sinensis.

At the same time, most identified N-glycoproteins (195, 74.4%) were targeted into the cell wall/extracellular/plasma membrane in C. sinensis, and thus they were designated as CWPs. Our results agreed well with those of the studies of tomato fruit [35] and B. distachyon leaf [43], in which 65 and 60% of N-glycoproteins were found to be located in the apoplast/cell wall/plasma membrane, respectively, demonstrating that N-glycoproteomics is a feasible method to identify and characterize CWPs.

It should be noted that 25 CWPs were newly identified ones in this study (Additional file 3: Table S3), and that more CWPs were identified through cell wall proteomics (501 CWPs) than N-glycoproteomics (195CWPs), indicating cell wall proteomics was more effective than N-glycoproteomics for CWP identification. However, the use of N-glycoproteomics as supplementation can further enhance CWP identification effect. Considering this, we propose that combined cell wall proteomics and N-glycoproteomics during CWP identification and characterization.

Possible functions of identified CWPs

Identified CWPs acting on cell wall polysaccharides

Glycoside hydrolases (GHs)

GHs are the overwhelming majority of identified CWPs with 19.0%. Possible substrates of most GH families are hemicelluloses (xyloglucan, xylans, glucomannans) and pectin (galactans, homogalacturonan). Of GHs identified in this study, GH16, GH29, GH31, and GH65 potentially act on xyloglucans; GH10 and GH51 show a possible action on xylans; and GH28 and GH35 can hydrolyze homogalacturonan and galactans, respectively [23, 44, 45] (Additional file 6: Table S5). Moreover, GH1, GH3, and GH5 possess broad substrates range, and their enzymes are reported to be involved in the modification and/or breakdown of cell wall hemicelluloses and pectins [46, 47], and to participate in lignification and secondary metabolism [48]. Identification of these GH families suggested that hemicelluloses and pectins might undergo important structural changes in the leaves of C. sinensis. Furthermore, GH127 (also known as DUF1680 domain protein), recently characterized as a novel β-L-arabinofuranosidase, might be involved in the degradation of cell wall polysaccharides and hydroxyproline-rich glycoproteins [49], and GH9 was known to catalyze the endohydrolysis of cellulose.

Some identified GHs might participate in defense against pathogens and various stresses. Chitin and β-1,3- or β-1,6-glucan are main components of cell walls of various fungi. GH17 acts as β-1,3-glucanase; GH18 and GH19 act as chitinases; and GH20 functions as key hydrolyzation enzyme of chitin, and these four GHs possess antifungal activity to degrade fungus cell walls and participate in defense against pathogens [45, 50]. Chitinase has been reported to respond to abiotic stress [42, 51]. GH37, a non-reducing sugar, was identified as a new CWP in this study and it has been found to be a universal stabiliser of protein conformation and probably contribute to various stress defense [52].

Several identified GHs including GH13, GH27, and GH32 might be implicated in the mobilization, allocation, and partitioning of storage reserves. GH13 is involved in the hydrolysis of starch and glycogen to yield glucose and maltose [53]. GH27 is one of three hydrolyzing enzymes of galactomannans (cell wall storage polysaccharide) [54], and GH32, as an invertase, is involved in long-distance nutrient allocation and carbohydrate partitioning [55, 56]. Additionally, a couple of GH enzymes including GH3, GH18, GH19, GH35, GH38, and GH79 are involved in post-translational modifications (PTMs) of glycoproteins [32, 45]. In this study, GH3, GH35, GH38, and GH79 were verified as N-glycoproteins.

Collectively, a large number of GHs associated with cell wall metabolism and defense were identified in this work, which is consistent with previous reports of sugarcane stems and leaves [26], B. distachyon grains [21], S. officinarum cell suspension [25]. Our data reveal the potential functions of identified GHs such as complex cell wall carbohydrate remodeling, pathogen and stress response, mobilization and allocation of storage reserves, and glycoprotein PTMs. Our results might be attributed to sustainable remodeling during plant growth and development and terrestrial habit of plants.

Other CWPs acting on polysaccharides

PMEs, PAEs, and PLs are pectin-modifying enzymes. PMEs catalyse the demethyl-esterification of homogalacturonan domain of pectin [57]. The degree of pectin methylation/demethylation affects cell wall stiffening and access to enzymes [58]. Demethyl-esterificated pectin favors the cleavage of the acidic polygalacturonic chains by GH28 and PLs. Likewise, PAEs can regulate pectin deacetyltation by cleaving the acetylester bond from pectin [59]. Overall, these enzymes play a major role in controlling cell wall plasticity/rheology by affecting pectin metabolism [60].

Trichome birefringence-like proteins and PNGase A are also two modification enzyme families in cell wall. The former is characterized as xylan acetyltransferases, and it is associated with the xylan O-acetylation mediation, secondary wall deposition, and pathogen resistance [61]. The latter is one of deglycosylation enzyme, and it is involved in the release of N-glycans from glycopeptides generated by the proteolysis of denatured glycoproteins [62].

Expansin, known as non-enzymatic protein and the most important structural protein, plays a central role in cell wall extension via their action on the cellulose-hemicellulose network, suggesting that expansin is essential for primary cell wall structure during plant growth- and development-related processes [63]. In addition, 4 cell wall GT families might be associated with the biosynthesis of cell wall polymers.

Identified CWPs functioning as proteases

Proteases are necessary for protein turnover, maturation of enzymes, and defense against pathogens [45, 64]. Consequently, plant proteases localized in the cell wall might be responsible for CWP degradation or maturation, and further they might play crucial roles in a striking variety of biological processes such as plant growth and branching, flower time regulation, and in defense responses.

Identified CWPs involved in redox

Class III PODs, a large multigene families, accounted for one half of OR functional class. Class III PODs are involved in lignin metabolism by catalyzing the oxidative polymerization of monolignols [65], stress responses, and signaling transduction via consuming hydrogen peroxide and generating reactive oxygen species [66]. Class III PODs can mediate cross-linking of cell wall compounds such as structural proteins, monolignols, and aromatic amino acids containing polysaccharides [67,68,69]. Like class III PODs, laccases are candidates for polymerizing monolignol unit into lignin, suggesting that laccases are essential for cell wall lignification [70, 71]. BBE-like proteins, as monolignol oxidoreductases, may participate in the mobilization and oxidation of monolignols required for polymerization processes [72]. Overall, three highly representative enzyme families in the redox class were considered to be involved in ligin production and subsequent the reinforcement of cell wall strength and rigidity, which supported plant defense against adverse environmental factors.

In addition, monocopper oxidase-like proteins (SKU5 and SKS1), blue copper proteins, and ascorbate oxidases were found, they might play a role in both cell wall loosening, expansion, and reticulation processes [24, 73].

Identified CWPs involved in signaling transduction

In this study, identified signaling transduction-related CWPs mainly consisted of FLAs and RLKs superfamily proteins. FLAs, heavily O-glycosylated CWPs, have been found to be correlated with cell wall formation [74], cell-to-cell adhesion and communication [75], and abiotic stress response [76]. RLKs, as primary cell wall “sensors”, are responsible for controlling diverse signaling events [77]. RLKs possess important functions in a wide variety of development- and defense-related processes, for example, they can recognize extracellular ligand to activate the intracellular kinase domain, resulting in downstream signaling transduction [78].

Identified CWPs related to lipid metabolism

CWPs in the functional class have been reported to be related to lipid metabolism [79,80,81,82,83]. LTPs are required for lipid export to the cell surface, and they are closely associated with cutin and wax formation [79]. One LTP has been found to be involved in cell wall extension by interacting with the cellulose/xyloglucan network [80]. GDSLs, a recently discovered subclass of lipolytic enzymes, possess multifunctional properties, and play important roles not only in the formation of surface cutin and epi-cuticular wax [81], but also in the tolerance to biotic and abiotic stresses [82, 83]. In summary, numerous LTPs and GDSLs might play important roles in cuticle assembly during the growth and development of C. sinensis leaf. The identification of CWPs related to PLMs is conducive to our understanding of leathery leaf of C. sinensis.

Structure proteins

Due to be still resistant to salt-extraction, structure proteins were eluted difficultly so far. In this study, five structure proteins were identified. LRR-EXTs have been reported to influence mechanical properties of cell wall by forming insolubilized, covalently crosslink with cell wall components [84], and they can perceive extracellular signals and indirectly relay them into the cytoplasm to regulate plant growth and salt tolerance, and consequent they are important for cell wall development, plant growth, and stress tolerance [85]. Non-classical AGPs have both proline-rich domain and non-proline-rich domain, may function in metal ion-binding, defense response, and they can interact with pectin [86, 87]. HPRG is an important structural components of plant cell walls, and are related to structural integrity, cell-cell interaction, and intercellular communication [88].

Identified CWPs related to other functions

Regarding several MPs-related CWPs, PAPs might be associated with the degradation of xyloglucan and oligosaccharides via dephosphorylating CWPs such as alpha xylosidase and β glucosidase [89]. DIRs are related to lignin polymerization [90, 91], and they play important roles in various stress responses and cell wall modification/reinforcement during cell wall integrity maintenance [92]. BCPs, GLPs, cupins, and Thaumatins have been previously reported to be associated with stress responses in plants [93,94,95,96].

Several enzymes of CWP inhibitor were also detected in this study, including PMEIs, PGIPs (polygalacturonase inhibitor-like), and Cys proteinase inhibitor. PMEIs partially inhibit the activity of PMEs and adjust the degree of pectin methyl-esterification. PGIPs specifically bind with polygalacturonases (GH28), thereby inhibiting the hydrolyzation of pectin and regulating pectin degradation, eventually triggering defense response against microbes and insects [97]. In summary, two couples (PMEIs and PME, PGIPs and PG) occurred coincidentally, and they modulate precisely pectin metabolism. Cys proteinase inhibitor exhibit inhibitory activities against specific Cys proteases, thus might function in preventing insect predation [98].

Roles of various CWPs in plant cell wall formation and defense response

Under dynamic environmental conditions, plants grow and develop continuously, and they always encounter various stresses and deleterious attack from insects and microbes. Plant cell walls, as the first barrier, change constantly to be adapted to environmental stresses. Doubtlessly, CWPs play central roles in altering cell wall properties. On basis of our results, a work model of identified CWPs mainly related to plant cell wall formation and defense response was proposed (Fig. 4).

Fig. 4
figure 4

A work model of CWPs identified from C. sinensis leaves. The figure is under copyright, each researcher can use and adapt it by citing our paper

To satisfy the requirement of normal growth and development, a large number of CWPs are activated to adjust cell wall structure. In this study, numerous CWPs related to PACs were identified, mainly including GH1, GH3, GH5, GH9, GH10, GH16, GH28, GH29, GH31, GH35, GH51, and GH65, and they might contribute to the rearrangement of cell wall structure, while expansin probably give rise to cell wall extension. Several CWPs associated with the formation and metabolism of secondary cell wall, such as class III PODs, BBEs, laccases, LTPs, GDSLs, and DIRs, may favor the reinforcement/modification of cell wall (Fig. 4).

Confronted with adverse environment, C. sinensis, a terrestrial plant, has no ability to escape. Therefore, it has evolved some mechanisms of defense responses such as altering cell wall properties. Many CWPs identified in this study are potentially involved in various defense. GH17, GH18, GH19, and GH20 have been reported to be involved mainly in response to pathogen stress as well as abiotic stress by hydrolyzing chitin. Class III PODs, monocopper oxidase-like proteins, blue copper proteins, and ascorbate oxidases are involved in response to various biotic and abiotic stresses by redox reaction. LTPs, GDSLs, and DIRs are also associated with defense response through the regulation of secondary cell wall. PGIPs and Cys proteinase inhibitor might function in improving protection against insects and pathogens [99] via inhibiting the activity of degradation enzymes of invaders. Likewise, BCPs, GLPs, cupins, and thaumatins also function in defense response (Fig. 4).

To sense dynamic environment and changing complex cell wall structures, plants have developed cell wall integrity-sensing pathway to transduce signals into cytoplasm. A number of sensors on the plasma membrane including RLKs and FLAs were identified in present study, which can mediate cross-talk between the cell wall and the cytoplasm in C. sinensis (Fig. 4).

Conclusions

This study combined cell wall proteomics and N-glycoproteomics to identify CWPs in C. sinensis. A total of 3880 proteins were identified by sequential salt extraction and UPLC-MS/MS. Meanwhile, 262 N-glycoproteins were identified by HILIC enrichment coupled with UPLC-MS/MS. Subsequently, 501 out of 3880 proteins and 195 out of 262 N-glycoproteins were designated as CWPs by multiple bioinformatics analysis. Of these designated CWPs, 118 were in common. In total, 578 CWPs were identified from C. sinensis leaves, 25 of which were determined as newly identified CWPs. This study was the first attempt of large-scale investigation of CWPs by cell wall proteomics and N-glycoproteomics in C. sinensis. It provides a reference for using a combined strategy of cell wall proteomics and N-glycoproteomics to improve CWP identification and characterization. Our findings promote the understanding of cell wall formation and defense response in C. sinensis.

Methods

Plant materials

From the tree top, the first to fifth leaves were collected from 20 uniform 2-year-old cutting seedlings of the Echa 1 variety (C. sinensis cv. ‘Echa 1’) in tea germplasm bank located in Wuhan city of Hubei province (China). The collected leaves were washed three times with Milli-Q water, and ground into fine power in liquid nitrogen immediately, and finally stored at − 80 °C for subsequent use.

Cell wall enrichment

Cell wall fractions were obtained from the leaves of C. sinensis using sequential washes as described by Printz et al. [100] with slight modification. Briefly, 5 g fine power of the leaves were homogenized with 3-fold volumes of 0.4 M sucrose buffer for 10 min, vortexed for 2 min, shaken overnight at 250 rpm at 4 °C, and then centrifuged. Subsequently, 0.6 M sucrose buffer was added into the precipitations, shaken for 30 min at 250 rpm at 4 °C and centrifuged. Afterwards, 1 M sucrose buffer was added into the precipitations again, suspended, and centrifuged. Finally, the precipitations were washed twice using 5 mM sodium acetate buffer. The final precipitations were cell wall fractions (pellet). Sucrose buffer contained 5 mM sodium acetate and 1% protease inhibitor cocktail (ApexBio). All the buffers (pH 4.6) were precooled at 4 °C, and the centrifugation was performed at 1000 rpm for 15 min at 4 °C.

Cell wall protein extraction

CWPs were extracted successively using CaCl2, EGTA, and LiCl according to the method reported by Printz et al. [100]. Briefly, 0.2 M CaCl2 buffer was firstly added into cell wall pellet, shaken for 30 min at 200 rpm at 4 °C, followed by centrifugation. Then the supernatants were collected. This process was repeated once, and the supernatants were pooled as CaCl2 fractions. Afterwards, cell wall pellets were mixed with 50 mM EGTA buffer, followed by shaking for 1 h at 300 rpm at 37 °C, centrifugation and supernatant collection. This procedure was repeated twice, and all the supernatants were collected as EGTA fractions. Cell wall pellets were finally re-suspended in 3 M LiCl buffer, homogenized overnight at 250 rpm, 4 °C. After centrifugation, the supernatants were collected once again. The CWPs were once again extracted from the cell wall pellets with 3 M LiCl buffer by shaking for 6 h at 250 rpm, 4 °C. The obtained supernatants were pooled and stored as LiCl fractions. Finally, CaCl2, EGTA, and LiCl fractions were combined as target CWPs (TCWPs) fractions. All extraction buffers were precooled at 4 °C, and the centrifugation was performed for 15 min at 10000 rpm at 4 °C.

Whole protein extraction

Whole proteins (WPs) were extracted from C. sinensis leaves according to several previous reports [42, 101, 102]. Briefly, about 0.5 g fine powder was firstly homogenized with 5 ml pre-cooled homogenization buffer [containing 20 mM Tris-HCl (pH 7.5), 250 mM sucrose, 10 mM EGTA, 1 mM PMSF, 1 mM DTT, and 1% (v/v) Triton], and then centrifuged at 12000 g for 20 min at 4 °C. The obtained supernatants were pooled and stored as whole protein fractions.

Protein precipitation and cleaning

According to our previous study [42], WP fractions and TCWP fractions were precipitated severally by Tris-phenol (pH ≥ 8.0) and ammonium acetate. In brief, the fractions were mixed with equal volume of Tris-phenol, and vortexed, followed by centrifugation at 12000 g for 20 min at 4 °C. Afterwards, the phenol phase was transferred carefully into other tubes, mixed thoroughly with 5 volumes of 0.1 M ammonium acetate in 100% methanol, and incubated at − 80 °C overnight. The precipitated proteins were washed twice with 0.1 M ammonium acetate and acetone, separately. The protein pellets were lyophilized and then dissolved into lysis buffers [containing 7 M urea, 2 M thiourea, 4% CHAPS, 250 mM DTT, and 0.2% (v/v) Bio-Lyte]. Protein concentration was determined with BCA kit according to the manufacturer’s instructions.

Protein digestion

Before trypsin digestion, WPs and TCWPs were reduced with 5 mM dithiothreitol for 30 min at 56 °C, and alkylated with 11 mM iodoacetamide for 15 min at room temperature in darkness, followed by urea dilution to concentration < 2 M through the addition of 100 mM triethylammonium bicarbonate. Afterwards, WPs and TCWPs were digested firstly by trypsin (1:50 trypsin/protein) overnight at 37 °C, and then by trypsin (1:100 trypsin/protein) for 4 h. Finally, tryptic peptides were desalted by Strata X C18 SPE column (Phenomenex, USA) and concentrated by centrifugal concentrator.

HPLC fractionation

After tryptic digestion, the peptides from WPs and TCWPs were fractionated severally by the use of high pH reversed-phase HPLC (high-performance liquid chromatography) with Agilent 300 Extend C18 column (5 μm particles, 4.6 mm inner diameter, and 250 mm length). Briefly, the digested peptides were first separated into 60 fractions with a gradient of 8 to 32% acetonitrile (pH 9.0) for more than 60 min. Subsequently, the peptides were pooled into 4 fractions and dried by vacuum centrifugation for further use.

Affinity enrichment of N-glycopeptides

To enrich N-glycosylation peptides, the dried peptides from WPs were firstly dissolved in 40 μL enrichment buffer (containing 80% acetonitrile and 1% trifluoroacetic acid), and then loaded into HILIC micro-column to separate glycopeptides from non-glycopeptides by centrifugation for 15 min at 4000 g. To remove non-specifically adsorbed peptides, HILIC micro-column was washed three times with enrichment buffers. Subsequently, the bound peptides were eluted from the micro-column with 10% acetonitrile, and then vacuum-dried. The lyophilized N-glycopeptides were reconstituted in 50 μL NH4CO3 buffer (50 mM) in heavy oxygen water and incubated with 2 μL PNGase F at 37 °C overnight. Finally, the resultant N-glycopeptides were desalted with C18 ZipTips (Millipore) according to the manufacturer’s instructions and lyophilized for LC-MS/MS analysis.

UPLC-MS/MS analysis

For LC-MS/MS analysis, the peptides were firstly dissolved in solvent A (containing 0.1% (v/v) formic acid and 2% acetonitrile), and then gradient-eluted in EASY-nLC 1000 UPLC system. Peptide separation was conducted with home-made reversed phase column (25 cm length, 100 μm ID). TCWP peptides were gradient eluted as follows: 450 nL/min constant flow; starting with 7% ~ 25% solvent B (containing 0.1% formic acid in 90% acetonitrile) for 0–40 min, followed by 25% ~ 35% solvent B for 40–52 min, 35% ~ 80% solvent B for 52–56 min, and 80% solvent B for 56–60 min. Deglycosylated peptides were gradient eluted with following procedures: 500 nL/min constant flow, starting with 4 to 20% solvent B for 0–24 min, 20 to 32% solvent B for 24–32 min, 32 to 80% solvent B for 32–36 min, and finally maintaining in 80% solvent B for 36–40 min.

Subsequently, the separated TCWP peptides and deglycosylated peptides were respectively injected into a nanoelectrospray ion source, followed by MS/MS analysis in Q Exactive™ and Orbitrap Fusion mass spectrometer (Thermo Fisher scientific). Briefly, the applied electrospray voltage was 2.0 kV, and the intact peptides and their secondary fragments were detected and analyzed by Orbitrap with a data-dependent acquisition mode automatically switching between MS scan and MS/MS scan.

TCWP peptides were fully scanned at a resolution of 70,000 with m/s scan range of 350–1800. Afterwards, the top 10 most intense parent ions per scan were selected for higher-energy collisional dissociation fragmentation (HCD) at 28% collision energy. The generated fragments were further analyzed at a resolution of 17,500 with a fixed first mass of 100 m/z. To increase the effective utilization rate of mass spectrometry, the related parameters were set as follows: automatic gain control of 5E4, 30 s dynamic exclusion, 100 ms maximum inject, and signal threshold of 20,000 ions/s. Likewise, deglycosylated peptides were fully scanned at a resolution of 60,000 with m/s scan range of 350–1550. The top 20 most intense parent ions per scan were selected for HCD at 35% collision energy, and then the resultant fragments were analyzed at a resolution of 15,000 with a fixed first mass of 100 m/z. Similarly, the related MS parameters were set as follows: automatic gain control of 5E4, 15 s dynamic exclusion, 200 ms maximum inject and signal threshold of 5000 ions/s were used.

Database search

The resultant raw MS/MS data were processed using MaxQuant search engine (v.1.5.2.8) with the following query parameters: (i) tea tree genome database (Camellia_sinensis_4442 with 53,512 sequences [103];) concatenated with reverse decoy database and mass spectrometry contaminants database for MS/MS search; (ii) Trypsin/P for enzyme cleavage and 2 missing cleavages; (iii) mass tolerance of 20 ppm and 5 ppm for peptide ions in first search and main research, respectively, and 0.02 Da for fragment ions; (iv) length of 7 amino acid residues as minimum peptide length, and 5 as maximum modification number in a peptide; (v) Cysteine alkylation as fixed modification; (vi) Variable modification: methionine oxidation and N-terminal acetylation for TCWPs, and methionine oxidation and deamidation (NQ), asparagine deamidation (18O)for N-glycoproteins; (vii) FDR ≤ 1% for protein identification and peptide-spectrum match identification.

Multiple bioinformatics analyses

CWPs were predicted and functionally categorized using WallProtDB database [104]. Glycoside hydrolases and carbohydrate esterase were grouped according to CAZy database [105]. N-terminal signal peptide of identified proteins was predicted using SignalP [106]. Transmembrane domain was evaluated by TMHMM server [107]. Subcellular localization predication was performed using TargetP [108], WoLF PSORT [109], Loctree 3 [110], and Plant-mPLoc [111]. ER retention signal was checked using Prosite [112].

Availability of data and materials

All generated or analyzed data were included in this published article, and all relevant raw data were deposited in ProteomeXchange Consortium with the dataset identifier PXD026772.

Abbreviations

CWPs:

Cell wall proteins

TCWPs:

Target cell wall proteins

WPs:

Whole proteins

GHs:

Glycoside hydrolases

PACs:

Proteins acting on polysaccharides

Ps:

Proteases

ORs:

Oxido-reductases.

PSs:

Proteins involved in signaling.

PLMs:

Proteins related to lipid metabolism.

PIDs:

Proteins with interaction domains.

PODs:

Peroxidases.

References

  1. Roberts K. How the cell wall acquired a cellular context. Plant Physiol. 2001;125(1):127–30. https://doi.org/10.1104/pp.125.1.127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Carpita NC, Gibeaut DM. Structural models of primary cell walls in flowering plants, consistency of molecular structure with the physical properties of the walls during growth. Plant J. 1993;3(1):1–30. https://doi.org/10.1111/j.1365-313X.1993.tb00007.x.

    Article  CAS  PubMed  Google Scholar 

  3. Cassab GI, Varner JE. Cell wall proteins. Annu Rev Plant Physiol Plant Mol Biol. 1992;39(4):321–53.

    Google Scholar 

  4. Albenne C, Canut H, Elisabeth J. Plant cell wall proteomics: the leadership of Arabidopsis thaliana. Front Plant Sci. 2013;4(4):111.

    PubMed  PubMed Central  Google Scholar 

  5. Zhu JM, Chen SX, Alvarez S, Asirvatham VS, Schachtman DP, Wu YJ, et al. Cell wall proteome in the maize primary root elongation zone. I. Extraction and identification of water-soluble and lightly ionically bound proteins. Plant Physiol. 2006;140(1):311–25. https://doi.org/10.1104/pp.105.070219.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Jamet E, Canut H, Boudart G, Pont-Lezica RF. Cell wall proteins: a new insight through proteomics. Trends Plant Sci. 2006;11(1):33–9. https://doi.org/10.1016/j.tplants.2005.11.006.

    Article  CAS  PubMed  Google Scholar 

  7. Sasidharan R, Voesenek LACJ, Pierik R. The regulation of cell wall extensibility during shade avoidance: a study using two contrasting ecotypes of Stellaria longipes. Plant Physiol. 2008;148(3):1557–69. https://doi.org/10.1104/pp.108.125518.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Krzesłowska M. The cell wall in plant cell response to trace metals: polysaccharide remodeling and its role in defense strategy. Acta Physiol Plant. 2011;33(1):35–51. https://doi.org/10.1007/s11738-010-0581-z.

    Article  CAS  Google Scholar 

  9. Wang X, Komatsu S. Plant subcellular proteomics: application for exploring optimal cell function in soybean. J Proteome. 2016;143(4):45–56. https://doi.org/10.1016/j.jprot.2016.01.011.

    Article  CAS  Google Scholar 

  10. Calderan-Rodrigues MJ, Jame E, Douché T, Bonassi MBR, Catald TR, Fonseca JG, et al. Cell wall proteome of sugarcane stems: comparison of a destructive and a nondestructive extraction method showed differences in glycoside hydrolases and peroxidases. BMC Plant Biol. 2016;16(1):14. https://doi.org/10.1186/s12870-015-0677-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Duruflé H, San Clemente H, Balliau T, Zivy M, Dunand C, Jamet E. Cell wall proteome analysis of Arabidopsis thaliana mature stems. Proteomics. 2017;17(8):1600449. https://doi.org/10.1002/pmic.201600449.

    Article  CAS  Google Scholar 

  12. Hervé V, Duruflé H, San Clemente H, Albenne C, Balliau T, Zivy M, et al. An enlarged cell wall proteome of Arabidopsis thaliana rosettes. Proteomics. 2016;16(24):3183–7. https://doi.org/10.1002/pmic.201600290.

    Article  CAS  PubMed  Google Scholar 

  13. Nguyen-Kim H, San Clemente H, Balliau T, Zivy M, Dunand C, Albenne C, et al. Arabidopsis thaliana root cell wall proteomics: increasing the proteome coverage using a combinatorial peptide ligand library and description of unexpected Hyp in peroxidase amino acid sequences. Proteomics. 2016;16(3):491–503. https://doi.org/10.1002/pmic.201500129.

    Article  CAS  PubMed  Google Scholar 

  14. Chivasa S, Ndimba BK, Simon WJ, Robertson D, Yu XL, Knox JP, et al. Proteomic analysis of the Arabidopsis thaliana cell wall. Electrophoresis. 2002;23(11):1754–65. https://doi.org/10.1002/1522-2683(200206)23:11<1754::AID-ELPS1754>3.0.CO;2-E.

    Article  CAS  PubMed  Google Scholar 

  15. Borderies G, Jamet E, Lafitte C, Rossignol M, Jauneau A, Boudart G, et al. Proteomics of loosely bound cell wall proteins of Arabidopsis thaliana cell suspension cultures: a critical analysis. Electrophoresis. 2003;24(19–20):3421–32. https://doi.org/10.1002/elps.200305608.

    Article  CAS  PubMed  Google Scholar 

  16. Feiz L, Irshad M, Pont-Lezica RF, Canut H, Jamet E. Evaluation of cell wall preparations for proteomics: a new procedure for purifying cell walls from Arabidopsis hypocotyls. Plant Methods. 2006;2(1):10. https://doi.org/10.1186/1746-4811-2-10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Irshad M, Canut H, Borderies G, Pont-Lezica R, Jamet E. A new picture of cell wall protein dynamics in elongating cells of Arabidopsis thaliana: confirmed actors and newcomers. BMC Plant Biol. 2008;8(1):e94.

    Article  Google Scholar 

  18. Bayer EM, Bottrill AR, Walshaw J, Vigouroux M, Naldrett MJ, Thomas CL, et al. Arabidopsis cell wall proteome defined using multidimensional protein identification technology. Proteomics. 2006;6(1):301–11. https://doi.org/10.1002/pmic.200500046.

    Article  CAS  PubMed  Google Scholar 

  19. Boudart G, Jamet E, Rossignol M, Lafitte C, Borderies G, Jauneau A, et al. Cell wall proteins in apoplastic fluids of Arabidopsis thaliana rosettes: identification by mass spectrometry and bioinformatics. Proteomics. 2005;5(1):212–21. https://doi.org/10.1002/pmic.200400882.

    Article  CAS  PubMed  Google Scholar 

  20. Douché T, San Clemente H, Burlat V, Roujol D, Valot B, Zivy M, et al. Brachypodium distachyon as a model plant toward improved biofuel crops: search for secreted proteins involved in biogenesis and disassembly of cell wall polymers. Proteomics. 2013;13(16):2438–54. https://doi.org/10.1002/pmic.201200507.

    Article  CAS  PubMed  Google Scholar 

  21. Francin-Allami M, Merah K, Albenne C, Rogniaux H, Pavlovic M, Lollier V, et al. Cell wall proteomic of Brachypodium distachyon grains: a focus on cell wall remodeling proteins. Proteomics. 2015;15(13):2296–306. https://doi.org/10.1002/pmic.201400485.

    Article  CAS  PubMed  Google Scholar 

  22. Francin-Allami M, Lollier V, Pavlovic M, San Clemente H, Rogniaux H, Jamet E, et al. Understanding the remodelling of cell walls during Brachypodium distachyon grain development through a sub-cellular quantitative proteomic approach. Proteomes. 2016;4(3):21. https://doi.org/10.3390/proteomes4030021.

    Article  CAS  PubMed Central  Google Scholar 

  23. Chabi M, Goulas E, Lille U. A cell wall proteome and targeted cell wall analyses provide novel information on hemicellulose metabolism in flax. Mol Cell Proteomics. 2017;16(9):1634–51. https://doi.org/10.1074/mcp.M116.063727.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Day A, Fénart S, Neutelings G, Hawkins S, Rolando C, Tokarski C. Identification of cell wall proteins in the flax (Linum usitatissimum) stem. Proteomics. 2013;13(5):812–25. https://doi.org/10.1002/pmic.201200257.

    Article  CAS  PubMed  Google Scholar 

  25. Calderan-Rodrigues MJ, Jamet E, Bonassi MBCR, Guidetti-Gonzalez S, Begossi AC, Setem L, et al. Cell wall proteomics of sugarcane cell suspension cultures were predicted to be localized in non-secretory pathway. Proteomics. 2014;14(6):738–49. https://doi.org/10.1002/pmic.201300132.

    Article  CAS  PubMed  Google Scholar 

  26. Fonseca JG, Calderan-Rodrigues MJ, de Moraes FE, Cataldi TR, Jamet E, Labate CA. Cell wall proteome of sugarcane young and mature leaves and stems. Proteomics. 2018;189(2):9853.

    Google Scholar 

  27. Chen XY, Kim ST, Cho WK, Cho WK, Rim Y, Kim S, et al. Proteomics of weakly bound cell wall proteins in rice calli. J Plant Physiol. 2009;166(7):675–85. https://doi.org/10.1016/j.jplph.2008.09.010.

    Article  CAS  PubMed  Google Scholar 

  28. Pandey A, Rajamani U, Verma J, Pratigya S, Navjyoti C, Asis D, et al. Identification of extracellular matrix proteins of rice (Oryza sativa L.) involved in dehydration-responsive network: a proteomic approach. J Proteome Res. 2010;9(7):3443–64. https://doi.org/10.1021/pr901098p.

    Article  CAS  PubMed  Google Scholar 

  29. Cho WK, Hyun TK, Kumar D, Rim Y, Chen XY, Jo Y, et al. Proteomic analysis to identify tightly-bound cell wall protein in rice calli. Mol Cells. 2015;38(8):685–96. https://doi.org/10.14348/molcells.2015.0033.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Strasser R. Biological significance of complex N-glycans in plants and their impact on plant physiology. Front Plant Sci. 2004;5:363.

    Google Scholar 

  31. Nguema-Ona E, Vicré-Gibouin M, Gotté M, Plancot B, Lerouge P, Bardor M, et al. Cell wall O-glycoproteins and N-glycoproteins: aspects of biosynthesis and function. Front Plant Sci. 2014;5:499.

    Article  Google Scholar 

  32. Minic Z, Jamet E, Négroni L. Arsene der Garabedian P, Zivy M, et al. a sub-proteome of Arabidopsis thaliana mature stems trapped on Concanavalin a is enriched in cell wall glycoside hydrolases. J Exp Bot. 2007;58(10):2503–12. https://doi.org/10.1093/jxb/erm082.

    Article  CAS  PubMed  Google Scholar 

  33. Catala C, Howe KJ, Hucko S, Rose JK, Thannhauser TW. Towards characterization of the glycoproteome of tomato (Solanum lycopersicum) fruit using concanavalin a lectin affinity chromatography and LC-MALDI-MS/MS analysis. Proteomics. 2011;11(8):1530–44. https://doi.org/10.1002/pmic.201000424.

    Article  CAS  PubMed  Google Scholar 

  34. Ruiz-May E, Kim SJ, Brandizzi F, Rose JKC. The secreted plant N-glycoproteome and associated secretory pathways. Front Plant Sci. 2012;3:117.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Ruiz-May E, Hucko S, Howe KJ, Zhang S, Sherwood RW, Thannhauser TW, et al. A comparative study of lectin affinity based plant N-glycoproteome profiling using tomato fruit as a model. Mol Cell Proteomics. 2013;13(2):566–79. https://doi.org/10.1074/mcp.M113.028969.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Barba-Espin G, Dedvisitsakul P, Hagglund P, Svensson B, Finnie C. Gibberellic acid-induced aleurone layers responding to heat shock or tunicamycin provide insight into the N-glycoproteome, protein secretion, and endoplasmic reticulum stress. Plant Physiol. 2014;164(2):951–65. https://doi.org/10.1104/pp.113.233163.

    Article  CAS  PubMed  Google Scholar 

  37. Dey S, Giri B. Fluoride fact on human health and health problems: a review. Med Clin Rev. 2016;2(1):11.

    Article  Google Scholar 

  38. Ruan J, Wong MH. Accumulation of fluoride and aluminium related to different varieties of tea plant. Environ Geochem Health. 2001;23(1):53–63. https://doi.org/10.1023/A:1011082608631.

    Article  CAS  Google Scholar 

  39. Xie ZL, Chen Z, Sun WT, Guo XJ, Yin B, Wang JH. Distribution of aluminum and fluoride in tea plant and soil of tea garden in central and Southwest China. Chin Geogr Sci. 2007;17(4):376–82. https://doi.org/10.1007/s11769-007-0376-3.

    Article  Google Scholar 

  40. Gao HJ, Zhao Q, Zhang XC, Wan XC, Mao JD. Localization of fluoride and aluminum in subcellular fractions of tea leaves and roots. J Agr Food Chem 2014;62(10):2313–2319, DOI: https://doi.org/10.1021/jf4038437.

    Article  CAS  Google Scholar 

  41. Cai HM, Peng CY, Chen J, Hou RY, Gao HJ, Wan XC. X-ray photoelectron spectroscopy surface analysis of fluoride stress in tea (Camellia sinensis (L.) O. Kuntze) leaves. J Fluor Chem. 2014;158:11–5. https://doi.org/10.1016/j.jfluchem.2013.11.012.

    Article  CAS  Google Scholar 

  42. Liu YL, Cao D, Ma LL, Jin XF, Yang PF, Ye F, et al. TMT-based quantitative proteomics analysis reveals the response of Camellia sinensis to fluoride. J Proteome. 2018;176:71–81. https://doi.org/10.1016/j.jprot.2018.02.001.

    Article  CAS  Google Scholar 

  43. Zhang M, Chen GX, Lv DW, Li XH, Yan YM. N-linked glycoproteome profiling of seedling leaf in Brachypodium distachyon L. J Proteome. 2015;14(4):1727–38. https://doi.org/10.1021/pr501080r.

    Article  CAS  Google Scholar 

  44. Augur C, Stiefel V, Darvill A, Albersheim P, Puigdomenech P. Molecular cloning and pattern of expression of an α-L-fucosidase gene from pea seedlings. J Biol Chem. 1995;270(42):24839–43. https://doi.org/10.1074/jbc.270.42.24839.

    Article  CAS  PubMed  Google Scholar 

  45. Boudart G, Minic Z, Albenne C, Canut H, Jamet E, Pont-Lezica RF. Cell Wall Proteome. In: Šamaj J, Thelen J, editors. Plant proteomics. Berlin: Springer; 2007. p. 169–85.

    Google Scholar 

  46. Minic Z. Physiological roles of plant glycoside hydrolases. Planta. 2008;227(4):723–40. https://doi.org/10.1007/s00425-007-0668-y.

    Article  CAS  PubMed  Google Scholar 

  47. Minic Z, Jouanin L. Plant glycoside hydrolases involved in cell wall polysaccharide degradation. Plant Physiol Biochem. 2006;44(7–9):435–49. https://doi.org/10.1016/j.plaphy.2006.08.001.

    Article  CAS  PubMed  Google Scholar 

  48. Ketudat Cairns JR, Esen A. β -glucosidases. Cell Mol Life Sci. 2010;67(20):3389–405. https://doi.org/10.1007/s00018-010-0399-2.

    Article  CAS  PubMed  Google Scholar 

  49. Fujita K, Takashi Y, Obuchi E, Kitahara K, Suganuma T. Characterization of a novel beta-L-arabinofuranosidase in Bifidobacterium longum: functional elucidation of a DUF1680 protein family member. J Biol Chem. 2014;289(8):5240–9. https://doi.org/10.1074/jbc.M113.528711.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Rast DM, Baumgartner D, Mayer C, Hollenstein GO. Cell wall-associated enzymes in fungi. Phytochemistry. 2003;64(2):339–66. https://doi.org/10.1016/S0031-9422(03)00350-9.

    Article  CAS  PubMed  Google Scholar 

  51. Hradilova J, Rehulka P, Rehulkova H, Vrbova M, Griga M, Brzobohaty B. Comparative analysis of proteomic changes in contrasting flax cultivars upon cadmium exposure. Electrophoresis. 2010;31(2):421–31. https://doi.org/10.1002/elps.200900477.

    Article  CAS  PubMed  Google Scholar 

  52. Frison M, Parrou JL, Guillaumot D, Masquelier D, François J, Chaumont F, et al. The Arabidopsis thaliana trehalase is a plasma membrane-bound enzyme with extracellular activity. FEBS Lett. 2007;581(21):4010–6. https://doi.org/10.1016/j.febslet.2007.07.036.

    Article  CAS  PubMed  Google Scholar 

  53. Calderan-Rodrigues MJ, Fonseca JG, Edgar de Moraes F, Setem LV, Begossi AC, Labate CA. Plant cell wall proteomics: a focus on monocot species, Brachypodium distachyon, Saccharum spp. and Oryza sativa. Int J Mol Sci. 2019;20(8):1975. https://doi.org/10.3390/ijms20081975.

    Article  CAS  PubMed Central  Google Scholar 

  54. DeMason DA, Madore MA, Chandra Sekhar KN, Harris MJ. Role of α-galactosidase in cell wall metabolism of date palm (Phoenix dactylifera) endosperm. Protoplasma. 1992;166(3-4):177–86. https://doi.org/10.1007/BF01322780.

  55. Roitsch T, González MC. Function and regulation of plant invertases: sweet sensations. Trends Plant Sci. 2004;9(12):606–13. https://doi.org/10.1016/j.tplants.2004.10.009.

    Article  CAS  PubMed  Google Scholar 

  56. Veillet F, Gaillard C, Coutos-Thévenot P, La Camera S. Targeting the AtCWIN1 gene to explore the role of invertases in sucrose transport in roots and during Botrytis cinerea infection. Front Plant Sci. 2016;7:1899.

    Article  Google Scholar 

  57. Wolf S, Mouille G, Pelloux J. Homogalacturonan methyl-esterification and plant development. Mol Plant. 2009;2(5):851–60. https://doi.org/10.1093/mp/ssp066.

    Article  CAS  PubMed  Google Scholar 

  58. Atmodjo MA, Hao Z, Mohnen D. Evolving views of pectin biosynthesis. Ann Rev Plant Biol. 2013;64(1):747–79. https://doi.org/10.1146/annurev-arplant-042811-105534.

    Article  CAS  Google Scholar 

  59. Philippe F, Pelloux J, Rayon C. Plant pectin acetylesterase structure and function: new insights from bioinformatics analysis. BMC Genomics. 2017;18(1):456. https://doi.org/10.1186/s12864-017-3833-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Gall HL, Philippe F, Domon JM, Gillet F, Pelloux J, Rayon C. Cell wall metabolism in response to abiotic stress. Plants. 2015;4(1):112–66. https://doi.org/10.3390/plants4010112.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Gao YP, He CW, Zhang DM, Liu XL, Xu ZP, Tian YB, et al. Two trichome birefringence-like proteins mediate xylan acetylation, which is essential for leaf blight resistance in rice. Plant Physiol. 2017;173(1):470–81. https://doi.org/10.1104/pp.16.01618.

    Article  CAS  PubMed  Google Scholar 

  62. Megumi M, Yoshinobu K. Structural features of free N-glycans occurring in plants and functional features of de-N-glycosylation enzymes, ENGase, and PNGase: the presence of unusual plant complex type N-glycans. Front Plant Sci. 2014;5:429.

    Google Scholar 

  63. Cosgrove DJ. Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005;6(11):850–61. https://doi.org/10.1038/nrm1746.

    Article  CAS  PubMed  Google Scholar 

  64. Van der Hoorn RA. Plant proteases: from phenotypes to molecular mechanisms. Ann Rev Plant Biol. 2008;59(1):191–223. https://doi.org/10.1146/annurev.arplant.59.032607.092835.

    Article  CAS  Google Scholar 

  65. Fagerstedt KV, Kukkola EM, Koistinen VVT, Takahashi J, Marjamaa K. Cell wall lignin is polymerised by class III secretable plant peroxidases in Norway spruce. J Integr Plant Biol. 2010;52(2):186–94. https://doi.org/10.1111/j.1744-7909.2010.00928.x.

    Article  CAS  PubMed  Google Scholar 

  66. Shigeto J, Tsutsumi Y. Diverse functions and reactions of class III peroxidases. New Phytol. 2006;209(4):1395–402.

    Article  Google Scholar 

  67. Passardi F, Penel C, Dunand C. Performing the paradoxical: how plant peroxidases modify the cell wall. Trends Plant Sci. 2004;9(11):534–40. https://doi.org/10.1016/j.tplants.2004.09.002.

    Article  CAS  PubMed  Google Scholar 

  68. Wakabayashi K, Soga K, Hoson T. Phenylalanine ammonia-lyase and cell wall peroxidase are cooperatively involved in the extensive formation of ferulate networks in cell walls of developing rice shoots. J Plant Physiol. 2012;169(3):262–7. https://doi.org/10.1016/j.jplph.2011.10.002.

    Article  CAS  PubMed  Google Scholar 

  69. Francoz E, Ranocha P, Nguyen-Kim H, Jamet E, Burlat V, Dunand C. Roles of cell wall peroxidases in plant development. Phytochemistry. 2015;112(1):15–21. https://doi.org/10.1016/j.phytochem.2014.07.020.

    Article  CAS  PubMed  Google Scholar 

  70. Claus H. Laccases: structure, reaction, distribution. Micron. 2004;35(1–2):93–6. https://doi.org/10.1016/j.micron.2003.10.029.

    Article  CAS  PubMed  Google Scholar 

  71. Wang Y, Bouchabke-Coussa O, Lebris P, Antelme S, Soulhat C, Gineau E, et al. LACCASE5 is required for lignification of the Brachypodium distachyon culm. Plant Physiol. 2015;168(1):192–204. https://doi.org/10.1104/pp.114.255489.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Daniel B, Pavkov-Keller T, Steiner B, Dordic A, Gutmann A, Nidetzky B, et al. Oxidation of monolignols by members of the berberine bridge enzyme family suggests a role in plant cell wall metabolism. J Biol Chem. 2015;290(30:18770–81.

    Article  Google Scholar 

  73. Sedbrook JC, Carroll KL, Hung KF, Masson PH, Somerville CR. The Arabidopsis SKU5 gene encodes an extracellular glycosyl phosphatidylinositol-anchored glycoprotein involved in directional root growth. Plant Cell. 2002;14(7):1635–48. https://doi.org/10.1105/tpc.002360.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. MacMillan C, Taylor L, Bi Y, Southerton S, Southerton SG, Evans R, et al. The fasciclin-like arabinogalactan protein family of Eucalyptus grandis contains members that impact wood biology and biomechanics. New Phytol. 2015;206(4):1314–27. https://doi.org/10.1111/nph.13320.

    Article  CAS  PubMed  Google Scholar 

  75. Groover A, Ribischon M. Developmental mechanisms regulating secondary growth in woody plants. Curr Opin Plant Biol. 2006;9(1):55–8. https://doi.org/10.1016/j.pbi.2005.11.013.

    Article  CAS  PubMed  Google Scholar 

  76. Faik A, Abouzouhair J, Sarhan F. Putative fasciclin-like arabinogalactan-proteins (FLA) in wheat (TriticuM aestivum) and rice (Oryza sativa): identification and bioinformatic analyses. Mole Genet Genomics. 2006;276(5):478–94. https://doi.org/10.1007/s00438-006-0159-z.

    Article  CAS  Google Scholar 

  77. Osakabe Y, Yamaguchi-Shinozaki K, Shinozaki K, Tran LP. Sensing the environment: key roles of membrane-localized kinases in plant perception and response to abiotic stress. J Exp Bot. 2013;64(2):445–58. https://doi.org/10.1093/jxb/ers354.

    Article  CAS  PubMed  Google Scholar 

  78. Torii KU. Leucine-rich repeat receptor kinases in plants: structure, function, and signal transduction pathways. Int Rev Cytol 2004;234(0):1–46, DOI: https://doi.org/10.1016/S0074-7696(04)34001-5.

  79. Debono A, Yeats T, Rose J, Bird D, Jetter R, Kunst L, et al. Arabidopsis LTPG is a glycosylphosphatidylinositol-anchored lipid transfer protein required for export of lipids to the plant surface. Plant Cell. 2009;21(4):1230–8. https://doi.org/10.1105/tpc.108.064451.

  80. Nieuwland J, Feron R, Huisman BA, Fasolino A, Hilbers CW, Derksen J, et al. Lipid transfer proteins enhance cell wall extension in tobacco. Plant Cell. 2005;17(7):2009–19. https://doi.org/10.1105/tpc.105.032094.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Li C, Chen G, Mishina K, Yamaji N, Ma JF, Yukuhiro F, et al. A GDSL-motif esterase/ acyltransferase/lipase is responsible for leaf water retention in barley. Plant Direct. 2017;1(5):e00025. https://doi.org/10.1002/pld3.25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Gao M, Yin X, Yang W, Lam S, Tong X, Liu J, et al. GDSL lipases modulate immunity through lipid homeostasis in rice. PLoS Pathog. 2017;13(11):e1006724. https://doi.org/10.1371/journal.ppat.1006724.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Zhang HH, Wang ML, Li YQ, Yan W, Chang ZY, Ni HL, et al. GDSL esterase/lipases OsGELP34 and OsGELP110/OsGELP115 are essential for rice pollen development. J Integr Plant Biol. 2020;62(10):1574–93. https://doi.org/10.1111/jipb.12919.

    Article  CAS  PubMed  Google Scholar 

  84. Draeger C, Fabrice TN, Gineau E, Mouille G, Kuhn BM, Moller I, et al. Arabidopsis leucine-rich repeat extension (LRX) proteins modify cell wall composition and influence plant growth. BMC Plant Biol. 2015;15(1):155. https://doi.org/10.1186/s12870-015-0548-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Zhao CZ, Zayed O, Yu ZP, Jiang W, Zhu PP, Hsu CC, et al. Leucine-rich repeat extensin proteins regulate plant salt tolerance in Arabidopsis. P Natl Aca Sci USA. 2018;115(51):13123–8. https://doi.org/10.1073/pnas.1816991115.

    Article  CAS  Google Scholar 

  86. Liu CG, Mehdy MC. A nonclassical arabinogalactan protein gene highly expressed in vascular tissues, AGP31, is transcriptionally repressed by methyl jasmonic acid in Arabidopsis. Plant Physiol. 2007;145(3):863–74. https://doi.org/10.1104/pp.107.102657.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Gong SY, Huang GQ, Sun X, Li P, Zhao LL, Zhang DJ, et al. GhAGP31, a cotton non-classical arabinogalactan protein, is involved in response to cold stress during early seedling development. Plant Biol. 2012;14(3):447–57. https://doi.org/10.1111/j.1438-8677.2011.00518.x.

    Article  CAS  PubMed  Google Scholar 

  88. Jiao HJ, Liu X, Sun SG, Wang P, Qiao X, Li JM, et al. The unique evolutionary pattern of the Hydroxyproline-rich glycoproteins superfamily in Chinese white pear (Pyrus bretschneideri). BMC Plant Biol. 2018;18(1):36. https://doi.org/10.1186/s12870-018-1252-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Kaida R, Serada S, Norioka N, Norioka S, Neumetzler L, Pauly M, et al. Potential role for purple acid phosphatase in the dephosphorylation of cell wall proteins in tobacco cells. Plant Physiol. 2010;153(2):603–10. https://doi.org/10.1104/pp.110.154138.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Gang DR, Costa MA, Fujita M, Dinkova-Kostova AT, Wang HB, Burlat V, et al. Regiochemical control of monolignol radical coupling: a new paradigm for lignin and lignan biosynthesis. Chem Biol. 1996;6(3):133–41.

    Google Scholar 

  91. Pickel B, Pfannstiel J, Steudle A, Lehmann A, Gerken U, Pleiss J, et al. A model of dirigent proteins derived from structural and functional similarities with allene oxide cyclase and lipocalins. FEBS J. 2012;279(11):1980–93. https://doi.org/10.1111/j.1742-4658.2012.08580.x.

    Article  CAS  PubMed  Google Scholar 

  92. Paniagua C, Bilkova A, Jackson P, Dabravolski S, Ribe W, Didi V, et al. Dirigent proteins in plants: modulating cell wall metabolism during abiotic and biotic stress exposure. J Exp Bot. 2017;68(13):3287–301. https://doi.org/10.1093/jxb/erx141.

    Article  CAS  PubMed  Google Scholar 

  93. Ruan XM, Luo F, Li DD, Zhang J, Liu ZH, Xu WL, et al. Cotton BCP genes encoding putative blue copper-binding proteins are functionally expressed in fiber development and involved in response to high-salinity and heavy metal stresses. Physiol Plantarum. 2001;141(1):71–83.

    Article  Google Scholar 

  94. de Jesús-Pires C, Ferreira-Neto JRC, Pacifico Bezerra-Neto J, Kido EA, de Oliveira Silva RL, Pandolfi V, et al. Plant thaumatin-like proteins: function, evolution and biotechnological applications. Curr Protein Pep Sc. 2020;21(16):36–51. https://doi.org/10.2174/1389203720666190318164905.

    Article  CAS  Google Scholar 

  95. Li L, Xu X, Chen C, Shen Z. Genome-wide characterization and overexpression analysis of the Germin-like protein family in rice and Arabidopsis. Int J Mol Sci. 2016;17(10):1622. https://doi.org/10.3390/ijms17101622.

    Article  CAS  PubMed Central  Google Scholar 

  96. Gabrišova D, Klubicova K, Danchenko M, Gomory D, Berezhna VV, Skultety L, et al. Do cupins have a function beyond being seed storage proteins? Front Plant Sci. 2016;6:1–9.

    Article  Google Scholar 

  97. Kirsch R, Vurmaz E, Schaefer C, Eberl F, Sporer T, Haeger W, et al. Plants use identical inhibitors to protect their cell wall pectin against microbes and insects. Ecol Evol. 2020;10(8):3814–24. https://doi.org/10.1002/ece3.6180.

    Article  PubMed  PubMed Central  Google Scholar 

  98. Zhao Y, Botella MA, Subramanian L, Niu X, Nielsen SS, Bressan RA, et al. Two wound-inducible soybean cysteine proteinase inhibitors have greater insect digestive proteinase inhibitory activities than a constitutive homolog. Plant Physiol. 1996;111(4):1299–306. https://doi.org/10.1104/pp.111.4.1299.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  99. Ryan C. Protease inhibitors in plants: genes for improving defenses against insects and pathogens. Annu Rev Phytopathol. 1990;28(1):425–49. https://doi.org/10.1146/annurev.py.28.090190.002233.

    Article  CAS  Google Scholar 

  100. Printz B, Morais RDS, Wienkoop S, Sergeant K, Lutts S, Hausman JF, et al. An improved protocol to study the plant cell wall proteome. Front Plant Sci. 2015;6(3):237.

    PubMed  PubMed Central  Google Scholar 

  101. Liu YL, Chaturvedi P, Fu JL, Cai QQ, Weckwerth W, Yang PF. Induction and quantitative proteomic analysis of cell dedifferentiation during callus formation of Lotus (Nelumbo nucifera Gaertn. Spp. baijianlian). J Proteome. 2016;131:61–70. https://doi.org/10.1016/j.jprot.2015.10.010.

    Article  CAS  Google Scholar 

  102. Yang P, Li X, Wang X, Chen H, Chen F, Shen S. Proteomic analysis of rice (Oryza sativa) seeds during germination. Proteomics. 2007;7(18):3358–68. https://doi.org/10.1002/pmic.200700207.

    Article  CAS  PubMed  Google Scholar 

  103. Wei CL, Yang H, Wang SB, Zhao J, Liu C, Gao LP, et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. P Natl Aca Sci USA. 2018;115(18):E4151–8. https://doi.org/10.1073/pnas.1719622115.

    Article  CAS  Google Scholar 

  104. San Clemente H, Jamet E. WallProtDB, a database resource for plant cell wall proteomics. Plant Methods. 2015;11(1):2. https://doi.org/10.1186/s13007-015-0045-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  105. Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B. The carbohydrate-active enzymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 2009;37(Database):D233–8. https://doi.org/10.1093/nar/gkn663.

    Article  CAS  PubMed  Google Scholar 

  106. Almagro AJJ, Tsirigos KD, Sønderby CK, et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019;37(4):420–3. https://doi.org/10.1038/s41587-019-0036-z.

    Article  CAS  Google Scholar 

  107. Moller S, Croning MDR, Apweiler R. Evaluation of methods for the prediction of membrane spanning regions. Bioinformatics. 2001;17(7):646–53. https://doi.org/10.1093/bioinformatics/17.7.646.

    Article  CAS  PubMed  Google Scholar 

  108. Emanuelsson O, Nielsen H, Brunak S, von Heijne G. Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol. 2000;300(4):1005–16. https://doi.org/10.1006/jmbi.2000.3903.

    Article  CAS  PubMed  Google Scholar 

  109. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, et al. WoLF PSRT: protein localization predictor. Nucleic Acids Res. 2007;35(Web Server issue):W585–7.

    Article  Google Scholar 

  110. Goldberg T, Hecht M, Hamp T, Karl T, Yachdav G, Ahmed N, et al. LocTree3 prediction of localization. Nucleic Acids Res. 2014;42(Web Server issue):W350–5.

    Article  CAS  Google Scholar 

  111. Chou KC, Shen HB. Plant-mPLoc: a top-down strategy to augment the power for predicting plant protein subcellular localization. PLoS One. 2010;5(6):e11335. https://doi.org/10.1371/journal.pone.0011335.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  112. Sigrist CJA, Cerutti L, Hulo N, Gattiker A, Falquet L, Pagni M, et al. PROSITE: a documented database using patterns and profiles as motif descriptors. Brief Bioinform. 2002;3(3):265–74. https://doi.org/10.1093/bib/3.3.265.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Thanks Jingjie PTM BioLab (HangZhou) Co., Ltd. for MS/MS analysis. Thanks Mrs. Ping Liu (Huazhong Agricultural University) for overall help in checking the manuscript and refining the language.

Funding

This work was financially supported by Natural science Foundation of Hubei province (2019CFB600), and Key Research and Development Program of Hubei province (2020BBA038). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

YLL and XFJ conceived and designed the experiment; YLL conducted the experiment, drafted the manuscript, and modified continually the manuscript. LLM and DC assisted in performing the part of the experiment and data acquisition; JF and HJH assisted in data analyses. ZMG helped in checking and revising the manuscript. All author had read and approved the final version of the manuscript.

Corresponding author

Correspondence to Xiaofang Jin.

Ethics declarations

Ethics approval and consent to participate

The collection of plant materials in this study complied with the guidelines by Ministry of Agriculture and Rural Affairs of the People’s Republic of China (certificate no: GPD tea plant (2019)420015).

Consent for publication

Not Applicable.

Competing interests

The authors have declared no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Proteins identified from C. sinensis leaves by cell wall proteomics.

Additional file 2: Table S2.

N-glycoproteins identified from C. sinensis leaves.

Additional file 3: Table S3.

578 CWPs identified from C. sinensis leaves by cell wall proteomics and N-glycoproteomics.

Additional file 4: Table S4.

Plasma membrane proteins identified by cell wall proteomics.

Additional file 5: Fig.S1.

Functional categories of CWPs identified from C. sinensis leaves, B. distachyon leaves, and A. thaliana rosettes.

Additional file 6: Table S5.

Possible substrates of GHs identified from C. sinensis leaves.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Ma, L., Cao, D. et al. Investigation of cell wall proteins of C. sinensis leaves by combining cell wall proteomics and N-glycoproteomics. BMC Plant Biol 21, 384 (2021). https://doi.org/10.1186/s12870-021-03166-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-021-03166-4

Keywords