BACKGROUND
the field of protein discovery through mass spectrometry  continues to grow rapidly but the number of species for which finished  whole genome sequence data are available is currently not keeping pace.
for a large number of laboratories worldwide studying proteomes in 'non-mainstream' organisms, annotations of tandem mass spectra data must rely on open reading frame  predictions from expressed sequence tag  data from their species of interest or a phylogenetically close relative.
ests generated from single pass sequencing reactions are frequently not full length and the reading frames are unknown.
proteolytic peptide sequence databases derived from multiple, truncated predicted orfs per each of thousands of ests can hamper the ability of search engines such as mascot  <cit>  and algorithms such as paragon in proteinpilot software  <cit>  to make statistically robust protein identifications from ms/ms spectrum data  <cit> .
protein identifications from ms/ms spectra may be further complicated when the est data that are used to build a peptide sequence database are created based on one genotype for a given species.
we report here on the development of scripts for the generation of a predicted tryptic peptide sequence database based on est data in grapevine.
our computational approach accounts for multiple open reading frames, truncated predicted orfs, and the presence of n-terminal signal peptides, and may be useful for ms/ms-based protein discovery in any species for which est data are available.
quantitative protein expression profiling analyses in plants have increasingly implemented stable isotopic labeling as an advance or complement to two dimensional gel electrophoresis  methods.
isotope coded affinity tagging  reagents are used to covalently label cysteine residues with heavy or light hydrogen or carbon in two complex peptide samples, for example, wild type versus mutant genotypes.
the icat chemistry is used to purify labeled peptides via affinity chromatography and then samples are mixed and subjected to lc-ms/ms  <cit> .
one of the first reports on an icat application in plants was in wheat  where relative expression in monosomic deletion mutants was used to begin to clarify the influence of ancestral genomes on differential seed protein expression for breeding applications  <cit> .
the icat technique is limited, however, by the tagging of cysteine residues only, as well as the need for affinity purification of labeled peptides; invariably, information is lost through these steps.
an improvement to the icat technique involves the labeling of amine groups using a set of four or more isobaric tags.
the advantages of this technique, isobaric tagging for relative and absolute quantitation , are that most peptides are labeled, no affinity purification step is required, and the isobaric nature of the tags allows co-elution of identical peptides that are differentially tagged, thereby enriching detection sensitivity and accuracy in comparison to icat  <cit> .
few reports of itraq implementation in plant proteome studies have been reported but pioneering work in this field has been successful, for example, in further defining the organellar proteome in arabidopsis  <cit> , characterizing pathogen defense mechanisms in arabidopsis  <cit> , and clarifying micronutrient stress responses in barley   <cit> .
we were interested in characterizing ripening initiation in grape berries at the level of differential protein expression in order to better define the molecular control of this important process for grape growers and wine makers.
grape berry ripening is non-climacteric and ethylene does not act as a major signal initiating this process, as it does in climacteric species such as tomato .
abscisic acid , hexoses, and brassinosteroids  have previously been implicated in non-climacteric ripening regulation but how these and potentially other signaling pathways interact to effect major changes in berry biochemistry at ripening initiation is poorly understood.
the tissues in the grape berry consist of the seeds, the mesocarp , and the exocarp ; the pericarp refers to the mesocarp and exocarp, collectively.
primary and secondary compounds important for grape and wine products begin to accumulate in the exocarp and/or the mesocarp at ripening initiation  <cit> , so we considered that it was important to evaluate changes in the berry proteome separately in these tissues.
to date, a limited number of reports on proteome profiling in grapevine and grape berries have been published in which 2dge was employed  <cit> .
we considered that the itraq technique could be useful in surmounting some technical limitations encountered with 2dge and allow us to detect a greater number of proteins per sample.
in this report, we demonstrate the application of our computational approach to tryptic peptide sequence database development from a large collection of grapevine est data and validate its usefulness by showing improved detection and annotations of ms/ms data derived from grape exocarp and mesocarp total protein extracts.
we further provide new quantitative information on differential protein expression during ripening initiation in grape berries.
this is the first report in which itraq has been used to study differential protein expression in any fruit.
methods
plant material
grape clusters were sampled from v. vinifera cv.
cabernet sauvignon clone  <dig> grafted on rootstock 101- <dig> in a commercial vineyard near osoyoos, british columbia, in the  <dig> and  <dig> seasons.
sampling dates during each season were focused on the developmental stages undergoing ripening initiation.
clusters were sampled on a single date in  <dig>  august 12th, which was the timing of approximately 50% ripening initiation based on a turning pink color phenotype.
for the  <dig> season, the ripening initiation stage was sampled over a longer period , since in this growing season, ripening advanced slowly due to lower atmospheric temperatures.
five clusters from five different vines were sampled in each season and snap frozen directly in liquid nitrogen in the vineyard and then transported on dry ice to ubc vancouver where they were stored at -80°c.
individual grapes from each of the  <dig> and  <dig> clusters were developmentally staged based on a visual pigmentation assessment and were segregated for each season into green, pink/turning, fully turned red, and fully turned purple phenotypic classes.
for the  <dig> samples, green grapes were only taken from clusters collected on august 10th, since for this date and august 12th there was no visible change in color present in any of the grape clusters.
thirty grapes of similar sizes per pigmentation class per year were segregated for experimentation.
prior to total protein extraction, individual grapes were partially thawed in gloved hands and then, using a forceps, the exocarp tissue was carefully peeled away from the mesocarp and placed immediately into liquid nitrogen.
seeds were then carefully removed while keeping the remaining mesocarp tissue frozen in liquid nitrogen.
exocarp and mesocarp samples were ground to a powder under liquid nitrogen and then used for total protein extractions.
tissue preparation for protein extraction
preparation of exocarp tissue samples for protein extraction was performed according to a previously described protocol for olive leaf  <cit>  with some modifications described here.
the procedure was carried out on ice and centrifugations were performed at 4°c.
throughout the procedure, each wash was done by complete resuspending of the tissue pellet.
four hundred mg of powdered exocarp tissue was placed in a  <dig> ml g-tube .
the tissue was suspended in  <dig>  ml of a cold  ethyl acetate:ethanol ) solution by vortexing for  <dig> s; the ethyl acetate:ethanol extraction was previously found to be useful for removing pectins as well as pigments such as chlorophylls  <cit> .
following centrifugation for  <dig> min at  <dig> × g, the supernatant was removed and the ethyl acetate:ethanol extraction and centrifugation steps were repeated on the remaining tissue.
the sample was next extracted twice with cold  100% acetone by vortexing and centrifuging, as before.
subsequently, the tissue with added acetone was transferred from the g-tube to a mortar using a  <dig> ml pipette with the tip end excised to increase diameter and then the acetone was evaporated from the tissue at room temperature.
after the addition of 1/ <dig> vol of white quartz sand  to the tissue, it was ground to an even finer powder.
the powder was transferred back to a clean  <dig> ml g-tube by suspending the tissue in  <dig>  ml of cold  tca:acetone ) and vigorously mixed and centrifuged, as before.
extraction with 10% tca:acetone was repeated five to seven times, or until no more anthocyanins  could be extracted from the tissue.
this was followed by three washes with chilled  10% tca in water by vigorous mixing and centrifugation, as before, to extract the pectins and remaining anthocyanins from the tissue.
after this, the tissue was washed twice with cold  80% acetone and centrifuged, as before.
protein extraction was performed after drying the tissue pellet to completion in a speed vacuum extractor .
for the preparation of the mesocarp tissue, the same procedure for the exocarp was used with the following modifications.
three g of starting material was used per sample and the first extractions up to the grinding step with white quartz were done in  <dig> ml oakridge tubes.
since some protein can be extracted from the mesocarp via tca:acetone extraction alone  <cit> , a  <dig> min incubation time at -20°c was introduced after the first 100% acetone step and included in the subsequent tca:acetone containing steps to ensure that all of the protein remained precipitated.
in the tca:h2o step, the  <dig> min incubation was done on ice.
since no anthocyanins are present in mesocarp, only two tca:acetone extractions were carried out for the mesocarp tissue.
total protein extraction
two hundred to  <dig> mg of pre-extracted and dried exocarp or mesocarp tissue contained in a  <dig> ml g-tube was extracted by resuspending the pellet in  <dig>  ml cold tris-buffered phenol, ph  <dig> .
then,  <dig>  ml of dense sds buffer  was added.
the mixture was vortexed for  <dig> s and incubated on ice for  <dig> min with intermittent vortexing.
the phenol phase containing the protein as the top phase was separated by centrifugation at  <dig> × g for  <dig> min and transferred into a clean  <dig> ml g-tube.
the remaining sds phase was re-extracted with another  <dig>  ml tris-buffered phenol and incubated for  <dig> min before centrifuging and subsequent transfer and combination of the two phenol phases.
protein was precipitated by adding a minimum of  <dig> vol cold methanol plus  <dig>  m ammonium acetate to the combined phenol phase.
precipitation was carried out at -20°c for  <dig> min or overnight.
after centrifugation at  <dig> g for  <dig> min, the pellet was washed twice with cold methanol containing  <dig>  m ammonium acetate and subsequently with 80% acetone twice.
pellets were next dissolved in 200– <dig> μl fresh buffer containing  <dig> m urea, 2% chaps,  <dig> mm edta, and  <dig> mm hepes, ph  <dig> , to obtain a concentration of approximately  <dig>  μg/μl.
careful sonication on ice was used to dissolve the samples.
protein quantitation was done using a bicinchoninic acid  absorption assay  and read in a victor v plate reader  equipped with a photometric filter of  <dig> nm and  <dig> nm bandwidth.
the quality of each protein sample was checked via sds-page; all samples were devoid of indications of degradation and showed good resolution with low background.
total protein samples were then shipped on dry ice to the university of victoria-genome bc proteomics centre in victoria, bc, for itraq analyses.
using a second bca assay, each protein sample was re-quantified just before aliquoting  <dig> μg of each sample for itraq labeling steps.
experimental design and labeling of peptides with itraq reagents
the experimental design consisted of the four developmental stages described earlier for each of exocarp  <dig>  mesocarp  <dig>  exocarp  <dig>  and mesocarp  <dig>  two biological replicates were employed for each stage and tissue for the  <dig> samples, whereas one  <dig> sample was used for each stage of mesocarp or exocarp.
an additional technical replicate was carried out for exocarp  <dig>  representing separate itraq labeling reactions and analyses starting from the same protein sample.
labeling of peptides with itraq reagents  was performed according to the manufacturer's recommendations as follows.
one hundred μg of each protein sample in a maximum volume of  <dig> μl was precipitated overnight using 100% acetone and dissolved in  <dig> μl of denaturing buffer containing  <dig> μl denaturant and  <dig> μl reducing reagent  as provided in the itraq kit, followed by vortexing and incubation at 60°c for  <dig> h. one μl of cysteine blocking solution  was then added to each sample, followed by incubation at room temperature for  <dig> min.
these protein samples were digested overnight with trypsin  at 37°c.
itraq labeling was carried out by adding itraq reagents  <dig>   <dig>   <dig>  and  <dig> to either the exocarp or the mesocarp samples representing the four developmental stages, green, pink/turning stage, red/fully turned, and purple, respectively.
subsequently, these four samples were mixed by vortexing and further incubated at room temperature for  <dig> h.
the four itraq-labeled peptide samples were pooled together, diluted 1: <dig> with cation exchange sample buffer  containing 25% acetonitrile in  <dig> mm kh2po <dig>  and then adjusted to ph  <dig>  using phosphoric acid.
because of this acidification step, it is important to remove pectins prior to total protein extractions; we found in previous trials that the pectins likely polymerized and precipitated out of solution, converting samples mostly to a gelatinous state unsuitable for further analyses .
the combined peptide mixture was fractionated by strong cation exchange  chromatography on a biocad workstation , using a  <dig>  mm ×  <dig> cm polysulfoethyl aspartamide column .
first, the mixed samples were loaded in buffer a at a flow rate of  <dig>  ml/min.
once completely loaded, the column was washed for  <dig> min with buffer a. peptides were eluted by a linear gradient of  <dig> to  <dig> mm kcl in buffer b .
sixty-nine fractions were collected over the course of  <dig> min at a flow rate of  <dig> ml/min.
of these fractions, only  <dig> fractions containing the eluted labeled peptides as measured by optical density monitoring at  <dig> nm were chosen for analysis on a  <dig> h lc-ms/ms program.
the fractionated samples were reduced to  <dig> μl in a speed-vac  and transferred to autosampler tubes .
liquid chromatography and mass spectrometry
the samples were analyzed for identification and quantitation on a qstar pulsar i hybrid tandem mass spectrometry  system , fitted with a nano-electrospray ionization source  using a  <dig> μm fused silica emitter tip  and interfaced with an integrated lc system consisting of a famos autosampler, switchos ii switching pump, and ultimate micropump .
individual fractions containing peptides were injected onto a  <dig> μm ×  <dig> cm c <dig> pepmap guard column , resolved using a  <dig> μm ×  <dig> mm analytical column , and eluted using an automated binary gradient  from 100% buffer a ,  <dig> % formic acid in h2o) to 40% buffer b  in  <dig> min, then from 40% to 80% buffer b for  <dig> min.
ms time of flight  scans were acquired from m/z  <dig> to  <dig> for one second with up to two precursors selected for ms/ms from m/z  <dig> to  <dig> using information-dependent acquisition at  <dig>  seconds per scan; rolling collision energy was used to promote fragmentation.
custom predicted tryptic peptide database
a schema showing the pipeline for production of the predicted peptide database in support of this subsection is shown in figure  <dig>  all publicly available est data for each vitis species , including those from all v. vinifera  cultivars, were downloaded in august  <dig> as fasta files from the national center for biotechnology information .
these data were parsed on the basis of reported vitis species of origin with the vast majority being from v. vinifera cultivars.
since we were specifically interested in studying the proteome in v. vinifera cv.
cabernet sauvignon  pericarp tissue, an additional, more rigorous approach to the parsing of the cs ests was carried out in order to reduce or eliminate the potential for subsequent assembly of paralogous cs sequences into invalid contigs, thereby striving to strengthen the validity of protein identification in our itraq experiments.
cs ests were obtained from the ncbi genbank database or from an in-house est project  <cit>  and subdivided into the following categories based upon the reported source tissues for the cdnas used for single pass sequencing: whole berry including seed , berry without seed , skin without seed or flesh , seed only , and other tissues  including leaf, flower, tendril, and root.
because the in-house ests were also present in the ncbi genbank database, the corresponding entries in genbank were removed since the genbank entries do not have sequence quality  scores.
the following files containing est data comprised each of the above mentioned groups: vv ; ws ; v. labrusca; v. pseudoreticulata; v. riparia; v. rotundifolia; v. shuttleworthii), cso , css , csp , cse , and csb .
sequences were processed using cross_match  and trim <dig>  in order to remove vector sequences as well as ambiguous nucleotides at the sequence ends.
to successfully perform the above cleanup analyses, phred quality scores were used where available; otherwise, 'place-holder' quality scores were generated for any sequences for which no phred scores were available, as was the case for most of the ests in genbank.
place-holder quality scores were also used later in the cluster assembly process as discussed in more detail, below.
following the cross_match and trim <dig> processing, the sequences were further trimmed using perl scripts designed in-house to eliminate known invalid sequences  and trim polya/t tails, if present in a given sequence.
polya/t stretches were limited to  <dig> bp in order to prevent subsequent chimeric contig assembly based on those repeats.
if polya was followed by a >  <dig> bp stretch of ac, at, gc, or gt repeats, the polya stretch was trimmed to  <dig> bp and all sequence 3' to this was discarded; if polyt was preceded by a >  <dig> bp stretch of ac, at, gc, or gt repeating sequence, the polyt stretch was trimmed to  <dig> bp and all sequence 5' to this was discarded.
if polya started at least two thirds of the est sequence length, it was trimmed to  <dig> bp; if polyt started at less than one third of the est sequence, it was trimmed to  <dig> bp.
any part of a sequence that started or ended with >  <dig> bp of repeats of ac, at, gc, or gt was deleted.
if a sequence started or ended with 'n's , the 'n's were deleted and the corresponding quality scores were also removed.
to better ensure that contig assemblies were based on high quality nucleotide sequence data, percent 'n'  content was determined for each sequence.
if the percentage was >  <dig>  , the flanking  <dig> bp regions where scanned for 'n's and, if present, were trimmed to exclude the 'n's, thereby lowering the total 'n' percentage.
sequences shorter than  <dig> bp were trimmed to the first and last occurrences of an 'n'.
for resulting sequences longer than  <dig> bp, the 'n' percentage was recalculated and, if still >  <dig> %, a record of the sequence was made.
each of these sequences was then compared with other sequences in a combined dataset using blastn to determine its uniqueness.
if a given sequence was already represented in the dataset by another sequence with a lower 'n' content, the sequence in question was eliminated.
the curated sequence datasets were next clustered using pcap software  <cit>  with parameters of 95% overlap identity and  <dig> bp overlap length  <cit> .
pcap was used instead of cap <dig> in order to take advantage of parallelized processing.
parallelization provided the ability to distribute each dataset assembly workload across  <dig> cpus for significantly faster processing time.
the pcap assembly program was modified and recompiled with est_flag set at  <dig> .
the pcap assembly step was followed by a series of post-assembly steps .
we performed two clustering permutations in order to test the effects of database design on peptide identification using our itraq data.
first, we clustered all sequences together to create the "as" database, including ws, vv, and all cs_ files; all sequences were weighted evenly.
second, csb, cse, csp, cso, css, ws, and vv  were clustered separately with higher weighting  placed on cs sequences in the vv build and the original phred scores retained for the in-house cs sequences.
weighting was accomplished by assigning higher quality scores such that when polymorphisms were encountered by pcap in an assembly, preference was given for selection of the cs nucleotide for the resulting contig.
following assemblies, the generated contigs and singletons were merged into one file for each dataset .
any sequences longer than  <dig> bp were suspected to be chimeric, so they were parsed to a separate file, translated in all  <dig> frames, and peptides with a minimal size of  <dig> amino acids before a predicted stop codon were submitted to a blastx search against the nr database.
the resulting multiple peptides predicted within long contigs were coded with "lc", as well as with "f" for the translational frame, with the frame number  and the peptide number designated from among the multiple peptides .
a blastx analysis was next performed on each contig and singleton sequence against the nr database in order to identify the best frame for subsequent in silico translation.
the frame identified via blastx analysis was then used to generate the predicted orf  for a given contig or singleton.
in order to further curate predicted orfs, each was subjected to in silico cleavage at any 'unknown' amino acid  or stop codon and then compared to a similarly generated list of peptides from the corresponding best scoring protein sequence identified in the blastx search.
the 'best peptide' was then identified in the translation frame as the peptide with an exact match to the blastx peptide.
if no such peptide could be identified, the longest peptide generated by in silico cleavage of the sequence at each occurrence of an unknown amino acid and/or stop codon was used.
all sequences which resulted in "no hit found" in the blastx results were subsequently translated in all six frames and appended to the end of the 'best peptide' file.
in all cases where a six-frame translation was applied, the resulting peptides  were cleaved in silico at every unknown amino acid and/or stop codon and only those sequences  <dig> amino acids or longer were kept.
the resulting list of 'best peptides' for each of the sequences was then subjected to blastp analysis using the uniprotkb database in order to determine the sequence identity.
the five highest blastp hits for each query sequence were aligned using an in-house perl script to identify putative n- and/or c-termini.
if no consensus site could determined for an n- or c-terminus via alignment with similar sequences, then these sequences were trimmed at the tryptic digestion site nearest to the ends of the predicted orf to eliminate potentially truncated predicted tryptic peptides from the database.
the parameters programmed into the scripts included: 1) blastp e-values, where e = h indicates a less significant hit  and e = l indicates a stronger hit , 2) the difference in length of the vitis query sequence versus each of the top five subject sequences, and 3) the length of the exact match of amino acids to the top hit.
these parameters improved automation of accurate predictions of methionine  sites and identification of likely full-length amino acid sequences without requiring manual inspection of the blastp results.
sequences identified with a predicted methionine at the n-terminus were coded with ''.
determinations of c termini were done based on a small range cutoff  of amino acids between the stop codon in each top hit and the predicted stop codon in each corresponding vitis orf; if the difference was greater than two amino acids and deemed unclear, the c-terminal end of the predicted orf was trimmed at the nearest upstream tryptic cleavage site.
the detailed process, above, was applied to each data set individually .
in preparation for the merger of the datasets, the cs tissue-specific sequences were analyzed for uniqueness based on comparing each sequence to every other sequence and discarding all shorter sequences for each exact match.
once all cs duplicate sequences were removed, this dataset was merged with the remaining two sets, vv and ws.
this final set consisting of all of the sequences was then subjected to yet another uniqueness test where each sequence was compared to every other sequence but this time cs sequences where intentionally not removed, even if an exact duplicate existed in either the vv or the ws set and was of greater length.
this allowed for preferential retention of the cs sequences in order to keep information about the tissue of origin of a detected protein.
out of a total of  <dig> sequences submitted to the uniqueness test,  <dig> were identified with cs duplicates present due to the preferential retention of those sequences.
from the resulting sequences, only those that started with a predicted methionine were then submitted for signalp analysis  and processing which allowed for the identification and ultimate trimming of signal peptides giving rise to predicted mature protein sequences.
those trimmed sequences that had a predicted cleavable targeting signal were coded ''.
after removal of the predicted signal peptides, a final uniqueness test was performed.
analysis of ms/ms data
itraq ms/ms data were analyzed using proteinpilot software v.  <dig> .
<dig>  for both tryptic peptide identification and quantification.
the peptides and corresponding relative abundances were obtained in proteinpilot using a confidence cutoff  of >  <dig>  .
database searching for each sample was done on predicted tryptic peptide sequence data using either the msdb database  or in-house databases .
annotations and annotated protein names indicated in proteinpilot output files were coded to indicate several parameters specific to the orf identified as well as the est or contig from which the orf sequence was predicted.
itraq data representing the four ripening initiation stages in each of the three exocarp samples  were combined into a single tab delimited file.
likewise, itraq data representing each of the three mesocarp samples  were combined into a second tab delimited file.
duplicate entries among exocarp or mesocarp files were identified using an in-house script in the r environment with 'custom orf id' as the search string.
then, ratiometric data at each of the three comparisons using 'green' as the reference stage  were averaged prior to export for cluster analyses.
entries with the same name but different template cdna sources were not averaged since these may represent isoforms from different source tissues and/or cultivars.
we chose to cluster all proteins detected in the exocarp or mesocarp in order to capture all information on expression patterns detected, without restricting our analyses to only those proteins that were replicated amongst the individual exocarp or mesocarp files.
k-means clustering into four partitions was carried out on ratiometric data for the exocarp and mesocarp files separately using multiexperiment viewer  software .
we used a  <dig> -fold threshold for biological significance which was validated by consistencies between trends in protein expression presented here as increasing or decreasing with corresponding patterns of gene expression identified in previous publications .
RESULTS
protein detection and annotation
in order to strengthen the ability of ms/ms spectra annotation software such as proteinpilot to accurately identify peptide sequences in complex total protein samples, we performed weighted clustering on a large est collection generated from vitis spp.
to create the most cultivar-specific peptide map database possible for cabernet sauvignon.
in addition, we created perl scripts in order to find the n and c termini in translated orfs or otherwise to carry out trimming at known tryptic digestion sites.
this was done to help decrease the number of incorrectly annotated, non-tryptic  n- and/or c-terminal peptides incorporated into protein abundance quantitations.
to further increase the number of identifiable peptides, our in-house scripts searched all of the predicted orfs designated as containing a likely start methionine using signalp to identify and remove putative n-terminal signal peptides, resulting in a database containing mature protein sequence data for the predicted full length proteins.
to assess whether the custom tryptic peptide database improved protein discovery, we examined the number of high confidence proteins and peptides identified by proteinpilot using this database versus a common approach to annotation, a search of the msdb.
we performed this analysis on two itraq data sets, mesocarp 2005- <dig> and exocarp 2005- <dig>  derived from four stages of cabernet sauvignon berries at ripening initiation.
at a confidence level of 95% ,  <dig> proteins were identified in the mesocarp using the custom database, whereas only  <dig> proteins were identified in a search of the same data set using the msdb .
at a confidence level of 95%,  <dig> proteins were identified in the exocarp using the custom database, whereas  <dig> proteins were identified in a search of the same data set using the msdb .
these results indicate that in these two itraq data sets, the use of the custom tryptic peptide database improved high confidence protein discovery by  <dig> % and  <dig> %.
the number of high confidence peptides detected using the custom database was  <dig> -fold and  <dig> -fold higher in mesocarp and exocarp, respectively, in comparison to searches of the same itraq data sets using the msdb .
the greater difference in the number of high confidence peptides versus proteins detected using the custom database in comparison to the msdb indicates that the most important impact of implementing the custom database was that it afforded the identification of more high confidence peptides per protein than could be achieved using the msdb.
arelease  <dig>  hammersmith campus of imperial college london
best clustering weighted higher scores for ests derived from cabernet sauvignon in order to avoid the potential loss of snps specific to this cultivar during consensus sequence determination by the pcap software
cpredicted tryptic peptides that were determined to likely be truncated were trimmed and removed from the custom peptide database
the effects of weighting and trimming on high confidence peptide and protein detection were analyzed using the  <dig> exocarp itraq data set.
we performed a second build of all vitis spp.
ests using pcap.
for the purposes of this comparative analysis, all ests in this second build were equally weighted for consensus sequence determination using an arbitrary phred score; no higher weighting of cabernet sauvignon-derived ests was used.
this unweighted vitis est database was then used to generate a second predicted tryptic peptide database, including trimming of end-truncated peptides, as before.
the number of high confidence peptides and proteins detected were similar using the orf databases created from the weighted versus the unweighted est databases , suggesting that the effect of simple nucleotide polymorphisms  among the vitis spp.
ests is ultimately negligible at the level of peptide identification via ms/ms in grapevine.
the effectiveness of trimming predicted end-truncated peptides, on the other hand, was evident in the increased number of peptides identified  by proteinpilot in searches of the trimmed  orf database versus a third, untrimmed custom vitis orf database that we created for this comparison by omitting the trimming step shown in figure 1; both of these databases were weighted for cabernet sauvignon ests in the original pcap builds.
these findings indicate that automated searching for the best predicted methionine start site and removing truncated n- and/or c-terminal tryptic peptide sequences in predicted orfs based upon est data increases the ability of ms/ms spectra annotation software to identify peptide sequences with high statistical confidence.
technical and biological replication of expression ratios
proteins identified in itraq data sets from two technical replicates were analyzed.
two  <dig> μg aliquots from the same total protein sample, exocarp  <dig>  were separately labeled with itraq reagents, subjected to nanolc-ms/ms, and the custom predicted tryptic peptide database was searched with the resulting spectra using proteinpilot.
only proteins that were detected in each of the four ripening initiation stages were retained for the comparative analysis.
in-house custom scripts written in r language  <cit>  were employed to detect protein entries that were duplicated in the two proteinpilot output files representing the technical replicates, using 'custom orf id' as the search string.
out of a total of  <dig> proteins detected in both itraq data sets,  <dig> or 29% of these proteins were detected in both files , whereas the remainder were unique to either technical replicate with  <dig> proteins in one file and  <dig> proteins in the other.
we were interested in estimating consistency in trends in expression data along ripening initiation for each replicated protein entry that we detected .
we determined that strict estimation of replicability via pearson's correlation coefficient or concordance correlation coefficient analyses of expression ratios did not provide relevant quantitative information , since even small deviations in expression ratios can result in poor correlative data.
we considered that identifying consistent trends in protein accumulation above an arbitrary cutoff such as  <dig> -fold was most relevant for the purposes of discovering protein candidates for hypothesis formulation related to ripening control.
as an alternate approach to correlation estimation, a replication scoring system was devised, whereby expression ratios identified by manual inspection as differing by less than log <dig> for each duplicate protein entry for at least two of the three ratios received a score of  <dig>  expression ratios differing by between log <dig> and log <dig> for at least two of the three ratios received a score of  <dig>  and expression ratios differing by log <dig> or greater for at least two of the three ratios received a score of  <dig>  the mean replication score for the  <dig> duplicated protein entries in the technical replicates was determined to be  <dig>  ±  <dig>  ; 17% of the protein entries received a score of  <dig>  these data indicate that if a  <dig> -fold change or greater is considered as biologically relevant, over 80% of duplicate protein entries but only approximately one quarter of all proteins detected in two technical replicates would be expected to be identified as significantly changing along ripening initiation with similar quantitative trends, based on the technical replicates analyzed here.
proteins identified in itraq analyses of two biological replicates, exocarp  <dig> and exocarp 2005- <dig> were next analyzed in an identical manner to the two technical replicates.
out of a total of  <dig> proteins detected in the two itraq data sets,  <dig> or 33% of these proteins were detected in both files , whereas the remainder were unique to either biological replicate with  <dig> proteins in one file and  <dig> proteins in the other.
the same replication scoring system was implemented for proteins detected in both itraq data sets  as was done for the technically replicated proteins , above.
the mean replication score for the  <dig> proteins common to both biological replicates was determined to be  <dig>  ±  <dig>  ; 12% of protein entries received a score of  <dig>
to further explore the biological significance of replicability analyses, proteins common to either the technical replicates  or the biological replicates  were compared and it was determined that  <dig> protein entries were common to the technical replicates and the biological replicates, i.e.
were detected in all three data sets analyzed for replication ; we considered this approach as valid since all data analyzed here were derived from the exocarp.
the technical and biological replication results indicate that the 'shotgun' lc-ms/ms approach employed here with itraq-labeled exocarp total proteins is capable of repeatedly identifying approximately one third of the total number of proteins, whether they are technical or biological replicates.
expression profiling of proteins along grape berry ripening initiation
four clusters were generated for the exocarp in which trends in protein accumulation were increasing strongly , increasing gradually , not changing significantly , or decreasing , respectively, from the green through the fully turned ripening initiation stages .
additional file  <dig> lists proteins by cluster number along with corresponding log2-transformed ratiometric data for each protein entry.
ratiometric data were calculated relative to the green stage for each of the pink/turning, red, and purple/fully turned stages.
several protein isoforms significantly increasing along ripening initiation were identified with annotated functions in anthocyanin flavonoid biosynthesis and storage , defense ), and cell expansion .
conversely, several components of the photosynthetic machinery were identified as decreasing greater than  <dig> -fold along ripening initiation, which is consistent with a reduction in photosynthesis at this period of berry development.
we mined the exocarp data  for proteins that were increasing in abundance relative to the green stage and annotated as enzyme or transporter components of pathways leading to hypothesized regulators of ripening initiation in grapes, aba, glucose, and br.
a putative lytb  protein  increased  <dig> -fold in abundance and is responsible for the last step of the plastidic pathway to isopentenyl diphosphate , leading in part to the production of the plant hormones, aba and gibberellic acid.
other proteins of this mep pathway were also detected but only as slightly increasing in abundance along ripening initiation.
an isopentenyl-diphosphate δ-isomerase i protein was detected as increasing 2-fold along ripening initiation ; this enzyme controls a major early step in isoprenoid biosynthesis and is likely localized to the chloroplast.
we identified one component specific to the aba biosynthetic pathway, a protein similar to violaxanthin de-epoxidase from tea , which was stably expressed along ripening initiation.
cytosolic idp, potentially also formed in the plastid and exported to the cytoplasm, is incorporated in the biosynthesis of brs.
a putative grapevine ortholog to a br biosynthetic protein from pea, pslkb, , was identified as increasing  <dig> -fold, peaking at the third  stage tested.
a grapevine hexose transporter  was identified as increasing  <dig> -fold, which is consistent with hexose accumulation during ripening initiation.
few proteins annotated with signal transduction functions were detected in any of the four clusters that could hypothetically be involved in aba, glucose, and/or br signaling.
abscisic stress response protein  was previously implicated in cross-talk between aba and glucose signaling  <cit> ; expression data were variable, with some isoforms increasing up to 4-fold along ripening initiation and others showing stable expression.
a putative pirin protein of unknown function increased  <dig> -fold along ripening initiation.
a putative ortholog to the malus spp.
ttg <dig> wd <dig> repeat protein , which has previously been implicated in arabidopsis in the regulation of anthocyanin biosynthesis  <cit> , was detected as increasing, indicating that vvttg <dig> could play a similar role in the grape exocarp; a caveat is that the confidence level for vvttg <dig> unique peptide detection was less than 95%.
curiously, we did not identify peptides representing vvmyba <dig>  which is highly expressed at the transcriptional level during ripening initiation and encodes a transcription factor previously demonstrated to positively regulate anthocyanin biosynthesis genes in grape  <cit> .
four clusters were generated for the mesocarp in which trends in protein accumulation were increasing strongly , increasing gradually , not changing significantly , or decreasing , respectively, from the green through the fully turned ripening initiation stages .
additional file  <dig> lists proteins by cluster number along with corresponding log2-transformed ratiometric data for each protein entry, relative to the green stage.
similar to the exocarp data, we identified > 2-fold accumulation of several protein isoforms annotated with functions in cell enlargement , fruit softening  and defense  concomitant with a significant reduction in proteins inhibitory to fruit ripening .
we mined the mesocarp data  for proteins annotated as enzyme or transporter components of pathways leading to aba, glucose, and brassinosteroid accumulation.
we detected vvnced <dig>  increasing over 2-fold along ripening initiation, which represents expression of a key committed step specifically in aba biosynthesis in the plastid  <cit> .
a protein annotated as aba glucosyltransferase increased  <dig> -fold along ripening initiation in the mesocarp.
as in the exocarp, an isopentenyl-diphosphate δ-isomerase i was detected as increasing 2-fold along ripening initiation , plus, a putative cytosolic isoform was detected as stably expressed   <cit> .
farnesyl diphosphate synthase, a cytoplasmic enzyme leading to the br biosynthetic pathway via squalene biosynthesis, was detected in the mesocarp as increasing  <dig> -fold along ripening initiation.
a putative grapevine ortholog to the br biosynthetic protein in pea, pslkb, , and in cotton, ghdwarf <dig> , was detected as increasing  <dig> -fold.
the hexose transporter, vvht <dig>  was detected as increasing 2-fold, similar to its expression pattern that was detected in the exocarp.
we mined the mesocarp data further to identify proteins annotated with signal transduction functions.
no plasma membrane receptor candidates were identified as increasing but several were identified as not changing or decreasing along ripening initiation, most notably a receptor-like kinase similar to perk1-like protein from rice  that decreased  <dig> -fold at the onset of pigment accumulation to  <dig> -fold at the purple stage; this perk1-like protein was also identified as significantly decreasing in the exocarp .
the same pirin protein identified in the exocarp was similarly detected as increasing in the mesocarp along ripening initiation.
inconsistent itraq data were obtained for vvasr in the mesocarp with isoforms showing an increase, stable expression, or strong down-regulation.
a custom script written in the r environment was used to search for proteins common to the combined exocarp and mesocarp files  using 'custom orf id' as the search string.
protein isoforms detected in common between exocarp and mesocarp  are shown in additional file  <dig>  protein isoforms detected only in exocarp  are shown in additional file  <dig>  protein isoforms detected only in mesocarp  are shown in additional file  <dig>  isoforms of anthocyanin flavonoid biosynthetic proteins were only detected in the exocarp, which is consistent with induction timing and tissue localization of these pigments during ripening initiation and demonstrates the efficacy of our separation of exocarp and mesocarp for proteomic analyses.
additional file  <dig> shows the source tissues and number of unique high confidence peptides per protein for those proteins indicated here in the results section and in figure  <dig>
discussion
the itraq data obtained with exocarp and mesocarp total proteins confirmed previous ratiometric transcript abundance data for key components of aba and br biosynthesis, as well as the influx and accumulation of sugars during ripening initiation .
we confirmed for the first time the strong accumulation of several components of anthocyanin biosynthesis at the protein level, including a putative anthocyanin o-methyltransferase.
the detection of vvnced <dig> accumulation confirmed previous real-time rt-pcr data  <cit>  and further supports a role specifically for this nced family member in aba biosynthesis in the mesocarp during berry ripening initiation.
based on the sensitivity limitations to the shotgun proteomic technique employed here, we cannot conclude that vvnced <dig> accumulation, activity, and, consequently, aba biosynthesis are localized to the mesocarp; vvnced <dig> may be expressed in the exocarp but we did not detect it.
the previous detection of ζ-carotene desaturase transcripts in the exocarp  <cit>  argues against a model in which aba is synthesized in the mesocarp and transported to the exocarp whereupon it activates, in part, anthocyanin biosynthesis.
if the primary site of aba production in the developing fruit is the seed, however, and aba transported from the seed to the surrounding pericarp is an early signaling event initiating ripening, as we have previously hypothesized  <cit> , it is reasonable that the first tissue in the pericarp in which aba is synthesized is the mesocarp, given its closer proximity to the seed than the exocarp.
it remains to be determined whether a gradient in aba biosynthesis during early ripening stages proceeds from the seed through the mesocarp and, finally, in the exocarp.
the presence and possible moderate accumulation of a putative aba glucosyltransferase indicates that an attenuating mechanism for the aba signal may operate during ripening initiation to control aba homeostasis  <cit> .
the accumulation of the pirin protein is also intriguing, given this protein's previously demonstrated interaction with g-protein alpha in arabidopsis seeds and a potential role in modifying aba action via negative feedback control  <cit> .
we used a relatively new quantitative ms/ms-based approach, itraq, to advance our understanding of differential protein expression underlying non-climacteric ripening initiation in grape berries.
the itraq approach offered several advantages over 2dge methods for protein discovery, including increased detection sensitivity based on our findings reported here in comparison to previous reports on grape berry proteomics.
by using strong cation exchange and reverse phase column microcapillary chromatography coupled with nanospray ms/ms detection with total protein extracts from grape berries, we were able to resolve three-fold or more proteins per sample than would be expected using 2dge  <cit> .
one current limitation to an ms-based proteomic approach with grapevine is that there are no finished genome sequence data for grapevine, although two projects  <cit>  are undertaking assembly and annotation from which a high quality orfeome database can eventually be derived.
although there are over  <dig>  vitis spp.
ests deposited in genbank, v. vinifera is a highly heterozygous species  <cit> ; therefore, we considered that it might be important to weight builds of these ests via manipulation of phred scores so as to favor sequence data corresponding to our cultivar of interest, cabernet sauvignon, when snps were encountered by pcap, where applicable for a given contig assembly.
although we determined that weighting ests to the genotype of interest for est assembly provided no clear advantages in this study, we conclude that producing a tryptic peptide database targeted to vitis sequences and including removal of predicted truncated peptides improved protein detection and annotation.
furthermore, our findings indicate that a tryptic peptide database based on finished pinot noir whole genome sequence data will be valid to implement for proteome studies with any v.vinifera cultivar or vitis species, with the exception of cases where deletions have occurred in the pinot noir homozygous line  <cit> .
we chose to not include genome sequence data available for v. vinifera cv.
pinot noir  <cit>  due to significant gaps in current assemblies and the potential for inaccurate automated gene predictions.
until grapevine genome sequence assembly and annotation are finished, we propose that the predicted orf database presented here will be of value to the grapevine community in two significant ways.
while gaps exist in the genome sequence assemblies, the protein database presented here may provide information for 'missing' proteins either not yet predicted from the pinot noir genome sequence data and/or other vitis spp.
not represented in the pinot noir genome sequence data, e.g.
due to chromosomal deletions.
further, protein predictions based upon expressed sequences  represent 'real' proteins, so our database could potentially be used to validate orf predictions based on whole genome sequence data alone.
while itraq labeling coupled with nanolc-ms/ms proved overall to be an advance over 2dge in sensitivity and quantitation, consistent detection of predicted protein sequences between technical or biological replicates from grape exocarp was limited.
consistency in trends in ratiometric data along ripening initiation for those proteins that were detected in replicate exocarp samples was further limited.
inconsistent ratiometric data for some proteins detected in both biological replicates may represent differences in expression due, for example, to variability in seasonal growing conditions and not technical variability.
nonetheless, limited replicable detection of proteins among biological samples for technical reasons is commonly encountered with liquid chromatography and ms-based proteomics  <cit>  and could arise from variation in preparation of total proteins and/or itraq-labeled peptides from sample to sample.
underlying the apparent variation arising from sample preparations are constraints on detection imposed by the mass spectrometer, both with respect to the dynamic range of the instrument and the nature of selections of peptides in the first ms by the ms software for export to the collision cell prior to amino acid detection via the second ms. furthermore, our finding that two thirds of the replicated exocarp proteins were detected both in technical and biological replicates may have reflected the higher probability of detecting these mainly abundant, housekeeping-type proteins in any given total protein sample using a shotgun approach.
over-sampling of abundant peptides in digested total protein samples is a limitation of the shotgun approach to quantitative proteomics  <cit>  and likely precluded our ability to discover more proteins annotated with signal transduction functions that could regulate aba, br, and hexose responses during ripening initiation.
in order to increase detection sensitivity of low abundance regulatory proteins using shotgun proteomics techniques, it will likely be helpful to isolate membrane and nuclear proteins separately from cytoplasmic proteins prior to digestions.
affinity chromatography of berry protein extracts using antibodies directed against abundant proteins such as thaumatin detected in both exocarp and mesocarp may also improve detection sensitivity for low abundance proteins by selectively removing these proteins prior to itraq labeling steps.
similarly, advances in detection sensitivity in devices such as the fourier transform ion cyclotron resonance ms   <cit>  should allow us to delve deeper into the grape berry proteome in order to better understand the molecular control of non-climacteric ripening in this species.
nonetheless, the searchable quantitative data presented here on over  <dig> proteins detected in two separate berry tissues along four ripening initiation stages can serve as a framework public resource contributing towards our understanding of the dynamic grapevine proteome.
CONCLUSIONS
we determined that a predicted peptide database can be derived from grapevine est data using advanced clustering and trimming approaches and successfully implemented for quantitative proteome profiling.
we demonstrated in grapevine that by implementing a predicted peptide database targeted to an organism of interest, a significant increase in the number of high confidence peptides identified and annotated from ms/ms data was gained in comparison to more commonly used msdb searches.
furthermore, by using a shotgun quantitative proteomics approach in combination with a targeted predicted peptide database, we showed that greater numbers of high confidence peptides were detected than would be expected using 2d gel electrophoresis techniques.
we verified for the first time at the protein level, quantitative expression patterns for components of isoprenoid and flavonoid metabolism important for grape compositional chemistry, as well as identified new signal transduction candidate proteins associated with non-climacteric ripening initiation in grape berries.
abbreviations
aba: abscisic acid; blast: basic local alignment and search tool; br: brassinosteroids; est: expressed sequence tag; itraq: isobaric tagging for relative and absolute quantitation; orf: open reading frame; pcap: parallel contig assembly program
authors' contributions
jl designed the experiments, carried out protein extractions and data analyses, assisted in designing the scripts for the predicted peptide database, and drafted the manuscript.
ml designed the scripts and carried out programming to produce the predicted peptide databases presented here.
ds carried out mass spectrometry and database searches using the itraq ms/ms data.
stl performed primary data analyses, assisted in drafting the manuscript, and was responsible for overseeing the project.
supplementary material
additional file 1
duplicated proteins and associated replication scores for exocarp  <dig> technically replicated itraq experiments.
ratiometric data are given in log <dig> scale.
click here for file
additional file 2
duplicated proteins and associated replication scores for exocarp  <dig> and exocarp 2005- <dig> biologically replicated itraq experiments.
ratiometric data are given in log <dig> scale.
click here for file
additional file 3
proteins identified in all of the three itraq data sets analyzed  and exocarp 2005-1) and tested for replication in expression trends along ripening initiation.
click here for file
additional file 4
exocarp proteins listed for each k-means cluster along with corresponding log2-transformed ratiometric data for each protein entry.
exocarp proteins from the  <dig>  2005- <dig>  and 2005- <dig> datasets were combined into a single file.
duplicate entries were identified using an in-house script in the r environment with 'custom orf id' as the search string.
then, ratiometric data at each of the three comparisons using 'green' as the reference stage were averaged prior to export for k-means cluster analyses.
click here for file
additional file 5
mesocarp proteins listed for each k-means cluster along with corresponding log2-transformed ratiometric data for each protein entry.
mesocarp proteins from the  <dig>  2005- <dig>  and 2005- <dig> datasets were combined into a single file.
duplicate entries were identified using an in-house script in the r environment with 'custom orf id' as the search string.
then, ratiometric data at each of the three comparisons using 'green' as the reference stage were averaged prior to export for k-means cluster analyses.
click here for file
additional file 6
protein isoforms detected in common between the exocarp and mesocarp .
a custom script written in the r environment was used to search for proteins common to the combined exocarp and mesocarp files using 'custom orf id' as the search string.
click here for file
additional file 7
protein isoforms detected only in exocarp samples.
proteins excluded from additional file  <dig> that were detected in the combined exocarp file  are shown.
click here for file
additional file 8
protein isoforms detected only in mesocarp samples.
proteins excluded from additional file  <dig> that were detected in the combined mesocarp file  are shown.
click here for file
additional file 9
protein source tissues and number of unique high confidence peptides per protein.
proteins presented in the results and discussion sections as well as in figure  <dig> are shown.
the source tissue for each protein hit shown in column a corresponds to a proteinpilot output file, which can be downloaded at  the number of unique > 95% confidence peptides determined by the proteinpilot software for each protein is shown in column e. ratiometric data shown in columns g, h, and i are relative to the green stage  and are not log2-transformed.
click here for file
additional file 10
reference key for coding terms used in column a  and column b  in additional files  <dig> through  <dig>  the key is provided as a printable rapid reference for data mining in additional files  <dig> through  <dig>  as well as in the protein database available online, as shown in the acknowledgments.
click here for file
acknowledgements
the authors wish to thank monica elliot for assistance with protein preparations and itraq labeling.
we thank dr. fred y. peng for assistance with itraq data processing and cluster analyses.
the authors gratefully acknowledge funding for this research from genome canada and project management support from genome british columbia as part of the genome canada-genoma españa collaborative research and development initiative.
the authors also wish to thank the province of british columbia and the bc wine industry for their continued support of research activities in the university of british columbia wine research centre.
resources and raw data developed through the research presented here can be downloaded freely at:  these are: 1) scripts used for the customized grapevine protein database creation, 2) the protein database, and 3) the proteinpilot output files for exocarp  <dig>  exocarp 2005- <dig>  exocarp 2005- <dig>  mesocarp  <dig>  mesocarp 2005- <dig>  and mesocarp 2005- <dig>  proteinpilot software can be obtained by contacting applied biosystems.
the scripts for protein database creation will need to be modified prior to implementation with a species different to v. vinifera; the corresponding author may be contacted for additional information.

 the field of protein discovery through mass spectrometry continues to grow rapidly but the number of species for which finished whole genome sequence data are available is currently not keeping pace. for a large number of laboratories worldwide studying proteomes in'non-mainstream' organisms, annotations of tandem mass spectra data must rely on open reading frame predictions from expressed sequence tag data from their species of interest or a phylogenetically close relative. ests generated from single pass sequencing reactions are frequently not full length and the reading frames are unknown. proteolytic peptide sequence databases derived from multiple, truncated predicted orfs per each of thousands of ests can hamper the ability of search engines such as mascot and algorithms such as paragon in proteinpilot software to make statistically robust protein identifications from ms/ms spectrum data. protein identifications from ms/ms spectra may be further complicated when the est data that are used to build a peptide sequence database are created based on one genotype for a given species. we report here on the development of scripts for the generation of a predicted tryptic peptide sequence database based on est data in grapevine. our computational approach accounts for multiple open reading frames, truncated predicted orfs, and the presence of n-terminal signal peptides, and may be useful for ms/ms-based protein discovery in any species for which est data are available. quantitative protein expression profiling analyses in plants have increasingly implemented stable isotopic labeling as an advance or complement to two dimensional gel electrophoresis methods. isotope coded affinity tagging reagents are used to covalently label cysteine residues with heavy or light hydrogen or carbon in two complex peptide samples, for example, wild type versus mutant genotypes. the icat chemistry is used to purify labeled peptides via affinity chromatography and then samples are mixed and subjected to lc-ms/ms. one of the first reports on an icat application in plants was in wheat where relative expression in monosomic deletion mutants was used to begin to clarify the influence of ancestral genomes on differential seed protein expression for breeding applications. the icat technique is limited, however, by the tagging of cysteine residues only, as well as the need for affinity purification of labeled peptides; invariably, information is lost through these steps. an improvement to the icat technique involves the labeling of amine groups using a set of four or more isobaric tags. the advantages of this technique, isobaric tagging for relative and absolute quantitation, are that most peptides are labeled, no affinity purification step is required, and the isobaric nature of the tags allows co-elution of identical peptides that are differentially tagged, thereby enriching detection sensitivity and accuracy in comparison to icat. few reports of itraq implementation in plant proteome studies have been reported but pioneering work in this field has been successful, for example, in further defining the organellar proteome in arabidopsis, characterizing pathogen defense mechanisms in arabidopsis, and clarifying micronutrient stress responses in barley. we were interested in characterizing ripening initiation in grape berries at the level of differential protein expression in order to better define the molecular control of this important process for grape growers and wine makers. grape berry ripening is non-climacteric and ethylene does not act as a major signal initiating this process, as it does in climacteric species such as tomato. abscisic acid, hexoses, and brassinosteroids have previously been implicated in non-climacteric ripening regulation but how these and potentially other signaling pathways interact to effect major changes in berry biochemistry at ripening initiation is poorly understood. the tissues in the grape berry consist of the seeds, the mesocarp, and the exocarp; the pericarp refers to the mesocarp and exocarp, collectively. primary and secondary compounds important for grape and wine products begin to accumulate in the exocarp and/or the mesocarp at ripening initiation, so we considered that it was important to evaluate changes in the berry proteome separately in these tissues. to date, a limited number of reports on proteome profiling in grapevine and grape berries have been published in which 2dge was employed. we considered that the itraq technique could be useful in surmounting some technical limitations encountered with 2dge and allow us to detect a greater number of proteins per sample. in this report, we demonstrate the application of our computational approach to tryptic peptide sequence database development from a large collection of grapevine est data and validate its usefulness by showing improved detection and annotations of ms/ms data derived from grape exocarp and mesocarp total protein extracts. we further provide new quantitative information on differential protein expression during ripening initiation in grape berries. this is the first report in which itraq has been used to study differential protein expression in any fruit. methods plant material grape clusters were sampled from v. vinifera cv. cabernet sauvignon clone grafted on rootstock 101- in a commercial vineyard near osoyoos, british columbia, in the and seasons. sampling dates during each season were focused on the developmental stages undergoing ripening initiation. clusters were sampled on a single date in august 12th, which was the timing of approximately 50% ripening initiation based on a turning pink color phenotype. for the season, the ripening initiation stage was sampled over a longer period, since in this growing season, ripening advanced slowly due to lower atmospheric temperatures. five clusters from five different vines were sampled in each season and snap frozen directly in liquid nitrogen in the vineyard and then transported on dry ice to ubc vancouver where they were stored at -80°c. individual grapes from each of the and clusters were developmentally staged based on a visual pigmentation assessment and were segregated for each season into green, pink/turning, fully turned red, and fully turned purple phenotypic classes. for the samples, green grapes were only taken from clusters collected on august 10th, since for this date and august 12th there was no visible change in color present in any of the grape clusters. thirty grapes of similar sizes per pigmentation class per year were segregated for experimentation. prior to total protein extraction, individual grapes were partially thawed in gloved hands and then, using a forceps, the exocarp tissue was carefully peeled away from the mesocarp and placed immediately into liquid nitrogen. seeds were then carefully removed while keeping the remaining mesocarp tissue frozen in liquid nitrogen. exocarp and mesocarp samples were ground to a powder under liquid nitrogen and then used for total protein extractions. tissue preparation for protein extraction preparation of exocarp tissue samples for protein extraction was performed according to a previously described protocol for olive leaf with some modifications described here. the procedure was carried out on ice and centrifugations were performed at 4°c. throughout the procedure, each wash was done by complete resuspending of the tissue pellet. four hundred mg of powdered exocarp tissue was placed in a ml g-tube. the tissue was suspended in ml of a cold ethyl acetate:ethanol ) solution by vortexing for s; the ethyl acetate:ethanol extraction was previously found to be useful for removing pectins as well as pigments such as chlorophylls. following centrifugation for min at × g, the supernatant was removed and the ethyl acetate:ethanol extraction and centrifugation steps were repeated on the remaining tissue. the sample was next extracted twice with cold 100% acetone by vortexing and centrifuging, as before. subsequently, the tissue with added acetone was transferred from the g-tube to a mortar using a ml pipette with the tip end excised to increase diameter and then the acetone was evaporated from the tissue at room temperature. after the addition of 1/ vol of white quartz sand to the tissue, it was ground to an even finer powder. the powder was transferred back to a clean ml g-tube by suspending the tissue in ml of cold tca:acetone ) and vigorously mixed and centrifuged, as before. extraction with 10% tca:acetone was repeated five to seven times, or until no more anthocyanins could be extracted from the tissue. this was followed by three washes with chilled 10% tca in water by vigorous mixing and centrifugation, as before, to extract the pectins and remaining anthocyanins from the tissue. after this, the tissue was washed twice with cold 80% acetone and centrifuged, as before. protein extraction was performed after drying the tissue pellet to completion in a speed vacuum extractor. for the preparation of the mesocarp tissue, the same procedure for the exocarp was used with the following modifications. three g of starting material was used per sample and the first extractions up to the grinding step with white quartz were done in ml oakridge tubes. since some protein can be extracted from the mesocarp via tca:acetone extraction alone, a min incubation time at -20°c was introduced after the first 100% acetone step and included in the subsequent tca:acetone containing steps to ensure that all of the protein remained precipitated. in the tca:h2o step, the min incubation was done on ice. since no anthocyanins are present in mesocarp, only two tca:acetone extractions were carried out for the mesocarp tissue. total protein extraction two hundred to mg of pre-extracted and dried exocarp or mesocarp tissue contained in a ml g-tube was extracted by resuspending the pellet in ml cold tris-buffered phenol, ph. then, ml of dense sds buffer was added. the mixture was vortexed for s and incubated on ice for min with intermittent vortexing. the phenol phase containing the protein as the top phase was separated by centrifugation at × g for min and transferred into a clean ml g-tube. the remaining sds phase was re-extracted with another ml tris-buffered phenol and incubated for min before centrifuging and subsequent transfer and combination of the two phenol phases. protein was precipitated by adding a minimum of vol cold methanol plus m ammonium acetate to the combined phenol phase. precipitation was carried out at -20°c for min or overnight. after centrifugation at g for min, the pellet was washed twice with cold methanol containing m ammonium acetate and subsequently with 80% acetone twice. pellets were next dissolved in 200– μl fresh buffer containing m urea, 2% chaps, mm edta, and mm hepes, ph, to obtain a concentration of approximately μg/μl. careful sonication on ice was used to dissolve the samples. protein quantitation was done using a bicinchoninic acid absorption assay and read in a victor v plate reader equipped with a photometric filter of nm and nm bandwidth. the quality of each protein sample was checked via sds-page; all samples were devoid of indications of degradation and showed good resolution with low background. total protein samples were then shipped on dry ice to the university of victoria-genome bc proteomics centre in victoria, bc, for itraq analyses. using a second bca assay, each protein sample was re-quantified just before aliquoting μg of each sample for itraq labeling steps. experimental design and labeling of peptides with itraq reagents the experimental design consisted of the four developmental stages described earlier for each of exocarp mesocarp exocarp and mesocarp two biological replicates were employed for each stage and tissue for the samples, whereas one sample was used for each stage of mesocarp or exocarp. an additional technical replicate was carried out for exocarp representing separate itraq labeling reactions and analyses starting from the same protein sample. labeling of peptides with itraq reagents was performed according to the manufacturer's recommendations as follows. one hundred μg of each protein sample in a maximum volume of μl was precipitated overnight using 100% acetone and dissolved in μl of denaturing buffer containing μl denaturant and μl reducing reagent as provided in the itraq kit, followed by vortexing and incubation at 60°c for h. one μl of cysteine blocking solution was then added to each sample, followed by incubation at room temperature for min. these protein samples were digested overnight with trypsin at 37°c. itraq labeling was carried out by adding itraq reagents and to either the exocarp or the mesocarp samples representing the four developmental stages, green, pink/turning stage, red/fully turned, and purple, respectively. subsequently, these four samples were mixed by vortexing and further incubated at room temperature for h. the four itraq-labeled peptide samples were pooled together, diluted 1: with cation exchange sample buffer containing 25% acetonitrile in mm kh2po and then adjusted to ph using phosphoric acid. because of this acidification step, it is important to remove pectins prior to total protein extractions; we found in previous trials that the pectins likely polymerized and precipitated out of solution, converting samples mostly to a gelatinous state unsuitable for further analyses. the combined peptide mixture was fractionated by strong cation exchange chromatography on a biocad workstation, using a mm × cm polysulfoethyl aspartamide column. first, the mixed samples were loaded in buffer a at a flow rate of ml/min. once completely loaded, the column was washed for min with buffer a. peptides were eluted by a linear gradient of to mm kcl in buffer b. sixty-nine fractions were collected over the course of min at a flow rate of ml/min. of these fractions, only fractions containing the eluted labeled peptides as measured by optical density monitoring at nm were chosen for analysis on a h lc-ms/ms program. the fractionated samples were reduced to μl in a speed-vac and transferred to autosampler tubes. liquid chromatography and mass spectrometry the samples were analyzed for identification and quantitation on a qstar pulsar i hybrid tandem mass spectrometry system, fitted with a nano-electrospray ionization source using a μm fused silica emitter tip and interfaced with an integrated lc system consisting of a famos autosampler, switchos ii switching pump, and ultimate micropump. individual fractions containing peptides were injected onto a μm × cm c pepmap guard column, resolved using a μm × mm analytical column, and eluted using an automated binary gradient from 100% buffer a, % formic acid in h2o) to 40% buffer b in min, then from 40% to 80% buffer b for min. ms time of flight scans were acquired from m/z to for one second with up to two precursors selected for ms/ms from m/z to using information-dependent acquisition at seconds per scan; rolling collision energy was used to promote fragmentation. custom predicted tryptic peptide database a schema showing the pipeline for production of the predicted peptide database in support of this subsection is shown in figure all publicly available est data for each vitis species, including those from all v. vinifera cultivars, were downloaded in august as fasta files from the national center for biotechnology information. these data were parsed on the basis of reported vitis species of origin with the vast majority being from v. vinifera cultivars. since we were specifically interested in studying the proteome in v. vinifera cv. cabernet sauvignon pericarp tissue, an additional, more rigorous approach to the parsing of the cs ests was carried out in order to reduce or eliminate the potential for subsequent assembly of paralogous cs sequences into invalid contigs, thereby striving to strengthen the validity of protein identification in our itraq experiments. cs ests were obtained from the ncbi genbank database or from an in-house est project and subdivided into the following categories based upon the reported source tissues for the cdnas used for single pass sequencing: whole berry including seed, berry without seed, skin without seed or flesh, seed only, and other tissues including leaf, flower, tendril, and root. because the in-house ests were also present in the ncbi genbank database, the corresponding entries in genbank were removed since the genbank entries do not have sequence quality scores. the following files containing est data comprised each of the above mentioned groups: vv; ws; v. labrusca; v. pseudoreticulata; v. riparia; v. rotundifolia; v. shuttleworthii), cso, css, csp, cse, and csb. sequences were processed using cross_match and trim in order to remove vector sequences as well as ambiguous nucleotides at the sequence ends. to successfully perform the above cleanup analyses, phred quality scores were used where available; otherwise,'place-holder' quality scores were generated for any sequences for which no phred scores were available, as was the case for most of the ests in genbank. place-holder quality scores were also used later in the cluster assembly process as discussed in more detail, below. following the cross_match and trim processing, the sequences were further trimmed using perl scripts designed in-house to eliminate known invalid sequences and trim polya/t tails, if present in a given sequence. polya/t stretches were limited to bp in order to prevent subsequent chimeric contig assembly based on those repeats. if polya was followed by a > bp stretch of ac, at, gc, or gt repeats, the polya stretch was trimmed to bp and all sequence 3' to this was discarded; if polyt was preceded by a > bp stretch of ac, at, gc, or gt repeating sequence, the polyt stretch was trimmed to bp and all sequence 5' to this was discarded. if polya started at least two thirds of the est sequence length, it was trimmed to bp; if polyt started at less than one third of the est sequence, it was trimmed to bp. any part of a sequence that started or ended with > bp of repeats of ac, at, gc, or gt was deleted. if a sequence started or ended with'n's, the'n's were deleted and the corresponding quality scores were also removed. to better ensure that contig assemblies were based on high quality nucleotide sequence data, percent'n' content was determined for each sequence. if the percentage was >, the flanking bp regions where scanned for'n's and, if present, were trimmed to exclude the'n's, thereby lowering the total'n' percentage. sequences shorter than bp were trimmed to the first and last occurrences of an'n'. for resulting sequences longer than bp, the'n' percentage was recalculated and, if still > %, a record of the sequence was made. each of these sequences was then compared with other sequences in a combined dataset using blastn to determine its uniqueness. if a given sequence was already represented in the dataset by another sequence with a lower'n' content, the sequence in question was eliminated. the curated sequence datasets were next clustered using pcap software with parameters of 95% overlap identity and bp overlap length. pcap was used instead of cap in order to take advantage of parallelized processing. parallelization provided the ability to distribute each dataset assembly workload across cpus for significantly faster processing time. the pcap assembly program was modified and recompiled with est_flag set at. the pcap assembly step was followed by a series of post-assembly steps. we performed two clustering permutations in order to test the effects of database design on peptide identification using our itraq data. first, we clustered all sequences together to create the "as" database, including ws, vv, and all cs_ files; all sequences were weighted evenly. second, csb, cse, csp, cso, css, ws, and vv were clustered separately with higher weighting placed on cs sequences in the vv build and the original phred scores retained for the in-house cs sequences. weighting was accomplished by assigning higher quality scores such that when polymorphisms were encountered by pcap in an assembly, preference was given for selection of the cs nucleotide for the resulting contig. following assemblies, the generated contigs and singletons were merged into one file for each dataset. any sequences longer than bp were suspected to be chimeric, so they were parsed to a separate file, translated in all frames, and peptides with a minimal size of amino acids before a predicted stop codon were submitted to a blastx search against the nr database. the resulting multiple peptides predicted within long contigs were coded with "lc", as well as with "f" for the translational frame, with the frame number and the peptide number designated from among the multiple peptides. a blastx analysis was next performed on each contig and singleton sequence against the nr database in order to identify the best frame for subsequent in silico translation. the frame identified via blastx analysis was then used to generate the predicted orf for a given contig or singleton. in order to further curate predicted orfs, each was subjected to in silico cleavage at any'unknown' amino acid or stop codon and then compared to a similarly generated list of peptides from the corresponding best scoring protein sequence identified in the blastx search. the'best peptide' was then identified in the translation frame as the peptide with an exact match to the blastx peptide. if no such peptide could be identified, the longest peptide generated by in silico cleavage of the sequence at each occurrence of an unknown amino acid and/or stop codon was used. all sequences which resulted in "no hit found" in the blastx results were subsequently translated in all six frames and appended to the end of the'best peptide' file. in all cases where a six-frame translation was applied, the resulting peptides were cleaved in silico at every unknown amino acid and/or stop codon and only those sequences amino acids or longer were kept. the resulting list of'best peptides' for each of the sequences was then subjected to blastp analysis using the uniprotkb database in order to determine the sequence identity. the five highest blastp hits for each query sequence were aligned using an in-house perl script to identify putative n- and/or c-termini. if no consensus site could determined for an n- or c-terminus via alignment with similar sequences, then these sequences were trimmed at the tryptic digestion site nearest to the ends of the predicted orf to eliminate potentially truncated predicted tryptic peptides from the database. the parameters programmed into the scripts included: 1) blastp e-values, where e = h indicates a less significant hit and e = l indicates a stronger hit, 2) the difference in length of the vitis query sequence versus each of the top five subject sequences, and 3) the length of the exact match of amino acids to the top hit. these parameters improved automation of accurate predictions of methionine sites and identification of likely full-length amino acid sequences without requiring manual inspection of the blastp results. sequences identified with a predicted methionine at the n-terminus were coded with''. determinations of c termini were done based on a small range cutoff of amino acids between the stop codon in each top hit and the predicted stop codon in each corresponding vitis orf; if the difference was greater than two amino acids and deemed unclear, the c-terminal end of the predicted orf was trimmed at the nearest upstream tryptic cleavage site. the detailed process, above, was applied to each data set individually. in preparation for the merger of the datasets, the cs tissue-specific sequences were analyzed for uniqueness based on comparing each sequence to every other sequence and discarding all shorter sequences for each exact match. once all cs duplicate sequences were removed, this dataset was merged with the remaining two sets, vv and ws. this final set consisting of all of the sequences was then subjected to yet another uniqueness test where each sequence was compared to every other sequence but this time cs sequences where intentionally not removed, even if an exact duplicate existed in either the vv or the ws set and was of greater length. this allowed for preferential retention of the cs sequences in order to keep information about the tissue of origin of a detected protein. out of a total of sequences submitted to the uniqueness test, were identified with cs duplicates present due to the preferential retention of those sequences. from the resulting sequences, only those that started with a predicted methionine were then submitted for signalp analysis and processing which allowed for the identification and ultimate trimming of signal peptides giving rise to predicted mature protein sequences. those trimmed sequences that had a predicted cleavable targeting signal were coded''. after removal of the predicted signal peptides, a final uniqueness test was performed. analysis of ms/ms data itraq ms/ms data were analyzed using proteinpilot software v.. for both tryptic peptide identification and quantification. the peptides and corresponding relative abundances were obtained in proteinpilot using a confidence cutoff of >. database searching for each sample was done on predicted tryptic peptide sequence data using either the msdb database or in-house databases. annotations and annotated protein names indicated in proteinpilot output files were coded to indicate several parameters specific to the orf identified as well as the est or contig from which the orf sequence was predicted. itraq data representing the four ripening initiation stages in each of the three exocarp samples were combined into a single tab delimited file. likewise, itraq data representing each of the three mesocarp samples were combined into a second tab delimited file. duplicate entries among exocarp or mesocarp files were identified using an in-house script in the r environment with'custom orf id' as the search string. then, ratiometric data at each of the three comparisons using'green' as the reference stage were averaged prior to export for cluster analyses. entries with the same name but different template cdna sources were not averaged since these may represent isoforms from different source tissues and/or cultivars. we chose to cluster all proteins detected in the exocarp or mesocarp in order to capture all information on expression patterns detected, without restricting our analyses to only those proteins that were replicated amongst the individual exocarp or mesocarp files. k-means clustering into four partitions was carried out on ratiometric data for the exocarp and mesocarp files separately using multiexperiment viewer software. we used a -fold threshold for biological significance which was validated by consistencies between trends in protein expression presented here as increasing or decreasing with corresponding patterns of gene expression identified in previous publications.  protein detection and annotation in order to strengthen the ability of ms/ms spectra annotation software such as proteinpilot to accurately identify peptide sequences in complex total protein samples, we performed weighted clustering on a large est collection generated from vitis spp. to create the most cultivar-specific peptide map database possible for cabernet sauvignon. in addition, we created perl scripts in order to find the n and c termini in translated orfs or otherwise to carry out trimming at known tryptic digestion sites. this was done to help decrease the number of incorrectly annotated, non-tryptic n- and/or c-terminal peptides incorporated into protein abundance quantitations. to further increase the number of identifiable peptides, our in-house scripts searched all of the predicted orfs designated as containing a likely start methionine using signalp to identify and remove putative n-terminal signal peptides, resulting in a database containing mature protein sequence data for the predicted full length proteins. to assess whether the custom tryptic peptide database improved protein discovery, we examined the number of high confidence proteins and peptides identified by proteinpilot using this database versus a common approach to annotation, a search of the msdb. we performed this analysis on two itraq data sets, mesocarp 2005- and exocarp 2005- derived from four stages of cabernet sauvignon berries at ripening initiation. at a confidence level of 95%, proteins were identified in the mesocarp using the custom database, whereas only proteins were identified in a search of the same data set using the msdb. at a confidence level of 95%, proteins were identified in the exocarp using the custom database, whereas proteins were identified in a search of the same data set using the msdb. these results indicate that in these two itraq data sets, the use of the custom tryptic peptide database improved high confidence protein discovery by % and %. the number of high confidence peptides detected using the custom database was -fold and -fold higher in mesocarp and exocarp, respectively, in comparison to searches of the same itraq data sets using the msdb. the greater difference in the number of high confidence peptides versus proteins detected using the custom database in comparison to the msdb indicates that the most important impact of implementing the custom database was that it afforded the identification of more high confidence peptides per protein than could be achieved using the msdb. arelease hammersmith campus of imperial college london best clustering weighted higher scores for ests derived from cabernet sauvignon in order to avoid the potential loss of snps specific to this cultivar during consensus sequence determination by the pcap software cpredicted tryptic peptides that were determined to likely be truncated were trimmed and removed from the custom peptide database the effects of weighting and trimming on high confidence peptide and protein detection were analyzed using the exocarp itraq data set. we performed a second build of all vitis spp. ests using pcap. for the purposes of this comparative analysis, all ests in this second build were equally weighted for consensus sequence determination using an arbitrary phred score; no higher weighting of cabernet sauvignon-derived ests was used. this unweighted vitis est database was then used to generate a second predicted tryptic peptide database, including trimming of end-truncated peptides, as before. the number of high confidence peptides and proteins detected were similar using the orf databases created from the weighted versus the unweighted est databases, suggesting that the effect of simple nucleotide polymorphisms among the vitis spp. ests is ultimately negligible at the level of peptide identification via ms/ms in grapevine. the effectiveness of trimming predicted end-truncated peptides, on the other hand, was evident in the increased number of peptides identified by proteinpilot in searches of the trimmed orf database versus a third, untrimmed custom vitis orf database that we created for this comparison by omitting the trimming step shown in figure 1; both of these databases were weighted for cabernet sauvignon ests in the original pcap builds. these findings indicate that automated searching for the best predicted methionine start site and removing truncated n- and/or c-terminal tryptic peptide sequences in predicted orfs based upon est data increases the ability of ms/ms spectra annotation software to identify peptide sequences with high statistical confidence. technical and biological replication of expression ratios proteins identified in itraq data sets from two technical replicates were analyzed. two μg aliquots from the same total protein sample, exocarp were separately labeled with itraq reagents, subjected to nanolc-ms/ms, and the custom predicted tryptic peptide database was searched with the resulting spectra using proteinpilot. only proteins that were detected in each of the four ripening initiation stages were retained for the comparative analysis. in-house custom scripts written in r language were employed to detect protein entries that were duplicated in the two proteinpilot output files representing the technical replicates, using'custom orf id' as the search string. out of a total of proteins detected in both itraq data sets, or 29% of these proteins were detected in both files, whereas the remainder were unique to either technical replicate with proteins in one file and proteins in the other. we were interested in estimating consistency in trends in expression data along ripening initiation for each replicated protein entry that we detected. we determined that strict estimation of replicability via pearson's correlation coefficient or concordance correlation coefficient analyses of expression ratios did not provide relevant quantitative information, since even small deviations in expression ratios can result in poor correlative data. we considered that identifying consistent trends in protein accumulation above an arbitrary cutoff such as -fold was most relevant for the purposes of discovering protein candidates for hypothesis formulation related to ripening control. as an alternate approach to correlation estimation, a replication scoring system was devised, whereby expression ratios identified by manual inspection as differing by less than log for each duplicate protein entry for at least two of the three ratios received a score of expression ratios differing by between log and log for at least two of the three ratios received a score of and expression ratios differing by log or greater for at least two of the three ratios received a score of the mean replication score for the duplicated protein entries in the technical replicates was determined to be ±; 17% of the protein entries received a score of these data indicate that if a -fold change or greater is considered as biologically relevant, over 80% of duplicate protein entries but only approximately one quarter of all proteins detected in two technical replicates would be expected to be identified as significantly changing along ripening initiation with similar quantitative trends, based on the technical replicates analyzed here. proteins identified in itraq analyses of two biological replicates, exocarp and exocarp 2005- were next analyzed in an identical manner to the two technical replicates. out of a total of proteins detected in the two itraq data sets, or 33% of these proteins were detected in both files, whereas the remainder were unique to either biological replicate with proteins in one file and proteins in the other. the same replication scoring system was implemented for proteins detected in both itraq data sets as was done for the technically replicated proteins, above. the mean replication score for the proteins common to both biological replicates was determined to be ±; 12% of protein entries received a score of to further explore the biological significance of replicability analyses, proteins common to either the technical replicates or the biological replicates were compared and it was determined that protein entries were common to the technical replicates and the biological replicates, i.e. were detected in all three data sets analyzed for replication; we considered this approach as valid since all data analyzed here were derived from the exocarp. the technical and biological replication results indicate that the'shotgun' lc-ms/ms approach employed here with itraq-labeled exocarp total proteins is capable of repeatedly identifying approximately one third of the total number of proteins, whether they are technical or biological replicates. expression profiling of proteins along grape berry ripening initiation four clusters were generated for the exocarp in which trends in protein accumulation were increasing strongly, increasing gradually, not changing significantly, or decreasing, respectively, from the green through the fully turned ripening initiation stages. additional file lists proteins by cluster number along with corresponding log2-transformed ratiometric data for each protein entry. ratiometric data were calculated relative to the green stage for each of the pink/turning, red, and purple/fully turned stages. several protein isoforms significantly increasing along ripening initiation were identified with annotated functions in anthocyanin flavonoid biosynthesis and storage, defense ), and cell expansion. conversely, several components of the photosynthetic machinery were identified as decreasing greater than -fold along ripening initiation, which is consistent with a reduction in photosynthesis at this period of berry development. we mined the exocarp data for proteins that were increasing in abundance relative to the green stage and annotated as enzyme or transporter components of pathways leading to hypothesized regulators of ripening initiation in grapes, aba, glucose, and br. a putative lytb protein increased -fold in abundance and is responsible for the last step of the plastidic pathway to isopentenyl diphosphate, leading in part to the production of the plant hormones, aba and gibberellic acid. other proteins of this mep pathway were also detected but only as slightly increasing in abundance along ripening initiation. an isopentenyl-diphosphate δ-isomerase i protein was detected as increasing 2-fold along ripening initiation; this enzyme controls a major early step in isoprenoid biosynthesis and is likely localized to the chloroplast. we identified one component specific to the aba biosynthetic pathway, a protein similar to violaxanthin de-epoxidase from tea, which was stably expressed along ripening initiation. cytosolic idp, potentially also formed in the plastid and exported to the cytoplasm, is incorporated in the biosynthesis of brs. a putative grapevine ortholog to a br biosynthetic protein from pea, pslkb,, was identified as increasing -fold, peaking at the third stage tested. a grapevine hexose transporter was identified as increasing -fold, which is consistent with hexose accumulation during ripening initiation. few proteins annotated with signal transduction functions were detected in any of the four clusters that could hypothetically be involved in aba, glucose, and/or br signaling. abscisic stress response protein was previously implicated in cross-talk between aba and glucose signaling; expression data were variable, with some isoforms increasing up to 4-fold along ripening initiation and others showing stable expression. a putative pirin protein of unknown function increased -fold along ripening initiation. a putative ortholog to the malus spp. ttg wd repeat protein, which has previously been implicated in arabidopsis in the regulation of anthocyanin biosynthesis, was detected as increasing, indicating that vvttg could play a similar role in the grape exocarp; a caveat is that the confidence level for vvttg unique peptide detection was less than 95%. curiously, we did not identify peptides representing vvmyba which is highly expressed at the transcriptional level during ripening initiation and encodes a transcription factor previously demonstrated to positively regulate anthocyanin biosynthesis genes in grape. four clusters were generated for the mesocarp in which trends in protein accumulation were increasing strongly, increasing gradually, not changing significantly, or decreasing, respectively, from the green through the fully turned ripening initiation stages. additional file lists proteins by cluster number along with corresponding log2-transformed ratiometric data for each protein entry, relative to the green stage. similar to the exocarp data, we identified > 2-fold accumulation of several protein isoforms annotated with functions in cell enlargement, fruit softening and defense concomitant with a significant reduction in proteins inhibitory to fruit ripening. we mined the mesocarp data for proteins annotated as enzyme or transporter components of pathways leading to aba, glucose, and brassinosteroid accumulation. we detected vvnced increasing over 2-fold along ripening initiation, which represents expression of a key committed step specifically in aba biosynthesis in the plastid. a protein annotated as aba glucosyltransferase increased -fold along ripening initiation in the mesocarp. as in the exocarp, an isopentenyl-diphosphate δ-isomerase i was detected as increasing 2-fold along ripening initiation, plus, a putative cytosolic isoform was detected as stably expressed. farnesyl diphosphate synthase, a cytoplasmic enzyme leading to the br biosynthetic pathway via squalene biosynthesis, was detected in the mesocarp as increasing -fold along ripening initiation. a putative grapevine ortholog to the br biosynthetic protein in pea, pslkb,, and in cotton, ghdwarf, was detected as increasing -fold. the hexose transporter, vvht was detected as increasing 2-fold, similar to its expression pattern that was detected in the exocarp. we mined the mesocarp data further to identify proteins annotated with signal transduction functions. no plasma membrane receptor candidates were identified as increasing but several were identified as not changing or decreasing along ripening initiation, most notably a receptor-like kinase similar to perk1-like protein from rice that decreased -fold at the onset of pigment accumulation to -fold at the purple stage; this perk1-like protein was also identified as significantly decreasing in the exocarp. the same pirin protein identified in the exocarp was similarly detected as increasing in the mesocarp along ripening initiation. inconsistent itraq data were obtained for vvasr in the mesocarp with isoforms showing an increase, stable expression, or strong down-regulation. a custom script written in the r environment was used to search for proteins common to the combined exocarp and mesocarp files using'custom orf id' as the search string. protein isoforms detected in common between exocarp and mesocarp are shown in additional file protein isoforms detected only in exocarp are shown in additional file protein isoforms detected only in mesocarp are shown in additional file isoforms of anthocyanin flavonoid biosynthetic proteins were only detected in the exocarp, which is consistent with induction timing and tissue localization of these pigments during ripening initiation and demonstrates the efficacy of our separation of exocarp and mesocarp for proteomic analyses. additional file shows the source tissues and number of unique high confidence peptides per protein for those proteins indicated here in the results section and in figure discussion the itraq data obtained with exocarp and mesocarp total proteins confirmed previous ratiometric transcript abundance data for key components of aba and br biosynthesis, as well as the influx and accumulation of sugars during ripening initiation. we confirmed for the first time the strong accumulation of several components of anthocyanin biosynthesis at the protein level, including a putative anthocyanin o-methyltransferase. the detection of vvnced accumulation confirmed previous real-time rt-pcr data and further supports a role specifically for this nced family member in aba biosynthesis in the mesocarp during berry ripening initiation. based on the sensitivity limitations to the shotgun proteomic technique employed here, we cannot conclude that vvnced accumulation, activity, and, consequently, aba biosynthesis are localized to the mesocarp; vvnced may be expressed in the exocarp but we did not detect it. the previous detection of ζ-carotene desaturase transcripts in the exocarp argues against a model in which aba is synthesized in the mesocarp and transported to the exocarp whereupon it activates, in part, anthocyanin biosynthesis. if the primary site of aba production in the developing fruit is the seed, however, and aba transported from the seed to the surrounding pericarp is an early signaling event initiating ripening, as we have previously hypothesized, it is reasonable that the first tissue in the pericarp in which aba is synthesized is the mesocarp, given its closer proximity to the seed than the exocarp. it remains to be determined whether a gradient in aba biosynthesis during early ripening stages proceeds from the seed through the mesocarp and, finally, in the exocarp. the presence and possible moderate accumulation of a putative aba glucosyltransferase indicates that an attenuating mechanism for the aba signal may operate during ripening initiation to control aba homeostasis. the accumulation of the pirin protein is also intriguing, given this protein's previously demonstrated interaction with g-protein alpha in arabidopsis seeds and a potential role in modifying aba action via negative feedback control. we used a relatively new quantitative ms/ms-based approach, itraq, to advance our understanding of differential protein expression underlying non-climacteric ripening initiation in grape berries. the itraq approach offered several advantages over 2dge methods for protein discovery, including increased detection sensitivity based on our findings reported here in comparison to previous reports on grape berry proteomics. by using strong cation exchange and reverse phase column microcapillary chromatography coupled with nanospray ms/ms detection with total protein extracts from grape berries, we were able to resolve three-fold or more proteins per sample than would be expected using 2dge. one current limitation to an ms-based proteomic approach with grapevine is that there are no finished genome sequence data for grapevine, although two projects are undertaking assembly and annotation from which a high quality orfeome database can eventually be derived. although there are over vitis spp. ests deposited in genbank, v. vinifera is a highly heterozygous species; therefore, we considered that it might be important to weight builds of these ests via manipulation of phred scores so as to favor sequence data corresponding to our cultivar of interest, cabernet sauvignon, when snps were encountered by pcap, where applicable for a given contig assembly. although we determined that weighting ests to the genotype of interest for est assembly provided no clear advantages in this study, we conclude that producing a tryptic peptide database targeted to vitis sequences and including removal of predicted truncated peptides improved protein detection and annotation. furthermore, our findings indicate that a tryptic peptide database based on finished pinot noir whole genome sequence data will be valid to implement for proteome studies with any v.vinifera cultivar or vitis species, with the exception of cases where deletions have occurred in the pinot noir homozygous line. we chose to not include genome sequence data available for v. vinifera cv. pinot noir due to significant gaps in current assemblies and the potential for inaccurate automated gene predictions. until grapevine genome sequence assembly and annotation are finished, we propose that the predicted orf database presented here will be of value to the grapevine community in two significant ways. while gaps exist in the genome sequence assemblies, the protein database presented here may provide information for'missing' proteins either not yet predicted from the pinot noir genome sequence data and/or other vitis spp. not represented in the pinot noir genome sequence data, e.g. due to chromosomal deletions. further, protein predictions based upon expressed sequences represent'real' proteins, so our database could potentially be used to validate orf predictions based on whole genome sequence data alone. while itraq labeling coupled with nanolc-ms/ms proved overall to be an advance over 2dge in sensitivity and quantitation, consistent detection of predicted protein sequences between technical or biological replicates from grape exocarp was limited. consistency in trends in ratiometric data along ripening initiation for those proteins that were detected in replicate exocarp samples was further limited. inconsistent ratiometric data for some proteins detected in both biological replicates may represent differences in expression due, for example, to variability in seasonal growing conditions and not technical variability. nonetheless, limited replicable detection of proteins among biological samples for technical reasons is commonly encountered with liquid chromatography and ms-based proteomics and could arise from variation in preparation of total proteins and/or itraq-labeled peptides from sample to sample. underlying the apparent variation arising from sample preparations are constraints on detection imposed by the mass spectrometer, both with respect to the dynamic range of the instrument and the nature of selections of peptides in the first ms by the ms software for export to the collision cell prior to amino acid detection via the second ms. furthermore, our finding that two thirds of the replicated exocarp proteins were detected both in technical and biological replicates may have reflected the higher probability of detecting these mainly abundant, housekeeping-type proteins in any given total protein sample using a shotgun approach. over-sampling of abundant peptides in digested total protein samples is a limitation of the shotgun approach to quantitative proteomics and likely precluded our ability to discover more proteins annotated with signal transduction functions that could regulate aba, br, and hexose responses during ripening initiation. in order to increase detection sensitivity of low abundance regulatory proteins using shotgun proteomics techniques, it will likely be helpful to isolate membrane and nuclear proteins separately from cytoplasmic proteins prior to digestions. affinity chromatography of berry protein extracts using antibodies directed against abundant proteins such as thaumatin detected in both exocarp and mesocarp may also improve detection sensitivity for low abundance proteins by selectively removing these proteins prior to itraq labeling steps. similarly, advances in detection sensitivity in devices such as the fourier transform ion cyclotron resonance ms should allow us to delve deeper into the grape berry proteome in order to better understand the molecular control of non-climacteric ripening in this species. nonetheless, the searchable quantitative data presented here on over proteins detected in two separate berry tissues along four ripening initiation stages can serve as a framework public resource contributing towards our understanding of the dynamic grapevine proteome.  we determined that a predicted peptide database can be derived from grapevine est data using advanced clustering and trimming approaches and successfully implemented for quantitative proteome profiling. we demonstrated in grapevine that by implementing a predicted peptide database targeted to an organism of interest, a significant increase in the number of high confidence peptides identified and annotated from ms/ms data was gained in comparison to more commonly used msdb searches. furthermore, by using a shotgun quantitative proteomics approach in combination with a targeted predicted peptide database, we showed that greater numbers of high confidence peptides were detected than would be expected using 2d gel electrophoresis techniques. we verified for the first time at the protein level, quantitative expression patterns for components of isoprenoid and flavonoid metabolism important for grape compositional chemistry, as well as identified new signal transduction candidate proteins associated with non-climacteric ripening initiation in grape berries. abbreviations aba: abscisic acid; blast: basic local alignment and search tool; br: brassinosteroids; est: expressed sequence tag; itraq: isobaric tagging for relative and absolute quantitation; orf: open reading frame; pcap: parallel contig assembly program authors' contributions jl designed the experiments, carried out protein extractions and data analyses, assisted in designing the scripts for the predicted peptide database, and drafted the manuscript. ml designed the scripts and carried out programming to produce the predicted peptide databases presented here. ds carried out mass spectrometry and database searches using the itraq ms/ms data. stl performed primary data analyses, assisted in drafting the manuscript, and was responsible for overseeing the project. supplementary material additional file 1 duplicated proteins and associated replication scores for exocarp technically replicated itraq experiments. ratiometric data are given in log scale. click here for file additional file 2 duplicated proteins and associated replication scores for exocarp and exocarp 2005- biologically replicated itraq experiments. ratiometric data are given in log scale. click here for file additional file 3 proteins identified in all of the three itraq data sets analyzed and exocarp 2005-1) and tested for replication in expression trends along ripening initiation. click here for file additional file 4 exocarp proteins listed for each k-means cluster along with corresponding log2-transformed ratiometric data for each protein entry. exocarp proteins from the 2005- and 2005- datasets were combined into a single file. duplicate entries were identified using an in-house script in the r environment with'custom orf id' as the search string. then, ratiometric data at each of the three comparisons using'green' as the reference stage were averaged prior to export for k-means cluster analyses. click here for file additional file 5 mesocarp proteins listed for each k-means cluster along with corresponding log2-transformed ratiometric data for each protein entry. mesocarp proteins from the 2005- and 2005- datasets were combined into a single file. duplicate entries were identified using an in-house script in the r environment with'custom orf id' as the search string. then, ratiometric data at each of the three comparisons using'green' as the reference stage were averaged prior to export for k-means cluster analyses. click here for file additional file 6 protein isoforms detected in common between the exocarp and mesocarp. a custom script written in the r environment was used to search for proteins common to the combined exocarp and mesocarp files using'custom orf id' as the search string. click here for file additional file 7 protein isoforms detected only in exocarp samples. proteins excluded from additional file that were detected in the combined exocarp file are shown. click here for file additional file 8 protein isoforms detected only in mesocarp samples. proteins excluded from additional file that were detected in the combined mesocarp file are shown. click here for file additional file 9 protein source tissues and number of unique high confidence peptides per protein. proteins presented in the results and discussion sections as well as in figure are shown. the source tissue for each protein hit shown in column a corresponds to a proteinpilot output file, which can be downloaded at the number of unique > 95% confidence peptides determined by the proteinpilot software for each protein is shown in column e. ratiometric data shown in columns g, h, and i are relative to the green stage and are not log2-transformed. click here for file additional file 10 reference key for coding terms used in column a and column b in additional files through the key is provided as a printable rapid reference for data mining in additional files through as well as in the protein database available online, as shown in the acknowledgments. click here for file acknowledgements the authors wish to thank monica elliot for assistance with protein preparations and itraq labeling. we thank dr. fred y. peng for assistance with itraq data processing and cluster analyses. the authors gratefully acknowledge funding for this research from genome canada and project management support from genome british columbia as part of the genome canada-genoma españa collaborative research and development initiative. the authors also wish to thank the province of british columbia and the bc wine industry for their continued support of research activities in the university of british columbia wine research centre. resources and raw data developed through the research presented here can be downloaded freely at: these are: 1) scripts used for the customized grapevine protein database creation, 2) the protein database, and 3) the proteinpilot output files for exocarp exocarp 2005- exocarp 2005- mesocarp mesocarp 2005- and mesocarp 2005- proteinpilot software can be obtained by contacting applied biosystems. the scripts for protein database creation will need to be modified prior to implementation with a species different to v. vinifera; the corresponding author may be contacted for additional information.
