Identification of novel alternative splice isoforms of circulating proteins in a mouse model of human pancreatic cancer

Physiology of malate dehydrogenase and how dysregulation leads to disease

Malate dehydrogenase (MDH) is pivotal in mammalian tissue metabolism, participating in various pathways beyond its classical roles and highlighting its adaptability to cellular demands. This enzyme is involved in maintaining redox balance, lipid synthesis, and glutamine metabolism and supports rapidly proliferating cells' energetic and biosynthetic needs. The involvement of MDH in glutamine metabolism underlines its significance in cell physiology. In contrast, its contribution to lipid metabolism highlights its role in essential biosynthetic processes necessary for cell maintenance and proliferation. The enzyme's regulatory mechanisms, such as post-translational modifications, underscore its complexity and importance in metabolic regulation, positioning MDH as a potential target in metabolic dysregulation. Furthermore, the association of MDH with various pathologies, including cancer and neurological disorders, suggests its involvement in disease progression. The overexpression of MDH isoforms MDH1 and MDH2 in cancers like breast, prostate, and pancreatic ductal adenocarcinoma, alongside structural modifications, implies their critical role in the metabolic adaptation of tumor cells. Additionally, mutations in MDH2 linked to pheochromocytomas, paragangliomas, and other metabolic diseases emphasize MDH's role in metabolic homeostasis. This review spotlights MDH's potential as a biomarker and therapeutic target, advocating for further research into its multifunctional roles and regulatory mechanisms in health and disease.

ASV-ID, a Proteogenomic Workflow To Predict Candidate Protein Isoforms on the Basis of Transcript Evidence

One of the goals of the Chromosome-Centric Human Proteome Project (C-HPP) is to map and characterize the functions of protein isoforms produced by alternative splicing of genes. However, identifying alternative splice variants (ASVs) via mass spectrometry remains a major challenge, because ASVs usually contain highly homologous peptide sequences. A routine protein sequence analysis suggests that more than half of the investigated proteins do not generate two or more uniquely mapping peptides that would enable their isoforms to be distinguished. Here, we develop a new proteogenomics method, named "ASV-ID" (alternative splicing variants identification), which enables identification of ASVs by using a cell type-specific protein sequence database that is supported by RNA-Seq data. Using this workflow, we identify 1935 distinct proteins under highly stringent conditions. In fact, transcript evidence on these 841 proteins helps us distinguish them from other isoforms, despite the fact that these proteins are not predicted to make 2 or more uniquely mapping peptides. We also demonstrate that ASV-ID enables detection of 19 differently expressed isoforms present in several cell lines. Thus, a new workflow using ASV-ID has the potential to map yet-to-be-identified difficult protein isoforms in a simple and robust way.

Splice variants as novel targets in pancreatic ductal adenocarcinoma

Despite a wealth of genomic information, a comprehensive alternative splicing (AS) analysis of pancreatic ductal adenocarcinoma (PDAC) has not been performed yet. In the present study, we assessed whole exome-based transcriptome and AS profiles of 43 pancreas tissues using Affymetrix exon array. The AS analysis of PDAC indicated on average two AS probe-sets (ranging from 1-28) in 1,354 significantly identified protein-coding genes, with skipped exon and alternative first exon being the most frequently utilised. In addition to overrepresented extracellular matrix (ECM)-receptor interaction and focal adhesion that were also seen in transcriptome differential expression (DE) analysis, Fc gamma receptor-mediated phagocytosis and axon guidance AS genes were also highly represented. Of note, the highest numbers of AS probe-sets were found in collagen genes, which encode the characteristically abundant stroma seen in PDAC. We also describe a set of 37 'hypersensitive' genes which were frequently targeted by somatic mutations, copy number alterations, DE and AS, indicating their propensity for multidimensional regulation. We provide the most comprehensive overview of the AS landscape in PDAC with underlying changes in the spliceosomal machinery. We also collate a set of AS and DE genes encoding cell surface proteins, which present promising diagnostic and therapeutic targets in PDAC.

Identification of Alternative Splice Variants Using Unique Tryptic Peptide Sequences for Database Searches

Alternative splicing is a mechanism in eukaryotes by which different forms of mRNAs are generated from the same gene. Identification of alternative splice variants requires the identification of peptides specific for alternative splice forms. For this purpose, we generated a human database that contains only unique tryptic peptides specific for alternative splice forms from Swiss-Prot entries. Using this database allows an easy access to splice variant-specific peptide sequences that match to MS data. Furthermore, we combined this database without alternative splice variant-1-specific peptides with human Swiss-Prot. This combined database can be used as a general database for searching of LC-MS data. LC-MS data derived from in-solution digests of two different cell lines (LNCaP, HeLa) and phosphoproteomics studies were analyzed using these two databases. Several nonalternative splice variant-1-specific peptides were found in both cell lines, and some of them seemed to be cell-line-specific. Control and apoptotic phosphoproteomes from Jurkat T cells revealed several nonalternative splice variant-1-specific peptides, and some of them showed clear quantitative differences between the two states.

Tutorial on Protein Ontology Resources

The Protein Ontology (PRO) is the reference ontology for proteins in the Open Biomedical Ontologies (OBO) foundry and consists of three sub-ontologies representing protein classes of homologous genes, proteoforms (e.g., splice isoforms, sequence variants, and post-translationally modified forms), and protein complexes. PRO defines classes of proteins and protein complexes, both species-specific and species nonspecific, and indicates their relationships in a hierarchical framework, supporting accurate protein annotation at the appropriate level of granularity, analyses of protein conservation across species, and semantic reasoning. In the first section of this chapter, we describe the PRO framework including categories of PRO terms and the relationship of PRO to other ontologies and protein resources. Next, we provide a tutorial about the PRO website ( proconsortium.org ) where users can browse and search the PRO hierarchy, view reports on individual PRO terms, and visualize relationships among PRO terms in a hierarchical table view, a multiple sequence alignment view, and a Cytoscape network view. Finally, we describe several examples illustrating the unique and rich information available in PRO.

Using Galaxy for Proteomics

The area of informatics for mass spectrometry (MS)-based proteomics data has steadily grown over the last two decades. Numerous, effective software programs now exist for various aspects of proteomic informatics. However, many researchers still have difficulties in using these software. These difficulties arise from problems with running and integrating disparate software programs, scalability issues when dealing with large data volumes, and lack of ability to share and reproduce workflows comprised of different software. The Galaxy framework for bioinformatics provides an attractive option for solving many of these current issues in proteomic informatics. Originally developed as a workbench to enable genomic data analysis, numerous researchers are now turning to Galaxy to implement software for MS-based proteomics applications. Here, we provide an introduction to Galaxy and its features, and describe how software tools are deployed, published and shared via the scalable framework. We also describe some of the existing tools in Galaxy for basic MS-based proteomics data analysis and informatics. Finally, we describe how proteomics tools in Galaxy can be combined with other existing tools for genomic and transcriptomic data analysis to enable powerful multi-omic data analysis applications.

