PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of mcpMolecular & Cellular Proteomics : MCP
 
Mol Cell Proteomics. 2012 June; 11(6): M111.015479.
Published online 2012 February 15. doi:  10.1074/mcp.M111.015479
PMCID: PMC3433920

Human Biomarker Discovery and Predictive Models for Disease Progression for Idiopathic Pneumonia Syndrome Following Allogeneic Stem Cell Transplantation*An external file that holds a picture, illustration, etc.
Object name is sbox.jpg

Abstract

Allogeneic hematopoietic stem cell transplantation (SCT) is the only curative therapy for many malignant and nonmalignant conditions. Idiopathic pneumonia syndrome (IPS) is a frequently fatal complication that limits successful outcomes. Preclinical models suggest that IPS represents an immune mediated attack on the lung involving elements of both the adaptive and the innate immune system. However, the etiology of IPS in humans is less well understood. To explore the disease pathway and uncover potential biomarkers of disease, we performed two separate label-free, proteomics experiments defining the plasma protein profiles of allogeneic SCT patients with IPS. Samples obtained from SCT recipients without complications served as controls. The initial discovery study, intended to explore the disease pathway in humans, identified a set of 81 IPS-associated proteins. These data revealed similarities between the known IPS pathways in mice and the condition in humans, in particular in the acute phase response. In addition, pattern recognition pathways were judged to be significant as a function of development of IPS, and from this pathway we chose the lipopolysaccaharide-binding protein (LBP) protein as a candidate molecular diagnostic for IPS, and verified its increase as a function of disease using an ELISA assay. In a separately designed study, we identified protein-based classifiers that could predict, at day 0 of SCT, patients who: 1) progress to IPS and 2) respond to cytokine neutralization therapy. Using cross-validation strategies, we built highly predictive classifier models of both disease progression and therapeutic response. In sum, data generated in this report confirm previous clinical and experimental findings, provide new insights into the pathophysiology of IPS, identify potential molecular classifiers of the condition, and uncover a set of markers potentially of interest for patient stratification as a basis for individualized therapy.

Allogeneic hematopoietic stem cell transplantation (SCT)1 is the only curative therapy for many individuals with malignant and non malignant conditions. Unfortunately, a number of complications limit successful outcomes. Diffuse lung injury occurs in 25–55% of SCT recipients and can account for ~50% of all SCT related mortality (1, 2). Lung injury following SCT can be infectious or noninfectious in nature. Idiopathic pneumonia syndrome (IPS) is defined as widespread alveolar injury following SCT without an active lower respiratory tract infection (1, 2). The incidence of IPS ranges from 5–15% depending upon donor type and conditioning regimen intensity. The median time to onset of IPS is ~18 days, and the mortality rate is >70% despite aggressive treatment and supportive care (2, 3).

Several laboratories have used rodent models to elucidate the molecular mechanisms at play in the development of IPS (2). These models recapitulate many aspects of the clinical histopathology and physiological observations that are associated with human disease (reviewed in (2)). For example, histochemical staining of mouse lung tissue reveals dense mononuclear cell infiltrates consisting of macrophages and lymphocytes around pulmonary vessels and bronchioles. Increases in these mononuclear infiltrates are accompanied by high levels of tumor necrosis factor alpha (TNFα) in both lung tissue and broncho-alveolar lavage (BAL) fluid (4). Neutralization of TNFα via a number of strategies reduces the progression and severity of experimental IPS demonstrating a casual role of this protein in the pathophysiology of disease (58). These experimental insights have been translated into clinical studies that are currently evaluating the use of etanercept (Enbrel®; Amgen, Thousand Oaks, CA) a soluble TNF binding molecule for the treatment of IPS in both pediatric and adult SCT recipients (9).

Preclinical studies have helped refocus IPS away from an idiopathic clinical syndrome and toward a SCT complication wherein the lung is a target of both innate and adaptive immune dysregulation (2). To date, the majority of clinical research efforts involving IPS have focused primarily on describing disease characteristics with a few reports measuring immuno-modulatory proteins in BAL fluid and/or plasma via enzyme linked immunosorbent assays (ELISA) (9, 10). Although these studies have provided important insights into understanding the possible contribution of soluble cytokines to the inflammation engendered during IPS, the approaches taken thus far have not been comprehensive (they target only a few known inflammatory proteins) and are limited with respect to expanding our understanding of the disease pathway. Nevertheless these focused studies provide a set of expected results that should be replicated in unbiased, comprehensive studies.

Discovery-based, quantitative, proteomics provides a powerful, unbiased approach for identification and quantification of hundreds of proteins, and multiple reliable proteomic techniques are currently available (1114). Label-free protein expression analysis is one such peptide based (“bottom up”) proteomic technique, which capitalizes on the highly reproducible chromatography and parts-per-million mass accuracy available in current liquid chromatography/mass spectrometry (LC-MS) systems. Peptide species in these approaches are quantified by ion intensity, while individual peptides are grouped across samples based on precise mass and retention time measurements (15, 16). Here we describe two label-free protein expression analysis studies where key samples were derived from a clinical trial for the treatment of IPS. The goals of these two studies were different and thus the study designs and verification strategies varied, although both relied on label-free mass spectrometry as the main experimental platform. In the first experimental study we identified disease progression-associated proteins by analyzing peptide intensity changes over two time points (day of SCT and day of diagnosis) for cases who developed IPS compared with control patients studied at two comparable time points (day of SCT and 14–21 days later). In this case, the specific goal was to explore and potentially expand the disease pathway using molecular network analysis. The peptides/proteins detected as significantly changing in this analysis were verified by a number of methods including: pathways analysis, comparison to literature results including mouse models and human data, and verification of a specific protein of interest (LBP). The second label-free expression experiment was performed in order to build a predictive model of disease and treatment. Thus, we focused on predictors of both disease and treatment response using additional samples collected only at Day 0 of SCT. In this case a cross-validation (leave-one-out) strategy was used to understand classification power of the predictive biomarkers. The results of the first study both confirm previous clinical and experimental findings and provide new insights into the pathophysiology of IPS. The results of the second study suggest a predictive model that can identify, at the time of SCT, the patients who will progress to IPS and subclassify those IPS patients into responders and nonresponders to cytokine neutralization therapy. These markers are suitable for further study in a wider clinical setting.

EXPERIMENTAL PROCEDURES

Patients and Patient Controls

Samples analyzed in the label-free protein expression experiments were collected from patients with IPS who were enrolled on a clinical study conducted at the University of Michigan and the Dana Farber Cancer Institute. Information regarding patient exclusion criteria, the procedures used to diagnosis IPS, and therapeutic intervention with corticosteroids and etanercept have been described in detail (9). In brief, study patients were ≥1 year in age, recipients of allogeneic SCT and within 100 days of stem cell infusion. All patients had evidence of IPS as defined by NIH working group criteria (1, 2). There were no restrictions based upon underlying disease, stem cell source, transplant conditioning regimen, or donor/recipient HLA match. SCT control subjects received myelo-ablative conditioning and subsequently exhibited no signs of acute graft versus host disease, hepatic veno-occlusive disease, significant infection, or pulmonary dysfunction through day 100 following SCT. All patients (or their surrogates) and controls gave written, informed consent in accordance with the Declaration of Helsinki and the trial was approved by the respective Institutional Review Boards of the University of Michigan and the Dana-Farber Cancer Center.

Label-free Expression

Plasma Sample Collection and Preparation

