Exome-based proteogenomics of HEK-293 human cell line: Coding genomic variants identified at the level of shotgun proteome

Proteogenomics in Cancer: Then and Now

For years, the paths of sequencing technologies and mass spectrometry have occurred in isolation, with each developing its own unique culture and expertise. These two technologies are crucial for inspecting complementary aspects of the molecular phenotype across the central dogma. Integrative multiomics strives to bridge the analysis gap among different fields to complete more comprehensive mechanisms of life events and diseases. Proteogenomics is one integrated multiomics field. Here in this review, we mainly summarize and discuss three aspects: workflow of proteogenomics, proteogenomics applications in cancer research, and the SWOT (Strengths, Weaknesses, Opportunities, Threats) analysis of proteogenomics in cancer research. In conclusion, proteogenomics has a promising future as it clarifies the functional consequences of many unannotated genomic abnormalities or noncanonical variants and identifies driver genes and novel therapeutic targets across cancers, which would substantially accelerate the development of precision oncology.

Finding haplotypic signatures in proteins

BACKGROUND

The nonrandom distribution of alleles of common genomic variants produces haplotypes, which are fundamental in medical and population genetic studies. Consequently, protein-coding genes with different co-occurring sets of alleles can encode different amino acid sequences: protein haplotypes. These protein haplotypes are present in biological samples and detectable by mass spectrometry, but they are not accounted for in proteomic searches. Consequently, the impact of haplotypic variation on the results of proteomic searches and the discoverability of peptides specific to haplotypes remain unknown.

FINDINGS

Here, we study how common genetic haplotypes influence the proteomic search space and investigate the possibility to match peptides containing multiple amino acid substitutions to a publicly available data set of mass spectra. We found that for 12.42% of the discoverable amino acid substitutions encoded by common haplotypes, 2 or more substitutions may co-occur in the same peptide after tryptic digestion of the protein haplotypes. We identified 352 spectra that matched to such multivariant peptides, and out of the 4,582 amino acid substitutions identified, 6.37% were covered by multivariant peptides. However, the evaluation of the reliability of these matches remains challenging, suggesting that refined error rate estimation procedures are needed for such complex proteomic searches.

CONCLUSIONS

As these procedures become available and the ability to analyze protein haplotypes increases, we anticipate that proteomics will provide new information on the consequences of common variation, across tissues and time.

Multiaspect Examinations of Possible Alternative Mappings of Identified Variant Peptides: A Case Study on the HEK293 Cell Line

Adopting proteogenomics approach to validate single nucleotide variation events by identifying corresponding single amino acid variant peptides from mass spectrometry (MS)-based proteomics data facilitates translational and clinical research. Although variant peptides are usually identified from MS data with a stringent false discovery rate (FDR), FDR control could fail to eliminate dubious results caused by several issues; thus, postexamination to eliminate dubious results is required. However, comprehensive postexaminations of identification results are still lacking. Therefore, we propose a framework of three bottom-up levels, peptide-spectrum match, peptide, and variant event levels, that consists of rigorous 11-aspect examinations from the MS perspective to further confirm the reliability of variant events. As a proof of concept and showing feasibility, we demonstrate 11 examinations on the identified variant peptides from an HEK293 cell line data set, where various database search strategies were applied to maximize the number of identified variant PSMs with an FDR <1% for postexaminations. The results showed that only FDR criterion is insufficient to validate identified variant peptides and the 11 postexaminations can reveal low-confidence variant events detected by shotgun proteomics experiments. Therefore, we suggest that postexaminations of identified variant events based on the proposed framework are necessary for proteogenomics studies.

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.

False discovery rate: the Achilles' heel of proteogenomics

Proteogenomics refers to the integrated analysis of the genome and proteome that leverages mass-spectrometry (MS)-based proteomics data to improve genome annotations, understand gene expression control through proteoforms and find sequence variants to develop novel insights for disease classification and therapeutic strategies. However, proteogenomic studies often suffer from reduced sensitivity and specificity due to inflated database size. To control the error rates, proteogenomics depends on the target-decoy search strategy, the de-facto method for false discovery rate (FDR) estimation in proteomics. The proteogenomic databases constructed from three- or six-frame nucleotide database translation not only increase the search space and compute-time but also violate the equivalence of target and decoy databases. These searches result in poorer separation between target and decoy scores, leading to stringent FDR thresholds. Understanding these factors and applying modified strategies such as two-pass database search or peptide-class-specific FDR can result in a better interpretation of MS data without introducing additional statistical biases. Based on these considerations, a user can interpret the proteogenomics results appropriately and control false positives and negatives in a more informed manner. In this review, first, we briefly discuss the proteogenomic workflows and limitations in database construction, followed by various considerations that can influence potential novel discoveries in a proteogenomic study. We conclude with suggestions to counter these challenges for better proteogenomic data interpretation.

