Massive Proteogenomic Reanalysis of Publicly Available Proteomic Datasets of Human Tissues in Search for Protein Recoding via Adenosine-to-Inosine RNA Editing.

The proteogenomic search pipeline developed in this work has been applied for reanalysis of 40 publicly available shotgun proteomic datasets from various human tissues comprising more than 8000 individual LC-MS/MS runs, of which 5442 .raw data files were processed in total. This reanalysis was focused on searching for ADAR-mediated RNA editing events, their clustering across samples of different origins, and classification. In total, 33 recoded protein sites were identified in 21 datasets. Of those, 18 sites were detected in at least two datasets, representing the core human protein editome. In agreement with prior artworks, neural and cancer tissues were found to be enriched with recoded proteins. Quantitative analysis indicated that recoding the rate of specific sites did not directly depend on the levels of ADAR enzymes or targeted proteins themselves, rather it was governed by differential and yet undescribed regulation of interaction of enzymes with mRNA. Nine recoding sites conservative between humans and rodents were validated by targeted proteomics using stable isotope standards in the murine brain cortex and cerebellum, and an additional one was validated in human cerebrospinal fluid. In addition to previous data of the same type from cancer proteomes, we provide a comprehensive catalog of recoding events caused by ADAR RNA editing in the human proteome.

Identification of Alternative Splicing in Proteomes of Human Melanoma Cell Lines without RNA Sequencing Data.

Alternative splicing is one of the main regulation pathways in living cells beyond simple changes in the level of protein expression. Most of the approaches proposed in proteomics for the identification of specific splicing isoforms require a preliminary deep transcriptomic analysis of the sample under study, which is not always available, especially in the case of the re-analysis of previously acquired data. Herein, we developed new algorithms for the identification and validation of protein splice isoforms in proteomic data in the absence of RNA sequencing of the samples under study. The bioinformatic approaches were tested on the results of proteome analysis of human melanoma cell lines, obtained earlier by high-resolution liquid chromatography and mass spectrometry (LC-MS). A search for alternative splicing events for each of the cell lines studied was performed against the database generated from all known transcripts (RefSeq) and the one composed of peptide sequences, which included all biologically possible combinations of exons. The identifications were filtered using the prediction of both retention times and relative intensities of fragment ions in the corresponding mass spectra. The fragmentation mass spectra corresponding to the discovered alternative splicing events were additionally examined for artifacts. Selected splicing events were further validated at the mRNA level by quantitative PCR.

Electrochemical biosensor for trypsin activity assay based on cleavage of immobilized tyrosine-containing peptide.

In this work, we proposed a biosensor for trypsin proteolytic activity assay using immobilization of model peptides on screen-printed electrodes (SPE) modified with gold nanoparticles (AuNPs) prepared by electrosynthetic method. Sensing of proteolytic activity was based on electrochemical oxidation of tyrosine residues of peptides. We designed peptides containing N-terminal cysteine residue for immobilization on an SPE, modified with gold nanoparticles, trypsin-specific cleavage site and tyrosine residue as a redox label. The peptides were immobilized on SPE by formation of chemical bonds between mercapto groups of the N-terminal cysteine residues and AuNPs. After the incubation with trypsin, time-dependent cleavage of the immobilized peptides was observed by decline in tyrosine electrochemical oxidation signal. The kinetic parameters of trypsin, such as the catalytic constant (kcat), the Michaelis constant (KM) and the catalytic efficiency (kcat/KM), toward the CGGGRYR peptide were determined as 0.33 ± 0.01 min-1, 198 ± 24 nM and 0.0016 min-1 nM-1, respectively. Using the developed biosensor, we demonstrated the possibility of analysis of trypsin specificity toward the peptides with amino acid residues disrupting proteolysis. Further, we designed the peptides with proline or glutamic acid residues after the cleavage site (CGGRPYR and CGGREYR), and trypsin had reduced activity toward both of them according to the existing knowledge of the enzyme specificity. The developed biosensor system allows one to perform a comparative analysis of the protease steady-state kinetic parameters and specificity toward model peptides with different amino acid sequences.

Blood Plasma Proteome: A Meta-Analysis of the Results of Protein Quantification in Human Blood by Targeted Mass Spectrometry.

