Construction and Validation of a Protein Prognostic Model for Lung Squamous Cell Carcinoma.

Lung squamous cell carcinoma (LUSCC), as the major type of lung cancer, has high morbidity and mortality rates. The prognostic markers for LUSCC are much fewer than lung adenocarcinoma. Besides, protein biomarkers have advantages of economy, accuracy and stability. The aim of this study was to construct a protein prognostic model for LUSCC. The protein expression data of LUSCC were downloaded from The Cancer Protein Atlas (TCPA) database. Clinical data of LUSCC patients were downloaded from The Cancer Genome Atlas (TCGA) database. A total of 237 proteins were identified from 325 cases of LUSCC patients based on the TCPA and TCGA database. According to Kaplan-Meier survival analysis, univariate and multivariate Cox analysis, a prognostic prediction model was established which was consisted of 6 proteins (CHK1_pS345, CHK2, IRS1, PAXILLIN, BRCA2 and BRAF_pS445). After calculating the risk values of each patient according to the coefficient of each protein in the risk model, the LUSCC patients were divided into high risk group and low risk group. The survival analysis demonstrated that there was significant difference between these two groups (p= 4.877e-05). The area under the curve (AUC) value of the receiver operating characteristic (ROC) curve was 0.699, which suggesting that the prognostic risk model could effectively predict the survival of LUSCC patients. Univariate and multivariate analysis indicated that this prognostic model could be used as independent prognosis factors for LUSCC patients. Proteins co-expression analysis showed that there were 21 proteins co-expressed with the proteins in the risk model. In conclusion, our study constructed a protein prognostic model, which could effectively predict the prognosis of LUSCC patients.

Multiclass cancer classification in fresh frozen and formalin-fixed paraffin-embedded tissue by DigiWest multiplex protein analysis.

Histomorphology and immunohistochemistry are the most common ways of cancer classification in routine cancer diagnostics, but often reach their limits in determining the organ origin in metastasis. These cancers of unknown primary, which are mostly adenocarcinomas or squamous cell carcinomas, therefore require more sophisticated methodologies of classification. Here, we report a multiplex protein profiling-based approach for the classification of fresh frozen and formalin-fixed paraffin-embedded (FFPE) cancer tissue samples using the digital western blot technique DigiWest. A DigiWest-compatible FFPE extraction protocol was developed, and a total of 634 antibodies were tested in an initial set of 16 FFPE samples covering tumors from different origins. Of the 303 detected antibodies, 102 yielded significant correlation of signals in 25 pairs of fresh frozen and FFPE primary tumor samples, including head and neck squamous cell carcinomas (HNSC), lung squamous cell carcinomas (LUSC), lung adenocarcinomas (LUAD), colorectal adenocarcinomas (COAD), and pancreatic adenocarcinomas (PAAD). For this signature of 102 analytes (covering 88 total proteins and 14 phosphoproteins), a support vector machine (SVM) algorithm was developed. This allowed for the classification of the tissue of origin for all five tumor types studied here with high overall accuracies in both fresh frozen (90.4%) and FFPE (77.6%) samples. In addition, the SVM classifier reached an overall accuracy of 88% in an independent validation cohort of 25 FFPE tumor samples. Our results indicate that DigiWest-based protein profiling represents a valuable method for cancer classification, yielding conclusive and decisive data not only from fresh frozen specimens but also FFPE samples, thus making this approach attractive for routine clinical applications.

Digital quantitative assessment of PD-L1 using digital spatial profiling.

