Combined thioredoxin reductase and glutaminase inhibition exerts synergistic anti-tumor activity in MYC-high high-grade serous ovarian carcinoma.

Approximately 50%-55% of high-grade serous ovarian carcinoma (HGSOC) patients have MYC oncogenic pathway activation. Because MYC is not directly targetable, we have analyzed molecular pathways enriched in MYC-high HGSOC tumors to identify potential therapeutic targets. Here, we report that MYC-high HGSOC tumors show enrichment in genes controlled by NRF2, an antioxidant signaling pathway, along with increased thioredoxin redox activity. Treatment of MYC-high HGSOC tumors cells with US Food and Drug Administration (FDA)-approved thioredoxin reductase 1 (TrxR1) inhibitor auranofin resulted in significant growth suppression and apoptosis in MYC-high HGSOC cells in vitro and also significantly reduced tumor growth in an MYC-high HGSOC patient-derived tumor xenograft. We found that auranofin treatment inhibited glycolysis in MYC-high cells via oxidation-induced GAPDH inhibition. Interestingly, in response to auranofin-induced glycolysis inhibition, MYC-high HGSOC cells switched to glutamine metabolism for survival. Depletion of glutamine with either glutamine starvation or glutaminase (GLS1) inhibitor CB-839 exerted synergistic anti-tumor activity with auranofin in HGSOC cells and OVCAR-8 cell line xenograft. These findings suggest that applying a combined therapy of GLS1 inhibitor and TrxR1 inhibitor could effectively treat MYC-high HGSOC patients.

Protein-coding potential of non-canonical open reading frames in human transcriptome.

In recent years, proteogenomics and ribosome profiling studies have identified a large number of proteins encoded by noncoding regions in the human genome. They are encoded by small open reading frames (sORFs) in the untranslated regions (UTRs) of mRNAs and long non-coding RNAs (lncRNAs). These sORF encoded proteins (SEPs) are often <150AA and show poor evolutionary conservation. A subset of them have been functionally characterized and shown to play an important role in fundamental biological processes including cardiac and muscle function, DNA repair, embryonic development and various human diseases. How many novel protein-coding regions exist in the human genome and what fraction of them are functionally important remains a mystery. In this review, we discuss current progress in unraveling SEPs, approaches used for their identification, their limitations and reliability of these identifications. We also discuss functionally characterized SEPs and their involvement in various biological processes and diseases. Lastly, we provide insights into their distinctive features compared to canonical proteins and challenges associated with annotating these in protein reference databases.

Protein kinase TgCDPK7 regulates vesicular trafficking and phospholipid synthesis in Toxoplasma gondii.

Apicomplexan parasites are causative agents of major human diseases. Calcium Dependent Protein Kinases (CDPKs) are crucial components for the intracellular development of apicomplexan parasites and are thus considered attractive drug targets. CDPK7 is an atypical member of this family, which initial characterization suggested to be critical for intracellular development of both Apicomplexa Plasmodium falciparum and Toxoplasma gondii. However, the mechanisms via which it regulates parasite replication have remained unknown. We performed quantitative phosphoproteomics of T. gondii lacking TgCDPK7 to identify its parasitic targets. Our analysis lead to the identification of several putative TgCDPK7 substrates implicated in critical processes like phospholipid (PL) synthesis and vesicular trafficking. Strikingly, phosphorylation of TgRab11a via TgCDPK7 was critical for parasite intracellular development and protein trafficking. Lipidomic analysis combined with biochemical and cellular studies confirmed that TgCDPK7 regulates phosphatidylethanolamine (PE) levels in T. gondii. These studies provide novel insights into the regulation of these processes that are critical for parasite development by TgCDPK7.

Signaling alterations in oral keratinocytes in response to shisha and crude tobacco extract.

