{'fraction_non_alphanumeric': 0.033891758798821996, 'fraction_numerical': 0.021966880702940183, 'mean_word_length': 5.250301750150875, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Genome-scale Proteome Quantification by DEEP SEQ Mass Spectrometry

Advances in chemistry and massively parallel detection underlie DNA sequencing platforms that are poised for application in personalized medicine. In stark contrast, systematic generation of protein-level data lags well-behind genomics in virtually every aspect: depth of coverage, throughput, ease of sample preparation, and experimental time. Here, to bridge this gap, we develop an approach based on simple detergent lysis and single-enzyme digest, extreme, orthogonal separation of peptides, and true nanoflow LC-MS/MS that provides high peak capacity and ionization efficiency. This automated, deep efficient peptide sequencing and quantification (DEEP SEQ) mass spectrometry platform provides genome-scale proteome coverage equivalent to RNA-seq ribosomal profiling and accurate quantification for multiplexed isotope labels. In a model of the embryonic to epiblast transition in murine stem cells, we unambiguously quantify 11,352 gene products that span 70% of Swiss-Prot and capture protein regulation across the full detectable range of high-throughput gene expression and protein translation.

The field of DNA and RNA sequencing has progressed at a remarkable rate since release of the initial human genome draft sequences 1, 2 . Genotype-phenotype association studies continue to elucidate large genomic regions that are amplified, deleted, or otherwise altered in the context of human disease, although identification of specific causal gene elements has proven to be a significant challenge 3 . Despite these and other successes 4 , primary DNA sequence data does not capture downstream regulatory information for protein translation, degradation, and post-translational modification status 5,6,7 . These data are critical components of studies designed to quantitatively monitor biological response to perturbation or build predictive models of cellular physiology 8 . As a result, there remains a clear and unmet need for methodologies that can provide systematic and scalable sequence characterization of proteins as an important functional complement to DNA and RNA data.
The wide dynamic range of protein expression and vast array of post-translational modifications, coupled with the lack of an analyte amplification strategy analogous to PCR, presents tremendous challenges for genome-wide protein characterization, particularly for signal transduction and other key regulatory factors that are often present in low abundance. As a result, the majority of shotgun sequencing approaches frequently rely on low throughput pre-fractionation (e.g., prior to LC-MS/MS analysis) of subcellular compartments or intact proteins to improve dynamic range. While these techniques are subject to continued improvements, they have not yet achieved genome-scale protein identification and quantification. Moreover in many cases the labor-intensive nature of prefractionation protocols hinders widespread adoption and standardization.
In this work we forgo cellular-or protein-level pre-fractionation altogether and instead rely on direct detergent-based protein solubilization, followed by single-enzyme trypsin digest and extensive, fully automated temporal separation of peptides through multiple physicochemically orthogonal stages: high-pH reversed phase (RP) and strong anion exchange (SAX) dimensions coupled in series and directly with a narrow-bore, extended length (25 μm × 100 cm) low-pH RP analytical column operated in a true nanoflow regime. The latter chromatographic stage provides high peak capacity separation in a third orthogonal dimension, along with increased electrospray ionization efficiency for improved detection. The figures of merit for these individual components combine to yield an automated deep efficient peptide sequencing and quantification (DEEP SEQ) mass spectrometry platform that provides unprecedented separation capacity, rapid sequencing speed, and quantification of proteins across the entire range of mammalian gene expression and protein translation. We utilized multiplexed stable isotope labels in the context of a model designed to profile early perturbations in the naïve, self-renewing ground state of murine embryonic stem cells to quantify 211,535 unique peptide sequences that mapped unambiguously to 11,352 gene products. These results span ∼70% of the highly curated Swiss-Prot database, capture a vast majority of known pluripotent factors, and provide a depth and scale of proteome coverage commensurate with genome-wide analysis of protein translation by RNA-seq ribosomal profiling 9 .

Peptide quantification by DEEP SEQ mass spectrometry
Multiplexed isotope labels have emerged as an enabling technology to increase the throughput of quantitative proteomic studies and provide increased flexibility with respect to experimental variables (time course, dose response, etc.) 10 . As these reagents have enjoyed wider use, it has become apparent that quantification accuracy degrades significantly as the complexity of target analyte mixtures increases. Simultaneous fragmentation of precursors that overlap in m/z leads to a compression of discrete peptide ratios towards the mean of all measured values 11,12,13 . We reasoned that the extreme temporal separation provided on our protein deep sequencing platform would mitigate deleterious ratio compression effects. To explore this hypothesis we created a mixed-species quantification model 12,13 in which the contaminant peptides were present at a mass ratio of ∼6and ∼32-fold greater as compared to the target species, respectively (Fig. 1a). We acquired data at a depth of 20 DEEP SEQ fractions (Methods) with an LTQ-Orbitrap XL and used multiplierz 14,15 to compile and analyze extracted ion chromatograms for all detected peptides. We observed an average chromatographic peak width measured at half-height of 27.8 seconds (Fig. 1b), with over 97% of all detected peptides spanning no more than two adjacent fractions (Fig. 1c), yielding an empirical peak capacity 16 of ∼1.3E4. Importantly, projection of peptide sequence identifications as a function of first dimension organic and second dimension salt concentrations (Fig. 1d) revealed that peptides were distributed throughout the entire separation space, confirming that our platform provided orthogonal and extensive peptide separation in these experiments. We next plotted iTRAQ ratios for species-specific peptides and compared these to ratios observed in a conventional, single-dimension LC-MS/MS analysis of the same sample mixture. The median ratio for the mixed-species channels (116:114) in the standard LC-MS/MS analysis was compressed by some 40% as compared to the single-species ratios (117:115); however, the equivalent analysis with DEEP SEQ mass spectrometry revealed a relative ratio compression of only 8% ( Fig. 1e and Supplementary Data 1). Inspection of individual peptide spectra (Fig. 1f) reveals the marked improvement in quantification data provided by our DEEP SEQ platform.

Genome-wide protein quantification comparable to ribo-seq
To further explore the performance of our deep protein sequencing platform we sought to quantify changes in mESC protein expression resulting from withdrawal of leukemia inhibitory factor (LIF), a cytokine required to maintain self-renewal in these cells. After iTRAQ labeling we acquired MS/MS data at a depth of 20 DEEP SEQ fractions (Methods) on a Triple TOF 5600 mass spectrometer (Fig. 2a). A single experiment which required approximately 1 day for sample preparation and another 8 days for data acquisition yielded nearly 2 million MS/MS spectra and 0.6 million peptide spectral matches (PSMs) based on a search against the Uniprot database which is composed of Swiss-Prot (16,502 mouse genes) plus Trembl (9,769 mouse genes) and contains 26,271 mouse genes in total. These PSMs corresponded to180,867 peptides (≤1%FDR, including chemical modifications). Following the convention described by Qeli et al. 17 we mapped the set of 128,513 non-redundant peptide sequences into the genome as follows: First, 107,722 peptides were uniquely assignable to 9,818 gene IDs; we refer to these as "Class I" peptides and genes, respectively ( Fig. 2b). Of the remaining 20,791 peptides, 2,819 sequences mapped to at least two of another 1,737 gene IDs ("Class II" peptides and genes, respectively, Supplementary Data 2); importantly Class II peptides do not map to Class I genes. The remaining 17,972 peptides ("Class III" peptides) map to multiple genes, including one or more Class I genes. Class III peptides are included in the count of total peptide detections but are not otherwise considered for purposes of protein identification. Cumulatively across three biological replicate experiments, comprising 24 days of data acquisition, we obtained ∼5.9 million MS/MS spectra, ∼1.8 million PSMs, and identified 211,535 non-redundant peptide sequences (178,167 Class I, 3,496 Class II, and 29,872 Class III peptides, respectively). This peptide set encompassed 13,075 and 2,824 Class I and II genes, respectively (Fig. 2c), with the former (Class I genes) spanning 50% of UniProt and nearly 70% of the manually curated, non-redundant Swiss-Prot mouse database (Fig. 2d,e).
Two recent reports 18,19 that sought to analyze global protein expression in mESC utilized pre-fractionation of sub-cellular compartments, proteins, or peptides (or some combination) prior to LC-MS/MS. We took the union of published data from these studies and found that results from our DEEP SEQ platform encompass ∼73% and ≥95% of the reported peptide and gene IDs, respectively, while adding a significant quantity of new sequence information ( Fig. 3a and Supplementary Data 3). Although the fractionation approaches, mass spectrometry instrumentation, and database search engines varied across these three studies we observed very similar physicochemical properties for peptides identified in each data set (Fig. 3b,c). To explore the dynamic range of proteins identified across these data we next overlaid Class I and II genes derived from each peptide set with mRNA expression data for mESC 20 (Fig. 3d and Supplementary Data 4). Class I genes identified in our data spanned the full dynamic range of gene expression, with new Class I peptides identified through additional replicate experiments generally mapping to genes expressed at lower levels. Moreover, peptides identified by the two previous studies 18,19 were biased to highexpression genes as compared to those in our data. Strikingly the addition of Class II genes from our DEEP SEQ analysis improved overall proteome coverage only incrementally as compared to the set of Class I genes alone, with a small bias towards low-expression genes. Given the limitations of microarray data as a surrogate for protein expression, we next sought to compare our proteomic data with that from ribosomal profiling by RNA-seq, a technique that monitors protein translation on a genome-wide scale 21 . In a recent ribosomal profiling study, Ingolia et al. 9 measured translation of 12,674 transcripts from 19,022 protein coding genes in mESC. Our set of Class I and II genes spanned an equivalent fraction of the total gene space (63%) and encompassed 81% of all translation events represented in the ribosomal profiling data (Fig 3e and Supplementary Data 5). Again we observed that our deep sequencing data spanned the entire dynamic range of protein translation as represented by ribosomal profiling, including significant coverage (∼42%) of low-frequency translation events (Fig. 3f). As was observed with the comparison to gene expression by microarray, our data exhibited significantly higher coverage of low-expression genes as compared to that from previous studies of the mESC proteome 18,19 , while inclusion of DEEP SEQ Class II genes provided only a negligible improvement in proteome coverage. The latter observation is particularly notable given that Class II genes are defined by shared peptides; as a result the data for these gene products is ambiguous with respect to identification and quantification. For these reasons we focused the remainder of our analyses exclusively on Class I genes. Collectively these results suggest that our protein deep sequencing platform provides accurate quantitative data for multiplexed stable isotope reagents while simultaneously maximizing the discovery potential in the analysis of complex mammalian proteomes. In fact, using very stringent criteria for unambiguous protein identification (Class I peptides) our DEEP SEQ mass spectrometry approach provides data that span the full dynamic range of mammalian protein expression.

Functional proteome coverage by DEEP SEQ mass spectrometry
Regulatory proteins such as transcription factors, kinases and other signal transduction factors are often underrepresented in mass spectrometry-based whole-proteome studies. Hence we next sought to evaluate the coverage of functional protein classes provided by our DEEP SEQ platform. We found that the set of Class I peptides mapped unambiguously to ≥40% of genes across Gene Ontology (GO) functional categories, regardless of the evidence filter used (Fig. 4a). We also observed that the coverage varied across categories, from nearly 100% for genes encoding ribosome and proteasome proteins to approximately 20% for transmembrane receptors; the same trend was observed for data derived from ribosomal profiling, suggestive of a general correlation between coverage and protein abundance ( Fig  4b). This hypothesis was further corroborated by mapping functional protein classes to gene expression in mESC (Fig. 4c). Consistent with these data we found that replicate experiments tended to augment representation of proteins encoded by genes expressed at low levels. Overall, our set of stringently defined Class I genes spanned 52% of the functional proteome, on par with coverage provided by ribosomal profiling (56% coverage); importantly, our DEEP SEQ mass spectrometry data captured ≥70% of key regulatory protein classes, including kinases, phosphatases, and transcription factors.

Quantification of LIF-dependent functional proteome in mESC
To explore the functional proteome characteristic of mESC, we next compiled two reference sets of pluripotent factors derived from genetic depletion screens 22,23 and proteins biochemically associated with the master regulatory transcription factors Oct4 and Nanog 24, 25, 26 (Supplementary Data 6). The set of Class I genes from our biological triplicate analysis encompassed ∼81% (Fig. 5a) and ∼90% (Fig. 5b) of these reference gene sets, respectively, suggesting that our DEEP SEQ mass spectrometry platform can capture regulatory events associated with key mediators of pluripotency in mESC. A plot of quantitative data (Fig. 5c) revealed numerous proteins (Table 1) whose expression level reproducibly increased or decreased in response to LIF withdrawal (Supplementary Data 7). Query of these proteins against the functional and biochemical reference sets demonstrated an enrichment (Fishers exact test, Pval = 2.1E-5) for pluripotent factors (Fig. 5d). Biochemical validation (Fig. 5e) of LIF-mediated expression for a sub-set of proteins confirms the power of DEEP SEQ mass spectrometry to provide accurate quantification for multiplexed isobaric labels in the context of genome-wide proteome profiling.

Discussion
Continued advances in our understanding of genome variation 2 , gene expression 27,28 , and translation 9 , driven in-part by the advent of next-generation DNA/RNA sequencing technologies 3 , has effectively re-kindled interest in experiments designed to maximize discovery potential and simultaneously quantify known genomic or other biomolecular events. Establishing a parallel trajectory for this paradigm in proteomics is complicated by the broad range of protein abundance, diverse repertoire and stoichiometry of posttranslational modifications, as well as the stochastic nature of discovery-driven or shotgun MS/MS acquisition methods. In fact today, nearly two decades after Marc Wilkens first coined the term, "proteome" 29 , the dynamic range of protein expression in mammalian systems has represented an insurmountable hurdle for genome-wide proteome quantification. Despite these obstacles, the functional content of protein-level data represents an important complement to genomic-based studies, and hence there remains significant motivation to develop scalable platforms for proteomic analyses.
In order to achieve true, genome-scale sequence coverage along with high-fidelity protein quantification, we sought to significantly improve the analytical figures of merit for each component of our DEEP SEQ mass spectrometry platform. First the use of physicochemically orthogonal high-pH reversed phase and strong anion exchange separation stages 16,30 coupled with a narrow-bore (25 μm I.D.), extended length (100 cm) low-pH analytical column 31 provides extreme temporal separation of peptides with an empirical peak capacity of ∼1.3E4. Our experience with online multidimensional separations 16,30 indicates that further optimization of first and second dimension eluent concentrations can yield an improved distribution of peptides across each separation dimension, ultimately providing deeper proteome coverage beyond the ∼70% achieved herein. Second, the final dimension column geometry maintains the integrity of chromatographic separation at ultralow effluent flow rates (∼5nL/min), thus maximizing electrospray ionization efficiency 32,33 . Third, all separation stages in the DEEP SEQ configuration are implemented in microcapillary format and coupled in series, with the final dimension interfaced directly to the mass spectrometer, providing for fully automated operation, along with efficient capture and transfer of peptides across all separation stages.
Based on a stringent peptide-to-gene I.D. criterion, our DEEP SEQ mass spectrometry analysis of mESC provided quantitative expression data for 11,352 out of 16,502 genes in SwissProt (∼69%, Fig. 2e), a depth of coverage equivalent to a recent study in which RNAseq ribosomal profiling was used to measure translation of 12,674 out of 19,022 (∼66%) protein coding genes contained in the UCSC mouse database 9 . Importantly our data provide significant coverage (∼70-85%) for key regulatory protein families, including kinases, ubiquitin ligases, and transcription factors (Fig. 4b). In fact, proteins quantified on our deep sequencing platform span the full dynamic range of corresponding gene expression and protein translation profiles (Fig. 3d-f). Not surprisingly, improvements in proteome coverage across triplicate experiments were most pronounced for low-expression protein families (Fig. 4c).
The absence of protein amplification strategies analogous to PCR, combined with the limited improvements in dynamic range offered by each generation of mass spectrometry hardware, places a disproportionate burden on separation techniques to achieve robust detection of low-abundance proteins. Despite numerous methodological studies, the choice of separation strategies that will provide the best combination of sample yield, experiment time, and ultimately, proteome coverage is unresolved. Our results refute the notion that analysis of tryptic peptides alone is insufficient to characterize proteins across a wide dynamic range of abundance. In practice, DEEP SEQ mass spectrometry achieves significantly higher proteome coverage (Fig. 3a) from only a fraction (<5%) of the input required by techniques that rely on protein-or cellular-level pre-fractionation 18,19 . In fact, while chromatographic and sub-cellular fractionation are slow and low-throughput compared to chemical amplification, the time required for our DEEP SEQ analysis (25 days for biological triplicates, including sample preparation) is not dramatically different from other recent attempts at deep protein sequencing by mass spectrometry (for example, 21 days for a study that relied on the use of protein level fractionation and multiple enzymes 34 ), or for that matter genome-wide RNA-seq ribosomal profiling (9-12 days as described in a recent reveiw 35 ). Finally it is worth noting that our DEEP SEQ approach provides the serendipitous benefit of a robust and streamlined sample preparation protocol. Simple detergent solubilization, single enzyme digestion, peptide desalting, and iTRAQ labeling is somewhat reminiscent of the commoditized kits used in conjunction with next-generation DNA/RNA sequencing.
The quality of peptide fractionation provided by our deep protein sequencing platform is also evident in the analysis of a mixed-species iTRAQ quantification model (Fig. 1a). Even with contaminant species present in >30-fold excess relative to the target peptides, our DEEP SEQ mass spectrometry platform provided sufficient separation peak capacity to yield accurate iTRAQ ratios (Fig. 1e). These data are important in light of discrepancies reported between the throughput advantages afforded by multiplexed reagents and the welldocumented limitations encountered when using these labels to quantify proteins in complex mixtures 11,12,13,36 . In fact, results from three of these studies 12,13,36 demonstrated unequivocally that two dimensions of peptide fractionation are insufficient to abrogate suppression of multiplexed ratios for analysis of samples intended to represent complex proteomes. Importantly, DEEP SEQ mass spectrometry provides a platform to reconcile these juxtaposed observations and fully leverage the advantages of isobaric reagents in studies designed to quantify proteome response to perturbation on a genome-wide scale. Moreover, recent reports 37, 38 suggest the possibility of significant increases in the degree of multiplexing for these stable isotope reagents. These advances, along with improvements in mass spectrometry instrumentation 39,40 in addition to new chromatographic 41,42,43 and electrophoretic 44,45,46 separation platforms, will yield a concomitant increase in the throughput of DEEP SEQ mass spectrometry analysis.
Finally, our data provide insight into the transition of mESC from a naïve ground state, similar to that observed in the embryo inner cell mass prior to implantation, to a postimplantation epiblast-like stage (mEpiSC) that is poised for directed differentiation 47 . LIF is a critical cytokine that supports self-renewal in mESC through the Jak-Stat and Pi3k-Akt signaling axes, while opposing pathways, including Wnt-Gsk3β and Mek-Erk, mediate transcription programs that degrade pluripotent potential 48 . How these and other exogenous stimuli coordinately influence the core pluripotent genes (Oct4 and Sox2) and other peripheral factors to enforce the transcriptional ground state or mediate lineage commitment is not yet fully resolved.
Strikingly we detect regulated protein expression for downstream targets linked to each of these discrete pathways (Klf4, Klf5 49 , Lef1 50 , Esrrb, Tcfcp2l1 51 and Tbx3 48 ). Proteins having roles in other developmental contexts (Otx2 52 , Pou3f1/Oct6 53 , and CD9 54, 55 ) along with epigenetic factors (Dnmt3a/b) that are critical for high-fidelity gene expression 56 , were also regulated in a LIF-dependent manner. In addition to these studies that targeted specific genes or pathways, our DEEP SEQ analysis also quantified the majority of pluripotent factors previously defined by high-throughput functional genetic 22,23 and biochemical interaction 24, 25, 26 assays (Fig. 5a,b). Indeed, we reproducibly observed regulated expression for 50 proteins ( Table 1) and found that these were enriched for pluripotent and developmental genes (Fig. 5d). Importantly, this set of putative LIF-dependent protein targets spans the full dynamic range of gene expression and protein translation in mESC ( Fig. 5f,g). Thus, though results presented herein represent a single time point in the LIFmediated transition between naïve and primed epiblast states, the depth and scale of these data provide compelling evidence that our DEEP SEQ mass spectrometry platform can capture the vast majority of the mESC functional proteome in the context of more complex experiments 57, 58 designed to decipher individual contributions of the above pathways to self-renewal and lineage commitment. More generally our DEEP SEQ mass spectrometry platform represents a powerful and scalable approach for genome-wide profiling of protein expression and post-translational modification status in mammalian systems.
Mouse embryonic stem cell (mESC) line J1 was generously provided by Dr. Stuart Orkin (Dana-Farber Cancer Institute and Children's Hospital, Boston, MA). Initially, 10cm Nunclon tissue culture dishes (ThermoFisher Scientific, Waltham, MA) were coated with 0.1% Gelatin (EMD Millipore, Billerica, MA) at room temperature for 30 minutes. Gelatin was aspirated and J1 mESC were plated at a density of ∼6 × 10 4 /cm 2 in Dulbecco's modified Eagle's medium with high glucose supplemented with L-glutamine (Gibco Life Technologies, Carlsbad, CA), 15% embryonic stem cell validated FCS (Stem Cell Technologies, Vancouver BC, Canada), 2-mercaptoethanol, nucleosides (EMD Millipore, Billerica, MA), nonessential amino acids (Gibco Life Technologies, Carlsbad, CA), and 10 ng/mL murine LIF6 (EMD Millipore, Billerica, MA). Perturbation experiments were performed by establishing J1 mESC under the above conditions and then removing LIF6 for 48 hours. At the time of harvest, plates were washed with cold PBS to remove serum proteins and adherent cells were lysed by the addition of a boiling SDS (50mM Tris-HCl, pH 7.5, 7.5% SDS, 5% glycerol, 50mM DTT, 5mM EDTA). Lysed cells were centrifuged and the supernatant extract was stored at −80 °C.

Sample preparation for DEEP SEQ analysis
Proteins were precipitated by adding six volumes of cold (−20°C) acetone and resolubilized in digestion buffer containing 8 M urea and 0.1 M NH 4 HCO 3 . Total protein levels were measured by BCA. Dithiothreitol (DTT) was added to a final concentration of 10 mM and incubated for 30 minutes at 60 °C, followed by addition of methyl methanethiosulfonate (MMTS) (ThermoFisher Scientific, Waltham, MA) to 20 mM. After 30 min. incubation in the dark at room temperature, excess MMTS was quenched by addition of DTT to a final concentration of 20 mM. Reduced and alkylated proteins were diluted in 0.1M ammonium bicarbonate, followed by addition of trypsin, with overnight digestion 37 °C and end-overend rotation. Digested sample solutions were loaded onto SepPak C18 reverse phase cartridges (Waters Corp., Milford, MA) to remove urea and other salts. Eluted peptides (45% acetonitrile in water with 0.1% trifluoroacetic acid) were lyophilized by vacuum centrifugation.
Peptides were labeled with 4-plex iTRAQ reagents (AB Sciex, Framingham, MA). Two aliquots of 0.2 μg yeast peptides were labeled with 114 and 115. Separately, two other aliquots of 1.0 μg yeast peptides were labeled with 116 and 117. In addition, two aliquots of 6.3 μg peptides from mESC were labeled with 114 and 116. For the mESC perturbation experiment, peptides from the cells incubated without LIF were labeled as technical replicates with iTRAQ116 and 117, while peptides from the cells incubated with LIF were labeled as technical replicates with iTRAQ 114 and 115. For each reaction, peptides were resuspended in 500mM triethylammonium bicarbonate and mixed with the appropriate iTRAQ reagent in ethanol. Labeling was allowed to proceed at room temperature for one hour; samples were then combined and dried by vacuum centrifugation. Fresh aliquots of mESC were processed as described above to provide biological triplicates.

DEEP SEQ multi-dimension separation
Three dimension peptide separation was performed on a modified Waters (Milford, MA) NanoAcquity UHPLC system with binary and isocratic pumps, along with an autosampler and additional 6-port  (Supplementary Data 8). The binary pump delivered 0.1% formic acid at 8μL/min. to dilute the organic content and acidify the 1 st /2 nd dimension effluent prior to the 3 rd dimension pre-column, or provided for gradient elution (2-50% B in 580 minutes, A = 0.1% formic acid, B = acetonitrile with 0.1% formic acid) of peptides from the 3 rd dimension reversed phase columns for LC-MS/MS analysis at a flow rate of ∼5nL/min. A Digital PicoView electrospray source platform (New Objective, Woburn, MA) was used on both the Orbitrap XL and 5600 Triple TOF mass spectrometers to automatically position the emitter tip at the source inlet during LC-MS/MS acquisition or beneath a gravity-driven drip station during injection of peptide samples or 1 st /2 nd dimension eluents.

DEEP SEQ sample capacity and experiment time
Based on our published and unpublished studies, we estimate that the total loading capacity of our first dimension RP column (200 μm I.D. × 20 cm) is currently ∼100 μg. The DEEP SEQ platform is easily tailored to a variety of sample types. Generally, first and second dimension eluent concentrations in a range of 7 to 55% acetonitrile and 10 to 300 mM KCl represent the useful boundaries for peptide elution. Importantly, our experience to date suggests that these conditions are robust with respect to biological input, obviating the need to run repeated pilot experiments for every sample. Total time of analysis is another important consideration when using multidimensional fractionation techniques. For example, DEEP SEQ experiments performed at a depth of 20 fractions required 8 days for data acquisition. System reliability is particularly important given that a single, DEEP SEQ mass spectrometry analysis will typically require several days of continuous instrument time. Importantly, all data presented herein were acquired using a single 25μm I.D. ×100 cm resolving column, demonstrating the robustness of our platform.

DEEP SEQ mass spectrometry data acquisition parameters
The LTQ-Orbitrap XL mass spectrometer (Thermo, Waltham, MA) was operated in data dependent mode, such that the top 10 most abundant precursors in each MS scan were subjected to MS/MS in both CAD and HCD mode (CAD in the linear trap, normalized collision energy = 35%, precursor isolation width = 1.9 Da, intensity threshold for precursor selection = 20,000 HCD in the orbitrap, normalized collision energy = 47%, precursor isolation width = 1.0 Da, intensity threshold for precursor selection = 20,000). Dynamic exclusion was enabled with a repeat count of 1 and exclusion duration set to 20 seconds. Electrospray voltage was 2.2 kV. We enabled the Lock Mass feature and selected m/z = 445.120025 ([Si(CH 3 ) 2 O] 6 ) as the internal calibrant.
The 5600 Triple TOF (AB Sciex, Framingham, MA) was operated in information dependent mode (IDA), with the top 50 precursors (charge state +2 to +5, >100 counts) in each MS scan (800 ms, scan range 350-1500 m/z) subjected to MS/MS (minimum time 140 ms, scan range 100-1400m/z). A dynamic exclusion window of 20 s was used with unit resolution for precursor isolation. Electrospray voltage was 2.2 kV.

DEEP SEQ mass spectrometry data processing
Our API-based multiplierz software framework was used to extract and format MS/MS data from the Orbitrap XL for subsequent search against the Uniprot mouse database (downloaded on 11/02/2011) and a S. cerevisiae database (downloaded from http:// www.yeastgenome.org/01/06/2010). MS/MS data was searched using Protein Pilot V4.4 (AB Sciex, Framingham, MA) with the following parameters: instrument ="Orbi/FT MS (sub-ppm), LTQ MS/MS" for data acquired with the Orbitrap XL, or "5600 TripleTOF" for data acquired with the 5600 Triple TOF mass spectrometer. A fixed modification of +42 corresponding to MMTS of cysteine was also included. The sample type was set to "iTRAQ 4-plex (peptides labeling)". All peptide spectral matches (PSM) from biotriplicate acquisitions were combined for the FDR assessment. Only those peptides with scores at or above a PSM FDR threshold of 1% were further considered. All peptides were mapped back to the genome without consideration of splice isoforms from the same gene. Following the procedure of Qeli et al. 17 , peptides passing the 1% FDR filter were classified as either "unique" or "shared" based on whether they could be assigned to only one or more than one gene, respectively. Class I genes were defined as those identified only by uniquely assignable peptides. Class II genes were identified based on shared peptides that could not be assigned to any Class I gene. The results reported herein are based on Class I and Class II genes; no Class III genes were included in our analyses. Multiplierz scripts were used to systematically extract XICs and calculate corresponding peak widths for all identified peptides. To estimate the orthogonality of peptide fractionation, the number of unique peptides identified in each third dimension LC-MS/MS run was represented as a circle of proportional diameter and projected onto a 2D plot at the corresponding 1st dimension acetonitrile (x-axis) and 2nd dimension salt (y-axis) concentrations used in the experiment. Linear regression was performed using the function stats. linregress from Numpy version 1.4.1, with each unique peptide sequence used as a separate data point.
In the mixed species quantification experiments the iTRAQ 116 /iTRAQ 114 and iTRAQ 117 / iTRAQ 115 ratios were first normalized based on all spectra identified (final normalization factor was 1.14). We used R scripts to create box-plots for all iTRAQ ratios. For the quantitative analyses of mESC, we summed iTRAQ channels in each biological condition (±LIF) for Class I peptides to generate a final ratio for each Class I gene. Across three biological triplicate experiments a protein was considered to be regulated in expression (Table 1) if the following criteria were met in two out of the three experiments: (i) the iTRAQ intensities exceeded 200 counts and (ii) the Log 2 iTRAQ ratio was ≥1 or ≤ −1. Enrichment of pluripotent factors within the set of genes whose expression was reproducibly regulated across biotriplicate experiments was estimated using a Fisher's exact test. A null distribution and two-sided p-value were calculated based on code freely available at: (http:// research.microsoft.com/en-us/um/redmond/projects/MSCompBio).

Comparisons to mESC microarray and ribosomal profiling data
Comparison of DEEP SEQ mass spectrometry data to that from ribosome profiling for mESC was performed as follows: The UCSC mouse genome database (mm9) (http:// hgdownload.soe.ucsc.edu/goldenPath/mm9/database/) was downloaded on 10/02/2012. We aligned the peptides identified in triplicate DEEP SEQ mass spectrometry analyses against sequences in the UCSC database. A similar alignment was performed for mESC ribosomal profiling data previously reported 9 . In 1,034 cases we were unable to reconcile gene names from the UCSC and UniProt databases; these genes were removed from further consideration. Average gene expression levels across three wild type mESC lines (B21, J1 and R1) was calculated using published datasets (accession codes GSM338369, GSM338371, GSM338373) 20 . We used the R/Bioconductor software environment to renormalize all microarray data based on the multi-array averaging method (RMA) and to remap all probe sequences to Entrez Gene IDs using a custom Chip Definition File (CDF) from the Michigan Microarray Lab (version 13). The average mRNA expression level (in Log 2 space) between the three arrays was calculated and the Entrez Gene IDs were converted to UniProt gene symbols. UniProt entries mapping to more than one Entrez Gene ID were excluded. The final list contained 15,705 entries of (average Log 2 expression value) | (UniProt gene symbol) pairs.

Gene ontology and reference gene sets for pluripotency
Positive reference sets (PRS) of pluripotent genes were created from RNAi screening data downloaded from two previous studies 22,23 . Similarly, pluripotency reference genes based on biochemical interactors of Nanog and Pou5f1 (Oct4) were created from three previous studies 24

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.
Science, Baltimore, MD) for helpful discussions related to mESC ribosomal profiling data. In addition the authors thank Christie Hunter and Sean Seymour of AB Sciex for support and technical assistance with the 5600 Triple TOF mass spectrometer and Protein Pilot search software.  Proteins were solubilized and processed directly (e.g., no sub-cellular or protein-level fractionation) for iTRAQ labeling (two replicates per condition), followed by DEEP SEQ mass spectrometry analysis. (b) A single, 20-fraction DEEP SEQ mass spectrometry analysis provided nearly 2 million MS/MS spectra, ∼610,000 peptide spectral matches (PSMs, <1%FDR), and 180,867 peptides corresponding to 128,513 unique peptide sequences. High-confidence (<1%FDR) peptide sequences were classified based on whether they mapped uniquely to mouse gene I.D.s (Class I peptide, Class I gene) or were shared across two or more mouse genes outside the set of Class I genes (Class II peptides, Class II genes). These data required 8 days of continuous acquisition time and yielded 9,818 Class I gene I.D.s. (blue). (c) Across biological triplicates, DEEP SEQ mass spectrometry analysis identified 13,075 Class I genes derived from ∼5.9 million MS/MS spectra acquired over 24 days. These data spanned (d) 50% of all protein-coding genes in the Uni-Prot database (26,271 entries). (e) The set of Class I peptides mapped unambiguously to nearly 70% of protein-coding genes in the manually annotated, non-redundant Swiss-Prot database (11,352 out of 16,502 total entries). (a) High-confidence peptides based on the union of published data from two previous mESC proteomic studies 18,19 were mapped to Class I and II gene I.D.s as described (see Methods) and compared with data from this study. High-confidence peptides identified in all three studies exhibited similar physicochemical properties, including (b) molecular weight distribution and (c) number of amino acids per peptide. Comparison of (d) Class I and II gene I.D.s with normalized mRNA levels in mESC demonstrates that DEEP SEQ mass spectrometry data spans the full range of gene expression. (e) The fraction of genes represented in the UCSC mouse database (grey) captured by ribosomal profiling (yellow) and DEEP SEQ mass spectrometry (red). The set of Class I and II genes identified herein encompassed 81% of protein translation events as represented by RNA-seq ribosomal profiling. (f) Overlay of ribosomal profiling data with that from DEEP SEQ mass spectrometry analysis demonstrates that the set of Class I genes spans the full dynamic range of protein translation in mESC. (a) Relative coverage provided by DEEP SEQ mass spectrometry and ribosomal profiling as a function of gene-ontology (GO) category and evidence filter. As defined by GO, "Manual All" contains all curator-reviewed experimental and computational annotations while "Manual Experiment" includes only curator-reviewed annotations derived from direct assay (IDA), physical interaction (IPI), mutant phenotype (IMP), genetic interaction (IGI), and expression pattern (IEP) evidence categories. (b) Relative coverage of GO protein families for DEEP SEQ mass spectrometry and ribosomal profiling analyses. Protein membership in each family is derived from the annotations within "Manual Experiment." (c) Relative mRNA expression level plotted in box-plot format (black and white, y-axis, right) as a function of GO protein families. Boxes encompass the interquartile range with respective median values indicated with stripes; whiskers represent 1.5× the interquartile range with outliers shown as open circles.(blue, y-axis, left) Relative increase in GO protein family representation across bio-triplicate DEEP SEQ mass spectrometry analyses for Class I genes. The set of Class I gene I.D.s from bio-triplicate, 20-fraction DEEP SEQ mass spectrometry experiments encompassed (a) ∼80% of murine pluripotent factors as defined by systematic genetic depletion assays 22,23 and (b) ∼90% of the Nanog and Oct4 transcription factor network as defined through analysis of biochemical interactions 24,25,26 . (c) Scatter plot of iTRAQ log-intensity versus log-ratio for 13,075 class I gene I.D.s identified across biotriplicate, 20-fraction DEEP SEQ mass spectrometry experiments. Individual iTRAQ signals for each high-confidence Class I peptide were summed across technical and biological replicates; these aggregate ratios were then combined to provide protein-level ratios. (d) Fisher's Exact Test confirms that the set of regulated proteins (Table 1) identified by DEEP SEQ mass spectrometry analysis is enriched (two-sided Pval = 2E-5) for pluripotent factors as defined by systematic loss-of-function and biochemical interaction assays. (e) Selected pluripotent and developmental factors were probed by western blot in mESC at 24 and 48 hrs after withdrawal of LIF. Left-most column indicates molecular weight markers. Overlay of regulated Class I gene products as detected by DEEP SEQ mass spectrometry with (f) relative gene expression as measured by microarray and (g) relative protein translation as measured by RNA-seq ribosomal profiling.