The assessment of programmed death 1 ligand 1 (PD-L1) expression by Immunohistochemistry (IHC) is the US Food and Drug Administration (FDA)-approved predictive marker to select responders to checkpoint blockade anti-PD-1/PD-L1 axis immunotherapies. Different PD-L1 immunohistochemistry (IHC) assays use different antibodies and different scoring methods in tumor cells and immune cells. Multiple studies have compared the performance of these assays with variable results. Here, we investigate an alternative method for assessment of PD-L1 using a new technology known as digital spatial profiling. We use a previously described standardization tissue microarray (TMA) to assess the accuracy of the method and compare digital spatial profiler (DSP) to each FDA-approved PD-L1 assays, one LDT assay and three quantitative fluorescence assays. The standardized cell line Index tissue microarray contains 10 isogenic cells lines in triplicates expressing various ranges of PD-L1. The dynamic range of PD-L1 digital counts was measured in the ten cell lines on the Index TMA using the GeoMx DSP assay and read on the nCounter platform. The digital method shows very high correlation with immunohistochemistry scored with quantitative software and with quantitative fluorescence. High correlation of PD-L1 digital DSP counts were seen between rows on the same Index TMA. Finally, experiments from two Index TMAs showed reproducibility of DSP counts were independent of variable slide storage time over a three-week period after antibody labeling but before collection of cleaved tags. In summary, DSP appears to have quantitative potential comparable to quantitative immunohistochemistry. It is possible that this technology could be used as a PD-L1 protein measurement system for companion diagnostic testing for immune therapy.

Bayesian data integration and variable selection for pan-cancer survival prediction using protein expression data.

Accurate prognostic prediction using molecular information is a challenging area of research, which is essential to develop precision medicine. In this paper, we develop translational models to identify major actionable proteins that are associated with clinical outcomes, like the survival time of patients. There are considerable statistical and computational challenges due to the large dimension of the problems. Furthermore, data are available for different tumor types; hence data integration for various tumors is desirable. Having censored survival outcomes escalates one more level of complexity in the inferential procedure. We develop Bayesian hierarchical survival models, which accommodate all the challenges mentioned here. We use the hierarchical Bayesian accelerated failure time model for survival regression. Furthermore, we assume sparse horseshoe prior distribution for the regression coefficients to identify the major proteomic drivers. We borrow strength across tumor groups by introducing a correlation structure among the prior distributions. The proposed methods have been used to analyze data from the recently curated "The Cancer Proteome Atlas" (TCPA), which contains reverse-phase protein arrays-based high-quality protein expression data as well as detailed clinical annotation, including survival times. Our simulation and the TCPA data analysis illustrate the efficacy of the proposed integrative model, which links different tumors with the correlated prior structures.

RPPAware: A software suite to preprocess, analyze and visualize reverse phase protein array data.

Reverse Phase Protein Arrays (RPPA) is a high-throughput technology used to profile levels of protein expression. Handling the large datasets generated by RPPA can be facilitated by appropriate software tools. Here, we describe RPPAware, a free and intuitive software suite that was developed specifically for analysis and visualization of RPPA data. RPPAware is a portable tool that requires no installation and was built using Java. Many modules of the tool invoke R to utilize the statistical features. To demonstrate the utility of RPPAware, data generated from screening brain regions of a mouse model of Down syndrome with 62 antibodies were used as a case study. The ease of use and efficiency of RPPAware can accelerate data analysis to facilitate biological discovery. RPPAware 1.0 is freely available under GNU General Public License from the project website at http://downsyndrome.ucdenver.edu/iddrc/rppaware/home.htm along with a full documentation of the tool.

Integration of pan-cancer transcriptomics with RPPA proteomics reveals mechanisms of epithelial-mesenchymal transition.

Integrating data from multiple regulatory layers across cancer types could elucidate additional mechanisms of oncogenesis. Using antibody-based protein profiling of 736 cancer cell lines, along with matching transcriptomic data, we show that pan-cancer bimodality in the amounts of mRNA, protein, and protein phosphorylation reveals mechanisms related to the epithelial-mesenchymal transition (EMT). Based on the bimodal expression of E-cadherin, we define an EMT signature consisting of 239 genes, many of which were not previously associated with EMT. By querying gene expression signatures collected from cancer cell lines after small-molecule perturbations, we identify enrichment for histone deacetylase (HDAC) inhibitors as inducers of EMT, and kinase inhibitors as mesenchymal-to-epithelial transition (MET) promoters. Causal modeling of protein-based signaling identifies putative drivers of EMT. In conclusion, integrative analysis of pan-cancer proteomic and transcriptomic data reveals key regulatory mechanisms of oncogenic transformation.

Identification of differentially expressed peptides in high-throughput proteomics data.

