Proteomic profiling of lung adenocarcinoma indicates heightened DNA repair, antioxidant mechanisms and identifies LASP1 as a potential negative predictor of survival
- Johannes F. Fahrmann†1, 9,
- Dmitry Grapov†2,
- Brett S. Phinney3,
- Carol Stroble4,
- Brian C. DeFelice1,
- William Rom5,
- David R. Gandara4,
- Yanhong Zhang6,
- Oliver Fiehn1, 7,
- Harvey Pass8 and
- Suzanne Miyamoto4Email author
© The Author(s) 2016
Received: 12 July 2016
Accepted: 12 October 2016
Published: 27 October 2016
Lung cancer is the leading cause of cancer mortality in the United States. Non-small cell lung cancer accounts for 85% of all lung cancers for which adenocarcinoma is the most common histological type. Management of lung cancer is hindered by high false-positive rates due to difficulty resolving between benign and malignant tumors. Better molecular analysis comparing malignant and non-malignant tissues will provide important evidence of the underlying biology contributing to tumorigenesis.
We utilized a proteomics approach to analyze 38 malignant and non-malignant paired tissue samples obtained from current or former smokers with early stage (Stage IA/IB) lung adenocarcinoma. Statistical mixed effects modeling and orthogonal partial least squares discriminant analysis were used to identify key cancer-associated perturbations in the adenocarcinoma proteome. Identified proteins were subsequently assessed against clinicopathological variables.
Top cancer-associated protein alterations were characterized by: (1) elevations in APEX1, HYOU1 and PDIA4, indicative of increased DNA repair machinery and heightened anti-oxidant defense mechanisms; (2) increased LRPPRC, STOML2, COPG1 and EPRS, suggesting altered tumor metabolism and inflammation; (3) reductions in SPTB, SPTA1 and ANK1 implying dysregulation of membrane integrity; and (4) decreased SLCA41 suggesting altered pH regulation. Increased protein levels of HYOU1, EPRS and LASP1 in NSCLC adenocarcinoma was independently validated by tissue microarray immunohistochemistry. Immunohistochemistry for HYOU1 and EPRS indicated AUCs of 0.952 and 0.841, respectively, for classifying tissue as malignant. Increased LASP1 correlated with poor overall survival (HR 3.66 per unit increase; CI 1.37–9.78; p = 0.01).
These results reveal distinct proteomic changes associated with early stage lung adenocarcinoma that may be useful prognostic indicators and therapeutic targets.
Lung cancer is a leading cause of cancer mortality in both men and women in the United States [1–3]. Non-small cell lung cancer (NSCLC) accounts for 85% of all lung cancer cases for which NSCLC adenocarcinoma is the most common histological type . While use of low dose computerized tomography (LDCT) for screening of persons at high risk for lung cancer can reduce cancer mortality, it is plagued by high false positive rates (96%)  because it is unable to adequately distinguish indolent (benign) solid pulmonary nodules (SPNs) from malignant SPNs. Increased knowledge of the molecular perturbations caused by tumorigenesis is needed to better understand the underlying biology, as well as potentially assisting with diagnosis, prognosis and identification of additional treatment targets.
Improved ‘Omic’ based analytical methods (e.g. genomics, transcriptomics, proteomics and metabolomics) gives us greater ability to monitor all biochemical processes associated with tumorigenesis with increasingly smaller amounts of difficult-to-obtain clinical specimens. Proteomics is particularly well suited to identify the underlying biology of lung cancer. Protein expression is the ultimate product of gene expression and is controlled through transcriptional, translational and post-translational regulations, all of which are highly dependent on cellular signaling . Unlike genomics (DNA and RNA), proteomic analysis is more complex due to the presence of sequence variations, isoforms, and post-translational modifications, yielding multiple protein isoforms of the same gene Proteomics can uncover powerful links between gene function and tumorigenesis, help discover clinically useful diagnostic and prognostic biomarkers for early stage NSCLC adenocarcinoma [7–9]. In the current study, a shotgun tandem liquid chromatography mass spectrometry (LC–MS/MS) approach was used to characterize proteomic differences between 38 matched malignant and non-malignant lung tissue pairs obtained from current or former smokers with early stage (IA-IB) NSLCL adenocarcinoma. Statistical analysis and multivariate modeling were used to identify the top 10% of all measured protein changes that best distinguished adenocarcinoma from control tissues. Identified proteins were additionally evaluated against clinicopathological variables. Three top candidate proteins, EPRS, HYOU1 and LASP1 from the discovery study were independently validated with a tissue microarray containing 40 pairs of malignant and non-malignant tissues from patients with early stage NSCLC adenocarcinoma. We hypothesize that identification of cancer induced cellular and tissue level protein changes will provide candidate tissue-specific prognostic markers for early stage adenocarcinoma that may eventually be used to better distinguish adenocarcinoma from benign tissues, help identify potential therapeutic targets for treatment of lung cancer and, importantly, improve our understanding of the mechanism(s) leading to lung cancer.
De-identified malignant and adjacent non-malignant lung tissue samples were harvested in the operating room from patients having resection or lobectomy for non-small cell lung cancer, none of whom received preoperative treatment. The matching control lung tissue was always taken from areas 8–10 cm removed from the cancer bed. All procedures were approved by institutional IRB protocols (NYU) with patient consent. Tissues were immediately frozen in liquid nitrogen and stored at −80 °C. Specimens were annotated for age, gender, race, diagnosis (including stage), smoking status and pack-years. Criteria used to select patient cases were: (a) current or former smokers; (b) diagnosis of NSCLC adenocarcinoma; (c) early stage IA or IB; and (d) understood and signed the IRB consent form. Of the 38 samples analyzed in the study, 14 (36.8%) were from patients that had recurrent cancer within 5 years. The inclusion of these 14 subjects was to identify potential prognostic biomarkers of early stage adenocarcinoma.
Preparation of samples for proteomic analysis
Preparation of tissues and N-glycan release followed by ethanol precipitation of tissue proteins has been previously described in Ruhaak et al. .
Trypsin digestion of samples
Protein pellets were solubilized in 100 µL of AMBIC (6 M urea, 50 mM ammonium bicarbonate) with dithiothreitol added (final concentration 5 mM) and incubated at 37 °C for 30 min. Iodoacetamide (IAA) was added (final concentration 15 mM), incubated for 30 min (RT) and DTT added to quench the IAA reaction. Lys-C/trypsin (Promega) was next added (1:25 enzyme:protein ratio), samples incubated (4 h at 37 °C) after which the urea concentration was lowered (<1 M) and samples further digested (overnight, 37 °C). Samples were desalted (C18 Macro Spin columns, Nest Group) and lyophilized.
Reconstituted samples were randomized into 8 blocks of 4 normal and 4 tumor samples. Triplicate LC–MS/MS analyses were acquired for each sample. LC separation was done on a Waters Nano Acquity UHPLC (Waters Corporation) with a Proxeon nanospray source. Mass spectra was collected on an Orbitrap Q Exactive Plus mass spectrometer (Thermo Fisher Scientific) in a data-dependent mode with one MS precursor scan followed by 15 MS/MS scans as previously described . Detailed information on instrument parameters and mass spectra collection is provided in the Additional file 1: Supplemental materials and methods (1a). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium  via the PRIDE partner repository with the dataset identifier PXD002612.
Immunohistochemistry of lung cancer LC003 TMA
Preparation of LC003 TMA
A tissue microarray (LC003) was prepared from FFPE tissue blocks containing non-malignant and malignant (tumor) tissues obtained from 40 patients diagnosed with early stage (Stage I and II) NSCLC adenocarcinoma after patient consent using an IRB approved protocol (IRB 293828, UC Davis Cancer Center Biorepository). Three 0.6 mm (diameter) by 0.4 mm (length) cylindrical cores of tumor and adjacent non-malignant (control) lung tissues were collected from each case and placed in the same block (Quick-Ray Manual Tissue Microarrayer). One TMA section was stained with H&E (hematoxylin and eosin) to confirm the presence of tumor and normal lung tissue. Clinical data (gender, age, smoking history, diagnosis and stage) information was provided after analysis.
Immunohistochemistry of LC003 TMA
Immunohistochemical staining was performed on 4-μm FFPE tissue sections from the LC003 TMA for three different protein targets, HYOU1, EPRS and LASP1. Details of the IHC procedure are provided in Additional file 1: Supplemental materials and methods (1c).
Immunostained slides were scored by a pathologist (YZ) blinded to clinicopathologic information. Staining for HYOU, EPRS and LASP1 was semi-quantified assessed using both intensity and percentage of positive cells. Staining intensity was graded as 0 = negative (no cells stained); 1 = weak; 2 = moderate; 3 = strong. The H score was calculated by multiplying the percentage of positive cells demonstrating each density (score 0 to 3) and adding the results. The average H score of 3 cores in a case was calculated. High expression was defined as H score 251–300; medium expression was H score 151–250; low expression was H score less than 150.
Detailed information on statistical analyses are provided in Additional file 1: Supplemental Materials and methods (1d). A conservative approach was used to overview data quality metrics in order to identify the most robust proteomic measurements for downstream statistical and multivariate analyses . A total of 799 high quality protein measurements (spectral counts) based on absence of abundant missing zeros, replicate precision, protein sequence coverage and annotation to known genes were selected for further statistical and multivariate analyses.
Statistical analyses were conducted on mean normalized count data. Statistical tests were conducted using generalized negative binomial mixed effects regression with the patients as the random term . Significantly altered proteins were identified based on the comparison of the full model (Chi squared test) including age + gender + packs + tumor/control to a reduced model excluding the tumor/normal labels. Model p-values were adjusted for multiple hypotheses tested .
Orthogonal partial least squares discriminant analysis (O-PLS-DA) multivariate models were used to identify the top 10% of all protein discriminants between tumor and control tissues as previously described .
A Gaussian graphical model protein–protein empirical network was calculated for O-PLS-DA selected top discriminants (n = 16) as previously described .
All subjects were included in the survival analysis with the exception of six subjects who did not die of lung cancer 23–2262 days from their surgery, and who, before their death, never had recurrent lung cancer. Kaplan–Meier survival curves were generated using Prism v5.0 (GraphPad Software, Inc). Significance was determined using log-rank Mantel-Cox test. Cox proportional hazard models were carried out in R statistical software.
LASP1 transcriptomic data
LASP1 transcriptomic data was obtained from the Okayama NSCLC study  using the Oncomine Database . The study set was chosen due its specific focus on early stage NSCLC adenocarcinoma and availability of clinicopathological variables. Only Stage I (IA/IB) adenocarcinoma subjects were considered. Subject characteristics are provided in Additional file 1: Table S5. LASP1 values were log2-median centered normalized. Cox proportional hazard models were used to evaluate the association between LASP1 mRNA expression and overall survival.
Lung cancer patients
Total sample size, N
Gender, N (males/female)
Age, mean ± SD (min, max)
70 ± 9 (56, 91)
Packs per year, mean ± SD (min, max)
33 ± 25 (0, 100)
Current smoker, N (%)
Stage IA, N (%)
Stage IB, N (%)
Developed lymphovascular invasion, N (%)
Developed pleural invasion, N (%)
Progressed, N (%)
Second primary, N
Top 10% discriminants of adenocarcinoma
Control (mean ± SD)a
Tumor (mean ± SD)a
Spectrin, beta, erythrocytic
10.7 ± 9.3
2.82 ± 2.6
Solute carrier family 4 (anion exchanger), member 1 (Diego blood group)
8.75 ± 5.2
3.8 ± 2.7
Spectrin, alpha, erythrocytic 1
13.3 ± 12
4.28 ± 3.4
Ankyrin 1, erythrocytic
6.28 ± 5.7
1.84 ± 1.8
1.24 ± 0.87
3.15 ± 1.7
Hemoglobin, gamma A
16.9 ± 9.7
7.44 ± 6.2
Hemoglobin, gamma G
17 ± 9.7
7.43 ± 6.2
Coatomer protein complex, subunit gamma 1
1.56 ± 1
2.67 ± 1.4
Hypoxia up-regulated 1
2.12 ± 1.9
4.74 ± 2.3
Polymerase I and transcript release factor
7.35 ± 2.8
3.81 ± 2.8
Protein disulfide isomerase family A, member 4
2.69 ± 2.2
6.23 ± 3.4
APEX nuclease (multifunctional DNA repair enzyme) 1
0.778 ± 0.83
2.05 ± 1.1
Stomatin (EPB72)-like 2
0.675 ± 0.39
1.38 ± 0.57
Leucine-rich pentatricopeptide repeat containing
1.31 ± 1.4
3.52 ± 2.5
N-acetylneuraminic acid synthase
0.61 ± 0.55
1.64 ± 1.1
47.3 ± 32
21.7 ± 18
NSCLC adenocarcinoma is characterized by alterations in DNA repair mechanisms, antioxidant defense capacity, altered membrane integrity and metabolism
SPTB was determined to be the single most discriminatory protein of adenocarcinoma displaying a 70% reduction in tumor tissue relative to control tissue (Table 2). The adenocarcinoma-dependent reduction in SPTB was also consistently observed in 82% of subjects. Tumor-associated reductions in SPTB were associated with similar reductions in SPTA1 (70%), SLC4A1 (60%), and ANK1 (70%) in tumor tissue when compared to control tissue (Fig. 1; Table 2). The reduction in SPTA1 was also directly correlated with similar reductions in the hemoglobin subunits HBD, HBG1 and HBG2, which were also found to be directly correlated among each other (Fig. 1). Reductions in SPTA1 and SPTB were indirectly associated with APEX1, which showed a 2.7-fold increase in adenocarcinoma compared to non-malignant tissue and was consistently elevated in 82% of subjects (Fig. 1; Table 2). The increase in APEX1 was correlated with similar increases in HYOU1 and NANS, which were also increased 2.2- and 2.7-fold in malignant compared to non-malignant tissue (Fig. 1; Table 2). HYOU1 and NANS were both positively correlated with adenocarcinoma-associated increases in EPRS, which was also positively associated with LRPPRC and COPG1 (Fig. 1). EPRS was repeatedly elevated in adenocarcinoma relative to control in 79% of cancer subjects, whereas both LRPPRC and COPG1 were consistently elevated in 76% of subjects (Table 2). Adenocarcinoma-dependent elevations in LRPPRC were associated with similar increases in STOML2 and PDIA4, which exhibited 2- and 2.3-fold increases in adenocarcinoma relative to control tissue, respectively (Fig. 1; Table 2). NANS, COPG1 and PDIA4 were all negatively associated with adenocarcinoma-dependent reductions in PTRF, which was generally decreased in adenocarcinoma compared to control tissue (89% of patients).
mRNA expression data from the Okayama et al. study on NSCLC adenocarcinoma  was used to strengthen the adenocarcinoma-associated proteomic perturbations. Only Stage I NSCLC adenocarcinomas were considered. Consistent with our proteomic data, mRNA expression of APEX1, HYOU1, PDIA4, NANS, LRPPRC, EPRS and COPG1 were significantly (Mann–Whitney U < 0.05) higher in adenocarcinoma compared to control whereas mRNA abundance of HBG1, HBG2, HBD, and PTRF were significantly (Mann–Whitney U < 0.05) lower (Additional file 1: Figure S1). No significant differences were observed for SCL4A1, SPTB, SPTA1 and ANK1. STOML2 was not detected in the Okayama Lung dataset.
Validation of HYOU1 and EPRS by immunohistochemistry of LC003 NSCLC adenocarcinoma TMA
Immunohistochemistry of a tissue microarrays (LC003 TMA) for NSCLC adenocarcinoma and matched controls (n = 40) was used to validate our proteomic findings. We focused on EPRS and HYOU1 due to both proteins being in the top 10% features that distinguish malignant from control tissue and due to their positive staining based on the Cancer Protein Atlas. All TMA results were blindly scored using s manual “H” scoring system (“Methods”).
Association between proteomic signatures and overall survival
Hazard ratios for LASP1 mRNA abundance and overall survival in Okayama Dataset
≤61 years of age
>61 years of age
In the current study, we evaluated the proteome of 38 malignant and matched control tissue of stage IA and IB lung adenocarcinoma. Differential analysis identified 436 differentially expressed proteins in adenocarcinoma compared to control tissue, of which 367 remained significant following false-discovery rate adjustment. Orthogonal partial least squares discriminant analysis identified the top 10% proteins that significantly differed between adenocarcinoma and control tissue. A Gaussian graphical model network was used to identify conditionally independent empirical protein–protein relationships between O-PLS-DA selected discriminants for adenocarcinoma. Of the 16 distinguishing proteins, 8 were significantly elevated in lung adenocarcinoma (COPG1, STOML2, HYOU1, PDIA4, EPRS, APEX1, LRPPRC and NANS) whereas 8 proteins were significantly decreased in lung adenocarcinoma (SPTB, SPTA1, ANK1, SLC4A1, HBG1 and HBG2). A subsequent sub-analysis was conducted to further evaluate the prognostic capacity of using all identified proteins.
Top discriminatory proteins of adenocarcinoma indicate heightened intrinsic defense mechanisms, altered metabolism and perturbed membrane integrity
NANS (N-acetylneuraminic acid synthase) exhibited the largest increased (2.7-fold) in lung adenocarcinoma relative to control and was consistently elevated in 82% of the subjects. The elevation in NANs suggests increased biosynthesis of sialic acid, which is known to be elevated in many cancers and may exhibit immune-modulatory and anti-apoptotic functions . Elevation in NANS was positively correlated with adenocarcinoma-associated elevations in APEX1 (apurinicapyrimidinic endonuclease), an important component of base excision repair (BER) pathway and transcriptional modulator of genes that protect against oxidative stress . Over expression of APEX1 has been described in NSCLC and other cancers , while its down-regulation may induce apoptosis . The elevation in APEX1 is particularly interesting since radiotherapy/chemotherapy induce DNA damage and promote production of cytotoxic reactive oxygen species (ROS) . An increase in repair machinery and anti-oxidant defense systems within a tumor cell could interfere with these treatments. APEX1 has been linked to chemotherapy/radiotherapy resistance .
The elevation in APEX1 was paralleled by an increase in HYOU1, a heat shock protein, which has important roles in hypoxia and angiogenesis and is linked to tumor prognosis . Increased HYOU1 also positively correlated with PDIA4, a known modulator of redox status, which has been shown to promote drug resistance to cisplatin in lung adenocarcinoma [28, 29]. PDIA4 also directly correlated to LRPPRC. LRPPRC acts as a regulator of mitochondrial DNA-encoded mRNAs and participates in glucose homeostasis, energy metabolism and nuclear receptor activation [30, 31] and is abundantly expressed in NSCLC adenocarcinoma and other cancers . Knockdown of LRPPRC promotes apoptosis and reduces tumor invasiveness . Elevation in LRPPR was paralleled by increases in STOML2 and COPG1. STOML2, a member of the stomatin family, is upregulated in numerous cancers and linked to tumor aggressiveness . COPG1 is a component of the COPI complex, an integral component in lipid homeostasis that promotes lipolysis through the association of PNPLA2, a triglyceride lipase, with lipid droplets . Lipid metabolism is known to be perturbed in various malignancies acting as energy sources, signaling moieties and for biosynthesis of structural membrane lipids facilitating increased proliferative potential . Cancer cells not only rely on de novo fatty acid synthesis, but also take fatty acids from their surrounding microenvironment . We previously reported that malignant cells exhibit perturbations in free fatty acid profiles with most fatty acids being decreased in tumorous tissue, with the exception of arachidonic acid that was elevated in tumor compared to non-tumor tissue . Free fatty acids, particularly arachidonic acid, serve as substrates for generation of inflammatory mediators. This may also provide a basis for the tumor-associated elevation in EPRS, which was positively correlated with LRPPRC and COPG1. EPRS (glutamyl-prolyl-tRNA synthetase) is a bifunctional enzyme thought to be a gatekeeper of inflammatory gene translation, modulator of angiogenesis and regulator of amino-acid stress responses . EPRS is also regulated by c-Myc, which is commonly amplified in lung adenocarcinoma and modulates tumor metabolism [36, 37]. Collectively, increases in LRPRRC, COPG1, EPRS and STOML2 all point towards alterations in energetics, metabolism and alterations in inflammatory responses that accompany transformation.
Adenocarcinoma-associated reductions in spectrins, SPTB and SPTA1, and ankrin (ANK1) all interact with each another to regulate cell shape and membrane integrity, so paralleled changes likely reflect changes in cell adherence, which is a known hallmark of metastatic cancer . Alterations in SPTB and ANK1 also indicate changes in rearrangement of transmembrane proteins, including ion channels, which may account for the observed positive correlation between SPTB and SLC4A1, a membrane bound anion exchange transporter. SLC4A1 plays central roles in pH homeostasis and has been linked to tumor aggressiveness in numerous cancer types . Since most cancer cells exhibit a metabolic shift towards acidic-producing pathways, reflective of both oncogenic signaling and hypoxia, upregulation of pH-regulatory transport proteins may be important . An alteration in SLC4A1 represents an attractive target for therapeutic intervention and may provide diagnostic value.
LASP1 as a potential negative prognostic indicator for overall survival
LIM and SH3 domain protein (LASP1) a dynamic protein involved in cell structure, physiological processes and cell signaling, is significantly expressed in various malignancies and associated with tumor aggressiveness  and is reported to be an independent prognostic factor in patient’s survival for gastric cancer  and hepatocellular carcinoma . We also found that LASP1 was associated with poor overall survival in our cohort. When considered as a continuous variable, LASP1 indicated a HR of 3.66 (CI 1.37–9.78; p 0.01) per unit increase for increased risk of death. Notably, adenocarcinoma abundance of LASP1 was only significantly higher relative to control tissue in those subjects who died (Additional file 1: Table S4). To strengthen the association between LASP1 and overall survival, we utilized a second independent dataset on LASP1 mRNA expression in NSCLC adenocarcinoma . LASP1 mRNA expression was a significant negative predictor of overall survival. LASP1 has been shown to promote invasion and metastasis; particularly, LASP1 overexpression in pancreatic ductal adenocarcinomas was found to be significantly associated with lymph node metastasis and poor overall survival . Therefore, it is plausible that tumor-associated elevations in LASP1 may contribute to LVI, which could explain why we no longer observed a significant hazard ratio when a LASP1 cutoff of 2.1 was used and instead only LVI was found to be a significant predictor of overall survival. Collectively, these findings highlight the potential of LASP1 as a prognostic indicator for NSCLC adenocarcinoma. Further studies in larger cohorts are required to fully validate these findings.
One limitation of the current study, and others like it, was the lack of assessment of tissue microheterogeneity at the sub-biopsy level. We acknowledge that changes in specific proteins, such as hemoglobin subunits, may reflect contamination from red blood cells. However, changes in erythrocyte-associated proteins may equally reflect altered tumor angiogenesis . Chen et al. used 2D-PAGE, MALDI MS to analyze Stage I and III lung adenocarcinomas and found proteins that were heavily implicated in antioxidant response systems and cellular metabolism  consistent with our findings. Zhou et al. used 2D-DIGE with followed-up by TMA IHC and blood studies of TyrRS and MACF-1 . However, Zhou et al. evaluated tissues from stage II and III lung adenocarcinoma and thus might not be directly compared to our results for stage IA or IB given the inherent heterogeneity that exists among different tumor stages. Kikuchi et al. conducted in-depth proteomic profiling on lung adenocarcinomas and squamous cell carcinomas and compared these against control tissue . Similar to our findings, Kikuchi et al. also found tumor-associated elevations in numerous mentioned above proteins including HYOU1, NANS and PDIA4 .
In conclusion, proteomic changes in early stage NSCLC adenocarcinoma tissues are consistent with known cancer-dependent alterations in repair machinery, redox status, energetics, and inflammation. The current study identified candidate markers that may help identify at-risk subjects and assist with treatment. This study also suggests that LASP1 might serve as a potential negative prognostic marker for overall survival. Further studies in larger cohorts are warranted to confirm and validate these findings.
JFF: wrote manuscript, interpreted data and generated figs/tables. DG: wrote manuscript, interpreted data, performed statistical analyses. BSP: generated proteomic data, provide expertise and knowledge of proteomic findings, aided in study design, critically reviewed manuscript. CS: perform sample preparation and provided knowledge. BCF: generated proteomic data. WNM: provided knowledge and expertise about lung cancer, aided in study design and critically reviewed manuscript. DRG: provided knowledge and expertise about lung cancer, provided funding and critically reviewed manuscript. YZ: analyzed the TMA IHC results and generated IHC data. She helped with the interpretation of the TMA results and provided the images for the TMA IHC figures. OF: generated metabolomic data, provided expertise and knowledge of metabolomic findings, aided in study design, critically reviewed manuscript. HP: helped lead the study, aided in study design, provided samples, provided knowledge and expertise, critically reviewed manuscript. SM: led the study, conceived study design, secured funding to support the study, provided knowledge and expertise, helped write the manuscript, critically reviewed manuscript.
We acknowledge the LUNGevity and Thomas G. Labrecque Foundations award 201118739 (Miyamoto), NIH 1 U24 DK097154 NIH West Coast Metabolomics Center award (Fiehn), and NCI EDRN 5U01CA111295-07 (Pass). We acknowledge the efforts of Chandra Goparaju and Ki Prokrym (from NYU) who assisted with patient sample collections; Camille Rodriguez (UC Davis) for assistance with sample distribution; Regina Gandour Edwards (MD) and Irmgard Feldman (Pathology, UC Davis) from the UC Davis Cancer Center Biorepository for helping with the acquisition of the lung cancer tissue microarray (LC003) that was used for validation of protein targets in this study. We also acknowledge the administrative support of Audrey Sorensen (NYU) whose efforts on this project were greatly appreciated.
The authors declare that they have no competing interests.
Availability of data and material
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
Lung cancer (malignant) and control frozen tissues were obtained from the NYU Biorepository (Harvey Pass, MD) following an approved NYU institution IRB protocol with patient consent (Study Title: NYU Lung Cancer Biomarker Center, Study# i8896, Jun-Chieh Tsay, MD, PI). The Lung Cancer Tissue Microarray (TMA 003) was prepared by the UC Davis Cancer Center Biorepository using archival tissue blocks obtained from the UC Davis Cancer Center Biorepository using IRB approved protocol #218204-24, “Cancer Center Biorepository (CCB) - Tissue, Blood, Urine and Other Biological Materials”, Regina Gandour Edwards, MD, Pathology, Principal Investigator.
We acknowledge the LUNGevity and Thomas G. Labrecque Foundations award 201118739 (Miyamoto), who provided funding for the study, NIH 1 U24 DK097154 NIH West Coast Metabolomics Center award (Fiehn) supported the equipment, resources and personnel to staff the WCMC, and NCI EDRN 5U01CA111295-07 (Pass), who supported the acquisition of the patient tissue samples.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Prevention CfDCa. National Center for Health Statistics. CDC WONDER On-line Database, compiled from Compressed Mortality File 1999–2012; 2014, Series 20 No. 2R.Google Scholar
- Society AC. Cancer facts and figures (http://www.cancer.org/research/cancerfactsstatistics/cancerfactsfigures2014/index). 2014. p. 1–60.
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2015. CA Cancer J Clin. 2015;65:5–29.View ArticlePubMedGoogle Scholar
- Breathnach OS, Freidlin B, Conley B, Green MR, Johnson DH, Gandara DR, O’Connell M, Shepherd FA, Johnson BE. Twenty-two years of phase III trials for patients with advanced non-small-cell lung cancer: sobering results. J Clin Oncol. 2001;19:1734–42.PubMedGoogle Scholar
- National Lung Screening Trial Research Team, Aberle DR, Adams AM, Berg CD, Black WC, Clapp JD, Fagerstrom RM, Gareen IF, Gatsonis C, Marcus PM, Sicks JD. Reduced lung-cancer mortality with low-dose computed tomographic screening. N Engl J Med. 2011;365:395–409.View ArticleGoogle Scholar
- van der Wekken AJ, Hiltermann TJ, Groen HJ. The value of proteomics in lung cancer. Ann Transl Med. 2015;3:29.PubMedPubMed CentralGoogle Scholar
- Hasan N, Kumar R, Kavuru MS. Lung cancer screening beyond low-dose computed tomography: the role of novel biomarkers. Lung. 2014;192:639–48.View ArticlePubMedGoogle Scholar
- Patz EF Jr, Campa MJ, Gottlin EB, Kusmartseva I, Guan XR, Herndon JE 2nd. Panel of serum biomarkers for the diagnosis of lung cancer. J Clin Oncol. 2007;25:5578–83.View ArticlePubMedGoogle Scholar
- Li XJ, Hayward C, Fong PY, Dominguez M, Hunsucker SW, Lee LW, McLean M, Law S, Butler H, Schirm M, et al. A blood-based proteomic classifier for the molecular characterization of pulmonary nodules. Sci Transl Med. 2013;5:207ra142.PubMedPubMed CentralGoogle Scholar
- Ruhaak LR, Taylor SL, Stroble C, Nguyen UT, Parker EA, Song T, Lebrilla CB, Rom WN, Pass H, Kim K, et al. Differential N-glycosylation patterns in lung adenocarcinoma tissue. J Proteome Res. 2015;14:4538–49.View ArticlePubMedGoogle Scholar
- Grapov D, Lemay DG, Weber D, Phinney BS, Azulay Chertok IR, Gho DS, German JB, Smilowitz JT. The human colostrum whey proteome is altered in gestational diabetes mellitus. J Proteome Res. 2015;14:512–20.View ArticlePubMedGoogle Scholar
- Vizcaino JA, Deutsch EW, Wang R, Csordas A, Reisinger F, Rios D, Dianes JA, Sun Z, Farrah T, Bandeira N, et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat Biotechnol. 2014;32:223–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang B, Wang J, Wang X, Zhu J, Liu Q, Shi Z, Chambers MC, Zimmerman LJ, Shaddox KF, Kim S, et al. Proteogenomic characterization of human colon and rectal cancer. Nature. 2014;513:382–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu J, McCleland M, Stawiski EW, Gnad F, Mayba O, Haverty PM, Durinck S, Chen YJ, Klijn C, Jhunjhunwala S, et al. Integrated exome and transcriptome sequencing reveals ZAK isoform usage in gastric cancer. Nat Commun. 2014;5:3830.PubMedPubMed CentralGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate—a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995; 289–300.Google Scholar
- Wikoff WR, Grapov D, Fahrmann JF, DeFelice B, Rom WN, Pass HI, Kim K, Nguyen U, Taylor SL, Gandara DR, et al. Metabolomic markers of altered nucleotide metabolism in early stage adenocarcinoma. Cancer Prev Res (Phila). 2015;8:410–8.View ArticleGoogle Scholar
- Okayama H, Kohno T, Ishii Y, Shimada Y, Shiraishi K, Iwakawa R, Furuta K, Tsuta K, Shibata T, Yamamoto S, et al. Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 2012;72:100–11.View ArticlePubMedGoogle Scholar
- Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6:1–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Higgins KA, Chino JP, Ready N, D’Amico TA, Berry MF, Sporn T, Boyd J, Kelsey CR. Lymphovascular invasion in non-small-cell lung cancer: implications for staging and adjuvant therapy. J Thorac Oncol. 2012;7:1141–7.View ArticlePubMedGoogle Scholar
- Al-Alao BS, Gately K, Nicholson S, McGovern E, Young VK, O’Byrne KJ. Prognostic impact of vascular and lymphovascular invasion in early lung cancer. Asian Cardiovasc Thorac Ann. 2014;22:55–64.View ArticlePubMedGoogle Scholar
- Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10:7252–9.View ArticlePubMedGoogle Scholar
- Henry-Mowatt J, Dive C, Martinou JC, James D. Role of mitochondrial membrane permeabilization in apoptosis and cancer. Oncogene. 2004;23:2850–60.View ArticlePubMedGoogle Scholar
- Yoo DG, Song YJ, Cho EJ, Lee SK, Park JB, Yu JH, Lim SP, Kim JM, Jeon BH. Alteration of APE1/ref-1 expression in non-small cell lung cancer: the implications of impaired extracellular superoxide dismutase and catalase antioxidant systems. Lung Cancer. 2008;60:277–84.View ArticlePubMedGoogle Scholar
- Zaky A, Busso C, Izumi T, Chattopadhyay R, Bassiouny A, Mitra S, Bhakat KK. Regulation of the human AP-endonuclease (APE1/Ref-1) expression by the tumor suppressor p53 in response to DNA damage. Nucleic Acids Res. 2008;36:1555–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Malhotra V, Perry MC. Classical chemotherapy: mechanisms, toxicities and the therapeutic window. Cancer Biol Ther. 2003;2:S2–4.View ArticlePubMedGoogle Scholar
- Gossage L, Perry C, Abbotts R, Madhusudan S. Base excision repair factors are promising prognostic and predictive markers in cancer. Curr Mol Pharmacol. 2012;5:115–24.View ArticlePubMedGoogle Scholar
- Stojadinovic A, Hooke JA, Shriver CD, Nissan A, Kovatich AJ, Kao TC, Ponniah S, Peoples GE, Moroni M. HYOU1/Orp150 expression in breast cancer. Med Sci Monit. 2007;13:231–9.Google Scholar
- Laurindo FR, Pescatore LA, Fernandes Dde C. Protein disulfide isomerase in redox cell signaling and homeostasis. Free Radic Biol Med. 2012;52:1954–69.View ArticlePubMedGoogle Scholar
- Tufo G, Jones AW, Wang Z, Hamelin J, Tajeddine N, Esposti DD, Martel C, Boursier C, Gallerne C, Migdal C, et al. The protein disulfide isomerases PDIA4 and PDIA6 mediate resistance to cisplatin-induced cell death in lung adenocarcinoma. Cell Death Differ. 2014;21:685–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Tian T, Ikeda J, Wang Y, Mamat S, Luo W, Aozasa K, Morii E. Role of leucine-rich pentatricopeptide repeat motif-containing protein (LRPPRC) for anti-apoptosis and tumourigenesis in cancers. Eur J Cancer. 2012;48:2462–73.View ArticlePubMedGoogle Scholar
- Cooper MP, Qu L, Rohas LM, Lin J, Yang W, Erdjument-Bromage H, Tempst P, Spiegelman BM. Defects in energy homeostasis in Leigh syndrome French Canadian variant through PGC-1alpha/LRP130 complex. Genes Dev. 2006;20:2996–3009.View ArticlePubMedPubMed CentralGoogle Scholar
- Song L, Liu L, Wu Z, Lin C, Dai T, Yu C, Wang X, Wu J, Li M, Li J. Knockdown of stomatin-like protein 2 (STOML2) reduces the invasive ability of glioma cells through inhibition of the NF-kappaB/MMP-9 pathway. J Pathol. 2012;226:534–43.View ArticlePubMedGoogle Scholar
- Beller M, Sztalryd C, Southall N, Bell M, Jackle H, Auld DS, Oliver B. COPI complex is a regulator of lipid homeostasis. PLoS Biol. 2008;6:e292.View ArticlePubMedPubMed CentralGoogle Scholar
- Currie E, Schulze A, Zechner R, Walther TC, Farese RV Jr. Cellular fatty acid metabolism and cancer. Cell Metab. 2013;18:153–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Keller TL, Zocco D, Sundrud MS, Hendrick M, Edenius M, Yum J, Kim YJ, Lee HK, Cortese JF, Wirth DF, et al. Halofuginone and other febrifugine derivatives inhibit prolyl-tRNA synthetase. Nat Chem Biol. 2012;8:311–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Seo AN, Yang JM, Kim H, Jheon S, Kim K, Lee CT, Jin Y, Yun S, Chung JH, Paik JH. Clinicopathologic and prognostic significance of c-MYC copy number gain in lung adenocarcinomas. Br J Cancer. 2014;110:2688–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller DM, Thomas SD, Islam A, Muench D, Sedoris K. c-Myc and cancer metabolism. Clin Cancer Res. 2012;18:5546–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Bennett V, Healy J. Membrane domains based on ankyrin and spectrin associated with cell-cell interactions. Cold Spring Harb Perspect Biol. 2009;1:a003012.View ArticlePubMedPubMed CentralGoogle Scholar
- Gorbatenko A, Olesen CW, Boedtkjer E, Pedersen SF. Regulation and roles of bicarbonate transporters in cancer. Front Physiol. 2014;5:130.View ArticlePubMedPubMed CentralGoogle Scholar
- Orth MF, Cazes A, Butt E, Grunewald TG. An update on the LIM and SH3 domain protein 1 (LASP1): a versatile structural, signaling, and biomarker protein. Oncotarget. 2015;6:26–42.PubMedGoogle Scholar
- Zheng J, Yu S, Qiao Y, Zhang H, Liang S, Wang H, Liu Y, Zhou F, Jiang J, Lu S. LASP-1 promotes tumor proliferation and metastasis and is an independent unfavorable prognostic factor in gastric cancer. J Cancer Res Clin Oncol. 2014;140:1891–9.View ArticlePubMedGoogle Scholar
- Wang H, Li W, Jin X, Cui S, Zhao L. LIM and SH3 protein 1, a promoter of cell proliferation and migration, is a novel independent prognostic indicator in hepatocellular carcinoma. Eur J Cancer. 2013;49:974–83.View ArticlePubMedGoogle Scholar
- Zhao T, Ren H, Li J, Chen J, Zhang H, Xin W, Sun Y, Sun L, Yang Y, Sun J, et al. LASP1 is a HIF1alpha target gene critical for metastasis of pancreatic cancer. Cancer Res. 2015;75:111–9.View ArticlePubMedGoogle Scholar
- Lind JS, Smit EF. Angiogenesis inhibitors in the treatment of non-small cell lung cancer. Ther Adv Med Oncol. 2009;1:95–107.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen G, Gharib TG, Huang CC, Thomas DG, Shedden KA, Taylor JM, Kardia SL, Misek DE, Giordano TJ, Iannettoni MD, et al. Proteomic analysis of lung adenocarcinoma: identification of a highly expressed set of proteins in tumors. Clin Cancer Res. 2002;8:2298–305.PubMedGoogle Scholar
- Zhou X, Xue L, Hao L, Liu S, Zhou F, Xiong H, Qi X, Lin D, Shao S. Proteomics-based identification of tumor relevant proteins in lung adenocarcinoma. Biomed Pharmacother. 2013;67:621–7.View ArticlePubMedGoogle Scholar
- Kikuchi T, Hassanein M, Amann JM, Liu Q, Slebos RJ, Rahman SM, Kaufman JM, Zhang X, Hoeksema MD, Harris BK, et al. In-depth proteomic analysis of nonsmall cell lung cancer to discover molecular targets and candidate biomarkers. Mol Cell Proteomics. 2012;11:916–32.View ArticlePubMedPubMed CentralGoogle Scholar