Proteomic and electron microscopy study of myogenic differentiation of alveolar mucosa multipotent mesenchymal stromal cells in three-dimensional culture

Myocyte differentiation is featured by adaptation processes, including mitochondria repopulation and cytoskeleton re-organization. The difference between monolayer and spheroid cultured cells at the proteomic level is uncertain. We cultivated alveolar mucosa multipotent mesenchymal stromal cells in spheroids in a myogenic way for the proper conditioning of ECM architecture and cell morphology, which induced spontaneous myogenic differentiation of cells within spheroids. Electron microscopy analysis was used for the morphometry of mitochondria biogenesis, and proteomic was used complementary to unveil events underlying differences between two-dimensional/three-dimensional myoblasts differentiation. The prevalence of elongated mitochondria with an average area of 0.097 μm² was attributed to monolayer cells 7 days after the passage. The population of small mitochondria with a round shape and area of 0.049 μm² (p < 0.05) was observed in spheroid cells cultured under three-dimensional conditions. Cells in spheroids were quantitatively enriched in proteins of mitochondria biogenesis (DNM1L, IDH2, SSBP1), respiratory chain (ACO2, ATP5I, COX5A), extracellular proteins (COL12A1, COL6A1, COL6A2), and cytoskeleton (MYL6, MYL12B, MYH10). Most of the Rab-related transducers were inhibited in spheroid culture. The proteomic assay demonstrated delicate mechanisms of mitochondria autophagy and repopulation, cytoskeleton assembling, and biogenesis. Differences in the ultrastructure of mitochondria indicate active biogenesis under three-dimensional conditions.

Proteomes Are of Proteoforms: Embracing the Complexity

Proteomes are complex-much more so than genomes or transcriptomes. Thus, simplifying their analysis does not simplify the issue. Proteomes are of proteoforms, not canonical proteins. While having a catalogue of amino acid sequences provides invaluable information, this is the Proteome-lite. To dissect biological mechanisms and identify critical biomarkers/drug targets, we must assess the myriad of proteoforms that arise at any point before, after, and between translation and transcription (e.g., isoforms, splice variants, and post-translational modifications [PTM]), as well as newly defined species. There are numerous analytical methods currently used to address proteome depth and here we critically evaluate these in terms of the current 'state-of-the-field'. We thus discuss both pros and cons of available approaches and where improvements or refinements are needed to quantitatively characterize proteomes. To enable a next-generation approach, we suggest that advances lie in transdisciplinarity via integration of current proteomic methods to yield a unified discipline that capitalizes on the strongest qualities of each. Such a necessary (if not revolutionary) shift cannot be accomplished by a continued primary focus on proteo-genomics/-transcriptomics. We must embrace the complexity. Yes, these are the hard questions, and this will not be easy…but where is the fun in easy?

Cytotoxic and Apoptotic Effects of Celecoxib and Topotecan on AGS and HEK 293 Cell Lines

PURPOSE

This study is aimed to assess the anti-cancer effects of Celecoxib and topotecan against Human Gastric cancer cell line (AGS) in comparison to the control in an in-vitro study.

METHODS

In this experimental study, Celecoxib and topotecan was prepared at concentrations of 500, 250, 125, 62.5, 31.2, 15.6 and 7.8 mg/ml. The effect of celecoxib and topotecan separately and in mixed form were investigated on AGS and normal HEK cells. To investigate the cell survival, MTT method was used to study the pathway of apoptosis using flowcytometry and Caspase kits based on colorimetric. Finally, one-way ANOVA and t-test were used to analyze the data.

RESULTS

The results of this study indicated that Celecoxib was cytotoxic against AGS and HEK cell lines. The topotecan indicated a significant cytotoxicity against AGS cells and was not toxic against HEK cell line. Our results indicated that Celecoxib and topotecan have synergist effects against AGS and HEK cell lines and were more effective than separate celecoxib or topotecan.

CONCLUSION

The mixture of clecoxib and topotecan was more effective than celecoxib and topotecan in separate form. Our results indicated that use mixed forms of treatments can cause excellent therapeutic effects and can cause less side effects.

MinProtMaxVP: Generating a minimized number of protein variant sequences containing all possible variant peptides for proteogenomic analysis

