Lung Cancer Diagnosis from Proteomic Analysis of Preinvasive Lesions

Informative biomarkers for early detection would be of great utility in improving the outcome of patients with lung cancer, the leading cause of cancer deaths in the United States and worldwide (1). Recent attempts at early detection strategies using low-dose chest CT are showing promise in reducing lung cancer-related mortality (2). Yet, continued research is needed in assessing individuals for their risk in developing and for having lung cancer (3, 4).

Although only a minority of preinvasive lung lesions progress to cancer, they are associated with an increased risk of having lung cancer (5–10). These lesions are likely to harbor molecular alterations indicative of lung cancer risk in the entire affected epithelium, including tumors of different histological subtypes (11, 12). However, this histology-based assessment of these lesions is limited in many ways. First, grading is semiquantitative at best and the agreement between 2 observers for the reporting of the histopathology according to the World Health Organization (WHO) criteria has not been studied extensively (13, 14). Distinguishing mild from moderate dysplasia and severe dysplasia from carcinoma in situ is difficult, for example. Consequently, most studies evaluate bronchial preinvasive lesions that are either low or high grade. Second, the effect of the structural integrity of the biopsy on the interpretability of the histology features remains a significant obstacle to this area of research. Thus, the development of a quantitative measure of probability of having lung cancer based on the molecular analysis of the bronchial epithelium, and independent of the histological grade, may have important implications. Our study was therefore designed to measure the probability of having a lung cancer based on the proteomic analysis of endobronchial lesions.

We recently identified patterns of protein expression in fresh frozen bronchial lesions obtained by autofluorescence bronchoscopy and resected lung tissue samples at different stages of tumor progression using matrix-assisted laser desorption ionization mass spectrometry (MALDI MS; ref. 15). Our proteomic analysis generated signatures discriminating subtypes of lesions in a continuum from normal lung to invasive lung tumor.

Here we took our original observation further and hypothesized that the bronchial epithelium expresses a subset of proteins that once increased confer an increased probability of having lung cancer. To assess the probability of having a lung cancer based on the analysis of bronchial lesions, we are moving away from a pathological scale based on the grade of the lesions to a quantitative scale based on the proteomic analysis of these bronchial specimens. We selected 9 MALDI MS features [mass-to-charge ratio (m/z)] as a signature from 2 previous studies (15, 16). The novel proteomic signature of 9 features gives rise to a score based on MALDI MS signal intensity that is used as a continuous variable to predict the diagnosis of lung cancer and circumvents many of the problems associated with a histology-based assessment of the lesions. We identified 8 of the 9 features and tested their expression by immunohistochemistry (IHC) in formalin fixed paraffin-embedded (FFPE) tissues and Western blotting.

The study design is summarized in Figure 1. To derive a bronchial epithelium-based signature predictive of lung cancer, results from 2 independent studies published from our laboratory were examined for top ranking discriminatory features (m/z) identified by MALDI MS. The first study led to features distinguishing normal lung from invasive lung tumors (n = 93; ref. 16), and the second study led to features distinguishing normal bronchial epithelium, preinvasive bronchial lesions, and invasive lung tumors (n = 51; ref. 15). Nine features (m/z) were selected from these 144 patients as candidates of a novel signature from the intercept of best classifiers from these 2 previous studies based on their statistical significance with false discovery rate-adjusted P < 0.01 and based on the expert visual confirmation of the characteristics of the peak. The m/z values of the selected features were 4,749, 4,965, 6,175, 8,451, 8,565, 9,955, 11,048, 12,275, and 12,345. The features (m/z) were identified and their expression was verified by biochemical methods. To evaluate the proteomic signature of bronchial lesions as a biomarker diagnostic of lung cancer, we built a prediction model using the Rahman and colleagues cohort (51 patients, 41 with and 10 without lung cancer; ref. 15) and validated it in a new data set (60 patients, 40 with and 20 without lung cancer). The intensity of the signature (9 m/z values) was transformed to a MALDI MS score to predict the diagnosis of having lung cancer (see Statistical Analysis).

The project was approved by the local Institutional Review Board and informed consent was obtained on all individuals at all institutions.

