HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Full Text (PDF) Supplementl Data All Versions of this Article: M700440-MCP200v1 7/5/927    most recent Submit a response Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Request Permissions Glossary Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Colzani, M. Articles by Quadroni, M. Search for Related Content PubMed PubMed Citation Articles by Colzani, M. Articles by Quadroni, M. Social Bookmarking                      What's this?

Molecular & Cellular Proteomics 7:927-937, 2008.
© 2008 by The American Society for Biochemistry and Molecular Biology, Inc.

---

Research

Relative Protein Quantification by Isobaric SILAC with Immonium Ion Splitting (ISIS)^*,S

Mara Colzani^,, Frédéric Schütz^¶, Alexandra Potts, Patrice Waridel and Manfredo Quadroni^,||

From the Center for Integrative Genomics, University of Lausanne, Quartier Sorge, 1015 Lausanne-Dorigny, Switzerland and ^¶ Bioinformatics Core Facility, Swiss Institute of Bioinformatics, Quartier Sorge, 1015 Lausanne, Switzerland

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Metabolic labeling techniques have recently become popular tools for the quantitative profiling of proteomes. Classical stable isotope labeling with amino acids in cell cultures (SILAC) uses pairs of heavy/light isotopic forms of amino acids to introduce predictable mass differences in protein samples to be compared. After proteolysis, pairs of cognate precursor peptides can be correlated, and their intensities can be used for mass spectrometry-based relative protein quantification. We present an alternative SILAC approach by which two cell cultures are grown in media containing isobaric forms of amino acids, labeled either with ¹³C on the carbonyl (C-1) carbon or ¹⁵N on backbone nitrogen. Labeled peptides from both samples have the same nominal mass and nearly identical MS/MS spectra but generate upon fragmentation distinct immonium ions separated by 1 amu. When labeled protein samples are mixed, the intensities of these immonium ions can be used for the relative quantification of the parent proteins. We validated the labeling of cellular proteins with valine, isoleucine, and leucine with coverage of 97% of all tryptic peptides. We improved the sensitivity for the detection of the quantification ions on a pulsing instrument by using a specific fast scan event. The analysis of a protein mixture with a known heavy/light ratio showed reliable quantification. Finally the application of the technique to the analysis of two melanoma cell lines yielded quantitative data consistent with those obtained by a classical two-dimensional DIGE analysis of the same samples. Our method combines the features of the SILAC technique with the advantages of isobaric labeling schemes like iTRAQ. We discuss advantages and disadvantages of isobaric SILAC with immonium ion splitting as well as possible ways to improve it.

Many chemical labeling schemes have been devised. In turn, they can be subdivided in two groups based on the approaches needed for quantification. Most techniques rely on measuring the intensity of light/heavy parent ions in MS survey scans to establish intensity ratios. One technical issue with such an approach is the need to establish an unambiguous correlation (based on m/z and retention times in liquid chromatography) between tandem MS spectra used for identification and their precursor peaks in MS survey scans. ICAT (8) and a multitude of other labeling schemes (9–12) belong to this class.

The second group is formed by techniques that exploit quantitative information embedded in tandem MS spectra rather than in survey scans (13, 14). Earlier examples of this principle used protease-mediated ¹⁶O/¹⁸O labeling (15–17). Such approaches, however, require co-fragmentation of the precursors, and this necessitates very close or identical masses, such as in the iTRAQ¹ scheme. The latter is a special case that uses isobaric tagging reagents that fragment to give, for every peptide, diagnostic low mass ions used for quantification. Compared with methods with distinct mass precursors, these techniques present the advantages that no precursor mass splitting results in higher signal intensity and data handling is easier because only MS/MS data are needed, thus eliminating the need to integrate peaks along chromatographic runs. Appropriate design of the reagent has allowed the reporter ion mass to be varied, and this has been used to implement 4- and 8-fold multiplexing work flows.

By contrast, classical SILAC (18) is a metabolic labeling method similar, from the point of view of the quantification procedure, to non-isobaric precursor methods. To achieve accurate quantification, it needs high resolution mass data with low noise levels together with a powerful software able to correlate MS/MS identifications with full scan MS data and extract integrated ion intensities for the precursors and cognate, but often unidentified labeled analogues.

We propose an alternative SILAC scheme that combines some useful features of iTRAQ-like work flows with the convenience of SILAC for studying samples from cultured cells. The technique is based on the incorporation of isobaric labeled amino acids that generate distinct residue immonium ion fragments. We show that reliable quantification is possible purely on MS/MS spectra using very simple data extraction tools.

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Cell Culture and Labeling—
Unlabeled cell extracts were obtained from BJAB human B cells grown in standard RPMI 1640 medium supplemented with 10% fetal bovine serum and penicillin-streptomycin (all from Invitrogen).

Metabolic labeling was performed on two human melanoma cell lines: SBCL2, derived from an radial growth phase melanoma (19) and kindly provided by Dr. Giovanella (Stehlin Foundation for Cancer Research, Houston, TX), and SKMel28, derived from skin metastasis (20) and obtained from the ATCC collection (number HTB-72, LGC Promochem). SBCL2 and SKmel28 cells were grown in high glucose Dulbecco's modified Eagle's medium deficient in L-valine, L-leucine, L-isoleucine, and L-glutamine (Cell Culture Technologies, Gravesano, Switzerland). The medium was reconstituted from concentrated stocks according to the manufacturer's instructions, sterile filtered, and stored at 4 °C. L-Isoleucine, L-leucine, and L-valine labeled either with 1-¹³C (99% 1-¹³C enrichment) or ¹⁵N (98% ¹⁵N enrichment) were obtained from Cambridge Isotope Laboratories (Andover, MA). The powdered amino acids were diluted in PBS to obtain stock solutions of 40 mg/ml for valine and 15 mg/ml for leucine and isoleucine that were sterile filtered and stored at –80 °C. SKMel28 and SBCL2 cells were grown in the base medium supplemented with 4 mM L-glutamine (Amimed, Allschwil, Switzerland), 10% dialyzed fetal bovine serum (Invitrogen), penicillin-streptomycin (Invitrogen), and the amino acids 1-¹³C- (for SKMel28) or ¹⁵N-labeled (for SBCL2 and SKMel28 as specified under "Results") L-valine, L-leucine, and L-isoleucine. The final concentrations of Val, Leu, and Ile were 94, 105, and 105 mg/liter, respectively, corresponding to the standard composition of the Dulbecco's modified Eagle's medium. Cell lines were grown for 3 weeks to achieve complete labeling.