Identifying single-amino-acid variants (SAVs) from mass spectrometry-based experiments is critical for validating single-nucleotide variants (SNVs) at the protein level to facilitate biomedical research. Currently, two approaches are usually applied to convert SNV annotations into SAV-harboring protein sequences. One approach generates one sequence containing exactly one SAV, and the other all SAVs. However, they may neglect the possibility of SAV combinations, e.g., haplotypes, existing in bio-samples. Therefore, it is necessary to consider all SAV combinations of a protein when generating SAV-harboring protein sequences. In this paper, we propose MinProtMaxVP, a novel approach which selects a minimized number of SAV-harboring protein sequences generated from the exhaustive approach, while still accommodating all possible variant peptides, by solving a classic set covering problem. Our study on known haplotype variations of TAS2R38 justifies the necessity for MinProtMaxVP to consider all combinations of SAVs. The performance of MinProtMaxVP is demonstrated by an in silico study on OR2T27 with five SAVs and real experimental data of the HEK293 cell line. Furthermore, assuming simulated somatic and germline variants of OR2T27 in tumor and normal tissues demonstrates that when adopting the appropriate somatic and germline SAV integration strategy, MinProtMaxVP is adaptable to labeling and label-free mass spectrometry-based experiments. SIGNIFICANCE: We present MinProtMaxVP, a novel approach to generate SAV-harboring protein sequences for constructing a customized protein sequence database, which is used in database searching for variant peptide identification. This approach outperforms the existing approaches in generating all possible variant peptides to be included in protein sequences and possibly leading to identification of more variant peptides in proteogenomic analysis.

Proteomics Pipeline for Identifying Variant Proteins in Plasmodium falciparum Parasites Isolated from Children Presenting with Malaria

Plasmodium falciparum variant antigens named erythrocyte membrane protein 1 (PfEMP1) are important targets for developing a protective immunity to malaria caused by P. falciparum. One of the major challenges in P. falciparum proteomics studies is identifying PfEMP1s at the protein level due to antigenic variation. To identify these PfEMP1s using shotgun proteomics, we developed a pipeline that searches high-resolution mass spectrometry spectra against a custom protein sequence database. A local alignment algorithm, LAX, was developed as a part of the pipeline that matches peptide sequences to the most similar PfEMP1 and calculates a weight value based on peptide's uniqueness used for PfEMP1 protein inference. The pipeline was first validated in the analysis of a laboratory strain with a known PfEMP1, then it was implemented on the analysis of parasite isolates from malaria-infected pregnant women and finally on the analysis of parasite isolates from malaria-infected children where there was an increase of PfEMP1s identified in 27 out of 31 isolates using the expanded database.

The Effect of Bacterial Peptide p28 on Viability and Apoptosis Status of P53-null HeLa Cells

Purpose: Despite all the efforts for discovery of efficient anti-cancer therapeutics, cancer is still a major health concern worldwide. p28 is a bacterial small peptide which has been widely investigated due to its preferential cell internalization and anti-cancer activities. Intracellularly, p28 offers its anti-cancer traits by impeding the degradation of tumor-suppressor protein "p53". In this study, we investigated the potency of p28 in inducing apoptosis or decreasing cell viability in p53-null "HeLa" cell line. Methods: The coding sequence for p28 peptide was obtained from Pseudomonas aeruginosa by PCR amplification of the p28 gene. The coding gene was cloned in pET-28a vector and transformed into E. coli bacterial host. Subsequently, the expressed peptide was purified using Ni-NTA chromatography system and introduced into the target cells. The anti-proliferative and apoptotic activity of p28 on HeLa and HEK-293 cells were investigated using MTT and PEAnnexin V Flowcytometry assays. Results: Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and Western blotting confirmed the expression of p28 peptide in the bacterial host. Bradford assay revealed a concentration of 0.05 mg/mL for the purified p28. MTT assay of cells treated with p28 at concentrations of 0, 0.5, 1, 2 and 2.5 µM indicated 24h viability values of 97%, 89%, 88%, 87% and 84% for HeLa cells, respectively. Data obtained from flowcytometry analyses revealed 24h apoptosis rate of 7.17%, 8.05%, 8.63% and 8.84% for HeLa cells treated with 0, 0.5, 1, and 2 µM p28, respectively. Conclusion: MTT and flowcytometry apoptosis assays suggest no statistically significant effect of p28 on the viability and apoptosis status of p53-null HeLa cells when results compared to data obtained from HEK-293 cells (P>0.05). These results imply that anti-cancer efficacy of p28 is directly dependent on the presence of p53, suggesting p28 as an inappropriate therapeutic agent for treatment of cancers with negative p53 status.

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.

Integration of large-scale multi-omic datasets: a protein-centric view