Proteogenomics: Integrating Next-Generation Sequencing and Mass Spectrometry to Characterize Human Proteomic Variation

Mass spectrometry-based proteomics has emerged as the leading method for detection, quantification, and characterization of proteins. Nearly all proteomic workflows rely on proteomic databases to identify peptides and proteins, but these databases typically contain a generic set of proteins that lack variations unique to a given sample, precluding their detection. Fortunately, proteogenomics enables the detection of such proteomic variations and can be defined, broadly, as the use of nucleotide sequences to generate candidate protein sequences for mass spectrometry database searching. Proteogenomics is experiencing heightened significance due to two developments: (a) advances in DNA sequencing technologies that have made complete sequencing of human genomes and transcriptomes routine, and (b) the unveiling of the tremendous complexity of the human proteome as expressed at the levels of genes, cells, tissues, individuals, and populations. We review here the field of human proteogenomics, with an emphasis on its history, current implementations, the types of proteomic variations it reveals, and several important applications.

Unveiling alterative splice diversity from human oligodendrocyte proteome data

Oligodendrocytes produce and maintain the myelin sheath of axons in the central nervous system. Because misassembled myelin sheaths have been associated with brain disorders such as multiple sclerosis and schizophrenia, recent advances have been made towards the description of the oligodendrocyte proteome. The identification of splice variants represented in the proteome is as important as determining the level of oligodendrocyte-associated proteins. Here, we used an oligodendrocyte proteome dataset deposited in ProteomeXchange to search against a customized protein sequence file containing computationally predicted splice variants. Our approach resulted in the identification of 39 splice variants, including one variant from the GTPase KRAS gene and another from the human glutaminase gene family. We also detected the mRNA expression of five selected splice variants and demonstrated that a fraction of these have their canonical proteins participating in direct protein-protein interactions. In conclusion, we believe our findings contribute to the molecular characterization of oligodendrocytes and may encourage other research groups working with central nervous system disorders to investigate the biological significance of these splice variants. The splice variants identified in this study may encode proteins that could be targeted in novel treatment strategies and diagnostic methods.

SIGNIFICANCE

Several disorders of the central nervous system (CNS) are associated with misassembled myelin sheaths, which are produced and maintained by oligodendrocytes (OL). Recently, the OL proteome has been explored to identify key proteins and molecular functions associated with CNS disorders. We developed an innovative approach to select, with a higher level of confidence, a relevant list of splice variants from a proteome dataset and detected the mRNA expression of five selected variants: EEF1D, KRAS, MFF, SDR39U1, and SUGT1. We also described splice variants extracted from OL proteome data. Among the splice variants identified, some are from genes previously linked to CNS and related disorders. Our findings may contribute to oligodendrocyte characterization and encourage other research groups to investigate the biological role of splice variants and to improve current treatments and diagnostic methods for CNS disorders.

Decreased Serum Thrombospondin-1 Levels in Pancreatic Cancer Patients Up to 24 Months Prior to Clinical Diagnosis: Association with Diabetes Mellitus

PURPOSE

Identification of serum biomarkers enabling earlier diagnosis of pancreatic ductal adenocarcinoma (PDAC) could improve outcome. Serum protein profiles in patients with preclinical disease and at diagnosis were investigated.

EXPERIMENTAL DESIGN

Serum from cases up to 4 years prior to PDAC diagnosis and controls (UKCTOCS,n= 174) were studied, alongside samples from patients diagnosed with PDAC, chronic pancreatitis, benign biliary disease, type 2 diabetes mellitus, and healthy subjects (n= 298). Isobaric tags for relative and absolute quantification (iTRAQ) enabled comparisons of pooled serum from a test set (n= 150). Validation was undertaken using multiple reaction monitoring (MRM) and/or Western blotting in all 472 human samples and samples from a KPC mouse model.

RESULTS

iTRAQ identified thrombospondin-1 (TSP-1) as reduced preclinically and in diagnosed samples. MRM confirmed significant reduction in levels of TSP-1 up to 24 months prior to diagnosis. A combination of TSP-1 and CA19-9 gave an AUC of 0.86, significantly outperforming both markers alone (0.69 and 0.77, respectively;P< 0.01). TSP-1 was also decreased in PDAC patients compared with healthy controls (P< 0.05) and patients with benign biliary obstruction (P< 0.01). Low levels of TSP-1 correlated with poorer survival, preclinically (P< 0.05) and at clinical diagnosis (P< 0.02). In PDAC patients, reduced TSP-1 levels were more frequently observed in those with confirmed diabetes mellitus (P< 0.01). Significantly lower levels were also observed in PDAC patients with diabetes compared with individuals with type 2 diabetes mellitus (P= 0.01).

CONCLUSIONS

Circulating TSP-1 levels decrease up to 24 months prior to diagnosis of PDAC and significantly enhance the diagnostic performance of CA19-9. The influence of diabetes mellitus on biomarker behavior should be considered in future studies.

Glucose metabolic phenotype of pancreatic cancer

AIM: To construct a global “metabolic phenotype” of pancreatic ductal adenocarcinoma (PDAC) reflecting tumour-related metabolic enzyme expression.

METHODS: A systematic review of the literature was performed using OvidSP and PubMed databases using keywords “pancreatic cancer” and individual glycolytic and mitochondrial oxidative phosphorylation (MOP) enzymes. Both human and animal studies investigating the oncological effect of enzyme expression changes and inhibitors in both an in vitro and in vivo setting were included in the review. Data reporting changes in enzyme expression and the effects on PDAC cells, such as survival and metastatic potential, were extracted to construct a metabolic phenotype.

RESULTS: Seven hundred and ten papers were initially retrieved, and were screened to meet the review inclusion criteria. 107 unique articles were identified as reporting data involving glycolytic enzymes, and 28 articles involving MOP enzymes in PDAC. Data extraction followed a pre-defined protocol. There is consistent over-expression of glycolytic enzymes and lactate dehydrogenase in keeping with the Warburg effect to facilitate rapid adenosine-triphosphate production from glycolysis. Certain isoforms of these enzymes were over-expressed specifically in PDAC. Altering expression levels of HK, PGI, FBA, enolase, PK-M2 and LDA-A with metabolic inhibitors have shown a favourable effect on PDAC, thus identifying these as potential therapeutic targets. However, the Warburg effect on MOP enzymes is less clear, with different expression levels at different points in the Krebs cycle resulting in a fundamental change of metabolite levels, suggesting that other essential anabolic pathways are being stimulated.

