Genomic variability and protein species — Improving sequence coverage for proteogenomics

Accurate Prediction of Protein Sequences for Proteogenomics Data Integration

This book chapter discusses proteogenomics data integration and provides an overview into the different omics layer involved in defining the proteome of a living organism. Various aspects of genome variability affecting either the sequence or abundance level of proteins are discussed in this book chapter, such as the effect of single-nucleotide variants or larger genomic structural variants on the proteome. Next, various sequencing technologies are introduced and discussed from a proteogenomics data integration perspective such as those providing short- and long-read sequencing and listing their respective advantages and shortcomings for accurate protein variant prediction using genomic/transcriptomics sequencing data. Finally, the various bioinformatics tools used to process and analyze DNA/RNA sequencing data are discussed with the ultimate goal of obtaining accurately predicted sample-specific protein sequences that can be used as a drop-in replacement in existing approaches for peptide and protein identification using popular database search engines such as MSFragger, SearchGUI/PeptideShaker.

Evaluation of bottom-up and top-down mass spectrum identifications with different customized protein sequences databases

MOTIVATION

Generally, bottom-up and top-down are two complementary approaches for proteoforms identification. The inference of proteoforms relies on searching mass spectra against an accurate proteoform sequence database. A customized protein sequence database derived by RNA-Seq data can be used to better identify the proteoform existed in a studied species. However, the quality of sequences in customized databases which constructed by different strategies affect the performances of mass spectrometry (MS) identification. Additionally, performances of identifications between bottom-up and top-down using customized databases are also needed to be evaluated.

RESULTS

Three customized databases were constructed with different strategies separately. Two of them were based on translating assembled transcripts with or without genomic annotation, and the third one is a variant-extending protein database. By testing with bottom-up and top-down MS data separately, a variant-extending protein database could identify not only the most number of spectra but also the alleles expressed at the same time in diploid cells. An assembled database could identify the spectrum missed in reference database and amino acid (AA) alterations existed in studied species.

AVAILABILITY AND IMPLEMENTATION

Experimental results demonstrated that the proteoform sequences in an annotated database are more suitable for identifying AA alterations and peptide sequences missed in reference database. An unannotated database instead of a reference proteome database gets an enough high sensitivity of identifying mass spectra. The variant-extending reference database is the most sensitive to identify mass spectra and single AA variants.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Application of Displacement Chromatography to Online Two-Dimensional Liquid Chromatography Coupled to Tandem Mass Spectrometry Improves Peptide Separation Efficiency and Detectability for the Analysis of Complex Proteomes

The complexity of mammalian proteomes is a challenge in bottom-up proteomics. For a comprehensive proteome analysis, multidimensional separation strategies are necessary. Online two-dimensional liquid chromatography–tandem mass spectrometry (2D-LC-MS/MS) combining strong cation exchange (SCX) in the first dimension with reversed-phase (RP) chromatography in the second dimension provides a powerful approach to analyze complex proteomes. Although the combination of SCX with RP chromatography provides a good orthogonality, only a moderate separation is achieved in the first dimension for peptides with two (+2) or three (+3) positive charges. The aim of this study was to improve the performance of online SCX-RP-MS/MS by applying displacement chromatography to the first separation dimension. Compared to gradient chromatography mode (GCM), displacement chromatography mode (DCM) was expected to improve the separation of +2-peptides and +3-peptides, thus reducing complexity and increasing ionization and detectability. The results show that DCM provided a separation of +2-peptides and +3-peptides in remarkably sharp zones with a low degree of coelution, thus providing fractions with significantly higher purities compared to GCM. In particular, +2-peptides were separated over several fractions, which was not possible to achieve in GCM. The better separation in DCM resulted in a higher reproducibility and significantly higher identification rates for both peptides and proteins including a 2.6-fold increase for +2-peptides. The higher number of identified peptides in DCM resulted in significantly higher protein sequence coverages and a considerably higher number of unique peptides per protein. Compared to conventionally used salt-based GCM, DCM increased the performance of online SCX-RP-MS/MS and enabled comprehensive proteome profiling in the low microgram range.

Proteogenomics: Key Driver for Clinical Discovery and Personalized Medicine

Proteogenomics is a multi-omics research field that has the aim to efficiently integrate genomics, transcriptomics and proteomics. With this approach it is possible to identify new patient-specific proteoforms that may have implications in disease development, specifically in cancer. Understanding the impact of a large number of mutations detected at the genomics level is needed to assess the effects at the proteome level. Proteogenomics data integration would help in identifying molecular changes that are persistent across multiple molecular layers and enable better interpretation of molecular mechanisms of disease, such as the causal relationship between single nucleotide polymorphisms (SNPs) and the expression of transcripts and translation of proteins compared to mainstream proteomics approaches. Identifying patient-specific protein forms and getting a better picture of molecular mechanisms of disease opens the avenue for precision and personalized medicine. Proteogenomics is, however, a challenging interdisciplinary science that requires the understanding of sample preparation, data acquisition and processing for genomics, transcriptomics and proteomics. This chapter aims to guide the reader through the technology and bioinformatics aspects of these multi-omics approaches, illustrated with proteogenomics applications having clinical or biological relevance.

Control of Gene Expression by RNA Binding Protein Action on Alternative Translation Initiation Sites

Transcript levels do not faithfully predict protein levels, due to post-transcriptional regulation of gene expression mediated by RNA binding proteins (RBPs) and non-coding RNAs. We developed a multivariate linear regression model integrating RBP levels and predicted RBP-mRNA regulatory interactions from matched transcript and protein datasets. RBPs significantly improved the accuracy in predicting protein abundance of a portion of the total modeled mRNAs in three panels of tissues and cells and for different methods employed in the detection of mRNA and protein. The presence of upstream translation initiation sites (uTISs) at the mRNA 5’ untranslated regions was strongly associated with improvement in predictive accuracy. On the basis of these observations, we propose that the recently discovered widespread uTISs in the human genome can be a previously unappreciated substrate of translational control mediated by RBPs.

Gene expression is a dynamic program by which the information stored in the genome is rendered functional by production and degradation of two types of macromolecules, RNAs and proteins. mRNAs are templates for proteins; therefore we expect correspondence between quantities of mRNAs and proteins. Genome-wide studies instead indicate a marked discrepancy between them, when considering their steady-state levels or their variations across different conditions. We employed linear regression approaches with paired mRNA/protein datasets in order to develop a model predicting the protein level of a gene from both the mRNA level and the protein levels of RBPs inferred to bind the mRNA untranslated regions. The results of our analyses restricted the utility of RBPs to improve accuracy of predicted protein abundance to a small fraction of the total modelled genes, and identified a novel association of the improvement induced by RBPs with the presence of upstream translation sites. This finding suggests a new avenue of experimental studies aimed at exploring the hypothesis that RBPs could influence protein abundance by changing the preference for certain translation initiation sites.