Innovative mass spectrometry-based proteomics has enabled routine measurements of protein abundance, localization, interactions, and modifications, covering unique aspects of gene expression regulation and function. It is now time to move from isolated analyses of these datasets toward true integration of proteomics with other data types to gain insights from the interactions and interdependencies of biomolecules. When combined with genomic or transcriptomic data, proteomics expands genome annotation to identify variant or missing genes. Dynamic proteomic measurements can move analysis from predominantly concentration-based framework to that of synthesis and degradation of proteins. Proteomic data from thousands of cancer patients can foster identification of novel pathogenic mutations via detection of protein sequence changes that lead to dysregulated pathways in various tumors. Such comprehensive efforts can exploit the synergy arising from large and complex datasets to advance virtually every field of biology.

Evaluating the Possibility of Detecting Variants in Shotgun Proteomics via LeTE-Fusion Analysis Pipeline

In proteogenomic studies, many genome-annotated events, for example, single amino acid variation (SAAV) and short INDEL, are often unobserved in shotgun proteomics. Therefore, we propose an analysis pipeline called LeTE-fusion (Le, peptide length; T, theoretical values; E, experimental data) to first investigate whether peptides with certain lengths are observed more often in mass spectrometry (MS)-based proteomics, which may hinder peptide identification causing difficulty in detecting genome-annotated events. By applying LeTE-fusion on different MS-based proteome data sets, we found peptides within 7-20 amino acids are more frequently identified, possibly attributed to MS-related factors instead of proteases. We then further extended the usage of LeTE-fusion on four variant-containing-sequence data sets (SAAV-only) with various sample complexity up to the whole human proteome scale, which yields theoretically ∼70% variants observable in an ideal shotgun proteomics. However, only ∼40% of variants might be detectable in real shotgun proteomic experiments when LeTE-fusion utilizes the experimentally observed variant-site-containing wild-type peptides in PeptideAtlas to estimate the expected observable coverage of variants. Finally, we conducted a case study on HEK293 cell line with variants reported at genomic level that were also identified in shotgun proteomics to demonstrate the efficacy of LeTE-fusion on estimating expected observable coverage of variants. To the best of our knowledge, this is the first study to systematically investigate the detection limits of genome-annotated events via shotgun proteomics using such analysis pipeline.

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.

Proteogenomic Investigation of Strain Variation in Clinical Mycobacterium tuberculosis Isolates

Mycobacterium tuberculosis consists of a large number of different strains that display unique virulence characteristics. Whole-genome sequencing has revealed substantial genetic diversity among clinical M. tuberculosis isolates, and elucidating the phenotypic variation encoded by this genetic diversity will be of the utmost importance to fully understand M. tuberculosis biology and pathogenicity. In this study, we integrated whole-genome sequencing and mass spectrometry (GeLC-MS/MS) to reveal strain-specific characteristics in the proteomes of two clinical M. tuberculosis Latin American-Mediterranean isolates. Using this approach, we identified 59 peptides containing single amino acid variants, which covered ∼9% of all coding nonsynonymous single nucleotide variants detected by whole-genome sequencing. Furthermore, we identified 29 distinct peptides that mapped to a hypothetical protein not present in the M. tuberculosis H37Rv reference proteome. Here, we provide evidence for the expression of this protein in the clinical M. tuberculosis SAWC3651 isolate. The strain-specific databases enabled confirmation of genomic differences (i.e., large genomic regions of difference and nonsynonymous single nucleotide variants) in these two clinical M. tuberculosis isolates and allowed strain differentiation at the proteome level. Our results contribute to the growing field of clinical microbial proteogenomics and can improve our understanding of phenotypic variation in clinical M. tuberculosis isolates.

ADAR-mediated messenger RNA editing: analysis at the proteome level

Post-transcriptional RNA editing by RNA specific adenosine deaminases (ADAR) was discovered more than two decades ago. It provides additional regulation of animal and human transcriptome. In most cases, it occurs in nervous tissue, where, as a result of the reaction, adenosine is converted to inosine in particular sites of RNA. In case of messenger RNA, during translation, inosine is recognized as guanine leading to amino acid substitutions. Those substitutions are shown to affect substantially the function of proteins, e.g. subunits of the glutamate receptor. Nevertheless, most of the works on RNA editing use analysis of nucleic acids, even those which deal with a coding RNA. In this review, we propose the use of shotgun proteomics based on high resolution liquid chromatography and mass spectrometry for investigation of the effects of RNA editing at the protein level. Recently developed methods of big data processing allow combining the results of various omics techniques, being referred to as proteogenomics. The proposed proteogenomic approach for the analysis of RNA editing at the protein level will directly conduct a qualitative and quantitative analysis of protein edited sequences in the scale of whole proteome.