A meta-analysis of the results of targeted quantitative screening of human blood plasma was performed to generate a reference standard kit that can be used for health analytics. The panel included 53 of the 296 proteins that form a "stable" part of the proteome of a healthy individual; these proteins were found in at least 70% of samples and were characterized by an interindividual coefficient of variation <40%. The concentration range of the selected proteins was 10−10−10−3 M and enrichment analysis revealed their association with rare familial diseases. The concentration of ceruloplasmin was reduced by approximately three orders of magnitude in patients with neurological disorders compared to healthy volunteers, and those of gelsolin isoform 1 and complement factor H were abruptly reduced in patients with lung adenocarcinoma. Absolute quantitative data of the individual proteome of a healthy and diseased individual can be used as the basis for personalized medicine and health monitoring. Storage over time allows us to identify individual biomarkers in the molecular landscape and prevent pathological conditions.

Proteomic Analysis of Zebrafish Protein Recoding via mRNA Editing by ADAR Enzymes.

RNA editing by adenosine deaminases of the ADAR family can lead to protein recoding, since inosine formed from adenosine in mRNA is complementary to cytosine; the resulting codon editing might introduce amino acid substitutions into translated proteins. Proteome recoding can have functional consequences which have been described in many animals including humans. Using protein recoding database derived from publicly available transcriptome data, we identified for the first time the recoding sites in the zebrafish shotgun proteomes. Out of more than a hundred predicted recoding events, ten substitutions were found in six used datasets. Seven of them were in the AMPA glutamate receptor subunits, whose recoding has been well described, and are conserved among vertebrates. Three sites were specific for zebrafish proteins and were found in the transmembrane receptors astrotactin 1 and neuregulin 3b (proteins involved in the neuronal adhesion and signaling) and in the rims2b gene product (presynaptic membrane protein participating in the neurotransmitter release), respectively. Further studies are needed to elucidate the role of recoding of the said three proteins in the zebrafish.

Vexitoxins: conotoxin-like venom peptides from predatory gastropods of the genus Vexillum.

Venoms of predatory marine cone snails are intensely studied because of the biomedical applications of the neuropeptides that they contain, termed conotoxins. Meanwhile some gastropod lineages have independently acquired secretory glands strikingly similar to the venom gland of cone snails, suggesting that they possess similar venoms. Here we focus on the most diversified of these clades, the genus Vexillum. Based on the analysis of a multi-species proteo-transcriptomic dataset, we show that Vexillum species indeed produce complex venoms dominated by highly diversified short cysteine-rich peptides, vexitoxins. Vexitoxins possess the same precursor organization, display overlapping cysteine frameworks and share several common post-translational modifications with conotoxins. Some vexitoxins show sequence similarity to conotoxins and adopt similar domain conformations, including a pharmacologically relevant inhibitory cysteine knot motif. The Vexillum envenomation gland (gL) is a notably more recent evolutionary novelty than the conoidean venom gland. Thus, we hypothesize lower divergence between vexitoxin genes, and their ancestral 'somatic' counterparts compared to that in conotoxins, and we find support for this hypothesis in the evolution of the vexitoxin cluster V027. We use this example to discuss how future studies on vexitoxins can inform the origin of conotoxins, and how they may help to address outstanding questions in venom evolution.

Interplay between A-to-I Editing and Splicing of RNA: A Potential Point of Application for Cancer Therapy.

Adenosine-to-inosine RNA editing is a system of post-transcriptional modification widely distributed in metazoans which is catalyzed by ADAR enzymes and occurs mostly in double-stranded RNA (dsRNA) before splicing. This type of RNA editing changes the genetic code, as inosine generally pairs with cytosine in contrast to adenosine, and this expectably modulates RNA splicing. We review the interconnections between RNA editing and splicing in the context of human cancer. The editing of transcripts may have various effects on splicing, and resultant alternatively spliced isoforms may be either tumor-suppressive or oncogenic. Dysregulated RNA splicing in cancer often causes the release of excess amounts of dsRNA into cytosol, where specific dsRNA sensors provoke antiviral-like responses, including type I interferon signaling. These responses may arrest cell division, causing apoptosis and, externally, stimulate antitumor immunity. Thus, small-molecule spliceosome inhibitors have been shown to facilitate the antiviral-like signaling and are considered to be potential cancer therapies. In turn, a cytoplasmic isoform of ADAR can deaminate dsRNA in cytosol, thereby decreasing its levels and diminishing antitumor innate immunity. We propose that complete or partial inhibition of ADAR may enhance the proapoptotic and cytotoxic effects of splicing inhibitors and that it may be considered a promising addition to cancer therapies targeting RNA splicing.

Multiomic Profiling Identified EGF Receptor Signaling as a Potential Inhibitor of Type I Interferon Response in Models of Oncolytic Therapy by Vesicular Stomatitis Virus.