CONCLUSION: Further characterisation of the PDAC metabolic phenotype is necessary as currently there are few clinical studies and no successful clinical trials targeting metabolic enzymes.

Overexpression of Rad51C splice variants in colorectal tumors

Functional alterations in Rad51C are the cause of the Fanconi anemia complementation group O (FANCO) gene disorder. We have identified novel splice variants of Rad51C mRNA in colorectal tumors and cells. The alternatively spliced transcript variants are formed either without exon-7 (variant 1), without exon 6 and 7 (variant 2) or without exon 7 and 8 (variant 3).

Real time PCR analysis of nine pair-matched colorectal tumors and non-tumors showed that variant 1 was overexpressed in tumors compared to matched non-tumors. Among 38 colorectal tumor RNA samples analyzed, 18 contained variant 1, 12 contained variant 2, 14 contained variant 3, and eight expressed full length Rad51C exclusively. Bisulfite DNA sequencing showed promoter methylation of Rad51C in tumor cells. 5-azacytidine treatment of LS-174T cells caused a 14 fold increase in variant 1, a 4.8 fold increase for variant 3 and 3.4 fold for variant 2 compared to 2.5 fold increase in WT.

Expression of Rad51C variants is associated with FANCD2 foci positive colorectal tumors and is associated with microsatellite stability in those tumors. Further investigation is needed to elucidate differential function of the Rad51C variants to evaluate potential effects in drug resistance and DNA repair.

Proteogenomics for understanding oncology: recent advances and future prospects

The concept of proteogenomics has emerged rapidly as a valuable approach to integrate mass spectrometry-derived proteomic data with genomic and transcriptomic data. It is used to harness the full potential of the former dataset in the discovery of potential biomarkers, therapeutic targets and novel proteins associated with various biological processes including diseases. Proteogenomic strategies have been successfully utilized to identify novel genes and redefine annotation of existing gene models in various genomes. In recent years, this approach has been extended to the field of cancer biology to unravel complexities in the tumor genomes and proteomes. Standard proteomics workflows employing translated cancer genomes and transcriptomes can potentially identify peptides from mutant proteins, splice variants and fusion proteins in the tumor proteome, which in addition to the currently available biomarker panels can serve as potential diagnostic and prognostic biomarkers, besides having therapeutic utility. This review focuses on the role of proteogenomics to understand cancer biology.

GenomewidePDB 2.0: A Newly Upgraded Versatile Proteogenomic Database for the Chromosome-Centric Human Proteome Project

Since the launch of the Chromosome-centric Human Proteome Project (C-HPP) in 2012, the number of "missing" proteins has fallen to 2932, down from ∼5932 since the number was first counted in 2011. We compared the characteristics of missing proteins with those of already annotated proteins with respect to transcriptional expression pattern and the time periods in which newly identified proteins were annotated. We learned that missing proteins commonly exhibit lower levels of transcriptional expression and less tissue-specific expression compared with already annotated proteins. This makes it more difficult to identify missing proteins as time goes on. One of the C-HPP goals is to identify alternative spliced product of proteins (ASPs), which are usually difficult to find by shot-gun proteomic methods due to their sequence similarities with the representative proteins. To resolve this problem, it may be necessary to use a targeted proteomics approach (e.g., selected and multiple reaction monitoring [S/MRM] assays) and an innovative bioinformatics platform that enables the selection of target peptides for rarely expressed missing proteins or ASPs. Given that the success of efforts to identify missing proteins may rely on more informative public databases, it was necessary to upgrade the available integrative databases. To this end, we attempted to improve the features and utility of GenomewidePDB by integrating transcriptomic information (e.g., alternatively spliced transcripts), annotated peptide information, and an advanced search interface that can find proteins of interest when applying a targeted proteomics strategy. This upgraded version of the database, GenomewidePDB 2.0, may not only expedite identification of the remaining missing proteins but also enhance the exchange of information among the proteome community. GenomewidePDB 2.0 is available publicly at http://genomewidepdb.proteomix.org/.

Alternatively Spliced Homologous Exons Have Ancient Origins and Are Highly Expressed at the Protein Level

Alternative splicing of messenger RNA can generate a wide variety of mature RNA transcripts, and these transcripts may produce protein isoforms with diverse cellular functions. While there is much supporting evidence for the expression of alternative transcripts, the same is not true for the alternatively spliced protein products. Large-scale mass spectroscopy experiments have identified evidence of alternative splicing at the protein level, but with conflicting results. Here we carried out a rigorous analysis of the peptide evidence from eight large-scale proteomics experiments to assess the scale of alternative splicing that is detectable by high-resolution mass spectroscopy. We find fewer splice events than would be expected: we identified peptides for almost 64% of human protein coding genes, but detected just 282 splice events. This data suggests that most genes have a single dominant isoform at the protein level. Many of the alternative isoforms that we could identify were only subtly different from the main splice isoform. Very few of the splice events identified at the protein level disrupted functional domains, in stark contrast to the two thirds of splice events annotated in the human genome that would lead to the loss or damage of functional domains. The most striking result was that more than 20% of the splice isoforms we identified were generated by substituting one homologous exon for another. This is significantly more than would be expected from the frequency of these events in the genome. These homologous exon substitution events were remarkably conserved—all the homologous exons we identified evolved over 460 million years ago—and eight of the fourteen tissue-specific splice isoforms we identified were generated from homologous exons. The combination of proteomics evidence, ancient origin and tissue-specific splicing indicates that isoforms generated from homologous exons may have important cellular roles.

Alternative splicing is thought to be one means for generating the protein diversity necessary for the whole range of cellular functions. While the presence of alternatively spliced transcripts in the cell has been amply demonstrated, the same cannot be said for alternatively spliced proteins. The quest for alternative protein isoforms has focused primarily on the analysis of peptides from large-scale mass spectroscopy experiments, but evidence for alternative isoforms has been patchy and contradictory. A careful analysis of the peptide evidence is needed to fully understand the scale of alternative splicing detectable at the protein level. Here we analysed peptides from eight large-scale data sets, identifying just 282 splice events among 12,716 genes. This suggests that most genes have a single dominant isoform. Many of the alternative isoforms that we identified were only subtly different from the main splice variant, and one in five was generated by substitution of homologous exons by swapping one related exon for another. Remarkably, the alternative isoforms generated from homologous exons were highly conserved, first appearing 460 million years ago, and several appear to have tissue-specific roles in the brain and heart. Our results suggest that these particular isoforms are likely to have important cellular roles.

Proteomics analysis of bodily fluids in pancreatic cancer