A total of 69 bronchial biopsies and surgical lung tissue specimens were collected from the 60 patients in the validation cohort (Table 1 and Supplementary Table S1) who participated between 2003 and 2008 at Vanderbilt University Medical Center (VUMC), the Nashville Veteran Affairs Medical Center (VAMC), or the University of Colorado Health Science Center (UCHSC). Multiple bronchial specimens were obtained at the time of bronchoscopy or surgical resection. Patients undergoing bronchoscopy for clinical suspicion of lung cancer were consented for this research study and agreed to undergo autofluorescence bronchoscopy and to provide bronchial specimens at predetermined normal sites (with normal fluorescence ratio) and at suspicious sites (with increased fluorescence ratio; ref. 17) in addition to the samples taken for clinical diagnosis. Patients with lung cancer in both the validation cohort and the Rahman and colleagues (15) cohort included different subtypes and all stages of lung cancer (Supplementary Tables S1 and S2). As shown in these supplementary tables, some patients in both cohorts provided more than 1 tissue specimen. Specimens collected for research were snap frozen and later cut, stained, and graded by the pathologists (W.A.F. and A.L.G.) according to the WHO nomenclature (14). Only bronchial epithelium specimens were considered in the analysis so that no alveolar tissue sample was used in either the training or test set. Normal bronchial epithelium specimens were obtained from surgical specimens at least 2 cm away from the tumor. Two pathologists (A.L.G. and W.A.F.) reviewed the hematoxylin and eosin (H&E) stained sections of each lesion.

Bronchial biopsy specimens were frozen in liquid nitrogen immediately after collection, embedded in Tissue-Tek Optimal Cutting Temperature compound (IMEB) and stored at −80°C. For MALDI MS profiling, a 12 μm section was cut and thaw-mounted directly onto a gold-coated MALDI plate followed by fixation in graded ethanol washes (70%, 90%, and 95% for 30 seconds each). An adjacent 7 μm section was also cut and stained with H&E following standard protocols. Photomicrographs of H&E stained tissue sections were obtained with a Mirax Scan digital microscope slide scanner (Mirax) at a pixel resolution of 0.23 μm. The serial 12 μm histology images were annotated using Mirax viewer software by a pathologist (A.L.G.) using a small circular shape (200 μm) to mark areas of interest for matrix spotting (Supplementary Fig. S1). Color coding was used to designate areas of different histology grade within the tissue section. Annotated histology images were then overlaid with the MALDI plate images using Adobe Photoshop. Fiducial points were placed into at least 4 distinctive areas that are visible in the MALDI plate image and served as landmarks for registering the MALDI plate to the robotic spotter. Using these fiducial points, coordinates of the histology annotations are transformed into a system recognized by the spotter. An acoustic robotic microspotter (LabCyte) was used to place crystalline matrix (25 mg/mL sinapinic acid in 1:1 acetonitrile: 0.2% trifluoroacetic acid) spots (∼120 pL) that were typically 180–220 μm in diameter (18). Subsequent to spotting, spot placement was evaluated for accuracy compared to the location designated by the pathologist. Spots that deviated significantly from their desired locations were eliminated from further analysis. Tissue profile spectra were acquired from the spotted areas using an Autoflex II (Bruker Daltonics) MALDI-TOF MS equipped with a SmartBeam Laser (Nd:YAG, 355 nm) Bruker Autoflex II. An accelerating voltage 20 kV was used with an extraction voltage of 18.65 kV and an Einzel lens voltage of 6.00 kV. Delayed extraction (350 nanoseconds) was optimized for resolution at m/z 12,000. Data were collected over the range of 2 to 40 kDa. A total of 400 laser shots were collected for each profile spectrum. Spectra were converted to ASCII files for preprocessing and statistical analysis.

Validation of the identified proteins was performed by Western immunoblotting using surgically resected specimens of squamous cell carcinomas of the lung, as well as adjacent histologically normal tissue obtained from the Vanderbilt lung Specialized Programs of Research Excellence (SPORE) tissue repository. Samples were snap-frozen and stored at −80°C until protein lysates were prepared. Protein concentration of the supernatant was determined by the bicinchoninic acid (BCA) protein assay (Thermo Scientific) and 40 μg of total protein from each sample was separated on 15% SDS-polyacrylamide gels [ubiquitin, macrophage migration inhibitory factor (MIF), thymosin β4, cytochrome c] or 4% to 12% Bis-Tris gels [Invitrogen; acyl-coA binding protein (ACBP) and cystatin A (CSTA)]. Proteins were transferred onto a polyvinylidene diflouride (PVDF) membrane and treated with blocking solution before treatment with primary antibodies for thymosin β4 (1:1,000, Santa Cruz Biotechnology), ubiquitin (1:1,000, R&D Systems), ACBP (1:250, Abcam), CSTA (1:2,500, R&D Systems), cytochrome c (1:1,000, Cell Signaling), and MIF (1:500, R&D Systems). The blots were treated with horseradish peroxidase (HRP)-conjugated secondary antibodies and developed with the enhanced chemiluminescence system (Perkin-Elmer). Immunoblotting of β-actin (Sigma) was used as a loading control.