Blood was obtained for label-free expression immediately after stem cell transplantation (Day 0) and at either the time of diagnosis (for patients with IPS) or between 14 and 21 days after transplant for patient controls. Specimens were collected in heparinzed tubes and plasma was separated and stored at −80 °C until analysis. Aliquots for individual plasma samples for both studies were thawed and subsequently depleted of the seven most abundant proteins using a 4.6 × 100 mm multiple affinity removal system (MARS Hu7, Agilent Technologies, Santa Clara, CA) according to manufacturer's instructions. Each depleted sample was concentrated (5,000 molecular weight cutoff, Millipore, Billerica, MA) and buffer exchanged with 50 mm Tris pH 8.8 to a final volume of ~100 μl. Total protein concentrations were determined by 2D Quant Kit as described by the manufacturer (GE Healthcare Piscataway, NJ) and 10 μg of flow through and bound fractions for each sample were loaded onto a one dimensional SDS-PAGE gel (4–20% Tris-HCL) as a quality control measure to confirm success of the depletion step. Subsequent to digestion, each sample was adjusted to 60 μg in 50 μl volume. Twenty microliters of 0.2% Rapigest (Waters, Milford, MA) was added followed by dithiotheritol to a final concentration of 5 mm. The samples were reduced at 80 °C for 15 min and cooled to room temperature prior to alkylation with iodacetamide at a final concentration of 10 mm for 30 min. Proteolytic digestion was performed with bovine trypsin (Promega, Madison, WI) with a final enzyme to protein ratio of 1:10 (w/w) for 18 h at 37 °C. Elisa methods for LBP were performed by the Immunology Core Laboratory at the University of Michigan as previously published (9) and per manufacturer's recommendations (Cell Sciences, Canton, MA).

Discovery Label-free Expression Data Set

The categorical factors under study consist of a Disease Factor (DF) crossed with a Time factor (TF). The two levels of the DF were the two groups of patients both of whom received allogeneic SCT and were either subsequently diagnosed with IPS (Disease - D) or experienced no complications (Controls - C). The two levels of the TF included the two time-points studied: day of SCT (Day-0) and either the day of IPS diagnosis (for IPS patients) or day 14 to 21 (Day-14) for controls. The later corresponds to the median time of onset of IPS. With this design, the goal of the first analysis was tightly focused on identifying proteins changing over both time and disease development (interaction). Eleven patients were included in this initial study. Specimens were collected at each time point, therefore a total of 22 observations (within 4 experimental groups) were analyzed with n = 5 to 6 biological replicates per group (see statistical section for power analysis as well as supplemental Fig. S1). Samples were randomized across the analysis (without blocking). There is no common reference sample. No technical replicates or sample pooling was performed.

Predictive Model Label-free Expression Data Set

The experimental design for establishing the disease predictive model included SCT patients who subsequently developed IPS (Disease-D) or those who experienced no complications (Controls- C). The experimental design for establishing the responsiveness predictive model included SCT patients who did not develop IPS (Controls-C), and those with IPS who ultimately responded to etanercept therapy (Responders–R) or not respond (Nonresponders - NR). Samples utilized in this model were those only collected at day of SCT with the goal established to define predictive and prognostic biomarker variables within the context of the aforementioned clinical trial. Technical replicates for 23 individual samples (14 from IPS patients; nine from controls) were utilized giving a total of 46 observations that were randomized across the analysis.

Liquid Chromatography and Mass Spectrometry

Six hundred nanograms of each sample were analyzed by LC/MS/MS using a Dionex Ultimate 3000 capillary LC system (Dionex Sunnyvale, CA) and a LTQ-Orbitrap XL (Thermo, Waltham, MA). The order of sample injections was randomized over all samples for each study. Separation and detection of peptides was performed as previously published (17). The instruments were mass calibrated immediately before the analysis using the instrument protocol. Raw LC/MS/MS data were processed via Proteomarker software version 1 (Infochromics, Toronto, Canada) for the discovery set and Rosetta Elucidator version 3.3.01 (Rosetta Biosoftware, Seattle, WA) for the predictive model set. The MS/MS peak lists were subsequently searched by Mascot version 2.2.0 (Matrix Science London, UK). The database used was the human International Protein Index (IPI) (68020 sequences). Search settings were as follows: trypsin enzyme specificity; mass accuracy window for precursor ion, 10 ppm; mass accuracy window for fragment ions, 0.8 Daltons; variable modification including carbamidomethylation of cysteines, 1 missed cleavage and oxidation of methionine. The criteria for peptide identification were a mass accuracy of ≤10 ppm and an expectation value of p ≤ 0.05 and an estimated False Discovery Rate (FDR) of 1% (18, 19). Proteins that had >2 peptides matching the above criteria were considered confirmed assignments whereas proteins identified with one peptide regardless of the Mascot score were highlighted as tentative assignments. Automated differential quantification of peptides in a set of samples was accomplished with Proteomarker or Rosetta Elucidator as previously described (17, 20, 21).

Data Quality Control

The discovery label-free experiment was processed for data quality control and dimension reduction in order to reduce missing value rates and error rates as previously described (17). Briefly, the data set was a matrix of n = 22 observations/individuals (rows) times p = 14,253 aligned m/z intensities (columns). From empirical evidence within the data, we assumed a bimodal distribution of consensus peptide sequence Mascot scores. A diffset for analysis was rejected if the consensus peptide sequence score was below the 85.7th percentile (≤ 25.35), which represents the lowest mode of distribution of Mascot scores (see supplemental Fig. S2). Out of these, 2024 diffsets were identified. An additional prefiltering step was carried out to select those diffsets for which the observed count of missing values per peptide satisfied two criteria simultaneously: (1) it had to be less than the total sample size n minus half of the minimum sample size per experimental group, (2) it had to maximize the difference between the total number of remaining peptides after selection and the total number of missing values. Here with n = 22 and an initial number of peptides p = 14253, taking the maximum mean value from the above criteria, the maximum number of missing events allowed per peptide was bounded to 19, and the final number of peptides retained was p = 1093. Despite the previous QC filtering steps, missing values still remained in the preselected data. These missing values represent a mixture of truly absent peptides whose expressions are to be taken as background “noise,” and truly low abundance peptides whose levels of expression are close to the limit of detection. To account for this nonignorable left-censoring issue (i.e. the nonrandom nature of the missingness mechanism at play) and its extent, we used a probability model adapted from Wang et al. (22), which describes “artifactual missing events” in this type of data. This model makes inferences on the missing values of one sample, based on the information from other “similar” samples (technical replicates or nearest neighbors). It substitutes a missing measurement of intensity with its expected value of the true intensity given that it is unobservable. The imputation model enabled quantification of low-abundance peptides and reduced the proportion of missing data to less than 10%. The remaining missing values after imputation represent truly absent peptides in the samples and were imputed by taking an estimate of the background noise. This amount of missing data compared very favorably to previously published label-free data sets (17).

The predictive model label-free experiment was processed using Rosetta Elucidator as previously described (20, 21, 23). The predictive data set was a matrix of n = 23 rows (observations/individuals) times p = 12,588 columns (aligned m/z intensities). Peptide and protein assignments were made using peptide and protein teller algorithms (24, 25) in Rosetta Elucidator with a false discovery rate of 1% (Rosetta Biosoftware, Seattle, WA). Quantification and imputation were performed only on those peptides that passed the above identification criteria, which yielded 1722 peptides.

Statistical Analysis

Data Transformation

To help remove sources of systematic variation (bias and variance) in the measured intensities of the above data sets because of experimental artifacts and to ensure that the usual assumptions for statistical inferences are met (normality), we first applied a log-transformation on the variables (peptides). In addition, because the homoscedasticity assumption in multi-group designs is also required e.g. in ANOVA models, we also applied our recently developed “joint adaptive mean-variance regularization” procedure as described in (26, 27), now available as an R package called “MVR” (28). Briefly, the joint adaptive regularization procedure simultaneously overcomes the lack of degrees of freedom and the variance-mean dependence issue in this type of data set where the number of variables hugely dominates the number of samples (26, 28). Although, the procedure is designed to stabilize the variance across variables (peptides), we observed that it also translates into good variance stabilization effect across sample groups in a multi-group design, as is the case in this study.