With the advent of high-throughput proteomics, the type and amount of data pose a significant challenge to statistical approaches used to validate current quantitative analysis. Whereas many studies focus on the analysis at the protein level, the analysis of peptide-level data provides insight into changes at the sub-protein level, including splice variants, isoforms and a range of post-translational modifications. Statistical evaluation of liquid chromatography-mass spectrometry/mass spectrometry peptide-based label-free differential data is most commonly performed using a t-test or analysis of variance, often after the application of data imputation to reduce the number of missing values. In high-throughput proteomics, statistical analysis methods and imputation techniques are difficult to evaluate, given the lack of gold standard data sets. Here, we use experimental and resampled data to evaluate the performance of four statistical analysis methods and the added value of imputation, for different numbers of biological replicates. We find that three or four replicates are the minimum requirement for high-throughput data analysis and confident assignment of significant changes. Data imputation does increase sensitivity in some cases, but leads to a much higher actual false discovery rate. Additionally, we find that empirical Bayes method (limma) achieves the highest sensitivity, and we thus recommend its use for performing differential expression analysis at the peptide level.

Protein network construction using reverse phase protein array data.

In this paper, we introduce a novel computational method for constructing protein networks based on reverse phase protein array (RPPA) data to identify complex patterns in protein signaling. The method is applied to phosphoproteomic profiles of basal expression and activation/phosphorylation of 76 key signaling proteins in three breast cancer cell lines (MCF7, LCC1, and LCC9). Temporal RPPA data are acquired at 48h, 96h, and 144h after knocking down four genes in separate experiments. These genes are selected from a previous study as important determinants for breast cancer survival. Interaction networks are constructed by analyzing the expression levels of protein pairs using a multivariate analysis of variance model. A new scoring criterion is introduced to determine relevant protein pairs. Through a network topology based analysis, we search for wiring patterns to identify key proteins that are associated with significant changes in expression levels across various experimental conditions.

Development of multigene expression signature maps at the protein level from digitized immunohistochemistry slides.

Molecular classification of diseases based on multigene expression signatures is increasingly used for diagnosis, prognosis, and prediction of response to therapy. Immunohistochemistry (IHC) is an optimal method for validating expression signatures obtained using high-throughput genomics techniques since IHC allows a pathologist to examine gene expression at the protein level within the context of histologically interpretable tissue sections. Additionally, validated IHC assays may be readily implemented as clinical tests since IHC is performed on routinely processed clinical tissue samples. However, methods have not been available for automated n-gene expression profiling at the protein level using IHC data. We have developed methods to compute expression level maps (signature maps) of multiple genes from IHC data digitized on a commercial whole slide imaging system. Areas of cancer for these expression level maps are defined by a pathologist on adjacent, co-registered H&E slides, allowing assessment of IHC statistics and heterogeneity within the diseased tissue. This novel way of representing multiple IHC assays as signature maps will allow the development of n-gene expression profiling databases in three dimensions throughout virtual whole organ reconstructions.

A time-series DDP for functional proteomics profiles.

Using a new type of array technology, the reverse phase protein array (RPPA), we measure time-course protein expression for a set of selected markers that are known to coregulate biological functions in a pathway structure. To accommodate the complex dependent nature of the data, including temporal correlation and pathway dependence for the protein markers, we propose a mixed effects model with temporal and protein-specific components. We develop a sequence of random probability measures (RPM) to account for the dependence in time of the protein expression measurements. Marginally, for each RPM we assume a Dirichlet process model. The dependence is introduced by defining multivariate beta distributions for the unnormalized weights of the stick-breaking representation. We also acknowledge the pathway dependence among proteins via a conditionally autoregressive model. Applying our model to the RPPA data, we reveal a pathway-dependent functional profile for the set of proteins as well as marginal expression profiles over time for individual markers.

A Bayesian hierarchical model for identifying epitopes in peptide microarray data.