BACKGROUND: Tobacco consumption in smoking and non-smoking forms has been consequential in the rise of oral cancer cases. Among different forms, epidemiological studies from Middle Eastern countries and rural parts of northern India have reported increasing association of oral cancer with waterpipe (hookah) smoking. However, molecular mechanisms and role played by waterpipe smoking in the onset of oral carcinogenesis remains unexplored. METHODS: In this study, immortalized normal human oral keratinocytes were chronically treated with extracts of two varieties of waterpipe tobacco-crude tobacco and processed shisha. Phenotypic changes and molecular aberrations were examined using cell culture-based assays and mass spectrometry-based quantitative proteomic analysis, respectively. Bioinformatics analysis was utilized to analyze proteomics data and identify dysregulated pathways. RESULTS: Our data indicate that chronic treatment with waterpipe tobacco extracts increased proliferation, invasion, migration, and significant dysregulation of protein expression in oral keratinocytes. Altered expression of proteins involved in interferon signaling pathway were observed with both varieties of tobacco. Overexpression of cholesterol metabolism and vesicle-mediated transport proteins were identified exclusively in cells treated with crude tobacco extract. Bioinformatics analyses revealed different oncogenic response in oral cells based on the type of waterpipe tobacco used. CONCLUSIONS: This study may serve as a useful resource in understanding the early onset of oral cancer attributed to waterpipe smoking.

Integrating transcriptomic and proteomic data for accurate assembly and annotation of genomes.

Complementing genome sequence with deep transcriptome and proteome data could enable more accurate assembly and annotation of newly sequenced genomes. Here, we provide a proof-of-concept of an integrated approach for analysis of the genome and proteome of Anopheles stephensi, which is one of the most important vectors of the malaria parasite. To achieve broad coverage of genes, we carried out transcriptome sequencing and deep proteome profiling of multiple anatomically distinct sites. Based on transcriptomic data alone, we identified and corrected 535 events of incomplete genome assembly involving 1196 scaffolds and 868 protein-coding gene models. This proteogenomic approach enabled us to add 365 genes that were missed during genome annotation and identify 917 gene correction events through discovery of 151 novel exons, 297 protein extensions, 231 exon extensions, 192 novel protein start sites, 19 novel translational frames, 28 events of joining of exons, and 76 events of joining of adjacent genes as a single gene. Incorporation of proteomic evidence allowed us to change the designation of more than 87 predicted "noncoding RNAs" to conventional mRNAs coded by protein-coding genes. Importantly, extension of the newly corrected genome assemblies and gene models to 15 other newly assembled Anopheline genomes led to the discovery of a large number of apparent discrepancies in assembly and annotation of these genomes. Our data provide a framework for how future genome sequencing efforts should incorporate transcriptomic and proteomic analysis in combination with simultaneous manual curation to achieve near complete assembly and accurate annotation of genomes.

A draft map of the human proteome.