Power Analysis and Sample Size Calculation

We determined a sample size per group that would guarantee to detect a minimum effect of interest, at a certain power, while simultaneously controlling the positive False Discovery Rate (pFDR) (29, 30), as described in Liu et al. (31). After variance stabilization and normalization of the data (26), the usual distributional assumptions of test statistics become satisfied. Under the above assumption, all peptides are assumed to have approximately a common standard deviation σ across samples/groups and a common effect size Δ/σ (31). In the case of a balanced design, we calculated the group sample size n required to detect e.g. at least a 2.3-fold change (effect size Δ/σ = log2(2.3)/σ = 2) in the interaction effect between the Disease and Time factors as a function of power 1 − β, the parameter π0, interpreted as the proportion of nondifferentially expressed peptides, the common standard deviation σ, and for a fixed level α of pFDR. Based on preliminary data in our hands, estimates of the parameters π0 and of the common standard deviation σ of the peptides across experiments were [pi]0 ≈ 0.47–0.77 and [sigma with hat] = 0.84–1.00, respectively, which is consistent with estimates found in other platforms in high dimensional data (31). Under above assumptions and estimations, the sample size per group, required to detect at least a 2.3-fold change in the interaction effect, providing at least 80% power (1 − β = 0.8), and controlling for at most 5% pFDR (pFDR ≤ α = 0.05) was n = 5 (supplemental Fig. S1). Therefore, a minimum number of individual subjects to enroll in our study had to be at least 2n = 10 (for a total of 4n = 20 observations).

Discovery Label-free Protein Expression—Unsupervised Analyses

Potential groups and outliers among the samples were checked by a principal component analysis (PCA) (32). Clustering analysis was also performed using complete linkage hierarchical clustering and gap statistics to estimate the real number of clusters in the data (33).

Discovery Label-free Differential Expression Profiles of Peptide—Linear Mixed Model and Empirical Bayes Estimators for Statistical Inference

The thousands of retained peptides were analyzed by well-established techniques that have been developed for analyzing high-dimensional data. If we let yjkl be the intensity signal on the original scale of i-th individual and the j-th peptide, and zjkl = t(Yjkl) be the outcome on the transformed scale, where t(.) is the appropriate transformation mentioned earlier (26), a linear mixed ANOVA model for each individual peptide j is fitted as follows:

equation image

where Time (T)jk and Disease (D)jl and their interaction (TD)jkl are taken as the fixed effects of interest, whereas a Patient effect (P)ij is taken as a random effect, and where μj represents the average signal intensity for that peptide across all factors and observations, and the error term εijkl is assumed to be independent from the random factor (P)ij, both assumed to be normally distributed with mean 0 and some (unknown) variance components: εijklN(0,Σ); and (P)ijN(0,σP2). In this experimental design, contrasts were built for each of the fixed effects of interest, with coefficients to be estimated. Variables were ranked in order of evidence of differential expression. Corresponding p values were adjusted for multiple testing using a recent extension of the standard Benjamini-Hochberg procedure, which controls the expected False Discovery Rate (FDR) (34). This error rate, called the positive FDR (denoted pFDR) (29, 30), resulted in a procedure less conservative than the conventional FDR. As well as statistical significance, we report the estimated log2 transformed Fold Change (log2FC or M = log2(FC)). This represents the log-ratio (log2FC) for individual peptide between two experimental conditions in the case of a main effect and to a difference in log-ratios in the case of an interaction effect. An estimated average log2FC of –1 and 1 correspond to a ½ and twofold change respectively. Moderated t-, and B- statistics are also reported, these represent different measures of statistical significance. The Moderated t-statistic corresponds to the usual t-statistic except that information has been borrowed across variables (peptides), whereas the B-statistic is the empirical Bayes log2 of the posterior odds that the peptide is differentially expressed (35). Finally raw and adjusted p values (with FDR correction) are listed. Note that in every list all the peptides are ranked by adjusted p value and then by B-statistic. Implementations and algorithms of were from the freely available consortium CRAN. For linear modeling and supervised inferences, we used the package “limma” (36). Finally, for the control of the positive FDR, we used the package “qvalue” (37).

Predictive Model Building

Our multi-class prediction model uses a combination of “Spike and Slab” linear regression and Random Forest multiclass prediction models (38). Briefly, a multiclass random forest set of predictor variables (peptides) was built by following Ishwaran and Rao's approach for categorical response variable, where a combination of Spike and Slab Variable Selection (3840) and Random Forest Multiclass Classifier (41) is used. A pseudo Y response vector of response values is generated using each of the top N principal components (n = 10) of the peptide expression matrix X. The Generalized Elastic Net (4244) obtained from using Spike and Slab linear regression (3840, 45) is used to select the peptides, one regression fit at a time for each of the pseudo Y response vector. The final combined set of peptides is passed unto a random forest classifier to construct a multiclass random forest set of predictor peptides. This procedure is repeated B = 100 times to get “Variable Signatures” and “Model Sizes,” and estimate total and within-class “Prediction Error Rates.” Each Monte-Carlo replicate is based on balanced K-fold cross-validation (K = 3) with the data randomly split into a (K−1)/K set used for training and 1/K set for testing. For classification, the Prediction Error Rate is measured by the Misclassification Error Rate computed over the hold-out data (1/K). Full cross-validation was performed by including all steps of peptide selection and model construction (classifier) in order to avoid risks of selection bias and Prediction Error Rates bias. Classification models were selected by simultaneously minimizing the misclassification error rate (MER) and maximizing the area under the ROC curve (AUC). Note that the performance of a model peaks at a certain model size (models III and V in our results, Fig. 5) as overly parsimonious models tend to underfit the data (excess bias), and overly large model tend to overfit the data (excess variance) because of a degradation of the signal to noise ratio when marginally relevant proteins (variables) are added to the model. Computations were performed using the recent implementation of Spike & Slab models in the R library called “SpikeSlab,” available from the CRAN repository.

Fig. 5.
Predictive models for disease and response to etanercept therapy. A, Performance versus model size plots of disease (IPS) and responsiveness predictive models. Classification models were optimized by simultaneously minimizing the misclassification error ...

Ingenuity Pathway Analysis (IPA)

Following statistical analysis, significant peptides (p < FDR = 0.1) and their corresponding estimated log-ratios were imported into Ingenuity Pathway Analysis (Ingenuity Systems, Redwood City, CA). Once imported, a single log ratio value was generated at the protein level by determining the mean for all peptides identified for a specific protein.

RESULTS

Label-free Expression Analyses

The sample depletion protocol was reproducible across individual samples and yielded sufficient protein for label-free expression. Reproducible protein patterns via 1D-SDS-PAGE were observed across all samples for both flow-through and bound fractions (data not shown). These depleted samples were subsequently digested and analyzed by liquid chromatography tandem MS (LC/MS/MS). Fig. 1 highlights the chromatographic reproducibility in terms of intensity and retention time among samples collected from controls on Day 0 of SCT and 14 and IPS patients on Day 0 of SCT and the drastic differences in the chromatographic profile between these 3 groups compared with samples obtained from IPS patients at diagnosis. In addition, correlation analysis was performed on the peptide intensities analyzed in the predictive model label-free study for both biological and technical replicates. Excellent reproducibility was observed in this analysis for both technical and biological samples with mean correlations among control patient peptides of 0.94 and mean correlations of 0.90 for IPS patient peptides whereas technical replicates had mean correlations of 0.99. Additional detail for these correlations can be found in supplemental Fig. S3.

Fig. 1.
Chromatographic peptide expression patterns after SCT. Chromatographic reproducibility was observed among samples collected from control patients on day of SCT, controls on day 14 and IPS patients on day of SCT in contrast to distinct differences observed ...