Preparation of Protein Samples—
Cells at about 80% confluence were washed twice with PBS to remove serum proteins. Cells were detached by trypsinization and washed twice in PBS. Pellets of cells were resuspended in ice-cold hypotonic buffer (10 mM HEPES, 1.5 mM MgCl₂, and 10 mM KCl) and sonicated for three cycles of 5 s each on ice. The lysates were centrifuged twice (10 min at 13,000 rpm at 4 °C) to pellet cellular debris; each time the supernatant was recovered. The protein concentration of the second supernatant was determined using the Bradford protein assay (Bio-Rad); the relative protein concentration of samples to be mixed was verified by densitometry on whole lanes of SDS-PAGE gels after Coomassie staining. For mixing experiments, lysates were combined in a known ratio equal to 1:3 for ¹⁵N:¹³C SKMel28 extracts or 1:1 for ¹³C SKMel28:¹⁵N SBCL2. For the co-elution experiment, ¹⁵N:¹³C SKMel28 lysates were combined in 1:2 ratio.

2D Gel Electrophoresis—
2D DIGE experiments were performed according to the manufacturer's instructions (GE Healthcare) (21). Briefly 300 µg of protein from ¹⁵N SBCL2 and ¹³C SKMel28 cell lysates were precipitated with acetone overnight at –20 °C; the pellet was dissolved in 25 µl of DIGE Cell Lysis Buffer (30 mM Tris/HCl, pH 8.5, 7 M urea, 2 M thiourea, 4% (w/v) CHAPS). 50 µg of protein from ¹⁵N SBCL2 and ¹³C SKMel28 lysates were minimally labeled with 800 pmol of reconstituted Cy3 and Cy5 dyes, respectively. 25 µg of the two samples were mixed together, labeled with Cy2, and used as internal standard. The Cy2-, Cy3-, and Cy5-labeled samples were combined, supplemented with 2% pH 3–11 ampholytes (IPG Buffer, GE Healthcare), and applied by cup-loading on a 13-cm, pH 3–11 non-linear Immobiline DryStrip (GE Healthcare) previously rehydrated in DIGE Cell Lysis Buffer overnight. Isoelectric focusing was performed on an Ettan IPGphor (GE Healthcare) system to attain a total of 22,500 V-h. Prior to the second dimension, strips were equilibrated for 15 min in a reducing buffer containing 6 M urea, 2% (w/v) SDS, 30% (v/v) glycerol, 32 mM DTT, 100 mM Tris, pH 8. This was followed by a 15-min alkylation in a buffer containing 6 M urea, 2% (w/v) SDS, 30% (v/v) glycerol, 240 mM iodoacetamide, 100 mM Tris, pH 8. Second dimension migration was carried out on a 13 x 9-cm 8–16% Criterion precast gel (Bio-Rad) at a constant voltage of 80 V for 4 h. The gel was rinsed in water and then scanned using a Molecular Imager FX (Bio-Rad) scanner. Gel images were analyzed using ImageMaster 2D Platinum DIGE Software version 5.0 (GE Healthcare). The gel was Coomassie-stained overnight and destained with 10% acetic acid. Spots of interest were manually excised from the gel.

1D Gel Electrophoresis for Protein Digestion—
For 1D gel electrophoresis, 60 µg of protein from either pure or mixed lysates were subjected to limited electrophoretic separation on a 10% SDS-PAGE minigel, i.e. the migration was stopped after the front had moved by about 2.5 cm into the separating gel at which point the bands of a prestained size marker were visible and separated in the 20–250-kDa range. After Coomassie staining, each lane was cut into four fractions corresponding to regions of different molecular weights.

In-gel Protein Digestion with Trypsin—
Slices excised from 1D or 2D gels were transferred to 96-well plates. In-gel proteolytic cleavage with sequencing grade trypsin (Promega, Madison, WI) was performed automatically in a ProGest robotic work station (Genomic Solutions, Ann Arbor, MI) according to a described protocol (22) The liquid supernatant from the digestion was recovered and concentrated by evaporation. The final volume was adjusted to 60 µl with 0.1% TFA in 2% acetonitrile.

LC-MS/MS—
Digests were desalted using C₁₈ StageTips (Proxeon, Odense, Denmark), which have an estimated loading capacity equal to 10 µg. The peptides were eluted from the microcolumns with 80% acetonitrile after which the eluate was dried and resuspended in 7 µl of solvent A (2% acetonitrile, 0.5% formic acid). 2.5 µl of purified sample were injected on a reversed-phase C₁₈ column (PepMap100, 3 µm, 100 Å, LC Packings) and separated by nanoflow liquid chromatography on an Ultimate (LC Packings) system on line with an electrospray quadrupole-time-of-flight mass spectrometer (API QSTAR Pulsar i, Applied Biosystems/SCIEX, Concord, Ontario, Canada). The gradient used for separation was from 2 to 40% acetonitrile in water with 0.5% formic acid at a flow rate of 200 nl/min.

The mass spectrometer was controlled by the Analyst QS 1.1 software set to operate in information-dependent acquisition mode to automatically switch between MS and collision-induced dissociation MS/MS. Survey full-scan MS spectra were acquired from 400 to 1200 m/z in 1 s after which the two most intense ions with charge 2+ to 4+ were isolated for fragmentation. In a second method (method B), the three most intense ions with 2+ to 4+ charge were fragmented.

Each parent ion underwent two different MS/MS acquisitions. The first scan (1-s duration) spanned the 50–1200 m/z mass range, used a collision energy proportional to precursor mass, and was intended to collect sequence information. The second MS/MS scan (0.2-s duration) was focused on (iso)leucine and valine immonium ions and therefore covered narrow ranges of masses (71–74 and 85–88 m/z for method "A" or a single scan at 64–95 m/z for method "B") to achieve a better sampling of the reporter ions. The collision energy for this narrow range scan was fixed at 60 eV to promote a higher degree of fragmentation, and the masses were acquired in the "enhanced" mode that allows a considerable gain in signal intensity by specifically pulsing ions of a narrow mass range (23). Three (method A) or two (method B) MS/MS spectra were cumulated for each selected peptide; former target peptides and their isotopes were dynamically excluded for 90 s (with a mass tolerance of 50 milli-mass units).

All runs were performed in triplicate for the ¹⁵N SKMel28:¹³C SKMel28 1:3 mixture fraction using method A. The ¹⁵N SBCL2:¹³C SKMel28 1:1 mixture was analyzed in quadruplicate using both methods A and B in duplicate.

For the analysis of 2D DIGE spots by LC-MS a shorter gradient from 2 to 40% acetonitrile in 45 min was used. All MS parameters were the same as those described above.

For testing peptide co-elution, a special acquisition method was used in which two precursors were fragmented for 10 consecutive scans, giving a total cycle time of 21 s. The elution profile of the peak was then followed to monitor immonium ion fragment intensities during the same cycle.