Peptide Microarray Immunoassay (PMI for brevity) is a novel technology that enables researchers to map a large number of proteomic measurements at a peptide level, providing information regarding the relationship between antibody response and clinical sensitivity. PMI studies aim at recognizing antigen-specific antibodies from serum samples and at detecting epitope regions of the protein antigen. PMI data present new challenges for statistical analysis mainly due to the structural dependence among peptides. A PMI is made of a complete library of consecutive peptides. They are synthesized by systematically shifting a window of a fixed number of amino acids through the finite sequence of amino acids of the antigen protein as ordered in the primary structure of the protein. This implies that consecutive peptides have a certain number of amino acids in common and hence are structurally dependent. We propose a new flexible Bayesian hierarchical model framework, which allows one to detect recognized peptides and bound epitope regions in a single framework, taking into account the structural dependence between peptides through a suitable latent Markov structure. The proposed model is illustrated using PMI data from a recent study about egg allergy. A simulation study shows that the proposed model is more powerful and robust in terms of epitope detection than simpler models overlooking some of the dependence structure.

Quantitative evaluation of protein expression as a function of tissue microarray core diameter: is a large (1.5 mm) core better than a small (0.6 mm) core?

CONTEXT: Tissue microarrays (TMAs) have emerged as a high-throughput technology for protein evaluation in large cohorts. This technique allows maximization of tissue resources by analysis of sections from 0.6-mm to 1.5-mm core "biopsies" of standard formalin-fixed, paraffin-embedded tissue blocks and by the processing of hundreds of cases arrayed on a single recipient block in an identical manner. OBJECTIVE: To assess the expression of a series of biomarkers as a function of core size. Although pathologists frequently feel better if larger core sizes are used, there is no evidence in the literature showing that large cores are better (or worse) than small cores for assessment of TMAs. DESIGN: Estrogen receptor, HER2/neu, epidermal growth factor receptor, STAT3, mTOR, and phospho-p70 S6 kinase were measured by immunofluorescence with automated quantitative analysis. One random 0.6-mm field (one 0.6-mm spot) was compared to 6 to 12 fields per spot, representing 1-mm and 1.5-mm cores, for 3 different tumor types. RESULTS: We show that measurement of a single random 0.6-mm spot was comparable to analysis of the whole 1-mm or 1.5-mm spot (Pearson R coefficient varying from 0.87-0.98) for all markers tested. CONCLUSIONS: Since TMA technology is now being used in all phases of biomarker development, this work shows that TMAs with 0.6-mm cores are as representative as those with any common larger core size for optimization of standardized experimental conditions. Given that a greater number of 0.6-cores can be arrayed in a single master block, use of this core size allows increased throughput and decreased cost.

Feature extraction and dimensionality reduction for mass spectrometry data.

Mass spectrometry is being used to generate protein profiles from human serum, and proteomic data obtained from mass spectrometry have attracted great interest for the detection of early stage cancer. However, high dimensional mass spectrometry data cause considerable challenges. In this paper we propose a feature extraction algorithm based on wavelet analysis for high dimensional mass spectrometry data. A set of wavelet detail coefficients at different scale is used to detect the transient changes of mass spectrometry data. The experiments are performed on 2 datasets. A highly competitive accuracy, compared with the best performance of other kinds of classification models, is achieved. Experimental results show that the wavelet detail coefficients are efficient way to characterize features of high dimensional mass spectra and reduce the dimensionality of high dimensional mass spectra.

A designed peptide chip: protein fingerprinting technology with a dry peptide array and statistical data mining.

There has recently been increased interest in the potential for microarray technologies to study protein networks in a whole cell system within a single experiment. Protein-detecting microarrays are composed of numerous agents immobilized within a tiny area on solid surfaces to capture targeted proteins and to detect interactions in a high-throughput fashion. In this chapter, in order to extend the usability of peptide microarrays, we describe a novel dry peptide microarray format to obtain protein fingerprint (PFP) data sets and a statistical PFP data manipulation technique to quantitatively analyze targeted proteins.

An advanced application of protein microarrays: cell-based assays for functional genomics.