Cancer cell lines responded differentially to type I interferon treatment in models of oncolytic therapy using vesicular stomatitis virus (VSV). Two opposite cases were considered in this study, glioblastoma DBTRG-05MG and osteosarcoma HOS cell lines exhibiting resistance and sensitivity to VSV after the treatment, respectively. Type I interferon responses were compared for these cell lines by integrative analysis of the transcriptome, proteome, and RNA editome to identify molecular factors determining differential effects observed. Adenosine-to-inosine RNA editing was equally induced in both cell lines. However, transcriptome analysis showed that the number of differentially expressed genes was much higher in DBTRG-05MG with a specific enrichment in inflammatory proteins. Further, it was found that two genes, EGFR and HER2, were overexpressed in HOS cells compared with DBTRG-05MG, supporting recent reports that EGF receptor signaling attenuates interferon responses via HER2 co-receptor activity. Accordingly, combined treatment of cells with EGF receptor inhibitors such as gefitinib and type I interferon increases the resistance of sensitive cell lines to VSV. Moreover, sensitive cell lines had increased levels of HER2 protein compared with non-sensitive DBTRG-05MG. Presumably, the level of this protein expression in tumor cells might be a predictive biomarker of their resistance to oncolytic viral therapy.

Validating Amino Acid Variants in Proteogenomics Using Sequence Coverage by Multiple Reads.

Mass spectrometry-based proteome analysis implies matching the mass spectra of proteolytic peptides to amino acid sequences predicted from genomic sequences. Reliability of peptide variant identification in proteogenomic studies is often lacking. We propose a way to interpret shotgun proteomics results, specifically in the data-dependent acquisition mode, as protein sequence coverage by multiple reads as it is done in nucleic acid sequencing for calling of single nucleotide variants. Multiple reads for each sequence position could be provided by overlapping distinct peptides, thus confirming the presence of certain amino acid residues in the overlapping stretch with a lower false discovery rate. Overlapping distinct peptides originate from miscleaved tryptic peptides in combination with their properly cleaved counterparts and from peptides generated by multiple proteases after the same specimen is subject to parallel digestion and analyzed separately. We illustrate this approach using publicly available multiprotease data sets and our own data generated for the HEK-293 cell line digests obtained using trypsin, LysC, and GluC proteases. Totally, up to 30% of the whole proteome was covered by tryptic peptides with up to 7% covered twofold and more. The proteogenomic analysis of the HEK-293 cell line revealed 36 single amino acid variants, seven of which were supported by multiple reads.

Multi-Omics Analysis of Glioblastoma Cells' Sensitivity to Oncolytic Viruses.

Oncolytic viruses have gained momentum in the last decades as a promising tool for cancer treatment. Despite the progress, only a fraction of patients show a positive response to viral therapy. One of the key variable factors contributing to therapy outcomes is interferon-dependent antiviral mechanisms in tumor cells. Here, we evaluated this factor using patient-derived glioblastoma multiforme (GBM) cultures. Cell response to the type I interferons' (IFNs) stimulation was characterized at mRNA and protein levels. Omics analysis revealed that GBM cells overexpress interferon-stimulated genes (ISGs) and upregulate their proteins, similar to the normal cells. A conserved molecular pattern unambiguously differentiates between the preserved and defective responses. Comparing ISGs' portraits with titration-based measurements of cell sensitivity to a panel of viruses, the "strength" of IFN-induced resistance acquired by GBM cells was ranked. The study demonstrates that suppressing a single ISG and encoding an essential antiviral protein, does not necessarily increase sensitivity to viruses. Conversely, silencing IFIT3 and PLSCR1 genes in tumor cells can negatively affect the internalization of vesicular stomatitis and Newcastle disease viruses. We present evidence of a complex relationship between the interferon response genes and other factors affecting the sensitivity of tumor cells to viruses.

Cysteine alkylation methods in shotgun proteomics and their possible effects on methionine residues.