Proteomics study of pancreatic cancer using bodily fluids emphasizes biomarker discovery and clinical application, presenting unique prospect and challenges. Depending on the physiological nature of the bodily fluid and its proximity to pancreatic cancer, the proteomes of bodily fluids, such as pancreatic juice, pancreatic cyst fluid, blood, bile, and urine, can be substantially different in terms of protein constitution and the dynamic range of protein concentration. Thus, a comprehensive discovery and specific detection of cancer-associated proteins within these varied fluids is a complex task, requiring rigorous experiment design and a concerted approach. While major challenges still remain, fluid proteomics studies in pancreatic cancer to date have provided a wealth of information in revealing proteome alterations associated with pancreatic cancer in various bodily fluids.

Identification of sialylated glycoproteins from metabolically oligosaccharide engineered pancreatic cells

In this study, we investigated the use of metabolic oligosaccharide engineering and bio-orthogonal ligation reactions combined with lectin microarray and mass spectrometry to analyze sialoglycoproteins in the SW1990 human pancreatic cancer line. Specifically, cells were treated with the azido N-acetylmannosamine analog, 1,3,4-Bu3ManNAz, to label sialoglycoproteins with azide-modified sialic acids. The metabolically labeled sialoglyproteins were then biotinylated via the Staudinger ligation, and sialoglycopeptides containing azido-sialic acid glycans were immobilized to a solid support. The peptides linked to metabolically labeled sialylated glycans were then released from sialoglycopeptides and analyzed by mass spectrometry; in parallel, the glycans from azido-sialoglycoproteins were characterized by lectin microarrays. This method identified 75 unique N-glycosite-containing peptides from 55 different metabolically labeled sialoglycoproteins of which 42 were previously linked to cancer in the literature. A comparison of two of these glycoproteins, LAMP1 and ORP150, in histological tumor samples showed overexpression of these proteins in the cancerous tissue demonstrating that our approach constitutes a viable strategy to identify and discover sialoglycoproteins associated with cancer, which can serve as biomarkers for cancer diagnosis or targets for therapy.

Metabolic rewiring of pancreatic ductal adenocarcinoma: New routes to follow within the maze

Pancreatic ductal adenocarcinoma (PDAC) is a debilitating and almost universally fatal malignancy. Despite advances in understanding of the oncogenetics of the disease, very few clinical benefits have been shown. One of the main characteristics of PDAC is the tumor architecture where tumor cells are surrounded by a firm desmoplasia. By reducing vascularization, thus both oxygen and nutrients delivery to the tumor, this stroma causes the appearance of hypoxic zones driving metabolic adaptation in surviving tumor cells in order to cope with challenging conditions. This metabolic reprogramming promoted by environmental constraints enhances PDAC aggressiveness. In this review, we provide a brief overview of previous works regarding the importance of glucose and glutamine addiction of PDAC cells. In particular we aim to highlight the need for exploring the impact of metabolites other than glucose and glutamine, such as non-essential amino acids and oncometabolites, to find new treatments. We also discuss the need for progress in methodology for metabolites detection. The overall purpose of our review is to emphasize the need to look beyond what is currently known, with a focus on amino acid availability, in order to improve our understanding of PDAC biology.

Proteogenomics: concepts, applications and computational strategies

Proteogenomics is an area of research at the interface of proteomics and genomics. In this approach, customized protein sequence databases generated using genomic and transcriptomic information are used to help identify novel peptides (not present in reference protein sequence databases) from mass spectrometry-based proteomic data; in turn, the proteomic data can be used to provide protein-level evidence of gene expression and to help refine gene models. In recent years, owing to the emergence of new sequencing technologies such as RNA-seq and dramatic improvements in the depth and throughput of mass spectrometry-based proteomics, the pace of proteogenomic research has greatly accelerated. Here I review the current state of proteogenomic methods and applications, including computational strategies for building and using customized protein sequence databases. I also draw attention to the challenge of false positive identifications in proteogenomics and provide guidelines for analyzing the data and reporting the results of proteogenomic studies.

Onco-proteogenomics: cancer proteomics joins forces with genomics

The complexities of tumor genomes are rapidly being uncovered, but how they are regulated into functional proteomes remains poorly understood. Standard proteomics workflows use databases of known proteins, but these databases do not capture the uniqueness of the cancer transcriptome, with its point mutations, unusual splice variants and gene fusions. Onco-proteogenomics integrates mass spectrometry-generated data with genomic information to identify tumor-specific peptides. Linking tumor-derived DNA, RNA and protein measurements into a central-dogma perspective has the potential to improve our understanding of cancer biology.

Mapping Splicing Quantitative Trait Loci in RNA-Seq

BACKGROUND

One of the major mechanisms of generating mRNA diversity is alternative splicing, a regulated process that allows for the flexibility of producing functionally different proteins from the same genomic sequences. This process is often altered in cancer cells to produce aberrant proteins that drive the progression of cancer. A better understanding of the misregulation of alternative splicing will shed light on the development of novel targets for pharmacological interventions of cancer.

METHODS

In this study, we evaluated three statistical methods, random effects meta-regression, beta regression, and generalized linear mixed effects model, for the analysis of splicing quantitative trait loci (sQTL) using RNA-Seq data. All the three methods use exon-inclusion levels estimated by the PennSeq algorithm, a statistical method that utilizes paired-end reads and accounts for non-uniform sequencing coverage.

RESULTS

Using both simulated and real RNA-Seq datasets, we compared these three methods with GLiMMPS, a recently developed method for sQTL analysis. Our results indicate that the most reliable and powerful method was the random effects meta-regression approach, which identified sQTLs at low false discovery rates but higher power when compared to GLiMMPS.

CONCLUSIONS

We have evaluated three statistical methods for the analysis of sQTLs in RNA-Seq. Results from our study will be instructive for researchers in selecting the appropriate statistical methods for sQTL analysis.

Transcriptomic analysis of Rongchang pig brains and livers

Recent developments in high-throughput RNA sequencing (RNA-seq) technology have led to a dramatic impact on our understanding of the structure and expression profiles of the mammalian transcriptome. To gain insights into the usefulness of swine production and biomedical model, the transcriptome profiling of Rongchang pig brains and livers was characterized using RNA-seq technology to uncover functional candidate molecules. In the study, total RNAs from brains and livers of Rongchang pig were sequenced and 8.6Gb sequencing data was obtained. This analysis revealed tissue specificity through the identification of 5575 and 4600 differentially expressed genes (DEGs) in brains and livers, respectively and the functional analysis of DEGs. Furthermore, 83 neuropeptide gene transcripts, 69 neuropeptide receptor gene transcripts, 10 pro-neuropeptide convertase gene transcripts and many other neuropeptide related protein gene transcripts were identified. Totally, the major characteristics of the transcriptional profiles of Rongchang pig brains and livers were present.