Protein Identification and Quantification—
The Mascot.dll script (version 1.6b21) supplied by Matrix Science was used to extract spectra of each run using the following parameters: 2+ to 4+ charge state, peak centroiding, no deisotoping, spectra rejected if containing less than eight peaks, no minimum peak intensity required, report peak area. For mixing experiments, Mascot generic format flat text (.mgf) files deriving from different fractions of the same mixture were pooled before performing the database search. Proteins were identified using Mascot version 2.1 (Matrix Science, London, UK) and searching the UniProt database, restricted to human taxonomy. Database releases used were 8.7 of September, 19 2006 for experiments with unlabeled samples (71,233 sequences after taxonomy filter) and release 11.0 of May, 29 2007 (70,004 sequences after taxonomy filter) for experiments with labeled samples. One trypsin missed cleavage was allowed, cysteine carbamidomethylation was set as fixed modification, and methionine oxidation was defined as variable modification. Unless otherwise specified, peptide mass tolerance was set to 1.2 Da (0.6 Da for unlabeled samples), and the fragment mass tolerance was set to 0.3 Da. We used a broad peptide mass tolerance (1.2 Da) to include in Mascot version 2.1 searches the +1-Da ¹³C precursor peptides that can be incorrectly chosen during peak detection. For labeled samples, ¹³C labeling and ¹⁵N labeling on valine, isoleucine, and leucine were also set as fixed peptide modifications unless specified. The Mascot multidimensional protein identification technology (MudPIT) scoring system was used for all analyses of mixtures (the unlabeled sample was analyzed with standard scoring) with significance threshold at p = 0.05 and minimum ion score equal to 14. The same .mgf files used for protein identification mentioned above were searched against a decoy database automatically generated by Mascot version 2.2 and containing random sequences having the same length and average amino acid composition of the forward "normal" database to compute the false discovery rate (FDR) of protein identification. Significance threshold was set at p = 0.05, and minimum ion score was set at 14 as in the previous search used for protein identification; the resulting FDR values are reported under "Results" for the different shotgun analyses. Proteins spanning the same set of peptides, or a subset, were collapsed into a single entry on the hit list according to the principle of parsimony. For these cases, only one representative protein is listed. Proteins with matched peptides in common with other sequences were validated only if matched by at least one unique ("bold red") discriminating peptide. No further attempt was made to discriminate protein isoforms that were indistinguishable with the available mass spectrometry data.

To obtain protein quantification, an in-house built Perl (version 5.2.2.0) script extracted the intensity values of the following masses from the Mascot .dat result file: 72.081 amu for [¹³C]valine, 73.078 amu for [¹⁵N]valine, 86.097 amu for [¹³C](iso)leucine, and 86.094 amu for [¹⁵N](iso)leucine with a tolerance of ±0.020 amu. For each identified protein, these intensities were associated to the corresponding valine- and (iso)leucine-containing peptides. In the case of peptides containing both valine and (iso)leucine, the quantification was considered as independent. Only Val-containing peptides with non-zero values on both 72.081 and 73.078 m/z and Ile/Leu-containing peptides with non-zero values on both 86.097 and 87.094 m/z were used for estimation of abundance ratios. R software (24) (2.3.1 version) was used for all statistical calculations. For each protein presenting at least four quantified peptides, a straight line was fitted by least squares through the points represented by the ¹³C and ¹⁵N immonium ion intensities for each peptide; the corresponding slope is an estimate of the abundance ratio (25). The standard error of the estimate, the t value of the test for significance of the slope, and the corresponding p value as well as the 95% confidence interval were also calculated from the linear regression. Outlier points, which could result either from experimental variability or from the presence of differentially expressed isoforms of the same protein matched by the same set of peptides, were not removed before linear regression. Thus no attempt was made to discriminate such protein isoforms. The regression line was fitted twice using in turn the ¹³C and ¹⁵N immonium ion intensities as response variable; the two estimations of the abundance ratio were extremely close, and the one that produced the smallest residuals was used. Other methods were considered in particular orthogonal regression, which produced extremely similar results (data not shown). Because the ratio ¹⁵N/¹³C should in theory remain constant for all peptides regardless of the intensities of the individual ions, the regression line was constrained to pass through the origin. Examination of the plots for proteins with a large number of peptides (the first 30 Mascot protein matches in the 1:1 ¹³C SKMel28:¹⁵N SBCL2 data set) suggested that this assumption is correct. Unless otherwise specified, the reported ratios are defined heavy/light (H/L).

	RESULTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Protein Labeling Strategy—
Isobaric SILAC with immonium ion splitting (ISIS) differential labeling is based on the incorporation of isobaric analogues of specific amino acids capable of providing distinct immonium ions fragments whose intensity can be relatively quantified in MS/MS spectra. Essential amino acids containing the stable 1-¹³C (culture A) or -¹⁵N (culture B) isotope are added to a cell culture medium lacking these same amino acids but containing all others in unlabeled form (Fig. 1A). Cells are then grown for a time sufficient to obtain a near complete (98% or more) protein labeling, i.e. six or more cell divisions. Every labeled amino acid residue (be it 1-¹³C- or -¹⁵N-labeled) in a peptide or protein increases its nominal mass by 1 amu. Corresponding peptides originating from the cultures A and B should therefore be isobaric, but they are expected to generate upon collision-induced fragmentation different immonium ions. These residue ions are produced by the simultaneous cleavage of both the peptide bond N-terminal to the residue of interest as well as of the C–CO bond within the residue (26). Residues labeled with ¹³C on the carbonyl carbon (position 1) thus produce "light" immonium ions identical to those generated by naturally occurring (unlabeled) residues because the labeled C-1 atom is lost during MS/MS fragmentation. On the other hand, amino acid residues labeled on the carbon or ¹⁵N-labeled on the backbone nitrogen are expected to produce immonium ions 1 amu heavier than the natural ones because the label is retained after CID.

View larger version (20K): [in this window] [in a new window]   FIG. 1. A schematic of ISIS strategy. a, ISIS labeling. Isobaric peptides containing either 1-13C- or 15N-labeled amino acids produce in MS/MS spectra light and heavy immonium ions, respectively. R = CH(CH3)2 (valine), CH2CH(CH3)2 (leucine), or CH(CH3)CH2CH3 (isoleucine). During immonium ion formation, 1-13C-labeled amino acids lose the label and produce an ion 1 amu lighter than the one deriving from 15N-amino acids. b, ISIS analytical steps. Cell lines are metabolically labeled using 1-13C and 15N isobaric variants of the three essential amino acids valine, leucine, and isoleucine. The samples are mixed together after lysis, digested, fractionated by SDS-PAGE, and analyzed by LC-MS/MS. After a survey scan, each precursor undergoes two different MS/MS scans, a broad m/z range scan for protein identification and a narrow range scan focused on the immonium ion region, to detect light and heavy peak intensities. After database search, reporter immonium ion intensities are extracted from CID spectra and used to calculate heavy/light protein ratios.