Discovery Principal Component Analysis, Statistical Inferences and FDR Corrections

Grouping and outlier detection among the samples were evaluated by a PCA and clustering analysis (Fig. 2). The PCA biplot illustrates the separation of data into two distinct groups: one that included samples collected from patient controls at Day 0 and Day 14 and the IPS group at Day 0 whereas the other group contained all samples from the IPS patients at diagnosis. In addition, two potential outliers were observed in the IPS groups from Day 0 and diagnosis. These two samples were identified as biological outliers as they were collected from the same patient, but they were not removed from the analysis.

Fig. 2.
PCA analysis plots for label-free analysis. Two-dimensional bi-plot that uses the first two principal components to display information about the variables (peptides) and the relationship between samples as indicated by intersample distances. The position ...

The goal of the discovery study was to observe changes of protein expressions over time in the disease versus the control group. We therefore fitted the same linear mixed-effect model of a two-way analysis of variance univariately to each individual peptide. In this design, the interaction effect is reported as this contrast evaluates the Disease effect over Time. supplemental Table S1 highlights the peptides that were determined as having a significant change. Using a pFDR correction cutoff of 10%, 290 peptides were up (160) or down (130) regulated and no more than 29 of these selected peptides are expected to be false positives across the interaction contrast. These peptides represent 81 proteins whose abundance significantly changes over time with disease (supplemental Table S1). The volcano plot in Fig. 3 summarizes the distribution of peptides that were found to be significantly down (green) or up (red) regulated in this interaction. Table I highlights the subset of 81 proteins (from the peptides in supplemental Table S1) that met our confirmation criteria for protein identification described in Methods. Specific examples of significant immune-related proteins seen in the interaction set include lipopolysaccharide binding protein (LBP), lumican (LUM), attractin (ATRN), and multiple negative and positive acute phase reactant proteins like complement, apolipoprotein H (APOH), and C-reactive protein (CRP).

Fig. 3.
Volcano plots of interaction effect. The t test volcano plot arranges peptides by statistical significance (vertical axis) and magnitude of change (horizontal axis). The horizontal axis represents the estimated log fold change (denoted log2(FC) or M) ...
Table I
Confirmed proteins which were significant in the discovery label-free analysis and predictive models. The columns highlight the t statistic, number of peptides identified for a protein, percent protein coverage from identified peptides, FDR adjusted p ...

Ingenuity Pathway Analysis (IPA)

All significant proteins (interaction) were analyzed by IPA to understand pathways associated with the development of IPS in human subjects. Of the 81 proteins from Table I, the software mapped 72 of the proteins. Table II highlights the top canonical pathways, the number of proteins identified in each pathway, and their significance in the analysis. Many pathways and dysregulated proteins relate to inflammation and immune-dysregulation with the acute phase response being the top ranked canonical pathway. Fig. 4 summarizes the major connections in the acute phase response (TNFα/IL-6) signaling pathway and illustrates the corresponding acute phase proteins identified in the discovery study as significantly changing. The data indicates strong correlations observed for changes in plasma abundance of both negative (e.g. those that decrease) and positive (e.g. those that increase) acute phase reactants as a function of disease progression. For example, six of the eight negative acute phase reactants identified as significant in the interaction effect were observed to decrease in abundance whereas the plasma expression for 16 of 28 positive acute phase reactants identified were increased. Notably, this pathway is known to be dysregulated in mouse models of IPS (see above), thus this network analysis confirms the relevance of this response in humans. Additional canonical pathways identified in this analysis that are associated with immune modulation included: complement and coagulation system and pattern recognition of bacteria and viruses (Table II). Overall, the data suggests coordinate dysregulation in these associated pathways occurs during the development of IPS.

Table II
Top seven significant canonical pathways identified by pathway analysis. p value represents the significance of this pathway in the discovery label-free expression study as determined by IPA. The number of proteins correlates to the number that were identified ...
Fig. 4.
The acute phase canonical pathway and corresponding acute phase reactants identified as significant for interaction effect. Proteins highlighted in green correspond to a decrease in abundance and proteins highlighted in red correspond to an increase in ...

Predictive Models for Disease Progression to IPS and Response to Etanercept Therapy

The label-free expression analysis to derive prognostic and predictive biomarkers included 1722 peptides that passed identical prefiltering steps. These peptide variables were subjected to predictive model building as outlined in Methods. A total of eight different models, in which the number of peptide classifiers was varied, were evaluated for each prediction (disease and responsiveness). MER and AUC values are reported in Fig. 5. The performance ranged from a maximum AUC of 97% with 48 proteins to a minimum AUC of 80% using just three proteins in the model. Model III had the highest overall performance for prognosis (AUC = 97%) and included 95 peptides corresponding to 48 different proteins. Model V had the highest predictive (drug responsiveness) performance (AUC = 95%) and utilized 24 peptides corresponding to nine proteins. Moreover, these models accurately classified all patients at day of SCT who would progress to IPS and identified the subset that would respond to etanercept treatment (no false negatives). Supplemental Tables S2 and S3 highlight the peptides that were used in both models. Table I highlights the proteins that were identified as significant in the discovery analysis and were also found in the predictive models for progression to IPS and/or response to therapy. Interestingly, a number of proteins in the acute phase response signaling (APRS) pathway, which was the top canonical pathway in the discovery study, also showed significant differences in abundances at day 0 when comparing samples from control versus IPS patients as well as responders versus nonresponder subgroups. For instance both negative (APOH and histidine-rich glycoprotein-HRG) and positive (serum amyloid P component-APCS, complement, CRP, and Alpha 1 antichymotrypsin-SERPINA3) acute phase reactants identified in the discovery study were found in the diseases and response prediction models.

Fig. 6 highlights the mean abundance differences for selected proteins analyzed in the predictive model experiments that were interesting from the standpoint of inflammatory biology. For these proteins, T-tests were performed on the mean intensities between control and IPS patients and the responder versus nonresponder subset. Acute phase proteins APCS, inter-alpha trypsin inhibitor heavy chain 2 (ITIH2), and transrthyretin (TTR) were very significantly different in the control versus IPS comparison (p < 9e−04). Moreover, C3, ATRN, and LUM proteins, which are involved in the innate immune response were significantly different when comparing the responder versus nonresponder subset (p < 1e−08). The significant changes in abundance of these proteins suggest that dysregulation of innate immunity (and specifically the APRS pathway) occurs in patients who develop IPS and can be observed as early as the day of SCT. Therefore, these proteins have significant clinical value as they are candidate predictive biomarkers for patients with respect to both the development of disease (IPS) and/or responsiveness to (future) treatment.

Fig. 6.
Box-whisker plots for protein abundance at day of SCT for selected proteins across the clinical groups. The x axis for each plot represents the different groups of patients: Cont (red) = controls, Resp (green) = IPS patients who responded to etanercept ...

DISCUSSION

IPS remains a leading cause of treatment-related mortality following allogeneic SCT, but our understanding of the immune dysregulation responsible for pulmonary inflammation engendered in this context remains limited. Murine IPS models suggest that the lung is a target of two distinct, but inter-related pathways of immune-mediated injury: one involving soluble inflammatory effectors including TNFα, and the other driven by antigen specific, donor T-cell effectors (46). Based directly upon this preclinical work, we recently completed an early phase clinical trial using a soluble TNFα binding protein, etanercept (Enbrel®, Amgen, Thousand Oaks, CA) for the treatment of IPS. Results from this trial showed that the combination of etanercept and corticosteroids was associated with high response rates, improved survival, and acceptable toxicity in SCT recipients with IPS (9). The current report extends these initial biologic findings and further characterizes and compares the plasma proteome of patients with IPS to control patients who have an uncomplicated course early post-SCT.