In order to optimize sample preparation for shotgun proteomics, we compared four cysteine alkylating agents: iodoacetamide, chloroacetamide, 4-vinylpyridine and methyl methanethiosulfonate, and estimated their effects on the results of proteome analysis. Because alkylation may result in methionine modification in vitro, proteomics data were searched for methionine to isothreonine conversions, which may mimic genomic methionine to threonine substitutions found in proteogenomic analyses. We found that chloroacetamide was superior to the other reagents in terms of the number of identified peptides and undesirable off-site reactions. Among the reagents evaluated, iodoacetamide increased the rate of methionine-to-isothreonine conversion, especially if the sample was prepared in gel. The presence of proline following methionine in a protein sequence increased the modification rate as well. Generally, the methionine-to-isothreonine conversion events were relatively rare, but should be taken into account in proteogenomic studies when searching for single nucleotide polymorphism events at the protein level. Additionally, we have evaluated other methionine modifications, such as oxidation and carbamidomethylation. We found that carbamidomethylation may affect up to 80% of peptides containing methionine under the condition of iodoacetamide alkylation. In this case, carbamidomethylation of methionine is more common than oxidation and should be accounted for as a variable modification during proteomic search. SIGNIFICANCE: One of the most trending questions in bottom-up proteomics is the depth of proteome profiling, in other words, the coverage of proteins by identified tryptic peptides. In proteogenomics, where the identification of a single peptide, e.g. bearing an amino acid substitution, may be of interest, high sequence coverage is especially important. Chemical modifications during sample preparation may mimic biologically significant coding mutations at the proteome level. A typical example of such modification is methionine to isothreonine conversion during alkylation, which mimics methionine to threonine substitution in protein sequences due to respective genomic mutations. Therefore, the studies on the proper selection of alkylating reagents which balance the cysteine alkylation efficiency and the extent of methionine conversion upon conventional proteomic sample preparation workflow are crucial for the outcome of proteogenomic analyses and should present a general interest for the proteomic community.

Identification-Free Control over the Precursor Isotopic Mass Misassignment in Orbitrap-Based Proteomics.

Selection of a precursor ion from a peptide isotopic cluster to obtain a fragmentation mass spectrum is a crucial step in data-dependent proteome analysis. However, the monoisotopic mass assignment performed in this step is often an issue confronted by the data acquisition software of hybrid Orbitrap FTMS that is most widely used in proteomics. To address the problem, many data processing tools, such as raw data converters and search engines, have optional accounting for the precursor mass shift due to the isotopic error. These solutions require additional data preprocessing steps and lead to an increase in the search space, thus making the analysis longer and/or less reliable. In this work, we processed 100 Orbitrap-based LC-MS/MS runs from 10 publicly available data sets to examine the rate of precursor isotope misassignment. The effect from taking the isotope error into account during the search on the number of identified peptides varied in a wide range from 0 to 33%. Thus, it may be tempting to spend extra time before or during a search to account for the mass assignment issue. Alternatively, this effect can be predicted a priori using an identification-free metric, which can be a part of data quality control software. Based on the results obtained in this work, we propose such a metric be further added into the visual and intuitive quality control software, viQC, developed previously and available at https://github.com/lisavetasol/viQC. It takes about a minute to calculate and plot nine quality metrics, including the proposed one for typical proteome analysis.

Analysis of candidate genes expected to be essential for melanoma surviving.

INTRODUCTION: Cancers may be treated by selective targeting of the genes vital for their survival. A number of attempts have led to discovery of several genes essential for surviving of tumor cells of different types. In this work, we tried to analyze genes that were previously predicted to be essential for melanoma surviving. Here we present the results of transient siRNA-mediated knockdown of the four of such genes, namely, UNC45A, STK11IP, RHPN2 and ZNFX1, in melanoma cell line A375, then assayed the cells for their viability, proliferation and ability to migrate in vitro. In our study, the knockdown of the genes predicted as essential for melanoma survival does not lead to statistically significant changes in cell viability. On the other hand, for each of the studied genes, mobility assays showed that the knockdown of each of the target genes accelerates the speed of cells migrating. Possible explanation for such counterintuitive results may include insufficiency of the predicting computational models or the necessity of a multiplex knockdown of the genes. AIMS: To examine the hypothesis of essentiality of hypomutated genes for melanoma surviving we have performed knockdown of several genes in melanoma cell line and analyzed cell viability and their ability to migrate. METHODS: Knockdown was performed by siRNAs transfected by Metafectene PRO. The levels of mRNAs before and after knockdown were evaluated by RT-qPCR analysis. Cell viability and proliferation were assessed by MTT assay. Cell migration was assessed by wound healing assay. RESULTS: The knockdown of the genes predicted as essential for melanoma survival does not lead to statistically significant changes in cell viability. On the other hand, for each of the studied genes, mobility assays showed that the knockdown of each of the target genes accelerates the speed of cells migrating. CONCLUSION: Our results do not confirm initial hypothesis that the genes predicted essential for melanoma survival as a matter of fact support the survival of melanoma cells.