cabernet sauvignon, and lacking truncated n- and c-terminal fragments. we determined that a predicted peptide database for ms/ms applications can be derived from est data using advanced clustering and trimming approaches and successfully implemented for quantitative proteome profiling. we used itraq ms/ms in conjunction with custom peptide db searches to quantitatively characterize several important pathway components for berry ripening previously described at the transcriptional level and confirmed expression patterns for these at the protein level. by searching itraq ms/ms data generated from berry exocarp and mesocarp samples at ripening initiation, we determined that implementation of the custom database afforded a large improvement in high confidence peptide annotation in comparison to the msdb. for species for which finished genome sequence data are not available, identification of proteins relies on similarity to proteins from other species using comprehensive peptide map databases such as the msdb. itraq is a proteomics technique that uses isobaric tags for relative and absolute quantitation of tryptic peptides. we developed a bioinformatic pipeline for processing est data in order to produce a predicted tryptic peptide database specifically targeted to the wine grape cultivar, vitis vinifera cv. we were interested in characterizing ripening initiation in grape berries at the protein level in order to better define the molecular control of this important process for grape growers and wine makers. quantitative shotgun proteome profiling holds great promise for characterizing biological processes such as fruit ripening initiation and may be further improved by employing preparative techniques and/or analytical equipment that increase peptide detection sensitivity via a shotgun approach.
 itraq is a proteomics technique that uses isobaric tags for relative and absolute quantitation of tryptic peptides. in proteomics experiments, the detection and high confidence annotation of proteins and the significance of corresponding expression differences can depend on the quality and the species specificity of the tryptic peptide map database used for analysis of the data. for species for which finished genome sequence data are not available, identification of proteins relies on similarity to proteins from other species using comprehensive peptide map databases such as the msdb.  we were interested in characterizing ripening initiation in grape berries at the protein level in order to better define the molecular control of this important process for grape growers and wine makers. we developed a bioinformatic pipeline for processing est data in order to produce a predicted tryptic peptide database specifically targeted to the wine grape cultivar, vitis vinifera cv. cabernet sauvignon, and lacking truncated n- and c-terminal fragments. by searching itraq ms/ms data generated from berry exocarp and mesocarp samples at ripening initiation, we determined that implementation of the custom database afforded a large improvement in high confidence peptide annotation in comparison to the msdb. we used itraq ms/ms in conjunction with custom peptide db searches to quantitatively characterize several important pathway components for berry ripening previously described at the transcriptional level and confirmed expression patterns for these at the protein level.  we determined that a predicted peptide database for ms/ms applications can be derived from est data using advanced clustering and trimming approaches and successfully implemented for quantitative proteome profiling. quantitative shotgun proteome profiling holds great promise for characterizing biological processes such as fruit ripening initiation and may be further improved by employing preparative techniques and/or analytical equipment that increase peptide detection sensitivity via a shotgun approach. 