One canonical pathway identified from the IPA network analysis as disease associated was the APRS pathway. The APRS pathway is a key component of the innate immune system that is triggered by disturbances in homeostasis resulting from injury and/or infection. The APRS pathway involves the altered expression of proteins synthesized in the liver whose plasma concentrations increase (positive acute phase reactants) or decrease (negative acute phase reactants) in response to circulating inflammatory cytokines such IL-6, IL-1, and TNFα (4749). For example, TNFα binds to its receptors, TNFRI and TNFRII, and initiates series of signaling cascades, which culminates in the activation of nuclear factor kappa B (NF-kB) a key transcription factor in the acute phase response. The activation of NF-kB alters the transcription of inflammatory molecules that is ultimately reflected by changes in their plasma concentrations (Fig. 4). Similarly, IL-6 can also serve as potent inducer of the APRS pathway. Upon binding to the IL-6 receptor, IL-6 stimulates hepatocytes to produce acute phase response proteins through a series of signaling cascades including the activation of STAT3. Thus, specific serum proteins can serve as “readouts” of the APRS pathway and provide insights into the status of key cytokines involved in its regulation (50). ELISA assays have shown significant increases in both BAL fluid and plasma concentrations of TNFα and IL-6 in IPS patients (9), and the proteome readout observed in this report provides additional evidence for broad-based activation of the APRS pathway as seen by the dysregulated expression of many relevant downstream molecules (Fig. 4).

Activation of the APRS pathway in the clinical setting as revealed here also reinforces experimental data establishing a causal role for TNFα in the development of IPS. TNFα levels are elevated in the BAL fluid of mice with IPS, and the administration of a soluble, dimeric, TNF binding protein (rhTNFR:Fc; Amgen Corp, Thousand Oaks, CA) from week 4 to week 6 after SCT reduced the progression of lung injury during this time period (4, 5). The use of TNFα deficient mice has shown that lung injury after allogeneic SCT is dependent upon donor, rather than host derived TNFα and that cytokine production from both donor accessory cells (macrophage/monocytes) and T cells significantly contributes to this toxicity (6). Donor derived TNFα serves as both an effector and facilitator of lung injury. TNFα contributes to pulmonary vascular endothelial cell apoptosis that accompanies pulmonary inflammation during the development of IPS (8). In addition, TNFα secreted by donor T cells regulates the chemokine milieu in the lung within the first 2 weeks after SCT, which directly contributes to the subsequent recruitment of monocytes and macrophages as lung injury progresses (6). Recent studies have revealed that receptor/ligand interactions between TNFα and TNFRII contribute to leukocyte recruitment during IPS by regulating pulmonary ICAM-1 expression (7).

The innate immune response is initiated in large part by antigen presenting cells (APCs) and myeloid-derived accessory cells (macrophages, monocytes, and neutrophils) both of which require activation via pattern recognition of evolutionary conserved structures on pathogens (51). These structures are also known as pathogen-associated molecular patterns (PAMPs) and are detected via pattern recognition molecules such as Toll like receptors (TLRs) and LBP (52, 53). Once detected, the engagement of PAMPs with TLRs on the cell surface triggers a cascade of signaling events that culminate in the production of proinflammatory cytokines (54). Examples of PAMPs include lipopolysaccharide (LPS), lipoteichoic acid, peptidoglycan, and mannan (51, 53). LPS, also known as endotoxin, is part of the outer membrane of Gram-negative bacteria and is a potent PAMP responsible for the inflammatory response during endotoxic shock (55). Circulating LPS liberated from Gram-negative bacteria binds to LBP in the blood stream and forms a complex with its co-receptor CD14 on the cell surface enabling transfer of LPS to MD2 and subsequently to TLR4, which is followed by TLR4 oligomerization and signaling (54). Experimental studies in animal models have suggested a critical role for LPS in cytokine amplification observed during the development of IPS. LPS levels are elevated in the BAL fluid of mice with IPS, and the intravenous administration of LPS to mice with advanced GVHD significantly increases the severity of lung injury (4, 5). The enhanced inflammation is associated with large increases of TNFα and LPS in the BAL fluid and the development of alveolar hemorrhage (5). Furthermore, direct antagonism of LPS early in the time course of SCT reduces systemic levels of TNFα and significantly decreases the severity of GVHD and IPS (56), a finding that is linked to the responsiveness of donor (rather than host) cells to LPS stimulation (57).

As a number of pattern recognition molecules were identified as significant (for interaction) in the label-free expression study (Table I), and pattern recognition of bacteria and viruses was a significant annotated pathway identified by the network analysis (Table II) we illustrate specific examples of protein level changes of pattern recognition molecules in Fig. 7. In particular, LBP levels are seen to have increased significantly for IPS patients. To verify the changes in this protein, we conducted ELISA experiments on LBP for a cohort of 24 patients from the clinical trial (9). These data in Table III show a roughly fivefold increase in LBP values at diagnosis versus either controls (day 14) or IPS patients at day 0. Although Fig. 7 also shows a significant increase in LBP as a function of IPS diagnosis, we wished to also evaluate the level of fold change in intensity from the proteomics data for more direct comparison to ELISA data. Thus, we summed the individual un-normalized intensities for the four LBP peptides to create a surrogate protein level intensity of LBP for the same three groups as examined for ELISA. These results are also shown in Table III and reflect a similar fold change to the ELISA data. Thus, by multiple means of measurement, LBP is seen to be significantly increased in IPS patients, and may represent a robust molecular diagnostic for the condition with a straightforward assay method already available. Biologically, LBP is a well characterized acute phase protein that binds endotoxin (58, 59). TLR4 activation and sensitivity to endotoxin relies on a series of interactions starting with LBP that enables optimal presentation of endotoxin to TLR4 (60, 61). In addition to bacterial recognition, LBP also can indirectly regulate the inflammatory response. Numerous studies in LBP knockout mice have shown that the absence of LBP leads to reduced LPS responsiveness and immune activation. Thus, it is an interesting candidate for future evaluation.

Fig. 7.
Interaction profile plots of top selected proteins. Interaction plots of four selected proteins for ATRN = Attractin, LBP = Lipopolysaccharide-binding protein, LUM = Lumican and PGLYRP2 Peptidoglycan recognition protein 2. The x axis for each plot represents ...
Table III
LBP ELISA and Proteomic Results. ELISA results were performed on 24 subjects and are reported in μg/ml. p values for control day 0 and IPS at Day 0 versus IPS Dx were 0.01. Proteomic results were derived from 11 subjects and are reported as unnormalized ...

Peptidoglycan recognition protein 2 (PGLYRP2) and lumican are pattern recognition proteins that had similar interaction profiles in our analysis; both showed increases in abundance in IPS patients at day 0 and a significant decrease at diagnosis when compared with controls at the same time point (Fig. 7). Whereas LBP has a well-established role in the innate immune response, the biologic functions of PGLYRP2 and lumican have only recently been identified. PGLYRPs participate in the innate immune response and contain a well conserved peptidoglycan binding type 2 amidase domain across most animals (6264). Mammalian PGLYRP2 is produced by the liver, and initial studies suggested that it functioned as an anti-inflammatory scavenger of peptidoglycan (65). However, recent studies in a peptidoglycan-induced arthritis animal model have demonstrated that PGLYRP2 has pro-inflammatory properties; PGLYRP2 was found to locally modulate the inflammatory response via cooperation with other pattern recognition molecules specifically Nod2 and TLR4 both of which are important regulators of inflammatory cytokine production (66). In our study, four peptides of PGLYRP2 were identified in the label-free analysis as significant for the interaction effect (Fig. 7).