Revisiting the identification of canonical splice isoforms through integration of functional genomics and proteomics evidence

Canonical isoforms in different databases have been defined as the most prevalent, most conserved, most expressed, longest, or the one with the clearest description of domains or posttranslational modifications. In this article, we revisit these definitions of canonical isoforms based on functional genomics and proteomics evidence, focusing on mouse data. We report a novel functional relationship network-based approach for identifying the highest connected isoforms (HCIs). We show that 46% of these HCIs are not the longest transcripts. In addition, this approach revealed many genes that have more than one highly connected isoforms. Averaged across 175 RNA-seq datasets covering diverse tissues and conditions, 65% of the HCIs show higher expression levels than nonhighest connected isoforms at the transcript level. At the protein level, these HCIs highly overlap with the expressed splice variants, based on proteomic data from eight different normal tissues. These results suggest that a more confident definition of canonical isoforms can be made through integration of multiple lines of evidence, including HCIs defined by biological processes and pathways, expression prevalence at the transcript level, and relative or absolute abundance at the protein level. This integrative proteogenomics approach can successfully identify principal isoforms that are responsible for the canonical functions of genes.

Profiling proteoforms: promising follow-up of proteomics for biomarker discovery

Today, proteomics usually compares clinical samples by use of bottom-up profiling with high resolution mass spectrometry, where all protein products of a single gene are considered as an integral whole. At the same time, proteomics of proteoforms, which considers the variety of protein species, offers the potential to discover valuable biomarkers. Proteoforms are protein species that arise as a consequence of genetic polymorphisms, alternative splicing, post-translational modifications and other less-explored molecular events. The comprehensive observation of proteoforms has been an exclusive privilege of top-down proteomics. Here, we review the possibilities of a bottom-up approach to address the microheterogeneity of the human proteome. Special focus is given to shotgun proteomics and structure-based bioinformatics as a source of hypothetical proteoforms, which can potentially be verified by targeted mass spectrometry to determine the relevance of proteoforms to diseases.

Mapping Splicing Quantitative Trait Loci in RNA-Seq

BACKGROUND

One of the major mechanisms of generating mRNA diversity is alternative splicing, a regulated process that allows for the flexibility of producing functionally different proteins from the same genomic sequences. This process is often altered in cancer cells to produce aberrant proteins that drive the progression of cancer. A better understanding of the misregulation of alternative splicing will shed light on the development of novel targets for pharmacological interventions of cancer.

METHODS

In this study, we evaluated three statistical methods, random effects meta-regression, beta regression, and generalized linear mixed effects model, for the analysis of splicing quantitative trait loci (sQTL) using RNA-Seq data. All the three methods use exon-inclusion levels estimated by the PennSeq algorithm, a statistical method that utilizes paired-end reads and accounts for non-uniform sequencing coverage.

RESULTS

Using both simulated and real RNA-Seq datasets, we compared these three methods with GLiMMPS, a recently developed method for sQTL analysis. Our results indicate that the most reliable and powerful method was the random effects meta-regression approach, which identified sQTLs at low false discovery rates but higher power when compared to GLiMMPS.

CONCLUSIONS

We have evaluated three statistical methods for the analysis of sQTLs in RNA-Seq. Results from our study will be instructive for researchers in selecting the appropriate statistical methods for sQTL analysis.

The emerging era of genomic data integration for analyzing splice isoform function

The vast majority of multi-exon genes in humans undergo alternative splicing, which greatly increases the functional diversity of protein species. Predicting functions at the isoform level is essential to further our understanding of developmental abnormalities and cancers, which frequently exhibit aberrant splicing and dysregulation of isoform expression. However, determination of isoform function is very difficult, and efforts to predict isoform function have been limited in the functional genomics field. Deep sequencing of RNA now provides an unprecedented amount of expression data at the transcript level. We describe here emerging computational approaches that integrate such large-scale whole-transcriptome sequencing (RNA-seq) data for predicting the functions of alternatively spliced isoforms, and we discuss their applications in developmental and cancer biology. We outline future directions for isoform function prediction, emphasizing the need for heterogeneous genomic data integration and tissue-specific, dynamic isoform-level network modeling, which will allow the field to realize its full potential.

Integrating genomic, transcriptomic, and interactome data to improve Peptide and protein identification in shotgun proteomics

Mass spectrometry (MS)-based shotgun proteomics is an effective technology for global proteome profiling. The ultimate goal is to assign tandem MS spectra to peptides and subsequently infer proteins and their abundance. In addition to database searching and protein assembly algorithms, computational approaches have been developed to integrate genomic, transcriptomic, and interactome information to improve peptide and protein identification. Earlier efforts focus primarily on making databases more comprehensive using publicly available genomic and transcriptomic data. More recently, with the increasing affordability of the Next Generation Sequencing (NGS) technologies, personalized protein databases derived from sample-specific genomic and transcriptomic data have emerged as an attractive strategy. In addition, incorporating interactome data not only improves protein identification but also puts identified proteins into their functional context and thus facilitates data interpretation. In this paper, we survey the major integrative bioinformatics approaches that have been developed during the past decade and discuss their merits and demerits.

A new class of protein cancer biomarker candidates: differentially expressed splice variants of ERBB2 (HER2/neu) and ERBB1 (EGFR) in breast cancer cell lines

Combined RNA-Seq and proteomics analyses reveal striking differential expression of splice isoforms of key proteins in important cancer pathways and networks. Even between primary tumor cell lines from histologically similar inflammatory breast cancers, we find striking differences in hormone receptor-negative cell lines that are ERBB2 (Her2/neu)-amplified versus ERBB1 (EGFR) over-expressed with low ERBB2 activity. We have related these findings to protein-protein interaction networks, signaling and metabolic pathways, and methods for predicting functional variants among multiple alternative isoforms. Understanding the upstream ligands and regulators and the downstream pathways and interaction networks for ERBB receptors is certain to be important for explanation and prediction of the variable levels of expression and therapeutic responses of ERBB+tumors in the breast and in other organ sites. Alternative splicing is a remarkable evolutionary development that increases protein diversity from multi-exonic genes without requiring expansion of the genome. It is no longer sufficient to report the up- or down-expression of genes and proteins without dissecting the complexity due to alternative splicing. This article is part of a Special Issue entitled: 20Years of Proteomics in memory of Viatliano Pallini. Guest Editors: Luca Bini , Juan J. Calvete, Natacha Turck, Denis Hochstrasser and Jean-Charles Sanchez.

Strategies for Genomic and Proteomic Profiling of Cancers

