Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection.

BACKGROUND: Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. METHODS: We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. RESULTS: We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. CONCLUSIONS: This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Mapping the proteogenomic landscape enables prediction of drug response in acute myeloid leukemia.

Acute myeloid leukemia is a poor-prognosis cancer commonly stratified by genetic aberrations, but these mutations are often heterogeneous and fail to consistently predict therapeutic response. Here, we combine transcriptomic, proteomic, and phosphoproteomic datasets with ex vivo drug sensitivity data to help understand the underlying pathophysiology of AML beyond mutations. We measure the proteome and phosphoproteome of 210 patients and combine them with genomic and transcriptomic measurements to identify four proteogenomic subtypes that complement existing genetic subtypes. We build a predictor to classify samples into subtypes and map them to a "landscape" that identifies specific drug response patterns. We then build a drug response prediction model to identify drugs that target distinct subtypes and validate our findings on cell lines representing various stages of quizartinib resistance. Our results show how multiomics data together with drug sensitivity data can inform therapy stratification and drug combinations in AML.

Mass Spectrometry-Based Proteogenomics: New Therapeutic Opportunities for Precision Medicine.

Proteogenomics refers to the integration of comprehensive genomic, transcriptomic, and proteomic measurements from the same samples with the goal of fully understanding the regulatory processes converting genotypes to phenotypes, often with an emphasis on gaining a deeper understanding of disease processes. Although specific genetic mutations have long been known to drive the development of multiple cancers, gene mutations alone do not always predict prognosis or response to targeted therapy. The benefit of proteogenomics research is that information obtained from proteins and their corresponding pathways provides insight into therapeutic targets that can complement genomic information by providing an additional dimension regarding the underlying mechanisms and pathophysiology of tumors. This review describes the novel insights into tumor biology and drug resistance derived from proteogenomic analysis while highlighting the clinical potential of proteogenomic observations and advances in technique and analysis tools.

Proteomic and phosphoproteomic measurements enhance ability to predict ex vivo drug response in AML.

Acute Myeloid Leukemia (AML) affects 20,000 patients in the US annually with a five-year survival rate of approximately 25%. One reason for the low survival rate is the high prevalence of clonal evolution that gives rise to heterogeneous sub-populations of leukemic cells with diverse mutation spectra, which eventually leads to disease relapse. This genetic heterogeneity drives the activation of complex signaling pathways that is reflected at the protein level. This diversity makes it difficult to treat AML with targeted therapy, requiring custom patient treatment protocols tailored to each individual's leukemia. Toward this end, the Beat AML research program prospectively collected genomic and transcriptomic data from over 1000 AML patients and carried out ex vivo drug sensitivity assays to identify genomic signatures that could predict patient-specific drug responses. However, there are inherent weaknesses in using only genetic and transcriptomic measurements as surrogates of drug response, particularly the absence of direct information about phosphorylation-mediated signal transduction. As a member of the Clinical Proteomic Tumor Analysis Consortium, we have extended the molecular characterization of this cohort by collecting proteomic and phosphoproteomic measurements from a subset of these patient samples (38 in total) to evaluate the hypothesis that proteomic signatures can improve the ability to predict response to 26 drugs in AML ex vivo samples. In this work we describe our systematic, multi-omic approach to evaluate proteomic signatures of drug response and compare protein levels to other markers of drug response such as mutational patterns. We explore the nuances of this approach using two drugs that target key pathways activated in AML: quizartinib (FLT3) and trametinib (Ras/MEK), and show how patient-derived signatures can be interpreted biologically and validated in cell lines. In conclusion, this pilot study demonstrates strong promise for proteomics-based patient stratification to assess drug sensitivity in AML.

The AML microenvironment catalyzes a stepwise evolution to gilteritinib resistance.