Lumican is a small, leucine rich, extracellular matrix glycoprotein that has recently been shown to be involved in the innate immune response (67, 68). A study that investigated the role of lumican in wound healing found that mice deficient in lumican (Lum−/−) had poor recruitment of macrophages to the site of injury and a significant decrease in the induction of inflammatory cytokines (69). When challenged with LPS, Lum−/− mice are resistant to septic shock and death and are poor inducers of TNFα and IL-6 (70). This study also showed that lumican expression is induced during the innate immune response; LPS binds to lumican and co-precipitates with CD14, which is a key regulator of LPS sensing. Taken together, these data suggest lumican may enhance host sensitivity to microbial products like LPS. In the label-free analysis, four peptides for lumican were identified as significant for the interaction effect. Five of the nine IPS patients had higher levels at day 0 than controls at the same time point. A ~50% increase in average peptide abundance (non log scale) was observed for these subjects when compared with other subjects analyzed in this study. Interestingly, these IPS subjects all responded to etanercept therapy.

In addition to PAMP related proteins, attractin, a dipeptidyl peptidase IV- (DPP IV) like enzyme was significant in this analysis and shared a similar interaction profile with lumican and PLGYRP2 (Fig. 7). Attractin is a CUB (complement C1r/C1s, Uegf, Bmp1) family protease that has activity similar to DPP IV, and its activity is conserved in both the secreted and transmembrane forms (71). Interestingly, attractin is not expressed on resting T cells but is highly expressed on activated T cell and peripheral blood monocytes (71, 72). The enzymatic activity of attractin plays an important role in the cellular responses to stimuli as well as serving as modulator of T-cell activation and extravasation (71, 73, 74). Attractin may also aide in the formation of immuno-regulation clusters by a direct interaction with monocytes and T cells or by regulating the levels of circulating chemotactic cytokines (71); recent work conducted on peripheral blood monocytes demonstrate that attractin assists in adhesion to fibronectin and cytokine production in these cells (72). Many of the proteins discussed above, including the pattern recognition molecules, were significant for interaction in the discovery analysis and also good classifiers for disease and/or response to therapy. In this context, attractin peptides ranked in the top 10% of all peptide classifiers for response to therapy and had significant intensity values when comparing responder versus non responders (Fig. 6). Taken together, this data suggest that early immune dysregulation (as early as day 0 of SCT) can be identified in patients who ultimately develop IPS and when evaluated simultaneously, changes in key proteins can be leveraged to classify patients for effective treatment.

In sum, preclinical and early phase clinical findings supporting a role for TNFα in the development of IPS have paved the way for the development of two multicenter trails that are currently accruing pediatric (phase II trial) and adult (phase III trial) patients. The data reported herein extend this ongoing translational research effort utilizing label-free proteome expression analysis. We have characterized the plasma proteome of patients with and without IPS on day 0 and at the day of diagnosis of IPS or day 14 following SCT for controls. We have successfully identified a number of protein changes associated with IPS that are involved in innate immunity, pattern recognition and the APRS pathway many of which are consistent with observations from preclinical models whereas others have not previously been biologically linked to this disease process. Collectively, our results confirm and broaden what is known about the pathophysiology of human disease and firmly support mechanistic insights fostered by animal models. Moreover, statistical modeling identified a panel of proteins that can identify patients on day 0 of SCT who will ultimately develop IPS and predict their responsiveness to TNFα neutralization strategies. Importantly, insights generated from the proteomic analysis of samples collected at the bedside can now be (1) brought back to the bench to evaluate the specific contribution of novel proteins to the inflammation engendered during IPS using established animal systems, (2) used to direct validation efforts in other SCT patient cohorts, including those under treatment by etanercept, and ultimately (3) incorporated into the intelligent design of future clinical trials wherein stratification by anticipated response to treatment can form the basis of individualized therapy and improved outcomes for SCT recipients and a multitude of other patients with inflammatory conditions. The latter highlights the potential broader implications of our results.

Acknowledgments

We thank Guarav S. J. B. Rana, Kathleen C. Lundberg, Jennifer Burgoyne, and Xiaolin Li for technical support.

Footnotes

* This work was supported by Burroughs Wellcome Fund, The Leukemia and Lymphoma Society, and the Case Western Reserve University/Cleveland Clinic CTSA (Grant Number UL1 R024989) from the National Center for Research Resources (NCRR), a component of the National Institutes of Health and NIH Roadmap for Medical Research. We also acknowledge support from the Case Comprehensive Cancer Center Core Grant (P30-CA04370).

An external file that holds a picture, illustration, etc.
Object name is sbox.jpg This article contains supplemental Figs. S1 to S3 and Tables S1 to S3.

Conflict of interest disclosure: KRC and GY are co-PIs on clinical trials for which Amgen (the producers of etanercept) provides study drug and or central pharmacy support.

1 The abbreviations used are:

SCT
Allogeneic hematopoietic stem cell transplantation
IPS
Idiopathic pneumonia syndrome
TNFα
tumor necrosis factor alpha
BAL
broncho-alveolar lavage
ELISA
enzyme linked immunosorbent assays
LC-MS
liquid chromatography/mass spectrometry
TF
Time-Factor
DF
Disease-Factor
PCA
principal component analysis
Mixed two way ANOVA
linear mixed-effect model of analysis of variance
LBP
lipopolysaccharide binding protein
APOH
apolipoprotein H
CRP
C-reactive protein
HRG
histidine-rich glycoprotein-HRG
APCS
serum amyloid P component
SERPINA3
Alpha 1 antichymotrypsin
ITIH2
inter-alpha trypsin inhibitor heavy chain 2
TTR
transrthyretin
C3
Complement 3
ATRN
attractin
LUM
lumican
APRS
acute phase response signaling
NF-kB
nuclear factor kappa B
APCs
antigen presenting cells
PAMPs
pathogen-associated molecular patterns
TLRs
Toll like receptors
GVHD
graft vs host disease
PGLYRP2
peptidoglycan recognition protein 2
DPP IV
dipeptidyl peptidase IV.

REFERENCES