Five micron tissue sections were cut from FFPE lung tissue blocks, deparaffinized by placing in 3 xylene baths, 10 minutes each, hydrated through a series of ethanol baths of decreasing concentration, and placed in a buffer bath of TBS. Endogenous peroxide was quenched using Dako Peroxidase Blocking Reagent and nonspecific staining blocked using Dako Serum-Free Protein Block according to manufacturer's instructions. Optimal staining for each antibody was obtained by using a citrate antigen retrieval and the following primary antibody dilutions: thymosin β4 (1:500, BioDesign), ubiquitin (1:400, Dako), ACBP (1:80, Abcam), CSTA (1:8,000, Sigma), ProteinTech Group, cytochrome c (1:1,000, Cell Signaling), and MIF (1:200, Santa Cruz Biotechnology). Primary antibodies were applied at room temperature for 1 hour. All further steps and chromogen development were carried out according to manufacturer's instructions (Dako EnVision⁺ System, HRP). Slides were counterstained using a Gills #1 hematoxylin, rinsed appropriately and dehydrated through a series of ethanol baths of increasing concentration into xylene. A nonaqueous coverslipping material was used (Cytoseal or Permount). Results of the IHC staining were analyzed by the pathologist (A.L.G.).

Tissue microarrays (TMA) of lung cancer subtypes from 59 patients were prepared from paraffin blocks. There were 21 squamous carcinomas, 20 adenocarcinomas, 4 bronchioalveolar carcinomas, 4 large cell carcinomas, 2 non-small cell lung cancer (NSCLC), 5 carcinoids, and 1 each of adenosquamous, sarcoma, and small cell lung cancer (SCLC) represented on the TMA. Archived tissue blocks from 1989 to 2002 were retrieved from the files of VUMC and the Nashville VAMC pathology departments. For all tissue blocks, H&E-stained sections were reviewed by our pathologist (A.L.G.). Areas to be punched for array production were carefully marked. Tissues from 26 patients had spots in triplicates, 13 of which also had single spot of corresponding normal tissue obtained at least 2 cm away from the tumor. Tissues from 33 patients had spots in duplicates. Cores 0.6 mm in diameter were taken from the selected area of each specimen and inserted into a recipient paraffin block. Five-micrometer sections were cut from the arrays and mounted onto charged slides. Every 15th section was stained by H&E to confirm the presence of the histological feature of interest (tumor or normal histology). Immunohistochemical staining of the TMA for thymosin β4, ACBP, CSTA, and MIF were was performed according the protocol described earlier. IHC staining was scored 0 to 3 based on intensity by ALG. Scoring for ubiquitin was not performed because of poor staining quality.