Proteome-Wide Analysis of ADAR-Mediated Messenger RNA Editing during Fruit Fly Ontogeny.

Adenosine-to-inosine RNA editing is an enzymatic post-transcriptional modification which modulates immunity and neural transmission in multicellular organisms. In particular, it involves editing of mRNA codons with the resulting amino acid substitutions. We identified such sites for developmental proteomes of Drosophila melanogaster at the protein level using available data for 15 stages of fruit fly development from egg to imago and 14 time points of embryogenesis. In total, 40 sites were obtained, each belonging to a unique protein, including four sites related to embryogenesis. The interactome analysis has revealed that the majority of the editing-recoded proteins were associated with synaptic vesicle trafficking and actomyosin organization. Quantitation data analysis suggested the existence of a phase-specific RNA editing regulation with yet unknown mechanisms. These findings supported the transcriptome analysis results, which showed that a burst in the RNA editing occurs during insect metamorphosis from pupa to imago. Finally, targeted proteomic analysis was performed to quantify editing-recoded and genomically encoded versions of five proteins in brains of larvae, pupae, and imago insects, which showed a clear tendency toward an increase in the editing rate for each of them. These results will allow a better understanding of the protein role in physiological effects of RNA editing.

Adenosine-to-Inosine RNA Editing in Mouse and Human Brain Proteomes.

Proteogenomics is based on the use of customized genome or RNA sequencing databases for interrogation of shotgun proteomics data in search for proteome-level evidence of genome variations or RNA editing. In this work, the products of adenosine-to-inosine RNA editing in human and murine brain proteomes are identified using publicly available brain proteome LC-MS/MS datasets and an RNA editome database compiled from several sources. After filtering of false-positive results, 20 and 37 sites of editing in proteins belonging to 14 and 32 genes are identified for murine and human brain proteomes, respectively. Eight sites of editing identified with high spectral counts overlapped between human and mouse brain samples. Some of these sites have been previously reported using orthogonal methods, such as α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate receptors, CYFIP2, coatomer alpha. Also, differential editing between neurons and microglia is demonstrated in this work for some of the proteins from primary murine brain cell cultures. Because many edited sites are still not characterized functionally at the protein level, the results provide a necessary background for their further analysis in normal and diseased cells and tissues using targeted proteomic approaches.

Method for Identification of Threonine Isoforms in Peptides by Ultraviolet Photofragmentation of Cold Ions.

Identification of isomeric amino acid residues in peptides and proteins is challenging but often highly desired in proteomics. One of the practically important cases that require isomeric assignments is that associated with single-nucleotide polymorphism substitutions of Met residues by Thr in cancer-related proteins. These genetically encoded substitutions can yet be confused with the chemical modifications, arising from protein alkylation by iodoacetamide, which is commonly used in the standard procedure of sample preparation for proteomic analysis. Similar to the genetically encoded mutations, the alkylation also induces a conversion of methionine residues, but to the iso-threonine form. Recognition of the mutations therefore requires isoform-sensitive detection techniques. Herein, we demonstrate an analytical method for reliable identification of isoforms of threonine residues in tryptic peptides. It is based on ultraviolet photodissociation mass spectrometry of cryogenically cooled ions and a machine-learning algorithm. The measured photodissociation mass spectra exhibit isoform-specific patterns, which are independent of the residues adjacent to threonine or iso-threonine in a peptide sequence. A comprehensive metric-based evaluation demonstrates that, being calibrated with a set of model peptides, the method allows for isomeric identification of threonine residues in peptides of arbitrary sequence.

Profiling modifications for glioblastoma proteome using ultra-tolerant database search: Are the peptide mass shifts biologically relevant or chemically induced?

Peptide mass shifts were profiled using ultra-tolerant database search strategy for shotgun proteomics data sets of human glioblastoma cell lines demonstrating strong response to the type I interferon (IFNα-2b) treatment. The main objective of this profiling was revealing the cell response to IFN treatment at the level of protein modifications. To achieve this objective, statistically significant changes in peptide mass shift profiles between IFN treated and untreated glioblastoma samples were analyzed. Detailed analysis of MS/MS spectra allowed further interpretation of the observed mass shifts and differentiation between post-translational and artifact modifications. SIGNIFICANCE: Malignant cells typically acquire increased sensitivity to viruses due to the deregulated antiviral mechanisms. Therefore, a viral therapy is considered as one of the promising approaches to treat cancer. However, recent studies have demonstrated that malignant cells can preserve intact antiviral mechanisms, e.g. interferon signaling, and develop resistance to virus infection in response to interferon treatment. Post translational modifications, e.g. tyrosine phosphorylation, are the interferon signaling drivers. Thus, comprehensive characterization of modifications is crucially important, yet, most challenging problem in cancer proteomics. Here, we report on the application of the recently introduced ultra-tolerant search strategy for profiling peptide modifications in the human glioblastoma cell lines demonstrating strong response to the type I interferon (IFNα-2b) treatment. The specific aim of the study was identification of statistically significant changes in peptide mass shift profiles between IFN treated and untreated glioblastoma samples, as well as determination of whether these shifts represent the biologically relevant modification.