Omics-based technology platforms have made new kinds of cancer profiling tests feasible. There are several valuable examples in clinical practice, and many more under development. A concerted, transparent process of discovery with lock-down of candidate assays and classifiers and clear specification of intended clinical use is essential. The Institute of Medicine has now proposed a three-stage scheme of confirming and validating analytical findings, validating performance on clinical specimens, and demonstrating explicit clinical utility for an approvable test (Micheel et al., Evolution of translational omics: lessons learned and path forward, 2012).

Moving from unsequenced to sequenced genome: reanalysis of the proteome of Leishmania donovani

The kinetoplastid protozoan parasite, Leishmania donovani, is the causative agent of kala azar or visceral leishmaniasis. Kala azar is a severe form of leishmaniasis that is fatal in the majority of untreated cases. Studies on proteomic analysis of L. donovani thus far have been carried out using homology-based identification based on related Leishmania species (L. infantum, L. major and L. braziliensis) whose genomes have been sequenced. Recently, the genome of L. donovani was fully sequenced and the data became publicly available. We took advantage of the availability of its genomic sequence to carry out a more accurate proteogenomic analysis of L. donovani proteome using our previously generated dataset. This resulted in identification of 17,504 unique peptides upon database-dependent search against the annotated proteins in L. donovani. These peptides were assigned to 3999 unique proteins in L. donovani. 2296 proteins were identified in both the life stages of L. donovani, while 613 and 1090 proteins were identified only from amastigote and promastigote stages, respectively. The proteomic data was also searched against six-frame translated L. donovani genome, which led to 255 genome search-specific peptides (GSSPs) resulting in identification of 20 novel genes and correction of 40 existing gene models in L. donovani.

BIOLOGICAL SIGNIFICANCE

Leishmania donovani genome sequencing was recently completed, which permitted us to use a proteogenomic approach to map its proteome and to carry out annotation of it genome. This resulted in mapping of 50% (3999 proteins) of L. donovani proteome. Our study identified 20 novel genes previously not predicted from the L. donovani genome in addition to correcting annotations of 40 existing gene models. The identified proteins may help in better understanding of stage-specific protein expression profiles in L. donovani and to identify novel stage-specific drug targets in L. donovani which could be used in the treatment of leishmaniasis. This article is part of a Special Issue entitled: Trends in Microbial Proteomics.

Resolving chromosome-centric human proteome with translating mRNA analysis: a strategic demonstration

Chromosome-centric human proteome project (C-HPP) aims at differentiating chromosome-based and tissue-specific protein compositions in terms of protein expression, quantification, and modification. We previously found that the analysis of translating mRNA (mRNA attached to ribosome-nascent chain complex, RNC-mRNA) can explain over 94% of mRNA-protein abundance. Therefore, we propose here to use full-length RNC-mRNA information to illustrate protein expression both qualitatively and quantitatively. We performed RNA-seq on RNC-mRNA (RNC-seq) and detected 12,758 and 14,113 translating genes in human normal bronchial epithelial (HBE) cells and human colorectal adenocarcinoma Caco-2 cells, respectively. We found that most of these genes were mapped with >80% of coding sequence coverage. In Caco-2 cells, we provided translating evidence on 4180 significant single-nucleotide variations. While using RNC-mRNA data as a standard for proteomic data integration, both translating and protein evidence of 7876 genes can be acquired from four interlaboratory data sets with different MS platforms. In addition, we detected 1397 noncoding mRNAs that were attached to ribosomes, suggesting a potential source of new protein explorations. By comparing the two cell lines, a total of 677 differentially translated genes were found to be nonevenly distributed across chromosomes. In addition, 2105 genes in Caco-2 and 750 genes in HBE cells are expressed in a cell-specific manner. These genes are significantly and specifically clustered on multiple chromosomes, such as chromosome 19. We conclude that HPP/C-HPP investigations can be considerably improved by integrating RNC-mRNA analysis with MS, bioinformatics, and antibody-based verifications.

Tools to covisualize and coanalyze proteomic data with genomes and transcriptomes: validation of genes and alternative mRNA splicing

Direct links between proteomic and genomic/transcriptomic data are not frequently made, partly because of lack of appropriate bioinformatics tools. To help address this, we have developed the PG Nexus pipeline. The PG Nexus allows users to covisualize peptides in the context of genomes or genomic contigs, along with RNA-seq reads. This is done in the Integrated Genome Viewer (IGV). A Results Analyzer reports the precise base position where LC-MS/MS-derived peptides cover genes or gene isoforms, on the chromosomes or contigs where this occurs. In prokaryotes, the PG Nexus pipeline facilitates the validation of genes, where annotation or gene prediction is available, or the discovery of genes using a "virtual protein"-based unbiased approach. We illustrate this with a comprehensive proteogenomics analysis of two strains of Campylobacter concisus . For higher eukaryotes, the PG Nexus facilitates gene validation and supports the identification of mRNA splice junction boundaries and splice variants that are protein-coding. This is illustrated with an analysis of splice junctions covered by human phosphopeptides, and other examples of relevance to the Chromosome-Centric Human Proteome Project. The PG Nexus is open-source and available from https://github.com/IntersectAustralia/ap11_Samifier. It has been integrated into Galaxy and made available in the Galaxy tool shed.

Plasma proteomics, the Human Proteome Project, and cancer-associated alternative splice variant proteins

This article addresses three inter-related subjects: the development of the Human Plasma Proteome Peptide Atlas, the launch of the Human Proteome Project, and the emergence of alternative splice variant transcripts and proteins as important features of evolution and pathogenesis. The current Plasma Peptide Atlas provides evidence on which peptides have been detected for every protein confidently identified in plasma; there are links to their spectra and their estimated abundance, facilitating the planning of targeted proteomics for biomarker studies. The Human Proteome Project (HPP) combines a chromosome-centric C-HPP with a biology and disease-driven B/D-HPP, upon a foundation of mass spectrometry, antibody, and knowledgebase resource pillars. The HPP aims to identify the approximately 7000 "missing proteins" and to characterize all proteins and their many isoforms. Success will enable the larger research community to utilize newly-available peptides, spectra, informative MS transitions, and databases for targeted analyses of priority proteins for each organ and disease. Among the isoforms of proteins, splice variants have the special feature of greatly enlarging protein diversity without enlarging the genome; evidence is accumulating of striking differential expression of splice variants in cancers. In this era of RNA-sequencing and advanced mass spectrometry, it is no longer sufficient to speak simply of increased or decreased expression of genes or proteins without carefully examining the splice variants in the protein mixture produced from each multi-exon gene. This article is part of a Special Issue entitled: Biomarkers: A Proteomic Challenge.

Distinct splice variants and pathway enrichment in the cell-line models of aggressive human breast cancer subtypes

