Quantitative Tyrosine Phosphoproteomics of Epidermal Growth Factor Receptor (EGFR) Tyrosine Kinase Inhibitor-treated Lung Adenocarcinoma Cells Reveals Potential Novel Biomarkers of Therapeutic Response

Mutations in the Epidermal growth factor receptor (EGFR) kinase domain, such as the L858R missense mutation and deletions spanning the conserved sequence ⁷⁴⁷LREA⁷⁵⁰, are sensitive to tyrosine kinase inhibitors (TKIs). The gatekeeper site residue mutation, T790M accounts for around 60% of acquired resistance to EGFR TKIs. The first generation EGFR TKIs, erlotinib and gefitinib, and the second generation inhibitor, afatinib are FDA approved for initial treatment of EGFR mutated lung adenocarcinoma. The predominant biomarker of EGFR TKI responsiveness is the presence of EGFR TKI-sensitizing mutations. However, 30-40% of patients with EGFR mutations exhibit primary resistance to these TKIs, underscoring the unmet need of identifying additional biomarkers of treatment response. Here, we sought to characterize the dynamics of tyrosine phosphorylation upon EGFR TKI treatment of mutant EGFR-driven human lung adenocarcinoma cell lines with varying sensitivity to EGFR TKIs, erlotinib and afatinib. We employed stable isotope labeling with amino acids in cell culture (SILAC)-based quantitative mass spectrometry to identify and quantify tyrosine phosphorylated peptides. The proportion of tyrosine phosphorylated sites that had reduced phosphorylation upon erlotinib or afatinib treatment correlated with the degree of TKI-sensitivity. Afatinib, an irreversible EGFR TKI, more effectively inhibited tyrosine phosphorylation of a majority of the substrates. The phosphosites with phosphorylation SILAC ratios that correlated with the TKI-sensitivity of the cell lines include sites on kinases, such as EGFR-Y1197 and MAPK7-Y221, and adaptor proteins, such as SHC1-Y349/350, ERRFI1-Y394, GAB1-Y689, STAT5A-Y694, DLG3-Y705, and DAPP1-Y139, suggesting these are potential biomarkers of TKI sensitivity. DAPP1, is a novel target of mutant EGFR signaling and Y-139 is the major site of DAPP1 tyrosine phosphorylation. We also uncovered several off-target effects of these TKIs, such as MST1R-Y1238/Y1239 and MET-Y1252/1253. This study provides unique insight into the TKI-mediated modulation of mutant EGFR signaling, which can be applied to the development of biomarkers of EGFR TKI response.

Calcium calmodulin dependent kinase kinase 2 - a novel therapeutic target for gastric adenocarcinoma

Gastric cancer is one of the most common gastrointestinal malignancies and is associated with poor prognosis. Exploring alterations in the proteomic landscape of gastric cancer is likely to provide potential biomarkers for early detection and molecules for targeted therapeutic intervention. Using iTRAQ-based quantitative proteomic analysis, we identified 22 proteins that were overexpressed and 17 proteins that were downregulated in gastric tumor tissues as compared to the adjacent normal tissue. Calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) was found to be 7-fold overexpressed in gastric tumor tissues. Immunohistochemical labeling of tumor tissue microarrays for validation of CAMKK2 overexpression revealed that it was indeed overexpressed in 94% (92 of 98) of gastric cancer cases. Silencing of CAMKK2 using siRNA significantly reduced cell proliferation, colony formation and invasion of gastric cancer cells. Our results demonstrate that CAMKK2 signals in gastric cancer through AMPK activation and suggest that CAMKK2 could be a novel therapeutic target in gastric cancer.

A multi-omic analysis of human naïve CD4+ T cells

Background

Cellular function and diversity are orchestrated by complex interactions of fundamental biomolecules including DNA, RNA and proteins. Technological advances in genomics, epigenomics, transcriptomics and proteomics have enabled massively parallel and unbiased measurements. Such high-throughput technologies have been extensively used to carry out broad, unbiased studies, particularly in the context of human diseases. Nevertheless, a unified analysis of the genome, epigenome, transcriptome and proteome of a single human cell type to obtain a coherent view of the complex interplay between various biomolecules has not yet been undertaken. Here, we report the first multi-omic analysis of human primary naïve CD4+ T cells isolated from a single individual.

Results

Integrating multi-omics datasets allowed us to investigate genome-wide methylation and its effect on mRNA/protein expression patterns, extent of RNA editing under normal physiological conditions and allele specific expression in naïve CD4+ T cells. In addition, we carried out a multi-omic comparative analysis of naïve with primary resting memory CD4+ T cells to identify molecular changes underlying T cell differentiation. This analysis provided mechanistic insights into how several molecules involved in T cell receptor signaling are regulated at the DNA, RNA and protein levels. Phosphoproteomics revealed downstream signaling events that regulate these two cellular states. Availability of multi-omics data from an identical genetic background also allowed us to employ novel proteogenomics approaches to identify individual-specific variants and putative novel protein coding regions in the human genome.

Conclusions

We utilized multiple high-throughput technologies to derive a comprehensive profile of two primary human cell types, naïve CD4+ T cells and memory CD4+ T cells, from a single donor. Through vertical as well as horizontal integration of whole genome sequencing, methylation arrays, RNA-Seq, miRNA-Seq, proteomics, and phosphoproteomics, we derived an integrated and comparative map of these two closely related immune cells and identified potential molecular effectors of immune cell differentiation following antigen encounter.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0225-4) contains supplementary material, which is available to authorized users.

Global phosphotyrosine survey in triple-negative breast cancer reveals activation of multiple tyrosine kinase signaling pathways