Our study details the stepwise evolution of gilteritinib resistance in FLT3-mutated acute myeloid leukemia (AML). Early resistance is mediated by the bone marrow microenvironment, which protects residual leukemia cells. Over time, leukemia cells evolve intrinsic mechanisms of resistance, or late resistance. We mechanistically define both early and late resistance by integrating whole-exome sequencing, CRISPR-Cas9, metabolomics, proteomics, and pharmacologic approaches. Early resistant cells undergo metabolic reprogramming, grow more slowly, and are dependent upon Aurora kinase B (AURKB). Late resistant cells are characterized by expansion of pre-existing NRAS mutant subclones and continued metabolic reprogramming. Our model closely mirrors the timing and mutations of AML patients treated with gilteritinib. Pharmacological inhibition of AURKB resensitizes both early resistant cell cultures and primary leukemia cells from gilteritinib-treated AML patients. These findings support a combinatorial strategy to target early resistant AML cells with AURKB inhibitors and gilteritinib before the expansion of pre-existing resistance mutations occurs.

Night shift schedule causes circadian dysregulation of DNA repair genes and elevated DNA damage in humans.

Circadian disruption has been identified as a risk factor for health disorders such as obesity, cardiovascular disease, and cancer. Although epidemiological studies suggest an increased risk of various cancers associated with circadian misalignment due to night shift work, the underlying mechanisms have yet to be elucidated. We sought to investigate the potential mechanistic role that circadian disruption of cancer hallmark pathway genes may play in the increased cancer risk in shift workers. In a controlled laboratory study, we investigated the circadian transcriptome of cancer hallmark pathway genes and associated biological pathways in circulating leukocytes obtained from healthy young adults during a 24-hour constant routine protocol following 3 days of simulated day shift or night shift. The simulated night shift schedule significantly altered the normal circadian rhythmicity of genes involved in cancer hallmark pathways. A DNA repair pathway showed significant enrichment of rhythmic genes following the simulated day shift schedule, but not following the simulated night shift schedule. In functional assessments, we demonstrated that there was an increased sensitivity to both endogenous and exogenous sources of DNA damage after exposure to simulated night shift. Our results suggest that circadian dysregulation of DNA repair may increase DNA damage and potentiate elevated cancer risk in night shift workers.

Proteogenomic Characterization of Ovarian HGSC Implicates Mitotic Kinases, Replication Stress in Observed Chromosomal Instability.

In the absence of a dominant driving mutation other than uniformly present TP53 mutations, deeper understanding of the biology driving ovarian high-grade serous cancer (HGSC) requires analysis at a functional level, including post-translational modifications. Comprehensive proteogenomic and phosphoproteomic characterization of 83 prospectively collected ovarian HGSC and appropriate normal precursor tissue samples (fallopian tube) under strict control of ischemia time reveals pathways that significantly differentiate between HGSC and relevant normal tissues in the context of homologous repair deficiency (HRD) status. In addition to confirming key features of HGSC from previous studies, including a potential survival-associated signature and histone acetylation as a marker of HRD, deep phosphoproteomics provides insights regarding the potential role of proliferation-induced replication stress in promoting the characteristic chromosomal instability of HGSC and suggests potential therapeutic targets for use in precision medicine trials.

Unified feature association networks through integration of transcriptomic and proteomic data.

High-throughput multi-omics studies and corresponding network analyses of multi-omic data have rapidly expanded their impact over the last 10 years. As biological features of different types (e.g. transcripts, proteins, metabolites) interact within cellular systems, the greatest amount of knowledge can be gained from networks that incorporate multiple types of -omic data. However, biological and technical sources of variation diminish the ability to detect cross-type associations, yielding networks dominated by communities comprised of nodes of the same type. We describe here network building methods that can maximize edges between nodes of different data types leading to integrated networks, networks that have a large number of edges that link nodes of different-omic types (transcripts, proteins, lipids etc). We systematically rank several network inference methods and demonstrate that, in many cases, using a random forest method, GENIE3, produces the most integrated networks. This increase in integration does not come at the cost of accuracy as GENIE3 produces networks of approximately the same quality as the other network inference methods tested here. Using GENIE3, we also infer networks representing antibody-mediated Dengue virus cell invasion and receptor-mediated Dengue virus invasion. A number of functional pathways showed centrality differences between the two networks including genes responding to both GM-CSF and IL-4, which had a higher centrality value in an antibody-mediated vs. receptor-mediated Dengue network. Because a biological system involves the interplay of many different types of molecules, incorporating multiple data types into networks will improve their use as models of biological systems. The methods explored here are some of the first to specifically highlight and address the challenges associated with how such multi-omic networks can be assembled and how the greatest number of interactions can be inferred from different data types. The resulting networks can lead to the discovery of new host response patterns and interactions during viral infection, generate new hypotheses of pathogenic mechanisms and confirm mechanisms of disease.