Preliminary Analysis on Unlabeled Sample—
Among the 10 amino acids essential for mammalian cells, we selected Val, Leu, and Ile for labeling. These three amino acids are among the most abundant in proteins, and their frequency in the Swiss-Prot database is 6.72, 9.66, and 5.88%, respectively. In addition, they all generate characteristic intense immonium ions upon CID. We choose to use the combination of these three amino acids to increase the likelihood of obtaining at least one labeled residue in each peptide generated by digestion. The immonium ions of these amino acids in their natural or 1-¹³C-labeled form are localized at 72.081 m/z for valine and at 86.097 m/z for leucine and isoleucine. ¹⁵N-labeled analogues of the same amino acids are expected to give signals at 73.078 and 87.094 m/z, respectively.

Before trying to exploit immonium ion intensities for quantification purposes, it was necessary to verify that (i) ions at 72.081 and 86.097 m/z are detected in MS/MS spectra of unlabeled peptides and are reliable diagnostic signals for the presence of Val, Ile, and Leu and that (ii) no naturally occurring residue produces strong signals at 73.078 and 87.094 amu so that these "channels" can be used for the specific quantification of heavy labeled peptides. To this end, an unlabeled complex mixture of proteins from human cells was digested and analyzed by LC-MS/MS, generating a collection of 4759 spectra. Immonium ion intensities were then extracted from the resulting centroided spectra with an in-house built peak matching script; proteins identified by at least two unique peptides were analyzed. Intense signals at 72.081 (Val) and 86.097 m/z (Ile/Leu) were found to correlate with the presence of the corresponding amino acid (supplemental Table 1 and supplemental Fig. 1). Much weaker peaks were present at 73.078 and 87.094 m/z, and their intensity was proportional to those at 72.081 and 86.097 m/z, respectively (supplemental Table 1 and supplemental Fig. 2). In fact, the correlation analysis by linear regression of the intensities of the peaks pairs localized at 72.081/73.078 m/z and at 86.097/87.094 m/z showed that, for each peptide, the heavier peaks were on average 5.6 and 6.5% (R² = 0.97) of the light peaks, respectively, thus very close to the percentage expected from the theoretical isotope distribution (4.8% for Val and 6.0% for Leu/Ile). This suggests that the peaks at 73.078 and 87.094 m/z were produced by the naturally occurring heavier (mostly ¹³C₁) isotopic forms of the Val and Ile/Leu immonium ions.

We also explored the eventuality that the natural immonium ion of asparagine (which has a calculated mass of 87.056 amu) could interfere with the quantification based on the [¹⁵N](iso)leucine immonium ion because two peaks are localized at close m/z values. To this end, we repeated the correlation analysis relaxing from ±0.02 to ±0.2 amu the tolerance for the extraction of intensities at 87.094 m/z to include Asn immonium ion signals. As a result, a linear regression with a poor correlation coefficient was obtained (R² = 0.63) (supplemental Fig. 3); when Asn-containing peptides were removed from the data set, the correlation improved dramatically (R² = 0.96), and the regression line was almost identical to the one obtained using a narrow mass tolerance (supplemental Fig. 2B). This demonstrates that, at least in some sequences, Asn can generate strong immonium ion signals that could bias the quantification; this was avoided using a reduced tolerance (±0.02 amu) for peak matching, which effectively filtered out Asn immonium ions. Analysis of a ¹⁵N-labeled sample showed that at our working resolution (4200 at this mass) two distinct peaks were detected at 87.055 and 87.100 m/z, confirming that it was possible to discriminate between these two species (Fig. 2B).

View larger version (12K): [in this window] [in a new window]   FIG. 2. MS spectra and reporting immonium ion regions in unlabeled and heavy and light labeled peptides. a, MS spectra of the peptides GYSFTTTAER and QEYDESGPSIVHR from unlabeled cells (Un) and cells labeled with [1-13C]Val, [1-13C]Ile, and [1-13C]Leu (L) or [15N]Val, [15N]Ile, and [15N]Leu (H). The peptide GYSFTTTAER (containing neither Val, Ile, nor Leu) shows no mass shift following labeling. The QEYDESGPSIVHR peptide contains instead two labeled amino acids; thus both 13C and 15N labeling produce a mass shift of +2 amu (+0.66 m/z for this triple charged precursor). b, cumulative MS/MS spectra (sum of about 300 spectra) of immonium ions regions for unlabeled cells (Un) and cells labeled with 1-[13C]Val, [13C]Ile, and [13C]Leu (L) or [15N]Val, [15N]Ile, and [15N]Leu (H).

Labeling Efficiency and Isobaricity—
To test protein labeling, SKMel28 and SBCL2 cells were grown in medium containing either the 1-¹³C- or the -¹⁵N-labeled forms of the three amino acids (Val, Ile, and Leu), respectively. The two cultures were grown in identical condition for 3 weeks, a time considered sufficient to achieve the complete incorporation of labeled analogues into proteins (18). A third culture was maintained in unlabeled medium as negative control. Lysates for labeled and unlabeled cells were fractionated by SDS-PAGE, digested, and analyzed by LC-MS/MS on a QQ-TOF instrument. Tandem mass spectra were then submitted to database search by Mascot. By comparing the molecular weights of identical peptides identified from labeled versus unlabeled cells, it was possible to verify that each 1-¹³C- or ¹⁵N-amino acid analogue added 1 amu to the original mass and that 1-¹³C- and ¹⁵N-labeled peptides were isobaric (Fig. 2A). However, we noticed a difference in the isotope distribution of peptide signals in ¹⁵N samples with a higher abundance of heavier isotope species; this phenomenon and its implications are addressed under "Discussion." No unlabeled peaks were noticed during manual inspection of MS spectra. To comprehensively check labeling efficiency, a database search was performed on a large data set (mixture of SKMel28 and SBCl2 cells, see below) without selecting ISIS labels as fixed modification, i.e. with native peptide masses and a precursor mass tolerance of 0.3 Da. Over a total of 23,784 spectra submitted to Mascot, only 267 were matched, identifying 54 proteins (FDR above identity threshold equal to 17%). The majority of these matched spectra (193, 72%) contained no Val, Ile, or Leu and were thus expected to have unmodified masses. All the remaining 74 matched spectra were either false positives or contaminants (keratin and trypsin). On the other hand, when ISIS labels were used as fixed modification in Mascot, 12,796 spectra were matched using the same search parameters, identifying 604 proteins (FDR above identity threshold equal to 2.9%); the percentage of sequences not containing Val, Ile, and Leu was in this case much smaller (only 2.7%). This confirmed that the labeling was highly efficient.