The availability of human genome sequence has transformed biomedical research over the past decade. However, an equivalent map for the human proteome with direct measurements of proteins and peptides does not exist yet. Here we present a draft map of the human proteome using high-resolution Fourier-transform mass spectrometry. In-depth proteomic profiling of 30 histologically normal human samples, including 17 adult tissues, 7 fetal tissues and 6 purified primary haematopoietic cells, resulted in identification of proteins encoded by 17,294 genes accounting for approximately 84% of the total annotated protein-coding genes in humans. A unique and comprehensive strategy for proteogenomic analysis enabled us to discover a number of novel protein-coding regions, which includes translated pseudogenes, non-coding RNAs and upstream open reading frames. This large human proteome catalogue (available as an interactive web-based resource at http://www.humanproteomemap.org) will complement available human genome and transcriptome data to accelerate biomedical research in health and disease.

Quantitative phosphoproteomic analysis reveals system-wide signaling pathways regulated by site-specific phosphorylation of Keratin-8 in skin squamous cell carcinoma derived cell line.

Keratin 8/18, a simple epithelia specific keratin pair, is often aberrantly expressed in squamous cell carcinomas (SCC) where its expression is correlated with increased invasion and poor prognosis. Majority of Keratin 8 (K8) functions are governed by its phosphorylation at Serine73 (head-domain) and Serine431 (tail-domain) residues. Although, deregulation of K8 phosphorylation is associated with progression of different carcinomas, its role in skin-SCC and the underlying mechanism is obscure. In this direction, we performed tandem mass tag-based quantitative phosphoproteomics by expressing K8 wild type, phosphodead, and phosphomimetic mutants in K8-deficient A431 cells. Further analysis of our phosphoproteomics data showed a significant proportion of total phosphoproteome associated with migratory, proliferative, and invasive potential of these cells to be differentially phosphorylated. Differential phosphorylation of CDK1T14,Y15 , EIF4EBP1T46,T50 , EIF4BS422 , AKT1S1T246,S247 , CTTN1T401,S405,Y421 , and CAP1S307/309 in K8-S73A/D mutant and CTTN1T401,S405,Y421 , BUB1BS1043 , and CARHSP1S30,S32 in K8-S431A/D mutants as well as some anonymous phosphosites including MYCS176 , ZYXS344 , and PNNS692 could be potential candidates associated with K8 phosphorylation mediated tumorigenicity. Biochemical validation followed by phenotypic analysis further confirmed our quantitative phosphoproteomics data. In conclusion, our study provides the first global picture of K8 site-specific phosphorylation function in neoplastic progression of A431 cells and suggests various potential starting points for further mechanistic studies.

Cerebrospinal Fluid from Sporadic Amyotrophic Lateral Sclerosis Patients Induces Mitochondrial and Lysosomal Dysfunction.

In our laboratory, we have developed (1) an in vitro model of sporadic Amyotrophic Lateral Sclerosis (sALS) involving exposure of motor neurons to cerebrospinal fluid (CSF) from sALS patients and (2) an in vivo model involving intrathecal injection of sALS-CSF into rat pups. In the current study, we observed that spinal cord extract from the in vivo sALS model displayed elevated reactive oxygen species (ROS) and mitochondrial dysfunction. Quantitative proteomic analysis of sub-cellular fractions from spinal cord of the in vivo sALS model revealed down-regulation of 35 mitochondrial proteins and 4 lysosomal proteins. Many of the down-regulated mitochondrial proteins contribute to alterations in respiratory chain complexes and organellar morphology. Down-regulated lysosomal proteins Hexosaminidase, Sialidase and Aryl sulfatase also displayed lowered enzyme activity, thus validating the mass spectrometry data. Proteomic analysis and validation by western blot indicated that sALS-CSF induced the over-expression of the pro-apoptotic mitochondrial protein BNIP3L. In the in vitro model, sALS-CSF induced neurotoxicity and elevated ROS, while it lowered the mitochondrial membrane potential in rat spinal cord mitochondria in the in vivo model. Ultra structural alterations were evident in mitochondria of cultured motor neurons exposed to ALS-CSF. These observations indicate the first line evidence that sALS-CSF mediated mitochondrial and lysosomal defects collectively contribute to the pathogenesis underlying sALS.

Annotation of the zebrafish genome through an integrated transcriptomic and proteomic analysis.

Accurate annotation of protein-coding genes is one of the primary tasks upon the completion of whole genome sequencing of any organism. In this study, we used an integrated transcriptomic and proteomic strategy to validate and improve the existing zebrafish genome annotation. We undertook high-resolution mass-spectrometry-based proteomic profiling of 10 adult organs, whole adult fish body, and two developmental stages of zebrafish (SAT line), in addition to transcriptomic profiling of six organs. More than 7,000 proteins were identified from proteomic analyses, and ∼ 69,000 high-confidence transcripts were assembled from the RNA sequencing data. Approximately 15% of the transcripts mapped to intergenic regions, the majority of which are likely long non-coding RNAs. These high-quality transcriptomic and proteomic data were used to manually reannotate the zebrafish genome. We report the identification of 157 novel protein-coding genes. In addition, our data led to modification of existing gene structures including novel exons, changes in exon coordinates, changes in frame of translation, translation in annotated UTRs, and joining of genes. Finally, we discovered four instances of genome assembly errors that were supported by both proteomic and transcriptomic data. Our study shows how an integrative analysis of the transcriptome and the proteome can extend our understanding of even well-annotated genomes.

Phosphotyrosine profiling identifies ephrin receptor A2 as a potential therapeutic target in esophageal squamous-cell carcinoma.

Esophageal squamous-cell carcinoma (ESCC) is one of the most common malignancies in Asia. Currently, surgical resection of early-stage tumor is the best available treatment. However, most patients present late when surgery is not an option. Data suggest that chemotherapy regimens are inadequate for clinical management of advanced cancer. Targeted therapy has emerged as one of the most promising approaches to treat several malignancies. A prerequisite for developing targeted therapy is prior knowledge of proteins and pathways that drive proliferation in malignancies. We carried out phosphotyrosine profiling across four different ESCC cell lines and compared it to non-neoplastic Het-1A cell line to identify activated tyrosine kinase signaling pathways in ESCC. A total of 278 unique phosphopeptides were identified across these cell lines. This included several tyrosine kinases and their substrates that were hyperphosphorylated in ESCC. Ephrin receptor A2 (EPHA2), a receptor tyrosine kinase, was hyperphosphorylated in all the ESCC cell lines used in the study. EPHA2 is reported to be oncogenic in several cancers and is also known to promote metastasis. Immunohistochemistry-based studies have revealed EPHA2 is overexpressed in nearly 50% of ESCC. We demonstrated EPHA2 as a potential therapeutic target in ESCC by carrying out siRNA-based knockdown studies. Knockdown of EPHA2 in ESCC cell line TE8 resulted in significant decrease in cell proliferation and invasion, suggesting it is a promising therapeutic target in ESCC that warrants further evaluation.

Identification and characterization of proteins encoded by chromosome 12 as part of chromosome-centric human proteome project.

Chromosome-centric human proteome project (C-HPP) is a global initiative to comprehensively characterize proteins encoded by genes across all human chromosomes by teams focusing on individual chromosomes. Here, we report mass spectrometry-based identification and characterization of proteins encoded by genes on chromosome 12. Our study is based on proteomic profiling of 30 different histologically normal human tissues and cell types using high-resolution mass spectrometry. In our analysis, we identified 1,535 proteins encoded by 836 genes on human chromosome 12. This includes 89 genes that are designated as "missing proteins" by "neXtProt" as they did not have any prior evidence either by mass spectrometry or by antibody-based detection methods. We identified several variant peptides that reflected coding SNPs annotated in dbSNP database. We also confirmed the start sites of ∼200 proteins by identifying protein N-terminal acetylated peptides. We also identified alternative start sites for 11 proteins that were not annotated in public databases until now. Most importantly, we identified 12 novel protein coding regions on chromosome 12 using our proteogenomics strategy. All of the 12 regions have been annotated as pseudogenes in public databases. This study demonstrates that there is scope for significantly improving annotation of protein coding genes in the human genome using mass-spectrometry-derived data. Individual efforts as part of C-HPP initiative should significantly contribute toward enriching human protein annotation. The data have been deposited to ProteomeXchange with identifier PXD000561.

Acute neuroinflammation promotes a metabolic shift that alters extracellular vesicle cargo in the mouse brain cortex.

Neuroinflammation is initiated through microglial activation and cytokine release which can be induced through lipopolysaccharide treatment (LPS) leading to a transcriptional cascade culminating in the differential expression of target proteins. These differentially expressed proteins can then be packaged into extracellular vesicles (EVs), a form of cellular communication, further propagating the neuroinflammatory response over long distances. Despite this, the EV proteome in the brain, following LPS treatment, has not been investigated. Brain tissue and brain derived EVs (BDEVs) isolated from the cortex of LPS-treated mice underwent thorough characterisation to meet the minimal information for studies of extracellular vesicles guidelines before undergoing mass spectrometry analysis to identify the differentially expressed proteins. Fourteen differentially expressed proteins were identified in the LPS brain tissue samples compared to the controls and 57 were identified in the BDEVs isolated from the LPS treated mice compared to the controls. This included proteins associated with the initiation of the inflammatory response, epigenetic regulation, and metabolism. These results allude to a potential link between small EV cargo and early inflammatory signalling.

Microglial activation induces nitric oxide signalling and alters protein S-nitrosylation patterns in extracellular vesicles.

Neuroinflammation is an underlying feature of neurodegenerative conditions, often appearing early in the aetiology of a disease. Microglial activation, a prominent initiator of neuroinflammation, can be induced through lipopolysaccharide (LPS) treatment resulting in expression of the inducible form of nitric oxide synthase (iNOS), which produces nitric oxide (NO). NO post-translationally modifies cysteine thiols through S-nitrosylation, which can alter function of the target protein. Furthermore, packaging of these NO-modified proteins into extracellular vesicles (EVs) allows for the exertion of NO signalling in distant locations, resulting in further propagation of the neuroinflammatory phenotype. Despite this, the NO-modified proteome of activated microglial EVs has not been investigated. This study aimed to identify the protein post-translational modifications NO signalling induces in neuroinflammation. EVs isolated from LPS-treated microglia underwent mass spectral surface imaging using time of flight-secondary ion mass spectrometry (ToF-SIMS), in addition to iodolabelling and comparative proteomic analysis to identify post-translation S-nitrosylation modifications. ToF-SIMS imaging successfully identified cysteine thiol side chains modified through NO signalling in the LPS treated microglial-derived EV proteins. In addition, the iodolabelling proteomic analysis revealed that the EVs from LPS-treated microglia carried S-nitrosylated proteins indicative of neuroinflammation. These included known NO-modified proteins and those associated with LPS-induced microglial activation that may play an essential role in neuroinflammatory communication. Together, these results show activated microglia can exert broad NO signalling changes through the selective packaging of EVs during neuroinflammation.

Phosphotyrosine Profiling Using SILAC.

Tyrosine phosphorylation on proteins is an important posttranslational modification that regulates various processes in cells. Mass spectrometry-based phosphotyrosine profiling can reveal tyrosine kinase signaling activity in cells. Using quantitative proteomics strategies such as stable isotope labeling with amino acids in cell culture (SILAC) allows comparison of tyrosine kinase signaling activity across two to -three different conditions. In this book chapter, we discuss the reagents required and a step-by-step protocol to carry out phosphotyrosine profiling using SILAC.

Comparative analysis of tangential flow filtration and ultracentrifugation, both combined with subsequent size exclusion chromatography, for the isolation of small extracellular vesicles.

Small extracellular vesicles (sEVs) provide major promise for advances in cancer diagnostics, prognostics, and therapeutics, ascribed to their distinctive cargo reflective of pathophysiological status, active involvement in intercellular communication, as well as their ubiquity and stability in bodily fluids. As a result, the field of sEV research has expanded exponentially. Nevertheless, there is a lack of standardisation in methods for sEV isolation from cells grown in serum-containing media. The majority of researchers use serum-containing media for sEV harvest and employ ultracentrifugation as the primary isolation method. Ultracentrifugation is inefficient as it is devoid of the capacity to isolate high sEV yields without contamination of non-sEV materials or disruption of sEV integrity. We comprehensively evaluated a protocol using tangential flow filtration and size exclusion chromatography to isolate sEVs from a variety of human and murine cancer cell lines, including HeLa, MDA-MB-231, EO771 and B16F10. We directly compared the performance of traditional ultracentrifugation and tangential flow filtration methods, that had undergone further purification by size exclusion chromatography, in their capacity to separate sEVs, and rigorously characterised sEV properties using multiple quantification devices, protein analyses and both image and nano-flow cytometry. Ultracentrifugation and tangential flow filtration both enrich consistent sEV populations, with similar size distributions of particles ranging up to 200 nm. However, tangential flow filtration exceeds ultracentrifugation in isolating significantly higher yields of sEVs, making it more suitable for large-scale research applications. Our results demonstrate that tangential flow filtration is a reliable and robust sEV isolation approach that surpasses ultracentrifugation in yield, reproducibility, time, costs and scalability. These advantages allow for implementation in comprehensive research applications and downstream investigations.

Protocol for purification and identification of MHC class I immunopeptidome from cancer cell lines.

Major histocompatibility complexes (MHC) play a critical role in immunity by presenting peptides on the cell surface for T cell recognition. Identification of these peptides can be valuable to develop vaccines or immunotherapeutic strategies for infectious diseases and cancers. Mass spectrometry is the only tool available for unbiased identification of the immunopeptidome. Here, we describe a protocol for purification and identification of MHC class I peptides, including in-house purification of anti-MHC-antibody from hybridoma cells and the LC-MS/MS analysis of MHC-I bound peptides.

Temporal Quantitative Proteomics Reveals Proteomic and Phosphoproteomic Alterations Associated with Adaptive Response to Hypoxia in Melanoma Cells.

Hypoxia is a common feature in various solid tumours, including melanoma. Cancer cells in hypoxic environments are resistant to both chemotherapy and radiation. Hypoxia is also associated with immune suppression. Identification of proteins and pathways that regulate cancer cell survival in hypoxic environments can reveal potential vulnerabilities that can be exploited to improve the efficacy of anticancer therapies. We carried out temporal proteomic and phosphoproteomic profiling in melanoma cell lines to identify hypoxia-induced protein expression and phosphorylation changes. By employing a TMT-based quantitative proteomics strategy, we report the identification and quantitation of >7000 proteins and >10,000 phosphosites in melanoma cell lines grown in hypoxia. Proteomics data show metabolic reprogramming as one of the prominent adaptive responses in hypoxia. We identify several novel hypoxia-mediated phosphorylation changes that have not been reported before. They reveal kinase signalling pathways that are potentially involved in modulating cellular response to hypoxia. In addition to known protein expression changes, we identify several novel proteomic alterations associated with adaptive response to hypoxia. We show that cancer cells require the ubiquitin-proteasome system to survive in both normoxia and hypoxia. Inhibition of proteasome activity affects cell survival and may provide a novel therapeutic avenue to target cancer cells in hypoxia. Our study can serve as a valuable resource to pursue novel candidates to target hypoxia in cancers and improve the efficacy of anticancer therapies.

An integrated approach for identification of a panel of candidate genes arbitrated for invasion and metastasis in oral squamous cell carcinoma.

Oral squamous cell carcinoma (OSCC) is known for its aggressiveness associated with poor prognosis. The molecular mechanisms underlying the invasion and metastasis are still poorly understood. An improved understanding of these mechanisms shall precede the development of new diagnostic tools and targeted therapies. We report an integrated approach using bioinformatics to predict candidate genes, coupled with proteomics and immunohistochemistry for validating their presence and involvement in OSCC pathways heralding invasion and metastasis. Four genes POSTN, TNC, CAV1 and FSCN1 were identified. A protein-protein interaction network analysis teamed with pathway analysis led us to propose the role of the identified genes in invasion and metastasis in OSCC. Further analyses of archived FFPE blocks of various grades of oral cancer was carried out using TMT-based mass spectrometry and immunohistochemistry. Results of this study expressed a strong communiqué and interrelationship between these candidate genes. This study emphasizes the significance of a molecular biomarker panel as a diagnostic tool and its correlation with the invasion and metastatic pathway of OSCC. An insight into the probable association of CAF's and these biomarkers in the evolution and malignant transformation of OSCC further magnifies the molecular-biological spectrum of OSCC tumour microenvironment.

Proteomic and phosphoproteomic profiling of shammah induced signaling in oral keratinocytes.

Shammah is a smokeless tobacco product often mixed with lime, ash, black pepper and flavorings. Exposure to shammah has been linked with dental diseases and oral squamous cell carcinoma. There is limited literature on the prevalence of shammah and its role in pathobiology of oral cancer. In this study, we developed a cellular model to understand the effect of chronic shammah exposure on oral keratinocytes. Chronic exposure to shammah resulted in increased proliferation and invasiveness of non-transformed oral keratinocytes. Quantitative proteomics of shammah treated cells compared to untreated cells led to quantification of 4712 proteins of which 402 were found to be significantly altered. In addition, phosphoproteomics analysis of shammah treated cells compared to untreated revealed hyperphosphorylation of 36 proteins and hypophosphorylation of 83 proteins (twofold, p-value ≤ 0.05). Bioinformatics analysis of significantly altered proteins showed enrichment of proteins involved in extracellular matrix interactions, necroptosis and peroxisome mediated fatty acid oxidation. Kinase-Substrate Enrichment Analysis showed significant increase in activity of kinases such as ROCK1, RAF1, PRKCE and HIPK2 in shammah treated cells. These results provide better understanding of how shammah transforms non-neoplastic cells and warrants additional studies that may assist in improved early diagnosis and treatment of shammah induced oral cancer.

Proteomic Alterations Associated with Oral Cancer Patients with Tobacco Using Habits.

Tobacco abuse is a major risk factor associated with the development of oral squamous cell carcinoma. Differences in molecular aberrations induced by tobacco exposure by chewing or smoking form are not well studied in case of oral cancer. We used tandem mass tag-based quantitative proteomic approach to delineate proteomic alterations in oral cancer patients based on their history of tobacco using habits (patients who chewed tobacco, patients who smoked tobacco, and those with no history of tobacco consumption). Our data identified distinct dysregulation of biological processes and pathways in each patient cohort. Bioinformatics analysis of dysregulated proteins identified in our proteomic study revealed dysregulation of collagen formation and antigen processing/presentation pathway in oral cancer patients who smoked tobacco, whereas proteins associated with the process of keratinization showed enrichment in patients who chewed tobacco. In addition, we identified overexpression of proteins involved in immune pathways and downregulation of muscle contraction-mediated signaling events in all three cohorts, irrespective of tobacco using habits. This study lays the groundwork for identification of protein markers that may aid in identification of high-risk patients for cancer development based on the history of tobacco exposure habits.