An Integrative Analysis of Tumor Proteomic and Phosphoproteomic Profiles to Examine the Relationships Between Kinase Activity and Phosphorylation.

Phosphorylation of proteins is a key way cells regulate function, both at the individual protein level and at the level of signaling pathways. Kinases are responsible for phosphorylation of substrates, generally on serine, threonine, or tyrosine residues. Though particular sequence patterns can be identified that dictate whether a residue will be phosphorylated by a specific kinase, these patterns are not highly predictive of phosphorylation. The availability of large scale proteomic and phosphoproteomic data sets generated using mass-spectrometry-based approaches provides an opportunity to study the important relationship between kinase activity, substrate specificity, and phosphorylation. In this study, we analyze relationships between protein abundance and phosphopeptide abundance across more than 150 tumor samples and show that phosphorylation at specific phosphosites is not well correlated with overall kinase abundance. However, individual kinases show a clear and statistically significant difference in correlation among known phosphosite targets for that kinase and randomly selected phosphosites. We further investigate relationships between phosphorylation of known activating or inhibitory sites on kinases and phosphorylation of their target phosphosites. Combined with motif-based analysis, this approach can predict novel kinase targets and show which subsets of a kinase's target repertoire are specifically active in one condition versus another.

Phosphoproteome Analysis Reveals Estrogen-ER Pathway as a Modulator of mTOR Activity Via DEPTOR.

ER-positive breast tumors represent ∼70% of all breast cancer cases. Although their treatment with endocrine therapies is effective in the adjuvant or recurrent settings, the development of resistance compromises their effectiveness. The binding of estrogen to ERα, a transcription factor, triggers the regulation of the target genes (genomic pathway). Additionally, a cytoplasmic fraction of estrogen-bound ERα activates oncogenic signaling pathways such as PI3K/AKT/mTOR (nongenomic pathway). The upregulation of the estrogenic and the PI3K/AKT/mTOR signaling pathways are frequently associated with a poor outcome. To better characterize the connection between these two pathways, we performed a phosphoproteome analysis of ER-positive MCF7 breast cancer cells treated with estrogen or estrogen and the mTORC1 inhibitor rapamycin. Many proteins were identified as estrogen-regulated mTORC1 targets and among them, DEPTOR was selected for further characterization. DEPTOR binds to mTOR and inhibits the kinase activity of both mTOR complexes mTORC1 and mTORC2, but mitogen-activated mTOR promotes phosphorylation-mediated DEPTOR degradation. Although estrogen enhances the phosphorylation of DEPTOR by mTORC1, DEPTOR levels increase in estrogen-stimulated cells. We demonstrated that DEPTOR accumulation is the result of estrogen-ERα-mediated transcriptional upregulation of DEPTOR expression. Consequently, the elevated levels of DEPTOR partially counterbalance the estrogen-induced activation of mTORC1 and mTORC2. These results underscore the critical role of estrogen-ERα as a modulator of the PI3K/AKT/mTOR signaling pathway in ER-positive breast cancer cells. Additionally, these studies provide evidence supporting the use of dual PI3K/mTOR or dual mTORC1/2 inhibitors in combination with endocrine therapies as a first-line treatment option for the patients with ER-positive advanced breast cancer.

Proteogenomic Analysis of Human Colon Cancer Reveals New Therapeutic Opportunities.