Microarrays have become common tools for approaching different experimental questions: DNA, protein and peptide arrays offer the power of multiplexing the assay and by means of miniaturization technology, the possibility to reduce cost and amount of samples and reagents. Recently, a novel technology for functional assays has been proposed. Sabatini and co-workers have shown a cell-based microarrays method (1) that relies on the deposition and immobilization of an array of cDNA plasmids on a slide where cells are subsequently plated; the cDNA is then internalized by "reverse transfection" and cells overexpress or downregulate in each single spot the genes of interest. This approach allows the screening of different phenotypes in living cells of many genes in parallel on a single slide. To overcome some relevant limitations of this approach, we have implemented the technology by means of viral immobilization (2) on a novel surface of cluster-assembled nanostructured TiO2 (3) previously functionalized with an array of a docking protein. In this work, we present the detailed development of the "reverse infection cell-microarray based technology" in U2OS cells on a novel coated slide that represents an advanced application of protein arrays.

Qualitative and quantitative analysis of peptide microarray binding experiments using SVM-PEPARRAY.

A main objective of analyzing peptide array-based binding experiments is to uncover the relationship between a peptide sequence and the binding outcome. Limited by the peptide array technologies available for applications, few attempts have been made to construct qualitative or quantitative models that depict the peptide sequence:binding strength relationships in peptide microarray-based binding studies. There has been a long history of similar modeling efforts based on low-throughput binding data in the areas of T-cell epitope screening and kinase substrate mapping, however. The keen needs in peptide array applications and the success of the modeling efforts in related fields have prompted us to develop SVM-PEPARRAY, a Web-based program capable of constructing qualitative and quantitative models based on peptide microarray binding datasets using support vector machine (SVM) modeling methods. We expect that such modeling analysis will allow researchers to quickly extract sequence-based biological information from improved peptide array binding results and provide more precise and accurate information about the biological systems investigated.

PASE: a web-based platform for peptide/protein microarray experiments.

Peptide microarray technology requires bioinformatics and statistical tools to manage, store, and analyze the large amount of data produced. To address these needs, we developed a system called protein array software environment (PASE) that provides an integrated framework to manage and analyze microarray information from polypeptide chip technologies.

Messina: a novel analysis tool to identify biologically relevant molecules in disease.

BACKGROUND: Morphologically similar cancers display heterogeneous patterns of molecular aberrations and follow substantially different clinical courses. This diversity has become the basis for the definition of molecular phenotypes, with significant implications for therapy. Microarray or proteomic expression profiling is conventionally employed to identify disease-associated genes, however, traditional approaches for the analysis of profiling experiments may miss molecular aberrations which define biologically relevant subtypes. METHODOLOGY/PRINCIPAL FINDINGS: Here we present Messina, a method that can identify those genes that only sometimes show aberrant expression in cancer. We demonstrate with simulated data that Messina is highly sensitive and specific when used to identify genes which are aberrantly expressed in only a proportion of cancers, and compare Messina to contemporary analysis techniques. We illustrate Messina by using it to detect the aberrant expression of a gene that may play an important role in pancreatic cancer. CONCLUSIONS/SIGNIFICANCE: Messina allows the detection of genes with profiles typical of markers of molecular subtype, and complements existing methods to assist the identification of such markers. Messina is applicable to any global expression profiling data, and to allow its easy application has been packaged into a freely-available stand-alone software package.

MALDI/SELDI protein profiling of serum for the identification of cancer biomarkers.

The ability to visualize the full depth of the serum proteome in a high-throughput manner is a major goal of clinical proteomics. Methodologies, which combine higher throughput with the ability to observe differential protein expression levels, have been applied to this goal. An example of such a system is the coupling of robotic sample processing to matrix-assisted laser desorption time of flight mass spectrometry (MALDI-TOF-MS). Within this paradigm is a modification of MALDI-TOF termed surface-enhanced laser desorption/ionization-TOF (SELDI-TOF). Both conventional MALDI and SELDI have been used to generate protein expression profiles reflective of potential peptide changes in serum. This information can be used to identify proteins, which may enable new diagnostic and therapeutic strategies.

Open-source platform for the analysis of liquid chromatography-mass spectrometry (LC-MS) data.

The analysis of protein mixtures by liquid chromatography-mass spectrometry (LCMS) requires tools for viewing and navigating LC-MS data, locating peptides in LC-MS data, and eliminating low-quality peptides. msInspect, an open source platform, can carry out these steps for single experiments and can align and normalize peptide features in comparative studies with multiple LC-MS runs. In addition, msInspect can analyze quantitative studies with and without isotopic labels to generate peptide arrays.