1. Clark J. G., Hansen J. A., Hertz M. I., Parkman R., Jensen L., Peavy H. H. (1993) NHLBI workshop summary. Idiopathic pneumonia syndrome after bone marrow transplantation. Am. Rev. Respir Dis. 147, 1601–1606 [PubMed]
2. Panoskaltsis-Mortari A., Griese M., Madtes D. K., Belperio J. A., Haddad I. Y., Folz R. J., Cooke K. R. (2011) An official American Thoracic Society research statement: noninfectious lung injury after hematopoietic stem cell transplantation: idiopathic pneumonia syndrome. Am. J. Respir Crit. Care Med. 183, 1262–1279 [PMC free article] [PubMed]
3. Kantrow S. P., Hackman R. C., Boeckh M., Myerson D., Crawford S. W. (1997) Idiopathic pneumonia syndrome: changing spectrum of lung injury after marrow transplantation. Transplantation 63, 1079–1086 [PubMed]
4. Cooke K. R., Kobzik L., Martin T. R., Brewer J., Delmonte J., Jr., Crawford J. M., Ferrara J. L. (1996) An experimental model of idiopathic pneumonia syndrome after bone marrow transplantation: I. The roles of minor H antigens and endotoxin. Blood 88, 3230–3239 [PubMed]
5. Cooke K. R., Hill G. R., Gerbitz A., Kobzik L., Martin T. R., Crawford J. M., Brewer J. P., Ferrara J. L. (2000) Tumor necrosis factor-alpha neutralization reduces lung injury after experimental allogeneic bone marrow transplantation. Transplantation 70, 272–279 [PubMed]
6. Hildebrandt G. C., Olkiewicz K. M., Corrion L. A., Chang Y., Clouthier S. G., Liu C., Cooke K. R. (2004) Donor-derived TNF-alpha regulates pulmonary chemokine expression and the development of idiopathic pneumonia syndrome after allogeneic bone marrow transplantation. Blood 104, 586–593 [PubMed]
7. Hildebrandt G. C., Olkiewicz K. M., Corrion L. A., Clouthier S. G., Pierce E., Liu C., Cooke K. R. (2008) A role for TNF receptor type II in leukocyte infiltration into the lung during experimental idiopathic pneumonia syndrome. Biol. Blood Marrow Transplant. 14, 385–396 [PMC free article] [PubMed]
8. Gerbitz A., Nickoloff B. J., Olkiewicz K., Willmarth N. E., Hildebrandt G., Liu C., Kobzik L., Eissner G., Holler E., Ferrara J. L., Cooke K. R. (2004) A role for tumor necrosis factor-alpha-mediated endothelial apoptosis in the development of experimental idiopathic pneumonia syndrome. Transplantation 78, 494–502 [PubMed]
9. Yanik G. A., Ho V. T., Levine J. E., White E. S., Braun T., Antin J. H., Whitfield J., Custer J., Jones D., Ferrara J. L., Cooke K. R.(2008) The impact of soluble tumor necrosis factor receptor etanercept on the treatment of idiopathic pneumonia syndrome after allogeneic hematopoietic stem cell transplantation. Blood 112, 3073–3081 [PubMed]
10. Clark J. G., Madtes D. K., Martin T. R., Hackman R. C., Farrand A. L., Crawford S. W. (1999) Idiopathic pneumonia after bone marrow transplantation: cytokine activation and lipopolysaccharide amplification in the bronchoalveolar compartment. Crit. Care Med. 27, 1800–1806 [PubMed]
11. Rifai N., Gillette M. A., Carr S. A. (2006) Protein biomarker discovery and validation: the long and uncertain path to clinical utility. Nat. Biotechnol. 24, 971–983 [PubMed]
12. Veenstra T. D., Conrads T. P., Hood B. L., Avellino A. M., Ellenbogen R. G., Morrison R. S. (2005) Biomarkers: mining the biofluid proteome. Mol. Cell Proteomics 4, 409–418 [PubMed]
13. Zhao Y., Lee W. N., Xiao G. G. (2009) Quantitative proteomics and biomarker discovery in human cancer. Expert Rev. Proteomics 6, 115–118 [PMC free article] [PubMed]
14. Paczesny S., Krijanovski O. I., Braun T. M., Choi S. W., Clouthier S. G., Kuick R., Misek D. E., Cooke K. R., Kitko C. L., Weyand A., Bickley D., Jones D., Whitfield J., Reddy P., Levine J. E., Hanash S. M., Ferrara J. L.(2009) A biomarker panel for acute graft-versus-host disease. Blood 113, 273–278 [PubMed]
15. Wang W., Zhou H., Lin H., Roy S., Shaler T. A., Hill L. R., Norton S., Kumar P., Anderle M., Becker C. H. (2003) Quantification of proteins and metabolites by mass spectrometry without isotopic labeling or spiked standards. Anal. Chem. 75, 4818–4826 [PubMed]
16. Chelius D., Bondarenko P. V. (2002) Quantitative profiling of proteins in complex mixtures using liquid chromatography and mass spectrometry. J. Proteome Res. 1, 317–323 [PubMed]
17. Schlatzer D. M., Dazard J. E., Dharsee M., Ewing R. M., Ilchenko S., Stewart I., Christ G., Chance M. R. (2009) Urinary protein profiles in a rat model for diabetic complications. Mol. Cell Proteomics 8, 2145–2158 [PMC free article] [PubMed]
18. Elias J. E., Haas W., Faherty B. K., Gygi S. P. (2005) Comparative evaluation of mass spectrometry platforms used in large-scale proteomics investigations. Nat. Methods 2, 667–675 [PubMed]
19. Brosch M., Yu L., Hubbard T., Choudhary J. (2009) Accurate and sensitive peptide identification with Mascot Percolator. J. Proteome Res. 8, 3176–3181 [PMC free article] [PubMed]
20. Chan E. Y., Sutton J. N., Jacobs J. M., Bondarenko A., Smith R. D., Katze M. G. (2009) Dynamic host energetics and cytoskeletal proteomes in human immunodeficiency virus type 1-infected human primary CD4 cells: analysis by multiplexed label-free mass spectrometry. J. Virol. 83, 9283–9295 [PMC free article] [PubMed]
21. Neubert H., Bonnert T. P., Rumpel K., Hunt B. T., Henle E. S., James I. T. (2008) Label-free detection of differential protein expression by LC/MALDI mass spectrometry. J. Proteome Res. 7, 2270–2279 [PubMed]
22. Wang P., Tang H., Zhang H., Whiteaker J., Paulovich A. G., McIntosh M. (2006) Normalization regarding non-random missing values in high-throughput mass spectrometry data. In Pac Symp Biocomput. Lihue, Hawaii: World Scientific Press; 315–326 [PubMed]
23. Weng L., Dai H., Zhan Y., He Y., Stepaniants S. B., Bassett D. E. (2006) Rosetta error model for gene expression analysis. Bioinformatics 22, 1111–1121 [PubMed]
24. Nesvizhskii A. I., Keller A., Kolker E., Aebersold R. (2003) A statistical model for identifying proteins by tandem mass spectrometry. Anal. Chem. 75, 4646–4658 [PubMed]
25. Keller A., Nesvizhskii A. I., Kolker E., Aebersold R. (2002) Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal. Chem. 74, 5383–5392 [PubMed]
26. Dazard J.-E., Rao J. S. (2012) Joint Adaptive Mean-Variance Regularization and Variance Stabilization of High Dimensional Data. Comput. Statist. Data Anal. (in press) [PMC free article] [PubMed]
27. Dazard J.-E., Xu H., Rao J. S. (2011) R package MVR for Joint Adaptive Mean-Variance Regularization and Variance Stabilization. In JSM Proceedings. Section for Statistical Programmers and Analysts Miami Beach, FL. U.S.A.: American Statistical Association; (in press)
28. Dazard J.-E., Xu H., Santana A., Rao J. S. (2011) R package MVR. In Comprehensive R Archive Network
29. Storey J. D. (2002) A direct approach to false discovery rates. J. R. Statist. Soc. 64, 479–498
30. Storey J. D. (2003) The positive false discovery rate: A Bayesian interpretation and the q-value. Ann. Statistics 31, 2013–2035
31. Liu P., Hwang J. T. (2007) Quick calculation for sample size while controlling false discovery rate with application to microarray analysis. Bioinformatics 23, 739–746 [PubMed]
32. Hotelling H. (1933) Analysis of a complex of statistical variables into principal components. J. Ed. Psychol. 24, 417–441
33. Tibshirani R. (2001) Estimating the number of clusters in a data set via the gap statistic. J. R. Statist. Soc. 63, 411–423
34. Benjamini Y., Hochberg Y. (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Statist. Soc. 57, 289–300
35. Smyth G. K. (2004) Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl. Genet Mol. Biol. 3, Article3 [PubMed]
36. Smyth G. K., Thorne N., Wettenhall J. (2004) Limma R package: linear models for microarray. Comprehensive R Archive Network
37. Dabney A., Storey J. D. (2003) Qvalue R package: Q-value estimation for false discovery rate control. Comprehensive R Archive Network
38. Ishwaran H., Rao J. S. (2011) Generalized ridge regression: geometry and computational solutions when p is larger than n. Miami, University of Miami; Pg. 1–39
39. Ishwaran H., Rao J. S. (2005) Spike and slab variable selection: frequentist and Bayesian strategies. Ann. Statistics 33, 730–773
40. Ishwaran H., Rao J. S. (2005) Spike and slab gene selection for multigroup microarray data. J. Amer. Stat Assoc. 100, 764–780
41. Breiman L. (2001) Random forests. Mach. Learn. 45, 5–32
42. Efron B., Hastie T., Johnstone I., Tibshirani R. (2004) Least angle regression. Ann. Statistics 32, 407–499
43. Zou H., Hastie T. (2005) Regularization and variable selection via the elastic net. J. R Statist Soc. 67, 301–320
44. Friedman J., Hastie T., Tibshirani R. (2010) Regularization Paths for Generalized Linear Models via Coordinate Descent. Journal of Statistical Software 33, 1–22 [PMC free article] [PubMed]
45. Ishwaran H., Rao J. S. (2003) Detecting differentially expressed genes in microarrays using Bayesian model selection. J. Amer. Stat. Assoc. 98, 438–455
46. Panoskaltsis-Mortari A., Griese M., Madtes D. K., Belperio J. A., Haddad I. Y., Folz R. J., Cooke K. R. (2011) An official American Thoracic Society research statement: noninfectious lung injury after hematopoietic stem cell transplantation: idiopathic pneumonia syndrome. Am. J. Respir Crit. Care Med. 183, 1262–1279 [PMC free article] [PubMed]
47. Ihnatko R., Kubes M. (2007) TNF signaling: early events and phosphorylation. Gen. Physiol. Biophys. 26, 159–167 [PubMed]
48. Heinrich P. C., Castell J. V., Andus T. (1990) Interleukin-6 and the acute phase response. Biochem. J. 265, 621–636 [PubMed]
49. Schuett H., Luchtefeld M., Grothusen C., Grote K., Schieffer B. (2009) How much is too much? Interleukin-6 and its signalling in atherosclerosis. Thromb. Haemost. 102, 215–222 [PubMed]
50. Chen Y., Lim B. K., Peh S. C., Abdul-Rahman P. S., Hashim O. H. (2008) Profiling of serum and tissue high abundance acute-phase proteins of patients with epithelial and germ line ovarian carcinoma. Proteome Sci. 6, 20. [PMC free article] [PubMed]
51. Medzhitov R., Janeway C., Jr. (2000) Innate immune recognition: mechanisms and pathways. Immunol. Rev. 173, 89–97 [PubMed]
52. Heumann D., Lauener R., Ryffel B. (2003) The dual role of LBP and CD14 in response to Gram-negative bacteria or Gram-negative compounds. J. Endotoxin Res. 9, 381–384 [PubMed]
53. Mogensen T. H. (2009) Pathogen recognition and inflammatory signaling in innate immune defenses. Clin. Microbiol. Rev. 22, 240–273, Table of Contents [PMC free article] [PubMed]
54. Akira S., Uematsu S., Takeuchi O. (2006) Pathogen recognition and innate immunity. Cell 124, 783–801 [PubMed]
55. Trent M. S., Stead C. M., Tran A. X., Hankins J. V. (2006) Diversity of endotoxin and its impact on pathogenesis. J. Endotoxin Res. 12, 205–223 [PubMed]
56. Cooke K. R., Gerbitz A., Crawford J. M., Teshima T., Hill G. R., Tesolin A., Rossignol D. P., Ferrara J. L. (2001) LPS antagonism reduces graft-versus-host disease and preserves graft-versus-leukemia activity after experimental bone marrow transplantation. J. Clin. Invest. 107, 1581–1589 [PMC free article] [PubMed]
57. Cooke K. R., Hill G. R., Crawford J. M., Bungard D., Brinson Y. S., Delmonte J., Jr., Ferrara J. L. (1998) Tumor necrosis factor- alpha production to lipopolysaccharide stimulation by donor cells predicts the severity of experimental acute graft-versus-host disease. J. Clin. Invest. 102, 1882–1891 [PMC free article] [PubMed]
58. Gioannini T. L., Teghanemt A., Zhang D., Levis E. N., Weiss J. P. (2005) Monomeric endotoxin:protein complexes are essential for TLR4-dependent cell activation. J. Endotoxin Res. 11, 117–123 [PubMed]
59. Jerala R. (2007) Structural biology of the LPS recognition. Int. J. Med. Microbiol. 297, 353–363 [PubMed]
60. Beutler B., Rietschel E. T. (2003) Innate immune sensing and its roots: the story of endotoxin. Nat. Rev. Immunol. 3, 169–176 [PubMed]
61. Ulevitch R. J., Tobias P. S. (1999) Recognition of gram-negative bacteria and endotoxin by the innate immune system. Curr. Opin. Immunol. 11, 19–22 [PubMed]
62. Kang D., Liu G., Lundström A., Gelius E., Steiner H. (1998) A peptidoglycan recognition protein in innate immunity conserved from insects to humans. Proc. Natl. Acad. Sci. U.S.A. 95, 10078–10082 [PubMed]
63. Dziarski R., Gupta D. (2006) The peptidoglycan recognition proteins (PGRPs). Genome Biol. 7, 232. [PMC free article] [PubMed]
64. Werner T., Liu G., Kang D., Ekengren S., Steiner H., Hultmark D. (2000) A family of peptidoglycan recognition proteins in the fruit fly Drosophila melanogaster. Proc. Natl. Acad. Sci. U.S.A. 97, 13772–13777 [PubMed]
65. Gelius E., Persson C., Karlsson J., Steiner H. (2003) A mammalian peptidoglycan recognition protein with N-acetylmuramoyl-L-alanine amidase activity. Biochem. Biophys. Res. Commun. 306, 988–994 [PubMed]
66. Saha S., Qi J., Wang S., Wang M., Li X., Kim Y. G., Nuñez G., Gupta D., Dziarski R. (2009) PGLYRP-2 and Nod2 are both required for peptidoglycan-induced arthritis and local inflammation. Cell Host Microbe. 5, 137–150 [PMC free article] [PubMed]
67. Chakravarti S. (2002) Functions of lumican and fibromodulin: lessons from knockout mice. Glycoconj J. 19, 287–293 [PubMed]
68. Iozzo R. V. (1999) The biology of the small leucine-rich proteoglycans. Functional network of interactive proteins. J. Biol. Chem. 274, 18843–18846 [PubMed]
69. Vij N., Roberts L., Joyce S., Chakravarti S. (2005) Lumican regulates corneal inflammatory responses by modulating Fas-Fas ligand signaling. Invest Ophthalmol Vis. Sci. 46, 88–95 [PubMed]
70. Wu F., Vij N., Roberts L., Lopez-Briones S., Joyce S., Chakravarti S. (2007) A novel role of the lumican core protein in bacterial lipopolysaccharide-induced innate immune response. J. Biol. Chem. 282, 26409–26417 [PubMed]
71. Duke-Cohan J. S., Tang W., Schlossman S. F. (2000) Attractin: a cub-family protease involved in T cell-monocyte/macrophage interactions. Adv. Exp. Med. Biol. 477, 173–185 [PubMed]
72. Wrenger S., Faust J., Friedrich D., Hoffmann T., Hartig R., Lendeckel U., Kähne T., Thielitz A., Neubert K., Reinhold D. (2006) Attractin, a dipeptidyl peptidase IV/CD26-like enzyme, is expressed on human peripheral blood monocytes and potentially influences monocyte function. J. Leukoc Biol. 80, 621–629 [PubMed]
73. Ohnuma K., Dang N. H., Morimoto C. (2008) Revisiting an old acquaintance: CD26 and its molecular mechanisms in T cell function. Trends Immunol. 29, 295–301 [PubMed]
74. Duke-Cohan J. S., Gu J., McLaughlin D. F., Xu Y., Freeman G. J., Schlossman S. F. (1998) Attractin (DPPT-L), a member of the CUB family of cell adhesion and guidance proteins, is secreted by activated human T lymphocytes and modulates immune cell interactions. Proc. Natl. Acad. Sci. U.S.A. 95, 11336–11341 [PubMed]

Articles from Molecular & Cellular Proteomics : MCP are provided here courtesy of American Society for Biochemistry and Molecular Biology