1-¹³C- and ¹⁵N-labeled peptides were isobaric within the limits of measurement error. Their MS/MS spectra showed identical series of b and y ions (supplemental Fig. 4). In fact, the nature of the labeling is such that the only fragments expected to be split in ESI spectra after CID fragmentation are immonium and a ions; as expected, split light/heavy a ions of significant intensity were observed only for peptides containing the labeled amino acids localized within the first two N-terminal residues. Finally the masses of peptides that did not contain Val, Ile, or Leu were unchanged after labeling (Fig. 2A).

To highlight occurrence of immonium ions in a complex peptide mixture, we cumulated data by adding several hundreds of MS/MS spectra obtained from the correspondent fraction of unlabeled and [¹³C]Val/Ile/Leu- and [¹⁵N]Val/Ile/Leu-labeled mixture of proteins (Fig. 2B). In the immonium ion region of MS/MS spectra, intense peaks at 72.081 and 86.097 m/z were visible in the 1-¹³C-labeled sample, whereas equally intense signals at 73.078 and 87.094 m/z were visible in ¹⁵N-labeled extracts (Fig. 2B). Two peaks were detected at m/z 87.055 and 87.100, corresponding to the natural immonium ion of Asn and the "heavy" immonium ion of Leu/Ile. Weaker peaks at 72.081 and 86.097 m/z were visible in the ¹⁵N-labeled extracts probably due to residual amounts of unlabeled Val, Ile, and Leu amino acids in the peptides.

The co-elution of 1-¹³C- and ¹⁵N-labeled precursor peptides, a prerequisite for the accurate estimate of light/heavy ratios, was confirmed using a 2:1 mixture of completely labeled SKMel28 lysates. Because the precursors are isobaric, their elution profiles are indistinguishable in MS mode and can only be discriminated upon fragmentation. Therefore the sample was analyzed using a modified LC-MS/MS method that acquired CID spectra of the same precursor for 20 s per cycle; this cycle time covered the typical chromatographic peak width, which was estimated to be between 20 and 30 s. Although this led to the identification of a limited number of peptides, the intensity profiles of the light and heavy immonium ions could be monitored and plotted. An example is shown in supplemental Fig. 5 and demonstrates co-elution for all precursors detected in this analysis.

Protein Quantification: Method Validation with a Normalized Mixture—
We then tested the ability of the method to reliably quantify proteins in a complex environment by analyzing a fraction composed of 1-¹³C- and ¹⁵N-labeled lysates of SKMel28 cells mixed in a known ratio (3:1) before SDS-PAGE fractionation. The fraction of proteins with a molecular mass between 40 and 80 kDa was selected and analyzed as described under "Materials and Methods." After data submission to Mascot, 127 proteins were identified with a total of 2035 spectra matched. Of these, 1978 (97%) were matched to sequences containing at least one of the three amino acids used for labeling. The false positive rate calculated for peptide matches above identity threshold was 3.6%. Intensity values for the 72.081/73.078 and 86.097/87.094 m/z peak pairs were extracted from the Mascot .dat file (containing the matched spectra in centroided form) and associated to the corresponding peptides using an in-house programmed script; intensity pairs where one of the values was zero were discarded as specified under "Material and Methods." This yielded a total of 1149 peptides; 537 were quantified on valine, 1040 were quantified on (iso)leucine, and 428 peptides were quantified on both residues. Only the proteins presenting a minimum of four measurements were quantified to increase the confidence level of data (supplemental Fig. 6). Over a total of 127 proteins, it was thus possible to compute 80 protein quantifications (63%) as shown in Fig. 3A. Using a linear model, the mean heavy/light estimated ratio was equal to 0.361, thus very close to the expected value of 0.333. Confidence intervals were particularly small for proteins characterized by a large number of peptides and thus showing higher Mascot scores. The ratio distribution was spread from a minimum value of 0.155 to a maximum of 0.472; 50% of the ratios were between 0.307 and 0.379. The deviation from the mean ratio could reflect both the biological variability in the protein content of the cell (the two samples originated from two distinct cultures) and the experimental error; the latter was likely larger for the proteins identified with a lower Mascot score, which usually correlates with low signal intensity. Ratio values were transformed to a log base 2 scale to express -fold changes and to normalize the magnitude of under- and overexpressed proteins (Fig. 3A). This transformation yielded an approximately normally distributed data set.

View larger version (17K): [in this window] [in a new window]   FIG. 3. Quantification of a protein mixture with known ratio. Shown are heavy/light protein ratios obtained from a 40–80-kDa fraction of proteins from 1-13C- and 15N-labeled SKMel28 cells mixed in a 3:1 ratio. A, the upper plot shows the H/L (15N:1-13C) ratios for the 80 quantified proteins with the relative error bars representing a confidence interval of 95%. In the lower graph, the distribution of protein ratios is visualized by histogram and kernel density in log2 scale to normalize the magnitude of under- and overexpressed proteins. The box plot at the top shows the interquartile range containing 50% of the ratios between 0.307 and 0.379 (–1.706 and –1.399 in log2 scale), the median at 0.343 (–1.542 in log2 scale), (I), and the mean at 0.340 (–1.577 in log2 scale) (x). B, the upper plot shows high confidence protein ratios for the 65 proteins of 80 presenting a t value >9 and p value <0.001; error bars represent a confidence interval of 95%. In the lower graph, the protein ratio distribution is visualized as described for a. The median and the mean are at 0.350 (–1.511 in log2 scale, |) and 0.346 (–1.546 in log2 scale, x), respectively.