We performed the first proteogenomic study on a prospectively collected colon cancer cohort. Comparative proteomic and phosphoproteomic analysis of paired tumor and normal adjacent tissues produced a catalog of colon cancer-associated proteins and phosphosites, including known and putative new biomarkers, drug targets, and cancer/testis antigens. Proteogenomic integration not only prioritized genomically inferred targets, such as copy-number drivers and mutation-derived neoantigens, but also yielded novel findings. Phosphoproteomics data associated Rb phosphorylation with increased proliferation and decreased apoptosis in colon cancer, which explains why this classical tumor suppressor is amplified in colon tumors and suggests a rationale for targeting Rb phosphorylation in colon cancer. Proteomics identified an association between decreased CD8 T cell infiltration and increased glycolysis in microsatellite instability-high (MSI-H) tumors, suggesting glycolysis as a potential target to overcome the resistance of MSI-H tumors to immune checkpoint blockade. Proteogenomics presents new avenues for biological discoveries and therapeutic development.

Proteogenomic Analysis of Surgically Resected Lung Adenocarcinoma.

INTRODUCTION: Despite apparently complete surgical resection, approximately half of resected early-stage lung cancer patients relapse and die of their disease. Adjuvant chemotherapy reduces this risk by only 5% to 8%. Thus, there is a need for better identifying who benefits from adjuvant therapy, the drivers of relapse, and novel targets in this setting. METHODS: RNA sequencing and liquid chromatography/liquid chromatography-mass spectrometry proteomics data were generated from 51 surgically resected non-small cell lung tumors with known recurrence status. RESULTS: We present a rationale and framework for the incorporation of high-content RNA and protein measurements into integrative biomarkers and show the potential of this approach for predicting risk of recurrence in a group of lung adenocarcinomas. In addition, we characterize the relationship between mRNA and protein measurements in lung adenocarcinoma and show that it is outcome specific. CONCLUSIONS: Our results suggest that mRNA and protein data possess independent biological and clinical importance, which can be leveraged to create higher-powered expression biomarkers.

Co-expression of CD39 and CD103 identifies tumor-reactive CD8 T cells in human solid tumors.

Identifying tumor antigen-specific T cells from cancer patients has important implications for immunotherapy diagnostics and therapeutics. Here, we show that CD103+CD39+ tumor-infiltrating CD8 T cells (CD8 TIL) are enriched for tumor-reactive cells both in primary and metastatic tumors. This CD8 TIL subset is found across six different malignancies and displays an exhausted tissue-resident memory phenotype. CD103+CD39+ CD8 TILs have a distinct T-cell receptor (TCR) repertoire, with T-cell clones expanded in the tumor but present at low frequencies in the periphery. CD103+CD39+ CD8 TILs also efficiently kill autologous tumor cells in a MHC-class I-dependent manner. Finally, higher frequencies of CD103+CD39+ CD8 TILs in patients with head and neck cancer are associated with better overall survival. Our data thus describe an approach for detecting tumor-reactive CD8 TILs that will help define mechanisms of existing immunotherapy treatments, and may lead to future adoptive T-cell cancer therapies.

Inhibition of interleukin-1 receptor-associated kinase-1 is a therapeutic strategy for acute myeloid leukemia subtypes.

Interleukin-1 receptor-associated kinase 1 (IRAK1), an essential mediator of innate immunity and inflammatory responses, is constitutively active in multiple cancers. We evaluated the role of IRAK1 in acute myeloid leukemia (AML) and assessed the inhibitory activity of multikinase inhibitor pacritinib on IRAK1 in AML. We demonstrated that IRAK1 is overexpressed in AML and provides a survival signal to AML cells. Genetic knockdown of IRAK1 in primary AML samples and xenograft model showed a significant reduction in leukemia burden. Kinase profiling indicated pacritinib has potent inhibitory activity against IRAK1. Computational modeling combined with site-directed mutagenesis demonstrated high-affinity binding of pacritinib to the IRAK1 kinase domain. Pacritinib exposure reduced IRAK1 phosphorylation in AML cells. A higher percentage of primary AML samples showed robust sensitivity to pacritinib, which inhibits FLT3, JAK2, and IRAK1, relative to FLT3 inhibitor quizartinib or JAK1/2 inhibitor ruxolitinib, demonstrating the importance of IRAK1 inhibition. Pacritinib inhibited the growth of AML cells harboring a variety of genetic abnormalities not limited to FLT3 and JAK2. Pacritinib treatment reduced AML progenitors in vitro and the leukemia burden in AML xenograft model. Overall, IRAK1 contributes to the survival of leukemic cells, and the suppression of IRAK1 may be beneficial among heterogeneous AML subtypes.