Breast cancer is the most prevalent cancer in women worldwide. About 15-20% of all breast cancers are triple negative breast cancer (TNBC) and are often highly aggressive when compared to other subtypes of breast cancers. To better characterize the biology that underlies the TNBC phenotype, we profiled the phosphotyrosine proteome of a panel of twenty-six TNBC cell lines using quantitative high resolution Fourier transform mass spectrometry. A heterogeneous pattern of tyrosine kinase activation was observed based on 1,789 tyrosine-phosphorylated peptides identified from 969 proteins. One of the tyrosine kinases, AXL, was found to be activated in a majority of aggressive TNBC cell lines and was accompanied by a higher level of AXL expression. High levels of AXL expression are correlated with a significant decrease in patient survival. Treatment of cells bearing activated AXL with a humanized AXL antibody inhibited cell proliferation and migration in vitro, and tumor growth in mice. Overall, our global phosphoproteomic analysis provided new insights into the heterogeneity in the activation status of tyrosine kinase pathways in TNBCs. Our approach presents an effective means of identifying important novel biomarkers and targets for therapy such as AXL in TNBC.

Chronic exposure to chewing tobacco selects for overexpression of stearoyl-CoA desaturase in normal oral keratinocytes

Chewing tobacco is a common practice in certain socio-economic sections of southern Asia, particularly in the Indian subcontinent and has been well associated with head and neck squamous cell carcinoma. The molecular mechanisms of chewing tobacco which leads to malignancy remains unclear. In large majority of studies, short-term exposure to tobacco has been evaluated. From a biological perspective, however, long-term (chronic) exposure to tobacco mimics the pathogenesis of oral cancer more closely. We developed a cell line model to investigate the chronic effects of chewing tobacco. Chronic exposure to tobacco resulted in higher cellular proliferation and invasive ability of the normal oral keratinocytes (OKF6/TERT1). We carried out quantitative proteomic analysis of OKF6/TERT1 cells chronically treated with chewing tobacco compared to the untreated cells. We identified a total of 3,636 proteins among which expression of 408 proteins were found to be significantly altered. Among the overexpressed proteins, stearoyl-CoA desaturase (SCD) was found to be 2.6-fold overexpressed in the tobacco treated cells. Silencing/inhibition of SCD using its specific siRNA or inhibitor led to a decrease in cellular proliferation, invasion and colony forming ability of not only the tobacco treated cells but also in a panel of head and neck cancer cell lines. These findings suggest that chronic exposure to chewing tobacco induced carcinogenesis in non-malignant oral epithelial cells and SCD plays an essential role in this process. The current study provides evidence that SCD can act as a potential therapeutic target in head and neck squamous cell carcinoma, especially in patients who are users of tobacco.

Heterogeneity of pancreatic cancer metastases in a single patient revealed by quantitative proteomics

Many patients with pancreatic cancer have metastases to distant organs at the time of initial presentation. Recent studies examining the evolution of pancreatic cancer at the genetic level have shown that clonal complexity of metastatic pancreatic cancer is already initiated within primary tumors, and organ-specific metastases are derived from different subclones. However, we do not yet understand to what extent the evolution of pancreatic cancer contributes to proteomic and signaling alterations. We hypothesized that genetic heterogeneity of metastatic pancreatic cancer results in heterogeneity at the proteome level. To address this, we employed a model system in which cells isolated from three sites of metastasis (liver, lung, and peritoneum) from a single patient were compared. We used a SILAC-based accurate quantitative proteomic strategy combined with high-resolution mass spectrometry to analyze the total proteome and tyrosine phosphoproteome of each of the distal metastases. Our data revealed distinct patterns of both overall proteome expression and tyrosine kinase activities across the three different metastatic lesions. This heterogeneity was significant because it led to differential sensitivity of the neoplastic cells to small molecule inhibitors targeting various kinases and other pathways. For example, R428, a tyrosine kinase inhibitor that targets Axl receptor tyrosine kinase, was able to inhibit cells derived from lung and liver metastases much more effectively than cells from the peritoneal metastasis. Finally, we confirmed that administration of R428 in mice bearing xenografts of cells derived from the three different metastatic sites significantly diminished tumors formed from liver- and lung-metastasis-derived cell lines as compared with tumors derived from the peritoneal metastasis cell line. Overall, our data provide proof-of-principle support that personalized therapy of multiple organ metastases in a single patient should involve the administration of a combination of agents, with each agent targeted to the features of different subclones.

Pancreatic Cancer Database: an integrative resource for pancreatic cancer