Differential Analysis—
The ISIS method was applied to the differential analysis of proteins in SKMel28 and SBCL2 human cell lines, both derived from epithelial melanoma but at different stages of tumor progression. Because of their similar background, we expected an overall similar pattern of protein expression. On the other hand, these two cell lines display some visible differences in morphology and in adhesion behavior, so we speculated that some differences should be clearly detectable in their proteomes. SBCL2 and SKMel28 cells were differentially labeled and lysed. Extracts were mixed in 1:1 ratio and fractionated by 1D PAGE into four fractions, which were digested by trypsin and analyzed in replicates with ISIS-specific scan cycles as described under "Materials and Methods." We identified 582 proteins (13,504 peptides) with a false positive rate of identification of 2.5% (calculated for peptide matches above Mascot identity threshold). Over 582 proteins, we obtained 527 heavy/light ratios (91%). After applying the statistical cutoff values described in the previous section, 502 high confidence estimates were found (Fig. 4). The mean value was 0.921, thus close to the expected value of 1.00, whereas 50% of the ratios were between 0.761 and 1.009, likely reflecting the normal biological variability in the protein content of the cells. The box plot highlights the presence of both high and low "outlier" ratios deviating from the expected normal distribution, indicating proteins differentially expressed in the two cell lines. We arbitrarily decided to consider as significant only the proteins at least 2-fold over- or underexpressed, therefore presenting heavy/light ratios larger than 1.858 (more abundant in SBCL2) or smaller than 0.465 (more abundant in SKMel28). With these thresholds we identified 13 proteins more abundant in SBCL2 and 13 other proteins more abundant in the SKMel28 cell line. These proteins are listed in Table I.

View larger version (11K): [in this window] [in a new window]   FIG. 4. Relative quantification of proteins in two melanoma cell lines by ISIS. Shown are a histogram and kernel density estimate for the 502 high confidence H/L protein ratios obtained from a 1:1 mixture of 1-13C SKMel28 and 15N SBCL2 cells. The protein ratios are visualized in log2 scale to express -fold change and to normalize the magnitude of under- and overexpressed proteins. The box plot at the top shows that 50% of the ratios are between 0.761 and 1.009 (–0.395 and +0.013 in log2 scale), the median is 0.873 (–0.195 in log2 scale, |), and the mean is 0.929 (–0.188 in log2 scale, x) as well as the presence of outlier ratios with a -fold change >2.0.

View larger version (9K): [in this window] [in a new window]   FIG. 5. DIGE versus ISIS quantification of 2D gel spots from two melanoma cell lines. Shown is a scatter plot of DIGE versus ISIS H/L ratios for 21 spots obtained from 2D PAGE (empty circles) and for eight well defined spots (filled circles) with the respective regression lines and coefficients of determination. Values are in log2 scale.

View larger version (12K): [in this window] [in a new window]   FIG. 6. Histogram of "shotgun" ISIS versus "spot" ISIS quantifications for proteins present in both shotgun and spot analyses. The * highlights proteins whose spot quantification differs less than 20% from the shotgun quantification. Cofilin-1 was detected in four spots on the 2D gel; all spot ISIS ratios are compared with the unique global ratio for cofilin-1 obtained by shotgun ISIS. PEBP, phosphatidylethanolamine-binding protein; TPIS, triosephosphate isomerase; HNRPD, heterogeneous nuclear ribonucleoprotein D.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

One intriguing phenomenon we observed was that peptides from [¹⁵N]Val/Ile/Leu-labeled cells consistently displayed isotope distribution with higher abundances of the +1, +2, and +3 isotope peaks. This fact was also observed in peptides not containing Val/Ile/Leu (Fig. 2A) but was not detected in samples labeled with 1-¹³C-amino acids. This suggests that the amino nitrogen from Val, Ile, or Leu can be utilized by cells to be incorporated (although probably with low frequency) in other amino acids. Essential amino acids like Val, Ile, and Leu cannot be synthesized by mammalian cells, which, however, possess enzymes for their degradation. The first step in the catabolism of these three amino acids is transamination that transfers the NH₂ group to glutamate, which in turn can be used to synthesize other amino acids, thus probably leading to a low rate of incorporation of ¹⁵N elsewhere in the peptide sequence. This phenomenon, however, should not affect quantifications obtained with the ISIS method because only fragments produced by the supplied essential amino acids are measured. In our experiments the masses of precursors were not significantly modified, and we found that the only precaution to be taken was to use a sufficiently large (at least 2.5 amu) isolation window to capture all isotopic peaks for MS/MS.

One of the critical issues during development of this method was sensitivity. Immonium ions are not always strong peaks in MS/MS spectra acquired under CID conditions set to fragment peptides for protein identification. As a comparison, the chemical structure of the iTRAQ reporter ions has been designed to provide efficient fragmentation and high ionization efficiency for a higher sensitivity of detection. However, we found that the intensity of the immonium ion peaks can be significantly increased using a fast "zoom" scan with higher collision energy (60 eV) and by exploiting the capabilities of a pulsing instrument (QSTAR) to maximize the transmission efficiency of small fragment ions. The quantification scans of narrow range obtained were then combined with the CID spectra of full range before database search and quantification.

The data analysis work flow was simple; after protein identification by Mascot, an in-house built script extracted immonium ion peak intensities from the Mascot .dat file and associated them to the corresponding peptides. The relative protein quantification together with the associated statistics was computed with the freeware software R, using a linear regression method, for the proteins presenting at least four quantified peptides. Because the quantification in ISIS work flow is in principle identical to that performed with iTRAQ data, both freeware and recent commercial software packages (e.g. iTRACKER and Mascot 2.2) can be easily adapted to process ISIS data.

The ISIS strategy was tested on a complex mixture of fully labeled proteins with known mixing ratio, showing good specificity and sensitivity and providing reliable protein quantification. 97% of the peptide sequences were found to contain at least one of the amino acids used for labeling, allowing the quantification of 63% of the identified proteins. The confidence intervals for the computed ratios were small, especially for the protein characterized by intense peaks or a high number of peptides. Using this data set, p and t cutoff values were set to detect high confidence ratios and discard the "low confidence" ratios. The implemented method was then used for a preliminary analysis of the proteomes of two different melanoma cell lines at different stages of tumor progression and allowed the identification and the quantification of differentially expressed proteins. Over a total of 582 proteins identified from replicate analyses, it was possible to obtain 502 high confidence ratios (86%). 26 proteins were found to be differentially expressed in these two cell lines. Interestingly a number of proteins involved in cytoskeletal modeling and organization (plectin-3, myosin-9, annexin-A1, calpain-2, and A- and B-filamins) were found to be preferentially expressed in SBCL2 cells. The differential expression of this group of proteins could at least partially explain some differences in morphology that can be observed between SBCL2 and SKMel28 cells.

To further validate our method, we performed a DIGE analysis of the same lysate mixtures, and we compared the protein ratios obtained using DIGE with those from ISIS quantification. Although some difficulties were experienced due to the identification of more proteins in one spot and to the uncertainty of spot boundaries in crowded regions of the gel, the two approaches generally provided similar results. Whenever a 2D gel spot contained more than one protein, the ISIS quantification could determine two distinct ratios, providing more accurate quantitative data. On the other hand, the ISIS shotgun experiment without two-dimensional separation provided information for the total level of identified proteins, without distinction for the possible isoforms produced by post-translational modifications, as usual for shotgun mass spectrometry-based techniques. The ISIS quantification obtained in the shotgun analysis detected a maximum difference in protein expression equal to 4.5-fold (Table I). On the other hand, the ratio values obtained from ISIS quantification of 2D gel spots spanned a broader range, detecting a maximum ratio value of 15.2 (Figs. 5 and 6). This indicates that the ISIS quantification method can potentially detect changes in protein expression of at least 15-fold.

