All relevant data are within the paper.

Introduction {#sec001}
============

Tea (*Camellia sinensis* (L.) O. Kuntze) is an important commercial crop that is grown mostly in the tropical and subtropical countries of Asia, Africa and to some extent Latin America \[[@pone.0177816.ref001]\]. In recent years, tremendous attention has been focused on tea plant due to its pleasant flavor and bioactive substances. Tea leaves are a valuable source of secondary metabolic products, including flavonoids, alkaloids, polysaccharides and theanine \[[@pone.0177816.ref002]\]. Among these secondary metabolites, flavonoids (which comprise polyphenols, flavones, flavanonols and anthocyanins) are considered to contribute numerous pharmacological beneficial effects on human health \[[@pone.0177816.ref003]--[@pone.0177816.ref004]\]. The health benefits of tea are thought to account for tea's protective role against cardiovascular disease \[[@pone.0177816.ref003]\], atherosclerosis \[[@pone.0177816.ref004]\], oxidant activity \[[@pone.0177816.ref005]\], and cancer \[[@pone.0177816.ref006]\].

The color of tea leaves has been diversified through long-time natural hybridization and artificial selection. Study on purple leaves of tea plant will positively promote the diversification of tea products and the fully utilization of different tea germplasm resources. Furthermore, tea products made from purple tea leaves are highly preferred by consumers. Compared with conventional tea, the anthocyanin-rich purple tea has multiple benefic functions including strong antioxidant activity \[[@pone.0177816.ref007]\], inhibition of colorectal carcinoma cell proliferation \[[@pone.0177816.ref008]\], and reinforcing brain antioxidant capacity \[[@pone.0177816.ref009]\]. In addition, anthocyanins can be used as commercial food colorants and have been reported to exhibit antioxidants, chemopreventive, anti-bacterial, antiangiogenic, anti-inflammatory and anti-atherosclerotic properties \[[@pone.0177816.ref010]\]. Therefore, to enhance the health potency of tea, purple-colored leaves have become one of the main qualitative attributes targeted in tea breeding programs. New tea cultivars with purple leaves have recently been developed and released in Kenya, China, Japan and India \[[@pone.0177816.ref011], [@pone.0177816.ref012]--[@pone.0177816.ref014]\].

The contents of total phenolic compounds and total anthocyanins, including catechins were higher in purple tea leaves than in green tea ones \[[@pone.0177816.ref015], [@pone.0177816.ref016]\]. Previous study has found that malvidin is the most abundant anthocyanidin in the tea products derived from the purple colored tea clones \[[@pone.0177816.ref011]\]. Anthocyanins often accumulate in young leaves of plants (including tea plants), and make leaves present a purple color, then the purple color in leaves gradually turn into green as the leaves mature \[[@pone.0177816.ref017]\]. Additionally, anthocyanin in young leaves correlates with tolerance to photoinhibition by the attenuation of visible radiation. Maayan \[[@pone.0177816.ref018]\] have reported that ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco), as well as chlorophyll, is not present in early development stages and only slowly accumulates throughout the developmental stages of the leaves. The lower abundance of rubisco is closely related to the rate of photosynthesis in young purple leaves \[[@pone.0177816.ref019]\]. In previous study, we used cDNA-AFLP to investigate the gene expression profiles associated with changes in leaf color during the developmental stages of tea leaves, and 148 differentially expressed genes, which were involved in cell respiration, nucleic acid and protein metabolism, photosynthesis, cell expansion, biosynthesis of flavonoids and anthocyanins, were examined \[[@pone.0177816.ref019]\]. However, the identification of differentially expressed genes at mRNA level is insufficient for the prediction of protein expression levels, and the available proteomics information on the changes in metabolites during shoot development in purple-leaf-tea cultivars is limited. Therefore, in this study, we used the newly developed tea cultivar Wuyiqizhong18 as source material, and we conducted 2-DE to separate differentially abundant proteins at the proteomic level and characterize the molecular mechanisms underlying changes in leaf color during the growth of tea plants with purple young shoots.