Pancreatic cancer is the fourth leading cause of cancer-related death in the world. The etiology of pancreatic cancer is heterogeneous with a wide range of alterations that have already been reported at the level of the genome, transcriptome, and proteome. The past decade has witnessed a large number of experimental studies using high-throughput technology platforms to identify genes whose expression at the transcript or protein levels is altered in pancreatic cancer. Based on expression studies, a number of molecules have also been proposed as potential biomarkers for diagnosis and prognosis of this deadly cancer. Currently, there are no repositories which provide an integrative view of multiple Omics data sets from published research on pancreatic cancer. Here, we describe the development of a web-based resource, Pancreatic Cancer Database (http://www.pancreaticcancerdatabase.org), as a unified platform for pancreatic cancer research. PCD contains manually curated information pertaining to quantitative alterations in miRNA, mRNA, and proteins obtained from small-scale as well as high-throughput studies of pancreatic cancer tissues and cell lines. We believe that PCD will serve as an integrative platform for scientific community involved in pancreatic cancer research.

Differential proteomic analysis of synovial fluid from rheumatoid arthritis and osteoarthritis patients

Background

Rheumatoid arthritis and osteoarthritis are two common musculoskeletal disorders that affect the joints. Despite high prevalence rates, etiological factors involved in these disorders remain largely unknown. Dissecting the molecular aspects of these disorders will significantly contribute to improving their diagnosis and clinical management. In order to identify proteins that are differentially expressed between these two conditions, a quantitative proteomic profiling of synovial fluid obtained from rheumatoid arthritis and osteoarthritis patients was carried out by using iTRAQ labeling followed by high resolution mass spectrometry analysis.

Results

We have identified 575 proteins out of which 135 proteins were found to be differentially expressed by ≥3-fold in the synovial fluid of rheumatoid arthritis and osteoarthritis patients. Proteins not previously reported to be associated with rheumatoid arthritis including, coronin-1A (CORO1A), fibrinogen like-2 (FGL2), and macrophage capping protein (CAPG) were found to be upregulated in rheumatoid arthritis. Proteins such as CD5 molecule-like protein (CD5L), soluble scavenger receptor cysteine-rich domain-containing protein (SSC5D), and TTK protein kinase (TTK) were found to be upregulated in the synovial fluid of osteoarthritis patients. We confirmed the upregulation of CAPG in rheumatoid arthritis synovial fluid by multiple reaction monitoring assay as well as by Western blot. Pathway analysis of differentially expressed proteins revealed a significant enrichment of genes involved in glycolytic pathway in rheumatoid arthritis.

Conclusions

We report here the largest identification of proteins from the synovial fluid of rheumatoid arthritis and osteoarthritis patients using a quantitative proteomics approach. The novel proteins identified from our study needs to be explored further for their role in the disease pathogenesis of rheumatoid arthritis and osteoarthritis.

Sartaj Ahmad and Raja Sekhar Nirujogi contributed equally to this article.

Plasma Proteome Database as a resource for proteomics research: 2014 update

Plasma Proteome Database (PPD; http://www.plasmaproteomedatabase.org/) was initially described in the year 2005 as a part of Human Proteome Organization's (HUPO's) pilot initiative on Human Plasma Proteome Project. Since then, improvements in proteomic technologies and increased throughput have led to identification of a large number of novel plasma proteins. To keep up with this increase in data, we have significantly enriched the proteomic information in PPD. This database currently contains information on 10,546 proteins detected in serum/plasma of which 3784 have been reported in two or more studies. The latest version of the database also incorporates mass spectrometry-derived data including experimentally verified proteotypic peptides used for multiple reaction monitoring assays. Other novel features include published plasma/serum concentrations for 1278 proteins along with a separate category of plasma-derived extracellular vesicle proteins. As plasma proteins have become a major thrust in the field of biomarkers, we have enabled a batch-based query designated Plasma Proteome Explorer, which will permit the users in screening a list of proteins or peptides against known plasma proteins to assess novelty of their data set. We believe that PPD will facilitate both clinical and basic research by serving as a comprehensive reference of plasma proteins in humans and accelerate biomarker discovery and translation efforts.

SILAC-based quantitative proteomic analysis of gastric cancer secretome

PURPOSE

Gastric cancer is a commonly occurring cancer in Asia and one of the leading causes of cancer deaths. However, there is no reliable blood-based screening test for this cancer. Identifying proteins secreted from tumor cells could lead to the discovery of clinically useful biomarkers for early detection of gastric cancer.

EXPERIMENTAL DESIGN

A SILAC-based quantitative proteomic approach was employed to identify secreted proteins that were differentially expressed between neoplastic and non-neoplastic gastric epithelial cells. Proteins from the secretome were subjected to SDS-PAGE and SCX-based fractionation, followed by mass spectrometric analysis on an LTQ-Orbitrap Velos mass spectrometer. Immunohistochemical labeling was employed to validate a subset of candidates using tissue microarrays.

RESULTS

We identified 2205 proteins in the gastric cancer secretome of which 263 proteins were overexpressed greater than fourfold in gastric cancer-derived cell lines as compared to non-neoplastic gastric epithelial cells. Three candidate proteins, proprotein convertase subtilisin/kexin type 9 (PCSK9), lectin mannose binding 2 (LMAN2), and PDGFA-associated protein 1 (PDAP1) were validated by immunohistochemical labeling.

CONCLUSIONS AND CLINICAL RELEVANCE

We report here the largest cancer secretome described to date. The novel biomarkers identified in the current study are excellent candidates for further testing as early detection biomarkers for gastric adenocarcinoma.

Gene Expression Profiling of Tuberculous Meningitis Co-infected with HIV

Tuberculous meningitis (TBM) is a fatal form of Mycobacterium tuberculosis infection of the central nervous system (CNS). The similarities in the clinical and radiological findings in TBM cases with or without HIV make the diagnosis very challenging. Identification of genes, which are differentially expressed in brain tissues of HIV positive and HIV negative TBM patients, would enable better understanding of the molecular aspects of the infection and would also serve as an initial platform to evaluate potential biomarkers. Here, we report the identification of 796 differentially regulated genes in brain tissues of TBM patients co-infected with HIV using oligonucleotide DNA microarrays. We also performed immunohistochemical validation and confirmed the abundance of four gene products-glial fibrillary acidic protein (GFAP), serpin peptidase inhibitor, clade A member 3 (SERPINA3), thymidine phosphorylase (TYMP/ECGF1) and heat shock 70 kDa protein 8 (HSPA8). Our study paves the way for understanding the mechanism of TBM in HIV positive patients and for further validation of potential disease biomarkers.

Overexpression of ribosome binding protein 1 (RRBP1) in breast cancer

The molecular events that lead to malignant transformation and subsequent metastasis of breast carcinoma include alterations in the cells at genome, transcriptome and proteome levels. In this study, we used publicly available gene expression databases to identify those candidate genes which are upregulated at the mRNA level in breast cancers but have not been systematically validated at the protein level. Based on an extensive literature search, we identified ribosome binding protein 1 (RRBP1) as a candidate that is upregulated at the mRNA level in five different studies but its protein expression had not been investigated. Immunohistochemical labeling of breast cancer tissue microarrays was carried out to determine the expression of RRBP1 in a large panel of breast cancers. We found that RRBP1 was overexpressed in 84% (177/219) of breast carcinoma cases tested. The subcellular localization of RRBP1 was mainly observed to be in the cytoplasm with intense staining in the perinuclear region. Our findings suggest that RRBP1 is an interesting molecule that can be further studied for its potential to serve as a breast cancer biomarker. This study also demonstrates how the integration of biological data from available resources in conjunction with systematic evaluation approaches can be successfully applied to clinical proteomics.

Overexpression of Kinesin Associated Protein 3 (KIFAP3) in Breast Cancer

Gene expression profiling studies on breast cancer have generated catalogs of differentially expressed genes. However, many of these genes have not been investigated for their expression at the protein level. It is possible to systematically evaluate such genes in a high-throughput fashion for their overexpression at the protein level using breast cancer tissue microarrays. Our strategy involved integration of information from publicly available repositories of gene expression to prepare a list of genes upregulated at the mRNA level in breast cancer followed by curation of the published literature to identify those genes that were not tested for overexpression at the protein level. We identified Kinesin Associated Protein 3 (KIFAP3) as one such molecule for further validation at the protein level. Immunohistochemical labeling of KIFAP3 using breast cancer tissue microarrays revealed overexpression of KIFAP3 protein in 84% (240/285) of breast cancers indicating the utility of our integrated approach of combining computational analysis with experimental biology.

Abstract 4797: Global phosphotyrosine survey reveals the evidence of activation of multiple tyrosine kinase signaling pathways in basal-like breast cancer cells

About 10-15% of all breast cancers are basal-like breast cancers (BLBCs) that are highly aggressive and lack effective targeted therapeutic options as they lack ER, PR and Her-2. Reversible phosphorylation of proteins is one the most important post-translational modifications and is involved in almost all kinds of cellular processes. Tyrosine phosphorylation (pY) accounts for a minority of total phosphorylation, less than 0.1% by overall abundance. However, pY plays a disproportionately large role in diseases especially in cancer. More than half of the 90 tyrosine kinases identified in the human proteome have been implicated in cancer through gain-of-function mutations, gene amplification, overexpression, or as tumor suppressors and become attractive therapeutic targets. To assess activated tyrosine kinase signaling pathways in TNBCs, we collected a panel of 28 publicly available triple-negative breast cancer cell lines and systematically analyzed the invasiveness and anchorage independent growth of this panel of cell lines. We carried out global quantitative phosphotyrosine profiling using high resolution Fourier transform mass spectrometry. In all, we identified 2,133 tyrosine-phosphorylated peptides from more than 800 proteins. A number of tyrosine kinases were found to be differentially phosphorylated in different sets of breast cancer cells. Pathway analysis revealed that phosphorylation level of certain oncogenic kinases including EGFR, MET, FAK1 and SRC was elevated in a majority of basal-like breast cancer cells. Notably, multiple members of EPH receptor family, EPHA1, EPHA2, EPHA3, EPHB3 and EPHB4, were hyper-phosphorylated in multiple basal-like breast cancer cells. Supervised clustering analysis of phosphorylated proteins identified specific phospho-signatures that correlated with invasiveness and anti-anoikis colony formation ability. One of these hyper-phosphorylated proteins was DYRK2, a dual specificity tyrosine kinase. We observed that DYRK2 had an effect on colony formation and invasive ability on a subset of TNBC cell lines. Overall, our global phosphoproteomic study confirms high heterogeneity in the activation status of tyrosine kinases across triple negative breast cancer cells and suggests that some of them are attractive candidates as therapeutic targets.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4797. doi:1538-7445.AM2012-4797

Abstract 1269: Activation of diverse signaling pathways in pancreatic cancer revealed by phosphoproteomics

Pancreatic cancer is the fourth leading cause of cancer-related deaths in the United States. A number of genes that have somatic mutations including CDKN2A, SMAD4, TP53 and KRAS have been linked to pancreatic cancer. Recently, a comprehensive genomic analysis was performed in a panel of 14 primary pancreatic cancer cell lines and 10 xenografts, in which ∼1,500 somatic alterations were identified and a core set of 12 activated signaling pathways were defined. We decided to use a phosphoproteomics approach to identify tyrosine kinase signaling pathways that could serve as attractive therapeutic targets. Based on immunoblotting with anti-phosphotyrosine antibodies, we observed a heterogeneous pattern of tyrosine phosphorylation in the 14 primary pancreatic cancer cell lines whose genomes have been sequenced previously, which indicated that different tyrosine kinases were likely activated in different cells. To identify these kinases, phosphotyrosine-containing peptides were enriched by immunoprecipitation by anti-phosphotyrosine antibodies and analyzed on a high resolution Fourier transform mass spectrometer coupled to nanoflow reversed phase liquid chromatography. W employed the SILAC methodology to quantitate tyrosine phosphorylated peptides across the panel of cell lines. Our preliminary data indicates that the approach used in this study is suitable for global profiling of tyrosine kinase pathways that are abnormally activated in cancers. The differential alteration of tyrosine kinases was confirmed from the accurate quantitative tyrosine phosphoproteomic analysis. For example, HPNE, a normal pancreatic cell line showed the lowest level of phosphotyrosine levels, while each pancreatic cancer cell line could be differentiated based on its own set of activated tyrosine kinases. Interestingly, some of tyrosine kinases were activated in a small number of cell lines (e.g. AXL, SYK, DDR2, RON and RET), while other tyrosine kinases were found to be activated more broadly (e.g. YES1, TYK2, EPHA2 and EPHB4). Overall, we conclude that these aberrantly activated tyrosine kinases could be targeted to develop newer therapeutic regimens.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 1269. doi:1538-7445.AM2012-1269

Transcriptomic Profiling of Medial Temporal Lobe Epilepsy

Epilepsy is one of the most prevalent neurological disorders affecting ~1% of the population. Medial temporal lobe epilepsy (MTLE) is the most frequent type of epilepsy observed in adults who do not respond to pharmacological treatment. The reason for intractability in these patients has not been systematically studied. Further, no markers are available that can predict the subset of patients who will not respond to pharmacotherapy. To identify potential biomarkers of epileptogenicity, we compared the mRNA profiles of surgically resected tissue from seizure zones with non-seizure zones from cases of intractable MTLE. We identified 413 genes that exhibited ≥2-fold change that were statistically significant across these two groups. Several of these differentially expressed genes have not been previously described in the context of MTLE including claudin 11 (CLDN11) and bone morphogenetic protein receptor, type IB (BMPR1B). In addition, we found significant downregulation of a subset of gamma-aminobutyric acid (GABA) associated genes. We also identified molecules such as BACH2 and ADAMTS15, which are already known to be associated with epilepsy. We validated one upregulated molecule, serine/threonine kinase 31 (STK31) and one downregulated molecule, SMARCA4, by immunohistochemical labeling of tissue sections. These molecules need to be further confirmed in large-scale studies to determine their potential use as diagnostic as well as prognostic markers in intractable MTLE.

Proteogenomic analysis of Mycobacterium tuberculosis by high resolution mass spectrometry

Mass spectrometric sequencing of low abundance, integral membrane proteins, particularly the transmembrane domains, presents challenges that span the multiple phases of sample preparation including solubilization, purification, enzymatic digestion, peptide extraction, and chromatographic separation. We describe a method through which we have obtained high peptide coverage for 12 γ-aminobutyric acid type A receptor (GABA_A receptor) subunits from 2 picomoles of affinity-purified GABA_A receptors from rat brain neocortex. Focusing on the α₁ subunit, we identified peptides covering 96% of the protein sequence from fragmentation spectra (MS2) using a database searching algorithm and deduced 80% of the amino acid residues in the protein from de novo sequencing of Orbitrap spectra. The workflow combined microscale membrane protein solubilization, protein delipidation, in-solution multi-enzyme digestion, multiple stationary phases for peptide extraction, and acquisition of high-resolution full scan and fragmentation spectra. For de novo sequencing of peptides containing the transmembrane domains, timed digestions with chymotrypsin were utilized to generate peptides with overlapping sequences that were then recovered by sequential solid phase extraction using a C4 followed by a porous graphitic carbon stationary phase. The specificity of peptide identifications and amino acid residue sequences was increased by high mass accuracy and charge state assignment to parent and fragment ions. Analysis of three separate brain samples demonstrated that 78% of the sequence of the α₁ subunit was observed in all three replicates with an additional 13% covered in two of the three replicates, indicating a high degree of sequence coverage reproducibility. Label-free quantitative analysis was applied to the three replicates to determine the relative abundances of 11 γ-aminobutyric acid type A receptor subunits. The deep sequence MS data also revealed two N-glycosylation sites on the α₁ subunit, confirmed two splice variants of the γ₂ subunit (γ_2L and γ_2S) and resolved a database discrepancy in the sequence of the α₅ subunit.

Proteogenomic analysis of Mycobacterium tuberculosis by high resolution mass spectrometry

The genome sequencing of H37Rv strain of Mycobacterium tuberculosis was completed in 1998 followed by the whole genome sequencing of a clinical isolate, CDC1551 in 2002. Since then, the genomic sequences of a number of other strains have become available making it one of the better studied pathogenic bacterial species at the genomic level. However, annotation of its genome remains challenging because of high GC content and dissimilarity to other model prokaryotes. To this end, we carried out an in-depth proteogenomic analysis of the M. tuberculosis H37Rv strain using Fourier transform mass spectrometry with high resolution at both MS and tandem MS levels. In all, we identified 3176 proteins from Mycobacterium tuberculosis representing ~80% of its total predicted gene count. In addition to protein database search, we carried out a genome database search, which led to identification of ~250 novel peptides. Based on these novel genome search-specific peptides, we discovered 41 novel protein coding genes in the H37Rv genome. Using peptide evidence and alternative gene prediction tools, we also corrected 79 gene models. Finally, mass spectrometric data from N terminus-derived peptides confirmed 727 existing annotations for translational start sites while correcting those for 33 proteins. We report creation of a high confidence set of protein coding regions in Mycobacterium tuberculosis genome obtained by high resolution tandem mass-spectrometry at both precursor and fragment detection steps for the first time. This proteogenomic approach should be generally applicable to other organisms whose genomes have already been sequenced for obtaining a more accurate catalogue of protein-coding genes.

NetSlim: high-confidence curated signaling maps

We previously developed NetPath as a resource for comprehensive manually curated signal transduction pathways. The pathways in NetPath contain a large number of molecules and reactions which can sometimes be difficult to visualize or interpret given their complexity. To overcome this potential limitation, we have developed a set of more stringent curation and inclusion criteria for pathway reactions to generate high-confidence signaling maps. NetSlim is a new resource that contains this ‘core’ subset of reactions for each pathway for easy visualization and manipulation. The pathways in NetSlim are freely available at http://www.netpath.org/netslim.

Database URL:www.netpath.org/netslim

Quantitative tissue proteomics of esophageal squamous cell carcinoma for novel biomarker discovery

Esophageal squamous cell carcinoma (ESCC) is among the top ten most frequent malignancies worldwide. In this study, our objective was to identify potential biomarkers for ESCC through a quantitative proteomic approach using the isobaric tags for relative and absolute quantitation (iTRAQ) approach. We compared the protein expression profiles of ESCC tumor tissues with the corresponding adjacent normal tissue from ten patients. LC-MS/MS analysis of strong cation exchange chromatography fractions was carried out on an Accurate Mass QTOF mass spectrometer, which led to the identification of 687 proteins. In all, 257 proteins were identified as differentially expressed in ESCC as compared to normal. We found several previously known protein biomarkers to be upregulated in ESCC including thrombospondin 1 (THBS1), periostin 1 (POSTN) and heat shock 70 kDa protein 9 (HSPA9) confirming the validity of our approach. In addition, several novel proteins that had not been reported previously were identified in our screen. These novel biomarker candidates included prosaposin (PSAP), plectin 1 (PLEC1) and protein disulfide isomerase A 4 (PDIA4) that were further validated to be overexpressed by immunohistochemical labeling using tissue microarrays. The success of our study shows that this mass spectrometric strategy can be applied to cancers in general to develop a panel of candidate biomarkers, which can then be validated by other techniques.

A comprehensive map of the human urinary proteome

The study of the human urinary proteome has the potential to offer significant insights into normal physiology as well as disease pathology. The information obtained from such studies could be applied to the diagnosis of various diseases. The high sensitivity, resolution, and mass accuracy of the latest generation of mass spectrometers provides an opportunity to accurately catalog the proteins present in human urine, including those present at low levels. To this end, we carried out a comprehensive analysis of human urinary proteome from healthy individuals using high-resolution Fourier transform mass spectrometry. Importantly, we used the Orbitrap for detecting ions in both MS (resolution 60 000) and MS/MS (resolution 15 000) modes. To increase the depth of our analysis, we characterized both unfractionated as well as lectin-enriched proteins in our experiments. In all, we identified 1,823 proteins with less than 1% false discovery rate, of which 671 proteins have not previously been reported as constituents of human urine. This data set should serve as a comprehensive reference list for future studies aimed at identification and characterization of urinary biomarkers for various diseases.

Gene Expression Profiling of Gastric Cancer

Gastric cancer is the second leading cause of cancer death worldwide, both in men and women. A genomewide gene expression analysis was carried out to identify differentially expressed genes in gastric adenocarcinoma tissues as compared to adjacent normal tissues. We used Agilent's whole human genome oligonucleotide microarray platform representing ~41,000 genes to carry out gene expression analysis. Two-color microarray analysis was employed to directly compare the expression of genes between tumor and normal tissues. Through this approach, we identified several previously known candidate genes along with a number of novel candidate genes in gastric cancer. Testican-1 (SPOCK1) was one of the novel molecules that was 10-fold upregulated in tumors. Using tissue microarrays, we validated the expression of testican-1 by immunohistochemical staining. It was overexpressed in 56% (160/282) of the cases tested. Pathway analysis led to the identification of several networks in which SPOCK1 was among the topmost networks of interacting genes. By gene enrichment analysis, we identified several genes involved in cell adhesion and cell proliferation to be significantly upregulated while those corresponding to metabolic pathways were significantly downregulated. The differentially expressed genes identified in this study are candidate biomarkers for gastric adenoacarcinoma.

The Axl receptor tyrosine kinase is an adverse prognostic factor and a therapeutic target in esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) arises in the backdrop of reflux-induced metaplastic phenomenon known as Barrett esophagus. The prognosis of advanced EAC is dismal, and there is an urgent need for identifying molecular targets for therapy. Serial Analysis of Gene Expression (SAGE) was performed on metachronous mucosal biopsies from a patient who underwent progression to EAC during endoscopic surveillance. SAGE confirmed significant upregulation of Axl "tags" during the multistep progression of Barrett esophagus to EAC. In a cohort of 92 surgically resected EACs, Axl overexpression was associated with shortened median survival on both univariate (p < 0.004) and multivariate (p < 0.036) analysis. Genetic knockdown of Axl receptor tyrosine kinase (RTK) function was enabled in two EAC lines (OE33 and JH-EsoAd1) using lentiviral short hairpin RNA (shRNA). Genetic knockdown of Axl in EAC cell lines inhibited invasion, migration, and in vivo engraftment, which was accompanied by downregulation in the activity of the Ral GTPase proteins (RalA and RalB). Restoration of Ral activation rescued the transformed phenotype of EAC cell lines, suggesting a novel effector mechanism for Axl in cancer cells. Pharmacological inhibition of Axl was enabled using a small molecule antagonist, R428 (Rigel Pharmaceuticals). Pharmacological inhibition of Axl with R428 in EAC cell lines significantly reduced anchorage-independent growth, invasion and migration. Blockade of Axl function abrogated phosphorylation of ERBB2 (Her-2/neu) at the Tyr877 residue, indicative of receptor crosstalk. Axl RTK is an adverse prognostic factor in EAC. The availability of small molecule inhibitors of Axl function provides a tractable strategy for molecular therapy of established EAC.

SILAC-based quantitative proteomic approach to identify potential biomarkers from the esophageal squamous cell carcinoma secretome

The identification of secreted proteins that are differentially expressed between non-neoplastic and esophageal squamous cell carcinoma (ESCC) cells can provide potential biomarkers of ESCC. We used a SILAC-based quantitative proteomic approach to compare the secretome of ESCC cells with that of non-neoplastic esophageal squamous epithelial cells. Proteins were resolved by SDS-PAGE, and tandem mass spectrometry analysis (LC-MS/MS) of in-gel trypsin-digested peptides was carried out on a high-accuracy qTOF mass spectrometer. In total, we identified 441 proteins in the combined secretomes, including 120 proteins with > 2-fold upregulation in the ESCC secretome vs. that of non-neoplastic esophageal squamous epithelial cells. In this study, several potential protein biomarkers previously known to be increased in ESCC including matrix metalloproteinase 1, transferrin receptor, and transforming growth factor beta-induced 68 kDa were identified as overexpressed in the ESCC-derived secretome. In addition, we identified several novel proteins that have not been previously reported to be associated with ESCC. Among the novel candidate proteins identified, protein disulfide isomerase family a member 3 (PDIA3), GDP dissociation inhibitor 2 (GDI2), and lectin galactoside binding soluble 3 binding protein (LGALS3BP) were further validated by immunoblot analysis and immunohistochemical labeling using tissue microarrays. This tissue microarray analysis showed overexpression of protein disulfide isomerase family a member 3, GDP dissociation inhibitor 2, and lectin galactoside binding soluble 3 binding protein in 93%, 93% and 87% of 137 ESCC cases, respectively. Hence, we conclude that these potential biomarkers are excellent candidates for further evaluation to test their role and efficacy in the early detection of ESCC.

Comparative proteomics of human embryonic stem cells and embryonal carcinoma cells

Pluripotent human embryonic stem cells (ESCs) can be differentiated in vitro into a variety of cells which hold promise for transplantation therapy. Human embryonal carcinoma cells (ECCs), stem cells of human teratocarcinomas, are considered a close but malignant counterpart to human ESCs. In this study, a comprehensive quantitative proteomic analysis of ESCs and ECCs was carried out using the iTRAQ method. Using two-dimensional LC and MS/MS analyses, we identified and quantitated approximately 1800 proteins. Among these are proteins associated with pluripotency and development as well as tight junction signaling and TGFbeta receptor pathway. Nearly approximately 200 proteins exhibit more than twofold difference in abundance between ESCs and ECCs. Examples of early developmental markers high in ESCs include beta-galactoside-binding lectin, undifferentiated embryonic cell transcription factor-1, DNA cytosine methyltransferase 3beta isoform-B, melanoma antigen family-A4, and interferon-induced transmembrane protein-1. In contrast, CD99-antigen (CD99), growth differentiation factor-3, cellular retinoic acid binding protein-2, and developmental pluripotency associated-4 were among the highly expressed proteins in ECCs. Several proteins that were highly expressed in ECCs such as heat shock 27 kDa protein-1, mitogen-activated protein kinase kinase-1, nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor like-2, and S100 calcium-binding protein-A4 have also been attributed to malignancy in other systems. Importantly, immunocytochemistry was used to validate the proteomic analyses for a subset of the proteins. In summary, this is the first large-scale quantitative proteomic study of human ESCs and ECCs, which provides critical information about the regulators of these two closely related, but developmentally distinct, stem cells.

Overexpression of periostin and lumican in esophageal squamous cell carcinoma

To identify biomarkers for early detection for esophageal squamous cell carcinoma (ESCC), we previously carried out a genome-wide gene expression profiling study using an oligonucleotide microarray platform. This analysis led to identification of several transcripts that were significantly upregulated in ESCC compared to the adjacent normal epithelium. In the current study, we performed immunohistochemical analyses of protein products for two candidates genes identified from the DNA microarray analysis, periostin (POSTN) and lumican (LUM), using tissue microarrays. Increased expression of both periostin and lumican was observed in 100% of 137 different ESCC samples arrayed on tissue microarrays. Increased expression of periostin and lumican was observed in carcinoma as well as in stromal cell in the large majority of cases. These findings suggest that these candidates can be investigated in the sera of ESCC patients using ELISA or multiple reaction monitoring (MRM) type assays to further explore their utility as biomarkers.

A compendium of potential biomarkers of pancreatic cancer

Akhilesh Pandey and colleagues describe a compendium of potential biomarkers that can be systematically validated by the pancreatic cancer community.

Temporal analysis of neural differentiation using quantitative proteomics

The ability to derive neural progenitors, differentiated neurons and glial cells from human embryonic stem cells (hESCs) with high efficiency holds promise for a number of clinical applications. However, investigating the temporal events is crucial for defining the underlying mechanisms that drive this process of differentiation along different lineages. We carried out quantitative proteomic profiling using a multiplexed approach capable of analyzing eight different samples simultaneously to monitor the temporal dynamics of protein abundance as human embryonic stem cells differentiate into motor neurons or astrocytes. With this approach, a catalog of approximately 1200 proteins along with their relative quantitative expression patterns was generated. The differential expression of the large majority of these proteins has not previously been reported or studied in the context of neural differentiation. As expected, two of the widely used markers of pluripotency, alkaline phosphatase (ALPL) and LIN28, were found to be downregulated during differentiation, while S-100 and tenascin C were upregulated in astrocytes. Neurofilament 3 protein, doublecortin and CAM kinase-like 1 and nestin proteins were upregulated during motor neuron differentiation. We identified a number of proteins whose expression was largely confined to specific cell types, embryonic stem cells, embryoid bodies and differentiating motor neurons. For example, glycogen phosphorylase (PYGL) and fatty acid binding protein 5 (FABP5) were enriched in ESCs, while beta spectrin (SPTBN5) was highly expressed in embryoid bodies. Karyopherin, heat shock 27 kDa protein 1 and cellular retinoic acid binding protein 2 (CRABP2) were upregulated in differentiating motor neurons but were downregulated in mature motor neurons. We validated some of the novel markers of the differentiation process using immunoblotting and immunocytochemical labeling. To our knowledge, this is the first large-scale temporal proteomic profiling of human stem cell differentiation into neural cell types highlighting proteins with limited or undefined roles in neural fate.

Genomewide mRNA profiling of esophageal squamous cell carcinoma for identification of cancer biomarkers

Cancer of the esophagus is of two main types, each with distinct etiological and pathological characteristics. Esophageal squamous cell carcinoma (ESCC) is predominant type of esophageal cancers worldwide comprising almost 95% of cases. While ESCC is prevalent in the developing world, esophageal adenocarcinoma is commonly seen in the developed country, usually in association with Barrett's esophagus. In spite of its higher prevalence, ESCC has not been studied as intensively as esophageal adenocarcinoma. ESCC and esophageal adenocarcinoma are common cancers worldwide with poor survival rate among patients mainly because both of these cancers lack early biomarkers of identification. Molecular mechanisms contributing to initiation and progression of esophageal squamous cell carcinoma are still poorly understood. Development of DNA microarray technology allows high-throughput identification of gene expression profiles in cancers. In order to identify molecules as candidates for early diagnosis and/or as therapeutic targets, we analyzed the mRNA expression profiles of 20 cases of ESCC using whole genome DNA microarrays. A total of 2,235 genes were differentially regulated in the tumors as compared to the corresponding adjacent normal epithelium of which 881 were significantly upregulated. We validated two molecules that were not previously reported to be overexpressed in ESCC, oral cancer overexpressed 2 (ORAOV2) and fibroblast activation protein (FAP), by immunohistochemical labeling of tissue microarrays and archival tissue sections and found that they were overexpressed in 98% (116/118) and 68% (79/116) of cases, respectively. By gene enrichment analysis, we identified significant downregulation of several genes in the arachidonic acid metabolic pathway. Overall, using this approach we have identified a number of promising novel candidates that can be validated further for their potential to serve as biomarkers for ESCC.

Human Protein Reference Database--2009 update

Human Protein Reference Database (HPRD--http://www.hprd.org/), initially described in 2003, is a database of curated proteomic information pertaining to human proteins. We have recently added a number of new features in HPRD. These include PhosphoMotif Finder, which allows users to find the presence of over 320 experimentally verified phosphorylation motifs in proteins of interest. Another new feature is a protein distributed annotation system--Human Proteinpedia (http://www.humanproteinpedia.org/)--through which laboratories can submit their data, which is mapped onto protein entries in HPRD. Over 75 laboratories involved in proteomics research have already participated in this effort by submitting data for over 15,000 human proteins. The submitted data includes mass spectrometry and protein microarray-derived data, among other data types. Finally, HPRD is also linked to a compendium of human signaling pathways developed by our group, NetPath (http://www.netpath.org/), which currently contains annotations for several cancer and immune signaling pathways. Since the last update, more than 5500 new protein sequences have been added, making HPRD a comprehensive resource for studying the human proteome.

Comparisons of tyrosine phosphorylated proteins in cells expressing lung cancer-specific alleles of EGFR and KRAS

We have used unbiased phosphoproteomic approaches, based on quantitative mass spectrometry using stable isotope labeling with amino acids in cell culture (SILAC), to identify tyrosine phosphorylated proteins in isogenic human bronchial epithelial cells (HBECs) and human lung adenocarcinoma cell lines, expressing either of the two mutant alleles of EGFR (L858R and Del E746-A750), or a mutant KRAS allele, which are common in human lung adenocarcinomas. Tyrosine phosphorylation of signaling molecules was greater in HBECs expressing the mutant EGFRs than in cells expressing WT EGFR or mutant KRAS. Receptor tyrosine kinases (such as EGFR, ERBB2, MET, and IGF1R), and Mig-6, an inhibitor of EGFR signaling, were more phosphorylated in HBECs expressing mutant EGFR than in cells expressing WT EGFR or mutant RAS. Phosphorylation of some proteins differed in the two EGFR mutant-expressing cells; for example, some cell junction proteins (beta-catenin, plakoglobin, and E-cadherin) were more phosphorylated in HBECs expressing L858R EGFR than in cells expressing Del EGFR. There were also differences in degree of phosphorylation at individual tyrosine sites within a protein; for example, a previously uncharacterized phosphorylation site in the nucleotide-binding loop of the kinase domains of EGFR (Y727), ERBB2 (Y735), or ERBB4 (Y733), is phosphorylated significantly more in HBECs expressing the deletion mutant than in cells expressing the wild type or L858R EGFR. Signaling molecules not previously implicated in ERBB signaling, such as polymerase transcript release factor (PTRF), were also phosphorylated in cells expressing mutant EGFR. Bayesian network analysis of these and other datasets revealed that PTRF might be a potentially important component of the ERBB signaling network.

The nuclear receptor corepressor (N-CoR) contains three isoleucine motifs (I/LXXII) that serve as receptor interaction domains (IDs)

Unliganded thyroid hormone receptors (TRs) repress transcription through recruitment of corepressors, including nuclear receptor corepressor (N-CoR). We find that N-CoR contains three interaction domains (IDs) that bind to TR, rather than the previously reported two. The hitherto unrecognized ID (ID3) serves as a fully functional TR binding site, both in vivo and in vitro, and may be the most important for TR binding. Each ID motif contains a conserved hydrophobic core (I/LXXII) that resembles the hydrophobic core of nuclear receptor boxes (LXXLL), which mediates p160 coactivator binding to liganded nuclear receptors. Although the integrity of the I/LXXII motif is required for ID function, substitution of ID isoleucines with leucines did not allow ID peptides to bind to liganded TR, and substitution of NR box leucines with isoleucines did not allow NR box peptides to bind unliganded TR. This indicates that the binding preferences of N-CoR for unliganded TR and p160s for liganded TR are not dictated solely by the identity of conserved hydrophobic residues within their TR binding motifs. Examination of sequence conservation between IDs, and mutational analysis of individual IDs, suggests that they are comprised of the central hydrophobic core and distinct adjacent sequences that may make unique contacts with the TR surface. Accordingly, a hybrid peptide that contains distinct adjacent sequences from ID3 and ID1 shows enhanced binding to TR.