Comparison with Available Techniques—
The ISIS methodology potentially combines some typical advantages of metabolic labeling schemes with the ease of MS/MS-based quantification techniques such as iTRAQ for studying samples from cultured cells. Metabolic labeling allows mixing the differentially labeled samples just after cell lysis, thus eliminating sample-to-sample variations due to subsequent manipulations. This is especially important whenever complex multistep protein purification must be carried out and explains why SILAC and similar techniques are widely applied for the study of subcellular compartments or protein complexes. Although the number of amino acid analogues required for ISIS is high (six for the scheme described here), the reagents are significantly cheaper (about half the price for the same volume of medium) than those commonly used for classical SILAC ([¹³C₆]lysine and [¹³C₆]arginine). One distinctive advantage of the ISIS scheme is that when Val, Ile, and Leu are used for labeling, a considerable fraction of all peptides (37%) can be quantified twice through both channels, thus significantly increasing the number of data points and the confidence of quantification. Also in comparison with classical SILAC, the ISIS scheme (like iTRAQ) has the advantage that the spectra complexity is not increased by doubling of all peaks, which also results in a gain in sensitivity. The isobaric nature of the ISIS tags allows MS/MS-based quantification and very simple data analysis independent of native instrument data formats. Finally in contrast with the commercial SILAC, ISIS labeling should be easily applicable to analyze non-tryptic peptides and could be a viable solution to study fractions enriched with membrane or acidic proteins, which necessitate cleavage with proteases different from trypsin or chemical cleavage.

One main disadvantage of the ISIS method relative to the iTRAQ technique is the difficulty to modify the technique to perform higher order (>2 samples) multiplexing. One general limitation is the need of an instrument that is able to produce high levels of precursor fragmentation and transmission of very low mass ions.

Perspectives—
Some features of ISIS method can be still improved; for example the data extraction procedure can be further explored to obtain a better determination of peak areas by optimizing centroiding parameters. Furthermore natural isotopic distributions (¹³C content of natural amino acids) and isotopic impurities of reagents should be taken into account during quantification, and this should lead to less biased calculations especially for high ratios. Moreover we observed that, as expected, the peptide ratios obtained with low ¹³C and ¹⁵N ion intensities were much more variable than ratios obtained from intense immonium peaks. This indicates that a weighted linear regression may be appropriate to down-weight the least precise measurements (the results described here were produced using an unweighted regression). Also the use of orthogonal regression, instead of the least squares, could produce a more accurate estimate of the ratios and the corresponding standard errors. Finally as for the iTRAQ technique, an implementation of the ISIS protocol could benefit from the use of the latest generation of instruments, such as Orbitrap with an additional octopole collision cell for accurate measurement of low mass ions or QQ-TOF with intelligent data acquisition routines for the optimization of data acquisition as a function of signal intensity and signal/noise ratio.

	FOOTNOTES

 
Received, September 14, 2007, and in revised form, November 9, 2007.

Published, MCP Papers in Press, December 28, 2007, DOI 10.1074/mcp.M700440-MCP200

¹ The abbreviations used are: iTRAQ, isobaric tags for relative and absolute quantitation; SILAC, stable isotope labeling with amino acids in cell cultures; ISIS, immonium ion splitting; 2D, two-dimensional; 1D, one-dimensional; FDR, false discovery rate; H/L, heavy/light; QQ-TOF, quadrupole/quadrupole TOF.

^* The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

^S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material.

To whom correspondence may be addressed. Tel.: 41-21-692-56-76; Fax: 41-21-692-57-05; E-mail: mara.colzani{at}unil.ch

^|| To whom correspondence may be addressed. Tel.: 41-21-692-56-76; Fax: 41-21-692-57-05; E-mail: manfredo.quadroni{at}unil.ch

	REFERENCES

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 

1. Tao, W. A., and Aebersold, R. (2003 ) Advances in quantitative proteomics via stable isotope tagging and mass spectrometry. Curr. Opin. Biotechnol. 14, 110 –118[CrossRef][Medline]

2. Chen, X., Sun, L., Yu, Y., Xue, Y., and Yang, P. (2007 ) Amino acid-coded tagging approaches in quantitative proteomics. Expert Rev. Proteomics 4, 25 –37[CrossRef][Medline]

3. Ong, S. E., and Mann, M. (2005 ) Mass spectrometry-based proteomics turns quantitative. Nat. Chem. Biol. 1, 252 –262[CrossRef][Medline]

4. Brand, M., Ranish, J. A., Kummer, N. T., Hamilton, J., Igarashi, K., Francastel, C., Chi, T. H., Crabtree, G. R., Aebersold, R., and Groudine, M. (2004 ) Dynamic changes in transcription factor complexes during erythroid differentiation revealed by quantitative proteomics. Nat. Struct. Mol. Biol. 11, 73 –80[CrossRef][Medline]

5. Blagoev, B., Kratchmarova, I., Ong, S. E., Nielsen, M., Foster, L. J., and Mann, M. (2003 ) A proteomics strategy to elucidate functional protein-protein interactions applied to EGF signaling. Nat. Biotechnol. 21, 315 –318[CrossRef][Medline]

6. Gruhler, A., Olsen, J. V., Mohammed, S., Mortensen, P., Faergeman, N. J., Mann, M., and Jensen, O. N. (2005 ) Quantitative phosphoproteomics applied to the yeast pheromone signaling pathway. Mol. Cell. Proteomics 4, 310 –327[Abstract/Free Full Text]

7. Sachon, E., Mohammed, S., Bache, N., and Jensen, O. N. (2006 ) Phosphopeptide quantitation using amine-reactive isobaric tagging reagents and tandem mass spectrometry: application to proteins isolated by gel electrophoresis. Rapid Commun. Mass Spectrom. 20, 1127 –1134[CrossRef][Medline]

8. Gygi, S. P., Rist, B., Gerber, S. A., Turecek, F., Gelb, M. H., and Aebersold, R. (1999 ) Quantitative analysis of complex protein mixtures using isotope-coded affinity tags. Nat. Biotechnol. 17, 994 –999[CrossRef][Medline]

9. Goodlett, D. R., Keller, A., Watts, J. D., Newitt, R., Yi, E. C., Purvine, S., Eng, J. K., von Haller, P., Aebersold, R., and Kolker, E. (2001 ) Differential stable isotope labelling of peptides for quantitation and de novo sequence derivation. Rapid Commun. Mass Spectrom. 15, 1214 –1221[CrossRef][Medline]