Before statistical analysis, the MALDI MS data were preprocessed for calibration, baseline correction, denoising, normalization, peak detection, and quantification by a preprocessing package: Wavespec software (19; The raw data can be found at http://www.vicc.org/biostatistics/supp.php). We selected 9 features from the intercept of best classifiers to distinguish between normal- or low-grade preinvasive lesions and high-grade or invasive cancer from 2 previously described datasets (15, 16). The selection was based on the statistical significance of the discriminatory features (m/z) with false discovery rate-adjusted P < 0.01 and based on the expert visual confirmation of the characteristics of the peak. With these predetermined 9 candidate features, we carried out the following analysis to access their ability to predict whether these bronchial samples were associated with a patient with a lung cancer. First, we generated a weighted compound score (MALDI score) based on the weighted flexible compound covariate method (WFCCM; refs. 16 and 20): the preselected 9 features were reduced to 1 “summarized” feature by building MALDI score for each patient (see Supplementary Materials). For easier interpretation and comparisons, the MALDI scores were rescaled to a range between 0 and 5. Second, a logistic regression model was applied to test the hypothesis that the summarized 9 features profile (MALDI score) and histology group are associated, thus, this MALDI score has predictive ability in differentiating cancer patients from normal controls. The model was derived and internally validated using the Rahman and colleagues cohort (15). A bias-corrected estimate of predictive accuracy was obtained by bootstrap method (21). The class prediction performance was assessed by a widely used measurement of diagnostic discrimination: the area under receiver operating characteristic (ROC) curve and it follows the standards for reporting of diagnosis accuracy (22). Third, we then applied this internally validated model to the new independent validation cohort described earlier by carrying the estimated coefficients of this model to the independent validation set. The prediction ability in terms of area under ROC curve and the predicted probability of having lung cancer with the change of MALDI score were provided for validation set. To complete our analysis, we applied a permutation analysis for an alternative interpretation on the same hypothesis we were testing, and the permutation test reconfirmed such association. From patients who provided more than 1 bronchial biopsy, including biopsies of histological normal and invasive cancer tissues, we used the average of these signal intensities to represent the patient in the analysis (see Supplementary Materials).

Characteristics of the 60 patients in the validation cohort are summarized in Table 1. Characteristics of individual patients of the validation cohort and the Rahman and colleagues cohort (15) are given in Supplementary Tables S1 and S2. Representative images of H&E stained tissue sections of each histological type and images with spots for automatic targeted MALDI matrix deposition are shown in Supplementary Figure S1.

The proteomic signature of 9 features (m/z) was first used to build a prediction model using Rahman and colleagues (15) cohort to assess the probability of having lung cancer. The resultant MALDI score obtained as a weighted sum (see Statistical Section and Supplementary Materials) was used as a continuous variable in the validation cohort. The model showed that the MALDI score in this independent data set is statistically significantly associated with the cancer status. Based on the model, the estimated probability of having lung cancer can be calculated for any patient in the cohort with a MALDI score value. The estimated probabilities for some values of MALDI score are summarized in Table 2. To understand the characteristics of MALDI score better, 2 graphical summaries are provided in Figure 2. The distribution of MALDI scores across patients with and without cancer in the validation cohort (n = 60) is presented in Figure 2A and MALDI score prediction performance through the prediction model, is presented by an area under the ROC curve in Figure 2B. By bootstrap method, the bias-corrected area under the ROC curve is 0.90 [with 95% CI (0.82, 0.98)] in the Rahman and colleagues (15) cohort and 0.77 [95% CI (0.66, 0.88)] in the validation cohort.

Out of 9 features we have identified 8 corresponding to 6 different proteins listed in Table 3. Two of the 6 proteins were newly identified in this study. One of those features at m/z 9,955 was identified ACBP following in-gel digestion and LC-MS/MS analysis (see Supplementary Methods and Supplementary Fig. S2). Based on the features identified, the m/z value was consistent with the loss of the initial methionine residue and an acetylation modification to the N-terminal serine residue. Examination of the 42 Da mass shift in the b-ion series in the MS/MS spectra for the acetylated peptide, SQAEFEK and its miscleavage product SQAEFEKAAEEVR, supports such a modification as shown in Supplementary Figure S3A. Another protein identified was CSTA. Analysis of the gel band corresponding to m/z 11,048 was found to contain an acetylated lysine residue in the peptide sequence NKDDELTGF. The modification was confirmed by examination of the MS/MS spectra containing the 42 Da mass shift in the b-ion series as shown in Supplementary Figure S3B. The unmodified NKDDELTGF peptide was also identified, consistent with the m/z 11,006 observed in the MALDI MS spectrum. All peptides identified for both ACBP and CSTA are listed in Supplementary Table S3. The other 4 proteins identified by our group are thymosin β4 (23), ubiquitin (16), cytochrome c (24), and MIF (25).

Overexpression of each of the 6 identified proteins was then confirmed by Western blotting in lung cancers as shown in Figure 3. Representative data of 3 surgical specimens of squamous cell carcinomas of the lung and adjacent histologically normal lung tissue are presented. Protein expression was further validated in FFPE squamous cell carcinoma of the lung and normal tissues by IHC as shown in Figure 4. Squamous cell carcinomas showed moderate to strong cytoplasmic and nuclear expression of all these proteins. Thymosin β4 and cytochrome c were predominantly expressed in the cytoplasm. In normal bronchial epithelium weak to moderate staining of MIF was observed, while staining was negative to weak for thymosin β4, ubiquitin, ACBP, CSTA, and cytochrome c. We also validated the expression of thymosin β4, ACBP, CSTA, and MIF in a TMA representing lung cancer tissues from 59 patients by IHC. The immunostaining signals of these proteins were localizing to the cytoplasm. The summary of the results (intensity score 0–3) are presented in Supplementary Table S4. Of interest, CSTA staining intensity was associated with prolonged survival (P = 0.03). We found no significant association between staining intensity scores and gender, smoking status, histology, or pathological stage.

In this study we report a proteomic signature derived from the bronchial mucosa that confers a quantifiable probability of having lung cancer. We identified all but 1 of the 9 features of this signature (Table 3). Overexpression of the identified proteins was confirmed by Western blotting (Fig. 3) and immunohistochemical staining (Fig. 4 and Supplementary Table S4) using normal lung and NSCLC. The bronchial biopsy proteomic signature provides a quantitative measure of the probability of having lung cancer that goes beyond the histological evaluation of preinvasive lesions.

Preinvasive lesions are prevalent in the airways of patients at increased risk for and with lung cancer (8–10, 26, 27). Our ultimate goal is to predict the probability of developing lung cancer and to identify the molecular determinants of lung cancer development. In this study, we made an important step toward this goal by providing a quantitative measure of the probability of having lung cancer in individuals presenting with endobronchial lesions (Fig. 2A). Importantly, our diagnostic biomarker signature goes beyond histologic interpretation. It is quantifiable and used as a continuous variable. It was formally validated in an independent data set of specimens from 3 institutions (Supplementary Table S1). The differential expression of these candidate biomarker proteins revealed by proteomic technologies also allows to detect changes in protein expression during lung tumor progression which may otherwise remain undetected by histology alone. Moreover, the development of a proteomic diagnostic signature for lung cancer provides an opportunity to further our understanding about the role of these proteins in various cellular interactions and in pathways of lung tumorigenesis. In Figure 2A, we showed the difference in the MALDI score between patients with and without lung cancer giving the existing overlap. Although such overlap between boxplots is common in the biomarker field, logistical regression analysis indicated that the difference is statistically significant. Yet, these results do not necessarily imply a clinically useful classifier. Further prospective validation of the biomarker is therefore needed to address clinical utility. Finally, this signature may lead to an improved understanding of the progression of preinvasive lesions to an invasive phenotype, potentially enabling earlier therapeutic intervention, and allowing us to further select patients to enroll in surveillance studies.

Scarcity and the small size of preinvasive lesions of the airways remain the major challenges for studies of bronchial mucosa in lung tumor development. To accomplish our goal, we obtained specimens from 3 institutions (VUMC and VAMC in Nashville, TN, and UCHSC in Denver, CO). We also took advantage of significant improvements in proteomic technologies, including automated spotting and related software, MS instrumentation, database resources, and statistical analysis tools (19). Histology directed spotting of matrix using Mirax viewer software, for example, allowed us to limit data acquisition specifically from areas containing greater than 70% of the epithelial cells (18, 28).

Technological advances have also allowed us to confidently identify the proteins represented by our signature. Out of the 6 identified proteins contained in the signature, CSTA and ACBP deserve special attention. Although the role of ACBP in lung tumorigenesis is unknown, ACBP plays a crucial role in lipid metabolism and acts as acyl-CoA transporter and intracellular pool former. In spite of the fact that ACBP is ubiquitously expressed (29), its differential expression in normal epithelium and invasive tissues was clearly demonstrated by MALDI MS. To our knowledge, this is the first report of overexpression of ACBP in lung tumorigenesis.

Cystatins are members of a protein family with endogenous inhibitors of cysteine proteases such as cathepsins B, H, and L. Studies demonstrated an inverse relationship between increased cathepsin B expression and decreased CSTA in a variety of human tumors (30). The importance of CSTA in the development of malignant tumors (31) and its overexpression in lung cancer have been reported (32). Lung tumor cells are suggested to produce both cysteine proteases and cystatins in vitro which are regulated differently in histologic subtypes of lung cancers; however, their diagnostic and prognostic role in large pathological materials has not been studied extensively. One study involving preinvasive bronchial lesions and various lung carcinomas suggests that CSTA expression is associated with malignant transformation and predicts a risk of tumor recurrence as well as poor survival (32). Our studies suggest that CSTA is one of the promising candidate biomarkers for lung cancer and deserves further investigation in lung tumorigenesis. Of interest also is that cathepsin B which is a cysteine protease cleaves 2 terminal glycine residues of ubiquitin generating des-ubiquitin (33). Both monomeric ubiquitin and des-ubiquitin corresponding to m/z values 8,565 and 8,451 are among the 9 candidates of the signature described in this study, and in our previous study in established tumors (16). Overexpression of cathepsin B inhibitor CSTA, and des-ubiquitin, a product of cathepsin B proteolytic activity on ubiquitin, suggests a delicate balance between cathepsin B and CSTA may somehow perturb cellular homeostasis contributing to the overexpression of des-ubiquitin.

MIF is expressed and secreted in response to mitogens and integrin-dependent cell adhesion. Extracellular MIF is required in the steady-state activation of Rho GTPase family members, leading to cell growth and migratory phenotypes. Previous reports demonstrated that MIF expression is increased in premalignant, malignant, and metastatic tumors and a potential therapeutic target in NSCLC (34, 35). By use of a siRNA knockdown of MIF, the migratory and invasive potential of human lung adenocarcinoma cells was abolished (36).

Thymosin β4 sequesters cytoplasmic monomeric actin (37) and is a potent regulator of actin polymerization. Overexpression of thymosin β4 in tumors has been suggested to stimulate lung tumor metastasis by activating cell migration and angiogenesis (38). The gene expression of thymosin β4 was also significantly associated with metastasis in NSCLC patients and it was suggested as a prognostic parameter for NSCLC (39). A combined proteomic indicator including thymosin β4, thymosin β10, and calmodulin was shown to be able to classify NSCLC patients according to good or poor prognosis (40).

Our study has limitations. The patient sample size, especially for noncancer cases, in the training and validation sets was not large. Although the demographics of our study population reflect well individuals who are at high risk for lung cancer, the results may not generalize to those with lower epidemiological risk. Our proteomic signature represents a direct measure of alterations of the airway epithelium in the field of carcinogenesis, but the signature has not yet been prospectively validated as a predictor of lung cancer development. Also, we could not evaluate the possibility that the presence of such a signature represents a true susceptibility factor rather than a tumor marker. Our study also does not address the diagnostic value of the proteomic analysis of alveolar atypical hyperplasia thought to be possible precursors of adenocarcinoma in situ (formerly called bronchioalveolar carcinoma or BAC). Among other study limitations, although the methodology used to assess the probability of having lung cancer may be found useful as an MS-based analysis, the signature may require further development in the form of an immune-based multiplex assay to reach broad clinical utility.

Looking ahead we see this novel signature as a significant advance in the search for early detection strategies. As opposed to histologic evaluation which provides a few discrete subtypes with potential significant interobserver variability (41, 42), the quantitative nature of our signature, after further validation, may be able to better assess the diagnosis of lung cancer from endobronchial biopsies. We will need to further validate the diagnostic potential of the signature prospectively in a cohort of high-risk individuals. Determining the probability of having lung cancer based on a novel proteomic signature of bronchial lesions may have important clinical implications. This signature may facilitate the diagnosis of lung cancer and the monitoring of high-risk individuals for lung cancer in surveillance and chemoprevention trials.

No potential conflicts of interest were disclosed.

We thank all patients who volunteered to participate in this study and provided tissue samples. We thank Candace Murphy and Anna Ostrander for assistance with sample collection, Harriet Davis and Gabe Garcia for clinical data management, Jamie Allen for assistance with sample processing for MALDI MS data acquisition and organization, and Salisha Hill and Kristen Cheek for their assistance with LC-MS/MS.

This work was supported by the National Institute of Health (CA102353 to P.P. Massion, the Lung SPORE P50 CA 90949 to D.P. Carbone, and the Lung SPORE P50 CA58187 to Y.E. Miller and W.A. Franklin), Vanderbilt-Ingram Cancer Center Core Support Grant (P30 CA68485 to P.P. Massion and R.M. Caprioli), Department of Defense (W81XWH-05-1-0179 to R.M. Caprioli) and NIH/NIGMS (5R01-GM58008 to R.M. Caprioli).

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