Proteome Profiling Outperforms Transcriptome Profiling for Coexpression Based Gene Function Prediction.

Coexpression of mRNAs under multiple conditions is commonly used to infer cofunctionality of their gene products despite well-known limitations of this "guilt-by-association" (GBA) approach. Recent advancements in mass spectrometry-based proteomic technologies have enabled global expression profiling at the protein level; however, whether proteome profiling data can outperform transcriptome profiling data for coexpression based gene function prediction has not been systematically investigated. Here, we address this question by constructing and analyzing mRNA and protein coexpression networks for three cancer types with matched mRNA and protein profiling data from The Cancer Genome Atlas (TCGA) and the Clinical Proteomic Tumor Analysis Consortium (CPTAC). Our analyses revealed a marked difference in wiring between the mRNA and protein coexpression networks. Whereas protein coexpression was driven primarily by functional similarity between coexpressed genes, mRNA coexpression was driven by both cofunction and chromosomal colocalization of the genes. Functionally coherent mRNA modules were more likely to have their edges preserved in corresponding protein networks than functionally incoherent mRNA modules. Proteomic data strengthened the link between gene expression and function for at least 75% of Gene Ontology (GO) biological processes and 90% of KEGG pathways. A web application Gene2Net (http://cptac.gene2net.org) developed based on the three protein coexpression networks revealed novel gene-function relationships, such as linking ERBB2 (HER2) to lipid biosynthetic process in breast cancer, identifying PLG as a new gene involved in complement activation, and identifying AEBP1 as a new epithelial-mesenchymal transition (EMT) marker. Our results demonstrate that proteome profiling outperforms transcriptome profiling for coexpression based gene function prediction. Proteomics should be integrated if not preferred in gene function and human disease studies.

The effect of inhibition of PP1 and TNFα signaling on pathogenesis of SARS coronavirus.

BACKGROUND: The complex interplay between viral replication and host immune response during infection remains poorly understood. While many viruses are known to employ anti-immune strategies to facilitate their replication, highly pathogenic virus infections can also cause an excessive immune response that exacerbates, rather than reduces pathogenicity. To investigate this dichotomy in severe acute respiratory syndrome coronavirus (SARS-CoV), we developed a transcriptional network model of SARS-CoV infection in mice and used the model to prioritize candidate regulatory targets for further investigation. RESULTS: We validated our predictions in 18 different knockout (KO) mouse strains, showing that network topology provides significant predictive power to identify genes that are important for viral infection. We identified a novel player in the immune response to virus infection, Kepi, an inhibitory subunit of the protein phosphatase 1 (PP1) complex, which protects against SARS-CoV pathogenesis. We also found that receptors for the proinflammatory cytokine tumor necrosis factor alpha (TNFα) promote pathogenesis, presumably through excessive inflammation. CONCLUSIONS: The current study provides validation of network modeling approaches for identifying important players in virus infection pathogenesis, and a step forward in understanding the host response to an important infectious disease. The results presented here suggest the role of Kepi in the host response to SARS-CoV, as well as inflammatory activity driving pathogenesis through TNFα signaling in SARS-CoV infections. Though we have reported the utility of this approach in bacterial and cell culture studies previously, this is the first comprehensive study to confirm that network topology can be used to predict phenotypes in mice with experimental validation.

Conservation of protein abundance patterns reveals the regulatory architecture of the EGFR-MAPK pathway.

Various genetic mutations associated with cancer are known to alter cell signaling, but it is not clear whether they dysregulate signaling pathways by altering the abundance of pathway proteins. Using a combination of RNA sequencing and ultrasensitive targeted proteomics, we defined the primary components-16 core proteins and 10 feedback regulators-of the epidermal growth factor receptor (EGFR)-mitogen-activated protein kinase (MAPK) pathway in normal human mammary epithelial cells and then quantified their absolute abundance across a panel of normal and breast cancer cell lines as well as fibroblasts. We found that core pathway proteins were present at very similar concentrations across all cell types, with a variance similar to that of proteins previously shown to display conserved abundances across species. In contrast, EGFR and transcriptionally controlled feedback regulators were present at highly variable concentrations. The absolute abundance of most core proteins was between 50,000 and 70,000 copies per cell, but the adaptors SOS1, SOS2, and GAB1 were found at far lower amounts (2000 to 5000 copies per cell). MAPK signaling showed saturation in all cells between 3000 and 10,000 occupied EGFRs, consistent with the idea that adaptors limit signaling. Our results suggest that the relative stoichiometry of core MAPK pathway proteins is very similar across different cell types, with cell-specific differences mostly restricted to variable amounts of feedback regulators and receptors. The low abundance of adaptors relative to EGFR could be responsible for previous observations that only a fraction of total cell surface EGFR is capable of rapid endocytosis, high-affinity binding, and mitogenic signaling.

Integrated Proteogenomic Characterization of Human High-Grade Serous Ovarian Cancer.

To provide a detailed analysis of the molecular components and underlying mechanisms associated with ovarian cancer, we performed a comprehensive mass-spectrometry-based proteomic characterization of 174 ovarian tumors previously analyzed by The Cancer Genome Atlas (TCGA), of which 169 were high-grade serous carcinomas (HGSCs). Integrating our proteomic measurements with the genomic data yielded a number of insights into disease, such as how different copy-number alternations influence the proteome, the proteins associated with chromosomal instability, the sets of signaling pathways that diverse genome rearrangements converge on, and the ones most associated with short overall survival. Specific protein acetylations associated with homologous recombination deficiency suggest a potential means for stratifying patients for therapy. In addition to providing a valuable resource, these findings provide a view of how the somatic genome drives the cancer proteome and associations between protein and post-translational modification levels and clinical outcomes in HGSC. VIDEO ABSTRACT.

The landscape of viral proteomics and its potential to impact human health.

INTRODUCTION: Advances in mass spectrometry-based proteomic technologies are enhancing studies of viral pathogenesis. Identification and quantification of host and viral proteins and modifications in cells and extracellular fluids during infection provides useful information about pathogenesis, and will be critical for directing clinical interventions and diagnostics. AREAS COVERED: Herein we review and discuss a broad range of global proteomic studies conducted during viral infection, including those of cellular responses, protein modifications, virion packaging, and serum proteomics. We focus on viruses that impact human health and focus on experimental designs that reveal disease processes and surrogate markers. Expert commentary: Global proteomics is an important component of systems-level studies that aim to define how the interaction of humans and viruses leads to disease. Viral-community resource centers and strategies from other fields (e.g., cancer) will facilitate data sharing and platform-integration for systems-level analyses, and should provide recommended standards and assays for experimental designs and validation.

Reproducibility of Differential Proteomic Technologies in CPTAC Fractionated Xenografts.

The NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) employed a pair of reference xenograft proteomes for initial platform validation and ongoing quality control of its data collection for The Cancer Genome Atlas (TCGA) tumors. These two xenografts, representing basal and luminal-B human breast cancer, were fractionated and analyzed on six mass spectrometers in a total of 46 replicates divided between iTRAQ and label-free technologies, spanning a total of 1095 LC-MS/MS experiments. These data represent a unique opportunity to evaluate the stability of proteomic differentiation by mass spectrometry over many months of time for individual instruments or across instruments running dissimilar workflows. We evaluated iTRAQ reporter ions, label-free spectral counts, and label-free extracted ion chromatograms as strategies for data interpretation (source code is available from http://homepages.uc.edu/~wang2x7/Research.htm ). From these assessments, we found that differential genes from a single replicate were confirmed by other replicates on the same instrument from 61 to 93% of the time. When comparing across different instruments and quantitative technologies, using multiple replicates, differential genes were reproduced by other data sets from 67 to 99% of the time. Projecting gene differences to biological pathways and networks increased the degree of similarity. These overlaps send an encouraging message about the maturity of technologies for proteomic differentiation.