10. Munchbach, M., Quadroni, M., Miotto, G., and James, P. (2000 ) Quantitation and facilitated de novo sequencing of proteins by isotopic N-terminal labelling of peptides with a fragmentation-directing moiety. Anal. Chem. 72, 4047 –4057[Medline]

11. Schmidt, A., Kellermann, J., and Lottspeich, F. (2005 ) A novel strategy for quantitative proteomics using isotope-coded protein labels. Proteomics 5, 4 –15[CrossRef][Medline]

12. Zhou, H., Ranish, J. A., Watts, J. D., and Aebersold, R. (2002 ) Quantitative proteome analysis by solid-phase isotope tagging and mass spectrometry. Nat. Biotechnol. 20, 512 –515[CrossRef][Medline]

13. Thompson, A., Schafer, J., Kuhn, K., Kienle, S., Schwarz, J., Schmidt, G., Neumann, T., Johnstone, R., Mohammed, A. K., and Hamon, C. (2003 ) Tandem mass tags: a novel quantification strategy for comparative analysis of complex protein mixtures by MS/MS. Anal. Chem. 75, 1895 –1904[Medline]

14. Ross, P. L., Huang, Y. N., Marchese, J. N., Williamson, B., Parker, K., Hattan, S., Khainovski, N., Pillai, S., Dey, S., Daniels, S., Purkayastha, S., Juhasz, P., Martin, S., Bartlet-Jones, M., He, F., Jacobson, A., and Pappin, D. J. (2004 ) Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol. Cell. Proteomics 3, 1154 –1169[Abstract/Free Full Text]

15. Heller, M., Mattou, H., Menzel, C., and Yao, X. (2003 ) Trypsin catalyzed ¹⁶O-to-¹⁸O exchange for comparative proteomics: tandem mass spectrometry comparison using MALDI-TOF, ESI-QTOF, and ESI-ion trap mass spectrometers. J. Am. Soc. Mass Spectrom. 14, 704 –718[CrossRef][Medline]

16. Wang, Y. K., Ma, Z., Quinn, D. F., and Fu, E. W. (2001 ) Inverse ¹⁸O labelling mass spectrometry for the rapid identification of marker/target proteins. Anal. Chem. 73, 3742 –3750[Medline]

17. Schnolzer, M., Jedrzejewski, P., and Lehmann, W. D. (1996 ) Protease-catalyzed incorporation of ¹⁸O into peptide fragments and its application for protein sequencing by electrospray and matrix-assisted laser desorption/ionization mass spectrometry. Electrophoresis 17, 945 –953[CrossRef][Medline]

18. Ong, S. E., Blagoev, B., Kratchmarova, I., Kristensen, D. B., Steen, H., Pandey, A., and Mann, M. (2002 ) Stable isotope labeling by amino acids in cell culture, SILAC, as a simple and accurate approach to expression proteomics. Mol. Cell. Proteomics 1, 376 –386[Abstract/Free Full Text]

19. Hsu, M. Y., Shih, D. T., Meier, F. E., Van Belle, P., Hsu, J. Y., Elder, D. E., Buck, C. A., and Herlyn, M. (1998 ) Adenoviral gene transfer of β3 integrin subunit induces conversion from radial to vertical growth phase in primary human melanoma. Am. J. Pathol. 153, 1435 –1442[Abstract/Free Full Text]

20. Satyamoorthy, K., DeJesus, E., Linnenbach, A. J., Kraj, B., Kornreich, D. L., Rendle, S., Elder, D. E., and Herlyn, M. (1997 ) Melanoma cell lines from different stages of progression and their biological and molecular analyses. Melanoma Res. 7, Suppl. 2, S35 –S42[Medline]

21. Unlu, M., Morgan, M. E., and Minden, J. S. (1997 ) Difference gel electrophoresis: a single gel method for detecting changes in protein extracts. Electrophoresis 18, 2071 –2077[CrossRef][Medline]

22. Wilm, M., Shevchenko, A., Houthaeve, T., Breit, S., Schweigerer, L., Fotsis, T., and Mann, M. (1996 ) Femtomole sequencing of proteins from polyacrylamide gels by nano-electrospray mass spectrometry. Nature 379, 466 –469[CrossRef][Medline]

23. Chernushevich, I. V., Loboda, A. V., and Thomson, B. A. (2001 ) An introduction to quadrupole-time-of-flight mass spectrometry. J. Mass Spectrom. 36, 849 –865[CrossRef][Medline]

24. R Development Core Team (2005 ) R: a Language and Environment for Statistical Computing , R Foundation for Statistical Computing, Vienna, Austria

25. Pan, C., Kora, G., Tabb, D. L., Pelletier, D. A., McDonald, W. H., Hurst, G. B., Hettich, R. L., and Samatova, N. F. (2006 ) Robust estimation of peptide abundance ratios and rigorous scoring of their variability and bias in quantitative shotgun proteomics. Anal. Chem. 78, 7110 –7120[Medline]

26. Falick, A. M., Hines, W. M., Medzihradszky, K. F., Baldwin, M. A., and Gibson, B. W. (1993 ) Low-mass ions produced from peptides by high-energy collision-induced dissociation in tandem mass spectrometry. J. Am. Soc. Mass Spectrom. 4, 882 –893[CrossRef]

27. Ong, S. E., Kratchmarova, I., and Mann, M. (2003 ) Properties of ¹³C-substituted arginine in stable isotope labelling by amino acids in cell culture (SILAC). J. Proteome Res. 2, 173 –181[CrossRef][Medline]

28. Lafaye, A., Labarre, J., Tabet, J. C., Ezan, E., and Junot, C. (2005 ) Liquid chromatography-mass spectrometry and ¹⁵N metabolic labelling for quantitative metabolic profiling. Anal. Chem. 77, 2026 –2033[Medline]

CiteULike

Complore

Connotea

Del.icio.us

Digg

Reddit

Technorati

This article has been cited by other articles:

				M. Wang, J. You, K. G. Bemis, T. J. Tegeler, and D. P. G. Brown Label-free mass spectrometry-based protein quantification technologies in proteomic analysis Brief Funct Genomic Proteomic, September 1, 2008; 7(5): 329 - 339. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) Supplementl Data All Versions of this Article: M700440-MCP200v1 7/5/927    most recent Submit a response Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Request Permissions Glossary Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Colzani, M. Articles by Quadroni, M. Search for Related Content PubMed PubMed Citation Articles by Colzani, M. Articles by Quadroni, M. Social Bookmarking                      What's this?