Materials and methods {#sec002}
=====================

Plant materials {#sec003}
---------------

Tea plants (*C*. *sinensis* (L.) O. Kuntze cv.Wuyiqizhong18) were grown in the Germplasm Repository of Tea Research Institute in Wuyi Mountain, Fujian province, China. Only the healthy, tender shoots (one bud and two leaves) were sampled on April 13, 2015, and then the purple and green leaves were packaged separately in aluminum foil. The plant characteristics in 'Wuyiqizhong18' tea cultivar were as follows: the bud were green, the color of the second leaf changes gradually to purple, and the leaves finally turned green again ([Fig 1](#pone.0177816.g001){ref-type="fig"}). The samples were immediately frozen in liquid nitrogen and stored at −80°C until use.

![The tender purple (TPL) and mature green leaves (MGL) of Wuyiqizhong18 tea plant.](pone.0177816.g001){#pone.0177816.g001}

Measurement of polyphenolic compounds {#sec004}
-------------------------------------

Anthocyanins were extracted with 1% HCl-methanol solution and measured using a spectrophotometer at 530 nm, according to method described by Hatlestad et al. \[[@pone.0177816.ref020]\]. The contents of total polyphenols were determined according to the International Organization for Standardization (ISO) 14502--1:2005 using Folin---Ciocalteu method \[[@pone.0177816.ref016]\]. The identification and quantification of catechins and caffeine were measured according to the method described by Wei et al. \[[@pone.0177816.ref021]\].

Protein extraction {#sec005}
------------------

Total protein was extracted from frozen leaf samples using a phenol extraction procedure according to Saravanan et al. \[[@pone.0177816.ref022]\] with some modifications. Briefly, about 1 g leaf sample was ground into fine powder in liquid N~2~ with a pestle and mortar. Approximately, 4 mL of ice-cooled extraction buffer \[100 mM Tris-HCl (pH 7.8), 50 mM L-ascorbic acid, 1% (v/v) Triton X-100, 1% (v/v) β-mercaptoethanol, 100 mM KCl, and 1 mM phenylmethylsulfonyl fluoride\] was added to the powder. The mixture was allowed to thaw slowly at room temperature and then gently pulverized for another two minutes. The homogenate was transferred to 10 ml centrifuge tube and centrifuged at 13,000 g for 15 min at 4°C after equal volume of Tris-saturated phenol (pH 8.0) was added. The phenolic phase was transferred to a 50 mL corning centrifuge tube with five volumes of methanol (with 100 mM ammonium acetate). The protein pellets were collected and re-suspended in 25 mL of ice-cooled pure methanol for 2 h at -20°C after centrifugation at 13,000 g for 10 min at 4°C. Protein pellets were collected by centrifugation at 13,000 g for 15 min at 4°C and washed twice with 25 mL of ice-cooled acetone, and then lyophilized at -20°C until use. About 10 mg protein pellet was dissolved in lysis buffer \[8 M urea, 1.5 M Thiourea, 4% (w/v) CHAPS, 2% (v/v) ampholyte (3--10), 15 mM DTT\]. The total protein concentration was measured by using Bradford protein assay \[[@pone.0177816.ref023]\].

2-DE and image analysis {#sec006}
-----------------------

Protein samples (about 1.5 mg) were rehydrated with a rehydration buffer containing 8 M urea, 2 M thiourea, 4% (w/v) 3-\[(3-cholamidopropyl)-dimethylammonio\] propanesulfonate (CHAPS, 0.5% (v/v) IPG buffer, 13 mM dithiothreitol (DTT), and 0.002% (w/v) bromophenol blue, then loaded onto a 24 cm, pH 4--7 linear immobilized pH gradient (IPG) strip (GE Healthcare, Uppsala, Sweden) on the rehydration cassette overnight. Isoelectric focusing electrophoresis (IEF) was conducted according to the methods described by You et al. \[[@pone.0177816.ref024]\].

Image analysis was carried out with PDQuest software version 8.0.1 (Bio-Rad, Hercules, CA, USA). Spot detection, matching and normalization were performed automatically, followed by manual verification. Only those candidate spots in all triplicate gels which had both a fold change of more than 1.5 and a *p* value of less than 0.05 by statistical test were considered differentially abundant proteins for mass spectrometry analysis.

Protein identification by MALDI-TOF/TOF-MS and bioinformatics analysis of proteins {#sec007}
----------------------------------------------------------------------------------

Protein spots were excised from the stained gels and digested with trypsin as described by You et al. \[[@pone.0177816.ref024]\]. Mass spectrometric analysis was conducted with a MALDI-TOF/TOF-MS 5800 (AB SCIEX, Shanghai, China). All acquired full-scan MS of samples were processed using TOF/TOF Explorer™ Software (AB SCIEX) in---default mode. The data were analyzed through the GPS Explorer software (Applied Biosystem, Foster City, CA, USA) using MASCOT software (Matrix Science Inc., MA, USA) with the database viridiplantae from NCBI website (1850050 sequences; 642453415 residues). The search criteria specified trypsin, allowing 1 missed cleavage, with a peptide tolerance of 100 ppm of and an MS/MS tolerance of 0.6 Da. Proteins exhibiting a score greater than 75 and statistical significance (*p* \< 0.05) were considered to be successfully identified proteins and were accepted.

Functional annotation of proteins was conducted using Gene Ontology (GO) annotation \[[@pone.0177816.ref025]\]. The differentially accumulated proteins were assigned to the Cluster of Orthologous Groups of proteins (COG) database and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database \[[@pone.0177816.ref026]\].

qPCR analysis {#sec008}
-------------

Total RNA extraction and qPCR were performed as previously described \[[@pone.0177816.ref019]\]. Specific primers were designed to investigate the gene expression at the transcript level of these representative differentially accumulated proteins. The sequences of the forward (F) and reverse (R) primers were listed in [S1 Table](#pone.0177816.s004){ref-type="supplementary-material"}. For the normalization of gene expression, *Camellia sinensis* 18S ribosomal RNA (AB120309.1) was used as a reference gene. Each treatment contained 3 biological replicates and each biological replicate had 3 technical replicates. The specificity of each primer pair was verified by melting curve analysis. The relative gene expression was calculated using the 2^-ΔΔCt^ algorithm below \[[@pone.0177816.ref027]\].

![](pone.0177816.e001.jpg){#pone.0177816.e001g}

Δ

Δ

C

t

=

(

C

t

g

e

n

e

o

f

i

n

t

e

r

e

s

t

−

C

t

i

n

t

e

r

n

a

l

c

o

n

t

r

o

l

)

t

r

e

a

t

m

e

n

t

−

(

C

t

g

e

n

e

o

f

i

n

t

e

r

e

s

t

−

C

t

i

n

t

e

r

n

a

l

c

o

n

t

r

o

l

)

c

o

n

t

r

o

l

Statistical analysis {#sec009}
--------------------

Significance tests were conducted using Student's t-test at *p* \< 0.05 level. The results are represented by the mean ± standard errors (SE).

To elucidate the overall proteomic responses of tea leaves, we performed a heatmap approach by using mean centered and standard deviation scaled with The Unscrambler (Version 10.0.1, Camo Software Inc., Woodbridge, NJ, USA) for the normalization of the values to give all variables an equal chance. After normalization of all the data sets, normalized values were used in the color scales system from Microsoft Office Excel 2010 (Microsoft Corporation, Redmond, WA, USA) to generate the heatmap.

The normalized data set was used for the principal component analysis (PCA) by performing normalized volume of the differentially abundant protein spots using the Unscrambler (version 10.0.1, Camo Software Inc., Woodbridge, NJ, USA). The PCA loading plots were used to determine the separation of the differentially expressed proteins based on the presence and absence between TPL and MGL samples. The PCA loading plots were performed in triplicate (n = 3).

Results {#sec010}
=======

Responses of polyphenolic compounds {#sec011}
-----------------------------------

In this study, we analyzed the contents of total anthocyanins (TA), total polyphenols (TP), catechins and caffeine in TPL and MGL tissues. As shown in [Fig 2](#pone.0177816.g002){ref-type="fig"}, the concentrations of TA and TP, (−)-epigallocatechin (EGC), (+)-catechin (C), (−)-epicatechin (EC), (−)-epigallocatechin 3-O-gallate (EGCG), (−)-epicatechin 3-O-gallate (ECG), total catechins and caffeine were significantly higher in TPL tissues than in MGL ones. These findings clearly indicate that the content of quality determinants in tea leaves decreases as tea leaves develop.

![Changes in the contents of total anthocyanins (TA), total polyphenols (TP), (−)-epigallocatechin (EGC), (+)-catechin (C), (−)-epicatechin (EC), (−)-epigallocatechin 3-O-gallate (EGCG), (−)-epicatechin 3-O-gallate (ECG), total catechins (TC) and caffeine (CAF) in mature green leaf (MGL) and tender purple leaf (TPL) tissues.\
Bars represent means ± SE (n = 4). Different letters above the bars indicate a significant difference at *p*\<0.05.](pone.0177816.g002){#pone.0177816.g002}

Identification of differentially expressed proteins {#sec012}
---------------------------------------------------

To explore the correlation between the proteomic and metabolite profiles of TPL and MGL tissues, 2-DE coupled with MS/MS technique was carried out with three biological replicates. More than 900 clear and reproducible spots were detected on each gel ([Fig 3](#pone.0177816.g003){ref-type="fig"}). After conducting analyses with PDQuest 8.0.1 software and MALDI-TOF/TOF-MS, we successfully identified 46 proteins ([S1 Fig](#pone.0177816.s001){ref-type="supplementary-material"}, [Table 1](#pone.0177816.t001){ref-type="table"}). Of these identified proteins, 26 proteins had high abundance in TPL tissues, while 20 proteins had low abundance compared with those in MGL tissues. The results from the database are presented in [Table 1](#pone.0177816.t001){ref-type="table"}.

![Representative 2-DE gels from *Camellia sinensis* of the development of young shoot purple-related tea plant.\
**A, mature green leaf (MGL) tissues; B, tender purple leaf (TPL) tissues**. Note: Differentially regulated proteins are numbered and indicated by arrows.](pone.0177816.g003){#pone.0177816.g003}

10.1371/journal.pone.0177816.t001

###### List of differentially expressed proteins identified using MALDI-TOF/TOF-MS in tender purple leaf (TPL) and mature green leaf (MGL) tissues of *Camellia sinensis*.

![](pone.0177816.t001){#pone.0177816.t001g}

  Spot No.[^a^](#t001fn001){ref-type="table-fn"}   Protein                                                         Accession No.[^b^](#t001fn002){ref-type="table-fn"}   Score   Expect     Species                               Mr      Pi     Ratio[^c^](#t001fn003){ref-type="table-fn"}±SD
  ------------------------------------------------ --------------------------------------------------------------- ----------------------------------------------------- ------- ---------- ------------------------------------- ------- ------ ------------------------------------------------
  ***Carbohydrate and energy metabolism***                                                                                                                                                                                                       
  2206                                             Oxygen evolving enhancer protein                                gi\|326467059                                         466     6.8e-041   *Litchi chinensis*                    39.25   5.35   0.31±0.02
  7209                                             Carbonic anhydrase                                              gi\|427199625                                         245     8.5e-019   *Vigna unguiculata*                   36.54   6.24   0
  8304                                             Cytosolic malate dehydrogenase                                  gi\|301015433                                         626     6.8e-057   *Camellia sinensis*                   48.71   6.4    100
  8306                                             Malate dehydrogenase                                            gi\|114479586                                         514     1.1e-045   *Citrus junos*                        46.03   6.45   100
  8308                                             Fructose-bisphosphate aldolase 3                                gi\|432139325                                         414     1.1e-035   *Camellia oleifera*                   47.34   6.47   0.49±0.05
  4304                                             Fructokinase                                                    gi\|402810391                                         290     2.7e-023   *Actinidia eriantha*                  47.12   5.62   2.27±0.18
  2409                                             Phosphoribulokinase, chloroplastic-like                         gi\|460400830                                         262     1.7e-020   *Solanum lycopersicum*                52.31   5.36   0.34±0.02
  6410                                             Phosphoglycerate kinase, chloroplastic-like                     gi\|356525742                                         741     2.1e-068   *Glycine max*                         53.4    6.06   0.37±0.01
  8702                                             Pyruvate decarboxylase family protein                           gi\|297319662                                         151     2.1e-009   *Arabidopsis lyrata subsp*.*lyrata*   69.17   6.33   2.39±0.01
  2403                                             Sedoheptulose-1,7-bisphosphatase, chloroplastic                 gi\|225466690                                         423     1.4e-036   *Vitis vinifera*                      49.92   5.21   0.37±0.04
  7705                                             Cytochrome b5 DIF-F                                             gi\|4204575                                           132     1.7e-007   *Petunia x hybrida*                   21.13   6.13   100
  3105                                             Cytochrome b6-f complex iron-sulfur subunit                     gi\|225461287                                         155     8.50E-10   *Vitis vinifera*                      25.52   5.48   0.31±0.01
  2608                                             ATP synthase CF1 alpha subunit (chloroplast)                    gi\|430728257                                         984     1.1e-092   *Camellia sinensis*                   65.8    5.32   0.44±0.01
  2603                                             ATP synthase CF1 alpha subunit (chloroplast)                    gi\|430728257                                         926     6.8e-087   *Camellia sinensis*                   66.56   5.2    0.22±0.03
  2605                                             ATP synthase CF1 alpha subunit (chloroplast)                    gi\|430728257                                         939     3.4e-088   *Camellia sinensis*                   65.7    5.26   0.36±0.02
  7311                                             ATP-dependent clp protease, putative                            gi\|223550217                                         395     8.50E-34   *Ricinus communis*                    47.61   6.14   0
  4401                                             Ribulose 1,5-bisphosphate carboxylase/oxygenase activase        gi\|359478916                                         320     2.7e-026   *Vitis vinifera*                      51.3    5.56   0.11±0.01
  2402                                             Ribulose 1,5-bisphosphate carboxylase/oxygenase activase        gi\|158726716                                         404     1.1e-034   *Flaveria bidentis*                   53.85   5.21   0.2±0.05
  8204                                             Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit   gi\|148590322                                         610     2.7e-055   *Parkia multijuga*                    39.21   6.46   0.2±0.02
  9502                                             Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit   gi\|18140571                                          662     1.7e-060   *Camellia sinensis*                   56.55   6.61   0.07±0.02
  8506                                             Ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit   gi\|18140563                                          578     4.30E-52   *Camellia albogigas*                  60.29   6.5    0.04±0.01
  7505                                             Ribulose 1,5-bisphosphate carboxylase                           gi\|17136048                                          683     1.4e-062   *Tetrameles nudiflora*                58.95   6.27   0.3±0.03
  5308                                             Carbon-nitrogen hydrolase family protein                        gi\|297317346                                         88      0.004      *Arabidopsis lyrata*                  48.08   5.86   3.33±0.25
  ***Secondary metabolism***                                                                                                                                                                                                                     
  4106                                             Chalcone isomerase                                              gi\|76152009                                          252     1.7e-019   *Camellia sinensis*                   33.98   5.69   2.02±0.17
  8406                                             Chalcone synthase                                               gi\|22086369                                          128     4.3e-007   *Rubus idaeus*                        52.87   6.47   3.13±0.25
  2405                                             Caffeine synthase                                               gi\|145952324                                         316     6.8e-026   *Camellia sinensis*                   51.14   5.26   2.53±0.1
  4412                                             Flavonol synthase                                               gi\|76786311                                          705     8.5e-065   *Camellia sinensis*                   50.17   5.71   2.43±0.19
  ***Protein metabolism***                                                                                                                                                                                                                       
  5601                                             Anthranilate N-benzoyltransferase protein, putative             gi\|223547579                                         105     8.5e-005   *Ricinus communis*                    63.79   5.74   0.33±0.03
  8701                                             TCP-1/cpn60 chaperonin family protein                           gi\|332006171                                         205     8.5e-015   *Arabidopsis thaliana*                69.81   6.31   2.43±0.1
  4407                                             Glutamine synthetase                                            gi\|42733460                                          570     2.70E-51   *Camellia sinensis*                   51.49   5.65   2.9±0.07
  6706                                             T-complex protein 1 subunit beta                                gi\|225459806                                         440     2.7e-038   *Vitis vinifera*                      66.93   6      3.08±0.25
  8202                                             proteasome subunit alpha type-6 isoform                         gi\|359482460                                         334     1.1e-027   *Vitis vinifera*                      35.66   6.44   2.02±0.14
  2607                                             alanine aminotransferase 2, mitochondrial-like, partial         gi\|449519802                                         150     2.7e-009   *Cucumis sativus*                     63.26   5.31   100
  801                                              Protein disulfide isomerase-like 1-4-like isoform 2             gi\|356539444                                         81      0.02       *Glycine max*                         76.99   4.85   2.22±0.64
  ***Nucleic acid metabolism***                                                                                                                                                                                                                  
  5109                                             Glycine-rich RNA-binding protein                                gi\|187373099                                         323     1.4e-026   *Nicotiana tabacum*                   19.97   5.86   100
  8704                                             Bifunctional purine biosynthesis protein, putative              gi\|223529724                                         162     1.70E-10   *Ricinus communis*                    67.22   6.43   5±0.48
  ***Cell wall*, *cytoskeleton and transport***                                                                                                                                                                                                  
  2504                                             alpha-tubulin                                                   gi\|385717688                                         892     1.7e-083   *Oryza sativa Japonica Group*         60.71   5.32   100
  1607                                             Tubulin, beta chain                                             gi\|222872170                                         721     2.1e-066   *Populus trichocarpa*                 62.8    5.09   2.39±0.28
  504                                              Bark storage protein A                                          gi\|225429975                                         116     6.8e-006   *Vitis vinifera*                      56.4    4.85   2.24±0.21
  ***Stress responses***                                                                                                                                                                                                                         
  5116                                             Copper/zinc superoxide dismutase 1                              gi\|381283808                                         373     1.4e-031   *Litchi chinensis*                    18.98   5.9    0.48±0.04
  6204                                             Thioredoxin-like protein CDSP32, chloroplastic                  gi\|225459760                                         324     1.1e-026   *Vitis vinifera*                      38.33   6      0.46±0.03
  4107                                             Cold shock domain protein                                       gi\|187609563                                         152     1.7e-009   *Oryza sativa Japonica Group*         24.74   5.69   2.72±0.04
  5106                                             Glutathione S-transferase                                       gi\|76365795                                          124     1.1e-006   *Citrus sinensis*                     31.77   5.82   100
  5114                                             Phospholipid hydroperoxide glutathione peroxidase               gi\|34786892                                          400     2.7e-034   *Citrus sinensis*                     27      5.82   2.05±0.07
  104                                              Translationally controlled tumor protein                        gi\|158120965                                         168     4.3e-011   *Salvia miltiorrhiza*                 27.22   4.71   3.58±0.74
  ***Other and unknown biological processes***                                                                                                                                                                                                   
  8109                                             predicted protein                                               gi\|222857852                                         63      1.1e-005   *Populus trichocarpa*                 26.83   6.39   100

^a^: Spot number corresponds to the 2-DE gel in [Fig 3](#pone.0177816.g003){ref-type="fig"}.

^b^: gi number is from NCBI database of matched protein.

^c^: Ratio means the ratio of TPL to MGL; 0 means protein spots were only detected in MGL tissues; 100 means protein spots were only detected in the TPL tissues.

According to biological functional properties, the 46 differentially abundant proteins were categorized into the following functional groups: carbohydrate and energy metabolism (50%), protein metabolism (15%), secondary metabolism (9%), stress responses (13%), nucleic acid metabolism (4%), cell wall, cytoskeleton and cell transport (7%), other and unknown biological processes (2%) ([Fig 4](#pone.0177816.g004){ref-type="fig"}).

![Functional classification of differentially expressed transcript derived fragments.](pone.0177816.g004){#pone.0177816.g004}

PCA and heatmap responses {#sec013}
-------------------------

PCA was performed with all 46 of the identified proteins. The PC1 and PC2 values, which accounted for 89% and 5% of the total X and Y variance, revealed that the differential protein responses between TPL and MGL were clearly separated by the first principal component ([Fig 5A](#pone.0177816.g005){ref-type="fig"}). Interestingly, the majority of differentially expressed proteins involved in carbohydrate and energy metabolism were highly clustered ([Fig 5B](#pone.0177816.g005){ref-type="fig"}). The heatmap of protein responses in TPL and MGL samples shows the correlation patterns of the differentially abundant proteins ([Fig 6](#pone.0177816.g006){ref-type="fig"}).

![Principal component analysis (PCA) of scores (A) and loading (B) plots of proteomic responses from mature green leaf (MGL) and tender purple leaf (TPL) tissues of tea plants.](pone.0177816.g005){#pone.0177816.g005}

![Heatmap of proteomic responses from mature green leaf (MGL) and tender purple leaf (TPL) tissues of tea plants.\
Proteomic levels correspond to the color scale. A color gradient from blue through to red represents a low level through to a high level of proteomic response. After calculating the mean of each metabolite, the normalized data were generated using the Unscrambler software.](pone.0177816.g006){#pone.0177816.g006}

Transcriptional expression analysis of differentially abundant proteins {#sec014}
-----------------------------------------------------------------------

To verify the proteomic results and evaluate the association between mRNA and protein levels, qPCR assays were conducted ([Fig 7](#pone.0177816.g007){ref-type="fig"}). From all identified proteins, fifteen gene products were selected to analyze the expression patterns of their mRNA levels: flavonol synthase (spot No. 412), chalcone synthase (spot No. 8406), ATP synthase CF1 alpha subunit (spot No. 2603), chalcone isomerase (spot No. 4106), malate dehydrogenase (spot No. 8306), chalcone synthase (spot No. 8406) and rubisco (spot No. 9502), fructokinase (spot No. 4304), cytosolic malate dehydrogenase (spot No. 8304), glutamine synthetase (spot No. 4407), alpha-tubulin (spot No. 2504), glutathione *S*-transferase (spot No. 5106), copper/zinc superoxide dismutase (spot No. 5116), alanine aminotransferase (spot No. 2607) and oxygen evolving enhancer protein (spot No. 2206). The expression levels of genes encoding flavonol synthase ([Fig 7A](#pone.0177816.g007){ref-type="fig"}), caffeine synthase ([Fig 7B](#pone.0177816.g007){ref-type="fig"}), chalcone isomerase ([Fig 7D](#pone.0177816.g007){ref-type="fig"}), malate dehydrogenase ([Fig 7E](#pone.0177816.g007){ref-type="fig"}), chalcone synthase ([Fig 7F](#pone.0177816.g007){ref-type="fig"}), fructokinase ([Fig 7H](#pone.0177816.g007){ref-type="fig"}), cytosolic malate dehydrogenase ([Fig 7I](#pone.0177816.g007){ref-type="fig"}), glutamine synthetase ([Fig 7J](#pone.0177816.g007){ref-type="fig"}), alpha-tubulin ([Fig 7K](#pone.0177816.g007){ref-type="fig"}), glutathione *S*-transferase ([Fig 7L](#pone.0177816.g007){ref-type="fig"}) and alanine aminotransferase ([Fig 7N](#pone.0177816.g007){ref-type="fig"}) were higher in TPL tissues than in MGL ones, while the expression levels of genes encoding the ATP synthase CF1 alpha subunit ([Fig 7C](#pone.0177816.g007){ref-type="fig"}), rubisco ([Fig 7G](#pone.0177816.g007){ref-type="fig"}), copper/zinc superoxide dismutase ([Fig 7M](#pone.0177816.g007){ref-type="fig"}) and oxygen evolving enhancer protein ([Fig 7O](#pone.0177816.g007){ref-type="fig"}) were lower in TPL tissues than in MGL ones. As shown in [Fig 7](#pone.0177816.g007){ref-type="fig"}, the transcription profiles of these genes matched the protein expression levels well, indicating that the levels of differentially abundant proteins were regulated at the transcriptional level. Only one gene encoding oxygen evolving enhancer protein exhibited different expression patterns under qPCR analysis compared with proteomic evaluation. The different results may be attributed to the posttranscriptional regulatory process and the different sensitivities of the analytical methods.

![Relative expression levels of genes encoding (A) flavonol synthase (S412), (B) caffeine synthase (S2405), (C) ATP synthase CF1 alpha subunit (S2603), (D) chalcone isomerase (S4106), (E) malate dehydrogenase (S8306), (F) chalcone synthase (S8406), (G) ribulose-1,5-bisphosphate carboxylase (S9502), (H) fructokinase, (I) cytosolic malate dehydrogenase, (J)glutamine synthetase, (K) alpha-tubulin, (L) glutathione *S*-transferase, (M) copper/zinc superoxide dismutase, (N) alanine aminotransferase, (O) oxygen evolving enhancer protein from *Camellia sinensis* leaves revealed by qPCR.\
Bars represent means ± SE (n = 9). Different letters above standard error bars indicate significant differences at *p* \< 0.05.](pone.0177816.g007){#pone.0177816.g007}

Discussion {#sec015}
==========

Anthocyanins often accumulate in young leaves, fruit or flowers, and then degrade as leaves develop \[[@pone.0177816.ref028]\]. This tendency has been reported on many plant species, including *Chrysobalanus icaco* \[[@pone.0177816.ref029]\], *Capsicum* spp \[[@pone.0177816.ref030]\] and *Brunfelsia calycina* \[[@pone.0177816.ref031]*\]*. A previous study by our group also showed that anthocyanin contents decreased with leaf growth and development, while the contents of chlorophyll *a* and chlorophyll *b* were higher in MGL samples. Thus, we used the cDNA-AFLP technique to reveal the regulatory mechanism of the alterations of leaf color during leaf development in Wuyiqizhong18 \[[@pone.0177816.ref019]\]. Nevertheless, due to post-transcriptional regulation, gene transcription data are insufficient to predict the expression levels of proteins. To further understanding the molecular mechanism underling the development of tea leaves, we used 2-DE technique to investigate the proteomic profile of different leaves and successfully identified 46 differentially abundant proteins in MGL and TPL samples in this study.

Proteins involved in carbohydrate and energy metabolism {#sec016}
-------------------------------------------------------

We identified 6 increased and 17 decreased proteins that were involved in carbohydrate and energy metabolisms from TPL tissues ([Table 1](#pone.0177816.t001){ref-type="table"}). In plants, the tricarboxylic acid cycle (TCA cycle) and the pentose phosphate pathway (PPP) provide not only energy and cofactors but also substrates for the synthesis of metabolites, or signals for feedback \[[@pone.0177816.ref032]\]. Interestingly, in the present study, we found that some key enzymes involved in carbon assimilation had high abundance in TPL, such as phosphoglycerate kinase (PGK, spot No. 6410), fructose-bisphosphate aldolase (spot No. 8308), sedoheptulose-1,7-bisphosphatase (spot No. 2403) and phosphoribulokinase (PRK, spot No. 2409) ([Fig 8](#pone.0177816.g008){ref-type="fig"}). Moreover, TPL sample has a lower abundance of rubisco (spot No. 4401, 2402, 8204, 9502, 8506 and 7505), a bi-functional enzyme that catalyzes the initial step of the Calvin cycle and the initial oxygenation reaction of photorespiration. This observation was well corresponded with our previous study, reporting that the transcriptional levels of the genes encoding some photosynthetic enzymes are decreased in TPL tissues \[[@pone.0177816.ref019]\]. The lower abundance of rubiscos would be responsible for the lower activity of CO~2~ assimilation in TPL tissues than in MGL ones. It is also possible that the lower abundance of rubiscos in TPL tissues due to the immature chloroplasts of TPL tissues \[[@pone.0177816.ref018]\].

![The differentially abundant proteins in the light reactions and Calvin cycle in tender leaves of *Camellia sinensis*.\
Note: Rubisco: Ribulose-1,5-bisphosphate carboxylase/oxygenase; FBA: Fructose-bisphosphate aldolase; SBPase: Sedoheptulose-1,7-bisphosphatase; PRK: Phosphoribulokinase; Arrows represent meansLow abundance.](pone.0177816.g008){#pone.0177816.g008}

Four enzymes (malate dehydrogenase, spot No. 8304 and 806; fructokinase, spot No. 4304; pyruvate decarboxylase, spot No. 8702) involved in organic acid metabolism and glycolysis have significantly higher abundance in TPL tissues than in MGL ones. It means that energy and secondary metabolism would be enhanced in TPL tissues. Previous research has suggested that anthocyanins in young leaves are likely associated with tolerance against photoinhibition by attenuating visible radiation \[[@pone.0177816.ref033], [@pone.0177816.ref034]\]. Since TPL tissues are rich in anthocyanins, light attenuation by anthocyanins may display shade acclimation of the photosynthetic machinery in purple leaves, which is compatible with a photoprotective function of anthocyanins \[[@pone.0177816.ref035]\]. The accumulation of anthocyanins in TPL tissues may serve as a filter of light and protect chloroplasts against photoinhibition \[[@pone.0177816.ref036]\]. Thus, anthocyanins likely play a key role in the photoprotection of young leaves, which is a hypothesis that has been addressed in a previous study \[[@pone.0177816.ref037]\].

Proteins involved in secondary metabolism {#sec017}
-----------------------------------------

Although the proportion of differentially abundance proteins related to secondary metabolism in our results was relatively low (9%) compared with the other categories, proteins related to anthocyanin and caffeine biosynthesis were found to be differentially abundances between TPL and MGL tissues ([Table 1](#pone.0177816.t001){ref-type="table"}). Tea leaves are a valuable source of secondary metabolites, including polyphenols, alkaloids, volatile oils, and others \[[@pone.0177816.ref002], [@pone.0177816.ref038]\]. Based on the proteomic results in this study, several proteins related to polyphenol biosynthesis were differentially expressed between the TPL and MGL tissues. CHS is the first committed enzyme involved in anthocyanin and flavonoid biosynthesis \[[@pone.0177816.ref039]\]. Gene silencing (via RNA interference or virus-induced gene silencing) of CHS in both of *Torenia hybrida* and petunia resulted in lower contents of total flavonoid and anthocyanins, and CHS enzyme activity of flower tissues \[[@pone.0177816.ref039], [@pone.0177816.ref040]\]. In the present study, the abundance of CHS (spot No. 8406) was increased in TPL tissues, indicating that the biosysnthesis of anthocyanin and flavonoid was enhanced. This outcome was consistent with our previous finding that CHS expression was a positive regulator of anthocyanin biosynthesis \[[@pone.0177816.ref019]\]. Furthermore, qPCR was employed to verify the higher gene expression level of CHS in TPL tissues ([Fig 7F](#pone.0177816.g007){ref-type="fig"}). CHI and FLS are two key enzymes involved in the flavonoid biosynthesis pathway \[[@pone.0177816.ref041]--[@pone.0177816.ref042]\]. In CHI-suppressed transgenic tobacco plants, they showed up to 25% reduction in anthocyanins content \[[@pone.0177816.ref043]\]. In the present study, the abundance of CHI (spot No. 4106) and FLS (spot No. 412) were higher in TPL tissues, which was consistent with the qPCR results, than in MGL ones ([Fig 7](#pone.0177816.g007){ref-type="fig"}). Thus, higher abundance of CHI, CHS and FLS may account for the production of the purple-shoot phenotype in Wuyiqizhong 18 and enhance the flavonoid biosynthesis. These results were also agreed with our metabolic data, which show that the contents of TA and TP were significantly higher in TPL tissues than in MGL ones ([Fig 2](#pone.0177816.g002){ref-type="fig"}).

Proteins involved in stress responses {#sec018}
-------------------------------------

The production of reactive oxidative species (ROS) is inevitable in plant grown in aerobic circumstance. As an adaptive strategy, plants have evolved several capabilities to scavenge ROS including endogenous antioxidant enzymes system and antioxidant products, such as superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPx), and glutathione reductase (GR). Due to different sub-cellular localization, the abundance of these anti-oxidative enzymes did not change in the same pattern in response to an environmental stimulus \[[@pone.0177816.ref044]\]. In our study, four proteins involved in stress responses have higher abundance in TPL tissues, while one protein named copper/zinc superoxide dismutase (Cu/Zn SOD, spot No. 5116) has lower abundance in TPL tissues than in MGL ones ([Table 1](#pone.0177816.t001){ref-type="table"}; [Fig 7M](#pone.0177816.g007){ref-type="fig"}). Cu/Zn SOD catalyzes the dismutation of ion superoxide into oxygen and hydrogen peroxide in plant chloroplast \[[@pone.0177816.ref045]\]. The lower abundance of Cu/Zn SOD in TPL tissues might be due to the well protection of anthocyanins. PHGPx, one of the protein isoforms of GPx, protects plants from pathogen attack and tolerance to other abiotic stresses \[[@pone.0177816.ref046]\]. A previous study suggested that only part of the absorbed light energy is utilized in photosynthesis, and the excess photon flux could lead to photo-damage and ROS production due to decreased CO~2~ assimilation in TPL tissues \[[@pone.0177816.ref019], [@pone.0177816.ref047]\]. As hypothesis, the increased abundance of PHGPx (spot No. 5114) and GST (spot No. 5106) observed in TPL tissues could quench the generated ROS by excess light and ameliorate the oxidative damage caused by ROS. Moreover, GPx and GSTs with GPx activity also decompose other toxic compound like alkyl hydroperoxides in addition to H~2~O~2~ generated by SODs activity \[[@pone.0177816.ref048]\]. Based on these results, we concluded that the accumulation of stress-related proteins in TPL tissues might dissipate the excess excitation energy and protect tender leaves against photo-damage.

Proteins involved in protein metabolism {#sec019}
---------------------------------------

Most of the proteins involved in protein metabolism have higher abundance in TPL tissues, compared with those in MGL tissues. In this study, the higher abundance of protein disulfide isomerase-like 1-4-like isoform 2 (PDIL, spot No.801), glutamine synthetase (GS, spot No.4407) and chaperonin family proteins (spot No.8701) indicated that there was an increased requirement for protein synthesis, correct assembly and turnover due to the development of tender leaves. Especially, GS (EC 6.3.1.2) is a key enzyme involved in the assimilation of ammonium into glutamine (Gln) \[[@pone.0177816.ref049], [@pone.0177816.ref050]\]. GS plays an important role in the synthesis of theanine, which is a unique amino acid found in tea plants and is the main component responsible for the taste of tea \[[@pone.0177816.ref051]\]. GS showed higher in tender leaves and lower abundance in older leaves, indicating that the synthesis of theanine in tender leaves was enhanced \[[@pone.0177816.ref019], [@pone.0177816.ref052]\].

Proteins involved in nucleic acid metabolism, the cellular cytoskeleton and transport {#sec020}
-------------------------------------------------------------------------------------

Glycine-rich RNA-binding proteins (GRPs) harbor an RNA-recognition motif (RRMs) at the N-terminus and a glycine-rich region at the C-terminus \[[@pone.0177816.ref053]\]. The expression patterns of GRPs are regulated by a number of external stimuli, such as cold stress conditions \[[@pone.0177816.ref054]\], high salinity and osmotic stress \[[@pone.0177816.ref055]\]. The high abundances of GRPs (spot No.5109) in TPL tissues could enhance stress resistance of tender tissue. Microtubules are dynamic heteropolymers of α- and β-tubulin, which are structural component of the cytoskeleton in eukaryote \[[@pone.0177816.ref056]\]. It plays key roles in many basic processes of eukaryotic cells, such as cell division, cell motility, intracellular transport and the control of cell shape \[[@pone.0177816.ref057]\]. The high abundance of α-tubulin (spot No. 2504) and β-tubulin (spot No. 1607) could maintain the stability of the cytoskeleton system, which is very important for the development of plants.

Conclusions {#sec021}
===========

Tea plants with purple coloration represent important tea germplasms due to their economic value and health benefits. The pronounced changes in the proteomic profile between TPL and MGL tissues were included: 1) the lower activity of proteins associated with CO~2~ assimilation, energy metabolism and photo flux efficiency and higher content of anthocyanins in TPL samples may protect tender leaves against photo-damage; 2) the higher abundance of CHS, CHI and FLS likely contributes to the synthesis of anthocyanins, catechins and flavonols in TPL tissues; 3) higher abundance of stress response proteins, such as GST and PHGPx, could enhance the tolerance of TPL tissues to adverse condition; and 4) the increased abundance of proteins related to protein synthesis, nucleic acids and cell wall proteins should be beneficial for the proliferation and expansion of leaf cell in TPL tissues. In summary, the results revealed that the higher abundance of two key proteins including CHS and CHI were responsible for the higher content of anthocyanins in TPL tissues compared with MGL ones. The higher abundance of GS protein related to the theanine biosynthesis may improve the flavor of tea products from TPL materials. This work enhances our understanding of the molecular mechanisms behind the changes in leaf color alteration between TPL and MGL tissues. Therefore, future studies should focus on the regulatory role of CHS and CHI in anthocyanin biosynthesis pathway of purple leaves of tea plants.

Supporting information {#sec022}
======================

###### Total ion chromatogram of protein spot 8406 excised from a 2-DE gel, trypsin-digested and analyzed by MALDI-TOF/TOF-MS.

\(A\) Peptide mass fingerprinting of protein spot 8406; (B) MS/MS spectra of ion 1799.8, the corresponding peptide sequence are shown.

(TIF)

###### 

Click here for additional data file.

###### Replicated gels of mature green leaf (MGL) and tender purple leaf (TPL) tissues.

(DOCX)

###### 

Click here for additional data file.

###### Master gel.

(PDF)

###### 

Click here for additional data file.

###### Primer sequences used for quantitative real-time PCR.

(DOCX)

###### 

Click here for additional data file.

We thank Irene Lavagi and Christian Craddock of Plant Bio Press Editing LLC for editing our manuscript.

[^1]: **Competing Interests:**The authors have declared that no competing interests exist.

[^2]: **Conceptualization:** QZ WS.**Data curation:** QZ JL XL.**Formal analysis:** WS XL.**Investigation:** QZ ZC.**Methodology:** QZ.**Project administration:** WS XL.**Resources:** WS XL.**Software:** QZ JL.**Supervision:** WS ZC.**Validation:** QZ.**Visualization:** QZ JL.**Writing -- original draft:** QZ JL.**Writing -- review & editing:** QZ XL WS.