This study was conducted as a part of the Chromosome-Centric Human Proteome Project (C-HPP) of the Human Proteome Organization. The United States team of C-HPP is focused on characterizing the protein-coding genes in chromosome 17. Despite its small size, chromosome 17 is rich in protein-coding genes; it contains many cancer-associated genes, including BRCA1, ERBB2, (Her2/neu), and TP53. The goal of this study was to examine the splice variants expressed in three ERBB2 expressed breast cancer cell-line models of hormone-receptor-negative breast cancers by integrating RNA-Seq and proteomic mass spectrometry data. The cell lines represent distinct phenotypic variations subtype: SKBR3 (ERBB2+ (overexpression)/ER-/PR-; adenocarcinoma), SUM190 (ERBB2+ (overexpression)/ER-/PR-; inflammatory breast cancer), and SUM149 (ERBB2 (low expression) ER-/PR-; inflammatory breast cancer). We identified more than one splice variant for 1167 genes expressed in at least one of the three cancer cell lines. We found multiple variants of genes that are in the signaling pathways downstream of ERBB2 along with variants specific to one cancer cell line compared with the other two cancer cell lines and with normal mammary cells. The overall transcript profiles based on read counts indicated more similarities between SKBR3 and SUM190. The top-ranking Gene Ontology and BioCarta pathways for the cell-line specific variants pointed to distinct key mechanisms including: amino sugar metabolism, caspase activity, and endocytosis in SKBR3; different aspects of metabolism, especially of lipids in SUM190; cell-to-cell adhesion, integrin, and ERK1/ERK2 signaling; and translational control in SUM149. The analyses indicated an enrichment in the electron transport chain processes in the ERBB2 overexpressed cell line models and an association of nucleotide binding, RNA splicing, and translation processes with the IBC models, SUM190 and SUM149. Detailed experimental studies on the distinct variants identified from each of these three breast cancer cell line models that may open opportunities for drug target discovery and help unveil their specific roles in cancer progression and metastasis.

Twenty-one proteins up-regulated in human H-ras oncogene transgenic rat pancreas cancers are up-regulated in human pancreas cancer

OBJECTIVES

We have established rat models of pancreatic ductal adenocarcinoma (PDAC) in which expression of a human H-ras(G12V) or K-ras(G12V) oncogene regulated by the Cre/lox system drives pancreatic carcinogenesis. Pancreatic ductal adenocarcinoma which develops in H-ras(G12V) and K-ras(G12V) transgenic rats is cytogenetically and histopathologically similar to human PDAC. The present study was designed to determine the feasibility of using the commercially available H-ras(G12V) transgenic rat to find diagnostic protein biomarkers for human pancreatic cancer.

METHODS

For an animal model to be useful for searching for protein biomarkers for a disease, it is essential that proteins that are up-regulated in the model are also up-regulated in humans. We used liquid chromatography-tandem mass spectrometry (LC-MS/MS) to compare H-ras(G12V) transgenic rat PDAC with surrounding normal pancreas tissue.

RESULTS

We identified 30 up-regulated proteins in the H-ras(G12V) transgenic rat PDAC lesions; importantly, 21 human homologs of these 30 rat proteins are up-regulated in human pancreatic cancer patients.

CONCLUSIONS

These results indicate that numerous proteins that are up-regulated in H-ras(G12V) transgenic rat PDAC are also up-regulated in human pancreatic cancer; therefore, this rat model can be used to search for diagnostic biomarkers for this disease.

Discovery and mass spectrometric analysis of novel splice-junction peptides using RNA-Seq

Human proteomic databases required for MS peptide identification are frequently updated and carefully curated, yet are still incomplete because it has been challenging to acquire every protein sequence from the diverse assemblage of proteoforms expressed in every tissue and cell type. In particular, alternative splicing has been shown to be a major source of this cell-specific proteomic variation. Many new alternative splice forms have been detected at the transcript level using next generation sequencing methods, especially RNA-Seq, but it is not known how many of these transcripts are being translated. Leveraging the unprecedented capabilities of next generation sequencing methods, we collected RNA-Seq and proteomics data from the same cell population (Jurkat cells) and created a bioinformatics pipeline that builds customized databases for the discovery of novel splice-junction peptides. Eighty million paired-end Illumina reads and ∼500,000 tandem mass spectra were used to identify 12,873 transcripts (19,320 including isoforms) and 6810 proteins. We developed a bioinformatics workflow to retrieve high-confidence, novel splice junction sequences from the RNA data, translate these sequences into the analogous polypeptide sequence, and create a customized splice junction database for MS searching. Based on the RefSeq gene models, we detected 136,123 annotated and 144,818 unannotated transcript junctions. Of those, 24,834 unannotated junctions passed various quality filters (e.g. minimum read depth) and these entries were translated into 33,589 polypeptide sequences and used for database searching. We discovered 57 splice junction peptides not present in the Uniprot-Trembl proteomic database comprising an array of different splicing events, including skipped exons, alternative donors and acceptors, and noncanonical transcriptional start sites. To our knowledge this is the first example of using sample-specific RNA-Seq data to create a splice-junction database and discover new peptides resulting from alternative splicing.

Innovations in proteomic profiling of cancers: alternative splice variants as a new class of cancer biomarker candidates and bridging of proteomics with structural biology

Alternative splicing allows a single gene to generate multiple RNA transcripts which can be translated into functionally diverse protein isoforms. Current knowledge of splicing is derived mainly from RNA transcripts, with very little known about the expression level, 3D structures, and functional differences of the proteins. Splicing is a remarkable phenomenon of molecular and biological evolution. Studies which simply report up-regulation or down-regulation of protein or mRNA expression are confounded by the effects of mixtures of these isoforms. Besides understanding the net biological effects of the mixtures, we may be able to develop biomarker tests based on the observable differential expression of particular splice variants or combinations of splice variants in specific disease states. Here we review our work on differential expression of splice variant proteins in cancers and the feasibility of integrating proteomic analysis with structure-based conformational predictions of the differences between such isoforms.

Whole human genome proteogenomic mapping for ENCODE cell line data: identifying protein-coding regions

BACKGROUND

Proteogenomic mapping is an approach that uses mass spectrometry data from proteins to directly map protein-coding genes and could aid in locating translational regions in the human genome. In concert with the ENcyclopedia of DNA Elements (ENCODE) project, we applied proteogenomic mapping to produce proteogenomic tracks for the UCSC Genome Browser, to explore which putative translational regions may be missing from the human genome.

RESULTS