Proteogenomics of Adenosine-to-Inosine RNA Editing in the Fruit Fly.

Adenosine-to-inosine RNA editing is one of the most common types of RNA editing, a posttranscriptional modification made by special enzymes. We present a proteomic study on this phenomenon for Drosophila melanogaster. Three proteome data sets were used in the study: two taken from public repository and the third one obtained here. A customized protein sequence database was generated using results of genome-wide adenosine-to-inosine RNA studies and applied for identifying the edited proteins. The total number of 68 edited peptides belonging to 59 proteins was identified in all data sets. Eight of them being shared between the whole insect, head, and brain proteomes. Seven edited sites belonging to synaptic vesicle and membrane trafficking proteins were selected for validation by orthogonal analysis by Multiple Reaction Monitoring. Five editing events in cpx, Syx1A, Cadps, CG4587, and EndoA were validated in fruit fly brain tissue at the proteome level using isotopically labeled standards. Ratios of unedited-to-edited proteoforms varied from 35:1 ( Syx1A) to 1:2 ( EndoA). Lys-137 to Glu editing of endophilin A may have functional consequences for its interaction to membrane. The work demonstrates the feasibility to identify the RNA editing event at the proteome level using shotgun proteomics and customized edited protein database.

FractionOptimizer: a method for optimal peptide fractionation in bottom-up proteomics.

Recent advances in mass spectrometry and separation technologies created the opportunities for deep proteome characterization using shotgun proteomics approaches. The "real world" sample complexity and high concentration range limit the sensitivity of this characterization. The common strategy for increasing the sensitivity is sample fractionation prior to analysis either at the protein or the peptide level. Typically, fractionation at the peptide level is performed using linear gradient high-performance liquid chromatography followed by uniform fraction collection. However, this way of peptide fractionation results in significantly suboptimal operation of the mass spectrometer due to the non-uniform distribution of peptides between the fractions. In this work, we propose an approach based on peptide retention time prediction allowing optimization of chromatographic conditions and fraction collection procedures. An open-source software implementing the approach called FractionOptimizer was developed and is available at http://hg.theorchromo.ru/FractionOptimizer . The performance of the developed tool was demonstrated for human embryonic kidney (HEK293) cell line lysate. In these experiments, we improved the uniformity of the peptides distribution between fractions. Moreover, in addition to 13,492 peptides, we found 6787 new peptides not identified in the experiments without fractionation and up to 800 new proteins (or 25%). Graphical abstract The analysis workflow employing FractionOptimizer software.

Proteogenomics of Malignant Melanoma Cell Lines: The Effect of Stringency of Exome Data Filtering on Variant Peptide Identification in Shotgun Proteomics.

The identification of genetically encoded variants at the proteome level is an important problem in cancer proteogenomics. The generation of customized protein databases from DNA or RNA sequencing data is a crucial stage of the identification workflow. Genomic data filtering applied at this stage may significantly modify variant search results, yet its effect is generally left out of the scope of proteogenomic studies. In this work, we focused on this impact using data of exome sequencing and LC-MS/MS analyses of six replicates for eight melanoma cell lines processed by a proteogenomics workflow. The main objectives were identifying variant peptides and revealing the role of the genomic data filtering in the variant identification. A series of six confidence thresholds for single nucleotide polymorphisms and indels from the exome data were applied to generate customized sequence databases of different stringency. In the searches against unfiltered databases, between 100 and 160 variant peptides were identified for each of the cell lines using X!Tandem and MS-GF+ search engines. The recovery rate for variant peptides was ∼1%, which is approximately three times lower than that of the wild-type peptides. Using unfiltered genomic databases for variant searches resulted in higher sensitivity and selectivity of the proteogenomic workflow and positively affected the ability to distinguish the cell lines based on variant peptide signatures.