We generated ~1 million high-resolution tandem mass (MS/MS) spectra for Tier 1 ENCODE cell lines K562 and GM12878 and mapped them against the UCSC hg19 human genome, and the GENCODE V7 annotated protein and transcript sets. We then compared the results from the three searches to identify the best-matching peptide for each MS/MS spectrum, thereby increasing the confidence of the putative new protein-coding regions found via the whole genome search. At a 1% false discovery rate, we identified 26,472, 24,406, and 13,128 peptides from the protein, transcript, and whole genome searches, respectively; of these, 481 were found solely via the whole genome search. The proteogenomic mapping data are available on the UCSC Genome Browser at http://genome.ucsc.edu/cgi-bin/hgTrackUi?db=hg19&g=wgEncodeUncBsuProt.

CONCLUSIONS

The whole genome search revealed that ~4% of the uniquely mapping identified peptides were located outside GENCODE V7 annotated exons. The comparison of the results from the disparate searches also identified 15% more spectra than would have been found solely from a protein database search. Therefore, whole genome proteogenomic mapping is a complementary method for genome annotation when performed in conjunction with other searches.

Comparative proteomics reveals a significant bias toward alternative protein isoforms with conserved structure and function

Advances in high-throughput mass spectrometry are making proteomics an increasingly important tool in genome annotation projects. Peptides detected in mass spectrometry experiments can be used to validate gene models and verify the translation of putative coding sequences (CDSs). Here, we have identified peptides that cover 35% of the genes annotated by the GENCODE consortium for the human genome as part of a comprehensive analysis of experimental spectra from two large publicly available mass spectrometry databases. We detected the translation to protein of “novel” and “putative” protein-coding transcripts as well as transcripts annotated as pseudogenes and nonsense-mediated decay targets.

We provide a detailed overview of the population of alternatively spliced protein isoforms that are detectable by peptide identification methods. We found that 150 genes expressed multiple alternative protein isoforms. This constitutes the largest set of reliably confirmed alternatively spliced proteins yet discovered. Three groups of genes were highly overrepresented. We detected alternative isoforms for 10 of the 25 possible heterogeneous nuclear ribonucleoproteins, proteins with a key role in the splicing process. Alternative isoforms generated from interchangeable homologous exons and from short indels were also significantly enriched, both in human experiments and in parallel analyses of mouse and Drosophila proteomics experiments. Our results show that a surprisingly high proportion (almost 25%) of the detected alternative isoforms are only subtly different from their constitutive counterparts.

Many of the alternative splicing events that give rise to these alternative isoforms are conserved in mouse. It was striking that very few of these conserved splicing events broke Pfam functional domains or would damage globular protein structures. This evidence of a strong bias toward subtle differences in CDS and likely conserved cellular function and structure is remarkable and strongly suggests that the translation of alternative transcripts may be subject to selective constraints.

Difference gel electrophoresis identifies differentially expressed proteins in endoscopically collected pancreatic fluid

Alterations in the pancreatic fluid proteome of individuals with chronic pancreatitis (CP) may offer insights into the development and progression of the disease. The endoscopic pancreatic function test (ePFT) can safely collect large volumes of pancreatic fluid that are potentially amenable to proteomic analyses using difference gel electrophoresis (DIGE) coupled with liquid chromatography-tandem mass spectrometry (LC-MS/MS). Pancreatic fluid was collected endoscopically using the ePFT method following secretin stimulation from three individuals with severe CP and three chronic abdominal pain (CAP) controls. The fluid was processed to minimize protein degradation and the protein profiles of each cohort, as determined by DIGE and LC-MS/MS, were compared. This DIGE-LC-MS/MS analysis reveals proteins that are differentially expressed in CP compared with CAP controls. Proteins with higher abundance in pancreatic fluid from CP individuals include: actin, desmoplankin, α-1-antitrypsin, SNC73, and serotransferrin. Those of relatively lower abundance include carboxypeptidase B, lipase, α-1-antichymotrypsin, α-2-macroglobulin, actin-related protein (Arp2/3) subunit 4, glyceraldehyde-3-phosphate dehydrogenase, and protein disulfide isomerase. Endoscopic collection (ePFT) in tandem with DIGE-LC-MS/MS is a suitable approach for pancreatic fluid proteome analysis; however, further optimization of our protocol, as outlined herein, may improve proteome coverage in future analyses.

Evaluating the potential of a novel oral lesion exudate collection method coupled with mass spectrometry-based proteomics for oral cancer biomarker discovery

Introduction

Early diagnosis of Oral Squamous Cell Carcinoma (OSCC) increases the survival rate of oral cancer. For early diagnosis, molecular biomarkers contained in samples collected non-invasively and directly from at-risk oral premalignant lesions (OPMLs) would be ideal.

Methods

In this pilot study we evaluated the potential of a novel method using commercial PerioPaper absorbent strips for non-invasive collection of oral lesion exudate material coupled with mass spectrometry-based proteomics for oral cancer biomarker discovery.

Results

Our evaluation focused on three core issues. First, using an "on-strip" processing method, we found that protein can be isolated from exudate samples in amounts compatible with large-scale mass spectrometry-based proteomic analysis. Second, we found that the OPML exudate proteome was distinct from that of whole saliva, while being similar to the OPML epithelial cell proteome, demonstrating the fidelity of our exudate collection method. Third, in a proof-of-principle study, we identified numerous, inflammation-associated proteins showing an expected increase in abundance in OPML exudates compared to healthy oral tissue exudates. These results demonstrate the feasibility of identifying differentially abundant proteins from exudate samples, which is essential for biomarker discovery studies.

Conclusions

Collectively, our findings demonstrate that our exudate collection method coupled with mass spectrometry-based proteomics has great potential for transforming OSCC biomarker discovery and clinical diagnostics assay development.

A proteogenomic analysis of Anopheles gambiae using high-resolution Fourier transform mass spectrometry

Anopheles gambiae is a major mosquito vector responsible for malaria transmission, whose genome sequence was reported in 2002. Genome annotation is a continuing effort, and many of the approximately 13,000 genes listed in VectorBase for Anopheles gambiae are predictions that have still not been validated by any other method. To identify protein-coding genes of An. gambiae based on its genomic sequence, we carried out a deep proteomic analysis using high-resolution Fourier transform mass spectrometry for both precursor and fragment ions. Based on peptide evidence, we were able to support or correct more than 6000 gene annotations including 80 novel gene structures and about 500 translational start sites. An additional validation by RT-PCR and cDNA sequencing was successfully performed for 105 selected genes. Our proteogenomic analysis led to the identification of 2682 genome search-specific peptides. Numerous cases of encoded proteins were documented in regions annotated as intergenic, introns, or untranslated regions. Using a database created to contain potential splice sites, we also identified 35 novel splice junctions. This is a first report to annotate the An. gambiae genome using high-accuracy mass spectrometry data as a complementary technology for genome annotation.